← Back to branch summary

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

« back to all changes in this revision

Viewing changes to tests/3material-droplet-p1dgp2-test-cty-cv-rhop1dg-3d/3material-droplet-p1dgp2-test-cty-cv-rhop1dg-3d.flml

  • Committer: Simon Mouradian
  • Date: 2012-03-30 11:11:47 UTC
  • mfrom: (3520.32.239 fluidity)
  • Revision ID: simon.mouradian06@imperial.ac.uk-20120330111147-y4n4y4th52x2dqcf
merge from lp:fluidity

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/3material-droplet-p1dgp2-test-cty-cv-rhop1dg-3d/3material-droplet-p1dgp2-test-cty-cv-rhop1dg-3d.flml
1
 
<?xml version='1.0' encoding='utf-8'?>
2
 
<fluidity_options>
3
 
  <simulation_name>
4
 
    <string_value lines="1">3material-droplet-p1dgp2-test-cty-cv-rhop1dg-3d</string_value>
5
 
  </simulation_name>
6
 
  <problem_type>
7
 
    <string_value lines="1">multimaterial</string_value>
8
 
  </problem_type>
9
 
  <geometry>
10
 
    <dimension>
11
 
      <integer_value rank="0">3</integer_value>
12
 
    </dimension>
13
 
    <mesh name="CoordinateMesh">
14
 
      <from_file file_name="cube">
15
 
        <format name="gmsh"/>
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">1</integer_value>
27
 
          </polynomial_degree>
28
 
        </mesh_shape>
29
 
        <mesh_continuity>
30
 
          <string_value>discontinuous</string_value>
31
 
        </mesh_continuity>
32
 
        <stat>
33
 
          <exclude_from_stat/>
34
 
        </stat>
35
 
      </from_mesh>
36
 
    </mesh>
37
 
    <mesh name="PressureMesh">
38
 
      <from_mesh>
39
 
        <mesh name="CoordinateMesh"/>
40
 
        <mesh_shape>
41
 
          <polynomial_degree>
42
 
            <integer_value rank="0">2</integer_value>
43
 
          </polynomial_degree>
44
 
        </mesh_shape>
45
 
        <stat>
46
 
          <exclude_from_stat/>
47
 
        </stat>
48
 
      </from_mesh>
49
 
    </mesh>
50
 
    <mesh name="OutputMesh">
51
 
      <from_mesh>
52
 
        <mesh name="CoordinateMesh"/>
53
 
        <mesh_shape>
54
 
          <polynomial_degree>
55
 
            <integer_value rank="0">2</integer_value>
56
 
          </polynomial_degree>
57
 
        </mesh_shape>
58
 
        <mesh_continuity>
59
 
          <string_value>discontinuous</string_value>
60
 
        </mesh_continuity>
61
 
        <stat>
62
 
          <exclude_from_stat/>
63
 
        </stat>
64
 
      </from_mesh>
65
 
    </mesh>
66
 
    <mesh name="DensityMesh">
67
 
      <from_mesh>
68
 
        <mesh name="CoordinateMesh"/>
69
 
        <mesh_shape>
70
 
          <polynomial_degree>
71
 
            <integer_value rank="0">1</integer_value>
72
 
          </polynomial_degree>
73
 
        </mesh_shape>
74
 
        <mesh_continuity>
75
 
          <string_value>discontinuous</string_value>
76
 
        </mesh_continuity>
77
 
        <stat>
78
 
          <exclude_from_stat/>
79
 
        </stat>
80
 
      </from_mesh>
81
 
    </mesh>
82
 
    <quadrature>
83
 
      <degree>
84
 
        <integer_value rank="0">5</integer_value>
85
 
      </degree>
86
 
      <controlvolume_surface_degree>
87
 
        <integer_value rank="0">4</integer_value>
88
 
      </controlvolume_surface_degree>
89
 
    </quadrature>
90
 
  </geometry>
91
 
  <io>
92
 
    <dump_format>
93
 
      <string_value>vtk</string_value>
94
 
    </dump_format>
95
 
    <dump_period_in_timesteps>
96
 
      <constant>
97
 
        <integer_value rank="0">10</integer_value>
98
 
      </constant>
99
 
    </dump_period_in_timesteps>
100
 
    <output_mesh name="OutputMesh"/>
101
 
    <stat>
102
 
      <output_at_start/>
103
 
    </stat>
104
 
  </io>
105
 
  <timestepping>
106
 
    <current_time>
107
 
      <real_value rank="0">0.0</real_value>
108
 
    </current_time>
109
 
    <timestep>
110
 
      <real_value rank="0">0.0025</real_value>
111
 
    </timestep>
112
 
    <finish_time>
113
 
      <real_value rank="0">200.0</real_value>
114
 
      <comment>200.0</comment>
115
 
    </finish_time>
116
 
    <final_timestep>
117
 
      <integer_value rank="0">50</integer_value>
118
 
    </final_timestep>
119
 
    <nonlinear_iterations>
120
 
      <integer_value rank="0">2</integer_value>
121
 
    </nonlinear_iterations>
122
 
    <adaptive_timestep>
123
 
      <requested_cfl>
124
 
        <real_value rank="0">1.0</real_value>
125
 
      </requested_cfl>
126
 
      <courant_number name="ControlVolumeCFLNumber">
127
 
        <mesh name="PressureMesh"/>
128
 
      </courant_number>
129
 
      <increase_tolerance>
130
 
        <real_value rank="0">1.05</real_value>
131
 
      </increase_tolerance>
132
 
    </adaptive_timestep>
133
 
  </timestepping>
134
 
  <physical_parameters>
135
 
    <gravity>
136
 
      <magnitude>
137
 
        <real_value rank="0">10.0</real_value>
138
 
      </magnitude>
139
 
      <vector_field name="GravityDirection" rank="1">
140
 
        <prescribed>
141
 
          <mesh name="CoordinateMesh"/>
142
 
          <value name="WholeMesh">
143
 
            <constant>
144
 
              <real_value shape="3" dim1="dim" rank="1">0.0 0.0 -1.0</real_value>
145
 
            </constant>
146
 
          </value>
147
 
          <output>
148
 
            <exclude_from_vtu/>
149
 
          </output>
150
 
          <stat>
151
 
            <exclude_from_stat/>
152
 
          </stat>
153
 
          <detectors>
154
 
            <exclude_from_detectors/>
155
 
          </detectors>
156
 
        </prescribed>
157
 
      </vector_field>
158
 
    </gravity>
159
 
  </physical_parameters>
160
 
  <material_phase name="Water">
161
 
    <equation_of_state>
162
 
      <fluids>
163
 
        <linear>
164
 
          <reference_density>
165
 
            <real_value rank="0">1.0e+03</real_value>
166
 
          </reference_density>
167
 
        </linear>
168
 
      </fluids>
169
 
    </equation_of_state>
170
 
    <scalar_field name="Pressure" rank="0">
171
 
      <prognostic>
172
 
        <mesh name="PressureMesh"/>
173
 
        <spatial_discretisation>
174
 
          <continuous_galerkin>
175
 
            <test_continuity_with_cv_dual/>
176
 
          </continuous_galerkin>
177
 
        </spatial_discretisation>
178
 
        <scheme>
179
 
          <poisson_pressure_solution>
180
 
            <string_value lines="1">never</string_value>
181
 
          </poisson_pressure_solution>
182
 
          <use_projection_method/>
183
 
        </scheme>
184
 
        <solver>
185
 
          <iterative_method name="gmres">
186
 
            <restart>
187
 
              <integer_value rank="0">30</integer_value>
188
 
            </restart>
189
 
          </iterative_method>
190
 
          <preconditioner name="eisenstat"/>
191
 
          <relative_error>
192
 
            <real_value rank="0">1.0E-10</real_value>
193
 
          </relative_error>
194
 
          <max_iterations>
195
 
            <integer_value rank="0">5000</integer_value>
196
 
          </max_iterations>
197
 
          <never_ignore_solver_failures/>
198
 
          <diagnostics>
199
 
            <print_norms/>
200
 
            <monitors/>
201
 
          </diagnostics>
202
 
        </solver>
203
 
        <output/>
204
 
        <stat/>
205
 
        <convergence>
206
 
          <include_in_convergence/>
207
 
        </convergence>
208
 
        <detectors>
209
 
          <exclude_from_detectors/>
210
 
        </detectors>
211
 
        <steady_state>
212
 
          <include_in_steady_state/>
213
 
        </steady_state>
214
 
        <consistent_interpolation/>
215
 
      </prognostic>
216
 
    </scalar_field>
217
 
    <scalar_field name="Density" rank="0">
218
 
      <diagnostic>
219
 
        <algorithm name="Internal" material_phase_support="multiple"/>
220
 
        <mesh name="DensityMesh"/>
221
 
        <output/>
222
 
        <stat/>
223
 
        <convergence>
224
 
          <include_in_convergence/>
225
 
        </convergence>
226
 
        <detectors>
227
 
          <include_in_detectors/>
228
 
        </detectors>
229
 
        <steady_state>
230
 
          <include_in_steady_state/>
231
 
        </steady_state>
232
 
      </diagnostic>
233
 
    </scalar_field>
234
 
    <vector_field name="Velocity" rank="1">
235
 
      <prognostic>
236
 
        <mesh name="VelocityMesh"/>
237
 
        <equation name="LinearMomentum"/>
238
 
        <spatial_discretisation>
239
 
          <discontinuous_galerkin>
240
 
            <viscosity_scheme>
241
 
              <compact_discontinuous_galerkin/>
242
 
            </viscosity_scheme>
243
 
            <advection_scheme>
244
 
              <upwind/>
245
 
              <integrate_advection_by_parts>
246
 
                <once/>
247
 
              </integrate_advection_by_parts>
248
 
              <integrate_conservation_term_by_parts/>
249
 
            </advection_scheme>
250
 
          </discontinuous_galerkin>
251
 
          <conservative_advection>
252
 
            <real_value rank="0">0.0</real_value>
253
 
          </conservative_advection>
254
 
        </spatial_discretisation>
255
 
        <temporal_discretisation>
256
 
          <theta>
257
 
            <real_value rank="0">0.5</real_value>
258
 
          </theta>
259
 
          <relaxation>
260
 
            <real_value rank="0">0.5</real_value>
261
 
          </relaxation>
262
 
        </temporal_discretisation>
263
 
        <solver>
264
 
          <iterative_method name="gmres">
265
 
            <restart>
266
 
              <integer_value rank="0">30</integer_value>
267
 
            </restart>
268
 
          </iterative_method>
269
 
          <preconditioner name="eisenstat"/>
270
 
          <relative_error>
271
 
            <real_value rank="0">1.0E-10</real_value>
272
 
          </relative_error>
273
 
          <max_iterations>
274
 
            <integer_value rank="0">5000</integer_value>
275
 
          </max_iterations>
276
 
          <never_ignore_solver_failures/>
277
 
          <diagnostics>
278
 
            <print_norms/>
279
 
            <monitors/>
280
 
          </diagnostics>
281
 
        </solver>
282
 
        <initial_condition name="WholeMesh">
283
 
          <constant>
284
 
            <real_value shape="3" dim1="dim" rank="1">0.0 0.0 0.0</real_value>
285
 
          </constant>
286
 
        </initial_condition>
287
 
        <boundary_conditions name="no_flow_out_base">
288
 
          <surface_ids>
289
 
            <integer_value shape="1" rank="1">1</integer_value>
290
 
          </surface_ids>
291
 
          <type name="no_normal_flow"/>
292
 
        </boundary_conditions>
293
 
        <boundary_conditions name="no_flow_out_sides">
294
 
          <surface_ids>
295
 
            <integer_value shape="4" rank="1">3 4 5 6</integer_value>
296
 
          </surface_ids>
297
 
          <type name="no_normal_flow"/>
298
 
        </boundary_conditions>
299
 
        <tensor_field name="Viscosity" rank="2">
300
 
          <diagnostic>
301
 
            <mesh name="VelocityMesh"/>
302
 
            <algorithm name="bulk_viscosity" material_phase_support="multiple "/>
303
 
            <output/>
304
 
            <stat>
305
 
              <include_in_stat/>
306
 
            </stat>
307
 
          </diagnostic>
308
 
        </tensor_field>
309
 
        <output/>
310
 
        <stat>
311
 
          <include_in_stat/>
312
 
          <previous_time_step>
313
 
            <exclude_from_stat/>
314
 
          </previous_time_step>
315
 
          <nonlinear_field>
316
 
            <exclude_from_stat/>
317
 
          </nonlinear_field>
318
 
        </stat>
319
 
        <convergence>
320
 
          <include_in_convergence/>
321
 
        </convergence>
322
 
        <detectors>
323
 
          <include_in_detectors/>
324
 
        </detectors>
325
 
        <steady_state>
326
 
          <include_in_steady_state/>
327
 
        </steady_state>
328
 
        <consistent_interpolation/>
329
 
      </prognostic>
330
 
    </vector_field>
331
 
    <scalar_field name="MaterialVolumeFraction" rank="0">
332
 
      <prognostic>
333
 
        <mesh name="PressureMesh"/>
334
 
        <equation name="AdvectionDiffusion"/>
335
 
        <spatial_discretisation>
336
 
          <coupled_cv>
337
 
            <face_value name="HyperC">
338
 
              <project_upwind_value_from_point>
339
 
                <reflect_off_domain_boundaries/>
340
 
                <bound_projected_value_locally/>
341
 
                <store_upwind_elements/>
342
 
              </project_upwind_value_from_point>
343
 
              <courant_number name="ControlVolumeCFLNumber"/>
344
 
            </face_value>
345
 
            <parent_sum>
346
 
              <target_maximum>
347
 
                <real_value rank="0">1.0</real_value>
348
 
              </target_maximum>
349
 
              <target_minimum>
350
 
                <real_value rank="0">0.0</real_value>
351
 
              </target_minimum>
352
 
            </parent_sum>
353
 
          </coupled_cv>
354
 
          <conservative_advection>
355
 
            <real_value rank="0">0.0</real_value>
356
 
          </conservative_advection>
357
 
        </spatial_discretisation>
358
 
        <temporal_discretisation>
359
 
          <theta>
360
 
            <real_value rank="0">0.0</real_value>
361
 
          </theta>
362
 
          <control_volumes>
363
 
            <maximum_courant_number_per_subcycle>
364
 
              <real_value rank="0">0.25</real_value>
365
 
              <courant_number name="ControlVolumeCFLNumber"/>
366
 
            </maximum_courant_number_per_subcycle>
367
 
            <pivot_theta>
368
 
              <real_value rank="0">0.0</real_value>
369
 
            </pivot_theta>
370
 
          </control_volumes>
371
 
        </temporal_discretisation>
372
 
        <explicit/>
373
 
        <initial_condition name="WholeMesh">
374
 
          <python>
375
 
            <string_value lines="20" type="code" language="python">def val(X,t):
376
 
   if (X[2]&lt;=2.0):
377
 
      return 1.0
378
 
   else:
379
 
      return 0.0</string_value>
380
 
          </python>
381
 
        </initial_condition>
382
 
        <boundary_conditions name="incoming_vfrac">
383
 
          <surface_ids>
384
 
            <integer_value shape="1" rank="1">2</integer_value>
385
 
          </surface_ids>
386
 
          <type name="dirichlet">
387
 
            <apply_weakly/>
388
 
            <constant>
389
 
              <real_value rank="0">0.0</real_value>
390
 
            </constant>
391
 
          </type>
392
 
        </boundary_conditions>
393
 
        <output/>
394
 
        <stat>
395
 
          <include_cv_stats/>
396
 
        </stat>
397
 
        <convergence>
398
 
          <include_in_convergence/>
399
 
        </convergence>
400
 
        <detectors>
401
 
          <include_in_detectors/>
402
 
        </detectors>
403
 
        <steady_state>
404
 
          <include_in_steady_state/>
405
 
        </steady_state>
406
 
        <consistent_interpolation/>
407
 
        <priority>
408
 
          <integer_value rank="0">2</integer_value>
409
 
        </priority>
410
 
      </prognostic>
411
 
    </scalar_field>
412
 
    <scalar_field name="ControlVolumeCFLNumber" rank="0">
413
 
      <diagnostic>
414
 
        <algorithm name="Internal" material_phase_support="multiple"/>
415
 
        <mesh name="PressureMesh"/>
416
 
        <output/>
417
 
        <stat/>
418
 
        <convergence>
419
 
          <include_in_convergence/>
420
 
        </convergence>
421
 
        <detectors>
422
 
          <include_in_detectors/>
423
 
        </detectors>
424
 
        <steady_state>
425
 
          <include_in_steady_state/>
426
 
        </steady_state>
427
 
      </diagnostic>
428
 
    </scalar_field>
429
 
    <scalar_field name="SumMaterialVolumeFractions" rank="0">
430
 
      <diagnostic>
431
 
        <algorithm name="Internal" material_phase_support="multiple"/>
432
 
        <mesh name="PressureMesh"/>
433
 
        <output/>
434
 
        <stat/>
435
 
        <convergence>
436
 
          <include_in_convergence/>
437
 
        </convergence>
438
 
        <detectors>
439
 
          <include_in_detectors/>
440
 
        </detectors>
441
 
        <steady_state>
442
 
          <include_in_steady_state/>
443
 
        </steady_state>
444
 
      </diagnostic>
445
 
    </scalar_field>
446
 
    <scalar_field name="ControlVolumeDivergence" rank="0">
447
 
      <diagnostic field_name="Velocity">
448
 
        <algorithm name="Internal" material_phase_support="multiple"/>
449
 
        <mesh name="PressureMesh"/>
450
 
        <output/>
451
 
        <stat/>
452
 
        <convergence>
453
 
          <include_in_convergence/>
454
 
        </convergence>
455
 
        <detectors>
456
 
          <include_in_detectors/>
457
 
        </detectors>
458
 
        <steady_state>
459
 
          <include_in_steady_state/>
460
 
        </steady_state>
461
 
      </diagnostic>
462
 
    </scalar_field>
463
 
    <scalar_field name="DG_CourantNumber" rank="0">
464
 
      <diagnostic>
465
 
        <algorithm name="Internal" material_phase_support="multiple"/>
466
 
        <mesh name="VelocityMesh"/>
467
 
        <output/>
468
 
        <stat/>
469
 
        <convergence>
470
 
          <include_in_convergence/>
471
 
        </convergence>
472
 
        <detectors>
473
 
          <include_in_detectors/>
474
 
        </detectors>
475
 
        <steady_state>
476
 
          <include_in_steady_state/>
477
 
        </steady_state>
478
 
      </diagnostic>
479
 
    </scalar_field>
480
 
    <scalar_field name="GridReynoldsNumber" rank="0">
481
 
      <diagnostic>
482
 
        <algorithm name="Internal" material_phase_support="multiple"/>
483
 
        <mesh name="VelocityMesh"/>
484
 
        <include_density_field/>
485
 
        <output/>
486
 
        <stat/>
487
 
        <convergence>
488
 
          <include_in_convergence/>
489
 
        </convergence>
490
 
        <detectors>
491
 
          <include_in_detectors/>
492
 
        </detectors>
493
 
        <steady_state>
494
 
          <include_in_steady_state/>
495
 
        </steady_state>
496
 
      </diagnostic>
497
 
    </scalar_field>
498
 
    <tensor_field name="MaterialViscosity" rank="2">
499
 
      <prescribed>
500
 
        <mesh name="VelocityMesh"/>
501
 
        <value name="WholeMesh">
502
 
          <isotropic>
503
 
            <constant>
504
 
              <real_value rank="0">1.0</real_value>
505
 
              <comment>should be 9.0e-04</comment>
506
 
            </constant>
507
 
          </isotropic>
508
 
        </value>
509
 
        <output/>
510
 
      </prescribed>
511
 
    </tensor_field>
512
 
  </material_phase>
513
 
  <material_phase name="Droplet">
514
 
    <equation_of_state>
515
 
      <fluids>
516
 
        <linear>
517
 
          <reference_density>
518
 
            <real_value rank="0">1.0e+04</real_value>
519
 
          </reference_density>
520
 
        </linear>
521
 
      </fluids>
522
 
    </equation_of_state>
523
 
    <scalar_field name="Pressure" rank="0">
524
 
      <aliased material_phase_name="Water" field_name="Pressure"/>
525
 
    </scalar_field>
526
 
    <scalar_field name="Density" rank="0">
527
 
      <aliased material_phase_name="Water" field_name="Density"/>
528
 
    </scalar_field>
529
 
    <vector_field name="Velocity" rank="1">
530
 
      <aliased material_phase_name="Water" field_name="Velocity"/>
531
 
    </vector_field>
532
 
    <scalar_field name="MaterialVolumeFraction" rank="0">
533
 
      <prognostic>
534
 
        <mesh name="PressureMesh"/>
535
 
        <equation name="AdvectionDiffusion"/>
536
 
        <spatial_discretisation>
537
 
          <coupled_cv>
538
 
            <face_value name="HyperC">
539
 
              <project_upwind_value_from_point>
540
 
                <reflect_off_domain_boundaries/>
541
 
                <bound_projected_value_locally/>
542
 
                <store_upwind_elements/>
543
 
              </project_upwind_value_from_point>
544
 
              <courant_number name="ControlVolumeCFLNumber"/>
545
 
            </face_value>
546
 
            <parent_sum>
547
 
              <target_maximum>
548
 
                <real_value rank="0">1.0</real_value>
549
 
              </target_maximum>
550
 
              <target_minimum>
551
 
                <real_value rank="0">0.0</real_value>
552
 
              </target_minimum>
553
 
            </parent_sum>
554
 
          </coupled_cv>
555
 
          <conservative_advection>
556
 
            <real_value rank="0">0.0</real_value>
557
 
          </conservative_advection>
558
 
        </spatial_discretisation>
559
 
        <temporal_discretisation>
560
 
          <theta>
561
 
            <real_value rank="0">0.0</real_value>
562
 
          </theta>
563
 
          <control_volumes>
564
 
            <maximum_courant_number_per_subcycle>
565
 
              <real_value rank="0">0.25</real_value>
566
 
              <courant_number name="ControlVolumeCFLNumber"/>
567
 
            </maximum_courant_number_per_subcycle>
568
 
            <pivot_theta>
569
 
              <real_value rank="0">0.0</real_value>
570
 
            </pivot_theta>
571
 
          </control_volumes>
572
 
        </temporal_discretisation>
573
 
        <explicit/>
574
 
        <initial_condition name="WholeMesh">
575
 
          <python>
576
 
            <string_value lines="20" type="code" language="python">def val(X,t):
577
 
   from math import sqrt
578
 
   dx1 = X[0]-0.0
579
 
   dx2 = X[1]-0.0
580
 
   dx3 = X[2]-2.7
581
 
   r=sqrt(dx1**2+dx2**2+dx3**2)
582
 
   if (r&lt;0.5):
583
 
      return 1.0
584
 
   else:
585
 
      return 0.0</string_value>
586
 
          </python>
587
 
        </initial_condition>
588
 
        <boundary_conditions name="incoming_vfrac">
589
 
          <surface_ids>
590
 
            <integer_value shape="1" rank="1">2</integer_value>
591
 
          </surface_ids>
592
 
          <type name="dirichlet">
593
 
            <apply_weakly/>
594
 
            <constant>
595
 
              <real_value rank="0">0.0</real_value>
596
 
            </constant>
597
 
          </type>
598
 
        </boundary_conditions>
599
 
        <output/>
600
 
        <stat>
601
 
          <include_cv_stats/>
602
 
        </stat>
603
 
        <convergence>
604
 
          <include_in_convergence/>
605
 
        </convergence>
606
 
        <detectors>
607
 
          <include_in_detectors/>
608
 
        </detectors>
609
 
        <steady_state>
610
 
          <include_in_steady_state/>
611
 
        </steady_state>
612
 
        <consistent_interpolation/>
613
 
        <priority>
614
 
          <integer_value rank="0">3</integer_value>
615
 
        </priority>
616
 
      </prognostic>
617
 
    </scalar_field>
618
 
    <tensor_field name="MaterialViscosity" rank="2">
619
 
      <prescribed>
620
 
        <mesh name="VelocityMesh"/>
621
 
        <value name="WholeMesh">
622
 
          <isotropic>
623
 
            <constant>
624
 
              <real_value rank="0">1.0</real_value>
625
 
              <comment>should be 5.0e-03</comment>
626
 
            </constant>
627
 
          </isotropic>
628
 
        </value>
629
 
        <output/>
630
 
      </prescribed>
631
 
    </tensor_field>
632
 
  </material_phase>
633
 
  <material_phase name="Air">
634
 
    <equation_of_state>
635
 
      <fluids>
636
 
        <linear>
637
 
          <reference_density>
638
 
            <real_value rank="0">1.0</real_value>
639
 
          </reference_density>
640
 
        </linear>
641
 
      </fluids>
642
 
    </equation_of_state>
643
 
    <scalar_field name="Pressure" rank="0">
644
 
      <aliased material_phase_name="Water" field_name="Pressure"/>
645
 
    </scalar_field>
646
 
    <scalar_field name="Density" rank="0">
647
 
      <aliased material_phase_name="Water" field_name="Density"/>
648
 
    </scalar_field>
649
 
    <vector_field name="Velocity" rank="1">
650
 
      <aliased material_phase_name="Water" field_name="Velocity"/>
651
 
    </vector_field>
652
 
    <scalar_field name="MaterialVolumeFraction" rank="0">
653
 
      <diagnostic>
654
 
        <algorithm name="Internal" material_phase_support="multiple"/>
655
 
        <mesh name="PressureMesh"/>
656
 
        <output/>
657
 
        <stat>
658
 
          <include_cv_stats/>
659
 
        </stat>
660
 
        <convergence>
661
 
          <include_in_convergence/>
662
 
        </convergence>
663
 
        <detectors>
664
 
          <include_in_detectors/>
665
 
        </detectors>
666
 
        <steady_state>
667
 
          <include_in_steady_state/>
668
 
        </steady_state>
669
 
      </diagnostic>
670
 
    </scalar_field>
671
 
    <tensor_field name="MaterialViscosity" rank="2">
672
 
      <prescribed>
673
 
        <mesh name="VelocityMesh"/>
674
 
        <value name="WholeMesh">
675
 
          <isotropic>
676
 
            <constant>
677
 
              <real_value rank="0">1.0</real_value>
678
 
              <comment>should be 2.0e-05</comment>
679
 
            </constant>
680
 
          </isotropic>
681
 
        </value>
682
 
        <output/>
683
 
      </prescribed>
684
 
    </tensor_field>
685
 
  </material_phase>
686
 
</fluidity_options>