fvOptions: limitVelocity: Fixed the correction factor

Resolves bug report https://bugs.openfoam.org/view.php?id=2634

fvOptions: limitVelocity: Fixed the correction factor
Resolves bug report https://bugs.openfoam.org/view.php?id=2634
bd4fcd16 · Will Bainbridge · Andrew Heather · 2e9fbcca · bd4fcd16
Commit bd4fcd16 authored 7 years ago by Will Bainbridge Committed by Andrew Heather 7 years ago
--- a/src/fvOptions/corrections/limitVelocity/limitVelocity.C
+++ b/src/fvOptions/corrections/limitVelocity/limitVelocity.C
@@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2016-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -94,7 +94,7 @@ void Foam::fv::limitVelocity::correct(volVectorField& U)

        if (magSqrUi > maxSqrU)
        {
-            Uif[celli] *= maxSqrU/magSqrUi;
+            Uif[celli] *= sqrt(maxSqrU/magSqrUi);
        }
    }

@@ -115,7 +115,7 @@ void Foam::fv::limitVelocity::correct(volVectorField& U)

                    if (magSqrUi > maxSqrU)
                    {
-                        Up[facei] *= maxSqrU/magSqrUi;
+                        Up[facei] *= sqrt(maxSqrU/magSqrUi);
                    }
                }
            }