zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Fixed a bug introduced by maintaining the coupled patch fields in

Browse Source 
sliced fields.
This commit is contained in:
henry
2008-06-24 17:01:01 +01:00
parent b2b5ace228
commit 7639571b4b
3 changed files with 30 additions and 16 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
View File

@ -55,6 +55,7 @@ void Foam::MULES::explicitSolve
Info<< "MULES: Solving for " << psi.name() << endl;
const fvMesh& mesh = psi.mesh();
psi.correctBoundaryConditions();
surfaceScalarField phiBD = upwind<scalar>(psi.mesh(), phi).flux(psi);
@ -76,7 +77,8 @@ void Foam::MULES::explicitSolve
),
mesh,
dimless,
allLambda
allLambda,
false // Use slices for the couples
);
limiter
@ -183,7 +185,8 @@ void Foam::MULES::implicitSolve
),
mesh,
dimless,
allCoLambda
allCoLambda,
false // Use slices for the couples
);
CoLambda == 1.0/max(CoCoeff*Cof, scalar(1));
@ -205,7 +208,8 @@ void Foam::MULES::implicitSolve
),
mesh,
dimless,
allLambda
allLambda,
false // Use slices for the couples
);
linear<scalar> CDs(mesh);
@ -347,7 +351,8 @@ void Foam::MULES::limiter
),
mesh,
dimless,
allLambda
allLambda,
false // Use slices for the couples
);
scalarField& lambdaIf = lambda;
@ -406,6 +411,9 @@ void Foam::MULES::limiter
{
label pfCelli = pFaceCells[pFacei];
psiMaxn[pfCelli] = max(psiMaxn[pfCelli], psiPf[pFacei]);
psiMinn[pfCelli] = min(psiMinn[pfCelli], psiPf[pFacei]);
psiMaxn[pfCelli] = max(psiMaxn[pfCelli], psiPNf[pFacei]);
psiMinn[pfCelli] = min(psiMinn[pfCelli], psiPNf[pFacei]);
}