[nicoela94]

Hello everyone,

I'm trying to run CFD simulation in order to get the lift coefficient and the drag coefficient of a NACA0018.

I use Pointwise as a meshing tool and the mesh I obtained did pass the checkMesh test. However I'm getting the following error message :

Create time

Create mesh for time = 0


SIMPLE: No convergence criteria found

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
bounding k, min: 0 max: 1000 average: 0.00375
RAS
{
model kOmegaSST;
turbulence on;
printCoeffs on;
alphaK1 0.85;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.856;
gamma1 0.555556;
gamma2 0.44;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
}

No MRF models present

No fvModels present
No fvConstraints present

Starting time loop

forces forces:
Not including porosity effects
forces forces write:
sum of forces:
pressure : (0 0 0)
viscous : (0 0 0)
porous : (0 0 0)
sum of moments:
pressure : (0 0 0)
viscous : (0 0 0)
porous : (0 0 0)

Time = 1

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.00271706, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.00541297, No Iterations 1
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#8 Foam::fvMatrix<double>::solve() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#9 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
Floating point exception

I guess this is due to bad boundary conditions but I cannot find where the error comes from.

Here are some more informations about my simulation:
- I use kOmegaSST model.
- y+ is about 1
- Reynolds number = 1e5
- I'm using a blunt trailing edge, the chord length is 0,1 m and the far field velocity is 10 m/s
- Supposing turbulence intensity equal to 0,5%, I computed :
> Turbulence kinetic energy = 0.00375 m²/s²
> Specific turbulence dissipation = 8.748177652797064 1/s


Thank you guys !