2 Chapter 2: Differential Equations of Motion
The need to understand the local detail of the flow drives the need for differential forms of the fluid flow equations. These differential forms are based on taking the limit of our previous conservation equations in the limit as the CV approaches an infinitesimal CV (or differential element). In general, these equations are derived from the evaluation of the conservation equations on an infinitesimal differential element level using vector calculus. More specifically, the approach gathers all the terms into a single volume integral that must apply everywhere, hence, provides a governing differential equation.
Starting with the mass conservation equation, Eq (1) in Chapter 1, we initiate by moving the time-differential term into the volume intetegral. Also, applying divergence or Guass’ thoerom (which analytically relates volume and surface integrals), the second (flux) term can be expressed as a volume integral. The results from this procedure yields the following
(1)
The volume integral terms can then be combined to provide the following application of CV applied to a differential element:
(2)
This equation must be satisfied throughout space and time and for any possible CV, which can only be achieved if the integrand also equals 0, hence, the conservation of mass in differential form given as
(3)
The first term ()can be described as the local time-rate-of-change of mass. The second term () can be described as changes associated with mass moving into and away from the local point (regarded as advection).
A second approach to derive the conservation of mass in differential form can be developed by utilizing a Taylor-series expansion about the center of the infinitesimal element as indicated in Fig. 1. The mass flux into the left side can be given as . The mass flux can then be extrapolated using the 1st term of the Taylor series to the right face as . The approach extends to all surfaces and yields the same governing equation, i.e. Eq. 3. See this video for the procedure:
The next aspect to consider is the limit as the flow is incompressible flow, i.e., is constant. This condition leads to a low that is Divergence Free. This literally implies that the conservation of volume and mass are the same. One can think of it as squishing in the x-direction is equally absorbed in the y and z directions. In this limit, conservation of mass is given as:
(4)
This simplified form drives analyses relevant to flow in liquids and low-speed applications, i.e., . Alternatively, high-speed flows and thermal flows that have variation in demand the conservative form given by Eq \label{eq:consMassDE2}.
(5)
In the condition when the flow is irrotational () and incompressible (), there is a connection with conservation of momentum through Bernoulli equation, Eq. 3 in Chapter 1, which links the fluid dynamics to provide pressures required for aerodynamic forces. Through these equations, combined with a variety of foundational, elementary flows, drive many aerodynamic theory and numerical predictions referred to as potential flow models. A nice tool and further discussion can be explored at potentialflow.com.
Now move to the momentum conservation equation,
(6)
Similar to before, the time-differential term is moved within the volume integral. And again, the second (flux) term is expressed as a volume integral using the Divergence Theorem. Focusing on the LHS of Eq (6), the procedure implies:
(9)
This equation must be satisfied throughout space and time and for any possible CV, which can only be achieved is the integrated also equals 0. This yields the conservation of mass in differential form given as:
(10)
The solution of this set of equations, in combination with conservation of mass, is very difficult due to instabilities in the equations that lead to turbulence (see box below). A subset of these equations are the Euler equations, which are simply the form when .
Turbulence, by definition, are fluctuations in the flow with a turbulence cascade. This cascade involves larger structures, breaking into smaller and smaller structures that eventually lead to very high viscous mechanisms and losses. A physical description of turbulence is described here:
The characterization and modeling of these high frequency turbulent fluctuations drive the Reynolds Averaged Navier Stokes Equations (RANS) which are briefly described in this video:
Lastly, similar approaches are used to obtain the differential form of the energy equation. The result in terms of a low-speed flow is given as
(12)
In general, the solutions to these equation sets drive fluid behavior and modeling associated with aerodynamics. These model equations serve to describe pressure, temperatures, and velocities throughout a fluid flow under a give set of both initial conditions (ICs) and boundary conditions (BCs). In steady aerodynamics we typically focus on BCs, which for a viscous calculation are:
- Free-stream velocity
- Free-steam pressure
- No slip condition on solid (aerodynamic) surfaces
In Potential Flow (inviscid) models, our BCs are typically
- Free-stream velocity
- Free-steam pressure
- Flow tangency or no penetration into solid (aerodynamic) surfaces
This flow tangency is a major consideration in the context of potential flow that drives much of the foundation of aerodynamics. And, in general, the goal of aerodynamics is to maximize lift while reducing drag. Therefore potential flow assumptions are often near valid and drive much of the most useful aerodynamic theory.
Example 1: Potential Flow Modeling of Airfoil (Boundary Element Method)
This is an example of a potential flow, boundary element method, applied to a cylinder.
Example 2: Potential Flow Model of Airfoil (XFOIL)
This is an example of the usage of Xfoil to analyze an airfoil
Example 3: CFD Model of Airfoil
This is an example of the application of StarCCM+ to an airfoil with a plain flap.
Key Skills: Consider these as you work through homework.
- Know how to derive these equations using both calculation and differential methods.
- Understand the basic concepts of derived models (potential flow and RANS) and what the methods mean and where they apply and do not apply.