← Back to branch summary

~fluidity-core/fluidity/excise-fldecomp

« back to all changes in this revision

Viewing changes to tests/mphase_dusty_gas_shock_tube/mphase_dusty_gas_shock_tube.flml

  • Committer: Mark Filipiak
  • Date: 2012-08-13 11:42:30 UTC
  • mfrom: (4003.1.23 dev-trunk)
  • Revision ID: mjf@staffmail.ed.ac.uk-20120813114230-wzoyf2gi4p4oxeh4
Merge in of the latest trunk.  To try to cure non-flredecomp tests that are passing at EPCC but failing in buildbot.

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/mphase_dusty_gas_shock_tube/mphase_dusty_gas_shock_tube.flml
 
1
<?xml version='1.0' encoding='utf-8'?>
 
2
<fluidity_options>
 
3
  <simulation_name>
 
4
    <string_value lines="1">mphase_dusty_gas_shock_tube</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">1</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
    <mesh name="DensityMesh">
 
48
      <from_mesh>
 
49
        <mesh name="CoordinateMesh"/>
 
50
        <mesh_shape>
 
51
          <polynomial_degree>
 
52
            <integer_value rank="0">1</integer_value>
 
53
          </polynomial_degree>
 
54
        </mesh_shape>
 
55
        <stat>
 
56
          <exclude_from_stat/>
 
57
        </stat>
 
58
      </from_mesh>
 
59
    </mesh>
 
60
    <quadrature>
 
61
      <degree>
 
62
        <integer_value rank="0">5</integer_value>
 
63
      </degree>
 
64
    </quadrature>
 
65
  </geometry>
 
66
  <io>
 
67
    <dump_format>
 
68
      <string_value>vtk</string_value>
 
69
    </dump_format>
 
70
    <dump_period>
 
71
      <constant>
 
72
        <real_value rank="0">0.0</real_value>
 
73
      </constant>
 
74
    </dump_period>
 
75
    <output_mesh name="VelocityMesh"/>
 
76
    <stat>
 
77
      <output_at_start/>
 
78
    </stat>
 
79
  </io>
 
80
  <timestepping>
 
81
    <current_time>
 
82
      <real_value rank="0">0.0</real_value>
 
83
    </current_time>
 
84
    <timestep>
 
85
      <real_value rank="0">0.000001</real_value>
 
86
    </timestep>
 
87
    <finish_time>
 
88
      <real_value rank="0">0.00037772998222</real_value>
 
89
      <comment>This corresponds to a non-dimensional time of tau = 4. Multiplying this by the reference time of l/u_ref gives 0.00037772998222 seconds.
 
90
 
 
91
l = (4/3)*(rho_p/rho_g)*d
 
92
u_ref = sqrt(p/rho_g)</comment>
 
93
    </finish_time>
 
94
    <nonlinear_iterations>
 
95
      <integer_value rank="0">3</integer_value>
 
96
      <tolerance>
 
97
        <real_value rank="0">1.0e-9</real_value>
 
98
        <infinity_norm/>
 
99
      </tolerance>
 
100
    </nonlinear_iterations>
 
101
  </timestepping>
 
102
  <material_phase name="Air">
 
103
    <equation_of_state>
 
104
      <compressible>
 
105
        <stiffened_gas>
 
106
          <ratio_specific_heats>
 
107
            <real_value rank="0">1.4</real_value>
 
108
          </ratio_specific_heats>
 
109
        </stiffened_gas>
 
110
      </compressible>
 
111
    </equation_of_state>
 
112
    <scalar_field name="Pressure" rank="0">
 
113
      <prognostic>
 
114
        <mesh name="PressureMesh"/>
 
115
        <spatial_discretisation>
 
116
          <continuous_galerkin>
 
117
            <remove_stabilisation_term/>
 
118
          </continuous_galerkin>
 
119
        </spatial_discretisation>
 
120
        <scheme>
 
121
          <poisson_pressure_solution>
 
122
            <string_value lines="1">never</string_value>
 
123
          </poisson_pressure_solution>
 
124
          <use_compressible_projection_method/>
 
125
        </scheme>
 
126
        <solver>
 
127
          <iterative_method name="preonly"/>
 
128
          <preconditioner name="lu"/>
 
129
          <relative_error>
 
130
            <real_value rank="0">1e-7</real_value>
 
131
          </relative_error>
 
132
          <max_iterations>
 
133
            <integer_value rank="0">10000</integer_value>
 
134
          </max_iterations>
 
135
          <never_ignore_solver_failures/>
 
136
          <diagnostics>
 
137
            <monitors/>
 
138
          </diagnostics>
 
139
        </solver>
 
140
        <output/>
 
141
        <stat/>
 
142
        <convergence>
 
143
          <include_in_convergence/>
 
144
        </convergence>
 
145
        <detectors>
 
146
          <exclude_from_detectors/>
 
147
        </detectors>
 
148
        <steady_state>
 
149
          <include_in_steady_state/>
 
150
        </steady_state>
 
151
        <no_interpolation/>
 
152
      </prognostic>
 
153
    </scalar_field>
 
154
    <scalar_field name="Density" rank="0">
 
155
      <prognostic>
 
156
        <mesh name="PressureMesh"/>
 
157
        <spatial_discretisation>
 
158
          <continuous_galerkin>
 
159
            <stabilisation>
 
160
              <no_stabilisation/>
 
161
            </stabilisation>
 
162
            <advection_terms/>
 
163
            <mass_terms/>
 
164
          </continuous_galerkin>
 
165
          <conservative_advection>
 
166
            <real_value rank="0">0.0</real_value>
 
167
          </conservative_advection>
 
168
        </spatial_discretisation>
 
169
        <temporal_discretisation>
 
170
          <theta>
 
171
            <real_value rank="0">1.0</real_value>
 
172
          </theta>
 
173
        </temporal_discretisation>
 
174
        <initial_condition name="WholeMesh">
 
175
          <python>
 
176
            <string_value lines="20" type="code" language="python">def val(X,t):
 
177
  if(X[0] &lt;= 0):
 
178
    return 12.3
 
179
  else:
 
180
    return 1.23</string_value>
 
181
          </python>
 
182
        </initial_condition>
 
183
        <output/>
 
184
        <stat/>
 
185
        <convergence>
 
186
          <include_in_convergence/>
 
187
        </convergence>
 
188
        <detectors>
 
189
          <include_in_detectors/>
 
190
        </detectors>
 
191
        <steady_state>
 
192
          <include_in_steady_state/>
 
193
        </steady_state>
 
194
        <consistent_interpolation/>
 
195
      </prognostic>
 
196
    </scalar_field>
 
197
    <vector_field name="Velocity" rank="1">
 
198
      <prognostic>
 
199
        <mesh name="VelocityMesh"/>
 
200
        <equation name="LinearMomentum"/>
 
201
        <spatial_discretisation>
 
202
          <continuous_galerkin>
 
203
            <stabilisation>
 
204
              <streamline_upwind>
 
205
                <nu_bar_optimal/>
 
206
                <nu_scale name="0.5">
 
207
                  <real_value shape="1" rank="0">0.5</real_value>
 
208
                </nu_scale>
 
209
              </streamline_upwind>
 
210
            </stabilisation>
 
211
            <mass_terms>
 
212
              <lump_mass_matrix/>
 
213
            </mass_terms>
 
214
            <advection_terms/>
 
215
            <stress_terms>
 
216
              <tensor_form/>
 
217
            </stress_terms>
 
218
          </continuous_galerkin>
 
219
          <conservative_advection>
 
220
            <real_value rank="0">0.0</real_value>
 
221
          </conservative_advection>
 
222
        </spatial_discretisation>
 
223
        <temporal_discretisation>
 
224
          <theta>
 
225
            <real_value rank="0">1.0</real_value>
 
226
          </theta>
 
227
          <relaxation>
 
228
            <real_value rank="0">0.5</real_value>
 
229
          </relaxation>
 
230
        </temporal_discretisation>
 
231
        <solver>
 
232
          <iterative_method name="preonly"/>
 
233
          <preconditioner name="lu"/>
 
234
          <relative_error>
 
235
            <real_value rank="0">1e-7</real_value>
 
236
          </relative_error>
 
237
          <max_iterations>
 
238
            <integer_value rank="0">10000</integer_value>
 
239
          </max_iterations>
 
240
          <never_ignore_solver_failures/>
 
241
          <diagnostics>
 
242
            <monitors/>
 
243
          </diagnostics>
 
244
        </solver>
 
245
        <initial_condition name="WholeMesh">
 
246
          <constant>
 
247
            <real_value shape="1" dim1="dim" rank="1">0.0</real_value>
 
248
          </constant>
 
249
        </initial_condition>
 
250
        <boundary_conditions name="Ends">
 
251
          <surface_ids>
 
252
            <integer_value shape="2" rank="1">1 2</integer_value>
 
253
          </surface_ids>
 
254
          <type name="dirichlet">
 
255
            <align_bc_with_cartesian>
 
256
              <x_component>
 
257
                <constant>
 
258
                  <real_value rank="0">0.0</real_value>
 
259
                </constant>
 
260
              </x_component>
 
261
            </align_bc_with_cartesian>
 
262
          </type>
 
263
        </boundary_conditions>
 
264
        <tensor_field name="Viscosity" rank="2">
 
265
          <prescribed>
 
266
            <value name="WholeMesh">
 
267
              <isotropic>
 
268
                <constant>
 
269
                  <real_value rank="0">1.83e-5</real_value>
 
270
                </constant>
 
271
              </isotropic>
 
272
            </value>
 
273
            <output/>
 
274
          </prescribed>
 
275
        </tensor_field>
 
276
        <output/>
 
277
        <stat>
 
278
          <include_in_stat/>
 
279
          <previous_time_step>
 
280
            <exclude_from_stat/>
 
281
          </previous_time_step>
 
282
          <nonlinear_field>
 
283
            <exclude_from_stat/>
 
284
          </nonlinear_field>
 
285
        </stat>
 
286
        <convergence>
 
287
          <include_in_convergence/>
 
288
        </convergence>
 
289
        <detectors>
 
290
          <include_in_detectors/>
 
291
        </detectors>
 
292
        <steady_state>
 
293
          <include_in_steady_state/>
 
294
        </steady_state>
 
295
        <consistent_interpolation/>
 
296
      </prognostic>
 
297
    </vector_field>
 
298
    <scalar_field name="PhaseVolumeFraction" rank="0">
 
299
      <diagnostic>
 
300
        <mesh name="CoordinateMesh"/>
 
301
        <algorithm name="Internal" material_phase_support="multiple"/>
 
302
        <output/>
 
303
        <stat/>
 
304
        <detectors>
 
305
          <include_in_detectors/>
 
306
        </detectors>
 
307
      </diagnostic>
 
308
    </scalar_field>
 
309
    <scalar_field name="CompressibleContinuityResidual" rank="0">
 
310
      <diagnostic>
 
311
        <mesh name="PressureMesh"/>
 
312
        <algorithm name="Internal" material_phase_support="multiple"/>
 
313
        <output/>
 
314
        <stat/>
 
315
        <detectors>
 
316
          <include_in_detectors/>
 
317
        </detectors>
 
318
        <solver>
 
319
          <iterative_method name="gmres">
 
320
            <restart>
 
321
              <integer_value rank="0">30</integer_value>
 
322
            </restart>
 
323
          </iterative_method>
 
324
          <preconditioner name="sor"/>
 
325
          <relative_error>
 
326
            <real_value rank="0">1.0e-7</real_value>
 
327
          </relative_error>
 
328
          <max_iterations>
 
329
            <integer_value rank="0">10000</integer_value>
 
330
          </max_iterations>
 
331
          <never_ignore_solver_failures/>
 
332
          <diagnostics>
 
333
            <monitors/>
 
334
          </diagnostics>
 
335
        </solver>
 
336
      </diagnostic>
 
337
    </scalar_field>
 
338
    <scalar_field name="InternalEnergy" rank="0">
 
339
      <prognostic>
 
340
        <mesh name="PressureMesh"/>
 
341
        <equation name="InternalEnergy">
 
342
          <density name="Density"/>
 
343
        </equation>
 
344
        <spatial_discretisation>
 
345
          <continuous_galerkin>
 
346
            <stabilisation>
 
347
              <streamline_upwind>
 
348
                <nu_bar_optimal/>
 
349
                <nu_scale name="0.5">
 
350
                  <real_value shape="1" rank="0">0.5</real_value>
 
351
                </nu_scale>
 
352
              </streamline_upwind>
 
353
            </stabilisation>
 
354
            <advection_terms/>
 
355
            <mass_terms/>
 
356
          </continuous_galerkin>
 
357
          <conservative_advection>
 
358
            <real_value rank="0">0.0</real_value>
 
359
          </conservative_advection>
 
360
        </spatial_discretisation>
 
361
        <temporal_discretisation>
 
362
          <theta>
 
363
            <real_value rank="0">1.0</real_value>
 
364
          </theta>
 
365
        </temporal_discretisation>
 
366
        <solver>
 
367
          <iterative_method name="preonly"/>
 
368
          <preconditioner name="lu"/>
 
369
          <relative_error>
 
370
            <real_value rank="0">1.0e-7</real_value>
 
371
          </relative_error>
 
372
          <max_iterations>
 
373
            <integer_value rank="0">10000</integer_value>
 
374
          </max_iterations>
 
375
          <never_ignore_solver_failures/>
 
376
          <diagnostics>
 
377
            <monitors/>
 
378
          </diagnostics>
 
379
        </solver>
 
380
        <initial_condition name="WholeMesh">
 
381
          <constant>
 
382
            <real_value rank="0">205894.308943</real_value>
 
383
            <comment>InternalEnergy from the formula: (3/2)*R*Temperature = (3/2)*287.04*300 K = 203252
 
384
InternalEnergy needed for a Pressure of 1.013e5 = 205894.308943</comment>
 
385
          </constant>
 
386
        </initial_condition>
 
387
        <output/>
 
388
        <stat/>
 
389
        <convergence>
 
390
          <include_in_convergence/>
 
391
        </convergence>
 
392
        <detectors>
 
393
          <include_in_detectors/>
 
394
        </detectors>
 
395
        <steady_state>
 
396
          <include_in_steady_state/>
 
397
        </steady_state>
 
398
        <consistent_interpolation/>
 
399
      </prognostic>
 
400
    </scalar_field>
 
401
    <scalar_field name="CFLNumber" rank="0">
 
402
      <diagnostic>
 
403
        <algorithm name="Internal" material_phase_support="multiple"/>
 
404
        <mesh name="VelocityMesh"/>
 
405
        <output/>
 
406
        <stat/>
 
407
        <convergence>
 
408
          <include_in_convergence/>
 
409
        </convergence>
 
410
        <detectors>
 
411
          <include_in_detectors/>
 
412
        </detectors>
 
413
        <steady_state>
 
414
          <include_in_steady_state/>
 
415
        </steady_state>
 
416
      </diagnostic>
 
417
    </scalar_field>
 
418
    <multiphase_properties>
 
419
      <effective_conductivity>
 
420
        <real_value rank="0">0.026</real_value>
 
421
      </effective_conductivity>
 
422
      <specific_heat>
 
423
        <real_value rank="0">758</real_value>
 
424
      </specific_heat>
 
425
    </multiphase_properties>
 
426
  </material_phase>
 
427
  <material_phase name="Dust">
 
428
    <equation_of_state>
 
429
      <fluids>
 
430
        <linear>
 
431
          <reference_density>
 
432
            <real_value rank="0">2500</real_value>
 
433
          </reference_density>
 
434
        </linear>
 
435
      </fluids>
 
436
    </equation_of_state>
 
437
    <scalar_field name="Pressure" rank="0">
 
438
      <aliased material_phase_name="Air" field_name="Pressure"/>
 
439
    </scalar_field>
 
440
    <scalar_field name="Density" rank="0">
 
441
      <diagnostic>
 
442
        <algorithm name="Internal" material_phase_support="multiple"/>
 
443
        <mesh name="PressureMesh"/>
 
444
        <output/>
 
445
        <stat/>
 
446
        <convergence>
 
447
          <include_in_convergence/>
 
448
        </convergence>
 
449
        <detectors>
 
450
          <include_in_detectors/>
 
451
        </detectors>
 
452
        <steady_state>
 
453
          <include_in_steady_state/>
 
454
        </steady_state>
 
455
      </diagnostic>
 
456
    </scalar_field>
 
457
    <vector_field name="Velocity" rank="1">
 
458
      <prognostic>
 
459
        <mesh name="VelocityMesh"/>
 
460
        <equation name="LinearMomentum"/>
 
461
        <spatial_discretisation>
 
462
          <continuous_galerkin>
 
463
            <stabilisation>
 
464
              <no_stabilisation/>
 
465
            </stabilisation>
 
466
            <mass_terms>
 
467
              <lump_mass_matrix/>
 
468
            </mass_terms>
 
469
            <advection_terms/>
 
470
            <stress_terms>
 
471
              <tensor_form/>
 
472
            </stress_terms>
 
473
          </continuous_galerkin>
 
474
          <conservative_advection>
 
475
            <real_value rank="0">0.0</real_value>
 
476
          </conservative_advection>
 
477
        </spatial_discretisation>
 
478
        <temporal_discretisation>
 
479
          <theta>
 
480
            <real_value rank="0">1.0</real_value>
 
481
          </theta>
 
482
          <relaxation>
 
483
            <real_value rank="0">0.5</real_value>
 
484
          </relaxation>
 
485
        </temporal_discretisation>
 
486
        <solver>
 
487
          <iterative_method name="preonly"/>
 
488
          <preconditioner name="lu"/>
 
489
          <relative_error>
 
490
            <real_value rank="0">1e-7</real_value>
 
491
          </relative_error>
 
492
          <max_iterations>
 
493
            <integer_value rank="0">10000</integer_value>
 
494
          </max_iterations>
 
495
          <never_ignore_solver_failures/>
 
496
          <diagnostics>
 
497
            <monitors/>
 
498
          </diagnostics>
 
499
        </solver>
 
500
        <initial_condition name="WholeMesh">
 
501
          <constant>
 
502
            <real_value shape="1" dim1="dim" rank="1">0.0</real_value>
 
503
          </constant>
 
504
        </initial_condition>
 
505
        <boundary_conditions name="Ends">
 
506
          <surface_ids>
 
507
            <integer_value shape="2" rank="1">1 2</integer_value>
 
508
          </surface_ids>
 
509
          <type name="dirichlet">
 
510
            <align_bc_with_cartesian>
 
511
              <x_component>
 
512
                <constant>
 
513
                  <real_value rank="0">0.0</real_value>
 
514
                </constant>
 
515
              </x_component>
 
516
            </align_bc_with_cartesian>
 
517
          </type>
 
518
        </boundary_conditions>
 
519
        <output/>
 
520
        <stat>
 
521
          <include_in_stat/>
 
522
          <previous_time_step>
 
523
            <exclude_from_stat/>
 
524
          </previous_time_step>
 
525
          <nonlinear_field>
 
526
            <exclude_from_stat/>
 
527
          </nonlinear_field>
 
528
        </stat>
 
529
        <convergence>
 
530
          <include_in_convergence/>
 
531
        </convergence>
 
532
        <detectors>
 
533
          <include_in_detectors/>
 
534
        </detectors>
 
535
        <steady_state>
 
536
          <include_in_steady_state/>
 
537
        </steady_state>
 
538
        <consistent_interpolation/>
 
539
      </prognostic>
 
540
    </vector_field>
 
541
    <scalar_field name="PhaseVolumeFraction" rank="0">
 
542
      <prognostic>
 
543
        <mesh name="CoordinateMesh"/>
 
544
        <equation name="AdvectionDiffusion"/>
 
545
        <spatial_discretisation>
 
546
          <control_volumes>
 
547
            <face_value name="FiniteElement">
 
548
              <limit_face_value>
 
549
                <limiter name="Sweby"/>
 
550
              </limit_face_value>
 
551
            </face_value>
 
552
            <diffusion_scheme name="BassiRebay"/>
 
553
          </control_volumes>
 
554
          <conservative_advection>
 
555
            <real_value rank="0">1.0</real_value>
 
556
          </conservative_advection>
 
557
        </spatial_discretisation>
 
558
        <temporal_discretisation>
 
559
          <theta>
 
560
            <real_value rank="0">1.0</real_value>
 
561
          </theta>
 
562
        </temporal_discretisation>
 
563
        <solver>
 
564
          <iterative_method name="preonly"/>
 
565
          <preconditioner name="lu"/>
 
566
          <relative_error>
 
567
            <real_value rank="0">1.0e-7</real_value>
 
568
          </relative_error>
 
569
          <max_iterations>
 
570
            <integer_value rank="0">10000</integer_value>
 
571
          </max_iterations>
 
572
          <never_ignore_solver_failures/>
 
573
          <diagnostics>
 
574
            <monitors/>
 
575
          </diagnostics>
 
576
        </solver>
 
577
        <initial_condition name="WholeMesh">
 
578
          <python>
 
579
            <string_value lines="20" type="code" language="python">def val(X,t):
 
580
  if(X[0] &lt;= 0):
 
581
    return 4.92e-8
 
582
  else:
 
583
    return 4.92e-4</string_value>
 
584
            <comment>alpha_p = 0.000492 in dusty gas region.
 
585
 
 
586
We know that the ratio of gas density (rho_g) to particle mass concentration (sigma_p) is assumed to be 1.0, so rho_g = sigma_p = 1.23 kg/m^3.
 
587
The normalised mass concentration is 1.0 in the dusty gas region, so we divide 1.23 by the particle density (2500 kg/m^3) get obtain the particle phase volume fraction.</comment>
 
588
          </python>
 
589
        </initial_condition>
 
590
        <boundary_conditions name="All">
 
591
          <surface_ids>
 
592
            <integer_value shape="2" rank="1">1 2</integer_value>
 
593
          </surface_ids>
 
594
          <type name="zero_flux"/>
 
595
        </boundary_conditions>
 
596
        <output/>
 
597
        <stat/>
 
598
        <convergence>
 
599
          <include_in_convergence/>
 
600
        </convergence>
 
601
        <detectors>
 
602
          <include_in_detectors/>
 
603
        </detectors>
 
604
        <steady_state>
 
605
          <include_in_steady_state/>
 
606
        </steady_state>
 
607
        <consistent_interpolation/>
 
608
      </prognostic>
 
609
    </scalar_field>
 
610
    <scalar_field name="ParticleReynoldsNumber" rank="0">
 
611
      <diagnostic>
 
612
        <algorithm name="particle_reynolds_number" material_phase_support="multiple">
 
613
          <depends>
 
614
            <string_value lines="1">Dust::Velocity,Air::Velocity,Air::PhaseVolumeFraction,Air::Density</string_value>
 
615
          </depends>
 
616
          <particle_diameter>
 
617
            <real_value rank="0">10e-6</real_value>
 
618
          </particle_diameter>
 
619
          <continuous_phase_name>Air</continuous_phase_name>
 
620
        </algorithm>
 
621
        <mesh name="VelocityMesh"/>
 
622
        <output/>
 
623
        <stat/>
 
624
        <convergence>
 
625
          <include_in_convergence/>
 
626
        </convergence>
 
627
        <detectors>
 
628
          <include_in_detectors/>
 
629
        </detectors>
 
630
        <steady_state>
 
631
          <include_in_steady_state/>
 
632
        </steady_state>
 
633
      </diagnostic>
 
634
    </scalar_field>
 
635
    <scalar_field name="CFLNumber" rank="0">
 
636
      <diagnostic>
 
637
        <algorithm name="Internal" material_phase_support="multiple"/>
 
638
        <mesh name="VelocityMesh"/>
 
639
        <output/>
 
640
        <stat/>
 
641
        <convergence>
 
642
          <include_in_convergence/>
 
643
        </convergence>
 
644
        <detectors>
 
645
          <include_in_detectors/>
 
646
        </detectors>
 
647
        <steady_state>
 
648
          <include_in_steady_state/>
 
649
        </steady_state>
 
650
      </diagnostic>
 
651
    </scalar_field>
 
652
    <scalar_field name="InternalEnergy" rank="0">
 
653
      <prognostic>
 
654
        <mesh name="PressureMesh"/>
 
655
        <equation name="InternalEnergy">
 
656
          <density name="Density"/>
 
657
        </equation>
 
658
        <spatial_discretisation>
 
659
          <continuous_galerkin>
 
660
            <stabilisation>
 
661
              <streamline_upwind>
 
662
                <nu_bar_optimal/>
 
663
                <nu_scale name="0.5">
 
664
                  <real_value shape="1" rank="0">0.5</real_value>
 
665
                </nu_scale>
 
666
              </streamline_upwind>
 
667
            </stabilisation>
 
668
            <advection_terms/>
 
669
            <mass_terms/>
 
670
          </continuous_galerkin>
 
671
          <conservative_advection>
 
672
            <real_value rank="0">0.0</real_value>
 
673
          </conservative_advection>
 
674
        </spatial_discretisation>
 
675
        <temporal_discretisation>
 
676
          <theta>
 
677
            <real_value rank="0">1.0</real_value>
 
678
          </theta>
 
679
        </temporal_discretisation>
 
680
        <solver>
 
681
          <iterative_method name="preonly"/>
 
682
          <preconditioner name="lu"/>
 
683
          <relative_error>
 
684
            <real_value rank="0">1.0e-7</real_value>
 
685
          </relative_error>
 
686
          <max_iterations>
 
687
            <integer_value rank="0">10000</integer_value>
 
688
          </max_iterations>
 
689
          <never_ignore_solver_failures/>
 
690
          <diagnostics>
 
691
            <monitors/>
 
692
          </diagnostics>
 
693
        </solver>
 
694
        <initial_condition name="WholeMesh">
 
695
          <constant>
 
696
            <real_value rank="0">205894.308943</real_value>
 
697
            <comment>InternalEnergy from the formula: (3/2)*R*Temperature = (3/2)*287.04*300 K = 203252
 
698
InternalEnergy needed for a Pressure of 1.013e5 = 205894.308943</comment>
 
699
          </constant>
 
700
        </initial_condition>
 
701
        <output/>
 
702
        <stat/>
 
703
        <convergence>
 
704
          <include_in_convergence/>
 
705
        </convergence>
 
706
        <detectors>
 
707
          <include_in_detectors/>
 
708
        </detectors>
 
709
        <steady_state>
 
710
          <include_in_steady_state/>
 
711
        </steady_state>
 
712
        <consistent_interpolation/>
 
713
      </prognostic>
 
714
    </scalar_field>
 
715
    <multiphase_properties>
 
716
      <particle_diameter>
 
717
        <real_value rank="0">10e-6</real_value>
 
718
      </particle_diameter>
 
719
      <specific_heat>
 
720
        <real_value rank="0">758</real_value>
 
721
        <comment>1406</comment>
 
722
      </specific_heat>
 
723
    </multiphase_properties>
 
724
  </material_phase>
 
725
  <multiphase_interaction>
 
726
    <fluid_particle_drag>
 
727
      <drag_correlation name="wen_yu"/>
 
728
    </fluid_particle_drag>
 
729
    <heat_transfer>
 
730
      <heat_transfer_coefficient name="gunn"/>
 
731
    </heat_transfer>
 
732
  </multiphase_interaction>
 
733
</fluidity_options>