[alberto]

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

[torvic]

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]

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

[alberto]

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

[torvic]

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

[msrinath80]

I just finished compiling Bernhard's setParabolicInlet.C on OpenFOAM 1.4.1. The missing header file (fvPatchFieldFields.H) is attached.

Attachments:
fvPatchFieldFields.H

[nzy102]

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

[hjasak]

Try this one (well, two):

parabolicVelocityFvPatchVectorField.H

parabolicVelocityFvPatchVectorField.C

Enjoy,

Hrv

[nzy102]

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]

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]

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

[hjasak]

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

[nzy102]

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

[msrinath80]

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.

[hjasak]

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

[nzy102]

Dr. Jasak:

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

Ning

[armin_h]

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]

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

[cedric_duprat]

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

[inddzen]

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