From 3eec2a49090195202da23a61da3adc2e6176cce3 Mon Sep 17 00:00:00 2001 From: David Blom Date: Sat, 30 Apr 2016 17:13:36 +0200 Subject: [PATCH] Bugfix myMovingWallVelocity boundary condition: the face centers of the (#273) current mesh should be used at each FSI outer iteration. Fixes #272 --- .../myMovingWallVelocityFvPatchVectorField.C | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/boundaryConditions/myMovingWallVelocity/myMovingWallVelocityFvPatchVectorField.C b/src/boundaryConditions/myMovingWallVelocity/myMovingWallVelocityFvPatchVectorField.C index 8a498bb5..f3bb21fc 100644 --- a/src/boundaryConditions/myMovingWallVelocity/myMovingWallVelocityFvPatchVectorField.C +++ b/src/boundaryConditions/myMovingWallVelocity/myMovingWallVelocityFvPatchVectorField.C @@ -99,10 +99,11 @@ namespace Foam { oldoldFc_ = oldFc_; oldFc_ = Fc_; - Fc_ = pp.faceCentres(); myTimeIndex_ = mesh.time().timeIndex(); } + Fc_ = pp.faceCentres(); + // const pointField& oldPoints = mesh.oldPoints(); const volVectorField & U = mesh.lookupObject( dimensionedInternalField().name() );