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

Power Law Velocity

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 10, 2013, 02:42
Default Power Law Velocity
  #1
Senior Member
 
starter
Join Date: Sep 2012
Posts: 109
Rep Power: 6
sihaqqi is on a distinguished road
Hi everyone

I am a beginner in openFoam. I have a geometry where two pipes are connected at 90degree. Flow is turbulent. I have made my geometry but now I want to have one seventh power law velocity profile set up in my boundary conditions. I have attached my files if anyone experienced can spare a few moments. Pipe diameter is 72.54mm. Mach number is 0.22 and medium is air. Fow enters vertically upwards in +Y and then turn in the bend at +X. Patches are inlet, outlet, upstream, downstream, symmetry.

Following is input of my "epsilon" file

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

internalField uniform 1;

boundaryField
{
inlet
{
type powerLawVelocity;
n (7 0 0);
y (0 0.7254 0);
maxValue 10;
value uniform (0 0 0);
}

outlet
{
type inletOutlet;
inletValue $internalField;
}

walls
{
type epsilonWallFunction;
value uniform 0;
}

symmetry
{
type inletOutlet;
}
}


Input for P

internalField uniform 0;

boundaryField
{
inlet
{
type powerLawVelocity;
n (7 0 0);
y (0 0.7254 0);
maxValue 10;
value uniform (0 0 0);
}

outlet
{
type fixedValue;
value uniform 0;
}

upstream
{
type wall;
}

downstream
{
type wall;
}

symmetry
{
type zeroGradient;
}

}
sihaqqi is offline   Reply With Quote

Old   March 10, 2013, 03:42
Default
  #2
Member
 
Michiel
Join Date: Oct 2010
Location: Delft, Netherlands
Posts: 97
Rep Power: 6
michielm is on a distinguished road
could you specify what you want us to answer/what the problem is you are experiencing?
michielm is offline   Reply With Quote

Old   March 10, 2013, 09:56
Default
  #3
Senior Member
 
starter
Join Date: Sep 2012
Posts: 109
Rep Power: 6
sihaqqi is on a distinguished road
Dear Michiel

Many thanks for your response. My problem is I want some piping tutorials with power law boundary conditions and the reason behind them. This is what mu Masters by research thesis is about. I want to know how to set up 1/7th power law velocity profile in all k, p, epsilon, u and Nut and NuTilda files so that my system works well. Patches are inlet, outlet, symmetry, upstream, downstream. Do you have any tutorial with some explanation on this type of problem with correct values of n and y. With my values, simpleFoam is giving errors which I do not understand at all.

Regards
sihaqqi is offline   Reply With Quote

Old   March 10, 2013, 10:53
Default
  #4
Member
 
Michiel
Join Date: Oct 2010
Location: Delft, Netherlands
Posts: 97
Rep Power: 6
michielm is on a distinguished road
I can't find a tutorial in the OpenFOAM tutorials folder (i tried: grep -r -i "powerlawvelocity" *), but you can find an example at the chalmers PhD course on CFD in OF:
Tutorial Power law velocity inlet, InterFoam

EDIT: just an addition, I noticed that the report of the chalmers course uses 1.5-dev so I guess the powerLawVelocity BC is also only available in the -dev versions of OpenFOAM.
michielm is offline   Reply With Quote

Reply

Tags
1/7th power law, velocity profile

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
Difficulty in calculating angular velocity of Savonius turbine simulation alfaruk CFX 8 December 3, 2013 17:51
plzzz help::energy equation using power law method sepidehkavousi Main CFD Forum 5 April 20, 2012 08:05
Using Power Law Model in Fluent kurne FLUENT 0 October 17, 2011 08:59
power law flow through porous media Francesca Chiusa FLUENT 1 November 3, 2006 15:30
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 09:11


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