




























































































Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
ment of Aerospace and Mechanical Engineering of the University of Notre ... Introduction to Physical Gas Dynamics, John Wiley, New York,.
Typology: Exams
1 / 298
This page cannot be seen from the preview
Don't miss anything!





























































































see Panton, Chapters 1-6, see Yih, Chapters 1-3, Appendix 1-2, see Aris.
We seek here to present an approach to fluid mechanics founded on the principles of rational continuum mechanics. There are many paths to understanding fluid mechanics, and good arguments can be made for each. A typical first undergraduate class will combine a mix of basic equations, coupled with strong physical motivations, and allows the student to develop a knowledge which is of great practical value and driven strongly by intuition. Such an approach works well within the confines of the intuition we develop in everyday life. It often fails when the engineer moves in to unfamiliar territory. For example, lack of fundamental understanding of high Mach number flows led to many aircraft and rocket failures in the 1950’s. In such cases, a return to the formalism of a careful theory, one which clearly exposes the strengths and weaknesses of all assumptions, is invaluable in both understanding the true fluid physics, and applying that knowledge to engineering design. Probably the most formal of approaches is that of the school of thought advocated most clearly by Truesdell,^1 sometimes known as Rational Continuum Mechanics. Truesdell de- veloped a broadly based theory which encompassed all materials which could be regarded as continua, including solids, liquids, and gases, in the limit when averaging volumes were sufficiently large so that the micro- and nanoscopic structure of these materials was unimpor- tant. For fluids (both liquid and gas), such length scales are often on the order of microns, while for solids, it may be somewhat smaller, depending on the type of crystalline structure. The difficulty of the Truesdellian approach is that it is burdened with a difficult notation
(^1) Clifford Ambrose Truesdell, III, 1919-2000, American continuum mechanician and natural philosopher. Taught at Indiana and Johns Hopkins Universities.
11
There are many subsets of mechanics, e.g. statistical mechanics, quantum mechanics, continuum mechanics, fluid mechanics, solid mechanics. Auto mechanics, while a legitimate topic for study, does not generally fall into the class of mechanics we consider here, though the intersection of the two sets is not the empty set.
Early mechanicians, such as Newton, dealt primarily with point mass and finite collections of particles. In one sense this is because such systems are the easiest to study, and it makes more sense to grasp the simple before the complex. External motivation was also present in the 18th century, which had a martial need to understand the motion of cannonballs and a theological need to understand the motion of planets. The discipline which considers systems of this type is often referred to as classical mechanics. Mathematically, such systems are generally characterized by a finite number of ordinary differential equations, and the properties of each particle (e.g. position, velocity) are taken to be functions of time only. Continuum mechanics, generally attributed to Euler,^6 considers instead an infinite num- ber of particles. In continuum mechanics every physical property (e.g. velocity, density, pressure) is taken to be functions of both time and space. There is an infinitesimal property variation from point to point in space. While variations are generally continuous, finite num- bers of surfaces of discontinuous property variation are allowed. This models, for example, the contact between one continuous body and another. Point discontinuities are not allowed, however. Finite valued material properties are required. Mathematically, such systems are characterized by a finite number of partial differential equations in which the properties of the continuum material are functions of both space and time. It is possible to show that a partial differential equation can be thought of as an infinite number of ordinary differen- tial equations, so this is consistent with our model of a continuum as an infinite number of particles.
The modifier “rational” was first applied by Truesdell to continuum mechanics to distin- guish the formal approach advocated by his school, from less formal, though mainly not irrational, approaches to continuum mechanics. Rational continuum mechanics is developed in a manner similar to that which Euclid^7 used for his geometry: formal definitions, axioms, and theorems, all accompanied by careful language and proofs. This course will generally
(^6) Leonhard Euler, 1707-1783, Swiss-born mathematician and physicist who served in the court of Cather- ine I of Russia in St. Petersburg, regarded by many as one of the greatest mechanicians. (^7) Euclid, Greek geometer of profound influence who taught in Alexandria, Egypt, during the reign of Ptolemy I Soter, who ruled 323-283 BC.
follow the less formal, albeit still rigorous, approach of Panton’s text, including the adoption of much of Panton’s notation.
The following are useful notions from Newtonian continuum mechanics. Here we use New- tonian to distinguish our mechanics from Einsteinian or relativistic mechanics.
x′^ = x − u (^) o t, (1.1) y ′^ = y − vo t, (1.2) z ′^ = z − w (^) o t, (1.3) t′^ = t. (1.4)
ρ = lim V → 0
i=1 mi V
Here V is the volume of the space considered, N is the number of particles contained within the volume, and mi is the mass of the ith particle. We can define a length scale
10 -8^10 -6^10 -4^ 0.01 1
x (m)
1
100
ρ (kg/m^3 )
variation on the continuum scale
variation on the sub-continuum molecular scale
Figure 1.1: Sketch of possible density variation of a gas near atmospheric pressure.
Example 1. Find the variation of mean free path with density for air.
We turn to Vincenti and Kruger for numerical parameter values, which are seen to be M =
λ =
(
) ( 10001 kmole mole^ ) √ 2 π
(
) ρ (3. 7 × 10 −^10 m) 2
, (1.7)
= 7.^8895 ×^10
− (^8) molecule mkg 2 ρ.^ (1.8) Note that the unit molecule is not really a dimension, but really is literally a “unit,” which may well be thought of as dimensionless. Thus, we can safely say
λ = 7.^8895 ×^10
− (^8) mkg 2 ρ.^ (1.9) A plot of the variation of mean free path λ as a function of ρ is given in Fig. 1.2. Vincenti and Kruger go on to consider an atmosphere with density of ρ = 1. 288 kg/m 3. For this density
λ = 7.^8895 ×^10
− (^8) mkg 2
, (1.10)
= 6. 125 × 10 −^8 m, (1.11) = 6. 125 × 10 −^2 μm. (1.12) Vincenti and Kruger also show the mean molecular speed under these conditions is roughly c = 500 m/s, so the mean time between collisions, τ , is
τ ∼ λ c =^6.^125 ×^10
− (^8) m 500 m s^ = 1.^225 ×^10
− (^10) s. (1.13)
0.01 0. 1 1 10
10 -
10 -
10 -
10 -
10 -80.
λ (m)
ρ (kg/m 3 )
ρatm ∼ 1.288 kg/m^3
λatm ~ 6.125 x 10-8^ m
Figure 1.2: Mean free path length, λ, as a function of density, ρ, for air.
Density is an example of a scalar property. We shall have more to say later about scalars. For now we say that a scalar property associates a single number with each point in time and space. We can think of this by writing the usual notation ρ(x, y, z, t), which indicates ρ has functional variation with position and time.
associates three scalars u, v, w with each point in space and time. We will see that a vector can be characterized as a scalar associated with a particular direction in space. Here we use a boldfaced notation for a vector. This is known as Gibbs^9 notation. We will soon study an alternate notation, developed by Einstein, and known as Cartesian^10 index notation.
Here we adopt a convention for the Gibbs notation, which we will find at times conflicts with other conventions, in which italics font (a) indicates a scalar, bold font (a) indicates a vector, upper case sans serif (A) indicates a second order tensor, over-lined upper case sans
seri (A) indicates a third order tensor, double over-lined upper case sans serif (A) indicates a fourth order tensor. In Cartesian index notation, their is no need to use anything except italics, as all terms are thought of as scalar components of a more expansive structure, with the structure indicated by the presence of subscripts. The essence of the Cartesian index notation is as follows. We can represent a three dimensional vector a as a linear combination of scalars and orthonormal basis vectors:
a = ax i + ay j + az k. (1.16)
We choose now to associate the subscript 1 with the x direction, the subscript 2 with the y direction, and the subscript 3 with the z direction. Further, we replace the orthonormal basis vectors i, j, and k, by e 1 , e 2 , and e 3. Then the vector a is represented by
a = a 1 e 1 + a 2 e 2 + a 3 e 3 =
i=
ai ei = ai ei = ai =
a 1 a 2 a 3
Following Einstein, we have adopted the convention that a summation is understood to exist when two indices, known as dummy indices, are repeated, and have further left the explicit representation of basis vectors out of our final version of the notation. We have also included a representation of a as a 3 × 1 column vector. We adopt the standard that all vectors can be thought of as column vectors. Often in matrix operations, we will need row vectors. They will be formed by taking the transpose, indicated by a superscript T , of a column vector. In the interest of clarity, full consistency with notions from matrix algebra, as well as transparent translation to the conventions of necessarily meticulous (as well as popular) software tools such as Matlab, we will scrupulously use the transpose notation. This comes at the expense of a more cluttered set of equations at times. We also note that most authors do not explicitly use the transpose notation, but its use is implicit.
1.2.1.2 Rotation of axes
The Cartesian index notation is developed to be valid under transformations from one Carte- sian coordinate system to another Cartesian coordinate system. It is not applicable to either general orthogonal systems (such as cylindrical or spherical) or non-orthogonal systems. It is straightforward, but tedious, to develop a more general system to handle generalized co- ordinate transformations, and Einstein did just that as well. For our purposes however, the simpler Cartesian index notation will suffice. We will consider a coordinate transformation which is a simple rotation of axes. This transformation preserves all angles; hence, right angles in the original Cartesian system will be right angles in the rotated, but still Cartesian system. We will require, ultimately, that whatever theory we develop must generate results in which physically relevant quantities such
β β β
α
α
α
P
Figure 1.3: Sketch of coordinate transformation which is a rotation of axes
as temperature, pressure, density, and velocity magnitude, are independent of the particular set of coordinates with which we choose to describe the system. To motivate this, let us consider a two-dimensional rotation from an unprimed system to a primed system. So, we seek a transformation which maps (x 1 , x 2 ) T^ → (x′ 1 , x′ 2 ) T^. We will rotate the unprimed system counterclockwise through an angle α to achieve the primed system. The rotation is sketched in Figure 1.3. Note that it is easy to show that the angle β = π/ 2 − α. Here a point P is identified by a particular set of coordinates (x∗ 1 , x∗ 2 ). One of the keys to all of continuum mechanics is realizing that while the location (or velocity, or stress, ...) of P may be represented differently in various coordinate systems, ultimately it must represent the same physical reality. Straightforward geometry shows the following relation between the primed and unprimed coordinate systems for x′ 1
x∗
′ 1 =^ x
∗ 1 cos^ α^ +^ x
∗ 2 cos^ β.^ (1.18)
More generally, we can say for an arbitrary point that
x′ 1 = x 1 cos α + x 2 cos β. (1.19)
We adopt the following notation