Solver for transonic flow?

Abhinay Kulkarni · February 22, 2012, 11:29

Hi Martin,

I am a beginer in openfoam as well as cfd. I wanted to somulate a stedy state compressible flow through a bent tube using the rhoSimpleFoam solver. My problem is the solver crashes after a few time time steps with the following error

Time = 186
smoothSolver: Solving for Ux, Initial residual = 0.00992166, Final residual = 0.000578669, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.00774996, Final residual = 0.000440937, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.0168676, Final residual = 0.000958473, No Iterations 4
DILUPBiCG: Solving for h, Initial residual = 0.0136817, Final residual = 8.63913e-05, No Iterations 2
#0 Foam::error:

rintStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::hPsiThermo<Foam:

ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:

erfectGas> > > > >::calculate() in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so"
#4 Foam::hPsiThermo<Foam:

ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:

erfectGas> > > > >::correct() in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so"
#5
in "/opt/openfoam201/platforms/linuxGccDPOpt/bin/rhoSimpleFoam"
#6 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#7
in "/opt/openfoam201/platforms/linuxGccDPOpt/bin/rhoSimpleFoam"
Gleitkomma-Ausnahme

Could you kindly help me out with understanding this error??

I dont know why but i am not able to upload my case directories on the forum. If you want to have alook at them i can send them to you via email.

Would be great to get some help from you.

Regards
Abhinay

MtnRunBeachBum · April 27, 2015, 11:02

I see that I am 6 years late to this discussion, however, thought I would add my 2 cents, coming from the experimental (wind tunnel) side of the house.

To the best of my knowledge, still, few CFD codes can reproduce transonic results well. The aero is just funny, to put it simply, in this regime - particularly right around 1.0M (say 0.93-1.05 in general). That being said, I would believe CFD, including OF, has come a long way in covering this regime since 2009, particularly with the density-based solvers.

Just some food for thought. Does anyone know of reports showing correlation of OF results (with the more recent versions) against older OF code or even wind tunnel data?

levvis · December 16, 2015, 04:59

Hi all

I'm also new to OpenFOAM, and obviously late to this discussion.

But anyone visiting this thread late, may want to look at another, more recent thread.

http://www.cfd-online.com/Forums/ope...ntralfoam.html

These guys have progressed a solver that can deal with transonic flow, and supersonic. Yet to fully try it, but I will be.

February 22, 2012, 11:29		#21
Abhinay Kulkarni New Member Join Date: Jan 2012 Posts: 23 Rep Power: 14	Hi Martin, I am a beginer in openfoam as well as cfd. I wanted to somulate a stedy state compressible flow through a bent tube using the rhoSimpleFoam solver. My problem is the solver crashes after a few time time steps with the following error Time = 186 smoothSolver: Solving for Ux, Initial residual = 0.00992166, Final residual = 0.000578669, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.00774996, Final residual = 0.000440937, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.0168676, Final residual = 0.000958473, No Iterations 4 DILUPBiCG: Solving for h, Initial residual = 0.0136817, Final residual = 8.63913e-05, No Iterations 2 #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #2 Uninterpreted: #3 Foam::hPsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::calculate() in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so" #4 Foam::hPsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::correct() in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so" #5 in "/opt/openfoam201/platforms/linuxGccDPOpt/bin/rhoSimpleFoam" #6 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6" #7 in "/opt/openfoam201/platforms/linuxGccDPOpt/bin/rhoSimpleFoam" Gleitkomma-Ausnahme Could you kindly help me out with understanding this error?? I dont know why but i am not able to upload my case directories on the forum. If you want to have alook at them i can send them to you via email. Would be great to get some help from you. Regards Abhinay

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Solver for an incompressible, turbulent flow with heat transfer	tH3f0rC3	OpenFOAM Running, Solving & CFD	9	June 17, 2019 06:12
How do I select solver options for external flow over an aircraft by fluent?	hadieliasi	FLUENT	5	May 2, 2011 03:54
pre-conditioning for low mach number compressible flow solver	Shenren_CN	Main CFD Forum	0	April 29, 2011 21:07
Inviscid Drag at subsonic, subcritical Mach #	Axel Rohde	Main CFD Forum	1	November 19, 2001 12:19
fluid flow fundas	ram	Main CFD Forum	5	June 17, 2000 21:31

April 27, 2015, 11:02		#22
MtnRunBeachBum New Member Jimmy Join Date: Apr 2015 Posts: 20 Rep Power: 11	I see that I am 6 years late to this discussion, however, thought I would add my 2 cents, coming from the experimental (wind tunnel) side of the house. To the best of my knowledge, still, few CFD codes can reproduce transonic results well. The aero is just funny, to put it simply, in this regime - particularly right around 1.0M (say 0.93-1.05 in general). That being said, I would believe CFD, including OF, has come a long way in covering this regime since 2009, particularly with the density-based solvers. Just some food for thought. Does anyone know of reports showing correlation of OF results (with the more recent versions) against older OF code or even wind tunnel data?

December 16, 2015, 04:59		#23
levvis New Member Lewis Join Date: Dec 2015 Posts: 14 Rep Power: 10	Hi all I'm also new to OpenFOAM, and obviously late to this discussion. But anyone visiting this thread late, may want to look at another, more recent thread. http://www.cfd-online.com/Forums/ope...ntralfoam.html These guys have progressed a solver that can deal with transonic flow, and supersonic. Yet to fully try it, but I will be.