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   March 13, 2009, 02:58
Default Hi Jose, it works perfectly
  #121
New Member
 
Johannes Alken
Join Date: Mar 2009
Posts: 5
Rep Power: 8
alki is on a distinguished road
Hi Jose,

it works perfectly. Thank you.


regards

Johannes
alki is offline   Reply With Quote

Old   June 7, 2009, 23:08
Default
  #122
New Member
 
vinu varghese
Join Date: Apr 2009
Posts: 9
Rep Power: 8
vinu is on a distinguished road
Hi

I am using OpenFOAM 1.4.1. I did as Bernhard Gschaider said. But got the following error message. Does anyone know what went wrong?.

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

Old   June 16, 2009, 09:13
Default
  #123
New Member
 
Brugiere Olivier
Join Date: Mar 2009
Posts: 27
Rep Power: 8
brugiere_olivier is on a distinguished road
Hi all,

I want to create an intersection between a patch and a plane. I've find this part of code and I put some lines but that it doesn't work.
// Create a patch set intesecting a plane
point pnt(4.5,2,2.25);
vector spanwiseDirection(0,0,1);
plane pl1(pnt,spanwiseDirection);
cuttingPlane cutPlane(mesh.boundaryMesh().findPatchID("mur_bas" ),pl1);
const labelList& cutCells = cutPlane.cells();
word setName("someCells");
cellSet currentSet(mesh, setName, cutCells);

the error's message is :

wallGradU.C: In function âint main(int, char**)â:
wallGradU.C:121: error: no matching function for call to âFoam::cuttingPlane::cuttingPlane(Foam::label, Foam:lane&)â
/home/users/brugiere/OpenFOAM/OpenFOAM-1.4.1/src/sampling/lnInclude/cuttingPlane.H:142: note: candidates are: Foam::cuttingPlane::cuttingPlane(const Foam:rimitiveMesh&, const Foam:lane&, const Foam::labelList&)
/home/users/brugiere/OpenFOAM/OpenFOAM-1.4.1/src/sampling/lnInclude/cuttingPlane.H:133: note: Foam::cuttingPlane::cuttingPlane(const Foam:rimitiveMesh&, const Foam:lane&)
/home/users/brugiere/OpenFOAM/OpenFOAM-1.4.1/src/sampling/lnInclude/cuttingPlane.H:112: note: Foam::cuttingPlane::cuttingPlane(const Foam:lane&)
/home/users/brugiere/OpenFOAM/OpenFOAM-1.4.1/src/sampling/lnInclude/cuttingPlane.H:58: note: Foam::cuttingPlane::cuttingPlane(const Foam::cuttingPlane&)
make: *** [Make/linux64GccDPOpt/wallGradU.o] Error 1

thank's for your messages

Olivier
brugiere_olivier is offline   Reply With Quote

Old   January 5, 2010, 14:02
Default
  #124
Senior Member
 
Join Date: Dec 2009
Posts: 112
Rep Power: 7
heavy_user is on a distinguished road
nvm................

Last edited by heavy_user; January 7, 2010 at 06:40.
heavy_user is offline   Reply With Quote

Old   February 10, 2010, 16:02
Smile Make setParabolicVelocity work for OF1.5
  #125
Member
 
Hamed Aghajani
Join Date: Mar 2009
Location: London, UK
Posts: 77
Rep Power: 8
haghajani is on a distinguished road
for those who may still interested in setParabolicVelocity

you should do the following two minor correction to make it work with OF-1.5
scalar maxVel(....args()[2])..
word ... args()[1]);

Best,
Hamed
haghajani is offline   Reply With Quote

Old   February 19, 2010, 12:18
Default Error Setting parabolicVelocity as BC in round Inlet
  #126
New Member
 
Join Date: Jan 2010
Posts: 9
Rep Power: 7
cfd_noob is on a distinguished road
Hi folks,

I am using OF 1.6 (stable) and have compiled the parabolicVelocity lib. I am having trouble with the simulation as such.
It inlet is a 2D pipe in the XY-Plane. The direction of the flow is in Z-Direction with an initial velocity U of 6.6 m/s.

The initial "0/U" file looks this way:

INJET
{
type parabolicVelocity;
n (0 0 1);
y (1 1 0);
maxValue 6.6;
value uniform (0 0 0);
}

However, the results I am getting are completely wrong! The U(z) distribution (Z-Plane => X=0 Y=0) at the INJET is not even circular (the INJET is a pipe). On the other hand, the distribution of the U(z) velocity in the Y=0 and X=0, not even not symmetrical but also incoherent.

I read on this thread that this code could also be used to set an initial parabolicVelocity profile for a 2D case, if the patch was 2D. This is my case.

Any help or any ideas why this might no be working?

best regards,
cfd_noob is offline   Reply With Quote

Old   February 23, 2010, 07:36
Default
  #127
New Member
 
Join Date: Jan 2010
Posts: 9
Rep Power: 7
cfd_noob is on a distinguished road
Anyone using this BC with openFoam 1.6?
cfd_noob is offline   Reply With Quote

Old   February 23, 2010, 07:53
Default
  #128
Senior Member
 
Join Date: Dec 2009
Posts: 112
Rep Power: 7
heavy_user is on a distinguished road
Hi Noob,

i have created a "tool" that writes an inlet profile.

In my case the inlet is in the x-y-plane, flowdirection is z.
You need to modify the equations and also change the name of the inlet.
(and pay attention to the flow direction and orientation of your case)

Then you compile it with wmake and you need to source it, then you should be able to execute the tool in the directory of your case.

I attached the files...

regards
Attached Files
File Type: zip setProfiles.zip (4.2 KB, 129 views)
heavy_user is offline   Reply With Quote

Old   February 28, 2010, 09:44
Default
  #129
New Member
 
Join Date: Jan 2010
Posts: 9
Rep Power: 7
cfd_noob is on a distinguished road
Hi Heavy_User,

Thank you, I will try it out!
cfd_noob is offline   Reply With Quote

Old   March 8, 2010, 12:41
Default Parabolic Profile and Lift&Drag Forces in OF-1.6
  #130
New Member
 
Evren
Join Date: Mar 2010
Posts: 20
Rep Power: 7
pbe_cfd is on a distinguished road
Hi,

I am working on benchmark problem "Flow Around Cylinder" and I am a freshman in OF world. I could obtain preliminary results for my geometry. Nevertheless, the results has to be improved by imposing a parabolic velocity on the inlet patch and the lift and drag forces has to be obtained for a quantitative comparison.

My main bottlenecks are
1) How to impose a parabolic profile in OF 1.6.
2) How to obtain lift and drag forces in OF 1.6.

Actually, bottleneck 1 is very well described by Bernhard in his several posts and I liked his way of implementation. Nevertheless, his code is good upto OF 1.4.1 (?). It gives following run time error when it's compiled and executed in OF 1.6.

////////////////////////////////////////////////////////////////////////////////////
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/user/ebayrakt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/home/user/ebayrakt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/user/ebayrakt/OpenFOAM/ebayrakt-1.6/applications/bin/linux64GccDPOpt/setParabolicInlet"
#4 __libc_start_main in "/lib64/libc.so.6"
#5 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/user/ebayrakt/OpenFOAM/ebayrakt-1.6/applications/bin/linux64GccDPOpt/setParabolicInlet"
Segmentation fault
////////////////////////////////////////////////////////////////////////////////////

I tried to debug but it seems it will take ages for me to fix it without external support. So, any suggestion is greatly appreciated.

About bottleneck 2, I've found a thread titled "Forces in v1.6" and I think that it is gonna solve the issue.

Cheers,
Evren

PS: May be, the code is all fine but naming convention or some libraries have been updated and it is just a matter of incompatibility of versions.
pbe_cfd is offline   Reply With Quote

Old   March 8, 2010, 12:46
Default
  #131
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,619
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi Evren

Make a search for groovyBC either here on the Forum or on the wiki, and you will find all you need.

Best regards,

Niels
ngj is offline   Reply With Quote

Old   March 8, 2010, 12:49
Default
  #132
New Member
 
Evren
Join Date: Mar 2010
Posts: 20
Rep Power: 7
pbe_cfd is on a distinguished road
Hi Niels,

I'm gonna check it out and post my feedback.

Thanks,
Evren
pbe_cfd is offline   Reply With Quote

Old   April 28, 2010, 05:11
Default
  #133
New Member
 
Evren
Join Date: Mar 2010
Posts: 20
Rep Power: 7
pbe_cfd is on a distinguished road
Hi all,

groovyBC seems to be a power tool to modify BCs, nevertheless it requires Bison to be installed. In my case, i just needed to implement parabolic BC at the inlet thus I preferred to follow the suggestions in the openfoamwiki (www.openfoamwiki.net). In wiki, under the title "Sig Turbomachinery / Tutorials", it is explained how to implement a parabolic BC in 2D. It's quite straight forward to modify it for 3D. And, if someone needs in 3D, I can send my source files.

Cheers,
evren
pbe_cfd is offline   Reply With Quote

Old   April 28, 2010, 16:27
Default
  #134
New Member
 
Fernando
Join Date: Feb 2010
Posts: 28
Rep Power: 7
fsalvucci is on a distinguished road
I have a problem with a customized 3D parabolic BC. I followed all the steps of the Sig Turbomachinery in wiki, but when i try to run the case, an error occurs. The write function is as follows:

void parabolicVelocityFvPatchVectorField::write(Ostream & os) const
{
fvPatchVectorField::write(os);
os.writeKeyword("maxValue")
<< maxValue_ << token::END_STATEMENT << nl;
os.writeKeyword("n")
<< n_ << token::END_STATEMENT << nl;
os.writeKeyword("y")
<< y_ << token::END_STATEMENT << nl;
os.writeKeyword("freq")
<< f_ << token::END_STATEMENT << nl;
os.writeKeyword("phi")
<< phi_ << token::END_STATEMENT << nl;
writeEntry("value", os);
}


In the 0/U files, i specify the boundary condition as:


inlet
{ type parabolicVelocity;
maxValue 0.452;
n (0 0 -1);
y (0 0 0);
freq 0;
phi 0;
}

And when i run the case, the following error occurs:

Cannot find 'value' entry on patch entrada of field U in file "./SmeriglioSilviaP/0/U"
which is required to set the values of the generic patch field.
(Actual type parabolicVelocity)
Please add the 'value' entry to the write function of the user-defined boundary-condition
or link the boundary-condition into libfoamUtil.so
file: ./SmeriglioSilviaP/0/U::boundaryField::entrada from line 35 to line 40.
From function genericFvPatchField<Type>::genericFvPatchField(con st fvPatch&, const Field<Type>&, const dictionary&)
in file fields/fvPatchFields/basic/generic/genericFvPatchField.C at line 72.

FOAM exiting

So, i add in 0/U an entry after phi:

value (0 0 0); As is sugested in wiki tutorial, and says:

FOAM FATAL ERROR:
gradientInternalCoeffs cannot be called for a genericFvPatchField (actual type parabolicVelocity)
on patch entrada of field U in file "./SmeriglioSilviaP/0/U" You are probably trying to solve for a field with a generic boundary condition.
From function genericFvPatchField<Type>::gradientInternalCoeffs( ) const
in file fields/fvPatchFields/basic/generic/genericFvPatchField.C at line 782.
FOAM exiting



Any ideas? I need help!







fsalvucci is offline   Reply With Quote

Old   April 29, 2010, 09:19
Default
  #135
New Member
 
Evren
Join Date: Mar 2010
Posts: 20
Rep Power: 7
pbe_cfd is on a distinguished road
Quote:
Originally Posted by fsalvucci View Post

FOAM FATAL ERROR:

gradientInternalCoeffs cannot be called for a genericFvPatchField (actual type parabolicVelocity)
on patch entrada of field U in file "./SmeriglioSilviaP/0/U" You are probably trying to solve for a field with a generic boundary condition.
From function genericFvPatchField<Type>::gradientInternalCoeffs( ) const
in file fields/fvPatchFields/basic/generic/genericFvPatchField.C at line 782.
FOAM exiting

Any ideas? I need help!
May it be that you forget to add this line in the controlDict, so to tell OF which library it should use for this type of BCs (<case_dir>/system/) file
libs ("libmyBCs.so");

If you build the dynamic library with another name then you should change "libmyBCs.so" with yours. You can check your own build libraries in this directory: $FOAM_USER_LIBBIN

I hope this fixes your problem
pbe_cfd is offline   Reply With Quote

Old   April 29, 2010, 10:32
Default
  #136
New Member
 
Fernando
Join Date: Feb 2010
Posts: 28
Rep Power: 7
fsalvucci is on a distinguished road
Thanks for the reply. Actually, i added that line in the controlDict file. And everything is as should be (or seems to), the fact is that the library libmyBCs.so isn´t anywhere. Should i create it before? i think that adding in the Make/files file in the BC folder this:

parabolicVelocity.......C
LIB = $(FOAM_USER_LIBBIN)/libmyBCs

Is all i needed to do and that by this, when i compile, it creates the library, is that right?

Thanks!
fsalvucci is offline   Reply With Quote

Old   April 30, 2010, 04:16
Default
  #137
New Member
 
Evren
Join Date: Mar 2010
Posts: 20
Rep Power: 7
pbe_cfd is on a distinguished road
Quote:
Originally Posted by fsalvucci View Post
And everything is as should be (or seems to), the fact is that the library libmyBCs.so isn´t anywhere.
Then, the necessary dynamic library is not build. Then, you have to compile and build it.

And if you compile the source code as it's described in the wiki page, it should create libmyBC.so in this directory $(FOAM_USER_LIBBIN).

Quote:
Originally Posted by fsalvucci View Post
Should i create it before? i think that adding in the Make/files file in the BC folder this:
The make file contains only these two lines:
///////////////////////////////////////////////////
parabolicVelocityFvPatchVectorField.C
LIB = $(FOAM_USER_LIBBIN)/libmyBCs
///////////////////////////////////////////////////
So libmyBCs will be created in the specified directory. Of course your source file might have another name then, you should write the proper file name. And this should do the work for you. You build the library, you tell the solver where to find this library. So, this should be all...

Good luck ;-)
Evren
pbe_cfd is offline   Reply With Quote

Old   April 30, 2010, 15:07
Default
  #138
New Member
 
Fernando
Join Date: Feb 2010
Posts: 28
Rep Power: 7
fsalvucci is on a distinguished road
Thanks very much!!!!!! it seems that there where some errors in the code of the BC that didnt allow the compliation and building to finish succesfully and that caused that the library woudlnt be linked. I could run succesfully now, ill see if it is doing what it was supposed to do.
Thanks!
fsalvucci is offline   Reply With Quote

Old   May 14, 2010, 05:15
Default general profile inlet bc
  #139
Member
 
Aldo Iannetti
Join Date: Feb 2010
Posts: 48
Rep Power: 7
aldo.iannetti is on a distinguished road
Hi
I'm a new OF 1.6 user, I'd like to know how to implement an interpolated inlet velocity profile taken from an external data file. Is this possibility already there in OF 1.6?
thanks







Quote:
Originally Posted by hjasak View Post
Hi, I did one of those a while back - it is called parabolicVelocityFvPatchVectorField and should be somewhere on the discussion group. If you'd like the one updated to version 1.3, send me and E-mail and I'll give you the files.

Enjoy,

Hrv
aldo.iannetti is offline   Reply With Quote

Old   August 7, 2010, 02:13
Default
  #140
Member
 
pradeep
Join Date: Jul 2010
Posts: 30
Rep Power: 7
pradeeppandeygbpec is on a distinguished road
hello...all..
i need to prescribe parabolic velocity profile at inlet to a 3-D circular pipe . The pipe is along the x-Axis...ie. x-axis as its centerline.
I know the parabolic velocity profile for cylindrical coordinates system but how to implement the same in Cartesian coordinates....i need "u" as a function of x,y,z at inlet.
plz help me out...
pradeeppandeygbpec 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:49.