zhyang-dev/OpenFOAM-6
1
0
Fork 0
mirror of https://github.com/OpenFOAM/OpenFOAM-6.git synced 2025-12-08 06:57:46 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

interRegionExplicitPorositySource: Corrected initialization of firstIter_

Browse Source 
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1978
This commit is contained in:
Henry Weller
2016-01-20 10:20:54 +00:00
parent 4426006d69
commit eb3cf2c880
1 changed files with 10 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
src/fvOptions/sources/interRegion/interRegionExplicitPorositySource/interRegionExplicitPorositySource.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2012-2015 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2012-2016 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -74,7 +74,7 @@ void Foam::fv::interRegionExplicitPorositySource::initialise()
new cellZone
(
zoneName,
nbrMesh.faceNeighbour(), // neighbour internal cells
nbrMesh.faceNeighbour(), // Neighbour internal cells
zoneID,
cellZones
)
@ -117,7 +117,7 @@ Foam::fv::interRegionExplicitPorositySource::interRegionExplicitPorositySource
:
interRegionOption(name, modelType, dict, mesh),
porosityPtr_(NULL),
firstIter_(-1),
firstIter_(false),
UName_(coeffs_.lookupOrDefault<word>("UName", "U")),
muName_(coeffs_.lookupOrDefault<word>("muName", "thermo:mu"))
{
@ -157,7 +157,7 @@ void Foam::fv::interRegionExplicitPorositySource::addSup
dimensionedVector("zero", U.dimensions(), vector::zero)
);
// map local velocity onto neighbour region
// Map local velocity onto neighbour region
meshInterp().mapSrcToTgt
(
U.internalField(),
@ -169,7 +169,7 @@ void Foam::fv::interRegionExplicitPorositySource::addSup
porosityPtr_->addResistance(nbrEqn);
// convert source from neighbour to local region
// Convert source from neighbour to local region
fvMatrix<vector> porosityEqn(U, eqn.dimensions());
scalarField& Udiag = porosityEqn.diag();
vectorField& Usource = porosityEqn.source();
@ -211,7 +211,7 @@ void Foam::fv::interRegionExplicitPorositySource::addSup
dimensionedVector("zero", U.dimensions(), vector::zero)
);
// map local velocity onto neighbour region
// Map local velocity onto neighbour region
meshInterp().mapSrcToTgt
(
U.internalField(),
@ -252,7 +252,7 @@ void Foam::fv::interRegionExplicitPorositySource::addSup
const volScalarField& mu =
mesh_.lookupObject<volScalarField>(muName_);
// map local rho onto neighbour region
// Map local rho onto neighbour region
meshInterp().mapSrcToTgt
(
rho.internalField(),
@ -260,7 +260,7 @@ void Foam::fv::interRegionExplicitPorositySource::addSup
rhoNbr.internalField()
);
// map local mu onto neighbour region
// Map local mu onto neighbour region
meshInterp().mapSrcToTgt
(
mu.internalField(),
@ -270,7 +270,7 @@ void Foam::fv::interRegionExplicitPorositySource::addSup
porosityPtr_->addResistance(nbrEqn, rhoNbr, muNbr);
// convert source from neighbour to local region
// Convert source from neighbour to local region
fvMatrix<vector> porosityEqn(U, eqn.dimensions());
scalarField& Udiag = porosityEqn.diag();
vectorField& Usource = porosityEqn.source();
@ -292,7 +292,7 @@ bool Foam::fv::interRegionExplicitPorositySource::read(const dictionary& dict)
coeffs_.readIfPresent("UName", UName_);
coeffs_.readIfPresent("muName", muName_);
// reset the porosity model?
// Reset the porosity model?
return true;
}