CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   RAMjet diffuser flow (http://www.cfd-online.com/Forums/openfoam/122584-ramjet-diffuser-flow.html)

devayabir August 22, 2013 11:22

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

devayabir August 22, 2013 12:00

body
 
1 Attachment(s)
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.

devayabir August 23, 2013 00:49

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

srahul November 19, 2013 12:33

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)

devayabir November 22, 2013 02:45

ya was able to sove it

srahul November 22, 2013 03:46

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

devayabir November 24, 2013 03:37

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.


All times are GMT -4. The time now is 06:52.