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

BUG: Correct division by zero when initial U = 0 in fvOption

Browse Source 
source directionalPressureGradientExplicitSource
This commit is contained in:
sergio
2018-12-06 09:04:08 -08:00
parent 942a7f8565
commit ff74469087
1 changed files with 2 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
src/fvOptions/sources/derived/directionalPressureGradientExplicitSource/directionalPressureGradientExplicitSource.C
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2015 OpenCFD Ltd. \\ / A nd | Copyright (C) 2015-2018 OpenCFD Ltd.
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -431,7 +431,7 @@ void Foam::fv::directionalPressureGradientExplicitSource::correct
{ {
label cellI = cells_[i]; label cellI = cells_[i];
const vector Ufnorm = UfCells[i]/mag(UfCells[i]); const vector Ufnorm = UfCells[i]/(mag(UfCells[i]) + SMALL);
const tensor D = rotationTensor(Ufnorm, flowDir_); const tensor D = rotationTensor(Ufnorm, flowDir_);
@ -452,7 +452,6 @@ void Foam::fv::directionalPressureGradientExplicitSource::correct
Info<< "UfCell:= " << UfCells[i] << "U : " << U[cellI] << endl; Info<< "UfCell:= " << UfCells[i] << "U : " << U[cellI] << endl;
} }
} }
writeProps(gradP0_ + dGradP_); writeProps(gradP0_ + dGradP_);
} }