Make the test pass first.

Author: zilonglicfd

src/adjoint/DASolver/DAPimpleDyMFoam/DAPimpleDyMFoam.C

@@ -67,15 +67,11 @@ label DAPimpleDyMFoam::solvePrimal()
 
     Foam::argList& args = argsPtr_();
 #include "createTime.H"
-//#include "createDynamicFvMesh.H"
-#include "createMesh.H"
-    mesh.moving(true);
+#include "createDynamicFvMesh.H"
 #include "initContinuityErrs.H"
-//#include "createDyMControls.H"
-#include "myCreateDyMControls.H"
+#include "createDyMControls.H"
 #include "createFieldsPimpleDyM.H"
-//#include "createUfIfPresent.H"
-#include "myCreateUfIfPresent.H"
+#include "createUfIfPresent.H"
 #include "CourantNo.H"
 #include "setInitialDeltaT.H"
 
@@ -118,7 +114,6 @@ label DAPimpleDyMFoam::solvePrimal()
                 );
 
                 mesh.movePoints(readPoints);
-                U.correctBoundaryConditions(); 
 
                 if (mesh.changing())
                 {

src/adjoint/DASolver/DAPimpleDyMFoam/DAPimpleDyMFoam.H

@@ -39,26 +39,9 @@
 #include "turbulentTransportModel.H"
 #include "pimpleControlDF.H"
 #include "pimpleControl.H"
-//#include "dynamicFvMesh.H"
+#include "dynamicFvMesh.H"
 #include "CorrectPhi.H"
 
-void myCorrectUf
-(
-    autoPtr<surfaceVectorField>& Uf,
-    const volVectorField& U,
-    const surfaceScalarField& phi
-)
-{
-    const fvMesh& mesh = U.mesh();
-
-    if (mesh.moving())
-    {
-        Uf() = fvc::interpolate(U);
-        surfaceVectorField n(mesh.Sf()/mesh.magSf());
-        Uf() += n*(phi/mesh.magSf() - (n & Uf()));
-    }
-}
-
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 namespace Foam

src/adjoint/DASolver/DAPimpleDyMFoam/myCreateDyMControls.H

@@ -58,8 +58,7 @@ U.correctBoundaryConditions();
 fvOptions.correct(U);
 
 // Correct Uf if the mesh is moving
-//fvc::correctUf(Uf, U, phi);
-myCorrectUf(Uf, U, phi);
+fvc::correctUf(Uf, U, phi);
 
 // Make the fluxes relative to the mesh motion
 fvc::makeRelative(phi, U);

src/adjoint/DASolver/DAPimpleDyMFoam/myCreateUfIfPresent.H

@@ -12,8 +12,6 @@
 gcomm = MPI.COMM_WORLD
 
 os.chdir("./reg_test_files-main/NACA0012DynamicMeshV4")
-if gcomm.rank == 0:
-    os.system("rm -rf processor* *.bin")
 
 # aero setup
 U0 = 10.0