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   August 14, 2007, 15:58
Default This library is the implementa
  #41
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
This library is the implementation of a parabolic velocity profile from Hrvoje which works under OF 1.4. It currently doesn't compile under 1.4.1.

parabolicVelocity.tar.gz

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   August 14, 2007, 16:26
Default Hi Alberto I do appreciate
  #42
Member
 
victor
Join Date: Mar 2009
Location: mexico city, MX
Posts: 50
Rep Power: 8
torvic is on a distinguished road
Hi Alberto

I do appreciate your fast response. thanks for both the fact about OF-1.4.1. and the file,

I'll try it and tell you,

thanks again

best

V
torvic is offline   Reply With Quote

Old   August 14, 2007, 21:26
Default Hi again Alberto, Now i can
  #43
Member
 
victor
Join Date: Mar 2009
Location: mexico city, MX
Posts: 50
Rep Power: 8
torvic is on a distinguished road
Hi again Alberto,

Now i can compile the BC. and added it to "U" in the dictionary. Thanks.

One last question, i want to make sure that the profile is taken into account in the computations. How can i "see" it. sorry if it is a naive question

thanks in advance

best

V
torvic is offline   Reply With Quote

Old   August 14, 2007, 22:50
Default Hi Victor, you can see it in
  #44
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
Hi Victor,
you can see it in two ways:

- use paraFoam to visualize the boundary condition. You should see the change in the velocity magnitude.

- extract data using the sample utility, by using a straight line across the BC. You can find how to use this utility in the User's manual.

P.S. Did you compile it using OF 1.4.1 or did you go back to 1.4?

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   August 15, 2007, 13:32
Default Hi Alberto, Thanks so much
  #45
Member
 
victor
Join Date: Mar 2009
Location: mexico city, MX
Posts: 50
Rep Power: 8
torvic is on a distinguished road
Hi Alberto,

Thanks so much again for your explanation. I will try those and tell you.

I compiled it under OF-1.4, since i followed your advice. However, i tried just in case, but it didn't work. I got these error messages (It's the final part, since it's a bit long)

-----------------
/home/foam-1.4.1/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFv PatchField.C:41: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fvPatch&, const Foam::DimensionedField<type,>&) [with Type = Foam::Vector<double>]
make: *** [Make/linuxGccDPOpt/parabolicVelocityFvPatchVectorField.o] Error 1
-----------------

Nevertheless, after compiled it in OF-1.4, i followed the manual and compiled it again under "foamUser". At this moment reactingFoam is running with no problem.

Thanks again Alberto, !

best

V
torvic is offline   Reply With Quote

Old   August 31, 2007, 14:32
Default I just finished compiling Bern
  #46
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
I just finished compiling Bernhard's setParabolicInlet.C on OpenFOAM 1.4.1. The missing header file (fvPatchFieldFields.H) is attached.

Attachments:
fvPatchFieldFields.H
msrinath80 is offline   Reply With Quote

Old   September 6, 2007, 15:29
Default Hi Alberto and Hrvoje: I tr
  #47
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 8
nzy102 is on a distinguished road
Hi Alberto and Hrvoje:

I tried to compile your parabolic velocity under openfoam 1.4.1 and got error message as follows:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[nzy102@lionxo foamUser]$ wmake libso
Making dependency list for source file parabolicVelocityFvPatchVectorField.C
SOURCE=parabolicVelocityFvPatchVectorField.C ; g++ -m64 -Dlinux64 -DDP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -march=opteron -O3 -DNoRepository -ftemplate-depth-40 -I/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude -IlnInclude -I. -I/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOpt/parabolicVelocityFvPatchVectorField.o
parabolicVelocityFvPatchVectorField.C: In constructor ‚Foam::parabolicVelocityFvPatchVectorField::parabo licVelocityFvPatchVectorField (const Foam::fvPatch&, const Foam::vectorField&)‚:
parabolicVelocityFvPatchVectorField.C:49: error: no matching function for call to ‚Foam::fixedValueFvPatchField<foam::vector<double> >::fixedValueFvPatchField(const Foam::fvPatch&, const Foam::Field<foam::vector<double> >&)‚
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:87: note: candidates are: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fixedValueFvPatchField<type>&, const Foam::DimensionedField<type,>&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:76: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fixedValueFvPatchField<type>&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:66: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fixedValueFvPatchField<type>&, const Foam::fvPatch&, const Foam::DimensionedField<type,>&, const Foam::fvPatchFieldMapper&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:53: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fvPatch&, const Foam::DimensionedField<type,>&, const Foam::dictionary&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:41: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fvPatch&, const Foam::DimensionedField<type,>&) [with Type = Foam::Vector<double>]
parabolicVelocityFvPatchVectorField.C: In constructor ‚Foam::parabolicVelocityFvPatchVectorField::parabo licVelocityFvPatchVectorField (const Foam::parabolicVelocityFvPatchVectorField&, const Foam::fvPatch&, const Foam::vectorField&, const Foam::fvPatchFieldMapper&)‚:
parabolicVelocityFvPatchVectorField.C:64: error: no matching function for call to ‚Foam::fixedValueFvPatchField<foam::vector<double> >::fixedValueFvPatchField(const Foam::parabolicVelocityFvPatchVectorField&, const Foam::fvPatch&, const Foam::Field<foam::vector<double> >&, const Foam::fvPatchFieldMapper&)‚
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:87: note: candidates are: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fixedValueFvPatchField<type>&, const Foam::DimensionedField<type,>&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:76: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fixedValueFvPatchField<type>&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:66: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fixedValueFvPatchField<type>&, const Foam::fvPatch&, const Foam::DimensionedField<type,>&, const Foam::fvPatchFieldMapper&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:53: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fvPatch&, const Foam::DimensionedField<type,>&, const Foam::dictionary&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:41: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fvPatch&, const Foam::DimensionedField<type,>&) [with Type = Foam::Vector<double>]
parabolicVelocityFvPatchVectorField.C: In constructor ‚Foam::parabolicVelocityFvPatchVectorField::parabo licVelocityFvPatchVectorField (const Foam::fvPatch&, const Foam::vectorField&, const Foam::dictionary&)‚:
parabolicVelocityFvPatchVectorField.C:78: error: no matching function for call to ‚Foam::fixedValueFvPatchField<foam::vector<double> >::fixedValueFvPatchField(const Foam::fvPatch&, const Foam::Field<foam::vector<double> >&)‚
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:87: note: candidates are: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fixedValueFvPatchField<type>&, const Foam::DimensionedField<type,>&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:76: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fixedValueFvPatchField<type>&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:66: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fixedValueFvPatchField<type>&, const Foam::fvPatch&, const Foam::DimensionedField<type,>&, const Foam::fvPatchFieldMapper&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:53: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fvPatch&, const Foam::DimensionedField<type,>&, const Foam::dictionary&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:41: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fvPatch&, const Foam::DimensionedField<type,>&) [with Type = Foam::Vector<double>]
parabolicVelocityFvPatchVectorField.C: In constructor ‚Foam::parabolicVelocityFvPatchVectorField::parabo licVelocityFvPatchVectorField (const Foam::parabolicVelocityFvPatchVectorField&, const Foam::vectorField&)‚:
parabolicVelocityFvPatchVectorField.C:103: error: no matching function for call to ‚Foam::fixedValueFvPatchField<foam::vector<double> >::fixedValueFvPatchField(const Foam::parabolicVelocityFvPatchVectorField&, const Foam::Field<foam::vector<double> >&)‚
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:87: note: candidates are: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fixedValueFvPatchField<type>&, const Foam::DimensionedField<type,>&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:76: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fixedValueFvPatchField<type>&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:66: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fixedValueFvPatchField<type>&, const Foam::fvPatch&, const Foam::DimensionedField<type,>&, const Foam::fvPatchFieldMapper&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:53: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fvPatch&, const Foam::DimensionedField<type,>&, const Foam::dictionary&) [with Type = Foam::Vector<double>]
/home1/nzy102/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fixedValueFvPat chField.C:41: note: Foam::fixedValueFvPatchField<type>::fixedValueFvPa tchField(const Foam::fvPatch&, const Foam::DimensionedField<type,>&) [with Type = Foam::Vector<double>]
make: *** [Make/linux64GccDPOpt/parabolicVelocityFvPatchVectorField.o] Error 1
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

can you give me some hints what is wrong? I really have to apply parabolic profile in OF 1.4.1. Thank you.

Ning
nzy102 is offline   Reply With Quote

Old   September 6, 2007, 18:20
Default Try this one (well, two):
  #48
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Try this one (well, two):

parabolicVelocityFvPatchVectorField.H

parabolicVelocityFvPatchVectorField.C

Enjoy,

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

Old   September 6, 2007, 21:13
Default Dr. Jasak: It compiles fine
  #49
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 8
nzy102 is on a distinguished road
Dr. Jasak:

It compiles fine now. When I ran my case, I got:

=================================================--> FOAM FATAL ERROR :
gradientInternalCoeffs cannot be called for a defaultFvPatchField (actual type parabolicVelocity)
on patch branch_right of field U in file "/home1/nzy102/OpenFOAM/nzy102-1.4.1/run/tutorials/simpleFoam/lvad_90_25_25_grid 2_corr1/0/U"
You are probably trying to solve for a field with a default boundary condition.

From function defaultFvPatchField<type>::gradientInternalCoeffs( ) const
in file fields/fvPatchFields/basic/default/defaultFvPatchField.C at line 694.
=================================================

I searched on the forum and one of your posts says: This means that the parabolic boundary condition is not picked up correctly. And you also mentioned you have written a long message on compiling and liking boundary conditions. I have searched on the forum very thoroughly but couldn't find it. Do you think you can post it again? Thank you.

Ning
nzy102 is offline   Reply With Quote

Old   September 6, 2007, 23:11
Default Dr. Jasak: One more questio
  #50
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 8
nzy102 is on a distinguished road
Dr. Jasak:

One more question: can you explain variables of "n" and "maxValue" in parabolicVelocityFvPatchVectorField.C? What format should I use to define n and maxValue in boundary condition as shown below?

inlet
{
type parabolicVelocity;
n ?;
maxValue ?;
}

Thanks a lot.

Ning
nzy102 is offline   Reply With Quote

Old   September 7, 2007, 00:34
Default I think I know n and maxValue
  #51
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 8
nzy102 is on a distinguished road
I think I know n and maxValue now. The question is about y. How to define y in boundary condition? I am working on a circular inlet which is not parallel to either x axis, or y axis, or z. I defined n as a unit vector perpendicular to the inlet surfce while maxValue is the value on the center of the surface, not quite sure about y. In your code, I saw the definition of coord as follows:

=================================================

scalarField coord = ((c - ctr) & y_)/((bb.max() - bb.min()) & y_);

=================================================

Coord has to be equal to r/R (R is radius of the inlet). If y_ is a constant vector, I don't see this happening. Can you explain this a little bit more? I really appreciate it.

Ning
nzy102 is offline   Reply With Quote

Old   September 7, 2007, 02:45
Default Run-time failure: you have to
  #52
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Run-time failure: you have to link the .o file coming out of the compilation with your executable and all will be well. Currently, the top-level code does nto recognise the name of your boundary condition and creates a default type instead.

Coordinate: OK, let's try to draw a parabola in space - it can be oriented in an arbitrary manner:

n - direction in which the velocity vector will point. That should be clear, right?

Now, in order to draw a parabola I have to go from value zero, up to max value and back to zero, OK? The RANGE over which I need to do thta is also clear: from one end of the patch that holds the b.c. up to the max value and down back to zero at the other end. Clear?

Finally, take a 2-D patch aligned with the x-y plane (i.e. its normal is the z-axis, meaning n = (0 0 1)). You can still draw a parabola in 2 ways: from zero to max and back to zero along the X-AXIS or the same along the Y-AXIS. Get it? The y tells you along which vector to walk to go from zero to max and back to zero (and the other axis will keep the same value).

Why don't you just try it out and see what happens?

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

Old   September 8, 2007, 23:55
Default Hey Hrvojie: I followed you
  #53
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 8
nzy102 is on a distinguished road
Hey Hrvojie:

I followed your suggestions and it turns out that the calculated flow rate based on the parabolic profile is almost two times as big as what it is supposed to be. So I made some changes in your parabolicVelocity files. Now the calculated flow rate is right. But the problem is that when I ran it in parallel with two nodes, it is getting even slower and the calcualted flow rate is pretty off from the actual value. So I am just wondering if there will be any problem with the code in the parallel mode. Thanks.





Ning
nzy102 is offline   Reply With Quote

Old   September 9, 2007, 02:29
Default Isn't the utility merely for s
  #54
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
Isn't the utility merely for setting a Boundary condition? In other words all it should do is enter the 0/ directory, open the U file, find the specified inlet patch and change the uniform value to a non-uniform list based on the coded parabola equation. If that is true then this should have no connection with parallel mode. Because once the 0/U file is modified, you use decomposePar and that should take care of processor assignment.
msrinath80 is offline   Reply With Quote

Old   September 9, 2007, 08:25
Default Well, it may have a connection
  #55
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Well, it may have a connection with parallelism because I need a global bounding box for a patch (and haven't done that). In any case, it's a Mickey Mouse modification and will be done at some stage (more important stuff to do now...)

I'm sure a competent user can do it without help.

As for the other comment (rubbisH), I explicitly say that the number in the boundary field is the peak velocity rather than a flow rate. If you want a flow rate, you either add more parameters or write your own b.c. which will give you some more interesting parallelisation issues.

In any case, I have expected a "thank you" but it does not seem to be forthcoming...

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

Old   September 10, 2007, 15:19
Default Dr. Jasak: I found a workar
  #56
Member
 
Ning Yang
Join Date: Mar 2009
Location: University Park, PA, USA
Posts: 85
Rep Power: 8
nzy102 is on a distinguished road
Dr. Jasak:

I found a workaround with my boundary conditions. Thank you for your help and I really appreciate it.

Ning
nzy102 is offline   Reply With Quote

Old   October 30, 2007, 23:45
Default Dear all I am trying to cre
  #57
New Member
 
Armin Hosseinian
Join Date: Mar 2009
Location: Perth, Western Australia, Australia
Posts: 17
Rep Power: 8
armin_h is on a distinguished road
Dear all

I am trying to create the velocity profile of the flow through the pipe.
The graphs which its appear and run with probe, is the same in all length of the pipe and the graph doesn't change in terms of the graphical shape.

I don't know if is there any other suggestion to make velocity profile which shows the velocity profile graph correctly.

I am using the OpenFoam 1.4.

Many thanks in advance
Armin
armin_h is offline   Reply With Quote

Old   December 10, 2007, 05:23
Default Hello all would you please
  #58
New Member
 
Armin Hosseinian
Join Date: Mar 2009
Location: Perth, Western Australia, Australia
Posts: 17
Rep Power: 8
armin_h is on a distinguished road
Hello all

would you please any one let me know how can i push flow into the pipe in openfoam?
I have created the geometry of the pipe and also the mesh.
I need to see the velocity profile which comes from the moving flow into the pipe.
I am new user of openfoam and i dont know how can i put the amounts for the u velocity and see the motion of the fluid inside the pipe.

I would appreciate any comments.

sorry if it is not an advance question.

Cheers
Armin
armin_h is offline   Reply With Quote

Old   December 10, 2007, 05:59
Default Hi Armin, welcome in Open-
  #59
Senior Member
 
Cedric DUPRAT
Join Date: Mar 2009
Location: Belgium
Posts: 179
Rep Power: 8
cedric_duprat is on a distinguished road
Hi Armin,

welcome in Open-Foam world
So, some basic before starting:
1- try, as possible, not to spam the forum by choosing the rignt thread for your question. It will be better for us also to give you the right answer.
2- try also not to send the same message 3 times (for example), because, we will receive all of them, don't worry.

well, after these basics, your answer.
try first to look into the tutorials which are very usefull to start and which explain how OF is working. In the simpleFoam, the .../tutorials/simpleFoam for example (RANS modelling) of the .../channelOodles/channel395 which is a periodic channel.

So you will get how to impose velocity at your inlet and how to make a run.

Then, as BANNARI said, the guides are also uselfull for beguiner. You won't lose your time if you work on that fisrt.

Hope it help

Cedric
cedric_duprat is offline   Reply With Quote

Old   December 11, 2007, 16:53
Default Hi all, I tried to implemen
  #60
New Member
 
inddzen
Join Date: Mar 2009
Posts: 4
Rep Power: 8
inddzen is on a distinguished road
Hi all,

I tried to implement the setParabolicInlet given by Bernhard and after all modifications since I'm using OF-1.4.1 I got the following error:

Making dependency list for source file parabolicInlet.C
SOURCE=parabolicInlet.C ; g++ -m64 -Dlinux64 -DDP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -march=opteron -O3 -DNoRepository -ftemplate-depth-40 -I/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/src/meshTools/lnInclude -I/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude -I/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude -IlnInclude -I. -I/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOpt/parabolicInlet.o
g++ -m64 -Dlinux64 -DDP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -march=opteron -O3 -DNoRepository -ftemplate-depth-40 -I/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/src/meshTools/lnInclude -I/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude -I/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude -IlnInclude -I. -I/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude -fPIC Make/linux64GccDPOpt/parabolicInlet.o -L/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt \
-lfiniteVolume -lOpenFOAM -liberty -ldl -lm -o /home/lhassa/OpenFOAM/lhassa-1.4.1/applications/bin/linux64GccDPOpt/parabolicInl et

/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so: undefined reference to `std::basic_ostream<char,> >& std::__ostream_insert<char,> >(std::basic_ostream<char,> >&, char const*, long)@GLIBCXX_3.4.9'

/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so: undefined reference to `std::basic_istream<char,> >& std::basic_istream<char,> >::_M_extract<float>(float&)@GLIBCXX_3.4.9'

/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so: undefined reference to `std::basic_ostream<char,> >& std::basic_ostream<char,> >::_M_insert<double>(double)@GLIBCXX_3.4.9'

/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so: undefined reference to `std::basic_istream<char,> >& std::basic_istream<char,> >::_M_extract<double>(double&)@GLIBCXX_3.4.9'

/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so: undefined reference to `std::basic_ostream<char,> >& std::basic_ostream<char,> >::_M_insert<long>(long)@GLIBCXX_3.4.9'

/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so: undefined reference to `std::basic_ostream<char,> >& std::basic_ostream<char,> >::_M_insert<unsigned>(unsigned long)@GLIBCXX_3.4.9'

/home/lhassa/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so: undefined reference to `std::basic_ostream<char,> >& std::basic_ostream<char,> >::_M_insert<long>(long long)@GLIBCXX_3.4.9'

collect2: ld returned 1 exit status

make: *** [/home/lhassa/OpenFOAM/lhassa-1.4.1/applications/bin/linux64GccDPOpt/parabolicIn let] Error 1

I'm working on Fedora8 X86_64, and I'm using the gcc-4.2.1.

Does anyone have any hints... and thanks in adavance
inddzen 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 14:10.