Equation of Motion in Streamline Coordinates

1 Equation of Motion in Streamline Coordinates Ain A. Sonin, MIT, 2004 Updated by Thomas Ober and Gareth McKinley, Oct. 2010 2.25 Advanced Fluid Mechanics Euler’s equation expresses the relationship between the velocity and the pressure fields in inviscid flow. Written in terms of streamline coordinates, this equation gives information about not only about the pressure-velocity relationship along a streamline (Bernoulli’s equation), but also about how these quantities are related as one moves in the direction transverse to the streamlines. The transverse relationship is often overlooked in textbooks, but is every bit as important for understanding many important flow phenomena, a good example being how lift is generated on wings. A streamline is a line drawn at a given instant in time so that its tangent is at every point in the direction of the local fluid velocity (Fig. 1). Streamlines indicate local flow direction, not speed, which usually varies along a streamline. In steady flow the streamline pattern remains fixed with time; in unsteady flow the streamline pattern may change from instant to instant. Fig. 1: Streamline coordinates In what follows, we simplify the exposition by considering only steady, inviscid flows with a conservative body forces (of which gravity is an example). A conservative force per unit mass G is one that may be expressed as the gradient of a time-invariant scalar function, G = -\U(r) , and the steady-state Euler equation reduces to (1) 2 � � 1 V · VV = - Vp - VU(r) . P (2) A uniform gravitational force per unit mass g pointing in the negative z direction is represented by the potential U = gz . (3) A streamline coordinate system is not chosen arbitrarily, but follows from the velocity field (which, we note, is not known à priori). Associated uniquely with any point r and time t in a flow field are (Fig. 2): the streamline that passes through the point (streamlines cannot cross), the streamline’s local radius of curvature R and center of curvature, and the following triad of orthogonal unit vectors: i1s : in the flow direction : in the normal direction, away from the local center of curvature !in y y y il : in the bi-normal direction, ( il = is x in ). The unit vectors define incremental distance ds measured along the streamline in the flow direction, dn measured in the normal direction, away from the center of curvature, and dl measured in the bi-normal direction. The radius of curvature R is defined as positive if in points away from the center of curvature, and negative if in points toward it. The unit vectors, the radius of curvature, and the center of curvature all change from point to point and in unsteady flows from time to time, depending on the velocity field. To transform Euler’s equation into streamline coordinates, we note that in those coordinates1, ! d ! d !d V = is + in +i ds dn dl (4) V = isV (5) and where V is the magnitude of the velocity vector V . From (4) and (5), ! d V ·V =V ds and thus 1 The gradient of a scalar function f (s,n,l ) from this definition and the expression is defined by Vf · dr � df (s, n, l ) = dr = is ds + in dn + il dl (6) !f !f !f ds + dn + dl . Equation !s !n !l (4) follows for an incremental displacement in streamline coordinates. 3 (V )V = V s (V is ) = is V2 i +V2 s . s 2 s (7) The unit vector in the last term of (7) changes orientation as one moves along the streamline. The change dis in is from s to s+ds is obtained with the construction shown in Fig. 2 as ! ! ! ds dis = -in d8 = -in R (8) Fig. 2: Incremental change in the streamwise unit vector from s to s+ds. from which we see that d is i =_ n ds R (9) Using (9) in (7), we obtain the convective acceleration as (V · V)V = is d (V 2 J V2 - in ds 2 R (10) The first term on the right is the convective acceleration in the direction of the velocity, and the second is the centripetal acceleration, toward the center of curvature. The pressure gradient in streamline coordinates is p = is p p p + in + il s n l (11) 4 Using (10) and (11) in (2), we obtain the equation of motion in streamline coordinates for steady, inviscid flow as s-direction: n-direction: d (V 2 l 1 dp dU =ds 2 P d s ds (12) V2 1 dp dU =R p dn dn (13) 1 dp dU P dl dl (14) - 0=- l-direction: In a uniform gravitational field U=gz and these equations read s-direction: ( s 1 2 ) V2 = 1 p s n-direction: V2 1 dp dz =-g R P dn dn l-direction: 0=- g 1 dp dz -g P dl dl z s (15) (16) (17) For constant-density flow in a uniform gravitational field, the equations simplify further to s-direction: d LV 2 ( p + Lgz + J=0 ds 2 (18) n-direction: d pV 2 ( p + pgz) = R dn (19) l-direction: l (p+ gz) = 0 (20) The s-direction equation (18) states Bernoulli’s theorem: the total pressure⎯the sum p + pgz + pV 2 2 of the static, gravitational, and dynamic pressures⎯remains invariant along a streamline. The n-direction equation (19) states that when there is flow and the streamlines curve, the sum p + pgz (which is constant in when the fluid is static) increases in the ndirection, that is, as one moves away from the local center of curvature. 5 The l-direction equation (20) states that p + pgz remains constant for small steps in the binormal direction, that is, the pressure distribution is quasi-hydrostatic distribution in the l-direction. EXAMPLE Consider the simple case of 2D, inviscid air flow over a smooth hill (Fig. 3). Far upstream of the hill the incident velocity is uniform at V= . The hill deflects the air around it, and a uniform flow is again established far downstream. Far upstream, above, and downstream of the hill, the pressure is constant at p= and the streamlines are straight (the hill does not perturb the flow at “infinity”). We shall assume that gravitational effects are negligible (the medium is air and the hill’s elevation is modest) and the free stream’s Mach number is small, so that and the density can be taken as constant. Based on the available equations, what can we say about the pressure and velocity distributions over the hill—where is the velocity higher than V= , for example, and where lower? Fig. 3: Sketch of streamlines in a 2D flow over a hill. To answer this question accurately we need to know the shapes of the streamlines throughout the flow field—or, at least, in the region that is perturbed by the hill. We don’t have this information, so we proceed by drawing a rough estimate of the streamline pattern, as shown in Fig. 3. The difference between the pressure at infinity and at the top of the hill, point (3), can be estimated by integrating equation (19) along the vertical path from (3) to ( 0 ). Since this path follows the local n-direction, R>0 everywhere along it. Neglecting the gravitational term, (19) gives dp pV 2 = dn R from which we see that (21) 6 p3 = p 3 V 2 dn >0 R (22) Thus p3 < p= , and according to Bernoulli’s equation (18), it follows that V3 > V= . Using similar arguments, we conclude that p1 = p= and V1 = V= , and p2 > p= and V2 < V= , etc. In principle, if R(n) and V(n) can be established or estimated, the integral in (21) can be evaluated. For example if we find that the flow perturbation caused by the hill is negligible at elevations greater than some multiple of the height of the hilltop, we might write for the path from (3) to ( 0 ) n R (23) Rhill e H where Rhill is the streamlines’ radius of curvature in the vicinity of the hilltop, n is measured from the top of the hill upward, and H = {h is some multiple β of the actual height h of the hilltop, the coefficient β being an empirical number. From Bernoulli’s equation (18) we also have that p+ V2 V2 , =p + 2 2 (24) Substituting for R and V into (21) from (23) and (24), respectively, we integrate (23) and obtain V 2 � R hill p3 = e 2 2H p � 1 (25) For a low hill such that 2H<