|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Anon
Join Date: Apr 2019
Posts: 4
Rep Power: 7 ![]() |
This is the error I'm getting when I solve my code, would appreciate any help
Code:
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 5.60524582e-05, Final residual = 5.976319661e-11, No Iterations 1000 smoothSolver: Solving for Uz, Initial residual = 3.327855114e-06, Final residual = 2.489972098e-15, No Iterations 2 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib64/libc.so.6" #3 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::eConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, bool) at ??:? #4 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::eConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:? #5 ? at ??:? #6 __libc_start_main in "/lib64/libc.so.6" #7 ? at ??:? /tmp/slurmd/job15209510/slurm_script: line 13: 12077 Floating point exceptionrhoCentralFoam |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 ![]() ![]() ![]() ![]() ![]() ![]() |
![]() ... Since you didn't follow the instructions given here: How to give enough info to get help ... I guess I will have to get my crystal ball out of the closet... Hold on... now where is it?.... rummage... rummage.... Oh, there it is... How did it get so dirty with dust? ... squicky sounds while cleaning... OK, now we're in business! Oh magic ball, magic ball... what could this error be?Oh, it's becoming clearer... ... something something... division by zero...Looks like you might have not correctly defined the boundary conditions on the pressure fields!? You're dealing with compressible flow, which means that a complete pressure and temperature formulation is needed, so that the fluid properties can be properly calculated back and forth. In other words: pressure defined with 0 or temperature defined with 0 will literally mean vacuum and absolute 0 K temperature, respectively. Not exactly what this solver can work with. If that wasn't the problem/solution, well... what does a crystal ball know about CFD anyway?
__________________
|
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Anon
Join Date: Apr 2019
Posts: 4
Rep Power: 7 ![]() |
I hope this gives more information. Sorry for the previous post as it does not have the most information. Hopefully this has more information as I checked the pressure temperature and velocity files, the parameters only have a zerogradient.
I just want to say we're pretty brand new with OpenFOAM, so we're struggling with concepts. Any help would be appreciated. Thanks. Code:
Reading thermophysical properties Selecting thermodynamics package hePsiThermo<pureMixture<sutherland<eConst<perfectGas<specie>>,sensibleInternalEnergy>>> Reading field U Creating turbulence model Selecting turbulence model type laminar Selecting laminar stress model Stokes fluxScheme: Kurganov Starting time loop --> FOAM Warning : From function void Foam::timeControl::read(const Foam::dictionary&) in file db/functionObjects/timeControl/timeControl.C at line 103 Reading "/fast/users/a1685400/projectTest2/system/controlDict.functions.forces" Using deprecated 'outputControl' Please use 'writeControl' with 'writeInterval' forces forces: --> FOAM Warning : From function Foam::labelHashSet Foam::polyBoundaryMesh::patchSet(const Foam::UList<Foam::wordRe>&, bool, bool) const in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 793 Cannot find any patch or group names matching Wall Not including porosity effects forceCoeffs forces: --> FOAM Warning : From function Foam::labelHashSet Foam::polyBoundaryMesh::patchSet(const Foam::UList<Foam::wordRe>&, bool, bool) const in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 793 Cannot find any patch or group names matching Wall Not including porosity effects --> FOAM Warning : From function void Foam::timeControl::read(const Foam::dictionary&) in file db/functionObjects/timeControl/timeControl.C at line 103 Reading "/fast/users/a1685400/projectTest2/system/controlDict.functions.probes" Using deprecated 'outputControl' Please use 'writeControl' with 'writeInterval' Mean and max Courant Numbers = 0.7855960745 171.2690737 deltaT = 5.838734162e-10 Time = 5.8387342e-10 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 0.999999951, Final residual = 2.409799906e-08, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 0.9999999999, Final residual = 2.437318018e-14, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 4.545764333e-08, Final residual = 1.88119813e-16, No Iterations 2 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 1.981686762e-07, Final residual = 5.972131121e-11, No Iterations 1000 ExecutionTime = 520.72 s ClockTime = 521 s forceCoeffs forces execute: Coefficients Cm : 0 (pressure: 0 viscous: 0) Cd : 0 (pressure: 0 viscous: 0) Cl : 0 (pressure: 0 viscous: 0) Cl(f) : 0 Cl(r) : 0 Mean and max Courant Numbers = 0.0004586886529 0.09984120229 deltaT = 5.848021723e-10 Time = 1.1686756e-09 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 1.446727025e-05, Final residual = 1.24969054e-09, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 2.307628967e-05, Final residual = 5.89913015e-13, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 1.901951443e-10, Final residual = 3.945358971e-18, No Iterations 2 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 1.031491159e-08, Final residual = 1.194967732e-11, No Iterations 1000 ExecutionTime = 935.25 s ClockTime = 935 s forceCoeffs forces execute: Coefficients Cm : 0 (pressure: 0 viscous: 0) Cd : 0 (pressure: 0 viscous: 0) Cl : 0 (pressure: 0 viscous: 0) Cl(f) : 0 Cl(r) : 0 |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 ![]() ![]() ![]() ![]() ![]() ![]() |
Quick answer: The output you've provided does not give any hints as to what exactly is the cause of the original report.
The only things I can deduce from the output you've provided is that:
What I or anyone else need in order to diagnose the problem:
If you want to diagnose the issue on your own, it's really simple: change in controlDict these lines to the values I've written: Code:
writeControl timeStep; writeInterval 1.0; |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
How to apply a turbulence model to the rhoCentralFoam solver? | arussell92 | OpenFOAM Pre-Processing | 18 | July 25, 2022 04:26 |
Modify rhoCentralFoam: other equations of state | fivos | OpenFOAM Programming & Development | 5 | July 29, 2020 13:17 |
Inviscid simulation in rhoCentralFoam | chencw | OpenFOAM Running, Solving & CFD | 0 | April 13, 2018 15:49 |
SigFpe when running ANY application in parallel | Pj. | OpenFOAM Running, Solving & CFD | 3 | April 23, 2015 14:53 |
rhoCentralFoam solver with Slip BCs fails in Parallel Only | JLight | OpenFOAM Running, Solving & CFD | 2 | October 11, 2012 21:08 |