# Solution of Navier-Stokes equations

(Difference between revisions)
 Revision as of 23:08, 20 December 2005 (view source)Jola (Talk | contribs) (removed extra spaces (we should not use extra line-feeds like this...))← Older edit Latest revision as of 10:01, 23 April 2009 (view source)Jola (Talk | contribs) m (Reverted edits by Diyaseb (Talk) to last version by Peter) (5 intermediate revisions not shown) Line 1: Line 1: - For the incompressible flows, the Navier-Stokes equation could be written in the form:
+ ==Introduction== + Most traditional CFD algorithms require the solution of the Navier-Stokes (N-S) equations.  There are two important issues that arise in the solution process: + # The N-S equations are nonlinear. + # The pressure-velocity coupling often requires special consideration. + + The details of the solution process depend upon the details of the flow to be solved.  The solution process for a incompressible flow can be very different than that for a compressible flow. + + ==Incompressible flow== + In situations in which the density is approximately constant, the flow may be termed [[Incompressible flow |incompressible]].  The Navier-Stokes equation may then be written as a continuity equation + :$:[itex] - \nabla \cdot \vec U = 0 + \nabla \cdot \vec U = 0$ - :${{\partial \vec U} \over {\partial t}} + \nabla \cdot \left( {\vec U\vec U} \right) - \nabla \cdot\left( {\nu \nabla \vec U} \right) = - \nabla p$ + + and as a momentum equation - There are two important issues regarding Navier-Stokes equations:
+ : {{\partial \vec U} \over {\partial t}} + \nabla  \cdot \left( {\vec U\vec U} \right) =  - \nabla p + \nabla  \cdot\left( {\nu \nabla \vec U} \right). [/itex] - # Non linearity of momentum equations + - # Pressure-velocity coupling + - ==Segregated Solver== + Here, we have written the nonlinear convective term (the second term on the right-hand side) in the so-called conservative form.  If a turbulence model is to be employed, then the equations will change.  If an eddy viscosity approach is to be used, then there are three likely modifications: - === The solution scheme === + *the flow variables will represent average (or filtered) quantities, + *the viscosity $\nu$ will actually be the sum of the fluid property and the calculated eddy viscosity (and will then be nonconstant, justifying keeping it inside the divergence operator), and + *the pressure will be modified to include normal-stress-like terms arising from the eddy viscosity assumption. - # Solve Momentum equations (u,v,w) + Thus, the equation as written will be valid for both turbulent or laminar flows, with some modification of the actual meaning of individual terms. - # Solve pressure correction equation + - ## Correct fluxes and velocities + - # Solve transport equations for other scalars + - == Coupled Solver == + It is somewhat difficult to categorize the solution methods for the incompressible N-S equations, as over the years many have been developed that are very similar.  However, the manner in which the difficulties associated with the pressure solution are handled (a process made difficult by the fact that pressure does not appear in the continuity equation) can be used to broadly categorize methods. - # Solve the Momentum equations- Pressure equation system in one go (u,v,w,p) + - # Solve transport equations for other scalars + - ---- + ==Compressible Flow== - Return to [[Numerical methods | Numerical Methods]] +

## Introduction

Most traditional CFD algorithms require the solution of the Navier-Stokes (N-S) equations. There are two important issues that arise in the solution process:

1. The N-S equations are nonlinear.
2. The pressure-velocity coupling often requires special consideration.

The details of the solution process depend upon the details of the flow to be solved. The solution process for a incompressible flow can be very different than that for a compressible flow.

## Incompressible flow

In situations in which the density is approximately constant, the flow may be termed incompressible. The Navier-Stokes equation may then be written as a continuity equation

$\nabla \cdot \vec U = 0$

and as a momentum equation

${{\partial \vec U} \over {\partial t}} + \nabla \cdot \left( {\vec U\vec U} \right) = - \nabla p + \nabla \cdot\left( {\nu \nabla \vec U} \right).$

Here, we have written the nonlinear convective term (the second term on the right-hand side) in the so-called conservative form. If a turbulence model is to be employed, then the equations will change. If an eddy viscosity approach is to be used, then there are three likely modifications:

• the flow variables will represent average (or filtered) quantities,
• the viscosity $\nu$ will actually be the sum of the fluid property and the calculated eddy viscosity (and will then be nonconstant, justifying keeping it inside the divergence operator), and
• the pressure will be modified to include normal-stress-like terms arising from the eddy viscosity assumption.

Thus, the equation as written will be valid for both turbulent or laminar flows, with some modification of the actual meaning of individual terms.

It is somewhat difficult to categorize the solution methods for the incompressible N-S equations, as over the years many have been developed that are very similar. However, the manner in which the difficulties associated with the pressure solution are handled (a process made difficult by the fact that pressure does not appear in the continuity equation) can be used to broadly categorize methods.