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

How to set fully developed flow aty inlet?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 25, 2015, 15:59
Default How to set fully developed flow aty inlet?
  #1
New Member
 
Fernando Soares Alves
Join Date: Oct 2013
Posts: 16
Rep Power: 5
FernandoSoares is on a distinguished road
Hi there Foamers,

I was wondering how to set properly a fully developed velocity profile WITH the pressure gradient required to push the fluid (newtonian). As far as I have checked mostly what has been adopted regards only mapping the velocity profile at the geometry inlet however with zeroGradient for the pressure field.

Is it just me or does these b.c. sound quiet off? Suppose I do know the pressure gradient in order to assure the development of the flow, which b.c. should I adopt at the outlet? I mean for a simple geometry, e.g. a backward facing step for instance.

I was trying something like

U:
Code:
    inlet
    {
        type            fixedValue;
        value           nonuniform List<vector>
     (developed velocity profile here ...)
     }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform (0 0 0);
        value           uniform (0 0 0);
     }
p:
Code:
    inlet
    {
        type            fixedGradient;
        gradient        uniform (~known Pressure Gradient);
    }
    outlet
    {
        type            outletInlet;
        outletValue     uniform 0;
        value           uniform 0;
    }
and I get no good convergence whatsover. (FYI I start the run with a mapped solution obtained from a previous result for a smaller geomtery). A plot of the residuals I am getting for simpleFoam with standard kEpsilon model is attached to this thread.

How is one suppose to set such a boundary condition in order to assure the velocity profile is in fact developed (i.e. with pressure gradient at the inlet)? I've already used zeroGradient for pressure, but it seems to me such condition doesn't agree with definition of developed flow.

Best regards to all and many thanks in advance for any advise.

F. Soares

PS: I found this thread on the subject emphasizing what I mean.
Attached Images
File Type: jpg Screenshot from 2015-09-25 16:46:53.jpg (91.3 KB, 21 views)
FernandoSoares is offline   Reply With Quote

Old   September 28, 2015, 02:57
Default
  #2
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,208
Blog Entries: 1
Rep Power: 17
nimasam is on a distinguished road
instead of using zeroGradient, employing fixedFluxPressure at inlet for pressure will consider your concerns. because fixedFluxPressure sets the pressure gradient based on velocity field
__________________
Telegram channel (https://telegram.me/openfoam4Iranian)
My Weblog (http://openfoam.blogfa.com/)
Training Course on OpenFOAM at (http://www.isme.ir/)
nimasam is offline   Reply With Quote

Old   September 30, 2015, 13:44
Thumbs up
  #3
New Member
 
Fernando Soares Alves
Join Date: Oct 2013
Posts: 16
Rep Power: 5
FernandoSoares is on a distinguished road
Sorry for the late reply, but thanks for your advise. I am trying to use the b.c. you mentioned but I am getting the following error message

Code:
--> FOAM FATAL ERROR: 
updateCoeffs(const scalarField& snGradp) MUST be called before updateCoeffs() or evaluate() to set the boundary gradient.

    From function fixedFluxPressureFvPatchScalarField::updateCoeffs()
    in file fields/fvPatchFields/derived/fixedFluxPressure/fixedFluxPressureFvPatchScalarField.C at line 151.

FOAM exiting
I must say I am not sure how to set it in my 0/p file. According to the damBreak tutorial

Code:
    leftWall
    {
        type            fixedFluxPressure;
        value           uniform 0;
    }

    rightWall
    {
        type            fixedFluxPressure;
        value           uniform 0;
    }

    lowerWall
    {
        type            fixedFluxPressure;
        value           uniform 0;
    }
although these are used at walls not at an inlet.

Are you familiar with the proper setting for fixedFluxPressure? What am I missing here?

Once again thank you for your help.

F. Soares
FernandoSoares is offline   Reply With Quote

Old   November 17, 2016, 05:18
Default
  #4
New Member
 
vidyadhar
Join Date: Jul 2016
Posts: 28
Rep Power: 2
vidyadhar is on a distinguished road
Hello Soares,

I am new to openfoam and I would like to simulate fully developed flow through a rectangular channel.

Since you have already tried, can you, please, help me in setting the velocity field for fully developed case.

Thanks & Regards,
vidyadhar
vidyadhar is offline   Reply With Quote

Old   November 17, 2016, 09:40
Default
  #5
New Member
 
Fernando Soares Alves
Join Date: Oct 2013
Posts: 16
Rep Power: 5
FernandoSoares is on a distinguished road
Quote:
Originally Posted by vidyadhar View Post
Hello Soares,

I am new to openfoam and I would like to simulate fully developed flow through a rectangular channel.

Since you have already tried, can you, please, help me in setting the velocity field for fully developed case.

Thanks & Regards,
vidyadhar
Hi there, Vidyadhar.

Perhaps what you need for such a geometry is simply the boundaryFoam solver. If I understood correctly, your problem is to find the velocity profile for a developed flow corresponding to a rectangular cross-sectioned pipe. Right? Check the tutorials in
Code:
$FOAM_TUTORIALS/incompressible/boundaryFoam
I believe you'll find what you're looking for in there!

Best regards,

Fernando Soares Alves.

PS: Check also this .pdf which summarizes the most used boundary conditions in the former OpenFOAM version 2.3.0. I hope it will also help you.
FernandoSoares is offline   Reply With Quote

Old   November 17, 2016, 10:29
Default
  #6
New Member
 
vidyadhar
Join Date: Jul 2016
Posts: 28
Rep Power: 2
vidyadhar is on a distinguished road
Hi Fernando,

thank you for the quick reply.

Actually, I want to give inlet condition as fully developed flow.

I want to study the flow behavior through a T junction (either with 1 inlet & 2 outlets or with 2 inlets & 1 outlet).

I am using Caelus6.04 version, which is similar to openfoam. I have openfoam 4.1 also.

I find in some forums, for defining fully developed flow at inlet: some are referring to groovyBCs; some are directly modifying their 0/U file.

As you have tried to use fully developed flow conditions at inlet for your geometry, can you please let me know how to implement for my case.

Thanks & Regards,
vidyadhar
vidyadhar is offline   Reply With Quote

Old   November 17, 2016, 13:02
Default
  #7
New Member
 
Fernando Soares Alves
Join Date: Oct 2013
Posts: 16
Rep Power: 5
FernandoSoares is on a distinguished road
Hi, Vidyadhar.

Actually I couldn't make it work. I've tried this fixedFluxPressure as mentioned above but kept getting the error message. Searching for this BC again today I've also found this blog entry. Ok, it's in japanese but Google translate is here for those sort of things. Besides you can also check the .pdf I mentioned to see the correct sintaxe for such BC.

Greetings and let me know if you could get it to work for your problem.

Fernando Soares Alves.
FernandoSoares is offline   Reply With Quote

Old   November 18, 2016, 01:16
Default
  #8
New Member
 
vidyadhar
Join Date: Jul 2016
Posts: 28
Rep Power: 2
vidyadhar is on a distinguished road
Hello Fernando,

Thank you very much.

I will let you know once I get any output.

Regards,
vidyadhar
vidyadhar 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
Couple outlet to inlet to achieve a fully developed flow. raghu.tejaswi Fluent UDF and Scheme Programming 8 July 18, 2013 09:04
fully developed flow simualtion with fluent bryant_k FLUENT 2 January 1, 2013 12:19
Nusselt numbers for fully developed laminar flow Danro CD-adapco 8 July 19, 2012 07:22
LRR Rstm pipe flow takes too long to get fully developed kjetil OpenFOAM Running, Solving & CFD 0 February 13, 2010 12:26
fully developed flow JAS FLUENT 4 February 12, 2007 07:01


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