← Back to branch summary

~fluidity-core/fluidity/sea-ice-branch

« back to all changes in this revision

Viewing changes to tests/driven_cavity_2d/driven_cavity_2d.flml

  • Committer: Simon Mouradian
  • Date: 2012-10-19 10:35:59 UTC
  • mfrom: (3520.32.371 fluidity)
  • Revision ID: simon.mouradian06@imperial.ac.uk-20121019103559-y36qa47phc69q8sc
merge from trunk

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/driven_cavity_2d/driven_cavity_2d.flml
4
4
    <string_value lines="1">driven_cavity_2d</string_value>
5
5
  </simulation_name>
6
6
  <problem_type>
7
 
    <string_value lines="1">oceans</string_value>
 
7
    <string_value lines="1">fluids</string_value>
8
8
  </problem_type>
9
9
  <geometry>
10
10
    <dimension>
36
36
    </mesh>
37
37
    <quadrature>
38
38
      <degree>
39
 
        <integer_value rank="0">4</integer_value>
 
39
        <integer_value rank="0">3</integer_value>
40
40
      </degree>
41
41
    </quadrature>
42
42
  </geometry>
46
46
    </dump_format>
47
47
    <dump_period>
48
48
      <constant>
49
 
        <real_value rank="0">10.0</real_value>
 
49
        <real_value rank="0">5.0</real_value>
50
50
      </constant>
51
51
    </dump_period>
52
52
    <output_mesh name="VelocityMesh"/>
62
62
    <finish_time>
63
63
      <real_value rank="0">60.0</real_value>
64
64
    </finish_time>
 
65
    <nonlinear_iterations>
 
66
      <integer_value rank="0">2</integer_value>
 
67
    </nonlinear_iterations>
65
68
  </timestepping>
66
69
  <physical_parameters/>
67
70
  <material_phase name="Fluid">
 
71
    <equation_of_state>
 
72
      <fluids>
 
73
        <linear>
 
74
          <reference_density>
 
75
            <real_value rank="0">1.0</real_value>
 
76
          </reference_density>
 
77
        </linear>
 
78
      </fluids>
 
79
    </equation_of_state>
68
80
    <scalar_field name="Pressure" rank="0">
69
81
      <prognostic>
70
82
        <mesh name="PressureMesh"/>
72
84
          <continuous_galerkin/>
73
85
        </spatial_discretisation>
74
86
        <reference_node>
75
 
          <integer_value rank="0">1</integer_value>
 
87
          <integer_value rank="0">10</integer_value>
76
88
        </reference_node>
77
89
        <scheme>
78
90
          <poisson_pressure_solution>
82
94
        </scheme>
83
95
        <solver>
84
96
          <iterative_method name="cg"/>
85
 
          <preconditioner name="sor"/>
 
97
          <preconditioner name="mg"/>
86
98
          <relative_error>
87
 
            <real_value rank="0">1.0e-7</real_value>
 
99
            <real_value rank="0">1.0e-6</real_value>
88
100
          </relative_error>
89
101
          <max_iterations>
90
102
            <integer_value rank="0">1000</integer_value>
94
106
            <monitors/>
95
107
          </diagnostics>
96
108
        </solver>
97
 
        <initial_condition name="WholeMesh">
98
 
          <constant>
99
 
            <real_value rank="0">0</real_value>
100
 
          </constant>
101
 
        </initial_condition>
102
109
        <output/>
103
110
        <stat/>
104
111
        <convergence>
108
115
          <exclude_from_detectors/>
109
116
        </detectors>
110
117
        <steady_state>
 
118
          <exclude_from_steady_state/>
 
119
        </steady_state>
 
120
        <no_interpolation/>
 
121
      </prognostic>
 
122
    </scalar_field>
 
123
    <scalar_field name="Density" rank="0">
 
124
      <diagnostic>
 
125
        <algorithm name="Internal" material_phase_support="multiple"/>
 
126
        <mesh name="VelocityMesh"/>
 
127
        <output/>
 
128
        <stat/>
 
129
        <convergence>
 
130
          <include_in_convergence/>
 
131
        </convergence>
 
132
        <detectors>
 
133
          <include_in_detectors/>
 
134
        </detectors>
 
135
        <steady_state>
111
136
          <include_in_steady_state/>
112
137
        </steady_state>
113
 
        <consistent_interpolation/>
114
 
      </prognostic>
 
138
      </diagnostic>
115
139
    </scalar_field>
116
140
    <vector_field name="Velocity" rank="1">
117
141
      <prognostic>
118
142
        <mesh name="VelocityMesh"/>
119
 
        <equation name="Boussinesq"/>
 
143
        <equation name="LinearMomentum"/>
120
144
        <spatial_discretisation>
121
145
          <continuous_galerkin>
122
146
            <stabilisation>
150
174
          </iterative_method>
151
175
          <preconditioner name="sor"/>
152
176
          <relative_error>
153
 
            <real_value rank="0">1.0e-7</real_value>
 
177
            <real_value rank="0">1.0e-6</real_value>
154
178
          </relative_error>
155
179
          <max_iterations>
156
180
            <integer_value rank="0">1000</integer_value>
172
196
          <type name="dirichlet">
173
197
            <align_bc_with_cartesian>
174
198
              <x_component>
175
 
                <constant>
176
 
                  <real_value rank="0">1.0</real_value>
177
 
                </constant>
 
199
                <python>
 
200
                  <string_value lines="20" type="code" language="python">def val(X,t):
 
201
  if(X[0]&lt;1e-3):
 
202
    return 0.0
 
203
  elif(X[0]&gt;1.-1e-3):
 
204
    return 0.0
 
205
  else:
 
206
    return 1.0</string_value>
 
207
                </python>
178
208
              </x_component>
179
209
              <y_component>
180
210
                <constant>
270
300
          <include_in_detectors/>
271
301
        </detectors>
272
302
        <steady_state>
273
 
          <include_in_steady_state/>
274
 
        </steady_state>
275
 
      </diagnostic>
276
 
    </scalar_field>
277
 
    <scalar_field name="StreamFunction" rank="0">
 
303
          <exclude_from_steady_state/>
 
304
        </steady_state>
 
305
      </diagnostic>
 
306
    </scalar_field>
 
307
    <scalar_field name="Vorticity2D" rank="0">
 
308
      <diagnostic>
 
309
        <algorithm source_field_type="vector" material_phase_support="single" name="curl_2d" source_field_name="Velocity">
 
310
          <lump_mass/>
 
311
        </algorithm>
 
312
        <mesh name="VelocityMesh"/>
 
313
        <output/>
 
314
        <stat/>
 
315
        <convergence>
 
316
          <include_in_convergence/>
 
317
        </convergence>
 
318
        <detectors>
 
319
          <include_in_detectors/>
 
320
        </detectors>
 
321
        <steady_state>
 
322
          <exclude_from_steady_state/>
 
323
        </steady_state>
 
324
      </diagnostic>
 
325
    </scalar_field>
 
326
    <scalar_field depends="Velocity" name="MultiplyConnectedStreamFunction" rank="0">
278
327
      <prognostic>
279
 
        <mesh name="PressureMesh"/>
 
328
        <mesh name="CoordinateMesh"/>
280
329
        <solver>
281
 
          <iterative_method name="gmres">
282
 
            <restart>
283
 
              <integer_value rank="0">30</integer_value>
284
 
            </restart>
285
 
          </iterative_method>
286
 
          <preconditioner name="sor"/>
 
330
          <iterative_method name="cg"/>
 
331
          <preconditioner name="mg"/>
287
332
          <relative_error>
288
 
            <real_value rank="0">1.0E-7</real_value>
 
333
            <real_value rank="0">1.0e-6</real_value>
289
334
          </relative_error>
290
335
          <max_iterations>
291
336
            <integer_value rank="0">1000</integer_value>
295
340
            <monitors/>
296
341
          </diagnostics>
297
342
        </solver>
 
343
        <boundary_conditions name="bottom">
 
344
          <surface_ids>
 
345
            <integer_value shape="1" rank="1">444</integer_value>
 
346
          </surface_ids>
 
347
          <primary_boundary/>
 
348
          <type name="dirichlet">
 
349
            <internally_calculated/>
 
350
          </type>
 
351
        </boundary_conditions>
 
352
        <boundary_conditions name="top">
 
353
          <surface_ids>
 
354
            <integer_value shape="1" rank="1">333</integer_value>
 
355
          </surface_ids>
 
356
          <secondary_boundary>
 
357
            <primary_point>
 
358
              <real_value shape="2" dim1="dim" rank="1">0.5 0.0</real_value>
 
359
            </primary_point>
 
360
            <secondary_point>
 
361
              <real_value shape="2" dim1="dim" rank="1">0.5 1.0</real_value>
 
362
            </secondary_point>
 
363
          </secondary_boundary>
 
364
          <type name="dirichlet">
 
365
            <internally_calculated/>
 
366
          </type>
 
367
        </boundary_conditions>
298
368
        <exclude_from_checkpointing/>
299
369
        <output/>
300
370
        <stat/>
301
371
        <detectors>
302
372
          <include_in_detectors/>
303
373
        </detectors>
304
 
        <consistent_interpolation/>
 
374
        <steady_state>
 
375
          <include_in_steady_state/>
 
376
        </steady_state>
 
377
        <no_interpolation/>
305
378
      </prognostic>
306
379
    </scalar_field>
307
380
  </material_phase>