diff --git a/dumux/freeflow/navierstokes/staggered/localresidual.hh b/dumux/freeflow/navierstokes/staggered/localresidual.hh
index 485204b8d4ae494a32267c55c14e63be81e8b707..73821e002df834c0299a82d8d60803f882e94aea 100644
--- a/dumux/freeflow/navierstokes/staggered/localresidual.hh
+++ b/dumux/freeflow/navierstokes/staggered/localresidual.hh
@@ -302,17 +302,18 @@ protected:
}
else if(bcTypes.isNeumann(Indices::velocity(scvf.directionIndex())))
{
- // set a given Neumann flux for the face on the boundary itself
+ // the source term has already been accounted for, here we
+ // add a given Neumann flux for the face on the boundary itself ...
const auto extrusionFactor = elemVolVars[scvf.insideScvIdx()].extrusionFactor();
- residual = problem.neumann(element, fvGeometry, elemVolVars, elemFaceVars, scvf)[Indices::velocity(scvf.directionIndex())]
+ residual += problem.neumann(element, fvGeometry, elemVolVars, elemFaceVars, scvf)[Indices::velocity(scvf.directionIndex())]
* extrusionFactor * scvf.area();
- // treat the remaining (frontal and lateral) faces of the staggered control volume
+ // ... and treat the fluxes of the remaining (frontal and lateral) faces of the staggered control volume
residual += computeFluxForFace(problem, element, scvf, fvGeometry, elemVolVars, elemFaceVars, elemFluxVarsCache);
}
else if(bcTypes.isSymmetry())
{
- // For symmetry boundary conditions, there is no flow accross the boundary and
+ // for symmetry boundary conditions, there is no flow accross the boundary and
// we therefore treat it like a Dirichlet boundary conditions with zero velocity
const Scalar velocity = elemFaceVars[scvf].velocitySelf();
const Scalar fixedValue = 0.0;
@@ -320,7 +321,8 @@ protected:
}
else if(bcTypes.isDirichlet(Indices::pressureIdx))
{
- // If none of the above conditions apply, we are at an "fixed pressure" boundary for which the velocity needs to be assembled
+ // if none of the above conditions apply, we are at an "fixed pressure" boundary for which the resdiual of the momentum balance needs to be assembled
+ // as if it where inside the domain and not on the boundary (source term has already been acounted for)
residual += computeFluxForFace(problem, element, scvf, fvGeometry, elemVolVars, elemFaceVars, elemFluxVarsCache);
}
else
diff --git a/dumux/multidomain/boundary/stokesdarcy/couplingdata.hh b/dumux/multidomain/boundary/stokesdarcy/couplingdata.hh
index 5846529cb7ae160a9c37a9cce145de35f5b81382..259d7da18e1e0406372e99e487246246b2ae190d 100644
--- a/dumux/multidomain/boundary/stokesdarcy/couplingdata.hh
+++ b/dumux/multidomain/boundary/stokesdarcy/couplingdata.hh
@@ -294,9 +294,14 @@ public:
}
else // use pressure reconstruction for single phase models
{
+ // v = -K/mu * (gradP + rho*g)
const Scalar velocity = stokesElemFaceVars[scvf].velocitySelf();
- const Scalar pressureInterFace = scvf.directionSign() * velocity * stokesContext.volVars.viscosity(darcyPhaseIdx)/darcyPermeability(scvf) * (stokesContext.element.geometry().center() - scvf.center()).two_norm() + darcyPressure;
- momentumFlux = pressureInterFace;
+ const Scalar mu = stokesContext.volVars.viscosity(darcyPhaseIdx);
+ const Scalar rho = stokesContext.volVars.density(darcyPhaseIdx);
+ const Scalar distance = (stokesContext.element.geometry().center() - scvf.center()).two_norm();
+ const Scalar g = couplingManager_.problem(darcyIdx).gravity()[scvf.directionIndex()];
+ const Scalar interfacePressure = ((scvf.directionSign() * velocity * (mu/darcyPermeability(scvf))) + rho * g) * distance + darcyPressure;
+ momentumFlux = interfacePressure;
}
// normalize pressure
diff --git a/dumux/multidomain/boundary/stokesdarcy/couplingmanager.hh b/dumux/multidomain/boundary/stokesdarcy/couplingmanager.hh
index d218acd81a0cb2ec7ddde18bccdbd9690b54a9cf..04b704f33e35d9634304f8240fa73d4d43b5a42a 100644
--- a/dumux/multidomain/boundary/stokesdarcy/couplingmanager.hh
+++ b/dumux/multidomain/boundary/stokesdarcy/couplingmanager.hh
@@ -30,6 +30,8 @@
#include <utility>
#include <memory>
+#include <dune/common/float_cmp.hh>
+#include <dune/common/exceptions.hh>
#include <dumux/common/properties.hh>
#include <dumux/multidomain/staggeredcouplingmanager.hh>
@@ -119,7 +121,6 @@ public:
using ParentType::couplingStencil;
using ParentType::updateCouplingContext;
- using DarcyGridVariables = GridVariables<2>;
using CouplingData = StokesDarcyCouplingData<MDTraits, StokesDarcyCouplingManager<MDTraits>>;
//! Constructor
@@ -137,6 +138,9 @@ public:
std::shared_ptr<const Problem<darcyIdx>> darcyProblem,
const SolutionVector& curSol)
{
+ if(Dune::FloatCmp::ne(stokesProblem->gravity(), darcyProblem->gravity()))
+ DUNE_THROW(Dune::InvalidStateException, "Both models must use the same gravity vector");
+
ParentType::init(std::make_tuple(stokesProblem, stokesProblem, darcyProblem));
this->curSol() = curSol;
couplingData_ = std::make_shared<CouplingData>(*this);
diff --git a/test/references/test_stokes1p2cnidarcy2p2cnihorizontal_darcy-reference.vtu b/test/references/test_stokes1p2cnidarcy2p2cnihorizontal_darcy-reference.vtu
index a4ebe346a5e3b202f1be4196c025c1cfa516ff19..8172835a1cca9873dba33103abe348141637ad2a 100644
--- a/test/references/test_stokes1p2cnidarcy2p2cnihorizontal_darcy-reference.vtu
+++ b/test/references/test_stokes1p2cnidarcy2p2cnihorizontal_darcy-reference.vtu
@@ -2,31 +2,31 @@
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
<UnstructuredGrid>
<Piece NumberOfCells="64" NumberOfPoints="81">
- <CellData Scalars="S_n" Vectors="velocity_liquid (m/s)">
- <DataArray type="Float32" Name="S_n" NumberOfComponents="1" format="ascii">
- 0.231386 0.231386 0.231387 0.231387 0.231387 0.231387 0.231387 0.231387 0.816991 0.816991 0.816993 0.816994
- 0.816995 0.816997 0.816997 0.816998 0.937394 0.937394 0.937393 0.937393 0.937393 0.937393 0.937393 0.937393
- 0.966181 0.966171 0.966161 0.966154 0.966149 0.966145 0.966142 0.966141 0.977391 0.977355 0.977326 0.977304
- 0.977287 0.977273 0.977264 0.977261 0.983867 0.983784 0.983724 0.983681 0.983647 0.983619 0.983599 0.983594
- 0.990183 0.989914 0.989734 0.989612 0.989514 0.989436 0.989376 0.989364 1 1 1 1
- 1 1 1 1
+ <CellData Scalars="porosity" Vectors="velocity_liquid (m/s)">
+ <DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii">
+ 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41
+ 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41
+ 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41
+ 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41
+ 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41
+ 0.41 0.41 0.41 0.41
+ </DataArray>
+ <DataArray type="Float32" Name="pc" NumberOfComponents="1" format="ascii">
+ 1346.28 1346.28 1346.28 1346.28 1346.28 1346.28 1346.28 1346.28 2053.9 2053.9 2053.9 2053.9
+ 2053.91 2053.91 2053.91 2053.91 2526.12 2526.11 2526.11 2526.11 2526.11 2526.11 2526.11 2526.11
+ 2849.26 2849.08 2848.92 2848.8 2848.71 2848.63 2848.58 2848.57 3100.58 3099.48 3098.61 3097.95
+ 3097.44 3097.04 3096.75 3096.67 3352.96 3348.7 3345.68 3343.5 3341.79 3340.41 3339.39 3339.16
+ 3721.5 3705.55 3694.95 3687.67 3681.92 3677.26 3673.76 3673.04 4302.49 4302.49 4302.49 4302.49
+ 4302.49 4302.49 4302.49 4302.49
</DataArray>
<DataArray type="Float32" Name="S_w" NumberOfComponents="1" format="ascii">
- 0.768614 0.768614 0.768613 0.768613 0.768613 0.768613 0.768613 0.768613 0.183009 0.183009 0.183007 0.183006
+ 0.768614 0.768614 0.768614 0.768613 0.768613 0.768613 0.768613 0.768613 0.183009 0.183009 0.183007 0.183006
0.183005 0.183004 0.183003 0.183002 0.0626056 0.0626064 0.0626069 0.062607 0.0626069 0.0626068 0.0626066 0.0626065
0.0338188 0.0338293 0.0338386 0.0338457 0.0338511 0.0338553 0.0338582 0.0338591 0.0226087 0.0226453 0.0226743 0.0226961
- 0.0227132 0.0227267 0.0227364 0.022739 0.016133 0.0162165 0.0162761 0.0163191 0.0163533 0.0163808 0.0164012 0.0164059
- 0.0098168 0.0100864 0.0102655 0.0103884 0.0104857 0.0105644 0.0106236 0.0106356 0 0 0 0
+ 0.0227132 0.0227267 0.0227364 0.022739 0.0161331 0.0162165 0.0162761 0.0163191 0.0163533 0.0163808 0.0164012 0.0164059
+ 0.00981681 0.0100864 0.0102655 0.0103884 0.0104857 0.0105644 0.0106236 0.0106356 0 0 0 0
0 0 0 0
</DataArray>
- <DataArray type="Float32" Name="p_n" NumberOfComponents="1" format="ascii">
- 100002 100002 100002 100002 100002 100002 100002 100002 100002 100002 100002 100002
- 100002 100002 100002 100002 100001 100001 100001 100001 100001 100001 100001 100001
- 100001 100001 100001 100001 100001 100001 100001 100001 100000 100000 100000 100000
- 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
- 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
- 100000 100000 100000 100000
- </DataArray>
<DataArray type="Float32" Name="p_w" NumberOfComponents="1" format="ascii">
98656.1 98656.1 98656.1 98656.1 98656.1 98656.1 98656.1 98656.1 97947.7 97947.7 97947.7 97947.7
97947.7 97947.7 97947.7 97947.7 97474.9 97474.9 97474.9 97474.9 97474.9 97474.9 97474.9 97474.9
@@ -35,14 +35,6 @@
96278.6 96294.5 96305.1 96312.4 96318.2 96322.8 96326.3 96327 95697.5 95697.5 95697.5 95697.5
95697.5 95697.5 95697.5 95697.5
</DataArray>
- <DataArray type="Float32" Name="pc" NumberOfComponents="1" format="ascii">
- 1346.28 1346.28 1346.28 1346.28 1346.28 1346.28 1346.28 1346.28 2053.9 2053.9 2053.9 2053.9
- 2053.91 2053.91 2053.91 2053.91 2526.12 2526.11 2526.11 2526.11 2526.11 2526.11 2526.11 2526.11
- 2849.26 2849.08 2848.92 2848.8 2848.71 2848.63 2848.58 2848.57 3100.58 3099.48 3098.61 3097.95
- 3097.44 3097.04 3096.75 3096.67 3352.96 3348.7 3345.68 3343.5 3341.79 3340.41 3339.39 3339.16
- 3721.5 3705.55 3694.95 3687.67 3681.92 3677.26 3673.76 3673.04 4302.49 4302.49 4302.49 4302.49
- 4302.49 4302.49 4302.49 4302.49
- </DataArray>
<DataArray type="Float32" Name="rho_w" NumberOfComponents="1" format="ascii">
999.385 999.387 999.389 999.392 999.395 999.397 999.399 999.4 999.352 999.355 999.36 999.367
999.373 999.379 999.383 999.385 999.308 999.315 999.327 999.34 999.352 999.363 999.37 999.373
@@ -51,14 +43,6 @@
999.123 999.192 999.25 999.289 999.321 999.346 999.364 999.37 271.631 328.146 369.108 400.775
427.636 450.858 469.655 473.181
</DataArray>
- <DataArray type="Float32" Name="rho_n" NumberOfComponents="1" format="ascii">
- 1.21091 1.21095 1.21104 1.21115 1.21127 1.21137 1.21145 1.21149 1.20961 1.20973 1.20993 1.21018
- 1.21043 1.21065 1.21081 1.21089 1.20787 1.20816 1.20862 1.20913 1.20961 1.21002 1.2103 1.21044
- 1.20607 1.20665 1.20746 1.20829 1.20904 1.20964 1.21005 1.21025 1.2044 1.2054 1.2066 1.20773
- 1.2087 1.20947 1.20999 1.21022 1.20295 1.20445 1.20601 1.20737 1.20851 1.20941 1.21002 1.21025
- 1.20176 1.20376 1.20561 1.20714 1.2084 1.20939 1.21008 1.21031 1.20645 1.20833 1.20996 1.2113
- 1.21238 1.21324 1.21382 1.21399
- </DataArray>
<DataArray type="Float32" Name="mob_w" NumberOfComponents="1" format="ascii">
343.952 343.852 343.673 343.445 343.206 342.991 342.83 342.744 4.84623 4.8427 4.83658 4.8292
4.82183 4.81545 4.81082 4.80841 0.196025 0.195688 0.195146 0.194538 0.193965 0.193492 0.19316 0.192991
@@ -67,36 +51,92 @@
0.000176098 0.000203354 0.000222036 0.000234824 0.000245161 0.000253738 0.000260397 0.000261626 0 0 0 0
0 0 0 0
</DataArray>
+ <DataArray type="Float32" Name="S_n" NumberOfComponents="1" format="ascii">
+ 0.231386 0.231386 0.231387 0.231387 0.231387 0.231387 0.231387 0.231387 0.816991 0.816991 0.816993 0.816994
+ 0.816995 0.816997 0.816997 0.816998 0.937394 0.937394 0.937393 0.937393 0.937393 0.937393 0.937393 0.937393
+ 0.966181 0.966171 0.966161 0.966154 0.966149 0.966145 0.966142 0.966141 0.977391 0.977355 0.977326 0.977304
+ 0.977287 0.977273 0.977264 0.977261 0.983867 0.983784 0.983724 0.983681 0.983647 0.983619 0.983599 0.983594
+ 0.990183 0.989914 0.989734 0.989612 0.989514 0.989436 0.989376 0.989364 1 1 1 1
+ 1 1 1 1
+ </DataArray>
+ <DataArray type="Float32" Name="p_n" NumberOfComponents="1" format="ascii">
+ 100002 100002 100002 100002 100002 100002 100002 100002 100002 100002 100002 100002
+ 100002 100002 100002 100002 100001 100001 100001 100001 100001 100001 100001 100001
+ 100001 100001 100001 100001 100001 100001 100001 100001 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000
+ </DataArray>
+ <DataArray type="Float32" Name="rho_n" NumberOfComponents="1" format="ascii">
+ 1.21091 1.21095 1.21104 1.21115 1.21127 1.21137 1.21145 1.21149 1.20961 1.20973 1.20993 1.21018
+ 1.21043 1.21065 1.21081 1.21089 1.20787 1.20816 1.20862 1.20913 1.20961 1.21002 1.2103 1.21044
+ 1.20607 1.20665 1.20746 1.20829 1.20904 1.20964 1.21005 1.21025 1.2044 1.2054 1.2066 1.20773
+ 1.2087 1.20947 1.20999 1.21022 1.20295 1.20445 1.20601 1.20737 1.20851 1.20941 1.21002 1.21025
+ 1.20176 1.20376 1.20561 1.20714 1.2084 1.20939 1.21008 1.21031 1.20645 1.20833 1.20996 1.2113
+ 1.21238 1.21324 1.21382 1.21399
+ </DataArray>
<DataArray type="Float32" Name="mob_n" NumberOfComponents="1" format="ascii">
- 3377.51 3377.59 3377.74 3377.93 3378.13 3378.31 3378.44 3378.51 41239.7 41242.2 41246.4 41251.6
+ 3377.51 3377.59 3377.74 3377.93 3378.13 3378.3 3378.44 3378.51 41239.7 41242.2 41246.4 41251.6
41256.8 41261.3 41264.6 41266.3 52669.1 52676.4 52688.3 52701.6 52714.1 52724.6 52732 52735.7
55280.4 55295.4 55317 55339.2 55359.2 55375.2 55386.3 55391.6 55915 55939.6 55971.4 56001.5
56027.5 56048 56062 56068.2 56139.2 56175.2 56215.5 56251.9 56282.3 56306.3 56322.6 56329.1
56257 56304.8 56351.1 56391.7 56425 56451.4 56469.4 56475.5 55904.6 56008.9 56091.7 56158.1
56212.1 56255.6 56287 56294.8
</DataArray>
+ <DataArray type="Float32" Name="x^H2O_liquid" NumberOfComponents="1" format="ascii">
+ 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983
+ 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983
+ 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983
+ 0.999983 0.999983 0.999983 0.999983 0.999984 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983
+ 0.999984 0.999984 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.271856 0.328394 0.369365 0.401039
+ 0.427905 0.45113 0.46993 0.473456
+ </DataArray>
+ <DataArray type="Float32" Name="x^Air_liquid" NumberOfComponents="1" format="ascii">
+ 1.70317e-05 1.70355e-05 1.70426e-05 1.70516e-05 1.7061e-05 1.70695e-05 1.70758e-05 1.70792e-05 1.69279e-05 1.69374e-05 1.69539e-05 1.6974e-05
+ 1.69941e-05 1.70115e-05 1.70243e-05 1.70309e-05 1.67902e-05 1.6813e-05 1.68493e-05 1.68901e-05 1.69288e-05 1.69611e-05 1.69839e-05 1.69955e-05
+ 1.66491e-05 1.66942e-05 1.67581e-05 1.6824e-05 1.68833e-05 1.69312e-05 1.69644e-05 1.69804e-05 1.65207e-05 1.65971e-05 1.66907e-05 1.67795e-05
+ 1.68563e-05 1.69176e-05 1.69595e-05 1.6978e-05 1.64113e-05 1.65243e-05 1.66447e-05 1.67513e-05 1.68412e-05 1.69128e-05 1.69619e-05 1.69809e-05
+ 1.63226e-05 1.64727e-05 1.66142e-05 1.67336e-05 1.68328e-05 1.6912e-05 1.6967e-05 1.69851e-05 1.64541e-05 1.66199e-05 1.6763e-05 1.68811e-05
+ 1.69784e-05 1.70562e-05 1.71105e-05 1.71257e-05
+ </DataArray>
+ <DataArray type="Float32" Name="x^H2O_gas" NumberOfComponents="1" format="ascii">
+ 0.0149404 0.0149307 0.014913 0.0148906 0.014867 0.0148459 0.01483 0.0148215 0.015202 0.0151778 0.0151359 0.0150853
+ 0.0150346 0.0149907 0.0149587 0.014942 0.0155538 0.0154953 0.0154021 0.0152978 0.0151993 0.0151176 0.01506 0.0150308
+ 0.0159207 0.015803 0.0156367 0.0154668 0.0153151 0.0151931 0.0151093 0.0150689 0.0162862 0.0160577 0.015812 0.0155814
+ 0.015384 0.0152278 0.0151215 0.0150747 0.0166416 0.0162745 0.0159323 0.0156543 0.0154226 0.0152399 0.0151154 0.0150673
+ 0.0169333 0.0164414 0.0160124 0.0157002 0.0154442 0.0152418 0.0151026 0.0150568 0.00467152 0.00544081 0.00593606 0.00630953
+ 0.00661476 0.00687553 0.00709052 0.0071241
+ </DataArray>
+ <DataArray type="Float32" Name="x^Air_gas" NumberOfComponents="1" format="ascii">
+ 0.98506 0.985069 0.985087 0.985109 0.985133 0.985154 0.98517 0.985178 0.984798 0.984822 0.984864 0.984915
+ 0.984965 0.985009 0.985041 0.985058 0.984446 0.984505 0.984598 0.984702 0.984801 0.984882 0.98494 0.984969
+ 0.984079 0.984197 0.984363 0.984533 0.984685 0.984807 0.984891 0.984931 0.983714 0.983942 0.984188 0.984419
+ 0.984616 0.984772 0.984878 0.984925 0.983358 0.983725 0.984068 0.984346 0.984577 0.98476 0.984885 0.984933
+ 0.983067 0.983559 0.983988 0.9843 0.984556 0.984758 0.984897 0.984943 0.995328 0.994559 0.994064 0.99369
+ 0.993385 0.993124 0.992909 0.992876
+ </DataArray>
<DataArray type="Float32" Name="rhoMolar_liquid" NumberOfComponents="1" format="ascii">
55474.1 55474.1 55474.2 55474.4 55474.6 55474.7 55474.8 55474.9 55472.2 55472.4 55472.7 55473
55473.4 55473.7 55473.9 55474 55469.8 55470.2 55470.8 55471.5 55472.2 55472.8 55473.2 55473.4
- 55467.2 55468 55469.2 55470.4 55471.4 55472.3 55472.8 55473.1 55464.5 55466.2 55467.9 55469.6
+ 55467.2 55468 55469.2 55470.4 55471.4 55472.3 55472.8 55473.1 55464.5 55466.2 55467.9 55469.5
55470.9 55472 55472.7 55473.1 55461.8 55464.6 55467.1 55469 55470.6 55471.9 55472.8 55473.1
55459.5 55463.3 55466.5 55468.7 55470.5 55471.9 55472.9 55473.2 55457.6 55462.3 55466.1 55468.5
55470.4 55471.9 55472.9 55473.2
</DataArray>
<DataArray type="Float32" Name="rhoMolar_gas" NumberOfComponents="1" format="ascii">
- 42.0507 42.0522 42.055 42.0584 42.0621 42.0654 42.0678 42.0692 42.0098 42.0136 42.0201 42.0279
- 42.0357 42.0426 42.0475 42.0501 41.9552 41.9642 41.9786 41.9947 42.01 42.0226 42.0315 42.036
- 41.8984 41.9166 41.9422 41.9685 41.9919 42.0108 42.0237 42.0299 41.8463 41.8772 41.9151 41.9507
- 41.9812 42.0053 42.0217 42.0289 41.8016 41.8477 41.8965 41.9394 41.9752 42.0034 42.0225 42.0299
- 41.765 41.8267 41.8841 41.9323 41.9718 42.0031 42.0244 42.0314 41.7329 41.8101 41.8745 41.9267
- 41.9691 42.0027 42.0261 42.0324
+ 42.0507 42.0522 42.055 42.0584 42.0621 42.0654 42.0679 42.0692 42.0098 42.0136 42.0201 42.0279
+ 42.0358 42.0426 42.0475 42.0501 41.9551 41.9642 41.9786 41.9947 42.01 42.0227 42.0316 42.0361
+ 41.8984 41.9165 41.9422 41.9684 41.9919 42.0108 42.0238 42.0301 41.8462 41.8772 41.915 41.9506
+ 41.9811 42.0053 42.0218 42.0291 41.8016 41.8476 41.8964 41.9393 41.9751 42.0034 42.0227 42.0301
+ 41.765 41.8266 41.884 41.9322 41.9717 42.003 42.0246 42.0317 41.7329 41.8101 41.8744 41.9265
+ 41.969 42.0027 42.0263 42.0327
</DataArray>
<DataArray type="Float32" Name="X^H2O_liquid" NumberOfComponents="1" format="ascii">
0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973
0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973
0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973
0.999973 0.999973 0.999973 0.999973 0.999974 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973
- 0.999974 0.999974 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.271845 0.328384 0.369355 0.401028
+ 0.999974 0.999974 0.999973 0.999973 0.999973 0.999973 0.999973 0.999973 0.271846 0.328384 0.369355 0.401028
0.427894 0.45112 0.46992 0.473446
</DataArray>
<DataArray type="Float32" Name="X^Air_liquid" NumberOfComponents="1" format="ascii">
@@ -108,12 +148,12 @@
2.72927e-05 2.74178e-05 2.7505e-05 2.75296e-05
</DataArray>
<DataArray type="Float32" Name="X^H2O_gas" NumberOfComponents="1" format="ascii">
- 0.0093468 0.00934066 0.00932955 0.00931544 0.00930062 0.00928729 0.00927731 0.00927198 0.00951134 0.00949612 0.00946977 0.00943795
- 0.00940603 0.00937839 0.00935828 0.00934775 0.00973279 0.00969598 0.0096373 0.00957168 0.00950969 0.00945827 0.00942204 0.00940366
- 0.00996376 0.00988963 0.00978498 0.00967805 0.00958255 0.0095058 0.00945301 0.0094276 0.0101939 0.01005 0.00989533 0.00975018
+ 0.00934679 0.00934066 0.00932955 0.00931544 0.00930062 0.00928729 0.00927731 0.00927198 0.00951134 0.00949612 0.00946977 0.00943795
+ 0.00940603 0.00937839 0.00935828 0.00934775 0.00973279 0.00969598 0.00963729 0.00957168 0.00950969 0.00945827 0.00942204 0.00940366
+ 0.00996376 0.00988963 0.00978498 0.00967804 0.00958255 0.0095058 0.00945301 0.0094276 0.0101939 0.01005 0.00989533 0.00975018
0.00962593 0.00952758 0.00946072 0.00943127 0.0104178 0.0101865 0.00997109 0.00979604 0.0096502 0.00953519 0.00945687 0.00942661
- 0.0106016 0.0102917 0.0100215 0.00982494 0.00966379 0.00953644 0.00944881 0.00942 0.00291115 0.00339155 0.00370095 0.00393436
- 0.00412517 0.00428821 0.00442266 0.00444366
+ 0.0106016 0.0102917 0.0100215 0.00982494 0.00966379 0.00953644 0.00944881 0.00942 0.00291116 0.00339155 0.00370095 0.00393436
+ 0.00412516 0.00428821 0.00442266 0.00444366
</DataArray>
<DataArray type="Float32" Name="X^Air_gas" NumberOfComponents="1" format="ascii">
0.990653 0.990659 0.99067 0.990685 0.990699 0.990713 0.990723 0.990728 0.990489 0.990504 0.99053 0.990562
@@ -123,47 +163,7 @@
0.989398 0.989708 0.989978 0.990175 0.990336 0.990464 0.990551 0.99058 0.997089 0.996608 0.996299 0.996066
0.995875 0.995712 0.995577 0.995556
</DataArray>
- <DataArray type="Float32" Name="x^H2O_liquid" NumberOfComponents="1" format="ascii">
- 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983
- 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983
- 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983
- 0.999983 0.999983 0.999983 0.999983 0.999984 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983
- 0.999984 0.999984 0.999983 0.999983 0.999983 0.999983 0.999983 0.999983 0.271855 0.328394 0.369365 0.401039
- 0.427905 0.45113 0.469931 0.473456
- </DataArray>
- <DataArray type="Float32" Name="x^Air_liquid" NumberOfComponents="1" format="ascii">
- 1.70317e-05 1.70355e-05 1.70426e-05 1.70516e-05 1.7061e-05 1.70695e-05 1.70758e-05 1.70792e-05 1.69279e-05 1.69374e-05 1.69539e-05 1.6974e-05
- 1.69941e-05 1.70115e-05 1.70243e-05 1.70309e-05 1.67902e-05 1.6813e-05 1.68493e-05 1.68901e-05 1.69288e-05 1.69611e-05 1.69839e-05 1.69955e-05
- 1.66491e-05 1.66941e-05 1.67581e-05 1.6824e-05 1.68833e-05 1.69312e-05 1.69644e-05 1.69804e-05 1.65207e-05 1.65971e-05 1.66907e-05 1.67795e-05
- 1.68563e-05 1.69176e-05 1.69595e-05 1.6978e-05 1.64113e-05 1.65243e-05 1.66447e-05 1.67513e-05 1.68412e-05 1.69128e-05 1.69619e-05 1.69809e-05
- 1.63226e-05 1.64727e-05 1.66142e-05 1.67336e-05 1.68328e-05 1.6912e-05 1.6967e-05 1.69851e-05 1.64541e-05 1.66199e-05 1.6763e-05 1.68811e-05
- 1.69784e-05 1.70562e-05 1.71105e-05 1.71257e-05
- </DataArray>
- <DataArray type="Float32" Name="x^H2O_gas" NumberOfComponents="1" format="ascii">
- 0.0149404 0.0149307 0.014913 0.0148906 0.014867 0.0148459 0.01483 0.0148215 0.015202 0.0151778 0.0151359 0.0150853
- 0.0150346 0.0149907 0.0149587 0.014942 0.0155538 0.0154953 0.0154021 0.0152978 0.0151993 0.0151176 0.01506 0.0150308
- 0.0159207 0.015803 0.0156367 0.0154668 0.0153151 0.0151931 0.0151093 0.0150689 0.0162862 0.0160577 0.015812 0.0155814
- 0.015384 0.0152278 0.0151215 0.0150747 0.0166416 0.0162745 0.0159323 0.0156543 0.0154226 0.0152399 0.0151154 0.0150673
- 0.0169333 0.0164414 0.0160124 0.0157002 0.0154442 0.0152418 0.0151026 0.0150568 0.00467151 0.00544081 0.00593606 0.00630953
- 0.00661476 0.00687553 0.00709052 0.0071241
- </DataArray>
- <DataArray type="Float32" Name="x^Air_gas" NumberOfComponents="1" format="ascii">
- 0.98506 0.985069 0.985087 0.985109 0.985133 0.985154 0.98517 0.985178 0.984798 0.984822 0.984864 0.984915
- 0.984965 0.985009 0.985041 0.985058 0.984446 0.984505 0.984598 0.984702 0.984801 0.984882 0.98494 0.984969
- 0.984079 0.984197 0.984363 0.984533 0.984685 0.984807 0.984891 0.984931 0.983714 0.983942 0.984188 0.984419
- 0.984616 0.984772 0.984878 0.984925 0.983358 0.983725 0.984068 0.984346 0.984577 0.98476 0.984885 0.984933
- 0.983067 0.983559 0.983988 0.9843 0.984556 0.984758 0.984897 0.984943 0.995328 0.994559 0.994064 0.99369
- 0.993385 0.993124 0.992909 0.992876
- </DataArray>
- <DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii">
- 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41
- 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41
- 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41
- 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41
- 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41
- 0.41 0.41 0.41 0.41
- </DataArray>
- <DataArray type="Float32" Name="phase presence" NumberOfComponents="1" format="ascii">
+ <DataArray type="Float32" Name="phasePresence" NumberOfComponents="1" format="ascii">
3 3 3 3 3 3 3 3 3 3 3 3
3 3 3 3 3 3 3 3 3 3 3 3
3 3 3 3 3 3 3 3 3 3 3 3
@@ -180,40 +180,40 @@
286.575 286.345 286.184 286.14
</DataArray>
<DataArray type="Float32" Name="velocity_liquid (m/s)" NumberOfComponents="3" format="ascii">
- -3.2019e-10 1.01795e-08 0 -4.27423e-10 9.88581e-09 0 -5.83728e-10 9.56262e-09 0 -6.3177e-10 9.28427e-09 0
- -5.88822e-10 9.0572e-09 0 -4.73142e-10 8.88043e-09 0 -3.04001e-10 8.76045e-09 0 -2.08623e-10 8.71084e-09 0
- 4.81065e-11 1.17697e-08 0 6.95072e-11 1.14803e-08 0 1.03203e-10 1.11729e-08 0 1.17437e-10 1.09126e-08 0
- 1.12316e-10 1.07011e-08 0 9.14173e-11 1.0536e-08 0 5.94309e-11 1.04234e-08 0 4.12821e-11 1.03781e-08 0
- -1.11461e-11 1.36937e-08 0 -8.97751e-12 1.33979e-08 0 -4.21085e-12 1.31044e-08 0 -2.44955e-13 1.28621e-08 0
- 1.48694e-12 1.26662e-08 0 1.76467e-12 1.25123e-08 0 1.70286e-12 1.24063e-08 0 1.72741e-12 1.23661e-08 0
- -4.16693e-11 1.41508e-08 0 -3.93484e-11 1.38211e-08 0 -3.24801e-11 1.35225e-08 0 -2.46601e-11 1.32804e-08 0
- -1.89465e-11 1.30844e-08 0 -1.39287e-11 1.29293e-08 0 -7.43816e-12 1.2821e-08 0 -3.52464e-12 1.27832e-08 0
- -6.43261e-11 1.43425e-08 0 -5.75798e-11 1.39652e-08 0 -4.43717e-11 1.3657e-08 0 -3.37402e-11 1.34154e-08 0
+ -3.2019e-10 1.01795e-08 0 -4.27422e-10 9.8858e-09 0 -5.83727e-10 9.56262e-09 0 -6.31769e-10 9.28427e-09 0
+ -5.88821e-10 9.0572e-09 0 -4.73141e-10 8.88043e-09 0 -3.04001e-10 8.76045e-09 0 -2.08623e-10 8.71084e-09 0
+ 4.81063e-11 1.17697e-08 0 6.95069e-11 1.14803e-08 0 1.03202e-10 1.11729e-08 0 1.17437e-10 1.09126e-08 0
+ 1.12315e-10 1.07011e-08 0 9.14168e-11 1.0536e-08 0 5.94306e-11 1.04234e-08 0 4.12819e-11 1.03781e-08 0
+ -1.1146e-11 1.36937e-08 0 -8.97748e-12 1.33979e-08 0 -4.21086e-12 1.31044e-08 0 -2.44983e-13 1.28621e-08 0
+ 1.48692e-12 1.26662e-08 0 1.76465e-12 1.25123e-08 0 1.70284e-12 1.24063e-08 0 1.72739e-12 1.23661e-08 0
+ -4.16691e-11 1.41508e-08 0 -3.93482e-11 1.38211e-08 0 -3.248e-11 1.35225e-08 0 -2.466e-11 1.32804e-08 0
+ -1.89465e-11 1.30844e-08 0 -1.39286e-11 1.29293e-08 0 -7.43815e-12 1.2821e-08 0 -3.52464e-12 1.27832e-08 0
+ -6.43257e-11 1.43425e-08 0 -5.75796e-11 1.39652e-08 0 -4.43715e-11 1.3657e-08 0 -3.37401e-11 1.34154e-08 0
-2.64787e-11 1.32193e-08 0 -2.00133e-11 1.30629e-08 0 -1.06096e-11 1.29513e-08 0 -4.57925e-12 1.29171e-08 0
- -6.91593e-11 1.44474e-08 0 -5.92703e-11 1.40313e-08 0 -4.24443e-11 1.37119e-08 0 -3.18076e-11 1.34704e-08 0
- -2.53837e-11 1.32742e-08 0 -1.96819e-11 1.31167e-08 0 -1.02857e-11 1.30019e-08 0 -3.86696e-12 1.29712e-08 0
- -2.75129e-11 9.23357e-09 0 -2.37367e-11 9.38685e-09 0 -1.72203e-11 9.48488e-09 0 -1.32237e-11 9.54575e-09 0
- -1.09944e-11 9.59605e-09 0 -8.87965e-12 9.64117e-09 0 -4.65884e-12 9.67951e-09 0 -1.57991e-12 9.68324e-09 0
- 0 3.98326e-09 0 0 4.72035e-09 0 0 5.24048e-09 0 0 5.60397e-09 0
- 0 5.90103e-09 0 0 6.14917e-09 0 0 6.3423e-09 0 0 6.37883e-09 0
+ -6.91588e-11 1.44474e-08 0 -5.927e-11 1.40313e-08 0 -4.24441e-11 1.37119e-08 0 -3.18076e-11 1.34704e-08 0
+ -2.53836e-11 1.32742e-08 0 -1.96819e-11 1.31167e-08 0 -1.02857e-11 1.30019e-08 0 -3.86697e-12 1.29712e-08 0
+ -2.75127e-11 9.23358e-09 0 -2.37365e-11 9.38685e-09 0 -1.72203e-11 9.48488e-09 0 -1.32237e-11 9.54575e-09 0
+ -1.09944e-11 9.59605e-09 0 -8.87962e-12 9.64117e-09 0 -4.65883e-12 9.67951e-09 0 -1.57991e-12 9.68324e-09 0
+ 0 3.98326e-09 0 0 4.72036e-09 0 0 5.24048e-09 0 0 5.60397e-09 0
+ 0 5.90103e-09 0 0 6.14917e-09 0 0 6.3423e-09 0 0 6.37882e-09 0
</DataArray>
<DataArray type="Float32" Name="velocity_gas (m/s)" NumberOfComponents="3" format="ascii">
- -4.2836e-09 -6.38625e-08 0 -6.33989e-09 -6.58319e-08 0 -9.49953e-09 -6.43835e-08 0 -1.07333e-08 -6.12961e-08 0
- -1.01596e-08 -5.79456e-08 0 -8.16051e-09 -5.50511e-08 0 -5.22289e-09 -5.297e-08 0 -3.58014e-09 -5.18524e-08 0
- -3.61274e-08 -1.08209e-08 0 -4.94709e-08 -1.99421e-08 0 -7.02719e-08 -3.1462e-08 0 -7.86167e-08 -4.50615e-08 0
- -7.42346e-08 -5.54879e-08 0 -5.93708e-08 -6.24027e-08 0 -3.7752e-08 -6.59554e-08 0 -2.57276e-08 -6.69315e-08 0
- -3.53123e-08 5.20477e-08 0 -4.68313e-08 2.95261e-08 0 -6.22718e-08 -2.4814e-09 0 -6.37441e-08 -3.70797e-08 0
- -5.49079e-08 -6.48175e-08 0 -3.98291e-08 -8.33114e-08 0 -2.24091e-08 -9.17955e-08 0 -1.36809e-08 -9.21732e-08 0
- -4.30799e-09 5.25882e-08 0 -6.91296e-09 2.96312e-08 0 -9.73744e-09 -1.04007e-08 0 -7.7358e-09 -5.00614e-08 0
- -1.56752e-09 -8.17479e-08 0 6.61684e-09 -1.01494e-07 0 1.21213e-08 -1.06343e-07 0 1.33886e-08 -9.72926e-08 0
- 6.20522e-08 1.15377e-08 0 6.24804e-08 1.46681e-08 0 6.15476e-08 -1.61043e-08 0 5.90905e-08 -5.57449e-08 0
- 5.77688e-08 -8.76408e-08 0 5.69774e-08 -1.06166e-07 0 5.11145e-08 -1.05147e-07 0 4.58173e-08 -7.9629e-08 0
- 1.56003e-07 -5.57782e-08 0 1.50879e-07 -5.17583e-09 0 1.35837e-07 -1.72799e-08 0 1.20042e-07 -5.57459e-08 0
- 1.09776e-07 -8.73602e-08 0 1.00922e-07 -1.04675e-07 0 8.7111e-08 -9.85917e-08 0 7.77637e-08 -5.25539e-08 0
- 2.57838e-07 6.87451e-06 0 2.39485e-07 6.19662e-06 0 2.00135e-07 5.62596e-06 0 1.68717e-07 5.17886e-06 0
- 1.50564e-07 4.81557e-06 0 1.35527e-07 4.52671e-06 0 1.16915e-07 4.32847e-06 0 1.05609e-07 4.35204e-06 0
- 2.21039e-08 1.38404e-05 0 3.09105e-08 1.24068e-05 0 4.29492e-08 1.12687e-05 0 4.84241e-08 1.04125e-05 0
- 5.3041e-08 9.71748e-06 0 5.97901e-08 9.15657e-06 0 8.65794e-08 8.75187e-06 0 1.08993e-07 8.74217e-06 0
+ -4.28519e-09 -6.37201e-08 0 -6.3454e-09 -6.5726e-08 0 -9.51115e-09 -6.43013e-08 0 -1.07478e-08 -6.12287e-08 0
+ -1.0174e-08 -5.78869e-08 0 -8.17248e-09 -5.49964e-08 0 -5.23064e-09 -5.29172e-08 0 -3.58546e-09 -5.18005e-08 0
+ -3.62312e-08 -1.06515e-08 0 -4.96079e-08 -1.98164e-08 0 -7.04587e-08 -3.13837e-08 0 -7.88211e-08 -4.50224e-08 0
+ -7.44261e-08 -5.54844e-08 0 -5.9524e-08 -6.24231e-08 0 -3.785e-08 -6.59901e-08 0 -2.57948e-08 -6.69722e-08 0
+ -3.55456e-08 5.23883e-08 0 -4.71267e-08 2.9748e-08 0 -6.26458e-08 -2.38662e-09 0 -6.41189e-08 -3.70894e-08 0
+ -5.52348e-08 -6.4911e-08 0 -4.00793e-08 -8.34638e-08 0 -2.25636e-08 -9.19839e-08 0 -1.3784e-08 -9.23733e-08 0
+ -4.75437e-09 5.3257e-08 0 -7.42588e-09 2.99846e-08 0 -1.03071e-08 -1.02938e-08 0 -8.25842e-09 -5.01222e-08 0
+ -2.00321e-09 -8.19241e-08 0 6.29231e-09 -1.01754e-07 0 1.19294e-08 -1.06656e-07 0 1.32679e-08 -9.76088e-08 0
+ 6.12922e-08 1.26006e-08 0 6.16994e-08 1.50889e-08 0 6.07981e-08 -1.60338e-08 0 5.84532e-08 -5.58689e-08 0
+ 5.72528e-08 -8.78916e-08 0 5.65958e-08 -1.06517e-07 0 5.09012e-08 -1.05571e-07 0 4.56992e-08 -8.0027e-08 0
+ 1.54857e-07 -5.42952e-08 0 1.49813e-07 -4.77175e-09 0 1.34949e-07 -1.72732e-08 0 1.19328e-07 -5.59301e-08 0
+ 1.09204e-07 -8.7669e-08 0 1.00495e-07 -1.05093e-07 0 8.68874e-08 -9.91135e-08 0 7.76658e-08 -5.30022e-08 0
+ 2.56288e-07 6.87638e-06 0 2.38155e-07 6.19694e-06 0 1.99155e-07 5.6259e-06 0 1.67954e-07 5.17863e-06 0
+ 1.49955e-07 4.81522e-06 0 1.35063e-07 4.52625e-06 0 1.16688e-07 4.32787e-06 0 1.05541e-07 4.35156e-06 0
+ 2.02008e-08 1.38424e-05 0 2.93705e-08 1.2407e-05 0 4.1922e-08 1.12686e-05 0 4.76368e-08 1.04123e-05 0
+ 5.24111e-08 9.71711e-06 0 5.92997e-08 9.15608e-06 0 8.63532e-08 8.75122e-06 0 1.08959e-07 8.74169e-06 0
</DataArray>
<DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
0 0 0 0 0 0 0 0 0 0 0 0