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

VoFPatchTransfer: Rationalised the names of the cached fields and their origin

Browse Source 
required rationalisation of the names of variables and functions in
surfaceFilmModels to clarify which relate to the film region.
This commit is contained in:
Henry Weller
2023-03-09 13:58:08 +00:00
parent 0b0957f03d
commit 5b9fe57c23
25 changed files with 311 additions and 263 deletions
Download Patch File Download Diff File Expand all files Collapse all files

79
applications/solvers/modules/compressibleVoF/fvModels/VoFSurfaceFilm/VoFPatchTransfer/VoFPatchTransfer.C
View File

@ -169,8 +169,6 @@ void VoFPatchTransfer::correct
const volScalarField& alphaVoF = thermo.alpha1();
const volScalarField& rhoVoF = thermo.thermo1().rho()();
const volScalarField& heVoF = thermo.thermo1().he();
const volScalarField& TVoF = thermo.thermo1().T();
const volScalarField CpVoF(thermo.thermo1().Cp());
forAll(patchIDs_, pidi)
{
@ -189,36 +187,44 @@ void VoFPatchTransfer::correct
if (primaryPatchi != -1)
{
scalarField deltaCoeffs
const scalarField deltaCoeffsVoF
(
film.primaryMesh().boundary()[primaryPatchi].deltaCoeffs()
film.toFilm
(
patchi,
film.primaryMesh().boundary()[primaryPatchi].deltaCoeffs()
)
);
film.toRegion(patchi, deltaCoeffs);
scalarField alphap(alphaVoF.boundaryField()[primaryPatchi]);
film.toRegion(patchi, alphap);
scalarField rhop(rhoVoF.boundaryField()[primaryPatchi]);
film.toRegion(patchi, rhop);
vectorField Up(UVoF.boundaryField()[primaryPatchi]);
film.toRegion(patchi, Up);
scalarField hp(heVoF.boundaryField()[primaryPatchi]);
film.toRegion(patchi, hp);
scalarField Tp(TVoF.boundaryField()[primaryPatchi]);
film.toRegion(patchi, Tp);
scalarField Cpp(CpVoF.boundaryField()[primaryPatchi]);
film.toRegion(patchi, Cpp);
scalarField Vp
const scalarField alphaVoFp
(
film.primaryMesh().boundary()[primaryPatchi]
.patchInternalField(film.primaryMesh().V())
film.toFilm(patchi, alphaVoF.boundaryField()[primaryPatchi])
);
const scalarField rhoVoFp
(
film.toFilm(patchi, rhoVoF.boundaryField()[primaryPatchi])
);
const vectorField UVoFp
(
film.toFilm(patchi, UVoF.boundaryField()[primaryPatchi])
);
const scalarField heVoFp
(
film.toFilm(patchi, heVoF.boundaryField()[primaryPatchi])
);
const scalarField VVoFp
(
film.toFilm
(
patchi,
film.primaryMesh().boundary()[primaryPatchi]
.patchInternalField(film.primaryMesh().V())
)
);
film.toRegion(patchi, Vp);
const polyPatch& pp = pbm[patchi];
const labelList& faceCells = pp.faceCells();
@ -235,10 +241,11 @@ void VoFPatchTransfer::correct
scalar dMass = 0;
// Film->VoF transfer
if
(
delta[celli] > 2*deltaFactorToVoF_/deltaCoeffs[facei]
|| alphap[facei] > alphaToVoF_
delta[celli] > 2*deltaFactorToVoF_/deltaCoeffsVoF[facei]
|| alphaVoFp[facei] > alphaToVoF_
)
{
dMass =
@ -249,19 +256,21 @@ void VoFPatchTransfer::correct
energyToTransfer[celli] += dMass*he[celli];
}
// VoF->film transfer
if
(
alphap[facei] > 0
&& delta[celli] < 2*deltaFactorToFilm_/deltaCoeffs[facei]
&& alphap[facei] < alphaToFilm_
alphaVoFp[facei] > 0
&& delta[celli] < 2*deltaFactorToFilm_/deltaCoeffsVoF[facei]
&& alphaVoFp[facei] < alphaToFilm_
)
{
dMass =
-transferRateCoeff_*alphap[facei]*rhop[facei]*Vp[facei];
-transferRateCoeff_
*alphaVoFp[facei]*rhoVoFp[facei]*VVoFp[facei];
massToTransfer[celli] += dMass;
momentumToTransfer[celli] += dMass*Up[facei];
energyToTransfer[celli] += dMass*hp[facei];
momentumToTransfer[celli] += dMass*UVoFp[facei];
energyToTransfer[celli] += dMass*heVoFp[facei];
}
availableMass[celli] -= dMass;