Skip to content
Snippets Groups Projects
Commit 919136bc authored by Anna Mareike Kostelecky's avatar Anna Mareike Kostelecky
Browse files

Merge branch...

Merge branch '1325-porousmediumflow-2pncmin-nonisothermal-use-current-permeability-for-neumann-gas-flux' into 'master'

[porousmediumflow][2pncmin][nonisothermal] Use current permeability for Neumann gas flux

Closes #1325

See merge request !3742
parents 36055b15 b8888a5a
No related branches found
No related tags found
1 merge request!3742[porousmediumflow][2pncmin][nonisothermal] Use current permeability for Neumann gas flux
Pipeline #40117 passed
+5
...@@ -217,7 +217,7 @@ public: ...@@ -217,7 +217,7 @@ public:
// get porous medium values: // get porous medium values:
const Scalar massFracH2OInside = volVars.massFraction(gasPhaseIdx, H2OIdx); const Scalar massFracH2OInside = volVars.massFraction(gasPhaseIdx, H2OIdx);
static const Scalar referencePermeability = getParam<Scalar>("SpatialParams.referencePermeability", 2.23e-14);
// calculate fluxes // calculate fluxes
// liquid phase // liquid phase
...@@ -248,7 +248,7 @@ public: ...@@ -248,7 +248,7 @@ public:
values[conti1EqIdx] = (volVars.pressure(gasPhaseIdx) - 1e5) values[conti1EqIdx] = (volVars.pressure(gasPhaseIdx) - 1e5)
/(globalPos - fvGeometry.scv(scvf.insideScvIdx()).center()).two_norm() /(globalPos - fvGeometry.scv(scvf.insideScvIdx()).center()).two_norm()
*volVars.mobility(gasPhaseIdx) *volVars.mobility(gasPhaseIdx)
*referencePermeability *volVars.permeability()
*volVars.molarDensity(gasPhaseIdx) *volVars.molarDensity(gasPhaseIdx)
*volVars.moleFraction(gasPhaseIdx, AirIdx); *volVars.moleFraction(gasPhaseIdx, AirIdx);
} }
...@@ -257,7 +257,7 @@ public: ...@@ -257,7 +257,7 @@ public:
values[conti1EqIdx] = (volVars.pressure(gasPhaseIdx) - 1e5) values[conti1EqIdx] = (volVars.pressure(gasPhaseIdx) - 1e5)
/(globalPos - fvGeometry.scv(scvf.insideScvIdx()).center()).two_norm() /(globalPos - fvGeometry.scv(scvf.insideScvIdx()).center()).two_norm()
*volVars.mobility(gasPhaseIdx) *volVars.mobility(gasPhaseIdx)
*referencePermeability *volVars.permeability()
*volVars.molarDensity(gasPhaseIdx) * (1-moleFracRefH2O); *volVars.molarDensity(gasPhaseIdx) * (1-moleFracRefH2O);
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment