← Back to branch summary

~f-milthaler/fluidity/fsi-model-stationary-solid-with-velocity

« back to all changes in this revision

Viewing changes to tests/mphase_inlet_velocity_bc_compressible/mphase_inlet_velocity_bc_compressible.flml

  • Committer: f.milthaler10 at uk
  • Date: 2013-11-06 13:43:56 UTC
  • mfrom: (3463.184.85 fluidity)
  • Revision ID: f.milthaler10@imperial.ac.ic.uk.-20131106134356-v3lw1dheesckywj0
merge from trunk

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/mphase_inlet_velocity_bc_compressible/mphase_inlet_velocity_bc_compressible.flml
 
1
<?xml version='1.0' encoding='utf-8'?>
 
2
<fluidity_options>
 
3
  <simulation_name>
 
4
    <string_value lines="1">mphase_inlet_velocity_bc_compressible</string_value>
 
5
  </simulation_name>
 
6
  <problem_type>
 
7
    <string_value lines="1">fluids</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">0</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">1</integer_value>
 
43
          </polynomial_degree>
 
44
        </mesh_shape>
 
45
        <stat>
 
46
          <exclude_from_stat/>
 
47
        </stat>
 
48
      </from_mesh>
 
49
    </mesh>
 
50
    <mesh name="DensityMesh">
 
51
      <from_mesh>
 
52
        <mesh name="CoordinateMesh"/>
 
53
        <mesh_shape>
 
54
          <polynomial_degree>
 
55
            <integer_value rank="0">1</integer_value>
 
56
          </polynomial_degree>
 
57
        </mesh_shape>
 
58
        <stat>
 
59
          <exclude_from_stat/>
 
60
        </stat>
 
61
      </from_mesh>
 
62
    </mesh>
 
63
    <quadrature>
 
64
      <degree>
 
65
        <integer_value rank="0">3</integer_value>
 
66
      </degree>
 
67
    </quadrature>
 
68
  </geometry>
 
69
  <io>
 
70
    <dump_format>
 
71
      <string_value>vtk</string_value>
 
72
    </dump_format>
 
73
    <dump_period>
 
74
      <constant>
 
75
        <real_value rank="0">0</real_value>
 
76
      </constant>
 
77
    </dump_period>
 
78
    <output_mesh name="PressureMesh"/>
 
79
    <stat>
 
80
      <output_at_start/>
 
81
    </stat>
 
82
  </io>
 
83
  <timestepping>
 
84
    <current_time>
 
85
      <real_value rank="0">0</real_value>
 
86
    </current_time>
 
87
    <timestep>
 
88
      <real_value rank="0">1.0</real_value>
 
89
    </timestep>
 
90
    <finish_time>
 
91
      <real_value rank="0">100.0</real_value>
 
92
    </finish_time>
 
93
    <nonlinear_iterations>
 
94
      <integer_value rank="0">50</integer_value>
 
95
      <tolerance>
 
96
        <real_value rank="0">1.0e-4</real_value>
 
97
        <infinity_norm/>
 
98
      </tolerance>
 
99
    </nonlinear_iterations>
 
100
  </timestepping>
 
101
  <physical_parameters>
 
102
    <gravity>
 
103
      <magnitude>
 
104
        <real_value rank="0">9.8</real_value>
 
105
      </magnitude>
 
106
      <vector_field name="GravityDirection" rank="1">
 
107
        <prescribed>
 
108
          <mesh name="CoordinateMesh"/>
 
109
          <value name="WholeMesh">
 
110
            <constant>
 
111
              <real_value shape="1" dim1="dim" rank="1">-1.0</real_value>
 
112
            </constant>
 
113
          </value>
 
114
          <output/>
 
115
          <stat>
 
116
            <include_in_stat/>
 
117
          </stat>
 
118
          <detectors>
 
119
            <exclude_from_detectors/>
 
120
          </detectors>
 
121
        </prescribed>
 
122
      </vector_field>
 
123
    </gravity>
 
124
  </physical_parameters>
 
125
  <material_phase name="Gas">
 
126
    <equation_of_state>
 
127
      <compressible>
 
128
        <stiffened_gas>
 
129
          <ratio_specific_heats>
 
130
            <real_value rank="0">1.33</real_value>
 
131
          </ratio_specific_heats>
 
132
        </stiffened_gas>
 
133
        <subtract_out_reference_profile/>
 
134
      </compressible>
 
135
    </equation_of_state>
 
136
    <scalar_field name="Pressure" rank="0">
 
137
      <prognostic>
 
138
        <mesh name="PressureMesh"/>
 
139
        <spatial_discretisation>
 
140
          <continuous_galerkin>
 
141
            <remove_stabilisation_term/>
 
142
            <integrate_continuity_by_parts/>
 
143
          </continuous_galerkin>
 
144
        </spatial_discretisation>
 
145
        <scheme>
 
146
          <poisson_pressure_solution>
 
147
            <string_value lines="1">never</string_value>
 
148
          </poisson_pressure_solution>
 
149
          <use_projection_method/>
 
150
        </scheme>
 
151
        <solver>
 
152
          <iterative_method name="gmres">
 
153
            <restart>
 
154
              <integer_value rank="0">30</integer_value>
 
155
            </restart>
 
156
          </iterative_method>
 
157
          <preconditioner name="sor"/>
 
158
          <relative_error>
 
159
            <real_value rank="0">1.0e-7</real_value>
 
160
          </relative_error>
 
161
          <max_iterations>
 
162
            <integer_value rank="0">1000</integer_value>
 
163
          </max_iterations>
 
164
          <never_ignore_solver_failures/>
 
165
          <diagnostics>
 
166
            <monitors/>
 
167
          </diagnostics>
 
168
        </solver>
 
169
        <initial_condition name="WholeMesh">
 
170
          <python>
 
171
            <string_value lines="20" type="code" language="python">def val(X,t):
 
172
  from math import exp
 
173
  T0 = 288.15 # K
 
174
  R = 287.058 # Specific gas constant in J/Kg
 
175
  g = 9.8 # m/s**2
 
176
  h = X[0] # m
 
177
  p0 = 101325
 
178
  e = T0*1406.0
 
179
  gamma = 1.33
 
180
  p = p0*exp(-(g*h)/((gamma-1.0)*e))
 
181
  return p</string_value>
 
182
          </python>
 
183
        </initial_condition>
 
184
        <boundary_conditions name="Right">
 
185
          <surface_ids>
 
186
            <integer_value shape="1" rank="1">2</integer_value>
 
187
          </surface_ids>
 
188
          <type name="dirichlet">
 
189
            <python>
 
190
              <string_value lines="20" type="code" language="python">def val(X,t):
 
191
  from math import exp
 
192
  T0 = 288.15 # K
 
193
  R = 287.058 # Specific gas constant in J/Kg
 
194
  g = 9.8 # m/s**2
 
195
  h = X[0] # m
 
196
  p0 = 101325
 
197
  e = T0*1406.0
 
198
  gamma = 1.33
 
199
  p = p0*exp(-(g*h)/((gamma-1.0)*e))
 
200
  return p</string_value>
 
201
            </python>
 
202
          </type>
 
203
        </boundary_conditions>
 
204
        <boundary_conditions name="Left">
 
205
          <surface_ids>
 
206
            <integer_value shape="1" rank="1">1</integer_value>
 
207
          </surface_ids>
 
208
          <type name="dirichlet">
 
209
            <constant>
 
210
              <real_value rank="0">101325</real_value>
 
211
            </constant>
 
212
          </type>
 
213
        </boundary_conditions>
 
214
        <output/>
 
215
        <stat/>
 
216
        <convergence>
 
217
          <include_in_convergence/>
 
218
        </convergence>
 
219
        <detectors>
 
220
          <exclude_from_detectors/>
 
221
        </detectors>
 
222
        <steady_state>
 
223
          <include_in_steady_state/>
 
224
        </steady_state>
 
225
        <no_interpolation/>
 
226
      </prognostic>
 
227
    </scalar_field>
 
228
    <scalar_field name="Density" rank="0">
 
229
      <prognostic>
 
230
        <mesh name="DensityMesh"/>
 
231
        <spatial_discretisation>
 
232
          <control_volumes>
 
233
            <face_value name="FirstOrderUpwind"/>
 
234
          </control_volumes>
 
235
          <conservative_advection>
 
236
            <real_value rank="0">1.0</real_value>
 
237
          </conservative_advection>
 
238
        </spatial_discretisation>
 
239
        <temporal_discretisation>
 
240
          <theta>
 
241
            <real_value rank="0">1.0</real_value>
 
242
          </theta>
 
243
        </temporal_discretisation>
 
244
        <initial_condition name="WholeMesh">
 
245
          <python>
 
246
            <string_value lines="20" type="code" language="python">def val(X,t):
 
247
  from math import exp
 
248
  T0 = 288.15 # K
 
249
  R = 287.058 # Specific gas constant in J/Kg
 
250
  g = 9.8 # m/s**2
 
251
  h = X[0] # m
 
252
  p0 = 101325
 
253
  e = T0*1406.0
 
254
  gamma = 1.33
 
255
  p = p0*exp(-(g*h)/((gamma-1.0)*e))
 
256
  rho = p/((gamma-1.0)*e)
 
257
  return rho</string_value>
 
258
          </python>
 
259
        </initial_condition>
 
260
        <boundary_conditions name="Left">
 
261
          <surface_ids>
 
262
            <integer_value shape="1" rank="1">1</integer_value>
 
263
          </surface_ids>
 
264
          <type name="dirichlet">
 
265
            <constant>
 
266
              <real_value rank="0">0.181985215</real_value>
 
267
            </constant>
 
268
          </type>
 
269
        </boundary_conditions>
 
270
        <output/>
 
271
        <stat/>
 
272
        <convergence>
 
273
          <include_in_convergence/>
 
274
        </convergence>
 
275
        <detectors>
 
276
          <include_in_detectors/>
 
277
        </detectors>
 
278
        <steady_state>
 
279
          <include_in_steady_state/>
 
280
        </steady_state>
 
281
        <consistent_interpolation/>
 
282
      </prognostic>
 
283
    </scalar_field>
 
284
    <vector_field name="Velocity" rank="1">
 
285
      <prognostic>
 
286
        <mesh name="VelocityMesh"/>
 
287
        <equation name="LinearMomentum"/>
 
288
        <spatial_discretisation>
 
289
          <discontinuous_galerkin>
 
290
            <viscosity_scheme>
 
291
              <compact_discontinuous_galerkin/>
 
292
              <tensor_form/>
 
293
            </viscosity_scheme>
 
294
            <advection_scheme>
 
295
              <upwind/>
 
296
              <integrate_advection_by_parts>
 
297
                <twice/>
 
298
              </integrate_advection_by_parts>
 
299
            </advection_scheme>
 
300
          </discontinuous_galerkin>
 
301
          <conservative_advection>
 
302
            <real_value rank="0">0</real_value>
 
303
          </conservative_advection>
 
304
        </spatial_discretisation>
 
305
        <temporal_discretisation>
 
306
          <theta>
 
307
            <real_value rank="0">1.0</real_value>
 
308
          </theta>
 
309
          <relaxation>
 
310
            <real_value rank="0">0.5</real_value>
 
311
          </relaxation>
 
312
        </temporal_discretisation>
 
313
        <solver>
 
314
          <iterative_method name="gmres">
 
315
            <restart>
 
316
              <integer_value rank="0">30</integer_value>
 
317
            </restart>
 
318
          </iterative_method>
 
319
          <preconditioner name="sor"/>
 
320
          <relative_error>
 
321
            <real_value rank="0">1.0e-7</real_value>
 
322
          </relative_error>
 
323
          <max_iterations>
 
324
            <integer_value rank="0">1000</integer_value>
 
325
          </max_iterations>
 
326
          <never_ignore_solver_failures/>
 
327
          <diagnostics>
 
328
            <monitors/>
 
329
          </diagnostics>
 
330
        </solver>
 
331
        <initial_condition name="WholeMesh">
 
332
          <constant>
 
333
            <real_value shape="1" dim1="dim" rank="1">0</real_value>
 
334
          </constant>
 
335
        </initial_condition>
 
336
        <boundary_conditions name="Left">
 
337
          <surface_ids>
 
338
            <integer_value shape="1" rank="1">1</integer_value>
 
339
          </surface_ids>
 
340
          <type name="dirichlet">
 
341
            <align_bc_with_cartesian>
 
342
              <x_component>
 
343
                <constant>
 
344
                  <real_value rank="0">300</real_value>
 
345
                </constant>
 
346
              </x_component>
 
347
            </align_bc_with_cartesian>
 
348
          </type>
 
349
        </boundary_conditions>
 
350
        <tensor_field name="Viscosity" rank="2">
 
351
          <prescribed>
 
352
            <value name="WholeMesh">
 
353
              <anisotropic_asymmetric>
 
354
                <constant>
 
355
                  <real_value symmetric="false" dim2="dim" shape="1 1" dim1="dim" rank="2">1.85e-5</real_value>
 
356
                </constant>
 
357
              </anisotropic_asymmetric>
 
358
            </value>
 
359
            <output/>
 
360
          </prescribed>
 
361
        </tensor_field>
 
362
        <output/>
 
363
        <stat>
 
364
          <include_in_stat/>
 
365
          <previous_time_step>
 
366
            <exclude_from_stat/>
 
367
          </previous_time_step>
 
368
          <nonlinear_field>
 
369
            <exclude_from_stat/>
 
370
          </nonlinear_field>
 
371
        </stat>
 
372
        <convergence>
 
373
          <include_in_convergence/>
 
374
        </convergence>
 
375
        <detectors>
 
376
          <include_in_detectors/>
 
377
        </detectors>
 
378
        <steady_state>
 
379
          <include_in_steady_state/>
 
380
        </steady_state>
 
381
        <consistent_interpolation/>
 
382
      </prognostic>
 
383
    </vector_field>
 
384
    <scalar_field name="CFLNumber" rank="0">
 
385
      <diagnostic>
 
386
        <algorithm name="Internal" material_phase_support="multiple"/>
 
387
        <mesh name="VelocityMesh"/>
 
388
        <output/>
 
389
        <stat/>
 
390
        <convergence>
 
391
          <include_in_convergence/>
 
392
        </convergence>
 
393
        <detectors>
 
394
          <include_in_detectors/>
 
395
        </detectors>
 
396
        <steady_state>
 
397
          <include_in_steady_state/>
 
398
        </steady_state>
 
399
      </diagnostic>
 
400
    </scalar_field>
 
401
    <scalar_field name="InternalEnergy" rank="0">
 
402
      <prognostic>
 
403
        <mesh name="PressureMesh"/>
 
404
        <equation name="InternalEnergy">
 
405
          <density name="Density">
 
406
            <discretisation_options>
 
407
              <spatial_discretisation>
 
408
                <control_volumes>
 
409
                  <face_value name="FirstOrderUpwind"/>
 
410
                </control_volumes>
 
411
              </spatial_discretisation>
 
412
              <temporal_discretisation>
 
413
                <theta>
 
414
                  <real_value rank="0">1.0</real_value>
 
415
                </theta>
 
416
                <control_volumes/>
 
417
              </temporal_discretisation>
 
418
            </discretisation_options>
 
419
          </density>
 
420
        </equation>
 
421
        <spatial_discretisation>
 
422
          <control_volumes>
 
423
            <face_value name="FirstOrderUpwind"/>
 
424
            <diffusion_scheme name="ElementGradient"/>
 
425
          </control_volumes>
 
426
          <conservative_advection>
 
427
            <real_value rank="0">0.0</real_value>
 
428
          </conservative_advection>
 
429
        </spatial_discretisation>
 
430
        <temporal_discretisation>
 
431
          <theta>
 
432
            <real_value rank="0">1.0</real_value>
 
433
          </theta>
 
434
        </temporal_discretisation>
 
435
        <solver>
 
436
          <iterative_method name="gmres">
 
437
            <restart>
 
438
              <integer_value rank="0">30</integer_value>
 
439
            </restart>
 
440
          </iterative_method>
 
441
          <preconditioner name="sor"/>
 
442
          <relative_error>
 
443
            <real_value rank="0">1.0e-7</real_value>
 
444
          </relative_error>
 
445
          <max_iterations>
 
446
            <integer_value rank="0">1000</integer_value>
 
447
          </max_iterations>
 
448
          <never_ignore_solver_failures/>
 
449
          <diagnostics>
 
450
            <monitors/>
 
451
          </diagnostics>
 
452
        </solver>
 
453
        <initial_condition name="WholeMesh">
 
454
          <constant>
 
455
            <real_value rank="0">405138.9</real_value>
 
456
          </constant>
 
457
        </initial_condition>
 
458
        <boundary_conditions name="Left">
 
459
          <surface_ids>
 
460
            <integer_value shape="1" rank="1">1</integer_value>
 
461
          </surface_ids>
 
462
          <type name="dirichlet">
 
463
            <constant>
 
464
              <real_value rank="0">1687200</real_value>
 
465
            </constant>
 
466
          </type>
 
467
        </boundary_conditions>
 
468
        <output/>
 
469
        <stat/>
 
470
        <convergence>
 
471
          <include_in_convergence/>
 
472
        </convergence>
 
473
        <detectors>
 
474
          <include_in_detectors/>
 
475
        </detectors>
 
476
        <steady_state>
 
477
          <include_in_steady_state/>
 
478
        </steady_state>
 
479
        <consistent_interpolation/>
 
480
      </prognostic>
 
481
    </scalar_field>
 
482
    <scalar_field name="HydrostaticReferencePressure" rank="0">
 
483
      <prescribed>
 
484
        <mesh name="PressureMesh"/>
 
485
        <value name="WholeMesh">
 
486
          <python>
 
487
            <string_value lines="20" type="code" language="python">def val(X,t):
 
488
  from math import exp
 
489
  T0 = 288.15 # K
 
490
  R = 287.058 # Specific gas constant in J/Kg
 
491
  g = 9.8 # m/s**2
 
492
  h = X[0] # m
 
493
  p0 = 101325
 
494
  e = T0*1406.0
 
495
  gamma = 1.33
 
496
  p = p0*exp(-(g*h)/((gamma-1.0)*e))
 
497
  return p</string_value>
 
498
          </python>
 
499
        </value>
 
500
        <output/>
 
501
        <stat/>
 
502
        <detectors>
 
503
          <exclude_from_detectors/>
 
504
        </detectors>
 
505
      </prescribed>
 
506
    </scalar_field>
 
507
    <scalar_field name="HydrostaticReferenceDensity" rank="0">
 
508
      <prescribed>
 
509
        <mesh name="PressureMesh"/>
 
510
        <value name="WholeMesh">
 
511
          <python>
 
512
            <string_value lines="20" type="code" language="python">def val(X,t):
 
513
  from math import exp
 
514
  T0 = 288.15 # K
 
515
  R = 287.058 # Specific gas constant in J/Kg
 
516
  g = 9.8 # m/s**2
 
517
  h = X[0] # m
 
518
  p0 = 101325
 
519
  e = T0*1406.0
 
520
  gamma = 1.33
 
521
  p = p0*exp(-(g*h)/((gamma-1.0)*e))
 
522
  rho = p/((gamma-1.0)*e)
 
523
  return rho</string_value>
 
524
          </python>
 
525
        </value>
 
526
        <output/>
 
527
        <stat/>
 
528
        <detectors>
 
529
          <exclude_from_detectors/>
 
530
        </detectors>
 
531
      </prescribed>
 
532
    </scalar_field>
 
533
    <scalar_field name="PhaseVolumeFraction" rank="0">
 
534
      <diagnostic>
 
535
        <mesh name="CoordinateMesh"/>
 
536
        <algorithm name="Internal" material_phase_support="multiple"/>
 
537
        <output/>
 
538
        <stat/>
 
539
        <detectors>
 
540
          <include_in_detectors/>
 
541
        </detectors>
 
542
      </diagnostic>
 
543
    </scalar_field>
 
544
    <multiphase_properties>
 
545
      <effective_conductivity>
 
546
        <real_value rank="0">0.5</real_value>
 
547
      </effective_conductivity>
 
548
      <specific_heat>
 
549
        <real_value rank="0">1406</real_value>
 
550
      </specific_heat>
 
551
    </multiphase_properties>
 
552
  </material_phase>
 
553
  <material_phase name="Particles">
 
554
    <equation_of_state>
 
555
      <fluids>
 
556
        <linear>
 
557
          <reference_density>
 
558
            <real_value rank="0">2400</real_value>
 
559
          </reference_density>
 
560
        </linear>
 
561
      </fluids>
 
562
    </equation_of_state>
 
563
    <scalar_field name="Pressure" rank="0">
 
564
      <aliased material_phase_name="Gas" field_name="Pressure"/>
 
565
    </scalar_field>
 
566
    <scalar_field name="Density" rank="0">
 
567
      <diagnostic>
 
568
        <algorithm name="Internal" material_phase_support="multiple"/>
 
569
        <mesh name="PressureMesh"/>
 
570
        <output/>
 
571
        <stat/>
 
572
        <convergence>
 
573
          <include_in_convergence/>
 
574
        </convergence>
 
575
        <detectors>
 
576
          <include_in_detectors/>
 
577
        </detectors>
 
578
        <steady_state>
 
579
          <include_in_steady_state/>
 
580
        </steady_state>
 
581
      </diagnostic>
 
582
    </scalar_field>
 
583
    <vector_field name="Velocity" rank="1">
 
584
      <prognostic>
 
585
        <mesh name="VelocityMesh"/>
 
586
        <equation name="LinearMomentum"/>
 
587
        <spatial_discretisation>
 
588
          <discontinuous_galerkin>
 
589
            <viscosity_scheme>
 
590
              <compact_discontinuous_galerkin/>
 
591
              <tensor_form/>
 
592
            </viscosity_scheme>
 
593
            <advection_scheme>
 
594
              <upwind/>
 
595
              <integrate_advection_by_parts>
 
596
                <twice/>
 
597
              </integrate_advection_by_parts>
 
598
            </advection_scheme>
 
599
          </discontinuous_galerkin>
 
600
          <conservative_advection>
 
601
            <real_value rank="0">0</real_value>
 
602
          </conservative_advection>
 
603
        </spatial_discretisation>
 
604
        <temporal_discretisation>
 
605
          <theta>
 
606
            <real_value rank="0">1.0</real_value>
 
607
          </theta>
 
608
          <relaxation>
 
609
            <real_value rank="0">0.5</real_value>
 
610
          </relaxation>
 
611
        </temporal_discretisation>
 
612
        <solver>
 
613
          <iterative_method name="gmres">
 
614
            <restart>
 
615
              <integer_value rank="0">30</integer_value>
 
616
            </restart>
 
617
          </iterative_method>
 
618
          <preconditioner name="sor"/>
 
619
          <relative_error>
 
620
            <real_value rank="0">1.0e-7</real_value>
 
621
          </relative_error>
 
622
          <max_iterations>
 
623
            <integer_value rank="0">1000</integer_value>
 
624
          </max_iterations>
 
625
          <never_ignore_solver_failures/>
 
626
          <diagnostics>
 
627
            <monitors/>
 
628
          </diagnostics>
 
629
        </solver>
 
630
        <initial_condition name="WholeMesh">
 
631
          <constant>
 
632
            <real_value shape="1" dim1="dim" rank="1">0</real_value>
 
633
          </constant>
 
634
        </initial_condition>
 
635
        <boundary_conditions name="Left">
 
636
          <surface_ids>
 
637
            <integer_value shape="1" rank="1">1</integer_value>
 
638
          </surface_ids>
 
639
          <type name="dirichlet">
 
640
            <align_bc_with_cartesian>
 
641
              <x_component>
 
642
                <constant>
 
643
                  <real_value rank="0">300</real_value>
 
644
                </constant>
 
645
              </x_component>
 
646
            </align_bc_with_cartesian>
 
647
          </type>
 
648
        </boundary_conditions>
 
649
        <tensor_field name="Viscosity" rank="2">
 
650
          <prescribed>
 
651
            <value name="WholeMesh">
 
652
              <isotropic>
 
653
                <constant>
 
654
                  <real_value rank="0">0.5</real_value>
 
655
                </constant>
 
656
              </isotropic>
 
657
            </value>
 
658
            <output/>
 
659
          </prescribed>
 
660
        </tensor_field>
 
661
        <output/>
 
662
        <stat>
 
663
          <include_in_stat/>
 
664
          <previous_time_step>
 
665
            <exclude_from_stat/>
 
666
          </previous_time_step>
 
667
          <nonlinear_field>
 
668
            <exclude_from_stat/>
 
669
          </nonlinear_field>
 
670
        </stat>
 
671
        <convergence>
 
672
          <include_in_convergence/>
 
673
        </convergence>
 
674
        <detectors>
 
675
          <include_in_detectors/>
 
676
        </detectors>
 
677
        <steady_state>
 
678
          <include_in_steady_state/>
 
679
        </steady_state>
 
680
        <consistent_interpolation/>
 
681
      </prognostic>
 
682
    </vector_field>
 
683
    <scalar_field name="CFLNumber" rank="0">
 
684
      <diagnostic>
 
685
        <algorithm name="Internal" material_phase_support="multiple"/>
 
686
        <mesh name="VelocityMesh"/>
 
687
        <output/>
 
688
        <stat/>
 
689
        <convergence>
 
690
          <include_in_convergence/>
 
691
        </convergence>
 
692
        <detectors>
 
693
          <include_in_detectors/>
 
694
        </detectors>
 
695
        <steady_state>
 
696
          <include_in_steady_state/>
 
697
        </steady_state>
 
698
      </diagnostic>
 
699
    </scalar_field>
 
700
    <scalar_field name="InternalEnergy" rank="0">
 
701
      <prognostic>
 
702
        <mesh name="PressureMesh"/>
 
703
        <equation name="InternalEnergy">
 
704
          <density name="Density">
 
705
            <discretisation_options>
 
706
              <spatial_discretisation>
 
707
                <control_volumes>
 
708
                  <face_value name="FirstOrderUpwind"/>
 
709
                </control_volumes>
 
710
              </spatial_discretisation>
 
711
              <temporal_discretisation>
 
712
                <theta>
 
713
                  <real_value rank="0">1.0</real_value>
 
714
                </theta>
 
715
                <control_volumes/>
 
716
              </temporal_discretisation>
 
717
            </discretisation_options>
 
718
          </density>
 
719
        </equation>
 
720
        <spatial_discretisation>
 
721
          <control_volumes>
 
722
            <face_value name="FirstOrderUpwind"/>
 
723
            <diffusion_scheme name="ElementGradient"/>
 
724
          </control_volumes>
 
725
          <conservative_advection>
 
726
            <real_value rank="0">0.0</real_value>
 
727
          </conservative_advection>
 
728
        </spatial_discretisation>
 
729
        <temporal_discretisation>
 
730
          <theta>
 
731
            <real_value rank="0">1.0</real_value>
 
732
          </theta>
 
733
        </temporal_discretisation>
 
734
        <solver>
 
735
          <iterative_method name="gmres">
 
736
            <restart>
 
737
              <integer_value rank="0">30</integer_value>
 
738
            </restart>
 
739
          </iterative_method>
 
740
          <preconditioner name="sor"/>
 
741
          <relative_error>
 
742
            <real_value rank="0">1.0e-7</real_value>
 
743
          </relative_error>
 
744
          <max_iterations>
 
745
            <integer_value rank="0">1000</integer_value>
 
746
          </max_iterations>
 
747
          <never_ignore_solver_failures/>
 
748
          <diagnostics>
 
749
            <monitors/>
 
750
          </diagnostics>
 
751
        </solver>
 
752
        <initial_condition name="WholeMesh">
 
753
          <constant>
 
754
            <real_value rank="0">274895.1</real_value>
 
755
          </constant>
 
756
        </initial_condition>
 
757
        <boundary_conditions name="Left">
 
758
          <surface_ids>
 
759
            <integer_value shape="1" rank="1">1</integer_value>
 
760
          </surface_ids>
 
761
          <type name="dirichlet">
 
762
            <constant>
 
763
              <real_value rank="0">1144800</real_value>
 
764
            </constant>
 
765
          </type>
 
766
        </boundary_conditions>
 
767
        <output/>
 
768
        <stat/>
 
769
        <convergence>
 
770
          <include_in_convergence/>
 
771
        </convergence>
 
772
        <detectors>
 
773
          <include_in_detectors/>
 
774
        </detectors>
 
775
        <steady_state>
 
776
          <include_in_steady_state/>
 
777
        </steady_state>
 
778
        <consistent_interpolation/>
 
779
      </prognostic>
 
780
    </scalar_field>
 
781
    <scalar_field name="PhaseVolumeFraction" rank="0">
 
782
      <prognostic>
 
783
        <mesh name="CoordinateMesh"/>
 
784
        <equation name="AdvectionDiffusion"/>
 
785
        <spatial_discretisation>
 
786
          <control_volumes>
 
787
            <face_value name="FirstOrderUpwind"/>
 
788
            <diffusion_scheme name="ElementGradient"/>
 
789
          </control_volumes>
 
790
          <conservative_advection>
 
791
            <real_value rank="0">0.0</real_value>
 
792
          </conservative_advection>
 
793
        </spatial_discretisation>
 
794
        <temporal_discretisation>
 
795
          <theta>
 
796
            <real_value rank="0">1.0</real_value>
 
797
          </theta>
 
798
        </temporal_discretisation>
 
799
        <solver>
 
800
          <iterative_method name="gmres">
 
801
            <restart>
 
802
              <integer_value rank="0">30</integer_value>
 
803
            </restart>
 
804
          </iterative_method>
 
805
          <preconditioner name="sor"/>
 
806
          <relative_error>
 
807
            <real_value rank="0">1.0e-7</real_value>
 
808
          </relative_error>
 
809
          <max_iterations>
 
810
            <integer_value rank="0">1000</integer_value>
 
811
          </max_iterations>
 
812
          <never_ignore_solver_failures/>
 
813
          <diagnostics>
 
814
            <monitors/>
 
815
          </diagnostics>
 
816
        </solver>
 
817
        <initial_condition name="WholeMesh">
 
818
          <constant>
 
819
            <real_value rank="0">1.0e-7</real_value>
 
820
          </constant>
 
821
        </initial_condition>
 
822
        <boundary_conditions name="Left">
 
823
          <surface_ids>
 
824
            <integer_value shape="1" rank="1">1</integer_value>
 
825
          </surface_ids>
 
826
          <type name="dirichlet">
 
827
            <constant>
 
828
              <real_value rank="0">1.0e-4</real_value>
 
829
            </constant>
 
830
          </type>
 
831
        </boundary_conditions>
 
832
        <output/>
 
833
        <stat/>
 
834
        <convergence>
 
835
          <include_in_convergence/>
 
836
        </convergence>
 
837
        <detectors>
 
838
          <include_in_detectors/>
 
839
        </detectors>
 
840
        <steady_state>
 
841
          <include_in_steady_state/>
 
842
        </steady_state>
 
843
        <consistent_interpolation/>
 
844
      </prognostic>
 
845
    </scalar_field>
 
846
    <multiphase_properties>
 
847
      <particle_diameter>
 
848
        <real_value rank="0">2e-4</real_value>
 
849
      </particle_diameter>
 
850
      <specific_heat>
 
851
        <real_value rank="0">954</real_value>
 
852
      </specific_heat>
 
853
    </multiphase_properties>
 
854
  </material_phase>
 
855
  <multiphase_interaction>
 
856
    <fluid_particle_drag>
 
857
      <drag_correlation name="wen_yu"/>
 
858
    </fluid_particle_drag>
 
859
    <heat_transfer>
 
860
      <heat_transfer_coefficient name="gunn"/>
 
861
    </heat_transfer>
 
862
  </multiphase_interaction>
 
863
</fluidity_options>