Scribd
Upload
49 views

Ppe PDF

The document discusses deriving the Poisson equation for pressure in incompressible fluid flow. It shows that by taking the divergence of the momentum equation and using the continuity equation, this results in a Poisson equation that couples pressure and velocity. In numerical simulations, the Poisson equation is used to solve for pressure and enforce continuity after solving the momentum equation, since numerical errors can cause the velocity field to not be fully divergence-free.

Uploaded by

abshaw
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
49 views

Ppe PDF

The document discusses deriving the Poisson equation for pressure in incompressible fluid flow. It shows that by taking the divergence of the momentum equation and using the continuity equation, this results in a Poisson equation that couples pressure and velocity. In numerical simulations, the Poisson equation is used to solve for pressure and enforce continuity after solving the momentum equation, since numerical errors can cause the velocity field to not be fully divergence-free.

Uploaded by

abshaw
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

1. Poisson Equation for Pressure — The Visual Room http://thevisualroom.com/poisson_for_pressure.

html

1. Poisson Equation for Pressure

For compressible flow, pressure and velocity can be coupled with the Equation of State. But for incompressible flow, there is no obvious way to
couple pressure and velocity.

So, take the divergence of the momentum equation and use the continuity equation to get a Poisson equation for pressure.

1.1. Poisson’s Equation in Continuous Domain

1.1.1. Navier Stokes Equations in Vector Notation

1.1.1.1. The Continuity Equation for an incompressible fluid:

Using the divergence operator (div):

Or using the nabla operator :

(1)

1.1.1.2. The Momentum Equation for an incompressible fluid:

Using the gradient (grad) and Laplacian operators ( ):

Using the nabla operator and Laplacian operators ( ):

(2)

Mass conservation (1) for is a kind of kinematic constraint to the momentum equation (2)
PROBLEM: There is no obvious way to couple the velocity and pressure (for compressible fluids we would have an equation of state, which
provides a relation between density and pressure)

1.1.1.3. Continuity Equation

1.1.1.4. Momentum Equations

1.1.2. The Divergence of the Momentum Equation in Partial Notation

Now take the divergence of the momentum equations, if is the “vector” of momentum equations, the divergence is:

1 of 3 3/27/2018, 2:22 PM
1. Poisson Equation for Pressure — The Visual Room http://thevisualroom.com/poisson_for_pressure.html

1.1.2.1. Add the LHS

Re-arrange:

Apply Continuity, so . Hence:

1.1.2.2. Add the RHS

Re-arrange:

Apply Continuity, so . Hence:

1.1.3. The Poisson Equation in Vector Notation

1.1.3.1. Equate LHS and RHS

1.1.3.2. In Vector Form

A Poisson Equation for pressure, which ensures that continuity is satisfied.


Now pressure and velocity are coupled in the continuous domain.

1.2. Poisson’s Equation in Numerical Domain

We have shown that the Poisson’s Equation is valid in the continuous domain, but in the numerical domain, we use discretisation

1.2.1. Momentum Equations in Vector Form

1.2.1.1. The Momentum Equation for an incompressible fluid:

2 of 3 3/27/2018, 2:22 PM
1. Poisson Equation for Pressure — The Visual Room http://thevisualroom.com/poisson_for_pressure.html

1.2.2. Discretised Momentum Equations in Vector Form

Discretise in time: FD in time, with pressure at time (the pressure that corresponds with the velocity at )

1.2.3. The Divergence of the Momentum Equations in Vector Form

Now take the divergence of the momentum equations, if is the “vector” of momentum equations, the divergence is:

In the numerical scheme, we want a divergence free velocity at the next step, i.e. to satisfy continuity
But at the current step we may have due to numerical errors
So we can’t cancel out the divergence of velocity in the numerical domain (although we could in the continuous domain)

This is a Fractional Step approach:

Solve the momentum equation for velocity, numerically (the divergence of the velocity might not be zero)
Solve the Poisson equation for pressure, forcing the divergence to be zero
Correct the velocity to satisfy continuity

1.2.4. The Poisson Equation in Vector Notation

Poisson Equation for at time and forcing

In the numerical domain the velocity field we are producing with the Navier Stokes equations is not completely divergence free
Think of the velocity obtained from the Navier Stokes as being at an intermediate step
And
We need so that continuity is satisfied

3 of 3 3/27/2018, 2:22 PM

You might also like