← Back to branch summary

~f-milthaler/fluidity/fsi-model-stationary-solid-with-velocity

« back to all changes in this revision

Viewing changes to tests/inlet_velocity_bc_compressible/inlet_velocity_bc_compressible_1d.flml

  • Committer: f.milthaler10 at uk
  • Date: 2013-11-06 13:43:56 UTC
  • mfrom: (3463.184.85 fluidity)
  • Revision ID: f.milthaler10@imperial.ac.ic.uk.-20131106134356-v3lw1dheesckywj0
merge from trunk

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/inlet_velocity_bc_compressible/inlet_velocity_bc_compressible_1d.flml
 
1
<?xml version='1.0' encoding='utf-8'?>
 
2
<fluidity_options>
 
3
  <simulation_name>
 
4
    <string_value lines="1">inlet_velocity_bc_compressible_1d</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">0</integer_value>
 
27
          </polynomial_degree>
 
28
        </mesh_shape>
 
29
        <mesh_continuity>
 
30
          <string_value>discontinuous</string_value>
 
31
        </mesh_continuity>
 
32
        <stat>
 
33
          <exclude_from_stat/>
 
34
        </stat>
 
35
      </from_mesh>
 
36
    </mesh>
 
37
    <mesh name="PressureMesh">
 
38
      <from_mesh>
 
39
        <mesh name="CoordinateMesh"/>
 
40
        <mesh_shape>
 
41
          <polynomial_degree>
 
42
            <integer_value rank="0">1</integer_value>
 
43
          </polynomial_degree>
 
44
        </mesh_shape>
 
45
        <stat>
 
46
          <exclude_from_stat/>
 
47
        </stat>
 
48
      </from_mesh>
 
49
    </mesh>
 
50
    <mesh name="DensityMesh">
 
51
      <from_mesh>
 
52
        <mesh name="CoordinateMesh"/>
 
53
        <mesh_shape>
 
54
          <polynomial_degree>
 
55
            <integer_value rank="0">1</integer_value>
 
56
          </polynomial_degree>
 
57
        </mesh_shape>
 
58
        <stat>
 
59
          <exclude_from_stat/>
 
60
        </stat>
 
61
      </from_mesh>
 
62
    </mesh>
 
63
    <mesh name="GeostrophicMesh">
 
64
      <from_mesh>
 
65
        <mesh name="CoordinateMesh"/>
 
66
        <mesh_shape>
 
67
          <polynomial_degree>
 
68
            <integer_value rank="0">1</integer_value>
 
69
          </polynomial_degree>
 
70
        </mesh_shape>
 
71
        <stat>
 
72
          <exclude_from_stat/>
 
73
        </stat>
 
74
      </from_mesh>
 
75
    </mesh>
 
76
    <quadrature>
 
77
      <degree>
 
78
        <integer_value rank="0">4</integer_value>
 
79
      </degree>
 
80
    </quadrature>
 
81
  </geometry>
 
82
  <io>
 
83
    <dump_format>
 
84
      <string_value>vtk</string_value>
 
85
    </dump_format>
 
86
    <dump_period>
 
87
      <constant>
 
88
        <real_value rank="0">0</real_value>
 
89
      </constant>
 
90
    </dump_period>
 
91
    <output_mesh name="PressureMesh"/>
 
92
    <stat>
 
93
      <output_at_start/>
 
94
    </stat>
 
95
  </io>
 
96
  <timestepping>
 
97
    <current_time>
 
98
      <real_value rank="0">0</real_value>
 
99
    </current_time>
 
100
    <timestep>
 
101
      <real_value rank="0">0.01</real_value>
 
102
    </timestep>
 
103
    <finish_time>
 
104
      <real_value rank="0">600.0</real_value>
 
105
    </finish_time>
 
106
    <nonlinear_iterations>
 
107
      <integer_value rank="0">2</integer_value>
 
108
      <tolerance>
 
109
        <real_value rank="0">1.0e-9</real_value>
 
110
        <infinity_norm/>
 
111
      </tolerance>
 
112
    </nonlinear_iterations>
 
113
    <adaptive_timestep>
 
114
      <requested_cfl>
 
115
        <real_value rank="0">0.5</real_value>
 
116
      </requested_cfl>
 
117
      <courant_number name="CFLNumber">
 
118
        <mesh name="VelocityMesh"/>
 
119
      </courant_number>
 
120
      <minimum_timestep>
 
121
        <real_value rank="0">0.001</real_value>
 
122
      </minimum_timestep>
 
123
      <increase_tolerance>
 
124
        <real_value rank="0">1.1</real_value>
 
125
      </increase_tolerance>
 
126
    </adaptive_timestep>
 
127
  </timestepping>
 
128
  <physical_parameters>
 
129
    <gravity>
 
130
      <magnitude>
 
131
        <real_value rank="0">9.8</real_value>
 
132
      </magnitude>
 
133
      <vector_field name="GravityDirection" rank="1">
 
134
        <prescribed>
 
135
          <mesh name="CoordinateMesh"/>
 
136
          <value name="WholeMesh">
 
137
            <constant>
 
138
              <real_value shape="1" dim1="dim" rank="1">-1</real_value>
 
139
            </constant>
 
140
          </value>
 
141
          <output/>
 
142
          <stat>
 
143
            <include_in_stat/>
 
144
          </stat>
 
145
          <detectors>
 
146
            <exclude_from_detectors/>
 
147
          </detectors>
 
148
        </prescribed>
 
149
      </vector_field>
 
150
    </gravity>
 
151
  </physical_parameters>
 
152
  <material_phase name="Gas">
 
153
    <equation_of_state>
 
154
      <compressible>
 
155
        <stiffened_gas>
 
156
          <ratio_specific_heats>
 
157
            <real_value rank="0">1.33</real_value>
 
158
          </ratio_specific_heats>
 
159
        </stiffened_gas>
 
160
        <subtract_out_reference_profile/>
 
161
      </compressible>
 
162
    </equation_of_state>
 
163
    <scalar_field name="Pressure" rank="0">
 
164
      <prognostic>
 
165
        <mesh name="PressureMesh"/>
 
166
        <spatial_discretisation>
 
167
          <continuous_galerkin>
 
168
            <remove_stabilisation_term/>
 
169
            <integrate_continuity_by_parts/>
 
170
          </continuous_galerkin>
 
171
        </spatial_discretisation>
 
172
        <scheme>
 
173
          <poisson_pressure_solution>
 
174
            <string_value lines="1">never</string_value>
 
175
          </poisson_pressure_solution>
 
176
          <use_projection_method/>
 
177
        </scheme>
 
178
        <solver>
 
179
          <iterative_method name="gmres">
 
180
            <restart>
 
181
              <integer_value rank="0">30</integer_value>
 
182
            </restart>
 
183
          </iterative_method>
 
184
          <preconditioner name="sor"/>
 
185
          <relative_error>
 
186
            <real_value rank="0">1.0e-7</real_value>
 
187
          </relative_error>
 
188
          <max_iterations>
 
189
            <integer_value rank="0">1000</integer_value>
 
190
          </max_iterations>
 
191
          <never_ignore_solver_failures/>
 
192
          <diagnostics>
 
193
            <monitors/>
 
194
          </diagnostics>
 
195
        </solver>
 
196
        <initial_condition name="WholeMesh">
 
197
          <python>
 
198
            <string_value lines="20" type="code" language="python">def val(X,t):
 
199
  from math import exp
 
200
  T0 = 288.15 # K
 
201
  R = 287.058 # Specific gas constant in J/Kg
 
202
  g = 9.8 # m/s**2
 
203
  h = X[0] # m
 
204
  p0 = 101325
 
205
  e = T0*1406.0
 
206
  gamma = 1.33
 
207
  p = p0*exp(-(g*h)/((gamma-1.0)*e))
 
208
  return p</string_value>
 
209
          </python>
 
210
        </initial_condition>
 
211
        <boundary_conditions name="Right">
 
212
          <surface_ids>
 
213
            <integer_value shape="1" rank="1">2</integer_value>
 
214
          </surface_ids>
 
215
          <type name="dirichlet">
 
216
            <python>
 
217
              <string_value lines="20" type="code" language="python">def val(X,t):
 
218
  from math import exp
 
219
  T0 = 288.15 # K
 
220
  R = 287.058 # Specific gas constant in J/Kg
 
221
  g = 9.8 # m/s**2
 
222
  h = X[0] # m
 
223
  p0 = 101325
 
224
  e = T0*1406.0
 
225
  gamma = 1.33
 
226
  p = p0*exp(-(g*h)/((gamma-1.0)*e))
 
227
  return p</string_value>
 
228
            </python>
 
229
          </type>
 
230
        </boundary_conditions>
 
231
        <output/>
 
232
        <stat/>
 
233
        <convergence>
 
234
          <include_in_convergence/>
 
235
        </convergence>
 
236
        <detectors>
 
237
          <exclude_from_detectors/>
 
238
        </detectors>
 
239
        <steady_state>
 
240
          <include_in_steady_state/>
 
241
        </steady_state>
 
242
        <no_interpolation/>
 
243
      </prognostic>
 
244
    </scalar_field>
 
245
    <scalar_field name="Density" rank="0">
 
246
      <prognostic>
 
247
        <mesh name="DensityMesh"/>
 
248
        <spatial_discretisation>
 
249
          <control_volumes>
 
250
            <face_value name="FirstOrderUpwind"/>
 
251
          </control_volumes>
 
252
          <conservative_advection>
 
253
            <real_value rank="0">1.0</real_value>
 
254
          </conservative_advection>
 
255
        </spatial_discretisation>
 
256
        <temporal_discretisation>
 
257
          <theta>
 
258
            <real_value rank="0">1.0</real_value>
 
259
          </theta>
 
260
        </temporal_discretisation>
 
261
        <initial_condition name="WholeMesh">
 
262
          <python>
 
263
            <string_value lines="20" type="code" language="python">def val(X,t):
 
264
  from math import exp
 
265
  T0 = 288.15 # K
 
266
  R = 287.058 # Specific gas constant in J/Kg
 
267
  g = 9.8 # m/s**2
 
268
  h = X[0] # m
 
269
  p0 = 101325
 
270
  e = T0*1406.0
 
271
  gamma = 1.33
 
272
  p = p0*exp(-(g*h)/((gamma-1.0)*e))
 
273
  rho = p/((gamma-1.0)*e)
 
274
  return rho</string_value>
 
275
          </python>
 
276
        </initial_condition>
 
277
        <boundary_conditions name="Left">
 
278
          <surface_ids>
 
279
            <integer_value shape="1" rank="1">1</integer_value>
 
280
          </surface_ids>
 
281
          <type name="dirichlet">
 
282
            <constant>
 
283
              <real_value rank="0">0.181985215</real_value>
 
284
            </constant>
 
285
          </type>
 
286
        </boundary_conditions>
 
287
        <output/>
 
288
        <stat/>
 
289
        <convergence>
 
290
          <include_in_convergence/>
 
291
        </convergence>
 
292
        <detectors>
 
293
          <include_in_detectors/>
 
294
        </detectors>
 
295
        <steady_state>
 
296
          <include_in_steady_state/>
 
297
        </steady_state>
 
298
        <consistent_interpolation/>
 
299
      </prognostic>
 
300
    </scalar_field>
 
301
    <vector_field name="Velocity" rank="1">
 
302
      <prognostic>
 
303
        <mesh name="VelocityMesh"/>
 
304
        <equation name="LinearMomentum"/>
 
305
        <spatial_discretisation>
 
306
          <discontinuous_galerkin>
 
307
            <viscosity_scheme>
 
308
              <compact_discontinuous_galerkin/>
 
309
              <tensor_form/>
 
310
            </viscosity_scheme>
 
311
            <advection_scheme>
 
312
              <upwind/>
 
313
              <integrate_advection_by_parts>
 
314
                <twice/>
 
315
              </integrate_advection_by_parts>
 
316
            </advection_scheme>
 
317
          </discontinuous_galerkin>
 
318
          <conservative_advection>
 
319
            <real_value rank="0">0</real_value>
 
320
          </conservative_advection>
 
321
        </spatial_discretisation>
 
322
        <temporal_discretisation>
 
323
          <theta>
 
324
            <real_value rank="0">1.0</real_value>
 
325
          </theta>
 
326
          <relaxation>
 
327
            <real_value rank="0">0.5</real_value>
 
328
          </relaxation>
 
329
        </temporal_discretisation>
 
330
        <solver>
 
331
          <iterative_method name="gmres">
 
332
            <restart>
 
333
              <integer_value rank="0">30</integer_value>
 
334
            </restart>
 
335
          </iterative_method>
 
336
          <preconditioner name="sor"/>
 
337
          <relative_error>
 
338
            <real_value rank="0">1.0e-7</real_value>
 
339
          </relative_error>
 
340
          <max_iterations>
 
341
            <integer_value rank="0">1000</integer_value>
 
342
          </max_iterations>
 
343
          <never_ignore_solver_failures/>
 
344
          <diagnostics>
 
345
            <monitors/>
 
346
          </diagnostics>
 
347
        </solver>
 
348
        <initial_condition name="WholeMesh">
 
349
          <constant>
 
350
            <real_value shape="1" dim1="dim" rank="1">0</real_value>
 
351
          </constant>
 
352
        </initial_condition>
 
353
        <boundary_conditions name="Left">
 
354
          <surface_ids>
 
355
            <integer_value shape="1" rank="1">1</integer_value>
 
356
          </surface_ids>
 
357
          <type name="dirichlet">
 
358
            <align_bc_with_cartesian>
 
359
              <x_component>
 
360
                <constant>
 
361
                  <real_value rank="0">300</real_value>
 
362
                </constant>
 
363
              </x_component>
 
364
            </align_bc_with_cartesian>
 
365
          </type>
 
366
        </boundary_conditions>
 
367
        <tensor_field name="Viscosity" rank="2">
 
368
          <prescribed>
 
369
            <value name="WholeMesh">
 
370
              <anisotropic_asymmetric>
 
371
                <constant>
 
372
                  <real_value symmetric="false" dim2="dim" shape="1 1" dim1="dim" rank="2">1.85e-5</real_value>
 
373
                </constant>
 
374
              </anisotropic_asymmetric>
 
375
            </value>
 
376
            <output/>
 
377
          </prescribed>
 
378
        </tensor_field>
 
379
        <output/>
 
380
        <stat>
 
381
          <include_in_stat/>
 
382
          <previous_time_step>
 
383
            <exclude_from_stat/>
 
384
          </previous_time_step>
 
385
          <nonlinear_field>
 
386
            <exclude_from_stat/>
 
387
          </nonlinear_field>
 
388
        </stat>
 
389
        <convergence>
 
390
          <include_in_convergence/>
 
391
        </convergence>
 
392
        <detectors>
 
393
          <include_in_detectors/>
 
394
        </detectors>
 
395
        <steady_state>
 
396
          <include_in_steady_state/>
 
397
        </steady_state>
 
398
        <consistent_interpolation/>
 
399
      </prognostic>
 
400
    </vector_field>
 
401
    <scalar_field name="CFLNumber" rank="0">
 
402
      <diagnostic>
 
403
        <algorithm name="Internal" material_phase_support="multiple"/>
 
404
        <mesh name="VelocityMesh"/>
 
405
        <output/>
 
406
        <stat/>
 
407
        <convergence>
 
408
          <include_in_convergence/>
 
409
        </convergence>
 
410
        <detectors>
 
411
          <include_in_detectors/>
 
412
        </detectors>
 
413
        <steady_state>
 
414
          <include_in_steady_state/>
 
415
        </steady_state>
 
416
      </diagnostic>
 
417
    </scalar_field>
 
418
    <scalar_field name="InternalEnergy" rank="0">
 
419
      <prognostic>
 
420
        <mesh name="PressureMesh"/>
 
421
        <equation name="InternalEnergy">
 
422
          <density name="Density"/>
 
423
        </equation>
 
424
        <spatial_discretisation>
 
425
          <control_volumes>
 
426
            <face_value name="FirstOrderUpwind"/>
 
427
            <diffusion_scheme name="BassiRebay"/>
 
428
          </control_volumes>
 
429
          <conservative_advection>
 
430
            <real_value rank="0">0.0</real_value>
 
431
          </conservative_advection>
 
432
        </spatial_discretisation>
 
433
        <temporal_discretisation>
 
434
          <theta>
 
435
            <real_value rank="0">1.0</real_value>
 
436
          </theta>
 
437
        </temporal_discretisation>
 
438
        <solver>
 
439
          <iterative_method name="gmres">
 
440
            <restart>
 
441
              <integer_value rank="0">30</integer_value>
 
442
            </restart>
 
443
          </iterative_method>
 
444
          <preconditioner name="sor"/>
 
445
          <relative_error>
 
446
            <real_value rank="0">1.0e-7</real_value>
 
447
          </relative_error>
 
448
          <max_iterations>
 
449
            <integer_value rank="0">1000</integer_value>
 
450
          </max_iterations>
 
451
          <never_ignore_solver_failures/>
 
452
          <diagnostics>
 
453
            <monitors/>
 
454
          </diagnostics>
 
455
        </solver>
 
456
        <initial_condition name="WholeMesh">
 
457
          <constant>
 
458
            <real_value rank="0">405138.9</real_value>
 
459
          </constant>
 
460
        </initial_condition>
 
461
        <boundary_conditions name="Left">
 
462
          <surface_ids>
 
463
            <integer_value shape="1" rank="1">1</integer_value>
 
464
          </surface_ids>
 
465
          <type name="dirichlet">
 
466
            <constant>
 
467
              <real_value rank="0">1687200</real_value>
 
468
            </constant>
 
469
          </type>
 
470
        </boundary_conditions>
 
471
        <output/>
 
472
        <stat/>
 
473
        <convergence>
 
474
          <include_in_convergence/>
 
475
        </convergence>
 
476
        <detectors>
 
477
          <include_in_detectors/>
 
478
        </detectors>
 
479
        <steady_state>
 
480
          <include_in_steady_state/>
 
481
        </steady_state>
 
482
        <consistent_interpolation/>
 
483
      </prognostic>
 
484
    </scalar_field>
 
485
    <scalar_field name="HydrostaticReferencePressure" rank="0">
 
486
      <prescribed>
 
487
        <mesh name="PressureMesh"/>
 
488
        <value name="WholeMesh">
 
489
          <python>
 
490
            <string_value lines="20" type="code" language="python">def val(X,t):
 
491
  from math import exp
 
492
  T0 = 288.15 # K
 
493
  R = 287.058 # Specific gas constant in J/Kg
 
494
  g = 9.8 # m/s**2
 
495
  h = X[0] # m
 
496
  p0 = 101325
 
497
  e = T0*1406.0
 
498
  gamma = 1.33
 
499
  p = p0*exp(-(g*h)/((gamma-1.0)*e))
 
500
  return p</string_value>
 
501
          </python>
 
502
        </value>
 
503
        <output/>
 
504
        <stat/>
 
505
        <detectors>
 
506
          <exclude_from_detectors/>
 
507
        </detectors>
 
508
      </prescribed>
 
509
    </scalar_field>
 
510
    <scalar_field name="HydrostaticReferenceDensity" rank="0">
 
511
      <prescribed>
 
512
        <mesh name="PressureMesh"/>
 
513
        <value name="WholeMesh">
 
514
          <python>
 
515
            <string_value lines="20" type="code" language="python">def val(X,t):
 
516
  from math import exp
 
517
  T0 = 288.15 # K
 
518
  R = 287.058 # Specific gas constant in J/Kg
 
519
  g = 9.8 # m/s**2
 
520
  h = X[0] # m
 
521
  p0 = 101325
 
522
  e = T0*1406.0
 
523
  gamma = 1.33
 
524
  p = p0*exp(-(g*h)/((gamma-1.0)*e))
 
525
  rho = p/((gamma-1.0)*e)
 
526
  return rho</string_value>
 
527
          </python>
 
528
        </value>
 
529
        <output/>
 
530
        <stat/>
 
531
        <detectors>
 
532
          <exclude_from_detectors/>
 
533
        </detectors>
 
534
      </prescribed>
 
535
    </scalar_field>
 
536
  </material_phase>
 
537
</fluidity_options>