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

ABL in OF2.0.x

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 22, 2011, 12:23
Default ABL in OF2.0.x
  #1
Senior Member
 
Andrea Pasquali
Join Date: Sep 2009
Location: Germany
Posts: 142
Rep Power: 7
andrea.pasquali is on a distinguished road
Hi,
I'm testing the ABL in OF 2.0.X.
I build a flat "2D" mesh with grading in Z.
In the picture attached you can see the velocity probed at 5 m from the inlet.
I found a very strange jump between the log function and the free-stream.
Can anyone help me to understand it? (I attached also the case)
Have anyone found problem with this BC?
Or maybe is there a bug?

Tahnks

Andrea
Attached Images
File Type: jpg ABL_OF_20x.jpg (34.1 KB, 192 views)
Attached Files
File Type: gz ABLtest.gz (5.1 KB, 77 views)
__________________
Andrea Pasquali
andrea.pasquali is offline   Reply With Quote

Old   July 26, 2011, 05:19
Default
  #2
New Member
 
Jochem
Join Date: May 2010
Posts: 28
Rep Power: 7
Jochem is on a distinguished road
Hello Andrea,

I have also struggled with this problem for quite some time.

I can probably guess what the problem is, but I am not sure. I can't seem to open your case.

The problem I think it is, is that in OF-2.0 the code uses 2 different equations for ABL. If you look at the code for ABL you can see that it uses all the parameters dictated in the ABLcondition file except Ustar. It uses all these parameters to calculate this Ustar with an equation Ustar=U_ref*kappa/log((Z_ref+Z0_)/min(Z0,0.001) So it will take the minimum of Z0 or 0.001.

A little bit further down the code it calculates the ABL inlet wind speed with equation Ustar=U_ref*kappa/log((Z_ref+Z0_)/(Z0). So there will be a shift in your Ustar and this can be the cause for your problem.

This difference in equation is new in OF-2.0, in OF-1.7.x there was no difference. I don't know the reason yet. You can solve this in two ways, you can adjust the code of you can take into account the difference when you choose the parameters.

I hope this helps.

I also have some problems with the nutkroughwallfunction for the turbineSiting case where this ABL is used. If someone reads this tread, maybe they can join the discussion. I am trying to simulate a flow over a flat plane and try to have always the same wind profile using the ABL condition and the wallfunction.

Regards,

Jochem
Jochem is offline   Reply With Quote

Old   August 9, 2011, 11:29
Default
  #3
Member
 
inginer's Avatar
 
Ovidiu Michiu
Join Date: Apr 2010
Location: Ingolstadt, Germany
Posts: 52
Rep Power: 7
inginer is on a distinguished road
Hello,

I'm using ABL to calculate the flow over a terrain. Now I'm doing the mesh with snappyHexMesh options, but still don't have a good surface inside the box yet.

At this moment I have a stupid question. I read but still it isn't so clear for me. So, here is the question:
How can you calculate
ABLConditions: Ustar, Uref, Href, z0, zGround,turbulentKE; and
initialConditions: flowVelocity, pressure, turbulentKE, turbulentEpsilon?

Best Regards,
Ovidiu
inginer is offline   Reply With Quote

Old   August 9, 2011, 19:12
Default
  #4
Senior Member
 
Julien de Charentenay
Join Date: Jun 2009
Location: Australia
Posts: 229
Rep Power: 9
julien.decharentenay is on a distinguished road
Send a message via Skype™ to julien.decharentenay
Could you replace the min(z0, 0.001) by max(z0,0.001)? This approach is traditionally used to prevent divide by zero problem. I think that it could like be a mistake.

Also the parenthesis do not match on the expression: Ustar=U_ref*kappa/log((Z_ref+Z0_)/(Z0).
__________________
---
Julien de Charentenay

Last edited by julien.decharentenay; August 9, 2011 at 19:13. Reason: Add another comment
julien.decharentenay is offline   Reply With Quote

Old   August 22, 2011, 05:56
Default
  #5
New Member
 
Jochem
Join Date: May 2010
Posts: 28
Rep Power: 7
Jochem is on a distinguished road
Hi Ovidi,

zGround is the lowest part in your domain, this is the easy one. The Abl condition and initial conditions can be calculated from winddata, I suspect that is the right way.

I use the simplified loglaw to calculate a surface roughness with 2 datapoints (same point of time, different hight). With this data I calculate Ustar, and fill in the surface roughness with z0. H_ref is the place where the free stream velocity U_ref is reached. H_ref is in my case the maximum height of the domain, because I suspect that the freestream velocity isn't reached within the height of my case. I use the formula to calculate U_ref.

Flow velocity and pressure can remain the same is in the tutorial. For turbulentKE and turbulent epsilon I have found a document on the forum, I don't know where so I've attached the document. You can see the formulas to calculate the constants for the initial condition.

Julien,

The way you've suggest is also a good way to prevent this mistake. And you are correct, I've made a mistake with the parenthesis.

If anyone wants to make additions or corrections to my explanation, please do so.

Regards,

Jochem
Attached Files
File Type: pdf Input.pdf (57.7 KB, 184 views)

Last edited by Jochem; August 30, 2011 at 09:01.
Jochem is offline   Reply With Quote

Old   August 29, 2011, 11:23
Default
  #6
Member
 
inginer's Avatar
 
Ovidiu Michiu
Join Date: Apr 2010
Location: Ingolstadt, Germany
Posts: 52
Rep Power: 7
inginer is on a distinguished road
Hi Jochem,

Thank you for useful information.
Also, can you attache again the document which calculates TurbulentKE and turbulent epsilon. Thks.

Best Regards,
Ovidiu,
inginer is offline   Reply With Quote

Old   August 30, 2011, 04:48
Default
  #7
Senior Member
 
Andrea Pasquali
Join Date: Sep 2009
Location: Germany
Posts: 142
Rep Power: 7
andrea.pasquali is on a distinguished road
Dear All,
I solved it using the 17x version where the ABL bc is ok.
For the zGround I think is better to use the lowest z of the inlet patch (Z_low_inlet).
If you have a model with a global lowest Z (Z_low_global) below the lowest Z inlet (Z_low_inlet) and you choice the first Z as zGround = Z_low_global, you make a mistake in your inlet log profile.
I think another problem is when your Z_low_inlet changes in Y direction.
Did anyone found this problem? Is in OpenFOAM something to take in account this?
I think the fast way to overcome this, it's to separate the inlet patch in different patches @ different Z_low_inlet.

Regards
__________________
Andrea Pasquali
andrea.pasquali is offline   Reply With Quote

Old   August 30, 2011, 07:25
Post
  #8
Member
 
inginer's Avatar
 
Ovidiu Michiu
Join Date: Apr 2010
Location: Ingolstadt, Germany
Posts: 52
Rep Power: 7
inginer is on a distinguished road
Hi,

I'm studding the wind over a complex terrain. I toke like example the simpleWindFoam tutorial. I start the simulations and I'm receiving this error after a few iterations:

#0 Foam::error:rintStack(Foam::Ostream&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linuxGccDPOpt/libOpenFOAM.so"
#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/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#6 Foam::incompressible::RASModels::kEpsilon::correct () in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#7 main in "/home/openfoam/OpenFOAM/openfoam-1.7.0/applications/bin/linuxGccDPOpt/simpleFoam"
#8 __libc_start_main in "/lib/libc.so.6"
#9 _start at /usr/src/packages/BUILD/glibc-2.10.1/csu/../sysdeps/i386/elf/start.S:122
Floating point exception


If you receive something similar can you explain why I'm receiving this kind of error?
Thank you in advance,

Best Regards,
Ovidiu
inginer is offline   Reply With Quote

Old   August 30, 2011, 08:48
Default
  #9
Senior Member
 
Andrea Pasquali
Join Date: Sep 2009
Location: Germany
Posts: 142
Rep Power: 7
andrea.pasquali is on a distinguished road
If I don't remember bad, int the tutorial you have in the file k:

value uniform 0.0;

So you divide per 0 and got floating point.
Change it as:

value $internalField;



Regards
__________________
Andrea Pasquali
andrea.pasquali is offline   Reply With Quote

Old   August 31, 2011, 06:06
Default
  #10
Member
 
inginer's Avatar
 
Ovidiu Michiu
Join Date: Apr 2010
Location: Ingolstadt, Germany
Posts: 52
Rep Power: 7
inginer is on a distinguished road
Hi,

Thank you guys for the information. It's working.
__________________
Best Regards,
Ovi
inginer is offline   Reply With Quote

Old   October 4, 2011, 11:28
Default
  #11
Member
 
inginer's Avatar
 
Ovidiu Michiu
Join Date: Apr 2010
Location: Ingolstadt, Germany
Posts: 52
Rep Power: 7
inginer is on a distinguished road
Hello,

I'm running a simulation with a wind passing on a terrain. The solution converge and I want to plot the velocity profile. I manage to do this but in excel (see the picture) and i want to plot in ParaView. How can I do it?

Thank you
Attached Images
File Type: jpg Vel_profile.JPG (16.8 KB, 74 views)
__________________
Best Regards,
Ovi
inginer is offline   Reply With Quote

Old   October 5, 2011, 08:29
Default
  #12
New Member
 
Jochem
Join Date: May 2010
Posts: 28
Rep Power: 7
Jochem is on a distinguished road
Hello,

You have to use the filter->Data Analysis -> plot over line. Type in the coordinates you want to use.

Regards,

Jochem
Jochem is offline   Reply With Quote

Old   August 28, 2012, 17:04
Default
  #13
New Member
 
Rafael Valenzuela Musura
Join Date: Feb 2012
Posts: 27
Rep Power: 5
rafamusura is on a distinguished road
Hi to everyone, I'm simulating the flow throw a delta wing in open sky, I need to generate the ABL but it doesn't work! I'm using the OF2.1.1 maybe that's the problem.

Please Ovidiu can you shre your 0 and patch files for taking a look at them?

thanks!
rafamusura is offline   Reply With Quote

Old   December 12, 2012, 11:51
Default
  #14
New Member
 
Burak
Join Date: Nov 2012
Posts: 5
Rep Power: 4
Burak_1984 is on a distinguished road
Hi there I am trying to run a different terrain in turbineSiting tutorial.But I ran into problems with that so I went back to the tutorial case and tried to edit the wind direction.I think my mesh is fine since it passes the checkmesh.The strange thing is when I set the ABL conditions.

Uref 8.16;
Href 70;
z0 uniform 0.1;
turbulentKE 1.3;
windDirection (1 0 0);
zDirection (0 0 1);
zGround uniform 0.0;

The problem is simple like in the Turbinesiting tutorial.From the west there is an inlet from east there is an outlet the terrain is defined like a wall and top and sides are no slip

With these initial Condition the solution converges in 66 iterations.So when I change the wind direction to (1 1 0) the solution crashes.But if I use (5 2 0) the solution still converges.

Create time

Create mesh for time = 1

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
bounding k, min: 0 max: 3.42645150226 average: 1.30628839732
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.11;
Prt 1;
}

Creating field source list from sourcesProperties


SIMPLE: convergence criteria
field p tolerance 0.01
field U tolerance 0.001
field "(k|epsilon)" tolerance 0.001


Starting time loop

Time = 2

smoothSolver: Solving for Ux, Initial residual = 0.473678546445, Final residual = 0.0253675867531, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.545508812572, Final residual = 0.0410704803125, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.490732100575, Final residual = 0.0369494371199, No Iterations 3
GAMG: Solving for p, Initial residual = 0.235186667396, Final residual = 0.0115677539951, No Iterations 4
time step continuity errors : sum local = 0.000123942210288, global = -1.07094754161e-05, cumulative = -1.07094754161e-05
smoothSolver: Solving for epsilon, Initial residual = 0.0142418185171, Final residual = 0.00129445182949, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.443069945282, Final residual = 0.0234141579111, No Iterations 3
ExecutionTime = 1.13 s ClockTime = 1 s



Time = 3

smoothSolver: Solving for Ux, Initial residual = 0.46653783383, Final residual = 0.0263371603412, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.309793353435, Final residual = 0.0170232186791, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.205484283561, Final residual = 0.0129615975667, No Iterations 3
GAMG: Solving for p, Initial residual = 0.16854903988, Final residual = 0.0155268143013, No Iterations 3
time step continuity errors : sum local = 0.000210952976415, global = 1.76511520803e-05, cumulative = 6.94167666417e-06
smoothSolver: Solving for epsilon, Initial residual = 0.00969959500694, Final residual = 0.000656979811585, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.344065799334, Final residual = 0.021449332508, No Iterations 3
ExecutionTime = 1.78 s ClockTime = 1 s

Time = 4

smoothSolver: Solving for Ux, Initial residual = 0.149131904737, Final residual = 0.00828018673637, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.186625072917, Final residual = 0.0100424579985, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.139523403237, Final residual = 0.00937017636232, No Iterations 3
GAMG: Solving for p, Initial residual = 0.20638215344, Final residual = 0.0181687040399, No Iterations 5
time step continuity errors : sum local = 0.000109995029081, global = -1.4898546506e-05, cumulative = -7.95686984185e-06
smoothSolver: Solving for epsilon, Initial residual = 0.00758607525273, Final residual = 0.000583556877572, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.215939251023, Final residual = 0.0141821403433, No Iterations 3
bounding k, min: 0 max: 3.32130962682 average: 1.32598953792
ExecutionTime = 2.47 s ClockTime = 2 s

Time = 5

smoothSolver: Solving for Ux, Initial residual = 0.231364300243, Final residual = 0.0131708756433, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.103479618454, Final residual = 0.00572050956008, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.278417864979, Final residual = 0.0194834887617, No Iterations 3
GAMG: Solving for p, Initial residual = 0.269803139772, Final residual = 0.0259705324425, No Iterations 2
time step continuity errors : sum local = 0.000175945244605, global = 5.86352385936e-05, cumulative = 5.06783687518e-05
#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"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#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 "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#5 at kEpsilon.C:0
#6 Foam::incompressible::RASModels::kEpsilon::correct () in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#7
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/simpleFoam"
#8 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#9
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/simpleFoam"
Floating point exception (core dumped)


As you can see there is a bounding at k just before it crashes.The last iteration I persume says there is a problem at Kepsilon model.Floating point expection usually means division by zero I figured that out but simply all the parameters are same and the solution just crashes.

I am using OpenFoam version 2.1 and I am using exactly the same wall conditions and solvers as proposed by the tutorial.

for Epsilon : AtmBoundaryLayerInletEpsilon and EpsilonWallFunction

for k: kqrWallFunction

for nut: nutkAtmRoughWallFunction

is used.

I also checked the previous errors where you had to change min(Zo,0.001) to max(Zo,001) but that doesnt exist.I also tried to change the relaxation values;no use.

I won't post my solution and other sheets since they are exactly the same with the tutorial ,since I didn't change anything.The terrain is also AskervereinHill.

If somebody could give me some insight that would be terrific

My Regards
Burak
Burak_1984 is offline   Reply With Quote

Old   December 17, 2012, 03:31
Default
  #15
New Member
 
Burak
Join Date: Nov 2012
Posts: 5
Rep Power: 4
Burak_1984 is on a distinguished road
Well I did some more digging and operation about my problem;

#2 Uninterpreted:


#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"


#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 "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"


#5 at kEpsilon.C:0


#6 Foam::incompressible::RASModels::kEpsilon::correct () in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"


Floating point exception (core dumped)



The highlighted parts indeed meant that in the model K or epsilon value was getting the value of "zero" and division was causing a floating point.




For my case for the modified "turbine siting tutorial" I changed the values of k at case/0/k where the highlighted parts were changed from 0.0 to 1.3 which is the internal field value that $turbulentKE calls from "initial conditions" .But I guess something other than zero would work since considering what I read it gets revised and updated no matter what anyway.


At the end it looks like this;


#include "include/initialConditions"

dimensions [0 2 -2 0 0 0 0];

internalField uniform $turbulentKE;

boundaryField
{
#include "include/ABLConditions"

outlet
{
type inletOutlet;
inletValue uniform 1.3;
value $internalField;
}
inlet
{
type uniformFixedValue;
value $turbulentKE;
uniformValue constant $turbulentKE;
}
"terrain_.*"
{
type kqRWallFunction;
value uniform 1.3;
}

ground
{
type zeroGradient;
}

#include "include/sideAndTopPatches"
}



This worked for me and the solution converged.It even worked for a vertical flow like wind direction ( -5 -5 0) Northeast-->SouthWest wind which was causing an immidiate crash for me.

I hope somebody finds this usefull;I was getting annoyed on the posts where the guy explained his situation and when he got the answer he wrote like "Okay I got the solution thanks" and didn't explain the outcome .

Cheers
Burak
Burak_1984 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
Simulating ABL flow cm_jubayer OpenFOAM 0 June 10, 2011 10:52
ABL shear pmarek OpenFOAM Paraview & paraFoam 0 March 9, 2011 06:35
Pressure BC's for ABL modelling jgpeder OpenFOAM 1 June 9, 2010 14:33


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