|September 11, 2012, 06:41||
Supersonic turbulent free jet LES
Join Date: Aug 2012
Posts: 35Rep Power: 4
I am doing an LES simulation of a turbulent supersonic free jet using the rhoPimpleFoam solver.
My solution domain is a converging diverging nozzle with the free-stream extending 30 nozzle exit diameters. and 10 diameters high. There is also a section of free-stream above the nozzle to capture the flow behavior at the nozzle lip. The computational domain is only a 1/8th section using cyclic side boundaries.
I have included pictures below.
velocity contour at exit.jpeg
velocity contour at crash.jpeg
This mesh is very coarse for LES but at this stage I am just trying to establish my BC's and get a stable run.
I am using the flowRateInletVelocity BC so the simulation is driven by a mass flow not a pressure difference (using zeroGradient pressure BC at inlet).
I am trying to avoid a 'numerical shock' so I am currently ramping this inlet BC but it is way too gradual, the table is:
Even at this very slight ramp the simulation is still crashing after only 0.003sec and at this time it should be barely moving yet there is a max velocity of 117m/s! So obviously this is completely unstable.
I also fear that my outlet is poorly defined. At this stage I have my inlet boundary, the nozzle wall, then all other boundaries are effectively an 'outlet', which I have named freestream. Is it OK to have more than one outlet or do I need to define an actual outlet (perhaps the rightmost face) and then the freestream separately using some far-field condition?
At the moment I am using a totalPressure condition for all the free-stream boundaries but will use waveTransmissive when I perform the actual run to dissipate the shocks.
Basically I keep getting this error message when I run the case (rho is diverging, enthalpy equation is crashing):
--> FOAM FATAL ERROR:
Maximum number of iterations exceeded
From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const
in file /home/opencfd/OpenFOAM/OpenFOAM-2.1.0/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69.
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Foam::specieThermo<Foam::hConstThermo<Foam:erfec tGas> >::T(double, double, double (Foam::specieThermo<Foam::hConstThermo<Foam:erfe ctGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam:erfe ctGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam:erfe ctGas> >::*)(double) const) const in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so"
#3 Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::calculate() in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so"
#4 Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::correct() in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so"
#6 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
I have done the following with no success:
- Running it in PISO mode with nOuterCorrectors = 3
- using backward ddt scheme
- have increased all rel tol by 10 times (kept rel tol final = 0)
- have decreased the relaxation factors to 0.7 (plan to increase it to 1 when flow is developed)
Does anyone have any ideas on how to increase the stability of my run? Should I just use a uniform velocity inlet condition?
A lot of info I know but any help would be really appreciated.
|January 20, 2013, 11:05||
Join Date: Nov 2012
Posts: 168Rep Power: 4
I had the same problem now. I am using rhoPimpleFoam (LES) to simulate a free turbulent jet. The jet exit velocity is 70 and the coflow velocity is 10. I used the standard rhoPimpleFoam solver and did not make any modification on it. However, I had the same error messege from the solver. I guess it is caused by the wrong calculations of enthalpy, rho. How did you sovle that problem half a year ago?
I really appreciate it if you can give me some suggestions.
|March 13, 2013, 11:03||
Join Date: Mar 2013
Posts: 12Rep Power: 4
I have the same problem. I try to simulate a supersonic LES flow with openFoam and this error about the thermoLibrary pop up every time.
|March 18, 2013, 02:17||
Join Date: Mar 2009
Posts: 59Rep Power: 9
A few suggestions, firstly, you might want to ensure that the 'transonic' switch is set to 'on' in fvSolutions (in the PISO subdict, not sure if that is mentioned above).....failing that, and secondly, modify the initial conditions (using
setFields) to allow the central column of gas to 'move'.
The thermo error you observe is merely a reflection of the fact that h (enthalpy) has overshot owing to a poor/unstable estimate of the convection
term in the hEqn. Once you stabilize U&phi the rest should follow.
I had a similar geometry some years back but after some consideration
was able to replace the downstream section with a pressure boundary condition (p fixed or waveTransmissive, I forget which) at the exit to the divergent section of the nozzle.
|July 1, 2013, 08:11||
Join Date: Feb 2013
Location: IIT Kanpur
Posts: 12Rep Power: 4
I am trying to solve supersonic flow over backstep using rhoCentralFoam and sonicFoambut I m faced with this strange problem :
The error is:
--> FOAM FATAL ERROR : Maximum number of iterations exceeded
From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const
in file /home/liu/OpenFOAM/OpenFOAM-1.2/src/thermophysicalModels/specie/lnInclude/specie ThermoI.H at line 83.
Looking forward to constructive suggestions.
|boundary condition, flowrateinletvelocity, inlet boundary condition|
|Thread||Thread Starter||Forum||Replies||Last Post|
|LES of jet injection giving pulsations||nsalpeter||Main CFD Forum||3||October 2, 2009 19:12|
|Help!simulation of free water jet||Annie||FLUENT||3||June 29, 2009 00:51|
|turbulent jet||ramo||Main CFD Forum||1||September 4, 2005 07:43|
|LES of jet in cross flow||dirk||FLUENT||0||April 25, 2005 08:39|
|Radial jet with free surface using FIDAP||Dr. Ahmed Hassaneen||Main CFD Forum||0||October 28, 1999 13:20|