From 9833e614da2e4a8e19ad21bba23e07f381ae7e7b Mon Sep 17 00:00:00 2001
From: Andrew Heather <>
Date: Fri, 3 Dec 2021 11:58:03 +0000
Subject: [PATCH] BUG: surfaceAlignedSBRStressFvMotionSolver - corrected units
(see #2265)
Note - this resolves the reported dimensions error, but it is not clear
if this code (has ever?) run or if this change is appropriate...
---
.../surfaceAlignedSBRStressFvMotionSolver.C | 10 ++++------
1 file changed, 4 insertions(+), 6 deletions(-)
diff --git a/src/fvMotionSolver/fvMotionSolvers/displacement/surfaceAlignedSBRStress/surfaceAlignedSBRStressFvMotionSolver.C b/src/fvMotionSolver/fvMotionSolvers/displacement/surfaceAlignedSBRStress/surfaceAlignedSBRStressFvMotionSolver.C
index 098f4fb7a4e..a888fd225df 100644
--- a/src/fvMotionSolver/fvMotionSolvers/displacement/surfaceAlignedSBRStress/surfaceAlignedSBRStressFvMotionSolver.C
+++ b/src/fvMotionSolver/fvMotionSolvers/displacement/surfaceAlignedSBRStress/surfaceAlignedSBRStressFvMotionSolver.C
@@ -384,11 +384,9 @@ void Foam::surfaceAlignedSBRStressFvMotionSolver::solve()
sigmaD_ = magNewSigmaD;
}
- const surfaceScalarField Df
- (
- dimensionedScalar("viscosity", dimViscosity, 1.0)
- *diffusivity().operator()()
- );
+ const dimensionedScalar oneViscosity("viscosity", dimViscosity, 1.0);
+
+ const surfaceScalarField Df(oneViscosity*diffusivity().operator()());
pointDisplacement_.boundaryFieldRef().updateCoeffs();
@@ -417,7 +415,7 @@ void Foam::surfaceAlignedSBRStressFvMotionSolver::solve()
)
)
==
- fvc::div(sigmaD_)
+ oneViscosity*fvc::div(sigmaD_)
+ fvOptions(cellDisp)
);
--
GitLab