~fluidity-core/fluidity/sea-ice-branch
» Revision
3645
: tests/3material-droplet-p1dgp2-test-cty-cv-rhop1dg/3material-droplet-p1dgp2-test-cty-cv-rhop1dg.flml
« back to all changes in this revision
Viewing changes to tests/3material-droplet-p1dgp2-test-cty-cv-rhop1dg/3material-droplet-p1dgp2-test-cty-cv-rhop1dg.flml
- browse files at revision 3645
- compare with another revision
- download diff
- download tarball
- view history from revision 3645
- Committer: Simon Mouradian
- Date: 2012-05-30 09:08:23 UTC
- mfrom: (3520.32.274 fluidity)
- Revision ID: simon.mouradian06@imperial.ac.uk-20120530090823-pq0fgemjpqn72j8d
merge from lp:fluidity
- files added:
- manual/fluidity_logo.svg
- tests/3material-droplet-p1dgp2-test-cty-cv-rhop1dg
- tests/3material-droplet-p1dgp2-test-cty-cv-rhop1dg-3d
- tests/3material-droplet-p1dgp2-test-cty-cv-rhop1dg/3material-droplet-p1dgp2-test-cty-cv-rhop1dg.flml
- tests/column-parallel
- tests/column-parallel/src
- files modified:
- assemble/Diagnostic_Fields_Matrices.F90
added
removed
tests/3material-droplet-p1dgp2-test-cty-cv-rhop1dg/3material-droplet-p1dgp2-test-cty-cv-rhop1dg.flml
1
<?xml version='1.0' encoding='utf-8'?>
2
<fluidity_options>
3
<simulation_name>
4
<string_value lines="1">3material-droplet-p1dgp2-test-cty-cv-rhop1dg</string_value>
5
</simulation_name>
6
<problem_type>
7
<string_value lines="1">multimaterial</string_value>
8
</problem_type>
9
<geometry>
10
<dimension>
11
<integer_value rank="0">2</integer_value>
12
</dimension>
13
<mesh name="CoordinateMesh">
14
<from_file file_name="square">
15
<format name="gmsh"/>
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
<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">2</integer_value>
43
</polynomial_degree>
44
</mesh_shape>
45
<stat>
46
<exclude_from_stat/>
47
</stat>
48
</from_mesh>
49
</mesh>
50
<mesh name="OutputMesh">
51
<from_mesh>
52
<mesh name="CoordinateMesh"/>
53
<mesh_shape>
54
<polynomial_degree>
55
<integer_value rank="0">2</integer_value>
56
</polynomial_degree>
57
</mesh_shape>
58
<mesh_continuity>
59
<string_value>discontinuous</string_value>
60
</mesh_continuity>
61
<stat>
62
<exclude_from_stat/>
63
</stat>
64
</from_mesh>
65
</mesh>
66
<mesh name="DensityMesh">
67
<from_mesh>
68
<mesh name="CoordinateMesh"/>
69
<mesh_shape>
70
<polynomial_degree>
71
<integer_value rank="0">1</integer_value>
72
</polynomial_degree>
73
</mesh_shape>
74
<mesh_continuity>
75
<string_value>discontinuous</string_value>
76
</mesh_continuity>
77
<stat>
78
<exclude_from_stat/>
79
</stat>
80
</from_mesh>
81
</mesh>
82
<quadrature>
83
<degree>
84
<integer_value rank="0">5</integer_value>
85
</degree>
86
<controlvolume_surface_degree>
87
<integer_value rank="0">4</integer_value>
88
</controlvolume_surface_degree>
89
</quadrature>
90
</geometry>
91
<io>
92
<dump_format>
93
<string_value>vtk</string_value>
94
</dump_format>
95
<dump_period_in_timesteps>
96
<constant>
97
<integer_value rank="0">10</integer_value>
98
</constant>
99
</dump_period_in_timesteps>
100
<output_mesh name="OutputMesh"/>
101
<stat>
102
<output_at_start/>
103
</stat>
104
</io>
105
<timestepping>
106
<current_time>
107
<real_value rank="0">0.0</real_value>
108
</current_time>
109
<timestep>
110
<real_value rank="0">0.0025</real_value>
111
</timestep>
112
<finish_time>
113
<real_value rank="0">200.0</real_value>
114
<comment>200.0</comment>
115
</finish_time>
116
<final_timestep>
117
<integer_value rank="0">100</integer_value>
118
</final_timestep>
119
<nonlinear_iterations>
120
<integer_value rank="0">2</integer_value>
121
</nonlinear_iterations>
122
<adaptive_timestep>
123
<requested_cfl>
124
<real_value rank="0">1.0</real_value>
125
</requested_cfl>
126
<courant_number name="ControlVolumeCFLNumber">
127
<mesh name="PressureMesh"/>
128
</courant_number>
129
<increase_tolerance>
130
<real_value rank="0">1.05</real_value>
131
</increase_tolerance>
132
</adaptive_timestep>
133
</timestepping>
134
<physical_parameters>
135
<gravity>
136
<magnitude>
137
<real_value rank="0">10.0</real_value>
138
</magnitude>
139
<vector_field name="GravityDirection" rank="1">
140
<prescribed>
141
<mesh name="CoordinateMesh"/>
142
<value name="WholeMesh">
143
<constant>
144
<real_value shape="2" dim1="dim" rank="1">0.0 -1.0</real_value>
145
</constant>
146
</value>
147
<output>
148
<exclude_from_vtu/>
149
</output>
150
<stat>
151
<exclude_from_stat/>
152
</stat>
153
<detectors>
154
<exclude_from_detectors/>
155
</detectors>
156
</prescribed>
157
</vector_field>
158
</gravity>
159
</physical_parameters>
160
<material_phase name="Water">
161
<equation_of_state>
162
<fluids>
163
<linear>
164
<reference_density>
165
<real_value rank="0">1.0e+03</real_value>
166
</reference_density>
167
</linear>
168
</fluids>
169
</equation_of_state>
170
<scalar_field name="Pressure" rank="0">
171
<prognostic>
172
<mesh name="PressureMesh"/>
173
<spatial_discretisation>
174
<continuous_galerkin>
175
<test_continuity_with_cv_dual/>
176
</continuous_galerkin>
177
</spatial_discretisation>
178
<scheme>
179
<poisson_pressure_solution>
180
<string_value lines="1">never</string_value>
181
</poisson_pressure_solution>
182
<use_projection_method/>
183
</scheme>
184
<solver>
185
<iterative_method name="gmres">
186
<restart>
187
<integer_value rank="0">30</integer_value>
188
</restart>
189
</iterative_method>
190
<preconditioner name="eisenstat"/>
191
<relative_error>
192
<real_value rank="0">1.0E-10</real_value>
193
</relative_error>
194
<max_iterations>
195
<integer_value rank="0">5000</integer_value>
196
</max_iterations>
197
<never_ignore_solver_failures/>
198
<diagnostics>
199
<print_norms/>
200
<monitors/>
201
</diagnostics>
202
</solver>
203
<output/>
204
<stat/>
205
<convergence>
206
<include_in_convergence/>
207
</convergence>
208
<detectors>
209
<exclude_from_detectors/>
210
</detectors>
211
<steady_state>
212
<include_in_steady_state/>
213
</steady_state>
214
<consistent_interpolation/>
215
</prognostic>
216
</scalar_field>
217
<scalar_field name="Density" rank="0">
218
<diagnostic>
219
<algorithm name="Internal" material_phase_support="multiple"/>
220
<mesh name="DensityMesh"/>
221
<output/>
222
<stat/>
223
<convergence>
224
<include_in_convergence/>
225
</convergence>
226
<detectors>
227
<include_in_detectors/>
228
</detectors>
229
<steady_state>
230
<include_in_steady_state/>
231
</steady_state>
232
</diagnostic>
233
</scalar_field>
234
<vector_field name="Velocity" rank="1">
235
<prognostic>
236
<mesh name="VelocityMesh"/>
237
<equation name="LinearMomentum"/>
238
<spatial_discretisation>
239
<discontinuous_galerkin>
240
<viscosity_scheme>
241
<compact_discontinuous_galerkin/>
242
</viscosity_scheme>
243
<advection_scheme>
244
<upwind/>
245
<integrate_advection_by_parts>
246
<once/>
247
</integrate_advection_by_parts>
248
<integrate_conservation_term_by_parts/>
249
</advection_scheme>
250
</discontinuous_galerkin>
251
<conservative_advection>
252
<real_value rank="0">0.0</real_value>
253
</conservative_advection>
254
</spatial_discretisation>
255
<temporal_discretisation>
256
<theta>
257
<real_value rank="0">0.5</real_value>
258
</theta>
259
<relaxation>
260
<real_value rank="0">0.5</real_value>
261
</relaxation>
262
</temporal_discretisation>
263
<solver>
264
<iterative_method name="gmres">
265
<restart>
266
<integer_value rank="0">30</integer_value>
267
</restart>
268
</iterative_method>
269
<preconditioner name="eisenstat"/>
270
<relative_error>
271
<real_value rank="0">1.0E-10</real_value>
272
</relative_error>
273
<max_iterations>
274
<integer_value rank="0">5000</integer_value>
275
</max_iterations>
276
<never_ignore_solver_failures/>
277
<diagnostics>
278
<print_norms/>
279
<monitors/>
280
</diagnostics>
281
</solver>
282
<initial_condition name="WholeMesh">
283
<constant>
284
<real_value shape="2" dim1="dim" rank="1">0.0 0.0</real_value>
285
</constant>
286
</initial_condition>
287
<boundary_conditions name="no_flow_out_base">
288
<surface_ids>
289
<integer_value shape="1" rank="1">7</integer_value>
290
</surface_ids>
291
<type name="no_normal_flow"/>
292
</boundary_conditions>
293
<boundary_conditions name="no_flow_out_sides">
294
<surface_ids>
295
<integer_value shape="2" rank="1">8 10</integer_value>
296
</surface_ids>
297
<type name="no_normal_flow"/>
298
</boundary_conditions>
299
<tensor_field name="Viscosity" rank="2">
300
<diagnostic>
301
<mesh name="VelocityMesh"/>
302
<algorithm name="bulk_viscosity" material_phase_support="multiple "/>
303
<output/>
304
<stat>
305
<include_in_stat/>
306
</stat>
307
</diagnostic>
308
</tensor_field>
309
<output/>
310
<stat>
311
<include_in_stat/>
312
<previous_time_step>
313
<exclude_from_stat/>
314
</previous_time_step>
315
<nonlinear_field>
316
<exclude_from_stat/>
317
</nonlinear_field>
318
</stat>
319
<convergence>
320
<include_in_convergence/>
321
</convergence>
322
<detectors>
323
<include_in_detectors/>
324
</detectors>
325
<steady_state>
326
<include_in_steady_state/>
327
</steady_state>
328
<consistent_interpolation/>
329
</prognostic>
330
</vector_field>
331
<scalar_field name="MaterialVolumeFraction" rank="0">
332
<prognostic>
333
<mesh name="PressureMesh"/>
334
<equation name="AdvectionDiffusion"/>
335
<spatial_discretisation>
336
<coupled_cv>
337
<face_value name="HyperC">
338
<project_upwind_value_from_point>
339
<reflect_off_domain_boundaries/>
340
<bound_projected_value_locally/>
341
<store_upwind_elements/>
342
</project_upwind_value_from_point>
343
<courant_number name="ControlVolumeCFLNumber"/>
344
</face_value>
345
<parent_sum>
346
<target_maximum>
347
<real_value rank="0">1.0</real_value>
348
</target_maximum>
349
<target_minimum>
350
<real_value rank="0">0.0</real_value>
351
</target_minimum>
352
</parent_sum>
353
</coupled_cv>
354
<conservative_advection>
355
<real_value rank="0">0.0</real_value>
356
</conservative_advection>
357
</spatial_discretisation>
358
<temporal_discretisation>
359
<theta>
360
<real_value rank="0">0.0</real_value>
361
</theta>
362
<control_volumes>
363
<maximum_courant_number_per_subcycle>
364
<real_value rank="0">0.25</real_value>
365
<courant_number name="ControlVolumeCFLNumber"/>
366
</maximum_courant_number_per_subcycle>
367
<pivot_theta>
368
<real_value rank="0">0.0</real_value>
369
</pivot_theta>
370
</control_volumes>
371
</temporal_discretisation>
372
<explicit/>
373
<initial_condition name="WholeMesh">
374
<python>
375
<string_value lines="20" type="code" language="python">def val(X,t):
376
if (X[1]<=0.0):
377
return 1.0
378
else:
379
return 0.0</string_value>
380
</python>
381
</initial_condition>
382
<boundary_conditions name="incoming_vfrac">
383
<surface_ids>
384
<integer_value shape="1" rank="1">9</integer_value>
385
</surface_ids>
386
<type name="dirichlet">
387
<apply_weakly/>
388
<constant>
389
<real_value rank="0">0.0</real_value>
390
</constant>
391
</type>
392
</boundary_conditions>
393
<output/>
394
<stat>
395
<include_cv_stats/>
396
</stat>
397
<convergence>
398
<include_in_convergence/>
399
</convergence>
400
<detectors>
401
<include_in_detectors/>
402
</detectors>
403
<steady_state>
404
<include_in_steady_state/>
405
</steady_state>
406
<consistent_interpolation/>
407
<priority>
408
<integer_value rank="0">2</integer_value>
409
</priority>
410
</prognostic>
411
</scalar_field>
412
<scalar_field name="ControlVolumeCFLNumber" rank="0">
413
<diagnostic>
414
<algorithm name="Internal" material_phase_support="multiple"/>
415
<mesh name="PressureMesh"/>
416
<output/>
417
<stat/>
418
<convergence>
419
<include_in_convergence/>
420
</convergence>
421
<detectors>
422
<include_in_detectors/>
423
</detectors>
424
<steady_state>
425
<include_in_steady_state/>
426
</steady_state>
427
</diagnostic>
428
</scalar_field>
429
<scalar_field name="SumMaterialVolumeFractions" rank="0">
430
<diagnostic>
431
<algorithm name="Internal" material_phase_support="multiple"/>
432
<mesh name="PressureMesh"/>
433
<output/>
434
<stat/>
435
<convergence>
436
<include_in_convergence/>
437
</convergence>
438
<detectors>
439
<include_in_detectors/>
440
</detectors>
441
<steady_state>
442
<include_in_steady_state/>
443
</steady_state>
444
</diagnostic>
445
</scalar_field>
446
<scalar_field name="ControlVolumeDivergence" rank="0">
447
<diagnostic field_name="Velocity">
448
<algorithm name="Internal" material_phase_support="multiple"/>
449
<mesh name="PressureMesh"/>
450
<output/>
451
<stat/>
452
<convergence>
453
<include_in_convergence/>
454
</convergence>
455
<detectors>
456
<include_in_detectors/>
457
</detectors>
458
<steady_state>
459
<include_in_steady_state/>
460
</steady_state>
461
</diagnostic>
462
</scalar_field>
463
<scalar_field name="DG_CourantNumber" rank="0">
464
<diagnostic>
465
<algorithm name="Internal" material_phase_support="multiple"/>
466
<mesh name="VelocityMesh"/>
467
<output/>
468
<stat/>
469
<convergence>
470
<include_in_convergence/>
471
</convergence>
472
<detectors>
473
<include_in_detectors/>
474
</detectors>
475
<steady_state>
476
<include_in_steady_state/>
477
</steady_state>
478
</diagnostic>
479
</scalar_field>
480
<scalar_field name="GridReynoldsNumber" rank="0">
481
<diagnostic>
482
<algorithm name="Internal" material_phase_support="multiple"/>
483
<mesh name="VelocityMesh"/>
484
<include_density_field/>
485
<output/>
486
<stat/>
487
<convergence>
488
<include_in_convergence/>
489
</convergence>
490
<detectors>
491
<include_in_detectors/>
492
</detectors>
493
<steady_state>
494
<include_in_steady_state/>
495
</steady_state>
496
</diagnostic>
497
</scalar_field>
498
<tensor_field name="MaterialViscosity" rank="2">
499
<prescribed>
500
<mesh name="VelocityMesh"/>
501
<value name="WholeMesh">
502
<isotropic>
503
<constant>
504
<real_value rank="0">1.0</real_value>
505
<comment>should be 9.0e-04</comment>
506
</constant>
507
</isotropic>
508
</value>
509
<output/>
510
</prescribed>
511
</tensor_field>
512
</material_phase>
513
<material_phase name="Droplet">
514
<equation_of_state>
515
<fluids>
516
<linear>
517
<reference_density>
518
<real_value rank="0">1.0e+04</real_value>
519
</reference_density>
520
</linear>
521
</fluids>
522
</equation_of_state>
523
<scalar_field name="Pressure" rank="0">
524
<aliased material_phase_name="Water" field_name="Pressure"/>
525
</scalar_field>
526
<scalar_field name="Density" rank="0">
527
<aliased material_phase_name="Water" field_name="Density"/>
528
</scalar_field>
529
<vector_field name="Velocity" rank="1">
530
<aliased material_phase_name="Water" field_name="Velocity"/>
531
</vector_field>
532
<scalar_field name="MaterialVolumeFraction" rank="0">
533
<prognostic>
534
<mesh name="PressureMesh"/>
535
<equation name="AdvectionDiffusion"/>
536
<spatial_discretisation>
537
<coupled_cv>
538
<face_value name="HyperC">
539
<project_upwind_value_from_point>
540
<reflect_off_domain_boundaries/>
541
<bound_projected_value_locally/>
542
<store_upwind_elements/>
543
</project_upwind_value_from_point>
544
<courant_number name="ControlVolumeCFLNumber"/>
545
</face_value>
546
<parent_sum>
547
<target_maximum>
548
<real_value rank="0">1.0</real_value>
549
</target_maximum>
550
<target_minimum>
551
<real_value rank="0">0.0</real_value>
552
</target_minimum>
553
</parent_sum>
554
</coupled_cv>
555
<conservative_advection>
556
<real_value rank="0">0.0</real_value>
557
</conservative_advection>
558
</spatial_discretisation>
559
<temporal_discretisation>
560
<theta>
561
<real_value rank="0">0.0</real_value>
562
</theta>
563
<control_volumes>
564
<maximum_courant_number_per_subcycle>
565
<real_value rank="0">0.25</real_value>
566
<courant_number name="ControlVolumeCFLNumber"/>
567
</maximum_courant_number_per_subcycle>
568
<pivot_theta>
569
<real_value rank="0">0.0</real_value>
570
</pivot_theta>
571
</control_volumes>
572
</temporal_discretisation>
573
<explicit/>
574
<initial_condition name="WholeMesh">
575
<python>
576
<string_value lines="20" type="code" language="python">def val(X,t):
577
from math import sqrt
578
dx1 = X[0]-0.0
579
dx2 = X[1]-0.7
580
r=sqrt(dx1**2+dx2**2)
581
if (r<0.5):
582
return 1.0
583
else:
584
return 0.0</string_value>
585
</python>
586
</initial_condition>
587
<boundary_conditions name="incoming_vfrac">
588
<surface_ids>
589
<integer_value shape="1" rank="1">9</integer_value>
590
</surface_ids>
591
<type name="dirichlet">
592
<apply_weakly/>
593
<constant>
594
<real_value rank="0">0.0</real_value>
595
</constant>
596
</type>
597
</boundary_conditions>
598
<output/>
599
<stat>
600
<include_cv_stats/>
601
</stat>
602
<convergence>
603
<include_in_convergence/>
604
</convergence>
605
<detectors>
606
<include_in_detectors/>
607
</detectors>
608
<steady_state>
609
<include_in_steady_state/>
610
</steady_state>
611
<consistent_interpolation/>
612
<priority>
613
<integer_value rank="0">3</integer_value>
614
</priority>
615
</prognostic>
616
</scalar_field>
617
<tensor_field name="MaterialViscosity" rank="2">
618
<prescribed>
619
<mesh name="VelocityMesh"/>
620
<value name="WholeMesh">
621
<isotropic>
622
<constant>
623
<real_value rank="0">1.0</real_value>
624
<comment>should be 5.0e-03</comment>
625
</constant>
626
</isotropic>
627
</value>
628
<output/>
629
</prescribed>
630
</tensor_field>
631
</material_phase>
632
<material_phase name="Air">
633
<equation_of_state>
634
<fluids>
635
<linear>
636
<reference_density>
637
<real_value rank="0">1.0</real_value>
638
</reference_density>
639
</linear>
640
</fluids>
641
</equation_of_state>
642
<scalar_field name="Pressure" rank="0">
643
<aliased material_phase_name="Water" field_name="Pressure"/>
644
</scalar_field>
645
<scalar_field name="Density" rank="0">
646
<aliased material_phase_name="Water" field_name="Density"/>
647
</scalar_field>
648
<vector_field name="Velocity" rank="1">
649
<aliased material_phase_name="Water" field_name="Velocity"/>
650
</vector_field>
651
<scalar_field name="MaterialVolumeFraction" rank="0">
652
<diagnostic>
653
<algorithm name="Internal" material_phase_support="multiple"/>
654
<mesh name="PressureMesh"/>
655
<output/>
656
<stat>
657
<include_cv_stats/>
658
</stat>
659
<convergence>
660
<include_in_convergence/>
661
</convergence>
662
<detectors>
663
<include_in_detectors/>
664
</detectors>
665
<steady_state>
666
<include_in_steady_state/>
667
</steady_state>
668
</diagnostic>
669
</scalar_field>
670
<tensor_field name="MaterialViscosity" rank="2">
671
<prescribed>
672
<mesh name="VelocityMesh"/>
673
<value name="WholeMesh">
674
<isotropic>
675
<constant>
676
<real_value rank="0">1.0</real_value>
677
<comment>should be 2.0e-05</comment>
678
</constant>
679
</isotropic>
680
</value>
681
<output/>
682
</prescribed>
683
</tensor_field>
684
</material_phase>
685
</fluidity_options>
Loggerhead is a web-based interface for Breezy
Version: 2.0.1