|
[Sponsors] |
sonicFoam for supersonic simulation with butane gas |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 20, 2017, 10:24 |
sonicFoam for supersonic simulation with butane gas
|
#1 |
Member
power
Join Date: Jun 2014
Posts: 86
Rep Power: 12 |
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) 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 |
|
April 21, 2017, 08:37 |
|
#2 |
Member
power
Join Date: Jun 2014
Posts: 86
Rep Power: 12 |
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 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 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 |
|
Tags |
butane, real gas model, sonicfoam, supersonic internal flows |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |