[Pyp]

This is the error I'm getting when I solve my code, would appreciate any help

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

[wyldckat]

... 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...
... something something... Google...
... something something... zero pressure?

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?

[Pyp]

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.

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

[wyldckat]

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:

1. You've incorrectly used the name "Wall" for the patch you've stated you wanted to monitor for forces, given that there is no patch with that name...
2. And that the example you've used for the forces function object is from an older version of OpenFOAM, hence the warnings about "outputControl"...
3. This:
   

   Mean and max Courant Numbers = 0.7855960745 171.2690737
   deltaT = 5.838734162e-10
   Time = 5.8387342e-10

   is indicative that the boundary conditions are utterly wrong, given that for a maximum Courant to be 171 and the time step be 1-e10, then you're dealing with an utterly massive implosion or explosion of fluid pressure being driven from cell A to cell B... worse case scenario, it may be on a scale similar to the speed of light, but I didn't do the math to double-check that...

I thought that the thread How to give enough info to get help was clearer on what is needed to diagnose a problem, but I guess it's not...


What I or anyone else need in order to diagnose the problem:

1. The files in the folder "0".
2. The file "constant/polyMesh/boundary", after the mesh is generated.
3. The files inside "constant", not the subfolders, because those describe how the models are set-up and one small mistake is enough to break everything.
4. The files inside "system".

Because as far as I can tell, the problem is that the pressure field is incorrectly initialized.




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:

writeControl      timeStep;

writeInterval 1.0;

And run the solver. It will write to disk all time steps and then you can see the results in ParaView to see how the fluid is behaving and how fast it's really going... also look at the pressure values...