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

Problem using atmBoundaryLayerInletVelocity

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 28, 2013, 10:21
Default Problem using atmBoundaryLayerInletVelocity
  #1
New Member
 
zako
Join Date: Apr 2012
Posts: 7
Rep Power: 5
bigzakk is on a distinguished road
Hi all,

I am trying to use atmBoundaryLayerInletVelocity for a very simple application however I am receiving this error i cant figure out. Thanks for any helps however small

Create time

Create mesh for time = 0

Reading field p

Reading field U

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 log in "/lib/x86_64-linux-gnu/libm.so.6"
#4 Foam::incompressible::atmBoundaryLayerInletVelocit yFvPatchVectorField::atmBoundaryLayerInletVelocity FvPatchVectorField(Foam::fvPatch const&, Foam:imensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#5 Foam::fvPatchField<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::incompres sible::atmBoundaryLayerInletVelocityFvPatchVectorF ield>::New(Foam::fvPatch const&, Foam:imensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#6 Foam::fvPatchField<Foam::Vector<double> >::New(Foam::fvPatch const&, Foam:imensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pisoFoam"
#7 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::GeometricB oundaryField(Foam::fvBoundaryMesh const&, Foam:imensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pisoFoam"
#8 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readField(Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pisoFoam"
#9 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readField(Foam::Istream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pisoFoam"
#10 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pisoFoam"
#11
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pisoFoam"
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pisoFoam"
Floating point exception (core dumped)
bigzakk is offline   Reply With Quote

Old   March 3, 2013, 17:43
Default
  #2
Senior Member
 
chegdan's Avatar
 
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 535
Rep Power: 18
chegdan will become famous soon enough
Zako,

"Floating point exception" is code for divide by zero or some other undefined numerical error. If you have something like a z0 = 0 and Href = 0 then there is going to be a log(0) in there for that BC. I tend to use 1e-20 instead of zero for small numbers that may induce this kind of error. Have a look at the BC itself in

Code:
$FOAM_SRC/turbulenceModels/incompressible/RAS/derivedFvPatchFields/atmBoundaryLayerInletVelocity
to track down exactly where this error is. Without seeing how you have defined the BC...I can't give you more info than this. Good Luck.
__________________
Dan

Find me on twitter @dancombest and LinkedIn
chegdan is offline   Reply With Quote

Old   March 4, 2013, 07:40
Default
  #3
New Member
 
zako
Join Date: Apr 2012
Posts: 7
Rep Power: 5
bigzakk is on a distinguished road
Thanks a million its works
bigzakk is offline   Reply With Quote

Old   March 9, 2013, 07:40
Default Problem with atmBoundaryLayerInletVelocity
  #4
New Member
 
Korichi Abdelkader
Join Date: Jan 2013
Posts: 6
Rep Power: 4
AbdelkaderDZ is on a distinguished road
Hi anyone,
I have error at the forth iteration when using atmBoundaryLayerInletVelocity for a simple case I use simpleFoam solver. Anyone can help me

ime = 4

DILUPBiCG: Solving for Ux, Initial residual = 0.0747753, Final residual = 0.00352462, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0842028, Final residual = 0.00496445, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.00994555, Final residual = 9.80074e-05, No Iterations 225
time step continuity errors : sum local = 0.0137186, global = 0.000340224, cumulative = -0.000110142
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
AbdelkaderDZ is offline   Reply With Quote

Old   March 9, 2013, 14:46
Default
  #5
Senior Member
 
Lieven
Join Date: Dec 2011
Location: Mol, Belgium
Posts: 295
Rep Power: 13
Lieven will become famous soon enough
Welcome to the forum AbdelkaderDZ.

For the looks of it, there is something wrong with your continuity. The cumulative error is simply way too big. If you are using a standard solver, this is often caused by incompatible boundary conditions. Could you post an overview of them of all fields (pressure, velocity, ...)?

Cheers,

L
Lieven is offline   Reply With Quote

Old   March 10, 2013, 03:10
Default Re:Re:Problem with atmBoundaryLayerInletVelocity
  #6
New Member
 
Korichi Abdelkader
Join Date: Jan 2013
Posts: 6
Rep Power: 4
AbdelkaderDZ is on a distinguished road
Thank you Lieven
I have used the standard simplefoam with the following boundary condition for u and p:
u boundary condition:

boundaryField
{
#include "include/ABLConditions"

outlet
{
type zeroGradient;
}

inlet
{
type atmBoundaryLayerInletVelocity;
Uref $Uref;
Href $Href;
n $windDirection;
z $zDirection;
z0 $z0;
value $internalField;
zGround $zGround;
}
upperWall
{
type symmetryPlane;
}
lowerWall
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
{
type empty;
}
p boundary condition
boundaryField
{
inlet
{
type zeroGradient;
}

outlet
{
type fixedValue;
value uniform 0;
}

upperWall
{
type symmetryPlane;
}

lowerWall
{
type zeroGradient;
}

frontAndBack
{
type empty;
}
}
AbdelkaderDZ is offline   Reply With Quote

Old   March 10, 2013, 05:31
Default
  #7
Senior Member
 
Lieven
Join Date: Dec 2011
Location: Mol, Belgium
Posts: 295
Rep Power: 13
Lieven will become famous soon enough
Nothing seems spectacular wrong with this so here are a few things you can do/check:
Does 'checkMesh' return 'mesh ok' (maybe post the final summary)?
Do you get the same problem when running with smaller time step?
One thing you can do is run your application using pimpleFoam instead of simpleFoam, which runs often more stable for larger time steps. Can you post the output of this here?
Which turbulence model are you using?

Chees,

L
Lieven 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
UDF compiling problem Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 04:43
Gambit - meshing over airfoil wrapping (?) problem JFDC FLUENT 1 July 11, 2011 05:59
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 06:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 19:13
Is this problem well posed? Thomas P. Abraham Main CFD Forum 5 September 8, 1999 14:52


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