exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Minor updates to oversetCourantNo.H and oversetAlphaCourantNo.H

Browse source 
Just to make it more self-sustainable
This commit is contained in:
Vuko Vukcevic 2019-05-15 17:47:07 +02:00 committed by Vuko Vukcevic
parent 189078e98e
commit f041b8e746
2 changed files with 16 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

17
src/overset/oversetMesh/include/oversetAlphaCourantNo.H
View file

@ -38,21 +38,28 @@ scalar maxAlphaCo
scalar alphaCoNum = 0.0; scalar alphaCoNum = 0.0;
scalar meanAlphaCoNum = 0.0; scalar meanAlphaCoNum = 0.0;
surfaceScalarField alpha1f = const surfaceScalarField alpha1f =
fvc::interpolate(min(max(alpha1, scalar(0)), scalar(1))); fvc::interpolate(min(max(alpha1, scalar(0)), scalar(1)));
const dimensionedScalar alphaOffset("alphaOffset", dimless, dAlpha); const dimensionedScalar alphaOffset
(
"alphaOffset",
dimless,
runTime.controlDict().lookupOrDefault("dAlpha", 0.01)
);
if (mesh.nInternalFaces()) if (mesh.nInternalFaces())
{ {
surfaceScalarField magAlphaPhi const oversetMesh& om = oversetMesh::New(mesh);
const surfaceScalarField magAlphaPhi
( (
pos(alpha1f - alphaOffset)* pos(alpha1f - alphaOffset)*
pos(scalar(1) - alphaOffset - alpha1f)* pos(scalar(1) - alphaOffset - alpha1f)*
mag(faceOversetMask*phi) mag(om.sGamma()*phi)
); );
surfaceScalarField SfUfbyDelta = const surfaceScalarField SfUfbyDelta =
mesh.surfaceInterpolation::deltaCoeffs()*magAlphaPhi; mesh.surfaceInterpolation::deltaCoeffs()*magAlphaPhi;
const scalar deltaT = runTime.deltaT().value(); const scalar deltaT = runTime.deltaT().value();

6
src/overset/oversetMesh/include/oversetCourantNo.H
View file

@ -36,9 +36,11 @@ scalar velMag = 0.0;
if (mesh.nInternalFaces()) if (mesh.nInternalFaces())
{ {
surfaceScalarField magPhi = mag(faceOversetMask*phi); const oversetMesh& om = oversetMesh::New(mesh);
surfaceScalarField SfUfbyDelta = const surfaceScalarField magPhi = mag(om.sGamma()*phi);
const surfaceScalarField SfUfbyDelta =
mesh.surfaceInterpolation::deltaCoeffs()*magPhi; mesh.surfaceInterpolation::deltaCoeffs()*magPhi;
CoNum = max(SfUfbyDelta/mesh.magSf()) CoNum = max(SfUfbyDelta/mesh.magSf())