« back to all changes in this revision
Viewing changes to examples/water_collapse/water_collapse.flml
- browse files at revision 1597
- compare with another revision
- download diff
- download tarball
- view history from revision 1597
- Committer: cwilson
- Date: 2010-07-23 15:06:43 UTC
- Revision ID: svn-v4:5bf5533e-7014-46e3-b1bb-cce4b9d03719:trunk:1622
Setting up an examples directory with the beginnings of the water collapse example in it for Gareth to play with. Have set this up like a test (with an xml file and a Makefile) but none of that bit of it runs properly yet as its just been copied from the water_collapse_2d test case and then had the mesh and input changed.
- files added:
- examples
- examples/water_collapse
- examples/water_collapse/pictures
added
removed
examples/water_collapse/water_collapse.flml
1
<?xml version='1.0' encoding='utf-8'?>
2
<fluidity_options>
3
<simulation_name>
4
<string_value lines="1">water_collapse</string_value>
5
</simulation_name>
6
<problem_type>
7
<string_value lines="1">multimaterial</string_value>
8
</problem_type>
9
<geometry>
10
<dimension replaces="NDIM">
11
<integer_value rank="0">2</integer_value>
12
</dimension>
13
<mesh name="CoordinateMesh">
14
<from_file file_name="2dtank">
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
<quadrature>
38
<degree replaces="NGI">
39
<integer_value rank="0">2</integer_value>
40
</degree>
41
<controlvolume_surface_degree>
42
<integer_value rank="0">1</integer_value>
43
</controlvolume_surface_degree>
44
</quadrature>
45
</geometry>
46
<io>
47
<dump_format>
48
<string_value>vtk</string_value>
49
</dump_format>
50
<dump_period>
51
<constant replaces="TIMDUM">
52
<real_value rank="0">0.05</real_value>
53
<comment>Should give approximately 50 dumps</comment>
54
</constant>
55
</dump_period>
56
<output_mesh name="CoordinateMesh"/>
57
<checkpointing>
58
<checkpoint_period_in_dumps>
59
<integer_value rank="0">50</integer_value>
60
</checkpoint_period_in_dumps>
61
<checkpoint_at_end/>
62
</checkpointing>
63
<stat/>
64
<detectors>
65
<static_detector name="P2">
66
<location>
67
<real_value shape="2" dim1="dim" rank="1">3.22 0.16</real_value>
68
</location>
69
</static_detector>
70
<static_detector name="P3">
71
<location>
72
<real_value shape="2" dim1="dim" rank="1">3.22 0.584</real_value>
73
</location>
74
</static_detector>
75
<static_detector name="P4">
76
<location>
77
<real_value shape="2" dim1="dim" rank="1">3.22 1.0</real_value>
78
</location>
79
</static_detector>
80
</detectors>
81
</io>
82
<timestepping>
83
<current_time replaces="ACCTIM">
84
<real_value rank="0">0.0</real_value>
85
</current_time>
86
<timestep replaces="DT">
87
<real_value rank="0">0.0002</real_value>
88
</timestep>
89
<finish_time replaces="LTIME">
90
<real_value rank="0">2.5</real_value>
91
</finish_time>
92
<nonlinear_iterations replaces="ITINOI">
93
<integer_value rank="0">2</integer_value>
94
</nonlinear_iterations>
95
<adaptive_timestep replaces="TIMAUT">
96
<requested_cfl>
97
<real_value rank="0">2.5</real_value>
98
</requested_cfl>
99
<courant_number name="ControlVolumeCFLNumber">
100
<mesh name="CoordinateMesh"/>
101
</courant_number>
102
<increase_tolerance>
103
<real_value rank="0">1.1</real_value>
104
</increase_tolerance>
105
</adaptive_timestep>
106
</timestepping>
107
<physical_parameters>
108
<gravity>
109
<magnitude>
110
<real_value rank="0">9.81</real_value>
111
</magnitude>
112
<vector_field name="GravityDirection" rank="1">
113
<prescribed>
114
<mesh name="CoordinateMesh"/>
115
<value name="WholeMesh">
116
<constant>
117
<real_value shape="2" dim1="dim" rank="1">0 -1</real_value>
118
</constant>
119
</value>
120
<output>
121
<exclude_from_vtu/>
122
</output>
123
<stat>
124
<exclude_from_stat/>
125
</stat>
126
<detectors>
127
<exclude_from_detectors/>
128
</detectors>
129
</prescribed>
130
</vector_field>
131
</gravity>
132
</physical_parameters>
133
<material_phase name="Air">
134
<equation_of_state>
135
<fluids>
136
<linear replaces="EQNSTA=0 or 1">
137
<reference_density replaces="DENINI">
138
<real_value rank="0">1.0</real_value>
139
</reference_density>
140
</linear>
141
</fluids>
142
</equation_of_state>
143
<vector_field name="Velocity" rank="1">
144
<aliased material_phase_name="Water" field_name="Velocity"/>
145
</vector_field>
146
<scalar_field name="MaterialVolumeFraction" rank="0">
147
<diagnostic>
148
<algorithm name="Internal" material_phase_support="multiple"/>
149
<mesh name="CoordinateMesh"/>
150
<output>
151
<exclude_from_vtu/>
152
</output>
153
<stat>
154
<exclude_from_stat/>
155
</stat>
156
<convergence>
157
<exclude_from_convergence/>
158
</convergence>
159
<detectors>
160
<exclude_from_detectors/>
161
</detectors>
162
<steady_state>
163
<exclude_from_steady_state/>
164
</steady_state>
165
</diagnostic>
166
</scalar_field>
167
</material_phase>
168
<material_phase name="Water">
169
<equation_of_state>
170
<fluids>
171
<linear replaces="EQNSTA=0 or 1">
172
<reference_density replaces="DENINI">
173
<real_value rank="0">1000.0</real_value>
174
</reference_density>
175
</linear>
176
</fluids>
177
</equation_of_state>
178
<scalar_field name="Pressure" rank="0">
179
<prognostic>
180
<mesh name="CoordinateMesh"/>
181
<spatial_discretisation>
182
<control_volumes replaces="NCOLOP = X1XXXXX -> INTERF, NOFILT"/>
183
</spatial_discretisation>
184
<scheme>
185
<poisson_pressure_solution replaces="POISON">
186
<string_value lines="1">never</string_value>
187
</poisson_pressure_solution>
188
<use_projection_method replaces="PROJEC"/>
189
</scheme>
190
<solver>
191
<iterative_method name="cg"/>
192
<preconditioner name="sor"/>
193
<relative_error>
194
<real_value rank="0">1e-10</real_value>
195
</relative_error>
196
<max_iterations>
197
<integer_value rank="0">10000</integer_value>
198
</max_iterations>
199
<never_ignore_solver_failures/>
200
<diagnostics>
201
<print_norms/>
202
<monitors/>
203
</diagnostics>
204
</solver>
205
<output/>
206
<stat/>
207
<convergence>
208
<exclude_from_convergence/>
209
</convergence>
210
<detectors>
211
<include_in_detectors/>
212
</detectors>
213
<steady_state>
214
<exclude_from_steady_state/>
215
</steady_state>
216
<consistent_interpolation/>
217
<enforce_discrete_properties>
218
<solenoidal_lagrange_update/>
219
</enforce_discrete_properties>
220
</prognostic>
221
</scalar_field>
222
<scalar_field name="Density" rank="0">
223
<diagnostic>
224
<algorithm name="Internal" material_phase_support="multiple"/>
225
<mesh name="CoordinateMesh"/>
226
<output/>
227
<stat>
228
<exclude_from_stat/>
229
</stat>
230
<convergence>
231
<exclude_from_convergence/>
232
</convergence>
233
<detectors>
234
<exclude_from_detectors/>
235
</detectors>
236
<steady_state>
237
<exclude_from_steady_state/>
238
</steady_state>
239
</diagnostic>
240
</scalar_field>
241
<vector_field name="Velocity" rank="1">
242
<prognostic>
243
<mesh name="VelocityMesh"/>
244
<equation name="LinearMomentum"/>
245
<spatial_discretisation>
246
<discontinuous_galerkin replaces="DISOPT">
247
<mass_terms>
248
<lump_mass_matrix/>
249
</mass_terms>
250
<viscosity_scheme>
251
<bassi_rebay/>
252
</viscosity_scheme>
253
<advection_scheme>
254
<upwind/>
255
<integrate_advection_by_parts>
256
<once/>
257
</integrate_advection_by_parts>
258
<integrate_conservation_term_by_parts/>
259
</advection_scheme>
260
</discontinuous_galerkin>
261
<conservative_advection replaces="BETA">
262
<real_value rank="0">0</real_value>
263
</conservative_advection>
264
</spatial_discretisation>
265
<temporal_discretisation>
266
<theta replaces="THETA">
267
<real_value rank="0">0.5</real_value>
268
</theta>
269
<relaxation replaces="ITHETA">
270
<real_value rank="0">0.5</real_value>
271
</relaxation>
272
</temporal_discretisation>
273
<solver>
274
<iterative_method name="gmres">
275
<restart>
276
<integer_value rank="0">30</integer_value>
277
</restart>
278
</iterative_method>
279
<preconditioner name="eisenstat"/>
280
<relative_error>
281
<real_value rank="0">1e-10</real_value>
282
</relative_error>
283
<max_iterations>
284
<integer_value rank="0">10000</integer_value>
285
</max_iterations>
286
<never_ignore_solver_failures/>
287
<diagnostics>
288
<print_norms/>
289
<monitors/>
290
</diagnostics>
291
</solver>
292
<initial_condition name="WholeMesh">
293
<constant>
294
<real_value shape="2" dim1="dim" rank="1">0.0 0.0</real_value>
295
</constant>
296
</initial_condition>
297
<boundary_conditions replaces="boundary, TMPER1 TMPER2 TMPERI" name="bottom">
298
<surface_ids>
299
<integer_value shape="1" rank="1">19</integer_value>
300
</surface_ids>
301
<type name="dirichlet">
302
<apply_weakly/>
303
<align_bc_with_cartesian>
304
<x_component>
305
<constant>
306
<real_value rank="0">0</real_value>
307
</constant>
308
</x_component>
309
<y_component>
310
<constant>
311
<real_value rank="0">0</real_value>
312
</constant>
313
</y_component>
314
</align_bc_with_cartesian>
315
</type>
316
</boundary_conditions>
317
<boundary_conditions replaces="boundary, TMPER1 TMPER2 TMPERI" name="leftside">
318
<surface_ids>
319
<integer_value shape="1" rank="1">22</integer_value>
320
</surface_ids>
321
<type name="dirichlet">
322
<apply_weakly/>
323
<align_bc_with_cartesian>
324
<x_component>
325
<constant>
326
<real_value rank="0">0</real_value>
327
</constant>
328
</x_component>
329
<y_component>
330
<constant>
331
<real_value rank="0">0</real_value>
332
</constant>
333
</y_component>
334
</align_bc_with_cartesian>
335
</type>
336
</boundary_conditions>
337
<boundary_conditions replaces="boundary, TMPER1 TMPER2 TMPERI" name="rightside">
338
<surface_ids>
339
<integer_value shape="1" rank="1">20</integer_value>
340
</surface_ids>
341
<type name="dirichlet">
342
<apply_weakly/>
343
<align_bc_with_cartesian>
344
<x_component>
345
<constant>
346
<real_value rank="0">0</real_value>
347
</constant>
348
</x_component>
349
<y_component>
350
<constant>
351
<real_value rank="0">0</real_value>
352
</constant>
353
</y_component>
354
</align_bc_with_cartesian>
355
</type>
356
</boundary_conditions>
357
<output/>
358
<stat>
359
<include_in_stat/>
360
<previous_time_step>
361
<exclude_from_stat/>
362
</previous_time_step>
363
<nonlinear_field>
364
<exclude_from_stat/>
365
</nonlinear_field>
366
</stat>
367
<convergence>
368
<exclude_from_convergence/>
369
</convergence>
370
<detectors>
371
<exclude_from_detectors/>
372
</detectors>
373
<steady_state>
374
<exclude_from_steady_state/>
375
</steady_state>
376
<galerkin_projection>
377
<discontinuous/>
378
</galerkin_projection>
379
<enforce_discrete_properties>
380
<solenoidal>
381
<interpolated_field>
382
<discontinuous>
383
<lump_mass_matrix/>
384
</discontinuous>
385
</interpolated_field>
386
<lagrange_multiplier>
387
<mesh name="CoordinateMesh"/>
388
<spatial_discretisation>
389
<control_volumes/>
390
</spatial_discretisation>
391
<update_scalar_field name="Pressure"/>
392
<solver>
393
<iterative_method name="cg"/>
394
<preconditioner name="sor"/>
395
<relative_error>
396
<real_value rank="0">1e-10</real_value>
397
</relative_error>
398
<max_iterations>
399
<integer_value rank="0">10000</integer_value>
400
</max_iterations>
401
<never_ignore_solver_failures/>
402
<diagnostics>
403
<monitors/>
404
</diagnostics>
405
</solver>
406
</lagrange_multiplier>
407
</solenoidal>
408
</enforce_discrete_properties>
409
</prognostic>
410
</vector_field>
411
<scalar_field name="ControlVolumeCFLNumber" rank="0">
412
<diagnostic>
413
<algorithm name="Internal" material_phase_support="multiple"/>
414
<mesh name="CoordinateMesh"/>
415
<output/>
416
<stat/>
417
<convergence>
418
<exclude_from_convergence/>
419
</convergence>
420
<detectors>
421
<exclude_from_detectors/>
422
</detectors>
423
<steady_state>
424
<exclude_from_steady_state/>
425
</steady_state>
426
</diagnostic>
427
</scalar_field>
428
<scalar_field name="ControlVolumeDivergence" rank="0">
429
<diagnostic field_name="Velocity">
430
<algorithm name="Internal" material_phase_support="multiple"/>
431
<mesh name="CoordinateMesh"/>
432
<output/>
433
<stat/>
434
<convergence>
435
<exclude_from_convergence/>
436
</convergence>
437
<detectors>
438
<exclude_from_detectors/>
439
</detectors>
440
<steady_state>
441
<exclude_from_steady_state/>
442
</steady_state>
443
</diagnostic>
444
</scalar_field>
445
<scalar_field name="MaterialVolumeFraction" rank="0">
446
<prognostic>
447
<mesh name="CoordinateMesh"/>
448
<equation name="AdvectionDiffusion"/>
449
<spatial_discretisation>
450
<coupled_cv>
451
<face_value replaces="MOD(INT(ABS(NDISOT)/1),10) = 8,9" name="HyperC">
452
<project_upwind_value_from_point replaces="NDISOT > 0 on simplex meshes">
453
<reflect_off_domain_boundaries/>
454
<bound_projected_value_locally/>
455
<store_upwind_elements replaces="ABS(NDISOT) >= 1000"/>
456
</project_upwind_value_from_point>
457
<courant_number name="ControlVolumeCFLNumber"/>
458
</face_value>
459
<parent_sum>
460
<target_maximum>
461
<real_value rank="0">1</real_value>
462
</target_maximum>
463
<target_minimum>
464
<real_value rank="0">0</real_value>
465
</target_minimum>
466
</parent_sum>
467
</coupled_cv>
468
<conservative_advection replaces="TBETA, DEFALT (TBETA = 0.0)">
469
<real_value rank="0">0</real_value>
470
</conservative_advection>
471
</spatial_discretisation>
472
<temporal_discretisation>
473
<theta replaces="TTHETA, DEFALT (TTHETA = 0.5)">
474
<real_value rank="0">0</real_value>
475
</theta>
476
<control_volumes>
477
<maximum_courant_number_per_subcycle>
478
<real_value rank="0">0.25</real_value>
479
<courant_number name="ControlVolumeCFLNumber"/>
480
</maximum_courant_number_per_subcycle>
481
<pivot_theta>
482
<real_value rank="0">0</real_value>
483
</pivot_theta>
484
</control_volumes>
485
</temporal_discretisation>
486
<explicit/>
487
<initial_condition name="WholeMesh">
488
<python>
489
<string_value lines="20" type="python">def val(X,t):
490
if((X[0]<=1.2) and (X[1]<=0.6)):
491
return 1.0
492
else:
493
return 0.0</string_value>
494
</python>
495
</initial_condition>
496
<output/>
497
<stat>
498
<include_mixing_stats replaces="heaviside.dat file" name="Mixing">
499
<control_volumes/>
500
<mixing_bin_bounds>
501
<real_value shape="9" rank="1">-1 -1e-08 1e-08 0.025 0.5 0.975 0.99999999 1.00000001 2</real_value>
502
</mixing_bin_bounds>
503
<tolerance>
504
<real_value rank="0">0</real_value>
505
</tolerance>
506
</include_mixing_stats>
507
</stat>
508
<convergence>
509
<exclude_from_convergence/>
510
</convergence>
511
<detectors>
512
<exclude_from_detectors/>
513
</detectors>
514
<steady_state>
515
<exclude_from_steady_state/>
516
</steady_state>
517
<adaptivity_options>
518
<absolute_measure replaces="ADOPTT = 0">
519
<scalar_field replaces="ADWEIT" name="InterpolationErrorBound" rank="0">
520
<prescribed>
521
<value name="WholeMesh">
522
<constant>
523
<real_value rank="0">0.075</real_value>
524
</constant>
525
</value>
526
<output>
527
<exclude_from_vtu/>
528
</output>
529
<stat>
530
<exclude_from_stat/>
531
</stat>
532
<detectors>
533
<exclude_from_detectors/>
534
</detectors>
535
</prescribed>
536
</scalar_field>
537
</absolute_measure>
538
</adaptivity_options>
539
<galerkin_projection>
540
<continuous>
541
<bounded name="Diffuse">
542
<boundedness_iterations>
543
<integer_value rank="0">10000</integer_value>
544
<tolerance>
545
<real_value rank="0">1.e-6</real_value>
546
</tolerance>
547
</boundedness_iterations>
548
<bounds>
549
<upper_bound>
550
<real_value rank="0">1.0</real_value>
551
<coupled/>
552
</upper_bound>
553
<lower_bound>
554
<real_value rank="0">0.0</real_value>
555
</lower_bound>
556
</bounds>
557
<repair_deviations>
558
<tolerance>
559
<real_value rank="0">1.e-10</real_value>
560
</tolerance>
561
</repair_deviations>
562
</bounded>
563
<solver>
564
<iterative_method name="cg"/>
565
<preconditioner name="sor"/>
566
<relative_error>
567
<real_value rank="0">1e-10</real_value>
568
</relative_error>
569
<max_iterations>
570
<integer_value rank="0">10000</integer_value>
571
</max_iterations>
572
<never_ignore_solver_failures/>
573
<diagnostics>
574
<monitors/>
575
</diagnostics>
576
</solver>
577
</continuous>
578
</galerkin_projection>
579
<priority>
580
<integer_value rank="0">2</integer_value>
581
</priority>
582
</prognostic>
583
</scalar_field>
584
</material_phase>
585
<mesh_adaptivity>
586
<hr_adaptivity replaces="ADMESH">
587
<period_in_timesteps>
588
<integer_value rank="0">10</integer_value>
589
</period_in_timesteps>
590
<maximum_number_of_nodes replaces="MXNODS">
591
<integer_value rank="0">19097</integer_value>
592
</maximum_number_of_nodes>
593
<anisotropic_gradation>
594
<tensor_field name="Gamma">
595
<anisotropic_symmetric>
596
<constant>
597
<real_value symmetric="true" dim2="dim" shape="2 2" dim1="dim" rank="2">1 0 0 1</real_value>
598
</constant>
599
</anisotropic_symmetric>
600
</tensor_field>
601
</anisotropic_gradation>
602
<metric_advection>
603
<spatial_discretisation>
604
<control_volumes>
605
<face_value replaces="MOD(INT(ABS(NDISOT)/1),10) = 0,1" name="FirstOrderUpwind"/>
606
</control_volumes>
607
<conservative_advection replaces="TBETA, DEFALT (TBETA = 0.0)">
608
<real_value rank="0">0</real_value>
609
</conservative_advection>
610
</spatial_discretisation>
611
<temporal_discretisation>
612
<theta replaces="TTHETA, DEFALT (TTHETA = 0.5)">
613
<real_value rank="0">0.5</real_value>
614
</theta>
615
<maximum_courant_number_per_subcycle>
616
<real_value rank="0">2.5</real_value>
617
<courant_number name="ControlVolumeCFLNumber"/>
618
</maximum_courant_number_per_subcycle>
619
<scale_advection_time>
620
<real_value rank="0">3</real_value>
621
</scale_advection_time>
622
<control_volumes>
623
<pivot_theta>
624
<real_value rank="0">0.5</real_value>
625
</pivot_theta>
626
</control_volumes>
627
</temporal_discretisation>
628
<solver>
629
<iterative_method name="gmres">
630
<restart>
631
<integer_value rank="0">30</integer_value>
632
</restart>
633
</iterative_method>
634
<preconditioner name="sor"/>
635
<relative_error>
636
<real_value rank="0">1e-07</real_value>
637
</relative_error>
638
<max_iterations>
639
<integer_value rank="0">10000</integer_value>
640
</max_iterations>
641
<never_ignore_solver_failures/>
642
<diagnostics>
643
<monitors/>
644
</diagnostics>
645
</solver>
646
<output/>
647
</metric_advection>
648
<tensor_field name="MinimumEdgeLengths">
649
<anisotropic_symmetric>
650
<constant>
651
<real_value symmetric="true" dim2="dim" shape="2 2" dim1="dim" rank="2">0.0033333333 0 0 0.0033333333</real_value>
652
</constant>
653
</anisotropic_symmetric>
654
</tensor_field>
655
<tensor_field name="MaximumEdgeLengths">
656
<anisotropic_symmetric>
657
<constant>
658
<real_value symmetric="true" dim2="dim" shape="2 2" dim1="dim" rank="2">1.61 0 0 1</real_value>
659
</constant>
660
</anisotropic_symmetric>
661
</tensor_field>
662
<adapt_at_first_timestep>
663
<number_of_adapts>
664
<integer_value rank="0">10</integer_value>
665
</number_of_adapts>
666
<output_adapted_mesh/>
667
</adapt_at_first_timestep>
668
</hr_adaptivity>
669
</mesh_adaptivity>
670
</fluidity_options>
Loggerhead is a web-based interface for Breezy
Version: 2.0.1