[frank79]

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:

/*--------------------------------*- 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?

Thanks in advance

Francesco

[frank79]

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]

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