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

sonicFoam for supersonic simulation with butane gas

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 20, 2017, 10:24
Exclamation sonicFoam for supersonic simulation with butane gas
  #1
Member
 
power
Join Date: Jun 2014
Posts: 86
Rep Power: 12
enginpower is on a distinguished road
Hello all,
I'm trying to simulate butane gas flowing through the geometry in the following picture sonicFoam.jpg which reach supersonic flow in some throats connecting in tube and out tube.
All patches not indicated in the picture, are signed like walls.

I want to simulate real gas, and I calculated all parameters to get it using Peng-Robinson equation of state as you can check in the files: shared.zip

You can also find inside them all applied boundary conditions. Actually it is the first time I'm using sonicFoam solver, so I'm not so expert about that.
Running the simulation I got this error, but it's not clear to me where is the problem:
Code:
Maximum number of iterations exceeded functions are : 

15
(
abort
coded
patchProbes
probes
psiReactionThermoMoleFractions
removeRegisteredObject
residuals
rhoReactionThermoMoleFractions
setTimeStep
sets
surfaces
systemCall
timeActivatedFileUpdate
writeDictionary
writeObjects
)



    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::janafThermo<Foam::PengRobinsonGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::janafThermo<Foam::PengRobinsonGas<Foam::specie> >, Foam::sensibleInternalEnergy>]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::PengRobinsonGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::PengRobinsonGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#4  ? at ??:?
#5  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6  ? at ??:?
Aborted (core dumped)
All other similar errors in the forum are different from that because it lists some functions but I've not understood where to set them. And what exeeds 15 iterations? Not clear this point.

Please give me some useful ways to get a good simulation of that, mainly on the boundary condition settings

Thanks

Last edited by enginpower; April 20, 2017 at 10:28. Reason: grammar
enginpower is offline   Reply With Quote

Old   April 21, 2017, 08:37
Default
  #2
Member
 
power
Join Date: Jun 2014
Posts: 86
Rep Power: 12
enginpower is on a distinguished road
I solved some errors.
Now it runs but I think it doesn't go over the first iteration. Here the log file is:
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 4.1
Exec   : sonicFoam
Date   : Apr 21 2017
Time   : 14:19:44
Host   : "engine-VirtualBox"
PID    : 10986
Case   : /home/engine/OpenFoam/2DdobleTube
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       sutherland;
    thermo          janaf;
    equationOfState PengRobinsonGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
}

Creating field kinetic energy K

No MRF models present

No finite volume options present


Starting time loop

Time = 5e-07

Courant Number mean: 0.576416 max: 108.333
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.0892473, No Iterations 1001
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.31323, No Iterations 1001
DILUPBiCG:  Solving for e, Initial residual = 1, Final residual = 9.86494e-10, No Iterations 189
It's not clear to me if it performs the first iteration or not. But what does it mean
Code:
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.0892473, No Iterations 1001
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.31323, No Iterations 1001
DILUPBiCG:  Solving for e, Initial residual = 1, Final residual = 9.86494e-10, No Iterations 189
in the last lines? Does it runs 189 iterations?
What can I do to understand what does not work in the simulation? Finally I chose a kepsilon model as starting model, but the boundary conditions are almost as set in the attached files (see above).
Hope it can help who wants to suggest me something interesting to learn new things

Thanks anyway

Last edited by enginpower; April 21, 2017 at 08:40. Reason: grammars
enginpower is offline   Reply With Quote

Reply

Tags
butane, real gas model, sonicfoam, supersonic internal 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
Simulation of radial gas turbine, Inlet velocity components dmaz CFX 17 May 19, 2019 09:45
Gas simulation setup anon_p CONVERGE 7 April 30, 2019 11:49
Ideal gas v.s. Real gas for wing simulation Anna Tian Main CFD Forum 1 December 16, 2013 11:39
SonicFoam simulation probleme bruxellois OpenFOAM 7 October 30, 2010 11:06
Gas pressure question Dan Moskal Main CFD Forum 0 October 24, 2002 22:02


All times are GMT -4. The time now is 16:22.