← Back to branch summary

~fluidity-core/fluidity/sea-ice-branch

« back to all changes in this revision

Viewing changes to tests/mphase_dusty_gas_shock_tube/single_phase_frozen_flow_test.flml

  • Committer: Simon Mouradian
  • Date: 2012-10-19 10:35:59 UTC
  • mfrom: (3520.32.371 fluidity)
  • Revision ID: simon.mouradian06@imperial.ac.uk-20121019103559-y36qa47phc69q8sc
merge from trunk

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/mphase_dusty_gas_shock_tube/single_phase_frozen_flow_test.flml
 
1
<?xml version='1.0' encoding='utf-8'?>
 
2
<fluidity_options>
 
3
  <simulation_name>
 
4
    <string_value lines="1">single_phase_frozen_flow_test</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">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
            <factorization_package name="petsc"/>
 
130
          </preconditioner>
 
131
          <relative_error>
 
132
            <real_value rank="0">1e-7</real_value>
 
133
          </relative_error>
 
134
          <max_iterations>
 
135
            <integer_value rank="0">10000</integer_value>
 
136
          </max_iterations>
 
137
          <never_ignore_solver_failures/>
 
138
          <diagnostics>
 
139
            <monitors/>
 
140
          </diagnostics>
 
141
        </solver>
 
142
        <output/>
 
143
        <stat/>
 
144
        <convergence>
 
145
          <include_in_convergence/>
 
146
        </convergence>
 
147
        <detectors>
 
148
          <exclude_from_detectors/>
 
149
        </detectors>
 
150
        <steady_state>
 
151
          <include_in_steady_state/>
 
152
        </steady_state>
 
153
        <no_interpolation/>
 
154
      </prognostic>
 
155
    </scalar_field>
 
156
    <scalar_field name="Density" rank="0">
 
157
      <prognostic>
 
158
        <mesh name="PressureMesh"/>
 
159
        <spatial_discretisation>
 
160
          <continuous_galerkin>
 
161
            <stabilisation>
 
162
              <no_stabilisation/>
 
163
            </stabilisation>
 
164
            <advection_terms/>
 
165
            <mass_terms/>
 
166
          </continuous_galerkin>
 
167
          <conservative_advection>
 
168
            <real_value rank="0">0.0</real_value>
 
169
          </conservative_advection>
 
170
        </spatial_discretisation>
 
171
        <temporal_discretisation>
 
172
          <theta>
 
173
            <real_value rank="0">1.0</real_value>
 
174
          </theta>
 
175
        </temporal_discretisation>
 
176
        <initial_condition name="WholeMesh">
 
177
          <python>
 
178
            <string_value lines="20" type="code" language="python">def val(X,t):
 
179
  if(X[0] &lt;= 0):
 
180
    return 12.3
 
181
  else:
 
182
    return 1.23</string_value>
 
183
          </python>
 
184
        </initial_condition>
 
185
        <output/>
 
186
        <stat/>
 
187
        <convergence>
 
188
          <include_in_convergence/>
 
189
        </convergence>
 
190
        <detectors>
 
191
          <include_in_detectors/>
 
192
        </detectors>
 
193
        <steady_state>
 
194
          <include_in_steady_state/>
 
195
        </steady_state>
 
196
        <consistent_interpolation/>
 
197
      </prognostic>
 
198
    </scalar_field>
 
199
    <vector_field name="Velocity" rank="1">
 
200
      <prognostic>
 
201
        <mesh name="VelocityMesh"/>
 
202
        <equation name="LinearMomentum"/>
 
203
        <spatial_discretisation>
 
204
          <continuous_galerkin>
 
205
            <stabilisation>
 
206
              <streamline_upwind>
 
207
                <nu_bar_optimal/>
 
208
                <nu_scale name="0.5">
 
209
                  <real_value shape="1" rank="0">0.5</real_value>
 
210
                </nu_scale>
 
211
              </streamline_upwind>
 
212
            </stabilisation>
 
213
            <mass_terms>
 
214
              <lump_mass_matrix/>
 
215
            </mass_terms>
 
216
            <advection_terms/>
 
217
            <stress_terms>
 
218
              <tensor_form/>
 
219
            </stress_terms>
 
220
          </continuous_galerkin>
 
221
          <conservative_advection>
 
222
            <real_value rank="0">0.0</real_value>
 
223
          </conservative_advection>
 
224
        </spatial_discretisation>
 
225
        <temporal_discretisation>
 
226
          <theta>
 
227
            <real_value rank="0">1.0</real_value>
 
228
          </theta>
 
229
          <relaxation>
 
230
            <real_value rank="0">0.5</real_value>
 
231
          </relaxation>
 
232
        </temporal_discretisation>
 
233
        <solver>
 
234
          <iterative_method name="preonly"/>
 
235
          <preconditioner name="lu">
 
236
            <factorization_package name="petsc"/>
 
237
          </preconditioner>
 
238
          <relative_error>
 
239
            <real_value rank="0">1e-7</real_value>
 
240
          </relative_error>
 
241
          <max_iterations>
 
242
            <integer_value rank="0">10000</integer_value>
 
243
          </max_iterations>
 
244
          <never_ignore_solver_failures/>
 
245
          <diagnostics>
 
246
            <monitors/>
 
247
          </diagnostics>
 
248
        </solver>
 
249
        <initial_condition name="WholeMesh">
 
250
          <constant>
 
251
            <real_value shape="1" dim1="dim" rank="1">0.0</real_value>
 
252
          </constant>
 
253
        </initial_condition>
 
254
        <boundary_conditions name="Ends">
 
255
          <surface_ids>
 
256
            <integer_value shape="2" rank="1">1 2</integer_value>
 
257
          </surface_ids>
 
258
          <type name="dirichlet">
 
259
            <align_bc_with_cartesian>
 
260
              <x_component>
 
261
                <constant>
 
262
                  <real_value rank="0">0.0</real_value>
 
263
                </constant>
 
264
              </x_component>
 
265
            </align_bc_with_cartesian>
 
266
          </type>
 
267
        </boundary_conditions>
 
268
        <tensor_field name="Viscosity" rank="2">
 
269
          <prescribed>
 
270
            <value name="WholeMesh">
 
271
              <isotropic>
 
272
                <constant>
 
273
                  <real_value rank="0">1.83e-5</real_value>
 
274
                </constant>
 
275
              </isotropic>
 
276
            </value>
 
277
            <output/>
 
278
          </prescribed>
 
279
        </tensor_field>
 
280
        <output/>
 
281
        <stat>
 
282
          <include_in_stat/>
 
283
          <previous_time_step>
 
284
            <exclude_from_stat/>
 
285
          </previous_time_step>
 
286
          <nonlinear_field>
 
287
            <exclude_from_stat/>
 
288
          </nonlinear_field>
 
289
        </stat>
 
290
        <convergence>
 
291
          <include_in_convergence/>
 
292
        </convergence>
 
293
        <detectors>
 
294
          <include_in_detectors/>
 
295
        </detectors>
 
296
        <steady_state>
 
297
          <include_in_steady_state/>
 
298
        </steady_state>
 
299
        <consistent_interpolation/>
 
300
      </prognostic>
 
301
    </vector_field>
 
302
    <scalar_field name="InternalEnergy" rank="0">
 
303
      <prognostic>
 
304
        <mesh name="PressureMesh"/>
 
305
        <equation name="InternalEnergy">
 
306
          <density name="Density"/>
 
307
        </equation>
 
308
        <spatial_discretisation>
 
309
          <continuous_galerkin>
 
310
            <stabilisation>
 
311
              <streamline_upwind>
 
312
                <nu_bar_optimal/>
 
313
                <nu_scale name="0.5">
 
314
                  <real_value shape="1" rank="0">0.5</real_value>
 
315
                </nu_scale>
 
316
              </streamline_upwind>
 
317
            </stabilisation>
 
318
            <advection_terms/>
 
319
            <mass_terms/>
 
320
          </continuous_galerkin>
 
321
          <conservative_advection>
 
322
            <real_value rank="0">0.0</real_value>
 
323
          </conservative_advection>
 
324
        </spatial_discretisation>
 
325
        <temporal_discretisation>
 
326
          <theta>
 
327
            <real_value rank="0">1.0</real_value>
 
328
          </theta>
 
329
        </temporal_discretisation>
 
330
        <solver>
 
331
          <iterative_method name="preonly"/>
 
332
          <preconditioner name="lu">
 
333
            <factorization_package name="petsc"/>
 
334
          </preconditioner>
 
335
          <relative_error>
 
336
            <real_value rank="0">1.0e-7</real_value>
 
337
          </relative_error>
 
338
          <max_iterations>
 
339
            <integer_value rank="0">10000</integer_value>
 
340
          </max_iterations>
 
341
          <never_ignore_solver_failures/>
 
342
          <diagnostics>
 
343
            <monitors/>
 
344
          </diagnostics>
 
345
        </solver>
 
346
        <initial_condition name="WholeMesh">
 
347
          <constant>
 
348
            <real_value rank="0">205894.308943</real_value>
 
349
            <comment>InternalEnergy from the formula: (3/2)*R*Temperature = (3/2)*287.04*300 K = 203252
 
350
InternalEnergy needed for a Pressure of 1.013e5 = 205894.308943</comment>
 
351
          </constant>
 
352
        </initial_condition>
 
353
        <output/>
 
354
        <stat/>
 
355
        <convergence>
 
356
          <include_in_convergence/>
 
357
        </convergence>
 
358
        <detectors>
 
359
          <include_in_detectors/>
 
360
        </detectors>
 
361
        <steady_state>
 
362
          <include_in_steady_state/>
 
363
        </steady_state>
 
364
        <consistent_interpolation/>
 
365
      </prognostic>
 
366
    </scalar_field>
 
367
    <scalar_field name="CFLNumber" rank="0">
 
368
      <diagnostic>
 
369
        <algorithm name="Internal" material_phase_support="multiple"/>
 
370
        <mesh name="VelocityMesh"/>
 
371
        <output/>
 
372
        <stat/>
 
373
        <convergence>
 
374
          <include_in_convergence/>
 
375
        </convergence>
 
376
        <detectors>
 
377
          <include_in_detectors/>
 
378
        </detectors>
 
379
        <steady_state>
 
380
          <include_in_steady_state/>
 
381
        </steady_state>
 
382
      </diagnostic>
 
383
    </scalar_field>
 
384
  </material_phase>
 
385
</fluidity_options>