Commit 0ccde397 authored by Gabi Seitz's avatar Gabi Seitz Committed by Timo Koch
Browse files

[vapor][fluidsystems] add method gasViscosityForMixtures

with contributions of Holger and Bernd

The following test results had to be adapted:
test_md_boundary_darcy2p2c_stokes1p2c_horizontal
test_md_boundary_darcy2p2cni_stokes1p2cni_horizontal
test_2p1cni_steaminjection_waterwet_box
test_2p1cni_steaminjection_waterwet_tpfa
test_2p2cni_waterair_box
test_2p2cni_waterair_buoyancy_box
test_dec2p2c
test_multiphysics2p2c
test_1pncmin_box
test_3p3cni_kuevette_box
test_3p3cni_kuevette_tpfa
test_3pwateroil_sagd_box
parent ac3b2fa4
......@@ -708,12 +708,14 @@ public:
* \param temperature temperature of component in \f$\mathrm{[K]}\f$
* \param pressure pressure of component in \f$\mathrm{[Pa]}\f$
*
* This method is only valid if pressure is below or at the vapor
* pressure of water.
* We assume pure water vapor here. For water in a mixture of other gaseous
* components, consider the free function h2oGasViscosityInMixture.
*
* See:
* We use the IAPWS Formulation, see:
* IAPWS: "Release on the IAPWS Formulation 2008 for the Viscosity
* of Ordinary Water Substance", http://www.iapws.org/relguide/visc.pdf \cite cooper2008
* This method is only valid if pressure is below or at the vapor
* pressure of water.
*/
static Scalar gasViscosity(Scalar temperature, Scalar pressure)
{
......@@ -723,6 +725,7 @@ public:
return Common::viscosity(temperature, rho);
}
/*!
* \brief The dynamic viscosity \f$\mathrm{[Pa*s]}\f$ of pure water.
*
......@@ -905,13 +908,88 @@ private:
Region2::dGamma_dPi(temperature, pressure) *
Rs * temperature / pressure;
}
}; // end class
};
template <class Scalar>
struct IsAqueous<H2O<Scalar>> : public std::true_type {};
} // end namespace Components
namespace FluidSystems::Detail {
// viscosity according to Reid, R.C.
template <class Scalar>
Scalar viscosityReid_(Scalar temperature)
{
constexpr Scalar tc = 647.3;
using std::max;
const Scalar tr = max(temperature/tc, 1e-8);
const Scalar fp0 = 1.0 + 0.221*(0.96 + 0.1*(tr - 0.7));
constexpr Scalar xi = 3.334e-3;
using std::pow;
using std::exp;
const Scalar eta_xi = (0.807*pow(tr, 0.618) - 0.357*exp((-0.449)*tr)
+ 0.34*exp((-4.058)*tr) + 0.018)*fp0;
return 1.0e-7*eta_xi/xi;
}
// viscosity according to Nagel, T. et al.
template <class Scalar>
Scalar viscosityNagel_(Scalar temperature)
{
constexpr Scalar a1 = -4.4189440e-6;
constexpr Scalar a2 = 4.6876380e-8;
constexpr Scalar a3 = -5.3894310e-12;
constexpr Scalar a4 = 3.2028560e-16;
constexpr Scalar a5 = 4.9191790e-22;
return a1 + a2*temperature + a3*temperature*temperature
+ a4*temperature*temperature*temperature
+ a5*temperature*temperature*temperature*temperature;
}
} // end FluidSystems::Detail
namespace FluidSystems {
/*!
* \brief The dynamic viscosity \f$\mathrm{[Pa*s]}\f$ of steam in a gas mixture.
*
* \param temperature temperature in \f$\mathrm{[K]}\f$
*
* We assume here that water is in mixture with other gaseous components.
* For pure water, use the gasViscosity function of Components::H2O.
*
* We apply two different laws depending on the gas temperature.
*
* For temperatures below 480 K see:
* "Reid, R.C., Prausnitz, J.M., Poling, B.E.: The Properties of
* Gases and Liquids (1987)"
* Lucas corresponding states method
* https://www.osti.gov/scitech/biblio/6504847} \cite{reid1987}
*
* For temperatures above 500 K see:
* Nagel, T. et al.: THC-Processes (2018)
* https://doi.org/10.1007/978-3-319-68225-9_12
*
* In the range 480 - 500 K, we interpolate between the two laws.
*/
template <class Scalar>
Scalar h2oGasViscosityInMixture(Scalar temperature, Scalar pressure)
{
if (temperature < 480.0)
return Detail::viscosityReid_(temperature);
else if (temperature > 500.0)
return Detail::viscosityNagel_(temperature);
else // interpolate
{
const Scalar op = (500.0 - temperature)/20.0;
return op*Detail::viscosityReid_(temperature)
+ (1.0 - op)*Detail::viscosityNagel_(temperature);
}
}
} // end namespace FluidSystems
} // end namespace Dumux
#endif
......@@ -481,8 +481,7 @@ public:
// Wilke method (Reid et al.):
Scalar muResult = 0;
const Scalar mu[numComponents] = {
H2O::gasViscosity(T,
H2O::vaporPressure(T)),
h2oGasViscosityInMixture(T, p),
Air::gasViscosity(T, p)
};
......
......@@ -349,7 +349,7 @@ public:
*/
Scalar muResult;
const Scalar mu[numComponents] = {
H2O::gasViscosity(fluidState.temperature(phaseIdx), H2O::vaporPressure(fluidState.temperature(phaseIdx))),
h2oGasViscosityInMixture(fluidState.temperature(phaseIdx), fluidState.pressure(phaseIdx)),
Air::gasViscosity(fluidState.temperature(phaseIdx), fluidState.pressure(phaseIdx)),
NAPL::gasViscosity(fluidState.temperature(phaseIdx), NAPL::vaporPressure(fluidState.temperature(phaseIdx)))
};
......
......@@ -348,7 +348,7 @@ public:
*/
Scalar muResult;
const Scalar mu[numComponents] = {
H2O::gasViscosity(fluidState.temperature(phaseIdx), H2O::vaporPressure(fluidState.temperature(phaseIdx))),
h2oGasViscosityInMixture(fluidState.temperature(phaseIdx),fluidState.pressure(phaseIdx)),
Air::gasViscosity(fluidState.temperature(phaseIdx), fluidState.pressure(phaseIdx)),
NAPL::gasViscosity(fluidState.temperature(phaseIdx), NAPL::vaporPressure(fluidState.temperature(phaseIdx)))
};
......
......@@ -337,7 +337,7 @@ public:
* -- compare e.g. with Promo Class p. 32/33
*/
const Scalar mu[numComponents] = {
H2O::gasViscosity(fluidState.temperature(phaseIdx), H2O::vaporPressure(fluidState.temperature(phaseIdx))),
h2oGasViscosityInMixture(fluidState.temperature(phaseIdx), fluidState.pressure(phaseIdx)),
HeavyOil::gasViscosity(fluidState.temperature(phaseIdx), HeavyOil::vaporPressure(fluidState.temperature(phaseIdx)))
};
......
......@@ -489,7 +489,7 @@ public:
// Wilke method (Reid et al.):
Scalar muResult = 0;
const Scalar mu[numComponents] = {
H2O::gasViscosity(T, H2O::vaporPressure(T)),
h2oGasViscosityInMixture(T, p),
N2::gasViscosity(T, p)
};
......
......@@ -563,7 +563,7 @@ public:
// Wilke method (Reid et al.):
Scalar muResult = 0;
const Scalar mu[numComponents] = {
H2O::gasViscosity(T, H2O::vaporPressure(T)),
h2oGasViscosityInMixture(T, p),
N2::gasViscosity(T, p),
O2::gasViscosity(T, p)
};
......
......@@ -7,10 +7,10 @@
200000 200000 200000 200000 200000 200000
</DataArray>
<DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
0.981775 0.981775 0.981775 0.981775 0.981775 0.981775
0.981785 0.981785 0.981785 0.981785 0.981785 0.981785
</DataArray>
<DataArray type="Float32" Name="mu" NumberOfComponents="1" format="ascii">
2.44123e-05 2.44123e-05 2.44123e-05 2.44123e-05 2.44123e-05 2.44123e-05
2.50048e-05 2.50048e-05 2.50048e-05 2.50048e-05 2.50048e-05 2.50048e-05
</DataArray>
<DataArray type="Float32" Name="delp" NumberOfComponents="1" format="ascii">
100000 100000 100000 100000 100000 100000
......@@ -27,11 +27,11 @@
<DataArray type="Float32" Name="X^H2O_gas" NumberOfComponents="1" format="ascii">
0.357618 0.357618 0.357618 0.357618 0.357618 0.357618
</DataArray>
<DataArray type="Float32" Name="precipVolFrac_CaO" NumberOfComponents="1" format="ascii">
0.105174 0.105174 0.105174 0.105174 0.105174 0.105174
<DataArray type="Float32" Name="precipitateVolumeFraction^CaO" NumberOfComponents="1" format="ascii">
0.105328 0.105328 0.105328 0.105328 0.105328 0.105328
</DataArray>
<DataArray type="Float32" Name="precipVolFrac_CaOH2" NumberOfComponents="1" format="ascii">
0.0943089 0.0943089 0.0943089 0.0943089 0.0943089 0.0943089
<DataArray type="Float32" Name="precipitateVolumeFraction^CaO2H2" NumberOfComponents="1" format="ascii">
0.0941565 0.0941565 0.0941565 0.0941565 0.0941565 0.0941565
</DataArray>
<DataArray type="Float32" Name="T" NumberOfComponents="1" format="ascii">
573.15 573.15 573.15 573.15 573.15 573.15
......@@ -40,10 +40,10 @@
8.53e-12 8.53e-12 8.53e-12 8.53e-12 8.53e-12 8.53e-12
</DataArray>
<DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii">
0.8 0.8 0.8 0.8 0.8 0.8
0.800516 0.800516 0.800516 0.800516 0.800516 0.800516
</DataArray>
<DataArray type="Float32" Name="reactionRate" NumberOfComponents="1" format="ascii">
7.84007 7.84007 7.84007 7.84007 7.84007 7.84007
7.8396 7.8396 7.8396 7.8396 7.8396 7.8396
</DataArray>
</PointData>
<CellData Scalars="process rank">
......
This source diff could not be displayed because it is too large. You can view the blob instead.
......@@ -4,18 +4,18 @@
<Piece NumberOfCells="64" NumberOfPoints="81">
<CellData Scalars="S_liq">
<DataArray type="Float32" Name="S_liq" NumberOfComponents="1" format="ascii">
0.614562 0.614562 0.614562 0.614562 0.614563 0.614563 0.614563 0.614563 0.125427 0.125428 0.125428 0.125428
0.125428 0.125428 0.125428 0.125428 0.0457042 0.0457048 0.0457055 0.045706 0.0457064 0.0457068 0.0457071 0.0457074
0.0257656 0.0257705 0.025775 0.0257786 0.0257815 0.0257841 0.0257863 0.025788 0.0172624 0.0172791 0.0172928 0.0173036
0.0173125 0.0173202 0.017327 0.0173325 0.0108362 0.010897 0.010943 0.0109788 0.0110083 0.0110337 0.0110561 0.0110753
0.614566 0.614566 0.614567 0.614567 0.614567 0.614567 0.614567 0.614567 0.125429 0.125429 0.125429 0.125429
0.125429 0.125429 0.125429 0.125429 0.0457046 0.0457052 0.0457058 0.0457064 0.0457068 0.0457072 0.0457075 0.0457077
0.0257658 0.0257707 0.0257751 0.0257787 0.0257817 0.0257843 0.0257865 0.0257882 0.0172625 0.0172792 0.0172929 0.0173037
0.0173127 0.0173203 0.0173271 0.0173326 0.0108364 0.0108971 0.0109431 0.0109789 0.0110085 0.0110339 0.0110563 0.0110754
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0
</DataArray>
<DataArray type="Float32" Name="p_liq" NumberOfComponents="1" format="ascii">
98496 98496 98496 98496 98496 98496 98496 98496 97789.1 97789.1 97789.1 97789.1
97789.1 97789.1 97789.1 97789.1 97316.9 97316.9 97316.9 97316.9 97316.9 97316.9 97316.9 97316.9
96986.3 96986.4 96986.6 96986.6 96986.7 96986.8 96986.8 96986.9 96702.2 96702.9 96703.5 96704
96704.4 96704.8 96705.1 96705.4 96339 96342.6 96345.3 96347.5 96349.2 96350.7 96352 96353.2
96986.3 96986.5 96986.6 96986.6 96986.7 96986.8 96986.8 96986.9 96702.2 96702.9 96703.6 96704
96704.5 96704.8 96705.1 96705.4 96339 96342.6 96345.3 96347.5 96349.2 96350.7 96352 96353.2
95697.6 95697.6 95697.6 95697.6 95697.6 95697.6 95697.6 95697.6 95697.5 95697.5 95697.5 95697.5
95697.5 95697.5 95697.5 95697.5
</DataArray>
......@@ -24,22 +24,22 @@
997.051 997.051 997.051 997.051 997.051 997.051 997.051 997.051 997.051 997.051 997.051 997.051
997.051 997.051 997.051 997.051 997.051 997.051 997.051 997.051 997.05 997.05 997.05 997.05
997.05 997.05 997.05 997.05 997.05 997.05 997.05 997.05 997.05 997.05 997.05 997.05
610.336 628.473 641.854 652.355 661.125 668.727 675.495 681.465 136.094 159.997 176.435 189.323
200.097 209.437 217.763 225.419
610.348 628.485 641.867 652.369 661.139 668.743 675.511 681.482 136.098 160.002 176.441 189.33
200.105 209.446 217.772 225.429
</DataArray>
<DataArray type="Float32" Name="mob_liq" NumberOfComponents="1" format="ascii">
233.927 233.927 233.927 233.927 233.927 233.928 233.928 233.928 2.13229 2.1323 2.1323 2.13231
2.13231 2.13231 2.13232 2.13232 0.0971432 0.0971476 0.0971519 0.0971555 0.0971585 0.097161 0.0971632 0.0971648
0.0143468 0.0143565 0.0143653 0.0143723 0.0143782 0.0143832 0.0143875 0.0143908 0.00321343 0.00322583 0.00323607 0.00324416
0.00325087 0.00325665 0.00326172 0.00326585 0.000390248 0.000401898 0.000410863 0.000417924 0.000423818 0.000428925 0.000433459 0.000437358
233.932 233.932 233.932 233.932 233.932 233.933 233.933 233.933 2.13236 2.13236 2.13236 2.13237
2.13237 2.13238 2.13238 2.13238 0.0971458 0.0971502 0.0971545 0.097158 0.097161 0.0971636 0.0971658 0.0971674
0.0143472 0.0143568 0.0143656 0.0143727 0.0143785 0.0143835 0.0143879 0.0143912 0.00321351 0.00322591 0.00323616 0.00324425
0.00325096 0.00325673 0.0032618 0.00326593 0.000390271 0.000401922 0.000410887 0.000417949 0.000423843 0.00042895 0.000433485 0.000437385
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0
</DataArray>
<DataArray type="Float32" Name="S_gas" NumberOfComponents="1" format="ascii">
0.385438 0.385438 0.385438 0.385438 0.385437 0.385437 0.385437 0.385437 0.874573 0.874572 0.874572 0.874572
0.874572 0.874572 0.874572 0.874572 0.954296 0.954295 0.954295 0.954294 0.954294 0.954293 0.954293 0.954293
0.974234 0.97423 0.974225 0.974221 0.974218 0.974216 0.974214 0.974212 0.982738 0.982721 0.982707 0.982696
0.982687 0.98268 0.982673 0.982668 0.989164 0.989103 0.989057 0.989021 0.988992 0.988966 0.988944 0.988925
0.385434 0.385434 0.385433 0.385433 0.385433 0.385433 0.385433 0.385433 0.874571 0.874571 0.874571 0.874571
0.874571 0.874571 0.874571 0.874571 0.954295 0.954295 0.954294 0.954294 0.954293 0.954293 0.954292 0.954292
0.974234 0.974229 0.974225 0.974221 0.974218 0.974216 0.974214 0.974212 0.982737 0.982721 0.982707 0.982696
0.982687 0.98268 0.982673 0.982667 0.989164 0.989103 0.989057 0.989021 0.988992 0.988966 0.988944 0.988925
1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1
</DataArray>
......@@ -60,18 +60,18 @@
1.16542 1.16529 1.16517 1.16507
</DataArray>
<DataArray type="Float32" Name="mob_gas" NumberOfComponents="1" format="ascii">
8028.89 8028.88 8028.88 8028.87 8028.86 8028.85 8028.85 8028.85 44360.2 44360.2 44360.2 44360.2
44360.2 44360.2 44360.1 44360.1 52869.7 52869.6 52869.6 52869.5 52869.5 52869.4 52869.4 52869.4
54506.3 54506 54505.7 54505.5 54505.3 54505.1 54505 54504.9 54935 54934.4 54933.9 54933.4
54933.1 54932.8 54932.5 54932.3 55123.7 55122.5 55121.6 55120.9 55120.3 55119.8 55119.4 55119
54834.1 54850.3 54862.3 54871.7 54879.5 54886.3 54892.4 54897.7 54413.3 54434.4 54448.9 54460.2
54469.8 54478 54485.4 54492.1
8018.61 8018.61 8018.6 8018.59 8018.59 8018.58 8018.57 8018.57 44304.4 44304.4 44304.4 44304.4
44304.4 44304.4 44304.4 44304.4 52803.3 52803.2 52803.2 52803.1 52803.1 52803.1 52803 52803
54437.9 54437.6 54437.3 54437.1 54436.9 54436.7 54436.6 54436.5 54866.1 54865.4 54864.9 54864.5
54864.1 54863.8 54863.6 54863.4 55054.5 55053.4 55052.4 55051.7 55051.1 55050.6 55050.2 55049.8
54792 54806.9 54818 54826.6 54833.8 54840.1 54845.7 54850.6 54404 54423.4 54436.8 54447.3
54456 54463.7 54470.4 54476.7
</DataArray>
<DataArray type="Float32" Name="pc" NumberOfComponents="1" format="ascii">
1506.34 1506.34 1506.34 1506.34 1506.34 1506.34 1506.34 1506.34 2212.36 2212.36 2212.36 2212.36
2212.36 2212.35 2212.35 2212.35 2684.11 2684.1 2684.09 2684.09 2684.08 2684.08 2684.07 2684.07
3014.27 3014.15 3014.04 3013.95 3013.87 3013.81 3013.76 3013.71 3298.2 3297.43 3296.81 3296.31
3295.9 3295.55 3295.24 3294.99 3661.17 3657.57 3654.85 3652.73 3650.99 3649.48 3648.16 3647.02
1506.34 1506.34 1506.34 1506.34 1506.34 1506.34 1506.34 1506.34 2212.35 2212.35 2212.35 2212.35
2212.35 2212.35 2212.35 2212.35 2684.1 2684.1 2684.09 2684.08 2684.08 2684.07 2684.07 2684.07
3014.27 3014.14 3014.03 3013.94 3013.87 3013.81 3013.75 3013.71 3298.19 3297.43 3296.8 3296.31
3295.9 3295.55 3295.24 3294.99 3661.16 3657.57 3654.84 3652.73 3650.98 3649.47 3648.15 3647.02
4302.49 4302.49 4302.49 4302.49 4302.49 4302.49 4302.49 4302.49 4302.49 4302.49 4302.49 4302.49
4302.49 4302.49 4302.49 4302.49
</DataArray>
......@@ -88,16 +88,16 @@
0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987
0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987
0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987
0.612125 0.630316 0.643736 0.654269 0.663064 0.67069 0.677478 0.683465 0.136476 0.160449 0.176937 0.189862
0.200668 0.210037 0.218387 0.226066
0.612137 0.630328 0.64375 0.654283 0.663079 0.670705 0.677494 0.683482 0.136479 0.160454 0.176943 0.189869
0.200676 0.210045 0.218396 0.226076
</DataArray>
<DataArray type="Float32" Name="X^H2O_liq" NumberOfComponents="1" format="ascii">
0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978
0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978
0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978
0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978
0.612117 0.630308 0.643728 0.65426 0.663056 0.670681 0.677469 0.683457 0.136467 0.160441 0.176928 0.189854
0.20066 0.210028 0.218378 0.226057
0.612128 0.63032 0.643741 0.654275 0.663071 0.670697 0.677486 0.683474 0.136471 0.160446 0.176934 0.189861
0.200668 0.210037 0.218388 0.226067
</DataArray>
<DataArray type="Float32" Name="x^Air_liq" NumberOfComponents="1" format="ascii">
1.34035e-05 1.34035e-05 1.34035e-05 1.34035e-05 1.34035e-05 1.34035e-05 1.34035e-05 1.34035e-05 1.34034e-05 1.34034e-05 1.34034e-05 1.34034e-05
......@@ -128,31 +128,31 @@
0.0317391 0.0317391 0.0317391 0.0317391 0.0317393 0.0317393 0.0317393 0.0317393 0.0317393 0.0317393 0.0317393 0.0317393
0.0317394 0.0317394 0.0317394 0.0317394 0.0317394 0.0317394 0.0317394 0.0317394 0.0317395 0.0317395 0.0317395 0.0317395
0.0317395 0.0317395 0.0317395 0.0317395 0.0317396 0.0317396 0.0317396 0.0317396 0.0317396 0.0317396 0.0317396 0.0317396
0.0194289 0.0200062 0.0204322 0.0207665 0.0210457 0.0212877 0.0215032 0.0216932 0.00433175 0.00509267 0.00561598 0.00602624
0.00636922 0.00666657 0.0069316 0.00717533
0.0194292 0.0200066 0.0204326 0.020767 0.0210461 0.0212882 0.0215037 0.0216937 0.00433186 0.00509282 0.00561617 0.00602646
0.00636947 0.00666685 0.0069319 0.00717566
</DataArray>
<DataArray type="Float32" Name="X^H2O_gas" NumberOfComponents="1" format="ascii">
0.0199835 0.0199835 0.0199835 0.0199835 0.0199835 0.0199835 0.0199835 0.0199835 0.0199837 0.0199837 0.0199837 0.0199837
0.0199837 0.0199837 0.0199837 0.0199837 0.0199838 0.0199838 0.0199838 0.0199838 0.0199838 0.0199838 0.0199838 0.0199838
0.0199839 0.0199839 0.0199839 0.0199839 0.0199839 0.0199839 0.0199839 0.0199839 0.0199839 0.0199839 0.0199839 0.0199839
0.0199839 0.0199839 0.0199839 0.0199839 0.019984 0.019984 0.019984 0.019984 0.019984 0.019984 0.019984 0.019984
0.0121755 0.0125401 0.0128092 0.0130204 0.0131969 0.0133499 0.0134861 0.0136063 0.00269907 0.00317411 0.00350097 0.00375731
0.00397167 0.00415756 0.00432327 0.0044757
0.0121758 0.0125404 0.0128095 0.0130207 0.0131972 0.0133502 0.0134864 0.0136066 0.00269914 0.00317421 0.00350109 0.00375745
0.00397183 0.00415773 0.00432346 0.00447591
</DataArray>
<DataArray type="Float32" Name="x^Air_gas" NumberOfComponents="1" format="ascii">
0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261
0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261
0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.968261 0.96826 0.96826 0.96826 0.96826
0.96826 0.96826 0.96826 0.96826 0.96826 0.96826 0.96826 0.96826 0.96826 0.96826 0.96826 0.96826
0.980571 0.979994 0.979568 0.979234 0.978954 0.978712 0.978497 0.978307 0.995668 0.994907 0.994384 0.993974
0.993631 0.993333 0.993068 0.992825
0.980571 0.979993 0.979567 0.979233 0.978954 0.978712 0.978496 0.978306 0.995668 0.994907 0.994384 0.993974
0.993631 0.993333 0.993068 0.992824
</DataArray>
<DataArray type="Float32" Name="X^Air_gas" NumberOfComponents="1" format="ascii">
0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016
0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016
0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016
0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016 0.980016
0.987824 0.98746 0.987191 0.98698 0.986803 0.98665 0.986514 0.986394 0.997301 0.996826 0.996499 0.996243
0.987824 0.98746 0.987191 0.986979 0.986803 0.98665 0.986514 0.986393 0.997301 0.996826 0.996499 0.996243
0.996028 0.995842 0.995677 0.995524
</DataArray>
<DataArray type="Float32" Name="rhoMolar_gas" NumberOfComponents="1" format="ascii">
......
......@@ -18,7 +18,7 @@
</DataArray>
<DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
1.16717 1.16683 1.16659 1.1664 1.16625 1.16612 1.166 1.16588 1.16725 1.1669 1.16667 1.16648
1.16633 1.16619 1.16607 1.16596 1.16736 1.16702 1.16678 1.16659 1.16644 1.1663 1.16618 1.16607
1.16632 1.16619 1.16607 1.16596 1.16736 1.16702 1.16678 1.16659 1.16644 1.1663 1.16618 1.16607
1.16752 1.16718 1.16695 1.16676 1.1666 1.16647 1.16635 1.16624 1.16773 1.16742 1.16719 1.167
1.16685 1.16671 1.16659 1.16648 1.16796 1.16771 1.16751 1.16734 1.16719 1.16706 1.16694 1.16683
1.16814 1.16801 1.16787 1.16774 1.16761 1.1675 1.16739 1.16729 1.16822 1.16818 1.16814 1.16808
......@@ -43,10 +43,10 @@
40.3385 40.3385 40.3385 40.3385 40.3385 40.3385 40.3385 40.3385
</DataArray>
<DataArray type="Float32" Name="X^Air_gas" NumberOfComponents="1" format="ascii">
0.998508 0.998017 0.997683 0.997421 0.997202 0.997012 0.996842 0.996685 0.998617 0.998124 0.99779 0.997527
0.997308 0.997117 0.996948 0.99679 0.998777 0.998284 0.997949 0.997687 0.997467 0.997276 0.997105 0.996947
0.999003 0.998519 0.998185 0.997922 0.997702 0.997511 0.99734 0.997182 0.999297 0.998848 0.998524 0.998265
0.998046 0.997856 0.997685 0.997526 0.999618 0.999267 0.998978 0.998736 0.998527 0.998342 0.998176 0.998019
0.998508 0.998016 0.997683 0.997421 0.997202 0.997011 0.996842 0.996685 0.998617 0.998124 0.99779 0.997527
0.997308 0.997117 0.996948 0.99679 0.998777 0.998284 0.997949 0.997686 0.997466 0.997276 0.997105 0.996947
0.999003 0.998519 0.998185 0.997922 0.997702 0.997511 0.99734 0.997181 0.999297 0.998848 0.998524 0.998264
0.998046 0.997855 0.997685 0.997526 0.999618 0.999267 0.998978 0.998735 0.998527 0.998342 0.998176 0.998018
0.999871 0.999686 0.99949 0.999303 0.999129 0.998969 0.99882 0.998674 0.99998 0.999934 0.999866 0.999784
0.999694 0.9996 0.999503 0.999398 0.999999 0.999995 0.999986 0.999973 0.999956 0.999934 0.999908 0.999874
1 1 0.999999 0.999999 0.999997 0.999996 0.999994 0.99999 1 1 1 1
......@@ -56,12 +56,12 @@
1 1 1 1 1 1 1 1
</DataArray>
<DataArray type="Float32" Name="x^Air_gas" NumberOfComponents="1" format="ascii">
0.997604 0.996815 0.99628 0.99586 0.995509 0.995205 0.994933 0.994682 0.997779 0.996988 0.996452 0.996031
0.995679 0.995374 0.995102 0.99485 0.998035 0.997245 0.996708 0.996286 0.995934 0.995628 0.995355 0.995102
0.997604 0.996815 0.99628 0.99586 0.995509 0.995205 0.994933 0.994682 0.997778 0.996988 0.996452 0.996031
0.995679 0.995374 0.995102 0.99485 0.998035 0.997245 0.996708 0.996286 0.995933 0.995628 0.995355 0.995102
0.998398 0.997621 0.997086 0.996664 0.996311 0.996005 0.995731 0.995477 0.99887 0.99815 0.997629 0.997213
0.996863 0.996557 0.996284 0.996029 0.999387 0.998822 0.998358 0.997969 0.997634 0.997338 0.997071 0.996819
0.999793 0.999495 0.999181 0.99888 0.998601 0.998344 0.998105 0.997871 0.999968 0.999894 0.999785 0.999653
0.999508 0.999357 0.999201 0.999033 0.999998 0.999991 0.999978 0.999957 0.999929 0.999894 0.999852 0.999798
0.996862 0.996557 0.996284 0.996029 0.999387 0.998822 0.998358 0.997969 0.997633 0.997338 0.997071 0.996818
0.999793 0.999495 0.99918 0.99888 0.998601 0.998344 0.998105 0.997871 0.999968 0.999894 0.999785 0.999653
0.999508 0.999356 0.999201 0.999033 0.999998 0.999991 0.999978 0.999957 0.999929 0.999894 0.999852 0.999798
1 1 0.999999 0.999998 0.999996 0.999993 0.99999 0.999984 1 1 1 1
1 1 1 0.999999 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1
......@@ -69,30 +69,30 @@
1 1 1 1 1 1 1 1
</DataArray>
<DataArray type="Float32" Name="X^H2O_gas" NumberOfComponents="1" format="ascii">
0.00149192 0.00198346 0.00231731 0.00257926 0.00279838 0.0029884 0.00315781 0.00331484 0.00138305 0.00187583 0.00221025 0.00247267
0.00269221 0.0028826 0.00305233 0.00320978 0.00122309 0.00171565 0.0020505 0.00231345 0.00253349 0.00272434 0.00289451 0.00305256
0.000997052 0.00148114 0.00181472 0.00207757 0.00229785 0.00248906 0.00265963 0.00281846 0.000703349 0.00115159 0.00147611 0.00173541
0.00195401 0.00214438 0.00231459 0.00247397 0.000381658 0.000732888 0.00102197 0.00126439 0.00147336 0.00165762 0.00182377 0.00198141
0.000129069 0.00031439 0.000509901 0.000697204 0.000870926 0.00103097 0.00117972 0.00132564 2.02138e-05 6.62085e-05 0.000133677 0.00021562
0.000305918 0.000400399 0.000497303 0.000601752 1.14304e-06 5.32299e-06 1.3785e-05 2.68527e-05 4.43156e-05 6.59113e-05 9.19227e-05 0.000125897
2.2131e-08 1.68022e-07 5.72937e-07 1.33036e-06 2.48681e-06 4.09501e-06 6.33481e-06 1.0012e-05 1.77603e-10 2.77914e-09 1.28338e-08 3.47577e-08
7.05942e-08 1.22727e-07 2.02435e-07 3.59451e-07 9.21173e-13 3.39729e-11 2.11376e-10 6.54976e-10 1.4099e-09 2.50852e-09 4.22372e-09 7.79319e-09
4.55821e-15 3.80956e-13 3.07831e-12 1.07179e-11 2.42397e-11 4.37258e-11 7.2934e-11 1.29275e-10 0 4.42008e-15 4.37864e-14 1.67272e-13
3.93802e-13 7.16329e-13 1.16607e-12 1.915e-12 0 0 6.90638e-16 2.62442e-15 6.49199e-15 1.17408e-14 1.85091e-14 2.8178e-14
0 0 1.38128e-16 0 1.38128e-16 1.38128e-16 2.76255e-16 5.5251e-16
0.00149197 0.00198354 0.00231742 0.00257938 0.00279852 0.00298856 0.00315799 0.00331503 0.0013831 0.00187591 0.00221035 0.00247279
0.00269235 0.00288275 0.0030525 0.00320996 0.00122313 0.00171573 0.0020506 0.00231357 0.00253362 0.0027245 0.00289468 0.00305275
0.000997094 0.00148121 0.00181482 0.00207769 0.00229798 0.00248921 0.0026598 0.00281864 0.000703385 0.00115165 0.0014762 0.00173552
0.00195414 0.00214453 0.00231475 0.00247414 0.000381685 0.00073294 0.00102204 0.00126449 0.00147348 0.00165775 0.00182392 0.00198157
0.000129082 0.000314423 0.000509955 0.000697278 0.000871019 0.00103108 0.00117984 0.00132578 2.02169e-05 6.62206e-05 0.0001337 0.000215658
0.00030597 0.000400465 0.000497383 0.000601846 1.14322e-06 5.32487e-06 1.37898e-05 2.68615e-05 4.43295e-05 6.59311e-05 9.1949e-05 0.000125931
2.21549e-08 1.68162e-07 5.73342e-07 1.33121e-06 2.4883e-06 4.09734e-06 6.33817e-06 1.00166e-05 1.78045e-10 2.78332e-09 1.28492e-08 3.47951e-08
7.06655e-08 1.22845e-07 2.02617e-07 3.59728e-07 9.24764e-13 3.40443e-11 2.1173e-10 6.5596e-10 1.41191e-09 2.51195e-09 4.22917e-09 7.80194e-09
4.55821e-15 3.81923e-13 3.08466e-12 1.07376e-11 2.42823e-11 4.38001e-11 7.30526e-11 1.29466e-10 1.38128e-16 4.42008e-15 4.39246e-14 1.67687e-13
3.94492e-13 7.17711e-13 1.16814e-12 1.91832e-12 0 1.38128e-16 6.90638e-16 2.62442e-15 6.49199e-15 1.1879e-14 1.85091e-14 2.83161e-14
0 -6.90638e-17 -6.90638e-17 0 1.38128e-16 1.38128e-16 2.76255e-16 5.5251e-16
</DataArray>
<DataArray type="Float32" Name="x^H2O_gas" NumberOfComponents="1" format="ascii">
0.00239614 0.00318465 0.00371992 0.00413975 0.00449085 0.00479525 0.00506657 0.00531801 0.00222143 0.00301203 0.00354829 0.00396894
0.00432074 0.00462577 0.00489764 0.00514978 0.00196469 0.0027551 0.00329215 0.00371372 0.0040664 0.00437223 0.00464486 0.00489802
0.00160182 0.00237884 0.00291402 0.00333556 0.00368872 0.0039952 0.00426855 0.00452302 0.00113017 0.00184992 0.00237077 0.0027868
0.00313742 0.00344268 0.00371556 0.00397102 0.000613385 0.00117762 0.00164183 0.00203099 0.00236636 0.002662 0.00292853 0.00318136
0.000207466 0.000505295 0.00081943 0.0011203 0.0013993 0.00165628 0.00189508 0.0021293 3.24939e-05 0.000106428 0.000214872 0.000346571
0.000491682 0.000643498 0.000799189 0.000966983 1.83748e-06 8.55685e-06 2.21597e-05 4.31659e-05 7.1237e-05 0.00010595 0.000147761 0.000202369
3.55763e-08 2.70101e-07 9.21014e-07 2.13859e-06 3.99762e-06 6.58285e-06 1.01834e-05 1.60944e-05 2.85502e-10 4.46756e-09 2.06307e-08 5.58742e-08
1.13483e-07 1.97287e-07 3.25421e-07 5.7783e-07 1.48082e-12 5.46125e-11 3.39795e-10 1.0529e-09 2.26646e-09 4.03252e-09 6.78977e-09 1.25278e-08
7.32747e-15 6.12399e-13 4.94849e-12 1.72293e-11 3.89662e-11 7.02907e-11 1.17244e-10 2.07815e-10 0 7.10543e-15 7.03881e-14 2.68896e-13
6.33049e-13 1.15152e-12 1.8745e-12 3.07843e-12 0 0 1.11022e-15 4.21885e-15 1.04361e-14 1.88738e-14 2.9754e-14 4.52971e-14
0 0 2.22045e-16 0 2.22045e-16 2.22045e-16 4.44089e-16 8.88178e-16
0.00239622 0.00318478 0.00372009 0.00413995 0.00449108 0.0047955 0.00506685 0.00531831 0.00222151 0.00301215 0.00354845 0.00396914
0.00432097 0.00462602 0.00489792 0.00515008 0.00196477 0.00275522 0.00329231 0.00371391 0.00406662 0.00437248 0.00464513 0.00489831
0.00160189 0.00237896 0.00291417 0.00333574 0.00368893 0.00399544 0.00426881 0.00452331 0.00113023 0.00185002 0.00237091 0.00278697
0.00313762 0.00344291 0.00371582 0.00397129 0.000613429 0.0011777 0.00164195 0.00203115 0.00236654 0.00266221 0.00292877 0.00318161
0.000207488 0.000505349 0.000819516 0.00112042 0.00139945 0.00165646 0.00189528 0.00212952 3.24989e-05 0.000106448 0.00021491 0.000346631
0.000491765 0.000643604 0.000799318 0.000967134 1.83776e-06 8.55988e-06 2.21673e-05 4.318e-05 7.12592e-05 0.000105982 0.000147803 0.000202423
3.56147e-08 2.70326e-07 9.21666e-07 2.13996e-06 4.00001e-06 6.58659e-06 1.01888e-05 1.6102e-05 2.86213e-10 4.47427e-09 2.06556e-08 5.59342e-08
1.13597e-07 1.97477e-07 3.25714e-07 5.78275e-07 1.48659e-12 5.47273e-11 3.40363e-10 1.05448e-09 2.26969e-09 4.03804e-09 6.79853e-09 1.25419e-08
7.32747e-15 6.13953e-13 4.9587e-12 1.72611e-11 3.90346e-11 7.04101e-11 1.17435e-10 2.0812e-10 2.22045e-16 7.10543e-15 7.06102e-14 2.69562e-13
6.34159e-13 1.15374e-12 1.87783e-12 3.08376e-12 0 2.22045e-16 1.11022e-15 4.21885e-15 1.04361e-14 1.90958e-14 2.9754e-14 4.55191e-14
0 -1.11022e-16 -1.11022e-16 0 2.22045e-16 2.22045e-16 4.44089e-16 8.88178e-16
</DataArray>
<DataArray type="Float32" Name="D^H2O_gas" NumberOfComponents="1" format="ascii">
2.49368e-05 2.49368e-05 2.49368e-05 2.49368e-05 2.49368e-05 2.49368e-05 2.49368e-05 2.49368e-05 2.49368e-05 2.49368e-05 2.49368e-05 2.49368e-05
......@@ -108,38 +108,38 @@
2.49374e-05 2.49374e-05 2.49374e-05 2.49374e-05 2.49374e-05 2.49374e-05 2.49374e-05 2.49374e-05
</DataArray>
<DataArray type="Float32" Name="velocity_gas (m/s)" NumberOfComponents="3" format="ascii">
0.00578835 8.32189e-06 0 0.0062025 1.14379e-05 0 0.006114 1.15037e-05 0 0.00602194 1.14602e-05 0
0.00594526 1.13929e-05 0 0.00588541 1.13473e-05 0 0.00583119 1.13753e-05 0 0.00569999 1.19281e-05 0
0.0191001 1.54974e-06 0 0.0194232 1.28385e-05 0 0.0191467 1.32181e-05 0 0.018859 1.29695e-05 0
0.0186191 1.25674e-05 0 0.0184318 1.22681e-05 0 0.0182622 1.23512e-05 0 0.017861 1.54334e-05 0
0.0390171 -8.21301e-06 0 0.0392015 1.74898e-05 0 0.0386451 1.89101e-05 0 0.0380657 1.79859e-05 0
0.0375825 1.64769e-05 0 0.0372048 1.53378e-05 0 0.0368627 1.56008e-05 0 0.0360786 2.66777e-05 0
0.0687786 -2.18959e-05 0 0.0687506 3.01327e-05 0 0.0677788 3.43886e-05 0 0.0667664 3.16422e-05 0
0.0659209 2.71308e-05 0 0.0652594 2.37107e-05 0 0.0646605 2.44507e-05 0 0.0633477 5.57893e-05 0
0.113166 -3.99757e-05 0 0.11281 6.16379e-05 0 0.111225 7.29947e-05 0 0.109572 6.576e-05 0
0.10819 5.37875e-05 0 0.107106 4.46876e-05 0 0.106127 4.65586e-05 0 0.104112 0.000123353 0
0.17918 -6.06228e-05 0 0.178315 0.000136415 0 0.175827 0.00016476 0 0.173233 0.00014706 0
0.171059 0.00011745 0 0.169351 9.4875e-05 0 0.16781 9.91745e-05 0 0.164914 0.000268105 0
0.276939 -7.38327e-05 0 0.275282 0.000308465 0 0.271476 0.000376046 0 0.267509 0.000334909 0
0.264174 0.000265022 0 0.261549 0.000211495 0 0.25919 0.00022029 0 0.255261 0.000556314 0
0.420775 -4.5111e-05 0 0.417903 0.00069675 0 0.412189 0.000851227 0 0.406241 0.000758985 0
0.401224 0.000599538 0 0.397259 0.000476579 0 0.393727 0.000491114 0 0.388657 0.00109051 0
0.630286 0.000123753 0 0.625597 0.00156115 0 0.617177 0.00190021 0 0.608434 0.00169872 0
0.601019 0.00134444 0 0.595133 0.00106819 0 0.589964 0.00107906 0 0.583644 0.00202226 0
0.930578 0.000684589 0 0.923242 0.00345355 0 0.911115 0.00416754 0 0.898574 0.00373923 0
0.887873 0.00297261 0 0.879334 0.00236336 0 0.871996 0.00231399 0 0.864257 0.00359109 0
1.34967 0.0022416 0 1.33861 0.00750669 0 1.32163 0.00892921 0 1.30422 0.00804006 0
1.2893 0.0064293 0 1.27736 0.00511023 0 1.26739 0.00479388 0 1.2581 0.00619226 0
1.90797 0.0061329 0 1.89197 0.0159432 0 1.86909 0.0185394 0 1.84605 0.016703 0
1.82633 0.0134208 0 1.81058 0.0106253 0 1.79792 0.0094639 0 1.78724 0.0104239 0
2.5874 0.0149906 0 2.56566 0.0327265 0 2.53666 0.0367774 0 2.50852 0.0328827 0
2.48483 0.026397 0 2.46622 0.0206793 0 2.45189 0.0174329 0 2.44068 0.0169605 0
3.2511 0.03314 0 3.22474 0.063064 0 3.19278 0.0673258 0 3.164 0.0587654 0
3.14103 0.0465244 0 3.1238 0.0357295 0 3.11124 0.0286375 0 3.10199 0.0256038 0
3.44725 0.0640355 0 3.42495 0.103159 0 3.40477 0.101114 0 3.39096 0.0837216 0
3.38247 0.063938 0 3.37753 0.0475677 0 3.37476 0.0366562 0 3.37315 0.0314175 0
1.96233 0.0414408 0 2.00107 0.0620094 0 2.04604 0.0580987 0 2.08525 0.0466093 0
2.11572 0.0347753 0 2.13824 0.025377 0 2.15494 0.0192078 0 2.16821 0.0162565 0
0.00578806 8.32416e-06 0 0.00620183 1.14388e-05 0 0.00611318 1.15043e-05 0 0.00602104 1.14606e-05 0
0.00594428 1.13932e-05 0 0.00588438 1.13476e-05 0 0.00583012 1.13754e-05 0 0.0056989 1.19281e-05 0
0.0190992 1.56204e-06 0 0.0194212 1.28423e-05 0 0.0191442 1.32205e-05 0 0.0188562 1.29712e-05 0
0.0186161 1.25688e-05 0 0.0184287 1.22692e-05 0 0.0182588 1.2352e-05 0 0.0178577 1.54335e-05 0
0.0390153 -8.16884e-06 0 0.0391975 1.75036e-05 0 0.0386402 1.89187e-05 0 0.0380603 1.79921e-05 0
0.0375766 1.64818e-05 0 0.0371986 1.53417e-05 0 0.0368562 1.56038e-05 0 0.036072 2.66783e-05 0
0.0687759 -2.17711e-05 0 0.0687441 3.01737e-05 0 0.0677708 3.4414e-05 0 0.0667573 3.16607e-05 0
0.065911 2.71454e-05 0 0.065249 2.37226e-05 0 0.0646495 2.44597e-05 0 0.0633364 5.57916e-05 0
0.113162 -3.96717e-05 0 0.112801 6.17464e-05 0 0.111213 7.30627e-05 0 0.109558 6.58097e-05 0
0.108175 5.38268e-05 0 0.10709 4.47194e-05 0 0.10611 4.65828e-05 0 0.104095 0.000123361 0
0.179175 -5.99637e-05 0 0.178302 0.000136679 0 0.175811 0.00016493 0 0.173214 0.000147184 0
0.171038 0.000117548 0 0.169329 9.49553e-05 0 0.167786 9.92359e-05 0 0.16489 0.000268128 0
0.276933 -7.25541e-05 0 0.275267 0.000309057 0 0.271456 0.00037644 0 0.267485 0.000335205 0
0.264148 0.000265258 0 0.26152 0.000211687 0 0.259159 0.000220439 0 0.255228 0.000556382 0
0.420769 -4.2897e-05 0 0.417889 0.000697956 0 0.412168 0.000852068 0 0.406216 0.000759631 0
0.401194 0.000600062 0 0.397225 0.000477011 0 0.39369 0.000491453 0 0.388619 0.00109068 0
0.630282 0.000127185 0 0.625584 0.00156336 0 0.617157 0.00190182 0 0.608409 0.0017 0
0.60099 0.00134549 0 0.595099 0.00106906 0 0.589927 0.00107975 0 0.583605 0.00202265 0
0.930575 0.000689382 0 0.923232 0.00345724 0 0.911099 0.00417035 0 0.898552 0.00374147 0
0.887847 0.00297447 0 0.879304 0.0023649 0 0.871962 0.00231522 0 0.864222 0.00359182 0
1.34967 0.00224763 0 1.3386 0.00751233 0 1.32162 0.00893372 0 1.30421 0.0080437 0
1.28928 0.0064323 0 1.27733 0.00511269 0 1.26736 0.00479579 0 1.25807 0.00619347 0
1.90797 0.00613968 0 1.89197 0.015951 0 1.86908 0.0185461 0 1.84604 0.0167085 0
1.82631 0.0134252 0 1.81057 0.0106289 0 1.7979 0.00946659 0 1.78722 0.0104257 0
2.5874 0.0149972 0 2.56566 0.0327359 0 2.53666 0.0367861 0 2.50852 0.0328901 0
2.48482 0.026403 0 2.46621 0.0206839 0 2.45188 0.0174363 0 2.44067 0.0169629 0
3.2511 0.0331455 0 3.22474 0.0630734 0 3.19278 0.0673354 0 3.164 0.0587738 0
3.14103 0.0465312 0 3.1238 0.0357348 0 3.11124 0.0286414 0 3.10199 0.0256066 0
3.44725 0.0640389 0 3.42496 0.103165 0 3.40477 0.101122 0 3.39097 0.0837285 0
3.38248 0.0639438 0 3.37754 0.0475722 0 3.37477 0.0366596 0 3.37316 0.0314201 0
1.96233 0.0414418 0 2.00107 0.0620118 0 2.04604 0.0581015 0 2.08525 0.046612 0
2.11573 0.0347776 0 2.13825 0.0253788 0 2.15495 0.0192092 0 2.16823 0.0162576 0
</DataArray>
<DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
0 0 0 0 0 0 0 0 0 0 0 0
......