zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

filmSurfaceVelocityFvPatchVectorField: Corrected to use the dynamic rather than kinematic viscosity

Browse Source
This commit is contained in:
Henry Weller
2023-05-05 21:55:03 +01:00
parent ab71ac6c1f
commit 6e9673e1e2
1 changed files with 31 additions and 16 deletions
Download Patch File Download Diff File Expand all files Collapse all files

47
applications/solvers/modules/isothermalFilm/derivedFvPatchFields/filmSurfaceVelocity/filmSurfaceVelocityFvPatchVectorField.C
View File

@ -24,7 +24,7 @@ License
\*---------------------------------------------------------------------------*/
#include "filmSurfaceVelocityFvPatchVectorField.H"
#include "momentumTransportModel.H"
#include "compressibleMomentumTransportModel.H"
#include "mappedPatchBase.H"
#include "addToRunTimeSelectionTable.H"
@ -122,6 +122,7 @@ void Foam::filmSurfaceVelocityFvPatchVectorField::updateCoeffs()
const fvPatch& patchNbr =
refCast<const fvMesh>(mpp.nbrMesh()).boundary()[patchiNbr];
// Neighbour patch face-cell velocity
const vectorField UpNbr
(
patchNbr.lookupPatchField<volVectorField, scalar>("U")
@ -136,43 +137,57 @@ void Foam::filmSurfaceVelocityFvPatchVectorField::updateCoeffs()
refValue() -= n*(n & refValue());
// Lookup the momentum transport model
const momentumTransportModel& transportModel =
db().lookupType<momentumTransportModel>();
const compressibleMomentumTransportModel& transportModel =
db().lookupType<compressibleMomentumTransportModel>();
// Get the patch laminar viscosity divided by delta
const tmp<scalarField> nuEffByDelta
// Patch density
const scalarField& rhop
(
transportModel.nuEff(patch().index())*patch().deltaCoeffs()
transportModel.rho().boundaryField()[patch().index()]
);
// Patch laminar dynamic viscosity divided by delta
const tmp<scalarField> muEffByDelta
(
rhop*transportModel.nuEff(patch().index())
*patch().deltaCoeffs()
);
// Lookup the neighbour momentum transport model
const compressibleMomentumTransportModel& transportModelNbr =
mpp.nbrMesh().lookupType<compressibleMomentumTransportModel>();
// Neighbour patch density
const scalarField& rhopNbr
(
transportModelNbr.rho().boundaryField()[patchiNbr]
);
if (Cs_ > 0)
{
// Calculate the drag coefficient from the drag constant
// and the magnitude of the velocity difference
const scalarField Ds(Cs_*mag(refValue() - *this));
const scalarField Ds(Cs_*rhopNbr*mag(refValue() - *this));
// Calculate the value-fraction from the balance between the
// external fluid drag and internal film stress
valueFraction() = Ds/(Ds + nuEffByDelta);
valueFraction() = Ds/(muEffByDelta + Ds);
}
else
{
// Lookup the neighbour momentum transport model
const momentumTransportModel& transportModelNbr =
mpp.nbrMesh().lookupType<momentumTransportModel>();
// Get the patch laminar viscosity
const tmp<scalarField> nuEffByDeltaNbr
// Get the neighbour patch laminar dynamic viscosity divided by delta
const tmp<scalarField> muEffByDeltaNbr
(
mpp.fromNeighbour
(
transportModelNbr.nuEff(patchiNbr)*patchNbr.deltaCoeffs()
rhopNbr*transportModelNbr.nuEff(patchiNbr)
*patchNbr.deltaCoeffs()
)
);
// Calculate the value-fraction from the balance between the
// external fluid and internal film stresses
valueFraction() = nuEffByDeltaNbr()/(nuEffByDelta + nuEffByDeltaNbr());
valueFraction() = muEffByDeltaNbr()/(muEffByDelta + muEffByDeltaNbr());
}
mixedFvPatchField<vector>::updateCoeffs();