|
[Sponsors] |
April 3, 2019, 15:41 |
floating point problem
|
#1 |
New Member
bren
Join Date: Apr 2019
Posts: 3
Rep Power: 7 |
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 |
|
|
|
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 05:12 |
MPI Error - simpleFoam - Floating Point Exception | scott | OpenFOAM Running, Solving & CFD | 3 | April 13, 2012 17:34 |
Ercoftac centrifugal pump case study bug, Openfoam 1.6-ext (Floating point exception) | ArianeJasmin | OpenFOAM Programming & Development | 1 | April 4, 2012 15:11 |
Problem with UDF compiling for kTkLW model | Wantami | FLUENT | 0 | July 18, 2011 06:11 |
having a problem whit floating point error and it's solution | m.beh | FLUENT | 2 | November 9, 2010 11:28 |