CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Shot Sleeve Wave Dynamics using InterFoam. Interface capturing at corner areas.

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 13, 2024, 04:56
Default Shot Sleeve Wave Dynamics using InterFoam. Interface capturing at corner areas.
  #1
Member
 
Michael Sukham
Join Date: Mar 2020
Location: India
Posts: 79
Rep Power: 6
2538sukham is on a distinguished road
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.

Code:
/*--------------------------------*- 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:

Code:
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.
Attached Images
File Type: jpg setFields.jpg (23.7 KB, 4 views)
File Type: png Screenshot 2023-12-04 170933.png (9.3 KB, 7 views)
2538sukham is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Wind turbine simulation Saturn CFX 58 July 3, 2020 01:13
Radiation interface hinca CFX 15 January 26, 2014 17:11
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 17:44
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 04:05
Replace periodic by inlet-outlet pair lego CFX 3 November 5, 2002 20:09


All times are GMT -4. The time now is 22:39.