From ff744690875284b9bd020d3e1a6bbffa80d2bdca Mon Sep 17 00:00:00 2001 From: sergio Date: Thu, 6 Dec 2018 09:04:08 -0800 Subject: [PATCH] BUG: Correct division by zero when initial U = 0 in fvOption source directionalPressureGradientExplicitSource --- .../directionalPressureGradientExplicitSource.C | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/src/fvOptions/sources/derived/directionalPressureGradientExplicitSource/directionalPressureGradientExplicitSource.C b/src/fvOptions/sources/derived/directionalPressureGradientExplicitSource/directionalPressureGradientExplicitSource.C index 9d6adf3027..91a3e495a1 100644 --- a/src/fvOptions/sources/derived/directionalPressureGradientExplicitSource/directionalPressureGradientExplicitSource.C +++ b/src/fvOptions/sources/derived/directionalPressureGradientExplicitSource/directionalPressureGradientExplicitSource.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2015 OpenCFD Ltd. + \\ / A nd | Copyright (C) 2015-2018 OpenCFD Ltd. \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -431,7 +431,7 @@ void Foam::fv::directionalPressureGradientExplicitSource::correct { 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_); @@ -452,7 +452,6 @@ void Foam::fv::directionalPressureGradientExplicitSource::correct Info<< "UfCell:= " << UfCells[i] << "U : " << U[cellI] << endl; } } - writeProps(gradP0_ + dGradP_); }