Simulation crash with dynamicRefineFvMesh and kOmegaSST - OF 2.3.x

nathanael · June 28, 2014, 14:32

Hello,

I have some trouble with dynamicRefineFvMesh, the simulation crash just before the computation of omega:

Code:

Selected 11492 cells for refinement out of 60106.
Refined from 60106 to 140550 cells.
Selected 0 split points out of a possible 11492.
Execution time for mesh.update() = 1.24 s
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 6.51074e-05, No Iterations 11
time step continuity errors : sum local = 1.53975e-10, global = 1.4405e-10, cumulative = -1.14609e-07
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.808133  Min(alpha1) = 0  Max(alpha1) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.808132  Min(alpha1) = 0  Max(alpha1) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.808132  Min(alpha1) = 0  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.000136642, Final residual = 9.66145e-07, No Iterations 6
time step continuity errors : sum local = 2.40051e-07, global = -1.43622e-07, cumulative = -2.58231e-07
GAMG:  Solving for p_rgh, Initial residual = 0.000844249, Final residual = 6.92963e-06, No Iterations 3
time step continuity errors : sum local = 2.3079e-07, global = -4.7405e-08, cumulative = -3.05636e-07
GAMG:  Solving for p_rgh, Initial residual = 0.000373485, Final residual = 6.69911e-09, No Iterations 17
time step continuity errors : sum local = 1.53759e-10, global = 2.86627e-11, cumulative = -3.05607e-07
#0  Foam::error::printStack(Foam::Ostream&) in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/lib/libOpenFOAM.so"
#2  
 at sigaction.c:0
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/lib/libOpenFOAM.so"
#4  void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/lib/libincompressibleRASModels.so"
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/lib/libincompressibleRASModels.so"
#6  Foam::incompressible::RASModels::kOmegaSST::correct() in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/lib/libincompressibleRASModels.so"
#7  
 in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/bin/interDyMFoam"
#8  __libc_start_main in "/lib64/libc.so.6"
#9  
 in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/bin/interDyMFoam"
Floating point exception

The computation run fine without dynamic mesh and with kOmegaSST as well as with dynamic mesh and no turbulence...

Testing with OF-2.2.x the solver complain about missing flux mapping table:

Code:

--> FOAM Warning : 
    From function dynamicRefineFvMesh::refine(const labelList&)
    in file dynamicRefineFvMesh/dynamicRefineFvMesh.C at line 313
    Cannot find surfaceScalarField rho*phi in user-provided flux mapping table 
6
(
rhoPhi none
phi U
phiAbs_0 U_0
ghf none
nHatf none
phiAbs U
)

After adding the required flux in dynamicMeshDict the solver solve the first time step and crash at second time step when solving omega for OF-2.3x and p_rgh for OF-2.2.x.

Any help would be appreciated

Nathanaël.

nimasam · June 29, 2014, 04:40

here the information, you should noticed

Quote:

Cannot find surfaceScalarField rho*phi in user-provided flux mapping table

so instead of rhoPhi, you should recall rho*phi in dynamicMeshRefinement dictionary

nathanael · June 29, 2014, 04:49

Hello,

Thanks for the repply, I already add the required flux, here is my dynamicMeshDict:

Code:

dynamicFvMesh   dynamicRefineFvMesh;

dynamicRefineFvMeshCoeffs
{
    // How often to refine
    refineInterval  1;
    
    // Field to be refinement on
    field           alpha.water;
    
    // Refine field inbetween lower..upper
    lowerRefineLevel 0.001;
    upperRefineLevel 0.999;
    
    // If value < unrefineLevel unrefine
    unrefineLevel   10;
    
    // Have slower than 2:1 refinement
    nBufferLayers   2;
    
    // Refine cells only up to maxRefinement levels
    maxRefinement   2;
    
    // Stop refinement if maxCells reached
    maxCells        2000000;
    
    // Flux field and corresponding velocity field. Fluxes on changed
    // faces get recalculated by interpolating the velocity. Use 'none'
    // on surfaceScalarFields that do not need to be reinterpolated.
    correctFluxes
    (
        (phi U)
        (phiAbs U)
        (phiAbs_0 U_0)
       // (phi none)
        (nHatf none)
        (rhoPhi none)
        (ghf none)
        
        (phi_0 none)
        (rho*phi none) 
        
    );
    // Write the refinement level as a volScalarField
    dumpLevel       false;
}

And here the log when running interDyMFoam:

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.x-896bd3fd1bcf
Exec   : interDyMFoam
Date   : Jun 29 2014
Time   : 09:47:51
Host   : "localhost.localdomain"
PID    : 3681
Case   : /mnt/hgfs/Data/Perso/Recherche/Simulation/CrashDynamicMesh
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicRefineFvMesh

PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Using dynamicCode for patch Inlet on field alpha.water at line 50 in "/mnt/hgfs/Data/Perso/Recherche/Simulation/CrashDynamicMesh/0/alpha.water.boundaryField.Inlet"
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
    alphaK1         0.85034;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.85616;
    gamma1          0.5532;
    gamma2          0.4403;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
}


Reading g
Calculating field g.h

No finite volume options present

GAMG:  Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Reading/calculating face velocity Uf

Courant Number mean: 0 max: 0

Starting time loop

Interface Courant Number mean: 0 max: 0
Courant Number mean: 0 max: 0
deltaT = 1.19976e-05
Time = 1.19976e-05

Selected 2874 cells for refinement out of 39988.
Refined from 39988 to 60106 cells.
Selected 0 split points out of a possible 2874.
Execution time for mesh.update() = 0.47 s
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 4.40804e-05, No Iterations 20
time step continuity errors : sum local = 3.82811e-11, global = 2.15865e-11, cumulative = 2.15865e-11
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.808134  Min(alpha1) = 0  Max(alpha1) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.808133  Min(alpha1) = 0  Max(alpha1) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.808133  Min(alpha1) = 0  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.00561127, No Iterations 6
time step continuity errors : sum local = 2.20727e-07, global = -1.2877e-07, cumulative = -1.28749e-07
GAMG:  Solving for p_rgh, Initial residual = 4.4142e-05, Final residual = 2.86487e-07, No Iterations 3
time step continuity errors : sum local = 1.05078e-07, global = 1.39143e-08, cumulative = -1.14834e-07
GAMG:  Solving for p_rgh, Initial residual = 9.41247e-06, Final residual = 3.9572e-09, No Iterations 7
time step continuity errors : sum local = 1.48954e-09, global = 8.13235e-11, cumulative = -1.14753e-07
smoothSolver:  Solving for omega, Initial residual = 0.00226246, Final residual = 7.42652e-09, No Iterations 7
bounding omega, min: -0.0235881 max: 12560.6 average: 8.45464
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 3.03672e-09, No Iterations 3
bounding k, min: -1.82238e-06 max: 0.000106311 average: 9.96208e-05
ExecutionTime = 3.23 s  ClockTime = 3 s

Interface Courant Number mean: 0 max: 0
Courant Number mean: 7.20312e-05 max: 0.0206844
deltaT = 1.43902e-05
Time = 2.63878e-05

Selected 11492 cells for refinement out of 60106.
Refined from 60106 to 140550 cells.
Selected 0 split points out of a possible 11492.
Execution time for mesh.update() = 1.19 s
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 6.51074e-05, No Iterations 11
time step continuity errors : sum local = 1.53975e-10, global = 1.4405e-10, cumulative = -1.14609e-07
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.808133  Min(alpha1) = 0  Max(alpha1) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.808132  Min(alpha1) = 0  Max(alpha1) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.808132  Min(alpha1) = 0  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.000136642, Final residual = 9.66145e-07, No Iterations 6
time step continuity errors : sum local = 2.40051e-07, global = -1.43622e-07, cumulative = -2.58231e-07
GAMG:  Solving for p_rgh, Initial residual = 0.000844249, Final residual = 6.92963e-06, No Iterations 3
time step continuity errors : sum local = 2.3079e-07, global = -4.7405e-08, cumulative = -3.05636e-07
GAMG:  Solving for p_rgh, Initial residual = 0.000373485, Final residual = 6.69911e-09, No Iterations 17
time step continuity errors : sum local = 1.53759e-10, global = 2.86627e-11, cumulative = -3.05607e-07
#0  Foam::error::printStack(Foam::Ostream&) in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/lib/libOpenFOAM.so"
#2  
 at sigaction.c:0
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/lib/libOpenFOAM.so"
#4  void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/lib/libincompressibleRASModels.so"
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/lib/libincompressibleRASModels.so"
#6  Foam::incompressible::RASModels::kOmegaSST::correct() in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/lib/libincompressibleRASModels.so"
#7  
 in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/bin/interDyMFoam"
#8  __libc_start_main in "/lib64/libc.so.6"
#9  
 in "/sw/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64Gcc48DPOpt/bin/interDyMFoam"
Floating point exception

nimasam · June 29, 2014, 07:08

it seems it returns to omega equation, it divides to zero

, check your BCs
can you run your problem with out mesh refinement?

nathanael · June 29, 2014, 18:02

Thanks for the interest in my problem.

Yes, I can run my case without refinement.

June 29, 2014, 07:08		#4
nimasam Senior Member Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,267 Blog Entries: 1 Rep Power: 25	it seems it returns to omega equation, it divides to zero , check your BCs can you run your problem with out mesh refinement? Aabadani likes this. __________________ My Personal Website (http://nimasamkhaniani.ir/) Telegram channel (https://t.me/cfd_foam)

June 29, 2014, 18:02		#5
nathanael New Member Nathanaël Geng Join Date: May 2012 Posts: 18 Rep Power: 14	Thanks for the interest in my problem. Yes, I can run my case without refinement.