CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   what this error comes from? (http://www.cfd-online.com/Forums/openfoam-solving/110152-what-error-comes.html)

immortality December 4, 2012 18:17

what this error comes from?
 
Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package ePsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::calculate() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#4 Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::ePsiThermo(Foam::fvMesh const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#5 Foam::basicPsiThermo::addfvMeshConstructorToTable< Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > > >::New(Foam::fvMesh const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#6 Foam::basicPsiThermo::New(Foam::fvMesh const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#7
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
Floating point exception

immortality December 5, 2012 03:08

I use a 64-bit version on a system by 2Gb Ram.could the problem arises from this point?or what?

Hanzo December 5, 2012 04:24

I would check the source code for

Code:

ePsiThermo<Foam::pureMixture<Foam::constTran  sport<Foam::specieThermo<Foam::hConstThermo<Foam::  perfectGas> > > > >::calculate()
There could be the chance that you somehow divide by a zero value like division by zero density or similar. I once had that case.

immortality December 5, 2012 06:40

may be this is because of initial condition.I put pressure and temperature zero due to empty channel that should become full.is it incorrect?

immortality December 5, 2012 08:26

for instance see this time step:


Mean and max Courant Numbers = 0.140537 0.296461
deltaT = 1.31579e-07
Time = 4.86842e-06

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
smoothSolver: Solving for Ux, Initial residual = 1.13168e-05, Final residual = 2.04904e-17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.00894171, Final residual = 1.97135e-17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 1.43136e-05, Final residual = 1.55745e-10, No Iterations 201
smoothSolver: Solving for omega, Initial residual = 0.000587127, Final residual = 1.92363e-11, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.00730572, Final residual = 1.34671e-10, No Iterations 2
ExecutionTime = 38.88 s ClockTime = 58 s

Hanzo December 5, 2012 23:07

Quote:

Originally Posted by immortality (Post 395886)
may be this is because of initial condition.I put pressure and temperature zero due to empty channel that should become full.is it incorrect?

Maybe you should post an overview what you are doing. Setting Temperature to zero may cause your error. I am not sure if there is a valid setting for such an initial condition.
Try to check the
Code:

calculate()
method if your error still occurs.

immortality December 22, 2012 11:57

thanks.how shuld I use calculate method?

Hanzo December 24, 2012 01:53

Did you try to find the source code of

Quote:

Originally Posted by Hanzo (Post 395858)

Code:

ePsiThermo<Foam::pureMixture<Foam::constTran  sport<Foam::specieThermo<Foam::hConstThermo<Foam::  perfectGas> > > > >::calculate()



All times are GMT -4. The time now is 16:26.