Commit 1f59360e authored by Will Bainbridge's avatar Will Bainbridge Committed by Andrew Heather
Browse files

fvOptions: limitVelocity: Fixed the correction factor

Resolves bug report https://bugs.openfoam.org/view.php?id=2634
parent 3b852844
...@@ -2,7 +2,7 @@ ...@@ -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) 2016 OpenFOAM Foundation \\ / A nd | Copyright (C) 2016-2017 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
...@@ -94,7 +94,7 @@ void Foam::fv::limitVelocity::correct(volVectorField& U) ...@@ -94,7 +94,7 @@ void Foam::fv::limitVelocity::correct(volVectorField& U)
if (magSqrUi > maxSqrU) if (magSqrUi > maxSqrU)
{ {
Uif[celli] *= maxSqrU/magSqrUi; Uif[celli] *= sqrt(maxSqrU/magSqrUi);
} }
} }
...@@ -115,7 +115,7 @@ void Foam::fv::limitVelocity::correct(volVectorField& U) ...@@ -115,7 +115,7 @@ void Foam::fv::limitVelocity::correct(volVectorField& U)
if (magSqrUi > maxSqrU) if (magSqrUi > maxSqrU)
{ {
Up[facei] *= maxSqrU/magSqrUi; Up[facei] *= sqrt(maxSqrU/magSqrUi);
} }
} }
} }
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment