CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User 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 Search this Thread 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: 17
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,912
Rep Power: 36
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 in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541)
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.
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: 17
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: 17
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: 434
Rep Power: 22
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: 17
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: 434
Rep Power: 22
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 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
energy in sonicFoam joern OpenFOAM Running, Solving & CFD 1 September 24, 2019 01:15
3-D Water Flow in Fluent to Evaluate Pressures ahammack FLUENT 2 September 1, 2009 23:41
Negative pressures in the flow domain mahut FLUENT 2 September 27, 2007 06: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 22:28


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