CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

simpleFoam Blows up with Parabolic Inlet Pipe Flow

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

LinkBack Thread Tools Search this Thread Display Modes
Old   December 6, 2010, 20:37
Default simpleFoam Blows up with Parabolic Inlet Pipe Flow
Stefano Wahono
Join Date: Aug 2010
Location: Melbourne, Australia
Posts: 42
Rep Power: 14
swahono is on a distinguished road
Hi All,

I am simulating a pipe with sudden expansion with a parabolic inlet velocity using simpleFoam.

The simulation started with laminar for 2000 iterations (residuals flat off after this) - solution attached.

Centreline Inlet Velocity is 0.178 m/s. Re_D=40,750.
GroovyBC was used to set the parabolic inlet velocity profile.
Then I restarted the simulation using RNGKEpsilon model.
The simulation crashed after 2 iterations! Consistently.
Error message indicate divide by zero problem for k and epsilon (See below).

Please NOTE that if I use NON-parabolic inlet (i.e. fixedValue), everything is fine - I got a good converged solution.

I used Upwind Div scheme with cellLimited interpolation scheme.
I have tried both PCG and GAMG for pressure, and PBiCG and GAMG for U|k|epsilon.
I have also tried varying URF.

I suspect the problem is in my BCs and bad starting solution for k and epsilon (since the flow is now parabolic from the inlet).
But, I can't figure out. I've also tried setting parabolic inlet profiles for k and epsilon to no avail.

Please Help anyone!.

the 0/epsilon:
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.05;
value uniform 1e-06;
the 0/k:
type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value uniform 0.0001;
Thank you very much.
I hope I do get any help this time.
Please.. It's quite difficult to get help from this forum.


Time = 2

DILUPBiCG:  Solving for Ux, Initial residual = 0.00855909, Final residual = 7.45812e-10, No Iterations 17
DILUPBiCG:  Solving for Uy, Initial residual = 0.0467024, Final residual = 6.71306e-10, No Iterations 17
DILUPBiCG:  Solving for Uz, Initial residual = 0.0490525, Final residual = 5.95911e-10, No Iterations 18
GAMG:  Solving for p, Initial residual = 0.192423, Final residual = 8.18271e-11, No Iterations 40
GAMG:  Solving for p, Initial residual = 0.0344647, Final residual = 6.74495e-11, No Iterations 35
GAMG:  Solving for p, Initial residual = 0.0101951, Final residual = 7.61577e-11, No Iterations 32
time step continuity errors : sum local = 8.98665e-13, global = 1.02105e-14, cumulative = 6.25679e-15
[4] #0  Foam::error::printStack(Foam::Ostream&) in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/"
[4] #1  Foam::sigFpe::sigFpeHandler(int) in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/"
[4] #2  __restore_rt at sigaction.c:0
[4] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/"
[4] #4  void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/"
[4] #5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/"
[4] #6  Foam::incompressible::RASModels::RNGkEpsilon::correct() in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/"
[4] #7  main in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/applications/bin/linux64GccDPOpt/simpleFoam"
[4] #8  __libc_start_main in "/lib64/"
[4] #9  Foam::regIOobject::writeObject(Foam::IOstream::streamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/applications/bin/linux64GccDPOpt/simpleFoam"
Attached Images
File Type: jpg n1200-laminar-Umag.jpg (46.0 KB, 24 views)
File Type: jpg ParabolaInletVelocity.jpg (19.3 KB, 16 views)
swahono is offline   Reply With Quote


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
Poiseuille flow (flow in a pipe) in SOLA-VOF. Parabolic distribution of velocities. Anton Main CFD Forum 9 April 29, 2009 11:31
Pulsating flow with non-reflective inlet boundary cc1000 CFX 6 April 27, 2009 09:10
Urgent: Modeling the drag on a stationary sphere in a parabolic Pres driven pipe flow farhan CFX 4 April 9, 2009 01:36
pulsating static pressure around zero at pipe flow David Kim FLUENT 0 May 16, 2006 14:47
flow in pipe with momentum and energy source Atit Koonsrisuk CFX 3 March 29, 2005 17:58

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