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

Parabolic inlet velocity profile

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

Like Tree5Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   December 14, 2006, 08:47
Default Hej, For my project I desid
  #21
mss
Guest
 
Posts: n/a
Hej,

For my project I desided to use the icoFoam case, but I need to
specify velocity inlet (it should be parabolic). Could u give me some
hint how I can do it?

Thank u,
Rita
  Reply With Quote

Old   December 14, 2006, 10:24
Default See Bernhard's responses above
  #22
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 698
Rep Power: 12
msrinath80 is on a distinguished road
See Bernhard's responses above. I used his utility without any problems for 2D cases. Haven't tried for 3D yet; but I guess there should be no problems.

Quoting from Bernhard:

Place the source in a directory. Create a directory called 'Make' in it in which you create a file named 'options' with this content

EXE_INC = \
-I$(LIB_SRC)/cfdTools/lnInclude \
-I$(LIB_SRC)/finiteVolume/lnInclude

EXE_LIBS = -lfiniteVolume

and a file called 'files' with this content

setParabolicInlet.C
EXE = $(FOAM_USER_APPBIN)/setParabolicInlet

After a wmake the executable is in your path and you can start using it.

BTW: you'll also need a creatFields.H:

Info<< "Vector field U\n" << endl;
volVectorField U
(
IOobject
(
"U",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
msrinath80 is offline   Reply With Quote

Old   December 14, 2006, 10:37
Default Thank u very much, I will try
  #23
mss
Guest
 
Posts: n/a
Thank u very much, I will try it.

Rita
  Reply With Quote

Old   January 17, 2007, 04:51
Default Hej, Sorry maybe it is a no
  #24
mss
Guest
 
Posts: n/a
Hej,

Sorry maybe it is a not good question about the parabolic inlet...But could u give me some hint where I should put exatly the setParabolicInlet.C Should I put it in the OpenFOAM/OpenFOAM-1.3/applications/solvers/incompressible/icoFoam or where? I am asking because as I understood that the way is very important in OpenFOAM.

Thank u,
Rita
  Reply With Quote

Old   January 17, 2007, 10:26
Default Not very critical to put it in
  #25
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 698
Rep Power: 12
msrinath80 is on a distinguished road
Not very critical to put it in a particular place. However, an organized arrangement always helps. I have mine put in a separate directory:

[username@hostname parabolic_velocity_inlet]$ pwd
/home/username/OpenFOAM/username-1.3/custom_utils/Bernhard_Gschaider/parabolic_v elocity_inlet
[username@hostname parabolic_velocity_inlet]$ ls -la
total 84K
drwxr-xr-x 3 username users 4.0K Nov 15 21:12 .
drwxr-xr-x 3 username users 4.0K Nov 14 16:58 ..
-rw-r--r-- 1 username users 156 Nov 14 16:58 createFields.H
drwxr-xr-x 3 username users 4.0K Nov 15 21:12 Make
-rw-r--r-- 1 username users 7.2K Nov 14 16:58 setParabolicInlet.C
-rw-r--r-- 1 username users 33K Nov 15 21:12 setParabolicInlet.dep
[username@hostname parabolic_velocity_inlet]$
msrinath80 is offline   Reply With Quote

Old   January 17, 2007, 18:00
Default I can't find the solution prop
  #26
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
I can't find the solution proposed by Hrvoje. Can someone post it here?

Regards,
Alberto
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   January 17, 2007, 18:07
Default Attached. http://www.cfd-o
  #27
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,761
Rep Power: 21
hjasak will become famous soon enough
Attached.

parabolicVelocity_HJ_17Jan2007.tgz

Enjoy,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   January 18, 2007, 09:55
Default Thanks! Alberto
  #28
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Thanks!

Alberto
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   January 23, 2007, 05:28
Default Hej pUl|, Thank you very
  #29
mss
Guest
 
Posts: n/a
Hej pUl|,

Thank you very much for advice. Now my ParabolicInlet is working
fine.


I created my own folder:
"OpenFOAM/margst-1.3/applications/solvers/Parabolic_velocity_inlet"
and compiled all files there.

The problem with this is that I can "see" my new case only within
OpenFOAM shell and ,for some reason, I can look at the results withy
Paraview.


For example, i create the class:

Case : IcoFoam,
Case root : OpenFOAM/margst-1.3/applications/solvers/Parabolic_velocity_inlet
Case name: test

But as soon as I close OpenFOAM, the case disappears for good from
the case brauser. How can one "save" the case? Plus, how can one use
ParaView for my newly created class?



Best regards,
Rita
  Reply With Quote

Old   January 24, 2007, 05:47
Default Hej I want to create a lin
  #30
mss
Guest
 
Posts: n/a
Hej

I want to create a linear profile of Temperature Boundary condition on
the one of the vertical walls. Before I was adding it
as a vector, but if I want to refine my mesh I have to change my vector
manualy.

Could you give me a hint what I should add to the file ParabolicInlet.C to
have both conditions (parabolic Inlet for velocity and linear profile for
temperature).

Thank you in advance,
Rita
  Reply With Quote

Old   April 17, 2007, 08:32
Default Hi Bernhard, stupid questio
  #31
New Member
 
Bummi
Join Date: Mar 2009
Posts: 3
Rep Power: 8
bummi is on a distinguished road
Hi Bernhard,

stupid question:

I would like to modify my boundary conditions by using your code. Everything works fine, but instead of a list with the data. Simply the maximum value of the velocity is responded in my boundary

So

Inlet
{
type fixedValue;
value uniform (0 0 0);
}
results in

Inlet
{
type fixedValue;
value uniform (1 0 0);
}

if my maximum velocity is 1. Any Idea????

Thanks Bummi
bummi is offline   Reply With Quote

Old   April 17, 2007, 17:27
Default Hi Bummi! You refering to m
  #32
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Hi Bummi!

You refering to my posting from April 26, 2006, are you? That code makes an assumption about the orientation of the patch: that it is parallel to the y-z-plane. Most propably your patch is parallel to xz or xy. Then the result should be the way you described it.
You'll have to modify the code to fit your geometry or rotate the geometry. Or write a more general utility.

regards
Bernhard
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   April 18, 2007, 02:34
Default Hi Bernhard, I'll try this!
  #33
New Member
 
Bummi
Join Date: Mar 2009
Posts: 3
Rep Power: 8
bummi is on a distinguished road
Hi Bernhard,

I'll try this! Thank you very much, especially for the fast response!!

regards

Bummi
bummi is offline   Reply With Quote

Old   April 18, 2007, 03:15
Default Hi Bernhard again, the dire
  #34
New Member
 
Bummi
Join Date: Mar 2009
Posts: 3
Rep Power: 8
bummi is on a distinguished road
Hi Bernhard again,

the direction was not the problem. It was correct, but I added some other parts to your code for modifying the inlet concentration and those caused the problem. By deleting them everything works.

Stupid me.

Regard

Bummi
bummi is offline   Reply With Quote

Old   April 20, 2007, 02:29
Default Hi Hvr, Could you re-post you
  #35
pnguyen1
Guest
 
Posts: n/a
Hi Hvr,
Could you re-post your "parabolicVelocity_HJ_17Jan2007.tgz", please ?
I could not download it !
Thank you very much!
Phuc-Danh
  Reply With Quote

Old   April 20, 2007, 02:34
Default The archive is good, I've just
  #36
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 647
Rep Power: 11
dmoroian is on a distinguished road
The archive is good, I've just tried it.
Be careful when you download it to change its name from something like parabolicVelocity_HJ_17Jan2007-3679.unk to parabolicVelocity_HJ_17Jan2007.tgz

Dragos
dmoroian is offline   Reply With Quote

Old   April 20, 2007, 02:50
Default Hi Hvr, Sorry, I can download
  #37
pnguyen1
Guest
 
Posts: n/a
Hi Hvr,
Sorry, I can download it now.
Thank you very much!
Phuc-Danh
  Reply With Quote

Old   April 20, 2007, 02:51
Default Hi Dragos, Thanks for your ad
  #38
pnguyen1
Guest
 
Posts: n/a
Hi Dragos,
Thanks for your additional precision !
Phuc-Danh
  Reply With Quote

Old   May 3, 2007, 07:54
Default Dragos, Once you download th
  #39
Member
 
Radu Mustata
Join Date: Mar 2009
Location: Zaragoza, Spain
Posts: 96
Rep Power: 8
r2d2 is on a distinguished road
Dragos,
Once you download the files, where do you have to
put them?
In:
~/OpenFOAM/OpenFOAM-1.3/src/finiteVolume/
fields/fvPatchFields/derivedFvPatchFields/

and

/home/radu/OpenFOAM/OpenFOAM-1.3/src/
finiteVolume/lnInclude/
???
What next? Do you have an example of U in one of the 0īs with the correct syntax?
Thanks,
Radu
r2d2 is offline   Reply With Quote

Old   August 14, 2007, 15:10
Default Hi to everyone, I'm trying
  #40
Member
 
victor
Join Date: Mar 2009
Location: mexico city, MX
Posts: 50
Rep Power: 8
torvic is on a distinguished road
Hi to everyone,

I'm trying to implement a parabolic profile in an inlet for reactingFoam, since I want to study an open flame and literature suggest to have such a developed profile before entering the crossflow plane. So, I'm trying two options: to compile the Bernhard's "SetParabolicInlet.C" and that of Dr. Jasak "ParabolicVelocity" under OF-1.4.1.

In the first one i get an error message concerning to fvPatchFieldFields.H. this is the message:

Making dependency list for source file parabolic.C
could not open file fvPatchFieldFields.H for source file parabolic.C
SOURCE=parabolic.C ; g++ -m32 -Dlinux -DDP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-40 -I/home/foam-1.4.1/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude -I/home/foam-1.4.1/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude -IlnInclude -I. -I/home/foam-1.4.1/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude -fPIC -pthread -c $SOURCE -o Make/linuxGccDPOpt/parabolic.o
parabolic.C:36:32: error: fvPatchFieldFields.H: No such file or directory
parabolic.C: In function 'int main(int, char**)':
parabolic.C:74: error: 'fvPatchVectorFieldField' was not declared in this scope
parabolic.C:74: error: 'Upatches' was not declared in this scope
make: *** [Make/linuxGccDPOpt/parabolic.o] Error 1

I looked for this file but didn't find it. I'm not sure if it is not anymore in OF-1.4.1, since I compiled it well under OF-1.3.

With respect to the one of Dr. Jasak, I modified the line he suggests (Alberto Passalacqua's post):

const vectorField& c = patch().Cf();

but I don't know if i have to compile as the one of bernhard's: create a Make, add file and options files.

I also tried to find the message Dr. Jasak refers to concerning a previous long message about compiling and linking new libraries (February 13 2007 answering to Thomas Jung) but i don't find it.

Please, I would thank so much if someone can give me a hint of how to proceed.

Thanks for reading,

best

V
torvic 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
2D air parabolic velocity profile ilker FLUENT 2 November 12, 2008 09:43
UDF Unsteady velocity parabolic profile Rashad FLUENT 2 February 29, 2008 03:36
parabolic velocity profile? bssdyl FLUENT 4 March 22, 2006 12:32
problem in 3d parabolic velocity profile Lokesh FLUENT 8 August 11, 2005 05:36
Parabolic temperature Inlet Profile in a tube majestywzh FLUENT 0 April 9, 2003 06:37


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