CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

floating point problem

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 3, 2019, 14:41
Post floating point problem
  #1
New Member
 
bren
Join Date: Apr 2019
Posts: 3
Rep Power: 7
flughafen is on a distinguished road
hello,
trying to do a simulation of transonic flow around a NACA 0012, these are the iterations working:

/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 6-e29811f5dff8
Exec : rhoSimpleFoam
Date : Apr 03 2019
Time : 19:33:42
Host : "fares-G500"
PID : 26089
I/O : uncollated
Case : /home/fares/OpenFOAM/fares-6/run/NACA0012_maillage_influence/50div
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


SIMPLE: Convergence criteria found
p: tolerance 0.0001
U: tolerance 0.0001
e: tolerance 0.0001
"(k|epsilon|omega)": tolerance 0.0001

Reading thermophysical properties

Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport sutherland;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
alphaEps 0.76923;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

No MRF models present

No finite volume options present

Starting time loop

fieldAverage fieldAverage1:
Starting averaging at time 0

forces forceCoeffs1:
Not including porosity effects
forceCoeffs forceCoeffs1:
Not including porosity effects
forces forces1:
Not including porosity effects
forces forceCoeffsAveraged1:
Not including porosity effects
forceCoeffs forceCoeffsAveraged1:
Not including porosity effects
forces forcesAveraged1:
Not including porosity effects
Time = 0.01

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.013686, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0148273, No Iterations 2
smoothSolver: Solving for e, Initial residual = 1, Final residual = 0.0120187, No Iterations 1
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0424933, No Iterations 9
time step continuity errors : sum local = 0.000205158, global = 8.41997e-06, cumulative = 8.41997e-06
smoothSolver: Solving for epsilon, Initial residual = 0.0105842, Final residual = 0.000108216, No Iterations 2
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0125213, No Iterations 2
ExecutionTime = 0.24 s ClockTime = 0 s

Reading/initialising field UMean
Reading/initialising field pMean
fieldAverage fieldAverage1 write:
Calculating averages

Writing average fields


forceCoeffs forceCoeffs1 write:
Cm = -5.1649e-20
Cd = 1.4965
Cl = 1.84411e-19
Cl(f) = 4.05567e-20
Cl(r) = 1.43855e-19

forces forces1 write:
sum of forces:
pressure : (40710 -108.964 5.01465e-15)
viscous : (47.5446 -0.00542527 7.86062e-18)
porous : (0 0 0)
sum of moments:
pressure : (-2.5818e-16 -1.40204e-15 55.203)
viscous : (-6.10257e-19 -4.63503e-18 0.00131471)
porous : (0 0 0)

forceCoeffs forceCoeffsAveraged1 write:
Cm = -5.1649e-20
Cd = 1.4965
Cl = 1.84411e-19
Cl(f) = 4.05567e-20
Cl(r) = 1.43855e-19

forces forcesAveraged1 write:
sum of forces:
pressure : (40710 -108.964 5.01465e-15)
viscous : (47.5446 -0.00542527 7.86062e-18)
porous : (0 0 0)
sum of moments:
pressure : (-2.5818e-16 -1.40204e-15 55.203)
viscous : (-6.10257e-19 -4.63503e-18 0.00131471)
porous : (0 0 0)

Time = 0.02

smoothSolver: Solving for Ux, Initial residual = 0.12776, Final residual = 0.00150867, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.363692, Final residual = 0.00411299, No Iterations 2
smoothSolver: Solving for e, Initial residual = 0.982097, Final residual = 0.00815914, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0312078, Final residual = 0.0012332, No Iterations 9
time step continuity errors : sum local = 0.000414326, global = 3.86681e-05, cumulative = 4.7088e-05
smoothSolver: Solving for epsilon, Initial residual = 0.0198776, Final residual = 0.000314401, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.278262, Final residual = 0.00435449, No Iterations 2
ExecutionTime = 0.31 s ClockTime = 0 s

fieldAverage fieldAverage1 write:
Calculating averages

Writing average fields


forceCoeffs forceCoeffs1 write:
Cm = -6.84639e-20
Cd = 2.07804
Cl = 2.37947e-19
Cl(f) = 5.05097e-20
Cl(r) = 1.87438e-19

forces forces1 write:
sum of forces:
pressure : (56521.4 -287.91 6.47175e-15)
viscous : (74.7537 -0.0236307 8.82175e-18)
porous : (0 0 0)
sum of moments:
pressure : (-3.33543e-16 -1.86143e-15 212.926)
viscous : (-6.6705e-19 -3.20608e-18 0.0177599)
porous : (0 0 0)

forceCoeffs forceCoeffsAveraged1 write:
Cm = -6.84639e-20
Cd = 2.07804
Cl = 2.37947e-19
Cl(f) = 5.05097e-20
Cl(r) = 1.87438e-19

forces forcesAveraged1 write:
sum of forces:
pressure : (56521.4 -287.91 6.47175e-15)
viscous : (74.7537 -0.0236307 8.82175e-18)
porous : (0 0 0)
sum of moments:
pressure : (-3.33543e-16 -1.86143e-15 212.926)
viscous : (-6.6705e-19 -3.20608e-18 0.0177599)
porous : (0 0 0)

Time = 0.03

smoothSolver: Solving for Ux, Initial residual = 0.112815, Final residual = 0.00143998, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.217492, Final residual = 0.00273693, No Iterations 2
smoothSolver: Solving for e, Initial residual = 0.544869, Final residual = 0.00365432, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0736136, Final residual = 0.0023271, No Iterations 6
time step continuity errors : sum local = 0.00072083, global = 8.21999e-05, cumulative = 0.000129288
smoothSolver: Solving for epsilon, Initial residual = 0.0599782, Final residual = 0.000729332, No Iterations 2
bounding epsilon, min: -26303.5 max: 414482 average: 2327.19
smoothSolver: Solving for k, Initial residual = 0.2211, Final residual = 0.00296538, No Iterations 2
bounding k, min: -232.538 max: 1175.73 average: 30.3552
ExecutionTime = 0.38 s ClockTime = 0 s

fieldAverage fieldAverage1 write:
Calculating averages

Writing average fields


forceCoeffs forceCoeffs1 write:
Cm = -7.45591e-20
Cd = 2.07067
Cl = 2.77413e-19
Cl(f) = 6.41474e-20
Cl(r) = 2.13266e-19

forces forces1 write:
sum of forces:
pressure : (56257.8 -206.193 7.54104e-15)
viscous : (137.617 -0.0885102 1.43959e-17)
porous : (0 0 0)
sum of moments:
pressure : (-3.80284e-16 -2.02561e-15 -37.2542)
viscous : (-1.03929e-18 -5.02648e-18 -0.00218632)
porous : (0 0 0)

forceCoeffs forceCoeffsAveraged1 write:
Cm = -7.45591e-20
Cd = 2.07067
Cl = 2.77413e-19
Cl(f) = 6.41474e-20
Cl(r) = 2.13266e-19

forces forcesAveraged1 write:
sum of forces:
pressure : (56257.8 -206.193 7.54104e-15)
viscous : (137.617 -0.0885102 1.43959e-17)
porous : (0 0 0)
sum of moments:
pressure : (-3.80284e-16 -2.02561e-15 -37.2542)
viscous : (-1.03929e-18 -5.02648e-18 -0.00218632)
porous : (0 0 0)

Time = 0.04

smoothSolver: Solving for Ux, Initial residual = 0.0440658, Final residual = 0.000444621, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.32332, Final residual = 0.00277778, No Iterations 2
smoothSolver: Solving for e, Initial residual = 0.0217963, Final residual = 0.000248265, No Iterations 1


/////////////////////////////////////////////////////////////////// and here is the error message :
#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::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"


i think that it is a pressure conditions related problem , please help me solve it
Attached Files
File Type: txt log.txt (8.9 KB, 1 views)
flughafen is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
OpenFoam problem with floating point exception Gabriela_ro OpenFOAM Running, Solving & CFD 0 February 26, 2018 04:12
MPI Error - simpleFoam - Floating Point Exception scott OpenFOAM Running, Solving & CFD 3 April 13, 2012 16:34
Ercoftac centrifugal pump case study bug, Openfoam 1.6-ext (Floating point exception) ArianeJasmin OpenFOAM Programming & Development 1 April 4, 2012 14:11
Problem with UDF compiling for kTkLW model Wantami FLUENT 0 July 18, 2011 05:11
having a problem whit floating point error and it's solution m.beh FLUENT 2 November 9, 2010 10:28


All times are GMT -4. The time now is 20:33.