The Navier-Stokes equations consists of a time-dependent continuity equation for conservation of mass, three time-dependent conservation of momentum equations and a time-dependent conservation of energy equation. There are four independent variables in the problem, the xy, and z spatial coordinates of some domain, and the time t. There are six dependent variables; the pressure p, density r, and temperature T (which is contained in the energy equation through the total energy Et) and three components of the velocity vector; the u component is in the x direction, the v component is in the y direction, and the w component is in the z direction, All of the dependent variables are functions of all four independent variables. The differential equations are therefore partial differential equations and not the ordinary differential equations that you study in a beginning calculus class.

You will notice that the differential symbol is different than the usual “d /dt” or “d /dx” that you see for ordinary differential equations. The symbol “partial” is is used to indicate partial derivatives. The symbol indicates that we are to hold all of the independent variables fixed, except the variable next to symbol, when computing a derivative. The set of equations are:

Continuity: partialr/partialt + partial(r * u)/partialx + partial(r * v)/partialy + partial(r * w)/partialz = 0

X – Momentum: partial(r * u)/partialt + partial(r * u^2)/partialx + partial(r * u * v)/partialy + partial(r * u * w)/partialz = – partialp/partialx

+ 1/Re * { partialtauxx/partialx + partialtauxy/partialy + partialtauxz/partialz}

Y – Momentum: partial(r * v)/partialt + partial(r * u * v)/partialx + partial(r * v^2)/partialy + partial(r * v * w)/partialz = – partialp/partialy

+ 1/Re * { partialtauxy/partialx + partialtauyy/partialy + partialtauyz/partialz}

Z – Momentum: partial(r * w)/partialt + partial(r * u * w)/partialx + partial(r * v * w)/partialy + partial(r * w^2)/partialz = – partialp/partialz

+ 1/Re * { partialtauxz/partialx + partialtauyz/partialy + partialtauzz/partialz}

Energy: partialEt/partialt + partial(u * Et)/partialx + partial(v * Et)/partialy + partial(w * Et)/partialz = – partial(r * u)/partialx – partial(r * v)/partialy – partial(r * w)/partialz

– 1/(Re*Pr) * { partialqx/partialx + partialqy/partialy + partialqz/partialz}

+ 1/Re * {partial(u * tauxx + v * tauxy + w * tauxz)/partialx + partial(u * tauxy + v * tauyy + w * tauyz)/partialy + partial(u * tauxz + v * tauyz + w * tauzz)/partialz}

where Re is the Reynolds number which is a similarity parameter that is the ratio of the scaling of the inertia of the flow to the viscous forces in the flow. The q variables are the heat flux components and Pr is the Prandtl number which is a similarity parameter that is the ratio of the viscous stresses to the thermal stresses. The tau variables are components of the stress tensor. A tensor is generated when you multiply two vectors in a certain way. Our velocity vector has three components; the stress tensor has nine components. Each component of the stress tensor is itself a second derivative of the velocity components.

The terms on the left hand side of the momentum equations are called the convection terms of the equations. Convection is a physical process that occurs in a flow of gas in which some property is transported by the ordered motion of the flow. The terms on the right hand side of the momentum equations that are multiplied by the inverse Reynolds number are called the diffusion termsDiffusion is a physical process that occurs in a flow of gas in which some property is transported by the random motion of the molecules of the gas. Diffusion is related to the stress tensor and to the viscosity of the gas. Turbulence, and the generation of boundary layers, are the result of diffusion in the flow. The Euler equations contain only the convection terms of the Navier-Stokes equations and can not, therefore, model boundary layers. There is a special simplification of the Navier-Stokes equations that describe boundary layer flows.

Notice that all of the dependent variables appear in each equation. To solve a flow problem, you have to solve all five equations simultaneously; that is why we call this a coupled system of equations. There are actually some other equation that are required to solve this system. We only show five equations for six unknowns. An equation of state relates the pressure, temperature, and density of the gas. And we need to specify all of the terms of the stress tensor. In CFD the stress tensor terms are often approximated by a turbulence model.



One response to “FIELDS: CHAPTER 8”

  1. This may take 3 minutes to read but apparently takes 100 years to understand. .