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

atmBoundaryLayerInletVelocity profiles

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 14, 2010, 09:51
Default atmBoundaryLayerInletVelocity profiles
  #1
Ivo
New Member
 
Ivo
Join Date: Jul 2010
Posts: 5
Rep Power: 8
Ivo is on a distinguished road
Hi,

I recently started using OpenFoam for modeling atmospheric processes. Since I am new to the program and OpenFoam version 1.7 includes a tutorial case, turbineSiting, which uses a simulated atmospheric boundary layer I decided to start by stuying this tutorial. The first problem I encountered was the fact that the Allrun script references the non-existing solver simpleWindFoam. However, since the case seems to run just fine using simpleFoam I used this instead.

Having run the case successfully I started to play with the input parameters and found out that the atmBoundaryLayerInletVelocity boundary condition displays some very strange behaviour. Defining the ABL by setting the velocity Uref to 15 m/s at a height Href of 300 meters and not changing any of the other parameters I get a velocity profile as shown in the attachment. The profile clearly doesn't reach a velocity of 15 m/s at the specified height of 300 meters, but around the 300 meter mark it does quickly increase to that value, after which it stays constant.

The formulas in the boundary condition source code seem valid enough to me, so perhaps there is a problem with one of the parameters which get passed to/from the routine? Or maybe simpleFoam is an improper solver for this problem after all?

Thank you very much in advance!
Attached Images
File Type: jpg inletVelocityProfile.jpg (13.7 KB, 292 views)
Ivo is offline   Reply With Quote

Old   July 15, 2010, 03:47
Default
  #2
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 14
henry is on a distinguished road
the simpleWindFoam demonstration solver is in the

OpenFOAM-1.7.?/tutorials/incompressible/simpleWindFoam/simpleWindFoam

directory. In the

OpenFOAM-1.7.?/tutorials/incompressible/simpleWindFoam

directory type

foamRunTutorials

which will compile simpleWindFoam and run it.

H
henry is offline   Reply With Quote

Old   July 15, 2010, 12:28
Default
  #3
Ivo
New Member
 
Ivo
Join Date: Jul 2010
Posts: 5
Rep Power: 8
Ivo is on a distinguished road
Hi Henry,

Thank you very much for your reply. I'm sorry I didn't know about the way to compile the simpleWindFoam solver. It tried your suggestion and I can now indeed compile and use it. However, the resulting velocity profile I get is unchanged from my original post, so it seems the problem is not with the solver. Do you perhaps have another suggestion what I could try?

Best wishes,

Ivo
Ivo is offline   Reply With Quote

Old   July 16, 2010, 05:35
Default atmBoundaryLayerInletVelocityFvPatchVectorField
  #4
New Member
 
Sergio Ferraris
Join Date: Jul 2010
Posts: 3
Rep Power: 8
SergioF is on a distinguished road
There was an error in the expression for the inlet velocity in the file
atmBoundaryLayerInletVelocityFvPatchVectorField.C

It was updated in our git repository for 1.7.x

Thanks

S
SergioF is offline   Reply With Quote

Old   July 16, 2010, 10:32
Default
  #5
Ivo
New Member
 
Ivo
Join Date: Jul 2010
Posts: 5
Rep Power: 8
Ivo is on a distinguished road
It now works fine, many thanks!!!

Ivo
Ivo is offline   Reply With Quote

Old   July 30, 2010, 11:05
Default
  #6
New Member
 
siri
Join Date: Apr 2010
Location: Norway
Posts: 16
Rep Power: 8
snippsnuske is on a distinguished road
Dear
I’m also working with the turbineSiting tutorial. It runs fine and I will now change the terrain file to a terrain file of my interest. I manage to use snappyHexMesh and a checkMesh gives me good result. But when I try to run simpleWndFoam I continuing to get this error massage. The only thing I change was the terrain. It looks like it reding the field p fine but it stops when rading U. I’m not used at reading different error massages and I would be very grateful if anyone can tell me what this message is indicating;
#13 __libc_start_main in "/lib/libc.so.6"#14 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/user/OpenFOAM/OpenFOAM-1.7.0/applications/bin/linuxGccDPOpt/simpleWindFoam"
snippsnuske is offline   Reply With Quote

Old   September 27, 2012, 04:54
Default
  #7
Member
 
Gregoire Junqua
Join Date: Jun 2011
Location: China
Posts: 58
Blog Entries: 1
Rep Power: 7
gregjunqua is on a distinguished road
Hi,
I am actualy trying to compile atmBoundaryLayerInletVelocity
the things is i get that error :


m::incompressible::atmBoundaryLayerInletVelocityFv PatchVectorField const&, Foam::fvPatch const&, Foam:imensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::fvPatchFieldMapper const&)'
Make/linuxGccDPOpt/simpleFoam.o:simpleFoam.C.text+0x3090): first defined here
Make/linuxGccDPOpt/atmBoundaryLayerInletVelocityFvPatchVectorField.o: (.bss+0x4): multiple definition of `Foam::incompressible::atmBoundaryLayerInletVeloci tyFvPatchVectorField::debug'
Make/linuxGccDPOpt/simpleFoam.o.bss+0x4): first defined here
Make/linuxGccDPOpt/atmBoundaryLayerInletVelocityFvPatchVectorField.o: (.bss+0x2): multiple definition of `Foam::incompressible::addatmBoundaryLayerInletVel ocityFvPatchVectorFieldpatchConstructorTofvPatchVe ctorFieldTable_'
Make/linuxGccDPOpt/simpleFoam.o.bss+0x2): first defined here
Make/linuxGccDPOpt/atmBoundaryLayerInletVelocityFvPatchVectorField.o: (.bss+0x1): multiple definition of `Foam::incompressible::addatmBoundaryLayerInletVel ocityFvPatchVectorFieldpatchMapperConstructorTofvP atchVectorFieldTable_'
Make/linuxGccDPOpt/simpleFoam.o.bss+0x1): first defined here
Make/linuxGccDPOpt/atmBoundaryLayerInletVelocityFvPatchVectorField.o: (.bss+0x0): multiple definition of `Foam::incompressible::addatmBoundaryLayerInletVel ocityFvPatchVectorFielddictionaryConstructorTofvPa tchVectorFieldTable_'
Make/linuxGccDPOpt/simpleFoam.o.bss+0x0): first defined here
collect2: ld a retourné 1 code d'état d'exécution
make: *** [/home/gregoire/OpenFOAM/OpenFOAM-2.1.1/platforms/linuxGccDPOpt/bin/bounSimpleFoam] Erreur 1

for come to here i had do :


• Add .atmBoundaryLayerInletVelocityFvPatchVectorField.C to the second line of Make/files, and
modify the final line to EXE = $(FOAM_USER_APPBIN)/
• Add #include "parabolicVelocityFvPatchVectorField.H" in the header of your
simpleFoam.C file, so that the solver knows the new boundary condition.
• wclean
rm -r Make/linux*
wmake

somebody can help me ?
or at least give me a good website for compile atmBoundaryLayerInlet
gregjunqua is offline   Reply With Quote

Reply

Tags
abl, boundary condition, turbinesiting

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
[ICEM] Reading/Writing Profiles in tilted geometry qascapri ANSYS Meshing & Geometry 0 August 14, 2009 11:04
Laminar wing profiles and drag gregorv OpenFOAM Running, Solving & CFD 4 December 4, 2007 14:25
Gas flow profiles in mixing tank srinivas FLUENT 0 November 8, 2006 05:52
Velocity profiles Rex FLUENT 4 July 23, 2002 19:35
component flow direction fieldnames for boundary profiles required. Ricky Wong FLUENT 1 May 12, 2000 10:36


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