|
[Sponsors] |
Can't simulate pressures over 80 bar in sonicFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
October 29, 2009, 08:24 |
Can't simulate pressures over 80 bar in sonicFoam
|
#1 |
Member
Dennis Rown
Join Date: Oct 2009
Posts: 51
Rep Power: 17 |
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. |
|
October 30, 2009, 00:18 |
|
#2 | |
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36 |
Quote:
__________________
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. |
||
October 30, 2009, 05:06 |
|
#3 |
Member
Dennis Rown
Join Date: Oct 2009
Posts: 51
Rep Power: 17 |
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 |
|
October 30, 2009, 09:33 |
|
#4 |
Member
Dennis Rown
Join Date: Oct 2009
Posts: 51
Rep Power: 17 |
I found out that with the solver "rhopSonicFoam" the problem does not occur!
|
|
October 30, 2009, 10:07 |
|
#5 |
Senior Member
Chris Sideroff
Join Date: Mar 2009
Location: Ottawa, ON, CAN
Posts: 434
Rep Power: 22 |
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?
|
|
October 30, 2009, 10:51 |
|
#6 |
Member
Dennis Rown
Join Date: Oct 2009
Posts: 51
Rep Power: 17 |
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; } |
|
October 30, 2009, 11:21 |
|
#7 |
Senior Member
Chris Sideroff
Join Date: Mar 2009
Location: Ottawa, ON, CAN
Posts: 434
Rep Power: 22 |
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.
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |