CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

sonicFoam "prism" tutorial not working when inlet velocity is increased

Register Blogs Members List Search Today's Posts Mark Forums Read

View Poll Results: what should be the first check to fix "negative initial temperature" error
0 folder 2 50.00%
constant folder 1 25.00%
system folder 1 25.00%
Voters: 4. You may not vote on this poll

Like Tree1Likes
  • 1 Post By hemanthgrylls

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   July 16, 2019, 05:08
Lightbulb sonicFoam "prism" tutorial not working when inlet velocity is increased
  #1
New Member
 
Dumbledore
Join Date: Jun 2019
Posts: 10
Rep Power: 6
hemanthgrylls is on a distinguished road
hello foamDudes

I took up a tutorial case $FOAM_TUTORIALS/compressible/sonicFoam/RAS/prism

the default inlet velocity boundary condition in the tutorial was 650(which is M =1.89)

In 0 folder i have changed velocity from (650 0 0) to (1030 0 0) at inlet, top and bottom wall boundaries to get to Mach 3 over the prism, i havent changed anything else and then this error shows up while i try to run it:

HTML Code:
Time = 3.25e-05

Courant Number mean: 0.056095 max: 0.45764
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.00361405, Final residual = 1.25973e-07, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0121551, Final residual = 5.68609e-07, No Iterations 2
smoothSolver:  Solving for e, Initial residual = 0.00566519, Final residual = 3.57166e-07, No Iterations 2


--> FOAM FATAL ERROR:
Negative initial temperature T0: -32.7705

    From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>]
    in file /home/pawan/OpenFOAM/OpenFOAM-v1812/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 54.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&)sh: 1: addr2line: not found
 addr2line failed
#1  Foam::error::abort()sh: 1: addr2line: not found
 addr2line failed
#2  Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::TEs(double, double, double) constsh: 1: addr2line: not found
 addr2line failed
#3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, bool)sh: 1: addr2line: not found
 addr2line failed
#4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct()sh: 1: addr2line: not found
 addr2line failed
#5  ?sh: 1: addr2line: not found
 addr2line failed
#6  __libc_start_mainsh: 1: addr2line: not found
 addr2line failed
#7  ?sh: 1: addr2line: not found
 addr2line failed
Aborted (core dumped)
please help me out in fixing this error to run the simulation at high mach numbers.
(i have later changed k and epsilon values in 0 folder but still the same story)

Last edited by hemanthgrylls; July 16, 2019 at 08:23. Reason: typos
hemanthgrylls is offline   Reply With Quote

 

Tags
sonicfoam, supersonic flows

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
serial udf to parallel udf radioss Fluent UDF and Scheme Programming 10 January 19, 2019 08:56
Defined inlet velocity is different with inlet velocity on CFX post jonpewpew CFX 2 November 2, 2017 16:40
Setting Density for Velocity Inlet Face arkie87 FLUENT 0 November 7, 2012 15:15
Velocity inlet boundary condition for porous medium Chander CFX 3 March 11, 2012 21:18
UDF paraboloid velocity inlet Ronak Shah FLUENT 0 June 4, 2003 09:44


All times are GMT -4. The time now is 04:36.