CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Why is fixedFluxPressure more stable ?

Register Blogs Community New Posts Updated Threads Search

Like Tree11Likes
  • 1 Post By vatavuk
  • 3 Post By Yann
  • 7 Post By hjasak

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 29, 2022, 01:26
Default Why is fixedFluxPressure more stable ?
  #1
New Member
 
HR
Join Date: Nov 2022
Posts: 17
Rep Power: 3
hr_nt is on a distinguished road
I recently an few multiphase scenarios using interFoam and found that fixedFluxPressure is more stable for the simulation than zeroGradient. What might be the issue for this ?
hr_nt is offline   Reply With Quote

Old   December 2, 2022, 18:17
Default
  #2
Senior Member
 
Paulo Vatavuk
Join Date: Mar 2009
Location: Campinas, Brasil
Posts: 196
Rep Power: 17
vatavuk is on a distinguished road
Hi HR,

At the 11th OpenFOAM Workshop there was a Training Course with the title "Learning how to use free surface flows in OpenFOAM 3.0" by Victoria Korchagova. This course is the only material that I found that explains this boundary condition. Is it for OF version 3, but I think that it also applies to newer versions. You can download the course material in this page:
https://openfoam-extend.sourceforge....s/courses.html
hr_nt likes this.
vatavuk is offline   Reply With Quote

Old   December 5, 2022, 03:48
Default
  #3
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 1,074
Rep Power: 26
Yann will become famous soon enough
The OpenFOAM user guide also mention it, despite not giving much details:

Quote:
5.2.3.3 Fixed flux pressure

In the above example, it can be seen that all the wall boundaries use a boundary condition named fixedFluxPressure. This boundary condition is used for pressure in situations where zeroGradient is generally used, but where body forces such as gravity and surface tension are present in the solution equations. The condition adjusts the gradient accordingly.
https://doc.cfd.direct/openfoam/user...v10/boundaries

There are also on old thread here discussing it: I need explanations about fixedFluxPressure

My personal experience with buoyantSimpleFoam is that I get a flowrate through the walls when using zeroGradient and it does not happen with fixedFluxPressure.

Yann
hjasak, acgnipper and hr_nt like this.
Yann is offline   Reply With Quote

Old   December 10, 2022, 09:09
Default
  #4
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,905
Rep Power: 33
hjasak will become famous soon enough
Hi,

This is because when you’re solving for p, zero gradient isn’t the correct boundary condition. Consider additional forces, such as gravity: if you zero gradient the pressure, such terms will give you flux through a wall.

Fixed flux pressure compensates for it to achieve zero flux: you get a pressure gradient on the wall to balance all other flux contributions.

Hrv
vatavuk, jherb, Yann and 4 others like this.
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Decomposing error simple method DanGode OpenFOAM Pre-Processing 1 July 27, 2021 06:58
fixedFluxPressure BC: updateCoeffs(const scalarField& snGradp) MUST be called ... MCrossover97 OpenFOAM Running, Solving & CFD 1 June 25, 2021 06:03
bugs of buoyantBoussinesqPisoFoam or fixedFluxPressure in foam-extend 3.1 and 3.2 Aaron_L OpenFOAM Bugs 5 July 30, 2016 07:11
fixedFluxPressure with PisoFoam me.ouda OpenFOAM Running, Solving & CFD 1 October 27, 2015 05:30
rhoPimpleFoam with fixedFluxPressure + flowRateInletVelocity --> p diverges tatu OpenFOAM Running, Solving & CFD 1 March 21, 2013 15:10


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