CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Can't simulate pressures over 80 bar in sonicFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 29, 2009, 08:24
Post Can't simulate pressures over 80 bar in sonicFoam
  #1
Member
 
Dennis Rown
Join Date: Oct 2009
Posts: 51
Rep Power: 7
shangzung is on a distinguished road
Hello!

(sorry, I meant 900 bar, threshold value is : 89984650Pa)

Somehow with the sonicFoam I can not set the pressure in the setFields over 80.000.000 Pa without the simulation crashing. Is this something of a border?

Did someone ever simulated pressures around 1kbar?

Last edited by shangzung; October 29, 2009 at 08:46.
shangzung is offline   Reply With Quote

Old   October 30, 2009, 00:18
Default
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by shangzung View Post
Hello!

(sorry, I meant 900 bar, threshold value is : 89984650Pa)

Somehow with the sonicFoam I can not set the pressure in the setFields over 80.000.000 Pa without the simulation crashing. Is this something of a border?

Did someone ever simulated pressures around 1kbar?
Error message?
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   October 30, 2009, 05:06
Default
  #3
Member
 
Dennis Rown
Join Date: Oct 2009
Posts: 51
Rep Power: 7
shangzung is on a distinguished road
It does not simulate even one timestep. I have to add that I modified the sonicFoam.C for the adjustable timestep to work.


I get this error:

Creating turbulence model

Selecting turbulence model type laminar

Starting time loop

Time = 1e-12

Courant Number mean: 0 max: 0
deltaT = 1.2e-12
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.50669e-16, No Iterations 1
DILUPBiCG: Solving for e, Initial residual = 6.66672e-16, Final residual = 6.66672e-16, No Iterations 0
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/simulation/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/simulation/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::LimitedScheme<double, Foam::limitedLinearLimiter<Foam::NVDTVD>, Foam::limitFuncs::magSqr>::limiter(Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/simulation/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam::limitedSurfaceInterpolationScheme<double>::w eights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/simulation/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::fv::gaussConvectionScheme<double>::fvmDiv(Fo am::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) const in "/home/simulation/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#6 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::word const&) in "/home/simulation/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/sonicFoam"
#7 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/home/simulation/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/sonicFoam"
#8 main in "/home/simulation/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/sonicFoam"
#9 __libc_start_main in "/lib64/libc.so.6"
#10 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Gleitkomma-Ausnahme
shangzung is offline   Reply With Quote

Old   October 30, 2009, 09:33
Default
  #4
Member
 
Dennis Rown
Join Date: Oct 2009
Posts: 51
Rep Power: 7
shangzung is on a distinguished road
I found out that with the solver "rhopSonicFoam" the problem does not occur!
shangzung is offline   Reply With Quote

Old   October 30, 2009, 10:07
Default
  #5
Senior Member
 
Chris Sideroff
Join Date: Mar 2009
Location: Ottawa, ON, CAN
Posts: 369
Rep Power: 12
cnsidero is on a distinguished road
Quote:
Originally Posted by shangzung View Post
It does not simulate even one timestep. I have to add that I modified the sonicFoam.C for the adjustable timestep to work.
Not an answer to your question but I am interested in getting adjustable timestep working in my solver. Would you be willing to share your adjustable timestep modifications to sonicFoam?
cnsidero is offline   Reply With Quote

Old   October 30, 2009, 10:51
Default
  #6
Member
 
Dennis Rown
Join Date: Oct 2009
Posts: 51
Rep Power: 7
shangzung is on a distinguished road
Yeah, here's the code: (I use OF 1.6)

#include "fvCFD.H"
#include "basicPsiThermo.H"
#include "turbulenceModel.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{
#include "setRootCase.H"
#include "createTime.H"
#include "createMesh.H"
#include "createFields.H"
#include "initContinuityErrs.H"
#include "readTimeControls.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Info<< "\nStarting time loop\n" << endl;

while (runTime.loop())
{
Info<< "Time = " << runTime.timeName() << nl << endl;

#include "readPISOControls.H"
#include "CourantNo.H"
#include "setDeltaT.H"
#include "rhoEqn.H"

#include "UEqn.H"

#include "eEqn.H"


// --- PISO loop

for (int corr=0; corr<nCorr; corr++)
{
#include "pEqn.H"
}

turbulence->correct();

rho = thermo.rho();

runTime.write();

Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
<< nl << endl;
}

Info<< "End\n" << endl;

return 0;
}
shangzung is offline   Reply With Quote

Old   October 30, 2009, 11:21
Default
  #7
Senior Member
 
Chris Sideroff
Join Date: Mar 2009
Location: Ottawa, ON, CAN
Posts: 369
Rep Power: 12
cnsidero is on a distinguished road
Oops sorry, I mistook "adjustable timestep" to be "local timstep". I am trying to implement a timestep that varies throughout the field dependent on the local Co.
cnsidero is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
3-D Water Flow in Fluent to Evaluate Pressures ahammack FLUENT 2 September 1, 2009 22:41
energy in sonicFoam joern OpenFOAM Running, Solving & CFD 0 June 16, 2009 15:53
Negative pressures in the flow domain mahut FLUENT 2 September 27, 2007 05:07
Guage vs. Absolute pressures in Fluent. Riaan FLUENT 1 January 21, 2005 15:22
What is best way to simulate big particle yueroo FLUENT 2 April 13, 2001 21:28


All times are GMT -4. The time now is 17:48.