^^WA-ITvi'B^^ it. Ha. NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS TECHNICAL MEMORANDUM 1398 FORMATION OF A VORTEX AT THE EDGE OF A PLATE By Leo Anton Translation of "Ausbildung eines Wirbels an der Kante einer Platte," Gb'ttinger Dissertation, Ingenieur-Archiv, vol. X, 1939. UNIVERSITY OF FLORIDA DOCUMENTS DEPARTMENT 120 MARSTON SCIENCE LIBRARY RO. BOX 11 7011 GAINESVILLE, FL 3261 1 -701 1 USA Washington March 1956 l^i, 1i J ^v^ 3?^V0^{. NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS TECHNICAL MEMORANDUM I598 FORMATION OF A VORTEX AT THE EDGE OF A PLATE^ By Leo Anton 1. INTRODUCTION If a plate is moved from the state of rest transversely to its width through a fluid at rest, a small vortex forms at every edge at the start of the motion; 'by continual added infliix of new parts of the fluid, that vortex becomes larger and grows into the flow. In case the flow is perpendicular to the plate, the position of the two vortices is symmetrical at the beginning of the motion of the fluid. In a later stage of development it is, as von Karman has proved, unstable. Then new vortices originate alternately at the upper and. lower end of the plate which group themselves behind the plate in a certain manner (von Karman 's vortex street). In what follows, we shall treat the initial flow just described for a plate perpendicular to the approaching flow, where, therefore, flow direction and plate form a right angle . Since in this case the vortices are in the positions of reflected images, our investigation is limited to the formation and the growth of one vortex at one plate edge. We visualize the plate as suddenly set in motion and then moved uniformly at constant velocity. For a system of axes moved with the plate, the velocity at infinity is therefore to remain invariable throughout the flow duration. Moreover, we assume the plate to be laterally extended to infinite length so that we deal with a plane nonstationary flow phenomenon. Solution of the problem is considerably facilitated by disregard of the friction. This is permissible since the friction of most fluids coming into question is very small. Thus we take a^ a basis an ideal fluid where, therefore, the viscosity is zero, and require furthermore that it be homogeneous and incompressible. In this fluid the flow described is then to be interpreted as a potential flow with free vor- tices distributed over a discontinuity surface, a so-called vortex sheet. The velocity field of such a flow is uniquely determined by the velocity at infinity and by the form and vortex distribution of this surface. A solution for the motion of a fluid about a plate, starting from the state of rest, has, so far, been achieved only for infinitely small angles of attack. This borderline case has been treated by H. Wagner Ausbildung eines Wirbels an der Kante einer Platte," Gottlnger Dissertation, Ingenieur-Archiv, vol. X, 1939^ PP- 411-^4-27. MCA TM 1398 (ref. 1). Here the free vortices lie in a plane vortex sheet which extends from the trailing edge. The solution of the flow problem consists in the calculation of the distribution density of the discontinuity sur- face. Wagner performed these calculations for accelerated and for uniform motion of the plate. For finite angles of attack, as in our case, the discontinuity surface forms a strongly curved vortex sheet which obviously is rolled up into a spiral vortex. Thus, besides the distribution density, the shape of the discontinuity surface also is unknown. We now make the following simplification. At the very beginning of the motion of the fluid, the magnitude of a vortex at the plate edge in proportion to the width of the plate is still very small. Consequently the vortices at both edges still lie outside of their mutual interference domain so that one deals practically with a flow about a plate of semi- infinite width. This means a great advantage. Since no unit length exists which would influence the process, the flow conditions must be independent of the magnitude of the vortex at the time . The growth of the vortex therefore consists in a similar increase of it. In this connection the treatises of L. Prandtl (ref. 2) and H. Kaden (ref. 5) must be men- tioned. Prandtl performed a general investigation concerning the flow around a corner for certain laws of acceleration. No solution for the plate flow resulted from the Prandtl formulation; however, one can derive from it the laws of similitude. In the treatise mentioned, Prandtl sur- mises that the rolled-up vortex sheet in its interior may be of the type of a spiral R = const. /cp"^ (R and cp signify polar coordinates counted from the center of the spiral, m is a number). Kaden treats the rolling up of an unstable discontinuity surface, unilaterally of infinite length, into a spiral vortex. As we shall see later, Kaden 's problem becomes identical with ours for the interior of the vortex. The same laws of similitude are valid for both cases. The solution for the vortex core / 2/5 R = const. /cp ' found by Kaden applies, therefore, also to the interior of our vortex sheet; the above-mentioned prediction of Prandtl is thereby confirmed. The outer part of the discontinuity surface, thus the transitional region from the spiral-shaped vortex core to the plate edge, we procure with the aid of graphical methods. The solution obtained is then used for finding the flow about the plate of finite width. 2. OUTLINE OF THE METHOD In this section we shall briefly describe the method by which one arrives at solution of the problem. We shall start from the plate of finite width, and shall then treat the plate of infinite width as a special case. NACA TM 1398 5 In the case of the vortex-free potential flow about the plate placed transversely, the velocity at the edges is known to become infinitely large; however, under actual conditions this is never possible. We must therefore demand also in this case that the velocity there always has finite values . We attain this by assuming in the fluid so-called Helmholtz discontinuity surfaces which extend from the plate edges. By this expression one understands vortex sheets consisting of free vortices and obeying Helmholtz' laws. In contrast to this is the plate, which also is a discontinuity surface but no longer obeys Helmholtz' laws. It forms a rigid surface and exerts, therefore, pressures on the flow. The plane in which the flow takes place will henceforth be designated as complex z-plane. The system of coordinates is selected so that the origin of the coordinates coincides with the center of the plate and the plate lies along an imaginary axis. The width of the plate is made equal to 2b . The movement of the plate is to take place from the left to the right, or, which means the same, the flow is to approach the plate from the positive side. The constant velocity at infinity - which, with oppo- site sign, may also be interpreted as speed of the plate traveling in a fluid at rest - is to be -v^ (fig- 1^ left). We conformally map the z-plane onto a ^-plane, by means of the transformation / 2 + b2 (1) with the plate, which represents, of course, a piece of the imaginary axis in the z-plane, becoming a piece of the real axis in the new plane (fig. 1, right). The edges z = +ib shift into the zero point ^ = 0; the point at z = transforms to the points ±b on the real ^-axis. Furthermore, in this transformation the region at infinity in the z-plane goes over, without change, into that at infinity in the ^-plane. There- fore the velocities at infinity in both planes are equal, thus v(z = Qo) = v(^ = 00) = -v^. As images of the discontinuity surfaces we obtain in the new plane again images which start from the point ^ = 0. The new vortex sheets lie symmetrically to the real axis whereby the latter becomes a streamline so that the image of the plate also remains a streamline-'-. The flow in the ^-plane corresponding to the plate flow is composed of the parallel flow and a flow caused by the free vortices of the discontinuity surfaces. We are now going to calculate its velocity field. For the case where the vortex sheets in the z-plane no longer lie symmetrical to the real axis, one must choose instead of (l) the transformation into a circle since the real axis in the ^-plane then is no longer a streamline. For our case the selected transformation is preferable to that into the circle because the graphical calculations become considerably simpler. MCA TM 1398 We denote by ^-p a variable point of the vortex sheet in the upper half plane and coordinate to it a circulation dP. The corresponding vortex point of the mirrored surface is then given by ^-p, and has the same vortex strength, but with opposite sign, thus -df. Both vortex points produce at an arbitrary point ^ of the plane a velocity the conjugate complex value of vhich is given by dv^(0 = dT 2.1^^ -^F ^_g If one puts dr= 7(^p)d^p wherein /(^p) represents the distribu- tion density at a point ^p, and d^„ a line element of the surface, one obtains by integration over all points of the surfaces the velocity induced by them (^p' signifies the temiinal point of the vortex sheet.) The parallel flow has the velocity -v^. For the total velocity of the superimposed flow one may now write =F (2) Between the velocities of both planes there exists the relationship H^ = HU^ (3) dz from which one obtains the velocities for the z-plane. The condition that the velocity should assume finite values at the plate edges remains to be satisfied. As one can see very quickly, this requirement is ful- filled when the expression (2) for C = disappears, because the magni- fication ratio d^/dz becomes infinite for • ^ = or z = tib . One obtains therewith the condition^ The part stemining from V]_(^) becomes real for ^ = 0. NACA TM 1398 which must always be satisfied. It is^ however, not sufficient for deter- mination of distribution density and shape of the discontinuity surface. As long as the vortices are still very smaH compared to the plate, the plate may be regarded as infinitely wide in comparison with the vor- tices. For this case, the lack of a fixed comparative length is the reason that the form of the discontinuity surface is independent of the magnitude of the vortex configuration. The development of the discon- tinuity surface consists, therefore, initially in a similar magnification. This similar magnification can be fulfilled only in the case of a certain shape of the discontinuity surface, and of a certain distribution of the vortices over it. On the basis of this condition, we can therefore cal- culate the form and vortex distribution for the beginning of the motion. Starting from this initial condition we can then find, by further obser- vation of the variation with time of the change in form, 'also the shapes and vortex distributions for the later times when the vortices are no longer small compared to the width of the plate . Before performing the limiting process to the plate of infinite width, we choose a new system of coordinates z' which arises from the former coordinate system by parallel displacement, and the zero point of which is shifted to the upper edge (z = +ib). We write z = z' + ib (5) and replace, in addition, z by z' = 5z or Z = ^ z' (6) b H where H represents a scale unit. Accordingly, equation (l) assumes the form ^ = l/fSz')^ + 2iHZ (7) The limiting process from the plate of finite width to that of infinite width (b — > 00) yields a new traas format ion I t^ ^ = \J2±EL or Z = - i |jj (8) 6 NACA TM 1598 This is the transformation of the plane with a slit along the negative imaginary axis onto the upper half plane of the ^-plane . If in the ^-plane v^^ continues denoting the free-stream velocity, the interpre- tation of v^ as traveling or, respectively, free-stream velocity, is lost in the Z-plane since in this plane the velocity at infinity tends ^0. The quantity v^^ appears in this plane as the velocity which prevails in the case of vortex- free flow at the point Z = iH. 2 All other considerations carried out so far concerning the plate of finite width may he directly transferred to the plate of infinite width. 5. THE INFINITELY WIDE PLATE (a) Laws of Similitude By the limiting process from a finitely wide to an infinitely wide plate which was just completed, a flow type was obtained for which the successive flow patterns are similar and, accordingly, the circulation of the vortex and the distribution density on the vortex sheet are similarly enlarged for this growth. We shall briefly derive here the laws which characterize this behavior, (As was remarked in the Introduction, the laws of similitude could be immediately derived from the quoted treatise of Prandtl.) We observe the distance Z of two points remaining in similar position during the similar enlargement of the discontinuity surface and assume it to be at the time t-, of the amount Z-|_ and at the time to of the amount Z2. We state as the time law for the enlargement Z2 \t2J where A is a number yet to be determined. The law for the circulation r about similarly situated regions may be obtained as follows. Since the free-stream velocity in the ^-plane is invariably Vqo, the same veloc- ity must prevail also for the similar enlargement at similarly located pointG5. This velocity is composed of the free-stream velocity v^o and of the fields of the individual vortex elements. The influence of a In the Z-plane this is not the case, since here (for vortex-free ii: 2 during the similar enlargement, thus does not remain in similar position. flow) the velocity is constant at the point 'Z =-— iH which is fixed NACA TM 1398 7 vortex element AF^ at the distance Zf from the point considered^ on the velocity is Ar/2nlj-. If during the gro-wth AF and Zf are modified, AF/Zf, due to the constant velocity, also must be constant. For the circulation F about similarly situated regions of the ^-plane, therefore, F2 z^ 2 is valid. In the Z-plane the corresponding distances axe ^1 _ hi"^' 2/ Z2 yZ^ Thus — becomes in the Z-plane r2 ^,t © If we designate by V the velocities which correspond to the law of similitude, we obtain for the velocity V(z) in the Z-plane which is pro- portional to r/z Vi_^iAi_ 1^ /tiN-^/2 Fg/Zg \j Z-]_ \t^i (11) If we replace in the ratios above the time t-|_ by the time 1 and the time tp by the time t, the laws of similitude read z = z^t^ F= r^t^s V = Vit"^/2 (12) MCA TM 139« wherein the number A remains to be determined. For calculation of the number A we form the velocity of growth of the distance I dZ , , . A-1 = 6 T At dt -^ Since this velocity, like all velocities, must vary proportionally -A/2 to t~ ' , there results A-l = -^ thus A^2 4 Therewith the laws of similitude read z = Zit2/5 r = rit^/5 V = Vit-1/5 (15) (b) Solution for the Vortex Core Our task now consists in finding a shape for the spiral which satisfies these laws. They are the same laws Kaden (ref. 5) obtained as a result. Thus one deals here, too, with a similar problem. In particular, both prob- lems become identical for the inner part of the vortex where in both cases a discontinuity surface consisting of free vortices is to be rolled up into a spiral vortex. It is therefore sufficient if we refer to Kaden 's paper and here only briefly mention the results. Let R and cp be polar coordinates of a system of coordinates the zero point of which coincides with the center of the spiral. For form and circulation F of the vortex core bounded by a circle of the radius R, Kaden obtained ^=m^/3, r = 2xv/S (11+) MCA TM 1598 The circumferential velocity^ for radius R is ^^ = ^ = ^ (15) Therein t is the time, counted from the beginning of the motion. The quantity X is a constant, still unknown for the time being. It has the dimension velocity times square root of a length. Thus we may also intro- duce for it the dimensionless constant k = -^ (16) (c) Solution for the Outer Loops of the Spiral As was shown before, the behavior of the spiral is characterized in /xt\2/3 the interior, that is, for large angles 9, by R = — . For small angles cp, in contrast, especially where the discontinuity surface adjoins the edge, considerable deviations from this form occur. This transitional region can be found as follows. We visualize the plate as lying parallel to the straight line cp = 0; this assumption is insignifi- cant for practical purposes since the spiral windings approach, in the inward direction, a circular form. We isolate, furthermore, at one point the inner part of the spiral, in which the form is prescribed with suffi- /xt^2/5 cient accuracy by R = — I . For calculation of the velocity field Vjtqpy outside of the spiral core, we visualize the latter as replaced by an isolated vortex at the center of gravity of the circulation. This is directly permissible since almost circular symmetry prevails in the interior of the spiral. The dissymmetry caused by the separation point brings it about that the center of gravity of the circulation does not coincide with the center of the spiral. For the calculations indicated later on, the separation point is placed at the point cp = 2.5fl:« The difference between the coordinates of the spiral center (a, h) and of the center of gravity (a;i_, h^) may be calculated, similarly as by Kaden, (ref. 3) to be approximately 10 NACA TM 1598 a^ - a = 0.006 R^ and h^ - h = O.O5I R^ (17) if R. is the radius of the core of the spiral. The outer loop which we now want to calculate, extends from the separation point (cp = 2.5rt) to the plate edge (cp = cp*). Its shape and vortex distribution 7 as well as the position and circulation F^ of the vortex core we assume, at first, arbitrarily, with a factor of simili- tude X common to both quantities and at present not yet determined, (equations (ik) and (I5)), and make therewith the transition into the ^-plane. Therein 7 is transformed into 7' = 7 ^ whereas the circulation rj_ remains unchanged. In this mapping plane (^-plane) the velocities stemming from the vortex may be calculated as the field of the vortex core and the distribution on the outer winding, likewise their mirrored images. By superimposltion of the undisturbed velocity v^ there results the velocity field v' in the ^-plane. First, one ascertains the velocity at the point ^ = 0. From the condi- tion that this velocity must be zero (equation (k)) results factor X, still undetermined at first, for the circulation f^ and the distribu- tion density 7. In the conformal mapping onto the Z-plane the velocities v' of the ^-plane are transformed into the velocities v of the Z-plane, according to the relationship - -, d^ V = V — dZ where v and v signify the conjugate values of the velocities. The velocity of the individual vortex concentrated at the core is obtained by omission of the field of this vortex. It must, however, be noted that in the transition from the ^-plane to the Z-plane the field of this vortex is deformed; hence an additional term appears in the conversion of the velocities. In the ^-plane the potential of the core vortex to be omitted is 1 NACA TM 1598 11 *k(?) —!-(?- to) In the Z-plane the potential of the vortex to be omitted is 23ti V ^0 / 2fli V'' ^0 2:ti V In order to calculate the velocity for the Z-plane we must^ therefore, subtract, besides the potential ^y-i^) ^ also the additional potential r. — — In K + t^ \ to which corresponds a velocity 2jti V 'J/ v' additional in the ^ -plane or, respectively. ^i 2iti U + C "1 rlr' 2n± K + U dZ additional in the Z-plane . Once one has calculated the velocities for individual points of the discontinuity surface one resolves them, most advantageously, into normal and tangential velocities v^^ and v^, and plots the latter as functions of the arc length s . By combination of the function values one obtains the velocities pertaining to each point of the vortex sheet; from them results the motion of the assiomed vortex sheet. A fluid particle at the point P of the spiral having the velocity V (fig. 2) moves by the distance PP-i = v in the unit time. Due to the 12 NACA TM 1598 similar enlargement, the point P is tranformed during this time into the similarly located point P2« The velocity at which the similarly situated points are displaced is assumed to be V so that the distance is PP2 = V. To have P-^ come to lie on the similarly enlarged spiral, the normal components V and Vj^ of the two velocities v and V must be equal or AVn = v^ - Vn = (19) We observe the point P on its path to the similarly situated point Pp and find that this path does not represent the motion of a fluid particle. Rather, the fluid flows, with the velocity v - V, through the point moving from P to P2. If V1J1 and Vrp signify the tangential components of the velocities V and V, the circulation in flowing with the fluid into the core inside of P or P2^ respectively, during the time At is Af = (vt - V^)7At Hence there results for the region inside of P, which is being similarly enlarged — = (vrp - V^)/ ot According to the laws of similitude (13)^ one must have for a region being similarly enlarged ^r ^ r 1 t-2/5 = _L St ^ 5 5t The two values of dP/dt must agree in order to have the similitude satisfied, therefore the equation MCA TM 1598 15 A = (Vrn dt ^ Vt)7 3t = (20) must apply. Since the distribution density is \ we may write instead of (20) also ht OS (Vrp - Vrp)7 ^ = 3t (21) With the aid of the equations (19) and (2l) it is now possible to find the shape and vortex distribution of the outer loop. We shall perform the calculation of these quantities for a time t t = « h5/2 X (20)' for which the equation (ik) of the spiral core assumes the simplified form R H ,2/5 (21)' 9' We select as the initial form for the entire spiral the form R = H/cp^/^ and place it as a first try so that It joins the plate at the point cp = cp-x- = jt/2 (fig. 5). The break originated at the edge is merely a flow which^ as one readily understands, disappears at the next moment. Except for the factor x^ the circulation 2XJR ^ (22) The duplication of the equation numbers 20 and 21 follows that of the original Germain docujnent. lU NACA TM 1598 is known to us. For the core which we shall assume to begin at cp = 2.5it, we have R = R = 1_^ (2.5rt) Thus r. becomes r, = 2x,. / 5 (25) V(2.5«)2/5 The distribution density is 7 - Sr _ ar dR - V R . , X ^21,) Bs Sr Ss \/9Ah5 + r5 \/icp2/3\|9/i^ + cp- In figure 4 this distribution density 7 is plotted as a function of the arc length s. The conforraal mapping (8), in combination with the reflection at the |-axis, leads to the double spiral indicated in figure 5> The distribution density in this plane (image plane) is represented in figure 6. When we calculate the velocity at the zero pointy we obtain, on the basis of equation {h^ y 8- condition for the still undetermined value X. There results <-L = 2.36\/Hv„ or k-[_ = 2.56 As an example for the determination of the velocity field we shall here briefly reproduce the calculation of the velocity at the center of gravity. In the ^-plane an element of the discontinuity surface ds induces at the center of gravity the velocity components . /(s ' ) , , /(s ' ) , dvt = sin \|/ ds , dv^, = cos ^ ds ^1 2jtRs' '1 2jtRs' MCA TM 1398 15 when Rg ' denotes the distance of the element from the center of gravity and ijr the angle of the radius vector from the center of gravity to the element with the ^-axis (fig. 5)« By integration over s' from the zero point to the core boundary (cp = 2.5jt, s' = 2.62 H) and over the corre- sponding curve (s') of the image one obtains vt ' and "v^ ' . In fig- ure 7 the values dv^ '/ds' and dv^ '/ds' are plotted against s' or s '■ respectively. The quantities ve ' and v^ ' are obtained by cir- cumscribing; with the aid of a planimeter^ the cross-hatched areas bounded by the curves and the abscissa axis These values represent the influence of the outer loop of the spiral. The image of the vortex core causes the velocities Ti = 0.22v_ v„ = ^2 2«2h' ^ ~ ^2 if h' signifies the distance of the center of gravity of the core from the |-axis. The vortex core itself does not contribute to the velocity of its own center of gravity in the ^ -plane. However^ it must be noted that, according to the explanations on pp. 10 and 11, in the conversion to the Z-plane a term stemming from the vortex core vt^' = 0.11 v_ Vr,, ' = 0.11 v„ ^13 - u.xi v„ v^^ appears , For the resultant velocity one obtains ^1 ' = - ^co + "^li + v^2 ^ "^^3 " -O.i^Ov^ Vti = v^^ + v^^ + v^^ = 0.11v„ l6 NACA TM 1598 whence there results, by conversion to the Z-plane, Vx = H —^ f- = - 0.29v<„ r + ^ ^y = H-^- -^= - 0.17v„ r + r\ According to 'the Law of similitude, Vx = - 0-57v<„ Vy = would have to be valid. The center of gravity thus shows for the assumptions made, relatively to its required motion, a wrong motion obliquely downward, toward the plate, with the velocity components ^^x = ^x = V^ = 0-08v^ ^V = ^y - Vy = - 0.l7v„ In the same manner the velocities for the points of the outer winding are calculated. Of course, one has to add a term which represents the influ- ence of the core; on the other hand, the additional term V:^ does not appear here. We find the calculated normal and tangential components plotted in figures 8a and 8c. As can be seen from figures 8b and 8d, equations (I9) and (21) are not satisfied. In other words, the shape and vortex distribution assumed deviate from the actual solution. We now modify the position of the core e.g. in the sense of Av^^ and Av„ and the shape and vortex distribution of the outer loop in the sense of the differences Av„ and A -^. The fonn of the spiral core remains "^ dt 'Xt\ 2/5 ^ H unchangedly R = (-i_i = . We now repeat the calculation procedure H cp5/2 NACA TM 1598 17 carried out so far. After several steps, the first two of which are represented in figures 9 "to 12, one arrives finally at a solution that satisfies the laws of similitude; it is plotted in figure 13 • In figure li<- is the radius of the spiral; in figure I5 its distri- bution density is plotted as functions of the center angle and of the arc length s, respectively. For more convenient further use, the most important constants necessary for characterization of the spiral vortex have heen briefly compiled below. For the solution indicated in figures 15 to I5, the calculation yields the time t = 1.42 i the constant X = 2.22H^/^ k = 2.22 00 the coordinates of the spiral center a = -0.53H h = -0.22H the total circulation of the spiral vortex r = lj-.20Hv„ For the assiomed core boundary at cp = 2.5rt the radius of the core is R^ = O.267H the circulation of the core Ti = 2.3Hv„ 18 NACA TM 1598 For an arbitrary time t one obtains a = -0.i+Ht2/3 h = -0.l65Ht2/5 R = r = 3.7Hv<„t = ^A. TT+2/5 1/3 a cp2/5 cp2/3 Ri = 0.2Ht^/5 = 0.5a r. = 2Hv^t^/5 ^ 5.i5v/I i4-. THE PLATE OF FINITE WIDTH With the solution found for the plate of infinite width, we know the initial flow conditions about a plate of finite width, where the vortices are still small compared to the plate width. We know that for this initial state the successive vortex images originate from one another by similar enlargement. However, as soon as the vortices assume a magnitude which is no longer negligible compared to the plate width, the presupposition for a similar growth no longer holds true. Yet, starting from the form and distribution density we fovmd for the still very small vortices, we can calculate the variation of this form and distribution density, using the values of Av^^ [(19)^ and A — (2l) which now, with growing vortices, more and more deviate from 0. We only must make use of the transformation by means of (l) instead of the trans- formation by means of (8). If we start from a known form and distribu- tion density at a time t-i , we obtain the deviation of the form and dis- tribution density from those which would result for similar enlargement according to the time laws (13) , at the time to as being t2 n tg by An = / AVj^dt and Aj = A — ^ dt ^^1 -^^1 bt I NACA TM 1398 19 The graphical determination of the required quantities is greatly hampered by the fact that the vortices are at the beginning very small and later very large. We can avoid this difficulty by visualizing the figure^ at every instant, enlarged or reduced in such a manner that the vortex would always remain the same if it would keep on growing according to the laws of similitude that are valid for the beginning. The figure in actual size we denote as figure I and the quantities valid for it by the subscript I, the enlarged figure as figure II and the corresponding quantities by the subscript II. At the initial stage, the enlarged vor- tex is to agree precisely with the one calculated for the plate of infinite width in the preceding section. Let us call the ratio of mag- nification in every case e. If we reduce in the magnified figure all velocities simultaneously at suitable points in the proportion l/v/e", the circulations F are magnified in the proportion \fe. The distribu- tion densities 7 are reduced in the proportion l/u^- For the potential flow about the plate of the width 2b in the plane I there results in the neighborhood of the edge the velocity TI = -pj=l= (25) ^^ For the plate of Infinite length treated in section 3 "the corresponding velocity was (26) wherein y' = b - y signifies the distance from the edge of the plate. In order to obtain for the beginning of the motion the same conditions as for the flow about the plate of infinite width treated before, we select the width of the plate 2b in such a manner that b = H For the above stipiilation, exactly the velocity V;„ prevails in the magnified figure, while e is still very large, at a point at a distance of b/2 from the edge. In the case of the plate of infinite length treated before, this point lay at a distance H/2 from the edge. Since we equated b = H, we obtain in the magnified figure precisely the flow treated in section 3. 20 NACA TM I598 According to equation (1^4-)^ for the growing vortex core of the plate of infinite width Since at the initial stage the flow about the edge of the plate of the finite width 2b = 2H coincides with the flow about the edge of the plate of infinite width^ this formula applies also to the initial stage in the case of a plate of finite width. By means of the magnification by the factor e, this initial vortex is to be transformed into the flow TT >, treated for which R = — ^j- = j—. Thus necessarily q)2/3 cp2/3' Re = e Xt\ 2/5 «^/ cp2/3 Hence there results . = (1)2/5 M-2/5 (^7) Since in the plane I^ at the beginning of the motion, the circulations 1/5 Fj grow, according to (l3), with t ' , we obtain in the plane II, where they appear magnified in the proportion ye, constant circulations and constant distribution densities. At the beginning of the motion e = <». The plate edges lie in the plane II at an infinite distance from one another so that we actually have the case of the plate of infinite width. In time, however, e attains smaller finite values, and we obtain in the plane II also a plate of finite width. Therewith the velocities become different, and we obtain deviations from the similar magnification. In order to calculate the velocities in the plane II, we transform this Zjj-plane onto a ^-plane which we shall characterize by the subscript III, by means of the function te = lP^- (28) NACA TM 1398 21 For very large e this mapping is transformed into the infinite plate according to (8) in the neighborhood of the edge. For finite e, however, other shapes result for the mapped vortices, and that is the very reason which causes the modifications of the velocities and there- with the deviations from similar vortex growth. As soon as, due to this deviation, the form and distribution density of the vortices have changed with respect to the similar growth, this also contributes to the variation of the velocities. Since, however, form and distribution density change, at first, only very slowly compared to the similar growth, one may in the plane II assume the form and distri- bution density as constant in turn through a large time interval, and need consider in this time interval only the modification of the trans- formation in the plane III, due to the modification of the value e. For such a time interval t-^ to t2 (the first starts with t = 0), one calculates for several intermediate times the normal and tangential components Vn^- and v-^ of the velocities of the plane II with the aid of the transformation onto the plane III in the same manner as in the case of the plate of infinite length. One forms furthermore the differ- ences with respect to the velocities of the similar magnification Vr and V+ ^11 'nil and obtains then by transfer to the plane I the values (l9) Av ^I t-'-ii ' ^f-ii = €Av - V, n 11/ (29) and corresponding to (21) 3t Ss-j- (-Ti - Vtj)7i 71 3t = -e Ssj I [ ^11 ^13 )/^-/ 3t 5/2 If we, finally, replace in the last term t by (b/e) ' Jt/x (equation (27)), we obtain B7 . A— 1 = -e' Bt =11 ■II ■^Ilj ^11 ^ X7ii 3jtb 3/2 (30) 22 NACA TM 1598 By graphical integration then results the displacement An of the points of the vortex sheet at right angles to it, and the modification of the distribution density compared to the similar magnification in the plane I as Z\ni ti ^ Av^^dt A7, = Jti A — i dt (31) (52) and the modification of form and distribution density in the plane II as Anjii = egAnj (53) ^''11 = ^ % (54) Therein £2 is the value of the magnification ratio e at the time t2. Due to the finite magnitude of the vortices compared with the plate ^7ti width, Av„ ^ and A f 0: also, the condition (k) , ''ll at that in the ^ -plane at the zero point the velocity must be zero, will no longer be satisfied. In the ascertainment of the vortex for the plate of infinite width we have been able to fulfill this condition by suitable definition of an as yet undetermined factor for the circulation. Due to this condition we found the quantity X or k, respectively. For the further development of the vortex, fonii and vortex strength and their variation with time are fixed. Only the strength of the vortices shed at the plate edge is still undetermined since we can here not form the differential quotient S/Ss occurring in (50). We must select it in such a manner that equation (h) is satisfied. We obtain therefore an additional modification of the distribution density starting from the plate edge which gradually is carried into the vortex by the flow. It NACA TM 1398 25 is true that it was shown in the quantitative calculation that the devia- tions from the condition (k) are extremely small, because the variation of the conformal representation [(l)] with e results in a positive velocity, whereas the modification of the form and distribution density according to (33) and. (3^) results in a negative velocity at the point ^ = and the two almost cancel one another. The calculation was carried out for the intervals e = 00 to e = 3> e = 3 to 6=2, aind e = 2 to e = 1. The results are compiled in the figures 16 to I9. True to expectation, the circulation F increases more slowly with time than it does according to the solution for the plate of infinite width (initial state, fig. 19) • In the final state it would perhaps approach a constant value which corresponds to a steady state of flow. However, according to experience the symmetrical vortex configuration becomes unstable from a certain magnitude onward, so that this steady state is not attained. 5. SIMVIARY The flow about the plate of infinite width may be represented as a potential flow with discontinuity surfaces which extend from the plate edges. For prescribed form and vortex distribution of the discontinuity surfaces, the velocity field may be calculated by means of a conformal representation. One condition is that the velocity at the plate edges must be finite. However, it is not sirfficient for determination of the form and vortex distribution of the surface. However, on the basis of a similitude requirement one succeeds in finding a solution of this problem for the plate of infinite width which is correct for the very beginning of the motion of the fluid. Starting from this solution, the further development of the vortex distribution and shape of the surface are observed in the case of a plate of finite width. Finally, I should like to express my special gratitude to Professor Betz for his suggestion of this investigation and his active support in carrying it out. Translated by Mary L. Mahler National Advisory Committee for Aeronautics 2k NACA TM 1398 REFERENCES 1. Wagner, H. : Z. angew. Math. Mech. 5, 1925, p. l?- 2. "Prandtl, L. : Uber die Entstehung von Wirbeln in der idealen Fliissigkeit, mit Anwendung auf die Tragfliigeltheorie und andere Aufgaben. Vortrage aus dem Gebiete der Hydro- und Aerodynamik, Innsbruck 1922, p. I8. 3. Kaden, H. : Ing.-Arch. 2, I93I, p. l40. NACA TM 1398 25 y O ^» z- Plane -V, 00 ^ ^00 U C- Plane -V, 00 -v^ ^ '00 Figure 1.- Plate with vortex sheets starting from the edges (left), and conEormal representation of the flow (right). Figure 2.- Actual velocity v of a fluid particle and displacement velo- city V corresponding to the similar magnification. Z=X+iY (Z=0) Figure 3.- Spiral according to the equation R = H/q)2/3_ 26 MCA TM 1598 0.5 r(s) K hI/2 0.4 \ N S^ 03 V. \ \ ^ 02 ,^ ^ - Ol / n-J l\ S 't \ I H 0.5 1.5 2.0 2.5 Figure 4.- Distribution density for the region tt/S = cp = 2, Sir pertaining to the spiral R = H/cp 2/3 NACA TM 1598 27 -e ' Figure 5.- The spiral and its image in the ^plane. 28 NACA TM 1398 i /(S'] '' K -H^ /2 04 / \ s 0.3 / \ s / \ 0.2 / \ / \ 01 / / s / H 0.5 1.0 1.5 2.0 2.5 Figure 6,- Distribution density 7(2') of the spiral in the ^-plane. NACA TM 1598 29 A 27r dvf, H 3/2 1 (^ ZTT - ds' ^3/2 K 1JD 1 0.5 s' H' i s' H Trrrn Tlllll TTTiiK H 1 1 <: > 1 > 1 f f \ : s 0.5 ' \ IF -10 Figure 7.- Influence of the distribution on the velocity components V. and V of the vortex center of gravity. 50 MCA TM 1598 vn 0.6 0.5 04 0.3 0.2 0,« -0.1 -0.2 -0.3 ) / \ / -\ Vn/ 1 1 n i\ J/ \ s y \ TT 1 M \ \ 3 \ -Vn '00 09 r 01 i c n / H 04 N / > n 9 'V J 0.9 08 Vqc ^^Jx_ ^1 \ 0.7 0.6 / / \^ /' \ 0.5 0,4 0.3 02 / r ^T 0.1 / \ / ' s H ni 1 \ ' _n 9 \ \ / -0.3 -04 \ J Figure 8.- The quantities which are decisive for the variation of form and distribution density for the initial spiral. NACA TM 1598 51 -X Figure 9.- First correction of the form. y(s) 1 4 Voo 0.9 0.8 0.7 0.6 0.5 04 N, ^ \ \ \ \ V ^ ^ 0.3 02 0.1 s H 1 02 0, 5 1.0 „ 1.5 2.0 2.5 f Voo i - 02 > at voo- 0.1 -0.1 -0.2 igure 10,- ] S '^* \ A s / 1- -U.l Lo \ /\ H \ 1 1 <2 ^ \ ^^ t r Dii 3tri b but] Lon der -0.2 LSity aj nd c ( iec] -siv e q la.n titles correction. 52 NACA TM 1598 Figure 11.- yis) 1.1 1.0 0.9 0.8 0.7 0,6 0.5 04 03 0.2 0.1 0" Second correction of the form. Voo S \, s \, V \ \ ^>^ ^"^^ \ s N s s H 0.5 0.2 0.1 -0.< -0.2 T^4 1.0 1.5 a 0.2 2.0 2.5 ■1- 0.1 -0.1 -0.2 A dyl dt H s N. ^ H \ i / I V / Figure 12.- Distribution density and decisive quantities after the second correction. MCA TM 1398 53 Figure 13.- Final form of the spiral for the plate of infinite width. i 0.8 0.7 06 0.5 04 Q3 02 0.1 R V A / /T* W^Outer loop \ tY' U 1 piral R=-4^ ^^ >::s >-a y,^/^ " ■■ ^==^ — - y * 1
\,
/
^
^^
^>
€=oo\
\"
\
^
\^
S
N^s
V
1
^
■
^
^
s.
s
b
0.5
1.0 1.5 2.0 2.5
Figiire 18.- Variation of the distribution density with time.
'b^oo
1 1 1 1 ' 1 1 1 I
Plate of infinite width — s^._
..
,.^
-j:SL
'J—
—
—
"""
^
-^
-^
P
late of finite widtf
1
/
t
'
b/v
03
0.5
1.0
1.5
Figure 19.- Growth of the circulation with time.
NACA - Langley Field. Va.
:^„-bS
b5«
< n! O
!>. o ■
s <= S •"
^ S "" 2
curias
;? ^ s ><
2 -a "" o
nl g xi >
D. O (U c
S <^ "' a
^ B J3 2
S * " o
rt m ? X!
1 rt
_ hD
j?
■a
0) o
T5
c
>
n.
a- -3
^ to
Cli ^
_o
a
>>^
o
o " J2 _ 5
1^- O ••^ "O -*-*
o''a ° I!
5 -S .a O
.§ 3 -g C
"O * .2
rt " £-S
° "> ■? n,
to T3 *J ■*-*
m '
cs
a -B o
" 0)
? QQ
f •a =a :§ ■
O rt O 0)
■S 13 S
to 0, ti
_ _ c ■" ^
"c •*-* i
oJ S a J:;
o c 0) o
tH .5 O O
3 3 b (U
■2-3^-^
1 ^
S o
(3 u
SJ 05 O — 1
00 Q) S ' t-
erod
enta
low,
o-g^ M-Hr..
< B li<
h, < Z £ > o.
(J
<
z
*J ^
< 0)
oiO
-go
SQ
C Ed
■"<:
a M
So
o>
^<=
°o 2 o
" t- Z
'-' T3 O
O X! K
< n) O
3
S bO
- to ia
H lo o
01 '-'""
01
S " £
3h^
I
. ^-*
bl, = 2 o.
^ o S oT
2 0) Ui CO
£1 <: "-■
3 O -
13^>
01 .S
< 2
a
o
5 >'
o -> <: 03
01 (U g '
(-1 3 O
° c < bo'-' ■
b. <: 2 « >
<
CJ
<
z
0) J
00
CO
<
c
m a,
t-. ■" ^
I ° o
^ M a
■^ T3 nl
c" =5 o
.2 C «,
t! 2 M
01 S 0)
-" T3 .-.
O 0) *
ix ^ I-
SP > m
c o
M o rt
o
u
<
2
CO
CO
<
U
<
z
0) 0) .
S " £
«; 3 ^
> « a,
i a* ^
"S I" ^
^ w a.
-S-
c " o
° c o
3 X! rt
•3 3 o
"> S a
£-2 c
■" tj ■-
o