RAMjet diffuser flow
 

Hi,
New to CFD trying to simulate flow in a ramjet diffuser with Openfoam
my geometry
Attachment 24734

Attachment 24735

Attachment 24736

Attachment 24737

Using rhoCentralFoam solver it keeps saying maxIter ireached limit specieThermo
Using sonicFoam now
simulation is running but rho=0 (isnt sonicFoam compressible solver)
also used setFields to set non uniform boundary conditions

defaultFieldValues
(
volScalarFieldValue p 85592
volScalarFieldValue epsilon 35008
volScalarFieldValue k 74
volVectorFieldValue U (703.2 0 0)
volScalarFieldValue T 279.06

);

regions
( cylinderAnnulusToCell
{
p1 ( 0.525985 0 0 ) ;
p2 ( 1.54485 0 0 );
outerRadius 0.140 ;
innerRadius 0.105 ;

fieldValues
(
volScalarFieldValue p 409540
volScalarFieldValue epsilon 35008
volScalarFieldValue k 74
volVectorFieldValue U (250 0 0)
volScalarFieldValue T 494

Plz advice

body
 

Attachment 24746
body geometry

used setFields to create a normal shock in th region of the cowl where it becomes straight as non uniform conditions to subsonic.Boundary Top is the semi cylindrical shell covering the region from inlet to cowl.

sonicFoam error
 

getting an error using sonicFoam solver
Time = 0.0237024

Courant Number mean: 0.0151219 max: 14.2721
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.000870388, Final residual = 1.50063e-06, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.00374136, Final residual = 1.09437e-09, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.00147596, Final residual = 1.95045e-06, No Iterations 1
DILUPBiCG: Solving for e, Initial residual = 0.000371632, Final residual = 3.91739e-06, No Iterations 1


--> 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/bir/OpenFOAM/OpenFOAM-2.1.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::specieThermo<Foam::eConstThermo<Foam::perfec tGas> >::T(double, double, double (Foam::specieThermo<Foam::eConstThermo<Foam::perfe ctGas> >::*)(double) const, double (Foam::specieThermo<Foam::eConstThermo<Foam::perfe ctGas> >::*)(double) const, double (Foam::specieThermo<Foam::eConstThermo<Foam::perfe ctGas> >::*)(double) const) const in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#3 Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::eConstThermo<Foam:: perfectGas> > > > >::calculate() in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#4 Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::eConstThermo<Foam:: perfectGas> > > > >::correct() in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#5
in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/sonicFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/sonicFoam"
Aborted (core dumped)

plz advice

Hi Devayan,

Is your problem solved.

--> 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/bir/OpenFOAM/OpenFOAM-2.1.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::specieThermo<Foam::eConstThermo<Foam::perfec tGas> >::T(double, double, double (Foam::specieThermo<Foam::eConstThermo<Foam::perfe ctGas> >::*)(double) const, double (Foam::specieThermo<Foam::eConstThermo<Foam::perfe ctGas> >::*)(double) const, double (Foam::specieThermo<Foam::eConstThermo<Foam::perfe ctGas> >::*)(double) const) const in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#3 Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::eConstThermo<Foam:: perfectGas> > > > >::calculate() in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#4 Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::eConstThermo<Foam:: perfectGas> > > > >::correct() in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#5
in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/sonicFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
in "/home/bir/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/sonicFoam"
Aborted (core dumped)

ya was able to sove it

Hi Devayan,

Many thankx for you reply,
can you offer me some suggestion how to eliminate this issue.

Looking forward to your reply.

Regards,
Rahul

I am a newbie at CFD. This issue I solved by using different boundary condition which is very problem specific. Please elaborate your issue maybe the problem lies there.