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

ENH: surface film - updates to help smooth restart behaviour

Browse Source
This commit is contained in:
andy
2014-04-30 16:55:03 +01:00
committed by Andrew Heather
parent e4bb5fa666
commit 17162a96cb
2 changed files with 36 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

23
src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -581,11 +581,11 @@ kinematicSingleLayer::kinematicSingleLayer
"phi",
time().timeName(),
regionMesh(),
IOobject::READ_IF_PRESENT,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
regionMesh(),
dimLength*dimMass/dimTime
dimensionedScalar("0", dimLength*dimMass/dimTime, 0.0)
),
primaryMassTrans_
@ -795,7 +795,22 @@ kinematicSingleLayer::kinematicSingleLayer
correctThermoFields();
deltaRho_ == delta_*rho_;
phi_ = fvc::interpolate(deltaRho_*U_) & regionMesh().Sf();
surfaceScalarField phi0
(
IOobject
(
"phi",
time().timeName(),
regionMesh(),
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE,
false
),
fvc::interpolate(deltaRho_*U_) & regionMesh().Sf()
);
phi_ == phi0;
}
}

18
src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C
View File

@ -554,8 +554,24 @@ thermoSingleLayer::thermoSingleLayer
// Update derived fields
hs_ == hs(T_);
deltaRho_ == delta_*rho_;
phi_ = fvc::interpolate(deltaRho_*U_) & regionMesh().Sf();
surfaceScalarField phi0
(
IOobject
(
"phi",
time().timeName(),
regionMesh(),
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE,
false
),
fvc::interpolate(deltaRho_*U_) & regionMesh().Sf()
);
phi_ == phi0;
// evaluate viscosity from user-model
viscosity_->correct(pPrimary_, T_);