[2538sukham]

Hello and greetings. I am simulating the wave dynamics in shot sleeve of high pressure die casting. I have the meshing done by SnappyHexMesh which produced a hex dominated mesh with some tet on the edge of the cylinder. Mesh motion is achieved by velocityComponentLaplacian since I have a moving plunger only in the x-direction. The motion is done by codedFixedValue in pointMotionUx.

My problem is when the liquid reached the corner of the cylinder at top. the interface by using MULES cAlpha 1 gives some oscillations.

Making cAlpha = 0 smears the interface but it does give a very nice interface by setting iso value at 0.5 in post-processing. I have some way in my mind

1. Refining the edge cells using a higher refinement of feature edge.

2. Using isoAdvector. (Although I am not sure of the settings with dynamic mesh)
3. Using cAlpha 0 and refining the feature edge cells.

I had some success with mapFields by extracting the stl file from latestTime using surfaceMeshExtract -latestTime and remeshing using snappy and mapping on the new mesh the field values. But again when ever the interface touches the upper corner, the problem persist. InterfaceCapturing on these corner points seems to give me the errors.

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2306                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       pointScalarField;
    object      pointMotionUx;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    movingWall
    {
        /*type            uniformFixedValue;
        uniformValue    constant -0.2;
        */
        type            codedFixedValue;
        value           uniform 0;      
        name            increaseToFixedValue;
        
        code 
        #{
          const scalar t = this->db().time().value();
          const scalar a = 1.8;
          const scalar b = 0.0254;
          operator==(a*b*(exp(a*t)-1));            
        #};
        
    }

    sleeveWall
    {
        type            slip;
        
    }

    axiSymmetric
    {
        type            slip;
        
    }

    endWall
    {
        type            uniformFixedValue;
        uniformValue    constant 0;
    }
    outletWall
    {
        type            uniformFixedValue;
        uniformValue    constant 0;
    }
    outlet
    {
        type            uniformFixedValue;
        uniformValue    constant 0;
    }
}

fvSchemes:

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    div(rhoPhi,U)   Gauss linearUpwind grad(U);                        // Gauss vanLeerV
    div(phi,alpha)  Gauss interfaceCompression limitedLinear01;     //Gauss vanLeer;
    div(phirb,alpha) Gauss interfaceCompression;            //Gauss vanLeer;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

If anyone has done some simulations, I would appreciate some tips on getting a academically viable simulation result.