← Back to branch summary

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

« back to all changes in this revision

Viewing changes to tests/mphase_inlet_velocity_bc_incompressible/mphase_inlet_velocity_bc_incompressible.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/mphase_inlet_velocity_bc_incompressible/mphase_inlet_velocity_bc_incompressible.flml
 
1
<?xml version='1.0' encoding='utf-8'?>
 
2
<fluidity_options>
 
3
  <simulation_name>
 
4
    <string_value lines="1">mphase_inlet_velocity_bc_incompressible</string_value>
 
5
  </simulation_name>
 
6
  <problem_type>
 
7
    <string_value lines="1">multiphase</string_value>
 
8
  </problem_type>
 
9
  <geometry>
 
10
    <dimension>
 
11
      <integer_value rank="0">1</integer_value>
 
12
    </dimension>
 
13
    <mesh name="CoordinateMesh">
 
14
      <from_file file_name="line">
 
15
        <format name="triangle"/>
 
16
        <stat>
 
17
          <include_in_stat/>
 
18
        </stat>
 
19
      </from_file>
 
20
    </mesh>
 
21
    <mesh name="VelocityMesh">
 
22
      <from_mesh>
 
23
        <mesh name="CoordinateMesh"/>
 
24
        <mesh_shape>
 
25
          <polynomial_degree>
 
26
            <integer_value rank="0">2</integer_value>
 
27
          </polynomial_degree>
 
28
        </mesh_shape>
 
29
        <stat>
 
30
          <exclude_from_stat/>
 
31
        </stat>
 
32
      </from_mesh>
 
33
    </mesh>
 
34
    <mesh name="PressureMesh">
 
35
      <from_mesh>
 
36
        <mesh name="CoordinateMesh"/>
 
37
        <mesh_shape>
 
38
          <polynomial_degree>
 
39
            <integer_value rank="0">1</integer_value>
 
40
          </polynomial_degree>
 
41
        </mesh_shape>
 
42
        <stat>
 
43
          <exclude_from_stat/>
 
44
        </stat>
 
45
      </from_mesh>
 
46
    </mesh>
 
47
    <quadrature>
 
48
      <degree>
 
49
        <integer_value rank="0">4</integer_value>
 
50
      </degree>
 
51
    </quadrature>
 
52
  </geometry>
 
53
  <io>
 
54
    <dump_format>
 
55
      <string_value>vtk</string_value>
 
56
    </dump_format>
 
57
    <dump_period>
 
58
      <constant>
 
59
        <real_value rank="0">0</real_value>
 
60
      </constant>
 
61
    </dump_period>
 
62
    <output_mesh name="VelocityMesh"/>
 
63
    <stat>
 
64
      <output_at_start/>
 
65
    </stat>
 
66
  </io>
 
67
  <timestepping>
 
68
    <current_time>
 
69
      <real_value rank="0">0</real_value>
 
70
    </current_time>
 
71
    <timestep>
 
72
      <real_value rank="0">0.01</real_value>
 
73
    </timestep>
 
74
    <finish_time>
 
75
      <real_value rank="0">10.0</real_value>
 
76
    </finish_time>
 
77
    <nonlinear_iterations>
 
78
      <integer_value rank="0">5</integer_value>
 
79
      <tolerance>
 
80
        <real_value rank="0">1.0e-9</real_value>
 
81
        <infinity_norm/>
 
82
      </tolerance>
 
83
    </nonlinear_iterations>
 
84
    <steady_state>
 
85
      <tolerance>
 
86
        <real_value rank="0">1.0e-9</real_value>
 
87
        <infinity_norm/>
 
88
      </tolerance>
 
89
    </steady_state>
 
90
  </timestepping>
 
91
  <material_phase name="Gas">
 
92
    <equation_of_state>
 
93
      <fluids>
 
94
        <linear>
 
95
          <reference_density>
 
96
            <real_value rank="0">1.23</real_value>
 
97
          </reference_density>
 
98
        </linear>
 
99
      </fluids>
 
100
    </equation_of_state>
 
101
    <scalar_field name="Pressure" rank="0">
 
102
      <prognostic>
 
103
        <mesh name="PressureMesh"/>
 
104
        <spatial_discretisation>
 
105
          <continuous_galerkin>
 
106
            <remove_stabilisation_term/>
 
107
            <integrate_continuity_by_parts/>
 
108
          </continuous_galerkin>
 
109
        </spatial_discretisation>
 
110
        <scheme>
 
111
          <poisson_pressure_solution>
 
112
            <string_value lines="1">never</string_value>
 
113
          </poisson_pressure_solution>
 
114
          <use_projection_method/>
 
115
        </scheme>
 
116
        <solver>
 
117
          <iterative_method name="preonly"/>
 
118
          <preconditioner name="lu">
 
119
            <factorization_package name="petsc"/>
 
120
          </preconditioner>
 
121
          <relative_error>
 
122
            <real_value rank="0">1.0e-7</real_value>
 
123
          </relative_error>
 
124
          <max_iterations>
 
125
            <integer_value rank="0">1000</integer_value>
 
126
          </max_iterations>
 
127
          <never_ignore_solver_failures/>
 
128
          <diagnostics>
 
129
            <monitors/>
 
130
          </diagnostics>
 
131
        </solver>
 
132
        <output/>
 
133
        <stat/>
 
134
        <convergence>
 
135
          <include_in_convergence/>
 
136
        </convergence>
 
137
        <detectors>
 
138
          <exclude_from_detectors/>
 
139
        </detectors>
 
140
        <steady_state>
 
141
          <include_in_steady_state/>
 
142
        </steady_state>
 
143
        <no_interpolation/>
 
144
      </prognostic>
 
145
    </scalar_field>
 
146
    <scalar_field name="Density" rank="0">
 
147
      <diagnostic>
 
148
        <algorithm name="Internal" material_phase_support="multiple"/>
 
149
        <mesh name="CoordinateMesh"/>
 
150
        <output/>
 
151
        <stat/>
 
152
        <convergence>
 
153
          <include_in_convergence/>
 
154
        </convergence>
 
155
        <detectors>
 
156
          <include_in_detectors/>
 
157
        </detectors>
 
158
        <steady_state>
 
159
          <include_in_steady_state/>
 
160
        </steady_state>
 
161
      </diagnostic>
 
162
    </scalar_field>
 
163
    <vector_field name="Velocity" rank="1">
 
164
      <prognostic>
 
165
        <mesh name="VelocityMesh"/>
 
166
        <equation name="LinearMomentum"/>
 
167
        <spatial_discretisation>
 
168
          <continuous_galerkin>
 
169
            <stabilisation>
 
170
              <no_stabilisation/>
 
171
            </stabilisation>
 
172
            <mass_terms>
 
173
              <lump_mass_matrix/>
 
174
            </mass_terms>
 
175
            <advection_terms/>
 
176
            <stress_terms>
 
177
              <tensor_form/>
 
178
            </stress_terms>
 
179
          </continuous_galerkin>
 
180
          <conservative_advection>
 
181
            <real_value rank="0">0</real_value>
 
182
          </conservative_advection>
 
183
        </spatial_discretisation>
 
184
        <temporal_discretisation>
 
185
          <theta>
 
186
            <real_value rank="0">1.0</real_value>
 
187
          </theta>
 
188
          <relaxation>
 
189
            <real_value rank="0">0.5</real_value>
 
190
          </relaxation>
 
191
        </temporal_discretisation>
 
192
        <solver>
 
193
          <iterative_method name="preonly"/>
 
194
          <preconditioner name="lu">
 
195
            <factorization_package name="petsc"/>
 
196
          </preconditioner>
 
197
          <relative_error>
 
198
            <real_value rank="0">1.0e-7</real_value>
 
199
          </relative_error>
 
200
          <max_iterations>
 
201
            <integer_value rank="0">1000</integer_value>
 
202
          </max_iterations>
 
203
          <never_ignore_solver_failures/>
 
204
          <diagnostics>
 
205
            <monitors/>
 
206
          </diagnostics>
 
207
        </solver>
 
208
        <initial_condition name="WholeMesh">
 
209
          <constant>
 
210
            <real_value shape="1" dim1="dim" rank="1">0.0</real_value>
 
211
          </constant>
 
212
        </initial_condition>
 
213
        <boundary_conditions name="Left">
 
214
          <surface_ids>
 
215
            <integer_value shape="1" rank="1">1</integer_value>
 
216
          </surface_ids>
 
217
          <type name="dirichlet">
 
218
            <align_bc_with_cartesian>
 
219
              <x_component>
 
220
                <constant>
 
221
                  <real_value rank="0">0.1</real_value>
 
222
                </constant>
 
223
              </x_component>
 
224
            </align_bc_with_cartesian>
 
225
          </type>
 
226
        </boundary_conditions>
 
227
        <boundary_conditions name="Right">
 
228
          <surface_ids>
 
229
            <integer_value shape="1" rank="1">2</integer_value>
 
230
          </surface_ids>
 
231
          <type name="dirichlet">
 
232
            <align_bc_with_cartesian>
 
233
              <x_component>
 
234
                <constant>
 
235
                  <real_value rank="0">0</real_value>
 
236
                </constant>
 
237
              </x_component>
 
238
            </align_bc_with_cartesian>
 
239
          </type>
 
240
        </boundary_conditions>
 
241
        <tensor_field name="Viscosity" rank="2">
 
242
          <prescribed>
 
243
            <value name="WholeMesh">
 
244
              <isotropic>
 
245
                <constant>
 
246
                  <real_value rank="0">1.78e-5</real_value>
 
247
                </constant>
 
248
              </isotropic>
 
249
            </value>
 
250
            <output/>
 
251
          </prescribed>
 
252
        </tensor_field>
 
253
        <output/>
 
254
        <stat>
 
255
          <include_in_stat/>
 
256
          <previous_time_step>
 
257
            <exclude_from_stat/>
 
258
          </previous_time_step>
 
259
          <nonlinear_field>
 
260
            <exclude_from_stat/>
 
261
          </nonlinear_field>
 
262
        </stat>
 
263
        <convergence>
 
264
          <include_in_convergence/>
 
265
        </convergence>
 
266
        <detectors>
 
267
          <include_in_detectors/>
 
268
        </detectors>
 
269
        <steady_state>
 
270
          <include_in_steady_state/>
 
271
        </steady_state>
 
272
        <consistent_interpolation/>
 
273
      </prognostic>
 
274
    </vector_field>
 
275
    <scalar_field name="PhaseVolumeFraction" rank="0">
 
276
      <prescribed>
 
277
        <mesh name="CoordinateMesh"/>
 
278
        <value name="WholeMesh">
 
279
          <constant>
 
280
            <real_value rank="0">0.5</real_value>
 
281
          </constant>
 
282
        </value>
 
283
        <output/>
 
284
        <stat/>
 
285
        <detectors>
 
286
          <exclude_from_detectors/>
 
287
        </detectors>
 
288
      </prescribed>
 
289
    </scalar_field>
 
290
    <scalar_field name="SumVelocityDivergence" rank="0">
 
291
      <diagnostic>
 
292
        <mesh name="PressureMesh"/>
 
293
        <algorithm name="Internal" material_phase_support="multiple"/>
 
294
        <integrate_divergence_by_parts/>
 
295
        <output/>
 
296
        <stat/>
 
297
        <detectors>
 
298
          <include_in_detectors/>
 
299
        </detectors>
 
300
        <solver>
 
301
          <iterative_method name="gmres">
 
302
            <restart>
 
303
              <integer_value rank="0">30</integer_value>
 
304
            </restart>
 
305
          </iterative_method>
 
306
          <preconditioner name="sor"/>
 
307
          <relative_error>
 
308
            <real_value rank="0">1.0e-7</real_value>
 
309
          </relative_error>
 
310
          <max_iterations>
 
311
            <integer_value rank="0">1000</integer_value>
 
312
          </max_iterations>
 
313
          <never_ignore_solver_failures/>
 
314
          <diagnostics>
 
315
            <monitors/>
 
316
          </diagnostics>
 
317
        </solver>
 
318
      </diagnostic>
 
319
    </scalar_field>
 
320
    <scalar_field name="FiniteElementDivergence" rank="0">
 
321
      <diagnostic field_name="Velocity">
 
322
        <algorithm legacy="true" name="Internal" material_phase_support="multiple"/>
 
323
        <mesh name="PressureMesh"/>
 
324
        <integrate_divergence_by_parts/>
 
325
        <solver>
 
326
          <iterative_method name="gmres">
 
327
            <restart>
 
328
              <integer_value rank="0">30</integer_value>
 
329
            </restart>
 
330
          </iterative_method>
 
331
          <preconditioner name="sor"/>
 
332
          <relative_error>
 
333
            <real_value rank="0">1.0e-7</real_value>
 
334
          </relative_error>
 
335
          <max_iterations>
 
336
            <integer_value rank="0">1000</integer_value>
 
337
          </max_iterations>
 
338
          <never_ignore_solver_failures/>
 
339
          <diagnostics>
 
340
            <monitors/>
 
341
          </diagnostics>
 
342
        </solver>
 
343
        <output/>
 
344
        <stat/>
 
345
        <convergence>
 
346
          <include_in_convergence/>
 
347
        </convergence>
 
348
        <detectors>
 
349
          <include_in_detectors/>
 
350
        </detectors>
 
351
        <steady_state>
 
352
          <include_in_steady_state/>
 
353
        </steady_state>
 
354
      </diagnostic>
 
355
    </scalar_field>
 
356
  </material_phase>
 
357
  <material_phase name="Particle">
 
358
    <equation_of_state>
 
359
      <fluids>
 
360
        <linear>
 
361
          <reference_density>
 
362
            <real_value rank="0">1.23</real_value>
 
363
          </reference_density>
 
364
        </linear>
 
365
      </fluids>
 
366
    </equation_of_state>
 
367
    <scalar_field name="Pressure" rank="0">
 
368
      <aliased material_phase_name="Gas" field_name="Pressure"/>
 
369
    </scalar_field>
 
370
    <scalar_field name="Density" rank="0">
 
371
      <diagnostic>
 
372
        <algorithm name="Internal" material_phase_support="multiple"/>
 
373
        <mesh name="CoordinateMesh"/>
 
374
        <output/>
 
375
        <stat/>
 
376
        <convergence>
 
377
          <include_in_convergence/>
 
378
        </convergence>
 
379
        <detectors>
 
380
          <include_in_detectors/>
 
381
        </detectors>
 
382
        <steady_state>
 
383
          <include_in_steady_state/>
 
384
        </steady_state>
 
385
      </diagnostic>
 
386
    </scalar_field>
 
387
    <vector_field name="Velocity" rank="1">
 
388
      <prognostic>
 
389
        <mesh name="VelocityMesh"/>
 
390
        <equation name="LinearMomentum"/>
 
391
        <spatial_discretisation>
 
392
          <continuous_galerkin>
 
393
            <stabilisation>
 
394
              <no_stabilisation/>
 
395
            </stabilisation>
 
396
            <mass_terms>
 
397
              <lump_mass_matrix/>
 
398
            </mass_terms>
 
399
            <advection_terms/>
 
400
            <stress_terms>
 
401
              <tensor_form/>
 
402
            </stress_terms>
 
403
          </continuous_galerkin>
 
404
          <conservative_advection>
 
405
            <real_value rank="0">0</real_value>
 
406
          </conservative_advection>
 
407
        </spatial_discretisation>
 
408
        <temporal_discretisation>
 
409
          <theta>
 
410
            <real_value rank="0">1.0</real_value>
 
411
          </theta>
 
412
          <relaxation>
 
413
            <real_value rank="0">0.5</real_value>
 
414
          </relaxation>
 
415
        </temporal_discretisation>
 
416
        <solver>
 
417
          <iterative_method name="preonly"/>
 
418
          <preconditioner name="lu">
 
419
            <factorization_package name="petsc"/>
 
420
          </preconditioner>
 
421
          <relative_error>
 
422
            <real_value rank="0">1.0e-7</real_value>
 
423
          </relative_error>
 
424
          <max_iterations>
 
425
            <integer_value rank="0">1000</integer_value>
 
426
          </max_iterations>
 
427
          <never_ignore_solver_failures/>
 
428
          <diagnostics>
 
429
            <monitors/>
 
430
          </diagnostics>
 
431
        </solver>
 
432
        <initial_condition name="WholeMesh">
 
433
          <constant>
 
434
            <real_value shape="1" dim1="dim" rank="1">0.0</real_value>
 
435
          </constant>
 
436
        </initial_condition>
 
437
        <boundary_conditions name="Right">
 
438
          <surface_ids>
 
439
            <integer_value shape="1" rank="1">2</integer_value>
 
440
          </surface_ids>
 
441
          <type name="dirichlet">
 
442
            <align_bc_with_cartesian>
 
443
              <x_component>
 
444
                <constant>
 
445
                  <real_value rank="0">0</real_value>
 
446
                </constant>
 
447
              </x_component>
 
448
            </align_bc_with_cartesian>
 
449
          </type>
 
450
        </boundary_conditions>
 
451
        <tensor_field name="Viscosity" rank="2">
 
452
          <prescribed>
 
453
            <value name="WholeMesh">
 
454
              <isotropic>
 
455
                <constant>
 
456
                  <real_value rank="0">1.78e-5</real_value>
 
457
                </constant>
 
458
              </isotropic>
 
459
            </value>
 
460
            <output/>
 
461
          </prescribed>
 
462
        </tensor_field>
 
463
        <output/>
 
464
        <stat>
 
465
          <include_in_stat/>
 
466
          <previous_time_step>
 
467
            <exclude_from_stat/>
 
468
          </previous_time_step>
 
469
          <nonlinear_field>
 
470
            <exclude_from_stat/>
 
471
          </nonlinear_field>
 
472
        </stat>
 
473
        <convergence>
 
474
          <include_in_convergence/>
 
475
        </convergence>
 
476
        <detectors>
 
477
          <include_in_detectors/>
 
478
        </detectors>
 
479
        <steady_state>
 
480
          <include_in_steady_state/>
 
481
        </steady_state>
 
482
        <consistent_interpolation/>
 
483
      </prognostic>
 
484
    </vector_field>
 
485
    <scalar_field name="PhaseVolumeFraction" rank="0">
 
486
      <diagnostic>
 
487
        <mesh name="CoordinateMesh"/>
 
488
        <algorithm name="Internal" material_phase_support="multiple"/>
 
489
        <output/>
 
490
        <stat/>
 
491
        <detectors>
 
492
          <include_in_detectors/>
 
493
        </detectors>
 
494
      </diagnostic>
 
495
    </scalar_field>
 
496
    <scalar_field name="FiniteElementDivergence" rank="0">
 
497
      <diagnostic field_name="Velocity">
 
498
        <algorithm legacy="true" name="Internal" material_phase_support="multiple"/>
 
499
        <mesh name="PressureMesh"/>
 
500
        <integrate_divergence_by_parts/>
 
501
        <solver>
 
502
          <iterative_method name="gmres">
 
503
            <restart>
 
504
              <integer_value rank="0">30</integer_value>
 
505
            </restart>
 
506
          </iterative_method>
 
507
          <preconditioner name="sor"/>
 
508
          <relative_error>
 
509
            <real_value rank="0">1.0e-7</real_value>
 
510
          </relative_error>
 
511
          <max_iterations>
 
512
            <integer_value rank="0">1000</integer_value>
 
513
          </max_iterations>
 
514
          <never_ignore_solver_failures/>
 
515
          <diagnostics>
 
516
            <monitors/>
 
517
          </diagnostics>
 
518
        </solver>
 
519
        <output/>
 
520
        <stat/>
 
521
        <convergence>
 
522
          <include_in_convergence/>
 
523
        </convergence>
 
524
        <detectors>
 
525
          <include_in_detectors/>
 
526
        </detectors>
 
527
        <steady_state>
 
528
          <include_in_steady_state/>
 
529
        </steady_state>
 
530
      </diagnostic>
 
531
    </scalar_field>
 
532
  </material_phase>
 
533
</fluidity_options>