[Arif14]

Dear CFD community,

I am trying to implement what seems to be a simple setup. I have an axisymmetric geometry (a section of a cylinder, as shown in the attached image). Inside this geometry, I have a mixture of water and another liquid-like phase, and I am using interIsoFoam for my simulation.
I would like to introduce a pre-conditioned state representing an increase in the diameter of the geometry. My approach was to first run a simulation using the moveDynamicMesh solver, applying a displacement to the top wall to increase the diameter. Then, I intended to save the updated coordinates and field data (velocity, air/liquid ratio, etc.) and use them as the initial condition for a second simulation with my original setup.

However, I have encountered an issue: while the top wall moves as expected, the side mesh does not properly follow the displacement, resulting in a highly distorted mesh (as seen in the attached figure). I searched for similar cases in this forum and OpenFOAM tutorials but could not find a satisfactory example.

Below you will find my dynamicMeshDict and pointDisplacement files. I would greatly appreciate any help or advice regarding my approach!

Thank you so much!

Best regards,

dynamicFvMesh   dynamicMotionSolverFvMesh;

motionSolverLibs ("libfvMotionSolvers.so");

solver displacementLaplacian;

diffusivity quadratic inverseDistance (topwall);

dimensions      [0 1 0 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type fixedValue;
        $internalField;
    }
    outlet
    {
        type fixedValue;
        $internalField;
    }
    topwall
    {
        type uniformFixedValue;
        uniformValue constant (0 0.002 0);  // Moves up by 0.002m
    }
    back
    {
        type wedge;
    }
    front
    {
        type wedge;
    }
    defaultFaces
    {
        type empty;
    }
}


// ************************************************************************* //