← Back to branch summary

~reducedmodelling/fluidity/ReducedModel_running

« back to all changes in this revision

Viewing changes to tests/mphase_sedimentation_2d/mphase_sedimentation_2d.xml

  • Committer: Christian Jacobs
  • Date: 2011-09-29 22:38:40 UTC
  • mfrom: (3408.1.160 ctjacobs-multiphase)
  • Revision ID: c.jacobs10@imperial.ac.uk-20110929223840-yiwz4nt5hfuraqey
Merging ctjacobs-multiphase branch revisions (from r3533 to r3568 inclusive) into trunk. Branch queue passes all unit, short and medium tests on buildbot.

Summary of changes:
===================

- Added multiphase support to the CV discretisation of the Pressure field. A multiphase MMS longtest for the P0-P1CV element pair has been added and shows first order convergence for Velocity and Pressure.

- When checking the field_bc_type in assemble_divergence_matrix_cv, we now compare against the more meaningful integer parameters BC_TYPE_WEAKDIRICHLET (=1), BC_TYPE_NO_NORMAL_FLOW (=2), etc.

- In Momentum_CG.F90: Using dnvfrac_t instead of du_t if PhaseVolumeFraction and Velocity are on different meshes.

In Divergence_Matrix_CG.F90: Using dnvfrac_t instead of dfield_t if nvfrac and field are on different meshes.

- A few minor changes to multiphase test cases and the manual.

- Added a multiphase diagnostic algorithm called apparent_density, which computes Density*PhaseVolumeFraction.

- Bug fix: The logical 'prognostic_velocity' was assigned a value of .true. regardless of whether the Velocity field was prognostic, prescribed or aliased. Now fixed.

- Bug fix for lp:854398: In Multimaterials.F90, we now only count the number of diagnostic MaterialVolumeFraction fields in the submaterials array instead of the entire state array.

- Bug fix for the particle_reynolds_number diagnostic algorithm: Including the continuous phase's PhaseVolumeFraction and using the magnitude of the relative velocity. Also updated the related comments in the schema.

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/mphase_sedimentation_2d/mphase_sedimentation_2d.xml
3
3
 
4
4
<testproblem>
5
5
   
6
 
  <name>mphase_sedimentation</name>
 
6
  <name>mphase_sedimentation_2d</name>
7
7
  <owner userid="ctj10"/>
8
8
  <tags>flml</tags>
9
9
  
20
20
  <variables>
21
21
    <variable name="f1_vfrac_max" language="python">
22
22
from fluidity_tools import stat_parser
23
 
s = stat_parser("mphase_sedimentation.stat")
 
23
s = stat_parser("mphase_sedimentation_2d.stat")
24
24
f1_vfrac_max = s["Fluid1"]["PhaseVolumeFraction"]["max"][-1]
25
25
    </variable>
26
26
    
27
27
    <variable name="f1_vfrac_min" language="python">
28
28
from fluidity_tools import stat_parser
29
 
s = stat_parser("mphase_sedimentation.stat")
 
29
s = stat_parser("mphase_sedimentation_2d.stat")
30
30
f1_vfrac_min = s["Fluid1"]["PhaseVolumeFraction"]["min"][-1]
31
31
    </variable>
32
32
    
33
33
    <variable name="f2_vfrac_max" language="python">
34
34
from fluidity_tools import stat_parser
35
 
s = stat_parser("mphase_sedimentation.stat")
 
35
s = stat_parser("mphase_sedimentation_2d.stat")
36
36
f2_vfrac_max = s["Fluid2"]["PhaseVolumeFraction"]["max"][-1]
37
37
    </variable>
38
38
    
39
39
    <variable name="f2_vfrac_min" language="python">
40
40
from fluidity_tools import stat_parser
41
 
s = stat_parser("mphase_sedimentation.stat")
 
41
s = stat_parser("mphase_sedimentation_2d.stat")
42
42
f2_vfrac_min = s["Fluid2"]["PhaseVolumeFraction"]["min"][-1]
43
43
    </variable>
44
44
    
45
45
    <variable name="f1_vfrac_integral" language="python">
46
46
from fluidity_tools import stat_parser
47
 
s = stat_parser("mphase_sedimentation.stat")
 
47
s = stat_parser("mphase_sedimentation_2d.stat")
48
48
f1_vfrac_integral = s["Fluid1"]["PhaseVolumeFraction"]["integral"][-1]
49
49
    </variable>
50
50
    
51
51
    <variable name="f2_vfrac_integral" language="python">
52
52
from fluidity_tools import stat_parser
53
 
s = stat_parser("mphase_sedimentation.stat")
 
53
s = stat_parser("mphase_sedimentation_2d.stat")
54
54
f2_vfrac_integral = s["Fluid2"]["PhaseVolumeFraction"]["integral"][-1]
55
55
    </variable>
56
56
    
57
57
    <variable name="divergence" language="python">
58
58
from fluidity_tools import stat_parser
59
 
s = stat_parser("mphase_sedimentation.stat")
 
59
s = stat_parser("mphase_sedimentation_2d.stat")
60
60
divergence = s["Fluid1"]["SumVelocityDivergence"]["max"][-1]
61
61
    </variable>
62
62
    
63
63
    <variable name="f1_vfrac" language="python">
64
64
import vtktools
65
65
 
66
 
file = vtktools.vtu('mphase_sedimentation_601.vtu')
 
66
file = vtktools.vtu('mphase_sedimentation_2d_601.vtu')
67
67
file.GetFieldNames()
68
68
f1_vfrac = file.GetScalarField('Fluid1::PhaseVolumeFraction')
69
69
    </variable>
71
71
    <variable name="f2_vfrac" language="python">
72
72
import vtktools
73
73
 
74
 
file = vtktools.vtu('mphase_sedimentation_601.vtu')
 
74
file = vtktools.vtu('mphase_sedimentation_2d_601.vtu')
75
75
file.GetFieldNames()
76
76
f2_vfrac = file.GetScalarField('Fluid2::PhaseVolumeFraction')
77
77
    </variable>
80
80
    <variable name="positions" language="python">
81
81
import vtktools
82
82
 
83
 
file = vtktools.vtu('mphase_sedimentation_601.vtu')
 
83
file = vtktools.vtu('mphase_sedimentation_2d_601.vtu')
84
84
file.GetFieldNames()
85
85
p = file.GetLocations()
86
86
positions = p[:,1]