CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Newbie: advantages/disadvantages of "penalty formulation" incompressible fluid flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 30, 2011, 15:32
Default Newbie: advantages/disadvantages of "penalty formulation" incompressible fluid flow
  #1
New Member
 
Join Date: Jan 2010
Posts: 28
Rep Power: 7
bzz77 is on a distinguished road
Hello all:

I'm a chemist, interested in getting into some fluid flow modeling (e.g., advection-diffusion). I'm a newbie, so please forgive any silly questions.

Today I came across something new (the "penalty formulation") and I'd really appreciate some basic help or directions to online references in understanding (1) how it works and (2) what the advantages/disadvantages are compared to the method I'm familiar with--using conservation of fluid mass to derive an expression for the fluid pressure p and then getting the velocity profile from that.

Background:
From what I've read, using the penalty method allows elimination of pressure from the momentum equation so that the mass and momentum equations can be combined. The result is a convection-diffusion equation where velocity can then be calculated directly.

This seems kind of weird to me though. Is it a bit hokey to do this? What approximations are involved? I deal mainly with multiphase phenomena, but I so want to work on some single-phase flow too... Could it be that the penalty formulation stuff might be better suited to single-phase flow?

The first step in using the penalty method is to replace the continuity equation by:
1/e + grad u ...

I don't understand what e is. Where does it come from? The book I have just calls it a "coefficient with large values."

I guess what I'm getting at is, what ARE the advantages of the penalty method? I see that someone I respect has used it and I'm wondering why he would have chosen it over the pressure calculation method. Any advice would be very much appreciated.
bzz77 is offline   Reply With Quote

Old   November 1, 2011, 00:39
Default
  #2
Super Moderator
 
praveen's Avatar
 
Praveen. C
Join Date: Mar 2009
Location: Bangalore
Posts: 244
Blog Entries: 6
Rep Power: 9
praveen is on a distinguished road
I have seen formulations in which

div(u) = 0

is replaced with

eps*p + div(u) = 0

Here eps is a small number. I think the advantage is that you dont have to solve a pressure poisson equation. You eliminate p from the momentum equation, solve for velocity, and then get p as

p = -(1/eps)*div(u)

For an analysis of the error introduced in this method see

http://www.math.purdue.edu/~shen/pub/Penalty.pdf

I am told that purists frown at such a method.
praveen is offline   Reply With Quote

Old   November 1, 2011, 04:32
Default
  #3
New Member
 
Join Date: Jan 2010
Posts: 28
Rep Power: 7
bzz77 is on a distinguished road
Hi Praveen:

Thanks a lot for your help!

It is very helpful to know that not solving the pressure equation is why this method is used, and that purists don't like it.

Thanks very much for the reference.
bzz77 is offline   Reply With Quote

Reply

Tags
diffusion, finite difference, finite element method, finite volume method, penalty formulation

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
Incompressible and compressible flow. Confused student. Main CFD Forum 22 January 30, 2014 07:32
Constant flow rate through a small area inside the fluid domain. robingilbert OpenFOAM 7 October 4, 2010 16:19
Turbulent flow of two incompressible fluids with the same physical properties ArneBab OpenFOAM Running, Solving & CFD 3 January 5, 2010 06:48
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
internal steady incompressible inviscid fluid flow Nina Shokina Main CFD Forum 0 October 21, 2003 06:11


All times are GMT -4. The time now is 22:47.