|February 7, 2012, 06:24||
Error, rhosimplefoam solver
Join Date: Feb 2012
Posts: 3Rep Power: 6
At first, great work you're doing here.
I'm completely new to openfoam, but i try to do some simulation of a gas burner.
when i start the solver, the following error-message appears:
Starting time loop
Time = 1
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.0290357, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0162366, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.014761, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.517934, Final residual = 0.00675935, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 9.90403e-05, No Iterations 598
DICPCG: Solving for p, Initial residual = 0.401541, Final residual = 3.90591e-05, No Iterations 214
DICPCG: Solving for p, Initial residual = 0.0953977, Final residual = 9.29665e-06, No Iterations 266
time step continuity errors : sum local = 13.4007, global = 0.347056, cumulative = 0.347056
rho max/min : 0.821852 0.001
DILUPBiCG: Solving for omega, Initial residual = 1, Final residual = 0.0530614, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.0229082, No Iterations 2
ExecutionTime = 712.65 s ClockTime = 724 s
Time = 2
DILUPBiCG: Solving for Ux, Initial residual = 0.400034, Final residual = 0.0251938, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.457358, Final residual = 0.0190353, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.421311, Final residual = 0.0167393, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.0557149, No Iterations 2
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 Foam::hPsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::calculate() in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#4 Foam::hPsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::correct() in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#6 __libc_start_main in "/lib/libc.so.6"
Floating point exception (core dumped)
Would be great if anyone could help me.
As i already mentioned, i'm completely new to openfoam, solvers,...
Thanks in advance for your help!
|March 6, 2012, 03:30||
Join Date: Aug 2009
Posts: 24Rep Power: 9
If you're still having this problem, can you post the content of your /constant/thermophysicalProperties dictionary file?
|July 17, 2012, 13:37||
I had exactly the same problem. I am trying to simulate a flow through a simple catalytic converter and I am using rhoPorousMRFSimpleFoam and I get the same error during the second iteration. I decreased the relaxation factors for p and rho to 0.05 each and that seems to work. I have decent results but unfortunately my timestep continuity error is quite high. It oscillates between 60 and 25
|July 17, 2012, 17:04||
check your fvsolution equation relaxation factors, make them bigger than 0.9.
you should have p and rho for field and the rest for equation relaxation factors.
also checking the thermo is a good idea, post it here maybe.
|July 18, 2012, 22:53||
Thanks for the reply. I am sorry but i don't understand your suggestion.
Do you want me to assign pressure and density relaxation to be greater than 0.9? Is there a reason why this might work?
Also, I setup the case almost exactly like the tutorial case ( rhoporousMRFsimpleFOAM - angledductexplicit) and I really don't understand the reason for the problem. I tried removing the porous zone and checking how rhoSimpleFoam would work and I still have the same problem. On the other hand, I don't have any problems with the incompressible case. This surely points towards the boundary and the initial conditions of the case.
Here are my boundary conditions
Inlet - zeroGradient
Outlet - fixedValue
Inlet - flowRateInletVelocity
outlet - inletOutlet
walls - fixedValue ( 0 for no-slip)
Inlet - fixedValue
Outlet - fixedvalue
walls - zeroGradient
This case is supposed to be a heat transfer case and hence I have to specific temperatures on both the inlet and outlet . I am not really confident if this over does the boundary conditions.
I am really sorry but the files are in my office and I don't think I have the permissions to share them. I will try my best to properly represent my case though.
|July 19, 2012, 09:58||
I was struggling with the same problem and I got one case to work by setting those values in fvsolution as follows.
Be advised I am using rhoSimplecFoam .
However, I am now finding out the mesh quality might play a role as well. I will keep you posted on my findings.
|July 19, 2012, 10:55||
Thanks a lot for the reply.
I actually was able to solve my problem.
The problem with my setup was that, the porous resistance was quite large and I somehow overlooked that and used an explicit porosity formulation which the solver didn't like. I changed the fvSolution dictionary to use the implicit porosity formulation and BAZINGA !!
But, I will definitely remember your advice for future problems.
Note - I had to use a pressure under relaxation value of 0.2 and density under relaxation of 0.05. It might have worked for higher values, but I didn't want to lose my almost perfect results :-)
|Thread||Thread Starter||Forum||Replies||Last Post|
|thobois class engineTopoChangerMesh error||Peter_600||OpenFOAM||4||August 2, 2014 09:52|
|A New Solver for Supersonic Combustion||nakul||OpenFOAM||6||March 30, 2014 03:00|
|A New Solver for Supersonic Combustion||nakul||OpenFOAM Announcements from Other Sources||18||February 19, 2013 08:48|
|Development of a Low mach PISO solver||nishant_hull||OpenFOAM Programming & Development||0||August 25, 2009 12:48|
|why the solver reject it? Anyone with experience?||bearcat||CFX||6||April 28, 2008 14:08|