CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   InterDyMFoam with rotating reference frame (https://www.cfd-online.com/Forums/openfoam-solving/98014-interdymfoam-rotating-reference-frame.html)

frank79 March 1, 2012 09:28

InterDyMFoam with rotating reference frame
 
Dear all,

I'm using InterDyMFoam to simulate a multiphase problem in which I need to refine dinamically the mesh according to the interface position as follows:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  1.6                                  |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    location    "constant";
    object      dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

//dynamicFvMesh  staticFvMesh;
dynamicFvMesh  dynamicRefineFvMesh;

dynamicRefineFvMeshCoeffs
{
    // How often to refine
    refineInterval  1;
  // Field to be refinement on
        field          alpha1;
    // Refine field inbetween lower..upper
    lowerRefineLevel 0.001;
    upperRefineLevel 0.999;
  // If value < unrefineLevel unrefine
    unrefineLevel  10;
  // Have slower than 2:1 refinement
    nBufferLayers  1;
        // Refine cells only up to maxRefinement levels
    maxRefinement  2;
        // Stop refinement if maxCells reached
    maxCells        2250000;
      // Flux field and corresponding velocity field. Fluxes on changed
    // faces get recalculated by interpolating the velocity.
        correctFluxes
            (
            (
                phi
            U
            )
                );
    // Write the refinement level as a volScalarField
      dumpLevel      true;
        }

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

Now I want to place my computational domain in a non inertial reference frame (i.e. rotation around z axis) keeping the dynamic mesh refinement due to the interface position.

Do you know how I can manage this issue?:confused:

Thanks in advance

Francesco:)

frank79 March 2, 2012 11:25

new code
 
I decided to hack the code of InterDyMFoam to add the Coriolis and the centrifugal term. At the moment I'm adding the SRFModel to the original InterDyMFoam. I'll let you know if it will work
bye

vonboett July 2, 2012 05:41

Hi Francesco,

I model a rotating drum partially filled with liquid, including internal faces using interDyMFoam, I use the solid-body motion function rotatingMotion and apply it to a cellZone containing the whole grid. Since the simulation always crushes after mesh.update() and I am stuck, I thought maybe you could post your fvSolution and fvSchemes dict because you seeem to be more successful?

Thanks,

Albrecht


All times are GMT -4. The time now is 18:25.