...
 
Commits (2)
...@@ -6,6 +6,7 @@ ...@@ -6,6 +6,7 @@
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
Copyright (C) 2011-2016 OpenFOAM Foundation Copyright (C) 2011-2016 OpenFOAM Foundation
Copyright (C) 2020 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
...@@ -177,7 +178,11 @@ Foam::LimitedScheme<Type, Limiter, LimitFunc>::limiter ...@@ -177,7 +178,11 @@ Foam::LimitedScheme<Type, Limiter, LimitFunc>::limiter
calcLimiter(phi, limiterField); calcLimiter(phi, limiterField);
return limiterField; return tmp<surfaceScalarField>::New
(
limiterFieldName,
limiterField
);
} }
else else
{ {
......
...@@ -6,6 +6,7 @@ ...@@ -6,6 +6,7 @@
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
Copyright (C) 2013-2019 OpenFOAM Foundation Copyright (C) 2013-2019 OpenFOAM Foundation
Copyright (C) 2020 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
...@@ -704,7 +705,7 @@ void Foam::multiphaseSystem::solve() ...@@ -704,7 +705,7 @@ void Foam::multiphaseSystem::solve()
phase.alphaRhoPhiRef() = phase.alphaRhoPhiRef() =
fvc::interpolate(phase.rho())*phase.alphaPhi(); fvc::interpolate(phase.rho())*phase.alphaPhi();
phase.clip(0, 1); phase.clip(SMALL, 1 - SMALL);
} }
calcAlphas(); calcAlphas();
......
...@@ -6,6 +6,7 @@ ...@@ -6,6 +6,7 @@
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
Copyright (C) 2013-2018 OpenFOAM Foundation Copyright (C) 2013-2018 OpenFOAM Foundation
Copyright (C) 2020 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
...@@ -339,7 +340,7 @@ void Foam::twoPhaseSystem::solve() ...@@ -339,7 +340,7 @@ void Foam::twoPhaseSystem::solve()
<< endl; << endl;
// Ensure the phase-fractions are bounded // Ensure the phase-fractions are bounded
alpha1.clip(0, 1); alpha1.clip(SMALL, 1 - SMALL);
// Update the phase-fraction of the other phase // Update the phase-fraction of the other phase
alpha2 = scalar(1) - alpha1; alpha2 = scalar(1) - alpha1;
......