~fluidity-core/fluidity/MultiFluids_Fast2
» Revision
4163
: legacy_reservoir_prototype/tests/Compositional_Model/CM_2D_2P2C_P1DGP1_10layers.mpml
« back to all changes in this revision
Viewing changes to legacy_reservoir_prototype/tests/Compositional_Model/CM_2D_2P2C_P1DGP1_10layers.mpml
- browse files at revision 4163
- compare with another revision
- download diff
- download tarball
- view history from revision 4163
- Committer: James Robert Percival
- Date: 2013-06-25 20:08:54 UTC
- mfrom: (4162.1.3 MultiFluids_Dev)
- Revision ID: j.percival@imperial.ac.uk-20130625200854-o2q4kqby28tk87pa
Merge back trunk
- files added:
- legacy_reservoir_prototype/tests/Compositional_Model
- legacy_reservoir_prototype/tests/Compositional_Model/src
- files modified:
- legacy_reservoir_prototype/src/Multiphase_TimeLoop.F90
added
removed
legacy_reservoir_prototype/tests/Compositional_Model/CM_2D_2P2C_P1DGP1_10layers.mpml
1
<?xml version='1.0' encoding='utf-8'?>
2
<multiphase_options>
3
<simulation_name>
4
<string_value lines="1">CM_2D_2P2C_P1DGP1_10layers</string_value>
5
</simulation_name>
6
<geometry>
7
<dimension>
8
<integer_value rank="0">2</integer_value>
9
</dimension>
10
<mesh name="CoordinateMesh">
11
<from_file file_name="2d_mesh_10layers">
12
<format name="triangle"/>
13
<stat>
14
<include_in_stat/>
15
</stat>
16
</from_file>
17
</mesh>
18
<mesh name="VelocityMesh">
19
<from_mesh>
20
<mesh name="CoordinateMesh"/>
21
<mesh_shape>
22
<polynomial_degree>
23
<integer_value rank="0">1</integer_value>
24
</polynomial_degree>
25
<element_type>
26
<string_value>overlapping</string_value>
27
</element_type>
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
<mesh_continuity>
46
<string_value>continuous</string_value>
47
</mesh_continuity>
48
<stat>
49
<exclude_from_stat/>
50
</stat>
51
</from_mesh>
52
</mesh>
53
<mesh name="P0DG">
54
<from_mesh>
55
<mesh name="CoordinateMesh"/>
56
<mesh_shape>
57
<polynomial_degree>
58
<integer_value rank="0">0</integer_value>
59
</polynomial_degree>
60
</mesh_shape>
61
<mesh_continuity>
62
<string_value>discontinuous</string_value>
63
</mesh_continuity>
64
<stat>
65
<exclude_from_stat/>
66
</stat>
67
</from_mesh>
68
</mesh>
69
<mesh name="VelocityMesh_Continuous">
70
<from_mesh>
71
<mesh name="CoordinateMesh"/>
72
<mesh_shape>
73
<polynomial_degree>
74
<integer_value rank="0">1</integer_value>
75
</polynomial_degree>
76
</mesh_shape>
77
<mesh_continuity>
78
<string_value>continuous</string_value>
79
</mesh_continuity>
80
<stat>
81
<exclude_from_stat/>
82
</stat>
83
</from_mesh>
84
</mesh>
85
<mesh name="PressureMesh_Continuous">
86
<from_mesh>
87
<mesh name="CoordinateMesh"/>
88
<mesh_shape>
89
<polynomial_degree>
90
<integer_value rank="0">1</integer_value>
91
</polynomial_degree>
92
</mesh_shape>
93
<mesh_continuity>
94
<string_value>continuous</string_value>
95
</mesh_continuity>
96
<stat>
97
<exclude_from_stat/>
98
</stat>
99
</from_mesh>
100
</mesh>
101
<mesh name="PressureMesh_Discontinuous">
102
<from_mesh>
103
<mesh name="CoordinateMesh"/>
104
<mesh_shape>
105
<polynomial_degree>
106
<integer_value rank="0">1</integer_value>
107
</polynomial_degree>
108
</mesh_shape>
109
<mesh_continuity>
110
<string_value>discontinuous</string_value>
111
</mesh_continuity>
112
<stat>
113
<exclude_from_stat/>
114
</stat>
115
</from_mesh>
116
</mesh>
117
<quadrature>
118
<degree>
119
<integer_value rank="0">5</integer_value>
120
</degree>
121
</quadrature>
122
</geometry>
123
<io>
124
<dump_format>
125
<string_value>vtk</string_value>
126
</dump_format>
127
<dump_period_in_timesteps>
128
<constant>
129
<integer_value rank="0">1</integer_value>
130
</constant>
131
</dump_period_in_timesteps>
132
<output_mesh name="P0DG"/>
133
<stat/>
134
</io>
135
<timestepping>
136
<current_time>
137
<real_value rank="0">0.0</real_value>
138
</current_time>
139
<timestep>
140
<real_value rank="0">1e-3</real_value>
141
</timestep>
142
<finish_time>
143
<real_value rank="0">0.5</real_value>
144
</finish_time>
145
<nonlinear_iterations>
146
<integer_value rank="0">3</integer_value>
147
</nonlinear_iterations>
148
</timestepping>
149
<physical_parameters>
150
<mobility>
151
<real_value rank="0">10.0</real_value>
152
</mobility>
153
</physical_parameters>
154
<material_phase name="phase1">
155
<scalar_field name="Pressure" rank="0">
156
<prognostic>
157
<mesh name="PressureMesh"/>
158
<spatial_discretisation>
159
<continuous_galerkin/>
160
</spatial_discretisation>
161
<scheme>
162
<poisson_pressure_solution>
163
<string_value lines="1">never</string_value>
164
</poisson_pressure_solution>
165
<use_projection_method/>
166
</scheme>
167
<solver>
168
<iterative_method name="gmres">
169
<restart>
170
<integer_value rank="0">30</integer_value>
171
</restart>
172
</iterative_method>
173
<preconditioner name="sor"/>
174
<relative_error>
175
<real_value rank="0">1e-8</real_value>
176
</relative_error>
177
<absolute_error>
178
<real_value rank="0">1e-10</real_value>
179
</absolute_error>
180
<max_iterations>
181
<integer_value rank="0">2500</integer_value>
182
</max_iterations>
183
<never_ignore_solver_failures/>
184
<diagnostics>
185
<monitors/>
186
</diagnostics>
187
</solver>
188
<initial_condition name="WholeMesh">
189
<constant>
190
<real_value rank="0">0.0</real_value>
191
</constant>
192
</initial_condition>
193
<boundary_conditions name="outflow">
194
<surface_ids>
195
<integer_value shape="1" rank="1">7</integer_value>
196
</surface_ids>
197
<type name="dirichlet">
198
<apply_weakly/>
199
<constant>
200
<real_value rank="0">0.0</real_value>
201
</constant>
202
</type>
203
</boundary_conditions>
204
<output/>
205
<stat/>
206
<convergence>
207
<include_in_convergence/>
208
</convergence>
209
<detectors>
210
<exclude_from_detectors/>
211
</detectors>
212
<steady_state>
213
<include_in_steady_state/>
214
</steady_state>
215
<no_interpolation/>
216
</prognostic>
217
</scalar_field>
218
<scalar_field name="Density" rank="0">
219
<prognostic>
220
<mesh name="PressureMesh"/>
221
<spatial_discretisation>
222
<control_volumes>
223
<face_value name="FirstOrderUpwind"/>
224
</control_volumes>
225
<conservative_advection>
226
<real_value rank="0">0.0</real_value>
227
</conservative_advection>
228
</spatial_discretisation>
229
<temporal_discretisation>
230
<theta>
231
<real_value rank="0">1.0</real_value>
232
</theta>
233
</temporal_discretisation>
234
<initial_condition name="WholeMesh">
235
<constant>
236
<real_value rank="0">1.0</real_value>
237
</constant>
238
</initial_condition>
239
<output/>
240
<stat/>
241
<convergence>
242
<include_in_convergence/>
243
</convergence>
244
<detectors>
245
<include_in_detectors/>
246
</detectors>
247
<steady_state>
248
<include_in_steady_state/>
249
</steady_state>
250
<consistent_interpolation/>
251
</prognostic>
252
</scalar_field>
253
<vector_field name="Velocity" rank="1">
254
<prognostic>
255
<mesh name="VelocityMesh"/>
256
<equation name="LinearMomentum"/>
257
<spatial_discretisation>
258
<discontinuous_galerkin>
259
<mass_terms/>
260
<viscosity_scheme>
261
<compact_discontinuous_galerkin/>
262
</viscosity_scheme>
263
<advection_scheme>
264
<upwind/>
265
<integrate_advection_by_parts>
266
<twice/>
267
</integrate_advection_by_parts>
268
</advection_scheme>
269
</discontinuous_galerkin>
270
<conservative_advection>
271
<real_value rank="0">1.0</real_value>
272
</conservative_advection>
273
</spatial_discretisation>
274
<temporal_discretisation>
275
<theta>
276
<real_value rank="0">1.0</real_value>
277
</theta>
278
<relaxation>
279
<real_value rank="0">1.0</real_value>
280
</relaxation>
281
</temporal_discretisation>
282
<solver>
283
<iterative_method name="gmres">
284
<restart>
285
<integer_value rank="0">30</integer_value>
286
</restart>
287
</iterative_method>
288
<preconditioner name="sor"/>
289
<relative_error>
290
<real_value rank="0">1e-7</real_value>
291
</relative_error>
292
<max_iterations>
293
<integer_value rank="0">500</integer_value>
294
</max_iterations>
295
<never_ignore_solver_failures/>
296
<diagnostics>
297
<monitors/>
298
</diagnostics>
299
</solver>
300
<initial_condition name="WholeMesh">
301
<constant>
302
<real_value shape="2" dim1="dim" rank="1">0.0 0.0</real_value>
303
</constant>
304
</initial_condition>
305
<boundary_conditions name="inflow">
306
<surface_ids>
307
<integer_value shape="1" rank="1">6</integer_value>
308
</surface_ids>
309
<type name="dirichlet">
310
<apply_weakly/>
311
<align_bc_with_cartesian>
312
<x_component>
313
<constant>
314
<real_value rank="0">1.0</real_value>
315
</constant>
316
</x_component>
317
<y_component>
318
<constant>
319
<real_value rank="0">0.0</real_value>
320
</constant>
321
</y_component>
322
</align_bc_with_cartesian>
323
</type>
324
</boundary_conditions>
325
<boundary_conditions name="sides">
326
<surface_ids>
327
<integer_value shape="1" rank="1">8</integer_value>
328
</surface_ids>
329
<type name="dirichlet">
330
<apply_weakly/>
331
<align_bc_with_cartesian>
332
<x_component>
333
<constant>
334
<real_value rank="0">0.0</real_value>
335
</constant>
336
</x_component>
337
<y_component>
338
<constant>
339
<real_value rank="0">0.0</real_value>
340
</constant>
341
</y_component>
342
</align_bc_with_cartesian>
343
</type>
344
</boundary_conditions>
345
<tensor_field name="Viscosity" rank="2">
346
<prescribed>
347
<mesh name="PressureMesh"/>
348
<value name="WholeMesh">
349
<isotropic>
350
<constant>
351
<real_value rank="0">1.0</real_value>
352
</constant>
353
</isotropic>
354
</value>
355
<output/>
356
</prescribed>
357
</tensor_field>
358
<vector_field name="Absorption" rank="1">
359
<diagnostic>
360
<mesh name="PressureMesh_Discontinuous"/>
361
<algorithm name="Internal" material_phase_support="multiple"/>
362
<output/>
363
<stat>
364
<include_in_stat/>
365
</stat>
366
<convergence>
367
<include_in_convergence/>
368
</convergence>
369
<detectors>
370
<include_in_detectors/>
371
</detectors>
372
<steady_state>
373
<include_in_steady_state/>
374
</steady_state>
375
</diagnostic>
376
<default_absorption/>
377
</vector_field>
378
<output/>
379
<stat>
380
<include_in_stat/>
381
<previous_time_step>
382
<exclude_from_stat/>
383
</previous_time_step>
384
<nonlinear_field>
385
<exclude_from_stat/>
386
</nonlinear_field>
387
</stat>
388
<convergence>
389
<include_in_convergence/>
390
</convergence>
391
<detectors>
392
<include_in_detectors/>
393
</detectors>
394
<steady_state>
395
<include_in_steady_state/>
396
</steady_state>
397
<consistent_interpolation/>
398
</prognostic>
399
</vector_field>
400
<scalar_field name="PhaseVolumeFraction" rank="0">
401
<prognostic>
402
<mesh name="PressureMesh"/>
403
<equation name="AdvectionDiffusion"/>
404
<spatial_discretisation>
405
<control_volumes>
406
<face_value name="FiniteElement">
407
<limit_face_value>
408
<limiter name="Sweby"/>
409
</limit_face_value>
410
</face_value>
411
<diffusion_scheme name="ElementGradient"/>
412
</control_volumes>
413
<conservative_advection>
414
<real_value rank="0">1.0</real_value>
415
</conservative_advection>
416
</spatial_discretisation>
417
<temporal_discretisation>
418
<theta>
419
<real_value rank="0">1.0</real_value>
420
</theta>
421
</temporal_discretisation>
422
<solver>
423
<iterative_method name="gmres">
424
<restart>
425
<integer_value rank="0">30</integer_value>
426
</restart>
427
</iterative_method>
428
<preconditioner name="sor"/>
429
<relative_error>
430
<real_value rank="0">1e-7</real_value>
431
</relative_error>
432
<max_iterations>
433
<integer_value rank="0">500</integer_value>
434
</max_iterations>
435
<never_ignore_solver_failures/>
436
<diagnostics>
437
<monitors/>
438
</diagnostics>
439
</solver>
440
<initial_condition name="WholeMesh">
441
<constant>
442
<real_value rank="0">0.0</real_value>
443
</constant>
444
</initial_condition>
445
<boundary_conditions name="inflow">
446
<surface_ids>
447
<integer_value shape="1" rank="1">6</integer_value>
448
</surface_ids>
449
<type name="dirichlet">
450
<apply_weakly/>
451
<constant>
452
<real_value rank="0">1.0</real_value>
453
</constant>
454
</type>
455
</boundary_conditions>
456
<output/>
457
<stat/>
458
<convergence>
459
<include_in_convergence/>
460
</convergence>
461
<detectors>
462
<include_in_detectors/>
463
</detectors>
464
<steady_state>
465
<include_in_steady_state/>
466
</steady_state>
467
<consistent_interpolation/>
468
</prognostic>
469
</scalar_field>
470
<multiphase_properties>
471
<relperm_type>
472
<Corey/>
473
</relperm_type>
474
</multiphase_properties>
475
</material_phase>
476
<material_phase name="phase2">
477
<scalar_field name="Pressure" rank="0">
478
<aliased material_phase_name="phase1" field_name="Pressure"/>
479
</scalar_field>
480
<scalar_field name="Density" rank="0">
481
<prognostic>
482
<mesh name="PressureMesh"/>
483
<spatial_discretisation>
484
<control_volumes>
485
<face_value name="FirstOrderUpwind"/>
486
</control_volumes>
487
<conservative_advection>
488
<real_value rank="0">0.0</real_value>
489
</conservative_advection>
490
</spatial_discretisation>
491
<temporal_discretisation>
492
<theta>
493
<real_value rank="0">1.0</real_value>
494
</theta>
495
</temporal_discretisation>
496
<initial_condition name="WholeMesh">
497
<constant>
498
<real_value rank="0">1.0</real_value>
499
</constant>
500
</initial_condition>
501
<output/>
502
<stat/>
503
<convergence>
504
<include_in_convergence/>
505
</convergence>
506
<detectors>
507
<include_in_detectors/>
508
</detectors>
509
<steady_state>
510
<include_in_steady_state/>
511
</steady_state>
512
<consistent_interpolation/>
513
</prognostic>
514
</scalar_field>
515
<vector_field name="Velocity" rank="1">
516
<prognostic>
517
<mesh name="VelocityMesh"/>
518
<equation name="LinearMomentum"/>
519
<spatial_discretisation>
520
<discontinuous_galerkin>
521
<viscosity_scheme>
522
<compact_discontinuous_galerkin/>
523
</viscosity_scheme>
524
<advection_scheme>
525
<upwind/>
526
<integrate_advection_by_parts>
527
<twice/>
528
</integrate_advection_by_parts>
529
</advection_scheme>
530
</discontinuous_galerkin>
531
<conservative_advection>
532
<real_value rank="0">1.0</real_value>
533
</conservative_advection>
534
</spatial_discretisation>
535
<temporal_discretisation>
536
<theta>
537
<real_value rank="0">1.0</real_value>
538
</theta>
539
<relaxation>
540
<real_value rank="0">1.0</real_value>
541
</relaxation>
542
</temporal_discretisation>
543
<solver>
544
<iterative_method name="gmres">
545
<restart>
546
<integer_value rank="0">30</integer_value>
547
</restart>
548
</iterative_method>
549
<preconditioner name="sor"/>
550
<relative_error>
551
<real_value rank="0">1e-7</real_value>
552
</relative_error>
553
<max_iterations>
554
<integer_value rank="0">500</integer_value>
555
</max_iterations>
556
<never_ignore_solver_failures/>
557
<diagnostics>
558
<monitors/>
559
</diagnostics>
560
</solver>
561
<initial_condition name="WholeMesh">
562
<constant>
563
<real_value shape="2" dim1="dim" rank="1">0.0 0.0</real_value>
564
</constant>
565
</initial_condition>
566
<boundary_conditions name="zero_inflow">
567
<surface_ids>
568
<integer_value shape="1" rank="1">6</integer_value>
569
</surface_ids>
570
<type name="dirichlet">
571
<apply_weakly/>
572
<align_bc_with_cartesian>
573
<x_component>
574
<constant>
575
<real_value rank="0">0.0</real_value>
576
</constant>
577
</x_component>
578
<y_component>
579
<constant>
580
<real_value rank="0">0.0</real_value>
581
</constant>
582
</y_component>
583
</align_bc_with_cartesian>
584
</type>
585
</boundary_conditions>
586
<boundary_conditions name="sides">
587
<surface_ids>
588
<integer_value shape="1" rank="1">8</integer_value>
589
</surface_ids>
590
<type name="dirichlet">
591
<apply_weakly/>
592
<align_bc_with_cartesian>
593
<x_component>
594
<constant>
595
<real_value rank="0">0.0</real_value>
596
</constant>
597
</x_component>
598
<y_component>
599
<constant>
600
<real_value rank="0">0.0</real_value>
601
</constant>
602
</y_component>
603
</align_bc_with_cartesian>
604
</type>
605
</boundary_conditions>
606
<tensor_field name="Viscosity" rank="2">
607
<prescribed>
608
<mesh name="PressureMesh"/>
609
<value name="WholeMesh">
610
<isotropic>
611
<constant>
612
<real_value rank="0">1.0</real_value>
613
</constant>
614
</isotropic>
615
</value>
616
<output/>
617
</prescribed>
618
</tensor_field>
619
<vector_field name="Absorption" rank="1">
620
<diagnostic>
621
<mesh name="PressureMesh_Discontinuous"/>
622
<algorithm name="Internal" material_phase_support="multiple"/>
623
<output/>
624
<stat>
625
<include_in_stat/>
626
</stat>
627
<convergence>
628
<include_in_convergence/>
629
</convergence>
630
<detectors>
631
<include_in_detectors/>
632
</detectors>
633
<steady_state>
634
<include_in_steady_state/>
635
</steady_state>
636
</diagnostic>
637
<default_absorption/>
638
</vector_field>
639
<output/>
640
<stat>
641
<include_in_stat/>
642
<previous_time_step>
643
<exclude_from_stat/>
644
</previous_time_step>
645
<nonlinear_field>
646
<exclude_from_stat/>
647
</nonlinear_field>
648
</stat>
649
<convergence>
650
<include_in_convergence/>
651
</convergence>
652
<detectors>
653
<include_in_detectors/>
654
</detectors>
655
<steady_state>
656
<include_in_steady_state/>
657
</steady_state>
658
<consistent_interpolation/>
659
</prognostic>
660
</vector_field>
661
<scalar_field name="PhaseVolumeFraction" rank="0">
662
<prognostic>
663
<mesh name="PressureMesh"/>
664
<equation name="AdvectionDiffusion"/>
665
<spatial_discretisation>
666
<control_volumes>
667
<face_value name="FiniteElement">
668
<limit_face_value>
669
<limiter name="Sweby"/>
670
</limit_face_value>
671
</face_value>
672
<diffusion_scheme name="ElementGradient"/>
673
</control_volumes>
674
<conservative_advection>
675
<real_value rank="0">0.0</real_value>
676
</conservative_advection>
677
</spatial_discretisation>
678
<temporal_discretisation>
679
<theta>
680
<real_value rank="0">1.0</real_value>
681
</theta>
682
</temporal_discretisation>
683
<solver>
684
<iterative_method name="gmres">
685
<restart>
686
<integer_value rank="0">30</integer_value>
687
</restart>
688
</iterative_method>
689
<preconditioner name="sor"/>
690
<relative_error>
691
<real_value rank="0">1e-7</real_value>
692
</relative_error>
693
<max_iterations>
694
<integer_value rank="0">500</integer_value>
695
</max_iterations>
696
<never_ignore_solver_failures/>
697
<diagnostics>
698
<monitors/>
699
</diagnostics>
700
</solver>
701
<initial_condition name="WholeMesh">
702
<constant>
703
<real_value rank="0">1.0</real_value>
704
</constant>
705
</initial_condition>
706
<boundary_conditions name="inflow">
707
<surface_ids>
708
<integer_value shape="1" rank="1">6</integer_value>
709
</surface_ids>
710
<type name="dirichlet">
711
<apply_weakly/>
712
<constant>
713
<real_value rank="0">0.0</real_value>
714
</constant>
715
</type>
716
</boundary_conditions>
717
<output/>
718
<stat/>
719
<convergence>
720
<include_in_convergence/>
721
</convergence>
722
<detectors>
723
<include_in_detectors/>
724
</detectors>
725
<steady_state>
726
<include_in_steady_state/>
727
</steady_state>
728
<consistent_interpolation/>
729
</prognostic>
730
</scalar_field>
731
<multiphase_properties>
732
<relperm_type>
733
<Corey/>
734
</relperm_type>
735
</multiphase_properties>
736
</material_phase>
737
<material_phase name="Component1">
738
<is_multiphase_component>
739
<alpha_beta>
740
<real_value rank="0">1.0</real_value>
741
</alpha_beta>
742
<KComp_Sigmoid>
743
<K_Comp>
744
<real_value rank="0">4.5</real_value>
745
</K_Comp>
746
</KComp_Sigmoid>
747
<Comp_Sum2One>
748
<Relaxation_Coefficient>
749
<real_value rank="0">1.0</real_value>
750
</Relaxation_Coefficient>
751
</Comp_Sum2One>
752
</is_multiphase_component>
753
<vector_field name="Velocity" rank="1">
754
<aliased material_phase_name="phase1" field_name="Velocity"/>
755
</vector_field>
756
<scalar_field material_phase_name="phase1" name="ComponentMassFractionPhase1" rank="0">
757
<prognostic>
758
<mesh name="PressureMesh"/>
759
<equation_of_state>
760
<incompressible>
761
<linear>
762
<all_equal>
763
<real_value rank="0">1.0</real_value>
764
</all_equal>
765
</linear>
766
</incompressible>
767
</equation_of_state>
768
<equation name="AdvectionDiffusion"/>
769
<spatial_discretisation>
770
<control_volumes>
771
<face_value name="FirstOrderUpwind"/>
772
<diffusion_scheme name="ElementGradient"/>
773
</control_volumes>
774
<conservative_advection>
775
<real_value rank="0">1.0</real_value>
776
</conservative_advection>
777
</spatial_discretisation>
778
<temporal_discretisation>
779
<theta>
780
<real_value rank="0">1.0</real_value>
781
</theta>
782
</temporal_discretisation>
783
<solver>
784
<iterative_method name="gmres">
785
<restart>
786
<integer_value rank="0">30</integer_value>
787
</restart>
788
</iterative_method>
789
<preconditioner name="sor"/>
790
<relative_error>
791
<real_value rank="0">1e-7</real_value>
792
</relative_error>
793
<max_iterations>
794
<integer_value rank="0">5000</integer_value>
795
</max_iterations>
796
<never_ignore_solver_failures/>
797
<diagnostics>
798
<monitors/>
799
</diagnostics>
800
</solver>
801
<initial_condition name="WholeMesh">
802
<constant>
803
<real_value rank="0">0.0</real_value>
804
</constant>
805
</initial_condition>
806
<boundary_conditions name="inflow_C1P1">
807
<surface_ids>
808
<integer_value shape="1" rank="1">6</integer_value>
809
</surface_ids>
810
<type name="dirichlet">
811
<apply_weakly/>
812
<constant>
813
<real_value rank="0">1.0</real_value>
814
</constant>
815
</type>
816
</boundary_conditions>
817
<output/>
818
<stat/>
819
<convergence>
820
<include_in_convergence/>
821
</convergence>
822
<detectors>
823
<include_in_detectors/>
824
</detectors>
825
<steady_state>
826
<include_in_steady_state/>
827
</steady_state>
828
<consistent_interpolation/>
829
</prognostic>
830
</scalar_field>
831
<scalar_field material_phase_name="phase2" name="ComponentMassFractionPhase2" rank="0">
832
<prognostic>
833
<mesh name="PressureMesh"/>
834
<equation_of_state>
835
<incompressible>
836
<linear>
837
<all_equal>
838
<real_value rank="0">1.0</real_value>
839
</all_equal>
840
</linear>
841
</incompressible>
842
</equation_of_state>
843
<equation name="AdvectionDiffusion"/>
844
<spatial_discretisation>
845
<control_volumes>
846
<face_value name="FirstOrderUpwind"/>
847
<diffusion_scheme name="ElementGradient"/>
848
</control_volumes>
849
<conservative_advection>
850
<real_value rank="0">1.0</real_value>
851
</conservative_advection>
852
</spatial_discretisation>
853
<temporal_discretisation>
854
<theta>
855
<real_value rank="0">1.0</real_value>
856
</theta>
857
</temporal_discretisation>
858
<solver>
859
<iterative_method name="gmres">
860
<restart>
861
<integer_value rank="0">30</integer_value>
862
</restart>
863
</iterative_method>
864
<preconditioner name="sor"/>
865
<relative_error>
866
<real_value rank="0">1e-7</real_value>
867
</relative_error>
868
<max_iterations>
869
<integer_value rank="0">1000</integer_value>
870
</max_iterations>
871
<never_ignore_solver_failures/>
872
<diagnostics>
873
<monitors/>
874
</diagnostics>
875
</solver>
876
<initial_condition name="WholeMesh">
877
<constant>
878
<real_value rank="0">0.0</real_value>
879
</constant>
880
</initial_condition>
881
<boundary_conditions name="inflow_C1P2">
882
<surface_ids>
883
<integer_value shape="1" rank="1">6</integer_value>
884
</surface_ids>
885
<type name="dirichlet">
886
<apply_weakly/>
887
<constant>
888
<real_value rank="0">0.0</real_value>
889
</constant>
890
</type>
891
</boundary_conditions>
892
<output/>
893
<stat/>
894
<convergence>
895
<include_in_convergence/>
896
</convergence>
897
<detectors>
898
<include_in_detectors/>
899
</detectors>
900
<steady_state>
901
<include_in_steady_state/>
902
</steady_state>
903
<consistent_interpolation/>
904
</prognostic>
905
</scalar_field>
906
</material_phase>
907
<material_phase name="Component2">
908
<is_multiphase_component>
909
<alpha_beta>
910
<real_value rank="0">1.0</real_value>
911
</alpha_beta>
912
<KComp_Sigmoid>
913
<K_Comp>
914
<real_value rank="0">0.2</real_value>
915
</K_Comp>
916
</KComp_Sigmoid>
917
<Comp_Sum2One>
918
<Relaxation_Coefficient>
919
<real_value rank="0">1.0</real_value>
920
</Relaxation_Coefficient>
921
</Comp_Sum2One>
922
</is_multiphase_component>
923
<vector_field name="Velocity" rank="1">
924
<aliased material_phase_name="phase2" field_name="Velocity"/>
925
</vector_field>
926
<scalar_field material_phase_name="phase1" name="ComponentMassFractionPhase1" rank="0">
927
<prognostic>
928
<mesh name="PressureMesh"/>
929
<equation_of_state>
930
<incompressible>
931
<linear>
932
<all_equal>
933
<real_value rank="0">1.0</real_value>
934
</all_equal>
935
</linear>
936
</incompressible>
937
</equation_of_state>
938
<equation name="AdvectionDiffusion"/>
939
<spatial_discretisation>
940
<control_volumes>
941
<face_value name="FirstOrderUpwind"/>
942
<diffusion_scheme name="ElementGradient"/>
943
</control_volumes>
944
<conservative_advection>
945
<real_value rank="0">1.0</real_value>
946
</conservative_advection>
947
</spatial_discretisation>
948
<temporal_discretisation>
949
<theta>
950
<real_value rank="0">1.0</real_value>
951
</theta>
952
</temporal_discretisation>
953
<solver>
954
<iterative_method name="gmres">
955
<restart>
956
<integer_value rank="0">30</integer_value>
957
</restart>
958
</iterative_method>
959
<preconditioner name="prometheus"/>
960
<relative_error>
961
<real_value rank="0">1e-7</real_value>
962
</relative_error>
963
<max_iterations>
964
<integer_value rank="0">1000</integer_value>
965
</max_iterations>
966
<never_ignore_solver_failures/>
967
<diagnostics>
968
<monitors/>
969
</diagnostics>
970
</solver>
971
<initial_condition name="WholeMesh">
972
<constant>
973
<real_value rank="0">0.0</real_value>
974
</constant>
975
</initial_condition>
976
<boundary_conditions name="inflow_C2P1">
977
<surface_ids>
978
<integer_value shape="1" rank="1">6</integer_value>
979
</surface_ids>
980
<type name="dirichlet">
981
<apply_weakly/>
982
<constant>
983
<real_value rank="0">0.0</real_value>
984
</constant>
985
</type>
986
</boundary_conditions>
987
<output/>
988
<stat/>
989
<convergence>
990
<include_in_convergence/>
991
</convergence>
992
<detectors>
993
<include_in_detectors/>
994
</detectors>
995
<steady_state>
996
<include_in_steady_state/>
997
</steady_state>
998
<consistent_interpolation/>
999
</prognostic>
1000
</scalar_field>
1001
<scalar_field material_phase_name="phase2" name="ComponentMassFractionPhase2" rank="0">
1002
<prognostic>
1003
<mesh name="PressureMesh"/>
1004
<equation_of_state>
1005
<incompressible>
1006
<linear>
1007
<all_equal>
1008
<real_value rank="0">1.0</real_value>
1009
</all_equal>
1010
</linear>
1011
</incompressible>
1012
</equation_of_state>
1013
<equation name="AdvectionDiffusion"/>
1014
<spatial_discretisation>
1015
<control_volumes>
1016
<face_value name="FirstOrderUpwind"/>
1017
<diffusion_scheme name="ElementGradient"/>
1018
</control_volumes>
1019
<conservative_advection>
1020
<real_value rank="0">1.0</real_value>
1021
</conservative_advection>
1022
</spatial_discretisation>
1023
<temporal_discretisation>
1024
<theta>
1025
<real_value rank="0">1.0</real_value>
1026
</theta>
1027
</temporal_discretisation>
1028
<solver>
1029
<iterative_method name="gmres">
1030
<restart>
1031
<integer_value rank="0">30</integer_value>
1032
</restart>
1033
</iterative_method>
1034
<preconditioner name="prometheus"/>
1035
<relative_error>
1036
<real_value rank="0">1e-7</real_value>
1037
</relative_error>
1038
<max_iterations>
1039
<integer_value rank="0">1000</integer_value>
1040
</max_iterations>
1041
<never_ignore_solver_failures/>
1042
<diagnostics>
1043
<monitors/>
1044
</diagnostics>
1045
</solver>
1046
<initial_condition name="WholeMesh">
1047
<constant>
1048
<real_value rank="0">1.0</real_value>
1049
</constant>
1050
</initial_condition>
1051
<boundary_conditions name="inflow_C2P2">
1052
<surface_ids>
1053
<integer_value shape="1" rank="1">6</integer_value>
1054
</surface_ids>
1055
<type name="dirichlet">
1056
<apply_weakly/>
1057
<constant>
1058
<real_value rank="0">0.0</real_value>
1059
</constant>
1060
</type>
1061
</boundary_conditions>
1062
<output/>
1063
<stat/>
1064
<convergence>
1065
<include_in_convergence/>
1066
</convergence>
1067
<detectors>
1068
<include_in_detectors/>
1069
</detectors>
1070
<steady_state>
1071
<include_in_steady_state/>
1072
</steady_state>
1073
<consistent_interpolation/>
1074
</prognostic>
1075
</scalar_field>
1076
</material_phase>
1077
<porous_media>
1078
<scalar_field name="Porosity">
1079
<prescribed>
1080
<mesh name="P0DG"/>
1081
<value name="WholeMesh">
1082
<constant>
1083
<real_value rank="0">0.5</real_value>
1084
</constant>
1085
</value>
1086
<output/>
1087
<stat/>
1088
<detectors>
1089
<exclude_from_detectors/>
1090
</detectors>
1091
</prescribed>
1092
</scalar_field>
1093
<scalar_field name="Permeability">
1094
<prescribed>
1095
<mesh name="P0DG"/>
1096
<value name="WholeMesh">
1097
<constant>
1098
<real_value rank="0">1.0</real_value>
1099
</constant>
1100
</value>
1101
<output/>
1102
<stat/>
1103
<detectors>
1104
<exclude_from_detectors/>
1105
</detectors>
1106
</prescribed>
1107
</scalar_field>
1108
</porous_media>
1109
</multiphase_options>
Loggerhead is a web-based interface for Breezy
Version: 2.0.1