[Daizong]

Dear foamers,

I am a new OpenFOAM user (using openfoam2212) and now simulating the irregular sphere settling process with overPimpleDyMFoam. The mass is 0.000631573kg, volume is 6.31569E-07 m3, and density of sphere is 1000.007029 kg/m3. The water density is 1000 kg/m3

In 0 folder I set the initial condition (from the tutorial https://gitlab.com/mAlletto/openfoam...ttelingGravity) but set the turbulence model as RAS in constant/turbulenceProperties as following. I also attached the dynamicFvMeshDict,transportproperties, fvsolution, and fvschemes followed the error report.

My problem:

I tried the laminar model at beginning because it did not include the initial conditions for turbulence, but the linear velocity increased remarkably and the deltaT decreased to pretty small value with the simulation so I used the RAS model then. as showed in error report, there were time step continuity errors. I have no idea how to deal with it. Any help would be highly appreciated. (I guess the possible reasons are from the initial condition of turbulence).

initial condition
RAS
{
RASModel RNGkEpsilon;
turbulence on;
printcoeffs on;
}

flowVelocity (10 0 0);
pressure 0;
turbulentKE 1.5;
turbulentEpsilon 0.88;



error report
Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicOversetFvMesh
Selecting motion solver: sixDoFRigidBodyMotion
Applying motion to entire mesh
Selecting sixDoFSolver Newmark
Translational constraint tensor (1 0 0 0 1 0 0 0 1)
Rotational constraint tensor (1 0 0 0 1 0 0 0 1)

PIMPLE: no residual control data found. Calculations will employ 10 corrector loops

Reading field p

Reading field U

Reading/calculating face flux field phi

Creating cellMask field to block out hole cells

--> FOAM Warning :
From bool Foam:versetPolyPatch::master() const
in file oversetPolyPatch/oversetPolyPatch.C at line 145
The master overset patch is not the first patch. Generally the first patch should be an overset patch to guarantee consistent operation.
Creating interpolatedCells field

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model RNGkEpsilon
Reading/calculating face velocity Uf

No MRF models present

No finite volume options present
Courant Number mean: 0 max: 0

Starting time loop

forces forces1:
rho: rhoInf
Freestream density (rhoInf) set to 1000
Not including porosity effects

timeActivatedFileUpdate fileUpdate1 output:
time vs file list:
-1 "<system>/maxCo0.01"
0.01 "<system>/maxCo0.02"
0.04 "<system>/maxCo0.03"
0.06 "<system>/maxCo0.04"
0.07 "<system>/maxCo0.05"
0.09 "<system>/maxCo0.06"
0.1 "<system>/maxCo0.07"
0.2 "<system>/maxCo0.08"
0.3 "<system>/maxCo0.09"
0.4 "<system>/maxCo0.10"

timeActivatedFileUpdate: copying file
from: "<system>/maxCo0.01"
to : "<system>/maxCo"


From virtual bool Foam::regIOobject::readIfModified()
in file db/regIOobject/regIOobjectRead.C at line 277
Re-reading object controlDict from file "/mnt/e/sphereSettelingGravity2212_MP/sphereAndBackground/system/controlDict" because of modified file "/mnt/e/sphereSettelingGravity2212_MP/sphereAndBackground/system/maxCo"
Courant Number mean: 0 max: 0
deltaT = 0.0030303030303
Time = 0.0030303

forces forces:
rho: rhoInf
Freestream density (rhoInf) set to 1000
Not including porosity effects

6-DoF rigid body motion
Centre of rotation: (0 0 -9.00826446281e-06)
Centre of mass: (0 0 -9.00826446281e-06)
Orientation: (1 0 0 0 1 0 0 0 1)
Linear velocity: (0 0 -0.00594545454545)
Angular velocity: (0 0 0)
inverseDistance : detected 2 mesh regions
zone:0 nCells:192000 voxels59 59 59) bb-0.020000132665 -0.020000132665 -0.100000132665) (0.020000132665 0.020000132665 0.020000132665)
zone:1 nCells:18692 voxels59 59 59) bb-0.0150000519615 -0.0150000519615 -0.015009060226) (0.0150000519615 0.0150000519615 0.0149910436971)

Number of layers : 4
Average volumetric ratio : 4(16.6782887577 13.0753043973 11.650122156 10.0437345926)
Number of holes cells : 4(707 730 761 791)
Using layer : 2

Overset analysis : nCells : 210692
calculated : 207884
interpolated : 1120 (from local:1120 mixed local/remote:0 remote:0)
hole : 1688

PIMPLE: iteration 1
DILUPBiCGStab: Solving for p, Initial residual = 1, Final residual = 0.00773777695833, No Iterations 27
time step continuity errors : sum local = 8.56084336729e-08, global = -2.65620218706e-08, cumulative = -2.65620218706e-08
smoothSolver: Solving for epsilon, Initial residual = 0.129768588398, Final residual = 0.00129051453533, No Iterations 333
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.082583708451, No Iterations 1000
PIMPLE: iteration 2
DILUPBiCGStab: Solving for p, Initial residual = 0.841210113625, Final residual = 0.00740776961625, No Iterations 17
time step continuity errors : sum local = 0.00327267442251, global = -8.38237268899e-05, cumulative = -8.38502889118e-05
smoothSolver: Solving for epsilon, Initial residual = 0.17715554218, Final residual = 0.000159153340508, No Iterations 1
bounding epsilon, min: -578.957245995 max: 1669.09409763 average: 27.6035774427
smoothSolver: Solving for k, Initial residual = 0.0108869370046, Final residual = 3.78940422443e-05, No Iterations 2
PIMPLE: iteration 3
DILUPBiCGStab: Solving for p, Initial residual = 0.894907825867, Final residual = 0.00885006331495, No Iterations 50
time step continuity errors : sum local = 0.0002380157875, global = -2.27954502573e-05, cumulative = -0.000106645739169
smoothSolver: Solving for epsilon, Initial residual = 0.549210596174, Final residual = 0.0153268561057, No Iterations 1000
smoothSolver: Solving for k, Initial residual = 0.918117326688, Final residual = 0.00908047618361, No Iterations 30
PIMPLE: iteration 4
DILUPBiCGStab: Solving for p, Initial residual = 0.838254217105, Final residual = 0.00663166455272, No Iterations 22
time step continuity errors : sum local = 0.0656097180809, global = -0.0338107754635, cumulative = -0.0339174212027
#0 Foam::error:rintStack(Foam::Ostream&) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#1 Foam::sigFpe::sigHandler(int) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#7 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#8 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#9 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/liboverset.so
#10 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/liboverset.so
#11 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so
#12 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so
#13 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/bin/overPimpleDyMFoam
#14 ? in /lib/x86_64-linux-gnu/libc.so.6
#15 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#16 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/bin/overPimpleDyMFoam
Floating point exception

dynamicFvMeshDict

dynamicFvMesh dynamicOversetFvMesh;

dynamicOversetFvMeshCoeffs
{
// layerRelax 0.3;
}

solver sixDoFRigidBodyMotion;
accelerationRelaxation 0.4;

sixDoFRigidBodyMotionCoeffs
{

patches (combine_69);
innerDistance 100.;
outerDistance 101.;

accelerationRelaxation 0.4;
accelerationDamping 0.4;

// mass reduced by Buoyancy: (rhosolid - rhofluid ) * VolumeSphere
mass 4.44E-09;
centreOfMass ( 0 0 0 );
momentOfInertia ( 5.30E-14 5.30E-14 5.30E-14);
rho rhoInf;
rhoInf 1000;
g (0 0 -9.81);
report on;
reportToFile on;

solver
{
type Newmark;
}

constraints
{

}
}


transportProperties

DT DT [ 0 2 -1 0 0 0 0 ] 1; //4e-05;

transportModel Newtonian;

nu nu [ 0 2 -1 0 0 0 0 ] 1e-06;
rho rho [ 1 -3 0 0 0 0 0 ] 1000.0;

CrossPowerLawCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
m m [ 0 0 1 0 0 0 0 ] 1;
n n [ 0 0 0 0 0 0 0 ] 1;
}

BirdCarreauCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
k k [ 0 0 1 0 0 0 0 ] 0;
n n [ 0 0 0 0 0 0 0 ] 1;
}


Fvsolution

solvers
{
cellDisplacement
{
solver PCG;
preconditioner DIC;

tolerance 1e-06;
relTol 0;
maxIter 100;
}

p
{
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-5;
relTol 0.01;
maxIter 10;
}

pFinal
{
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-7;
relTol 0.01;
}

pcorr
{
$pFinal;
solver PCG;
preconditioner DIC;
tolerance 1e-7;
relTol 0.02;
}

pcorrFinal
{
$pcorr;
relTol 0.1;
tolerance 1e-7;
}

"(U|k|epsilon)"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-6;
relTol 0.01;
}

"(U|k|epsilon)Final"
{
$U;
tolerance 1e-6;
relTol 0;
}
}

PIMPLE
{
momentumPredictor no;
oversetAdjustPhi true;
nOuterCorrectors 10;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
massFluxInterpolation true;
ddtCorr false;
turbOnFinalIterOnly false;
moveMeshOuterCorrectors yes;
}

relaxationFactors
{
fields
{
p 1;
}
equations
{
"(U|k|omega)" 1;
".*Final" 1.0;
}
}

fvschemes

ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss limitedLinearV 1;
div(phi,k) Gauss limitedLinear 1;;
div(phi,epsilon) Gauss limitedLinear 1;;
div(phi,R) Gauss limitedLinear 1;;
div(R) Gauss limitedLinear 1;;
div(phi,omega) Gauss limitedLinear 1;
div(phid,p) Gauss limitedLinear 1;
div(phi,K) Gauss limitedLinear 1;
div(phi,e) Gauss limitedLinear 1;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected ;
}

oversetInterpolation
{
method inverseDistance;
holeLayers 4;
useLayer 2;
}

fluxRequired
{
default no;
pcorr ;
p ;
}