1
<?xml version='1.0' encoding='utf-8'?>
4
<string_value lines="1">backward_facing_step_2d_kim_kepsilon</string_value>
7
<string_value lines="1">fluids</string_value>
11
<integer_value rank="0">2</integer_value>
13
<mesh name="CoordinateMesh">
14
<from_file file_name="step2d-kim-med">
15
<format name="triangle"/>
21
<mesh name="VelocityMesh">
23
<mesh name="CoordinateMesh"/>
29
<mesh name="PressureMesh">
31
<mesh name="CoordinateMesh"/>
39
<integer_value rank="0">3</integer_value>
45
<string_value>vtk</string_value>
47
<dump_period_in_timesteps>
49
<integer_value rank="0">25</integer_value>
51
</dump_period_in_timesteps>
52
<output_mesh name="CoordinateMesh"/>
57
<real_value rank="0">0.0</real_value>
60
<real_value rank="0">0.01</real_value>
63
<real_value rank="0">300.0</real_value>
65
<nonlinear_iterations>
66
<integer_value rank="0">2</integer_value>
67
</nonlinear_iterations>
70
<real_value rank="0">1.e-4</real_value>
78
<material_phase name="Fluid">
79
<subgridscale_parameterisations>
81
<scalar_field name="TurbulentKineticEnergy" rank="0">
83
<mesh name="VelocityMesh"/>
84
<equation name="KEpsilon"/>
85
<spatial_discretisation>
87
<face_value name="FiniteElement">
89
<limiter name="Sweby">
90
<project_upwind_value_from_point>
91
<store_upwind_elements>
92
<store_upwind_quadrature/>
93
</store_upwind_elements>
94
</project_upwind_value_from_point>
98
<diffusion_scheme name="ElementGradient"/>
100
<conservative_advection>
101
<real_value rank="0">0</real_value>
102
</conservative_advection>
103
</spatial_discretisation>
104
<temporal_discretisation>
106
<real_value rank="0">1.0</real_value>
109
<number_advection_iterations>
110
<integer_value rank="0">3</integer_value>
111
</number_advection_iterations>
114
</temporal_discretisation>
116
<iterative_method name="cg"/>
117
<preconditioner name="hypre">
118
<hypre_type name="boomeramg"/>
121
<real_value rank="0">1.e-7</real_value>
124
<integer_value rank="0">1000</integer_value>
126
<never_ignore_solver_failures/>
131
<initial_condition name="WholeMesh">
133
<real_value rank="0">1.e-4</real_value>
136
<boundary_conditions name="inflow">
138
<integer_value shape="1" rank="1">9</integer_value>
140
<type name="dirichlet">
142
<real_value rank="0">1.e-4</real_value>
145
</boundary_conditions>
146
<boundary_conditions name="walls">
148
<integer_value shape="2" rank="1">10 12</integer_value>
150
<type name="k_epsilon">
151
<string_value>high_Re</string_value>
153
</boundary_conditions>
154
<tensor_field name="Diffusivity" rank="2">
156
<algorithm name="Internal" material_phase_support="multiple"/>
163
<scalar_field name="Source" rank="0">
165
<algorithm name="Internal" material_phase_support="multiple"/>
169
<include_in_detectors/>
173
<scalar_field name="Absorption" rank="0">
175
<algorithm name="Internal" material_phase_support="multiple"/>
179
<include_in_detectors/>
186
<include_in_convergence/>
189
<include_in_detectors/>
192
<include_in_steady_state/>
194
<consistent_interpolation/>
197
<scalar_field name="TurbulentDissipation" rank="0">
199
<mesh name="VelocityMesh"/>
200
<equation name="KEpsilon"/>
201
<spatial_discretisation>
203
<face_value name="FiniteElement">
205
<limiter name="Sweby">
206
<project_upwind_value_from_point>
207
<store_upwind_elements>
208
<store_upwind_quadrature/>
209
</store_upwind_elements>
210
</project_upwind_value_from_point>
214
<diffusion_scheme name="ElementGradient"/>
216
<conservative_advection>
217
<real_value rank="0">0</real_value>
218
</conservative_advection>
219
</spatial_discretisation>
220
<temporal_discretisation>
222
<real_value rank="0">1.0</real_value>
225
<number_advection_iterations>
226
<integer_value rank="0">3</integer_value>
227
</number_advection_iterations>
230
</temporal_discretisation>
232
<iterative_method name="cg"/>
233
<preconditioner name="hypre">
234
<hypre_type name="boomeramg"/>
237
<real_value rank="0">1.e-7</real_value>
240
<integer_value rank="0">1000</integer_value>
242
<never_ignore_solver_failures/>
247
<initial_condition name="WholeMesh">
249
<real_value rank="0">1.e-5</real_value>
252
<boundary_conditions name="inflow">
254
<integer_value shape="1" rank="1">9</integer_value>
256
<type name="dirichlet">
258
<real_value rank="0">1.e-5</real_value>
261
</boundary_conditions>
262
<boundary_conditions name="walls">
264
<integer_value shape="2" rank="1">10 12</integer_value>
266
<type name="k_epsilon">
267
<string_value>high_Re</string_value>
269
</boundary_conditions>
270
<tensor_field name="Diffusivity" rank="2">
272
<algorithm name="Internal" material_phase_support="multiple"/>
279
<scalar_field name="Source" rank="0">
281
<algorithm name="Internal" material_phase_support="multiple"/>
285
<include_in_detectors/>
289
<scalar_field name="Absorption" rank="0">
291
<algorithm name="Internal" material_phase_support="multiple"/>
295
<include_in_detectors/>
302
<include_in_convergence/>
305
<include_in_detectors/>
308
<include_in_steady_state/>
310
<consistent_interpolation/>
313
<tensor_field name="BackgroundViscosity" rank="2">
315
<mesh name="VelocityMesh"/>
316
<value name="WholeMesh">
317
<anisotropic_symmetric>
319
<real_value symmetric="true" dim2="dim" shape="2 2" dim1="dim" rank="2">1.74e-5 1.74e-5 1.74e-5 1.74e-5</real_value>
321
</anisotropic_symmetric>
326
<tensor_field name="EddyViscosity" rank="2">
328
<algorithm name="Internal" material_phase_support="multiple"/>
329
<mesh name="VelocityMesh"/>
336
<scalar_field name="ScalarEddyViscosity" rank="0">
338
<algorithm name="Internal" material_phase_support="multiple"/>
339
<mesh name="VelocityMesh"/>
343
<include_in_convergence/>
346
<include_in_detectors/>
349
<include_in_steady_state/>
353
<scalar_field name="LengthScale" rank="0">
355
<algorithm name="Internal" material_phase_support="multiple"/>
356
<mesh name="VelocityMesh"/>
360
<include_in_convergence/>
363
<include_in_detectors/>
366
<include_in_steady_state/>
370
<scalar_field name="f_1" rank="0">
372
<algorithm name="Internal" material_phase_support="multiple"/>
373
<mesh name="VelocityMesh"/>
377
<include_in_convergence/>
380
<include_in_detectors/>
383
<include_in_steady_state/>
387
<scalar_field name="f_2" rank="0">
389
<algorithm name="Internal" material_phase_support="multiple"/>
390
<mesh name="VelocityMesh"/>
394
<include_in_convergence/>
397
<include_in_detectors/>
400
<include_in_steady_state/>
404
<scalar_field name="f_mu" rank="0">
406
<algorithm name="Internal" material_phase_support="multiple"/>
407
<mesh name="VelocityMesh"/>
411
<include_in_convergence/>
414
<include_in_detectors/>
417
<include_in_steady_state/>
422
<real_value rank="0">10.0</real_value>
425
<real_value rank="0">.09</real_value>
428
<real_value rank="0">1.44</real_value>
431
<real_value rank="0">1.92</real_value>
434
<real_value rank="0">1</real_value>
437
<real_value rank="0">1.3</real_value>
440
<real_value rank="0">1.0</real_value>
442
<time_discretisation>
444
<real_value rank="0">0.5</real_value>
446
<source_term_implementation>
447
<production_term>source</production_term>
448
<destruction_term>absorbtion</destruction_term>
449
<buoyancy_term>source</buoyancy_term>
450
</source_term_implementation>
451
</time_discretisation>
456
</subgridscale_parameterisations>
457
<scalar_field name="Pressure" rank="0">
459
<mesh name="PressureMesh"/>
460
<spatial_discretisation>
461
<continuous_galerkin>
462
<integrate_continuity_by_parts/>
463
</continuous_galerkin>
464
</spatial_discretisation>
466
<poisson_pressure_solution>
467
<string_value lines="1">only first timestep</string_value>
468
</poisson_pressure_solution>
469
<use_projection_method/>
472
<iterative_method name="cg"/>
473
<preconditioner name="mg"/>
475
<real_value rank="0">1.0E-7</real_value>
478
<integer_value rank="0">10000</integer_value>
480
<never_ignore_solver_failures/>
488
<include_in_convergence/>
491
<exclude_from_detectors/>
494
<include_in_steady_state/>
499
<vector_field name="Velocity" rank="1">
501
<mesh name="VelocityMesh"/>
502
<equation name="Boussinesq"/>
503
<spatial_discretisation>
504
<continuous_galerkin>
513
<partial_stress_form/>
515
</continuous_galerkin>
516
<conservative_advection>
517
<real_value rank="0">0.0</real_value>
518
</conservative_advection>
519
</spatial_discretisation>
520
<temporal_discretisation>
522
<real_value rank="0">0.5</real_value>
525
<real_value rank="0">0.5</real_value>
527
</temporal_discretisation>
529
<iterative_method name="gmres">
531
<integer_value rank="0">30</integer_value>
534
<preconditioner name="sor"/>
536
<real_value rank="0">1.0E-7</real_value>
539
<integer_value rank="0">1000</integer_value>
541
<never_ignore_solver_failures/>
546
<initial_condition name="WholeMesh">
548
<real_value shape="2" dim1="dim" rank="1">0.0 0.0</real_value>
551
<boundary_conditions name="Inflow">
553
<integer_value shape="1" rank="1">9</integer_value>
555
<type name="dirichlet">
556
<align_bc_with_cartesian>
559
<string_value lines="20" type="code" language="python">def val(X, t):
569
return u</string_value>
570
<comment>max velocity = 2.3</comment>
575
<real_value rank="0">0.</real_value>
578
</align_bc_with_cartesian>
580
</boundary_conditions>
581
<boundary_conditions name="walls">
583
<integer_value shape="2" rank="1">10 12</integer_value>
585
<type name="log_law_of_wall">
587
<real_value rank="0">1.e-6</real_value>
590
</boundary_conditions>
591
<boundary_conditions name="no_normal">
593
<integer_value shape="2" rank="1">10 12</integer_value>
595
<type name="dirichlet">
596
<align_bc_with_surface>
599
<real_value rank="0">0.</real_value>
602
</align_bc_with_surface>
604
</boundary_conditions>
605
<tensor_field name="Viscosity" rank="2">
607
<algorithm name="Internal" material_phase_support="multiple"/>
616
<vector_field name="Source" rank="1">
618
<algorithm name="Internal" material_phase_support="multiple"/>
626
<include_in_convergence/>
629
<include_in_detectors/>
632
<include_in_steady_state/>
641
</previous_time_step>
647
<include_in_convergence/>
650
<include_in_detectors/>
653
<include_in_steady_state/>
655
<consistent_interpolation/>
658
<scalar_field name="CFLNumber" rank="0">
660
<algorithm name="Internal" material_phase_support="multiple"/>
661
<mesh name="VelocityMesh"/>
665
<include_in_convergence/>
668
<include_in_detectors/>
671
<exclude_from_steady_state/>