← Back to branch summary

~fluidity-core/fluidity/sediment

« back to all changes in this revision

Viewing changes to tests/k_epsilon_scalar_diffusivity/k_epsilon_scalar_diffusivity.flml

  • Committer: Sam Parkinson
  • Date: 2011-12-01 15:18:13 UTC
  • mfrom: (3413.1.30 sediment_restructure)
  • Revision ID: s.parkinson11@imperial.ac.uk-20111201151813-9vygndld22tdzepx
merge with sediment_restructure branch

Show diffs side-by-side

added added

removed removed

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