1
<?xml version='1.0' encoding='utf-8'?>
4
<string_value lines="1">MMS_B</string_value>
7
<string_value lines="1">fluids</string_value>
11
<integer_value rank="0">2</integer_value>
13
<mesh name="CoordinateMesh">
15
<mesh name="InputMesh"/>
21
<mesh name="VelocityMesh">
23
<mesh name="InputMesh"/>
26
<integer_value rank="0">2</integer_value>
34
<mesh name="InputMesh">
35
<from_file file_name="MMS_A">
36
<format name="triangle"/>
44
<integer_value rank="0">5</integer_value>
50
<string_value>vtk</string_value>
54
<real_value rank="0">25</real_value>
57
<output_mesh name="VelocityMesh"/>
62
<real_value rank="0">0.0</real_value>
65
<real_value rank="0">0.01</real_value>
68
<real_value rank="0">1</real_value>
70
<nonlinear_iterations>
71
<integer_value rank="0">2</integer_value>
72
</nonlinear_iterations>
75
<real_value rank="0">1.e-8</real_value>
80
<material_phase name="NS">
81
<vector_field name="Velocity" rank="1">
83
<mesh name="VelocityMesh"/>
84
<value name="WholeMesh">
86
<real_value shape="2" dim1="dim" rank="1">0. 0.</real_value>
94
<exclude_from_detectors/>
98
<scalar_field name="Tracer" rank="0">
100
<mesh name="VelocityMesh"/>
101
<equation name="AdvectionDiffusion"/>
102
<spatial_discretisation>
103
<continuous_galerkin>
109
</continuous_galerkin>
110
<conservative_advection>
111
<real_value rank="0">0.</real_value>
112
</conservative_advection>
113
</spatial_discretisation>
114
<temporal_discretisation>
116
<real_value rank="0">0.5</real_value>
118
</temporal_discretisation>
120
<iterative_method name="gmres">
122
<integer_value rank="0">30</integer_value>
125
<preconditioner name="sor"/>
127
<real_value rank="0">1.e-7</real_value>
130
<integer_value rank="0">1000</integer_value>
132
<never_ignore_solver_failures/>
137
<initial_condition name="WholeMesh">
139
<string_value lines="20" type="code" language="python">def val(X,t):
140
from math import pi, exp
145
return u</string_value>
148
<boundary_conditions name="sides">
150
<integer_value shape="4" rank="1">7 8 9 10</integer_value>
152
<type name="dirichlet">
154
<string_value lines="20" type="code" language="python">def val(X,t):
155
from math import pi, exp
160
return u</string_value>
163
</boundary_conditions>
167
<include_in_convergence/>
170
<include_in_detectors/>
173
<include_in_steady_state/>
175
<consistent_interpolation/>
178
<scalar_field name="FilteredTracer" rank="0">
180
<algorithm source_field_type="component" material_phase_support="single" name="helmholtz_smoothed_scalar" source_field_name="Tracer">
181
<smoothing_length_scale>
182
<real_value symmetric="true" dim2="dim" shape="2 2" dim1="dim" rank="2">0.005 0 0 0.005</real_value>
183
</smoothing_length_scale>
185
<iterative_method name="cg"/>
186
<preconditioner name="sor"/>
188
<real_value rank="0">1.e-7</real_value>
191
<integer_value rank="0">1000</integer_value>
193
<never_ignore_solver_failures/>
199
<mesh name="VelocityMesh"/>
203
<include_in_convergence/>
206
<include_in_detectors/>
209
<include_in_steady_state/>
213
<scalar_field name="Difference" rank="0">
215
<algorithm source_field_2_type="scalar" name="scalar_difference" source_field_1_name="Tracer" source_field_2_name="FilteredTracer" material_phase_support="single" source_field_1_type="scalar"/>
216
<mesh name="VelocityMesh"/>
220
<include_in_convergence/>
223
<include_in_detectors/>
226
<include_in_steady_state/>