Bugfix: set boundary velocity on a moving IB surface

src/immersedBoundary/immersedBoundaryDynamicMesh/movingImmersedBoundary/movingImmersedBoundary.C

@@ -77,6 +77,23 @@ void Foam::movingImmersedBoundary::movePoints() const
             << abort(FatalError);
     }
 
+    // Get non-const reference to velocity field
+    volVectorField& U = const_cast<volVectorField&>
+    (
+        mesh().lookupObject<volVectorField>("U")
+    );
+
+    // Get non-const reference to patch field
+    immersedBoundaryFvPatchVectorField& ibPatchField =
+        refCast<immersedBoundaryFvPatchVectorField>
+        (
+            U.boundaryField()[patchID]    
+        );
+
+    // Set refValue_ to moving boundary velocity
+    ibPatchField.refValue() = 
+        transform(sbmfPtr_->velocity(), refIbSurface_.points());
+
     const immersedBoundaryPolyPatch& cibPatch =
         refCast<const immersedBoundaryPolyPatch>
         (

src/immersedBoundary/immersedBoundaryDynamicMesh/movingImmersedBoundary/movingImmersedBoundary.H

@@ -39,6 +39,8 @@ SourceFiles
 #include "fvMesh.H"
 #include "solidBodyMotionFunction.H"
 #include "triSurfaceMesh.H"
+#include "fvPatchFields.H"
+#include "immersedBoundaryFvPatchFields.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //