1
<?xml version='1.0' encoding='utf-8'?>
4
<string_value lines="1">MMS_D_cv</string_value>
7
<string_value lines="1">fluids</string_value>
11
<integer_value rank="0">3</integer_value>
13
<mesh name="CoordinateMesh">
14
<from_file file_name="MMS_D">
21
<mesh name="VelocityMesh">
23
<mesh name="CoordinateMesh"/>
26
<integer_value rank="0">1</integer_value>
36
<integer_value rank="0">5</integer_value>
38
<controlvolume_surface_degree>
39
<integer_value rank="0">5</integer_value>
40
</controlvolume_surface_degree>
45
<string_value>vtk</string_value>
49
<real_value rank="0">1000.0</real_value>
52
<output_mesh name="VelocityMesh"/>
57
<real_value rank="0">0.0</real_value>
60
<real_value rank="0">0.0010676037854987386</real_value>
61
<comment>0.0010676037854987386 gives a cv cfl number of 0.5</comment>
64
<real_value rank="0">5.0</real_value>
68
<real_value rank="0">1.E-10</real_value>
74
<material_phase name="Burgers">
75
<vector_field name="Velocity" rank="1">
77
<mesh name="VelocityMesh"/>
78
<value name="WholeMesh">
80
<string_value lines="20" type="python">def val(XX, t):
81
from math import sin,cos
89
return (u, v, w)</string_value>
97
<exclude_from_detectors/>
101
<scalar_field name="NumericalSolution" rank="0">
103
<mesh name="VelocityMesh"/>
104
<equation name="AdvectionDiffusion"/>
105
<spatial_discretisation>
108
<exclude_mass_terms/>
110
<face_value name="FirstOrderUpwind"/>
111
<diffusion_scheme name="BassiRebay"/>
113
<conservative_advection>
114
<real_value rank="0">0.0</real_value>
115
</conservative_advection>
116
</spatial_discretisation>
117
<temporal_discretisation>
119
<real_value rank="0">1.0</real_value>
121
</temporal_discretisation>
123
<iterative_method name="gmres">
125
<integer_value rank="0">30</integer_value>
128
<preconditioner name="sor"/>
130
<real_value rank="0">1.E-10</real_value>
133
<integer_value rank="0">350</integer_value>
135
<never_ignore_solver_failures/>
140
<initial_condition name="WholeMesh">
142
<real_value rank="0">0.0</real_value>
145
<boundary_conditions name="sides">
147
<integer_value shape="2" rank="1">7 10</integer_value>
149
<type name="dirichlet">
151
<string_value lines="20" type="python">def val(XX, t):
152
from math import sin,cos,sqrt
158
u = sin(25*x*y+omega*t)-2*y/(sqrt(x));
159
return u</string_value>
162
</boundary_conditions>
163
<scalar_field name="Source" rank="0">
165
<value name="WholeMesh">
167
<string_value lines="20" type="python">def val(XX, t):
168
from math import sin,cos,sqrt
176
S = (25*y*cos(25*x*y) + 1.0*y/x1p5)*sin(5*(y2 + x2)) + (25*x*cos(25*x*y) - 2/xp5)*cos(3*(x2 - y2));
177
return S</string_value>
183
<exclude_from_detectors/>
188
<include_previous_time_step/>
192
<include_in_convergence/>
195
<include_in_detectors/>
198
<include_in_steady_state/>
200
<consistent_interpolation/>
203
<scalar_field name="AnalyticalSolution" rank="0">
205
<mesh name="VelocityMesh"/>
206
<value name="WholeMesh">
208
<string_value lines="20" type="python">def val(XX, t):
209
from math import sin,cos,sqrt
213
u = sin(25*x*y + omega*t) - 2*y/(sqrt(x));
214
return u</string_value>
220
<exclude_from_detectors/>
224
<scalar_field name="AbsoluteDifference" rank="0">
225
<diagnostic field_name_b="NumericalSolution" field_name_a="AnalyticalSolution">
226
<algorithm name="Internal" material_phase_support="multiple"/>
227
<mesh name="VelocityMesh"/>
233
<include_in_convergence/>
236
<include_in_detectors/>
239
<include_in_steady_state/>