Inertial Instability of Arbitrarily Meandering Currents Governed by the Eccentrically Cyclogeostrophic Equation

Size: px
Start display at page:

Download "Inertial Instability of Arbitrarily Meandering Currents Governed by the Eccentrically Cyclogeostrophic Equation"

Transcription

1 Jornal of Oceanography, Vol. 59, pp. 163 to 17, 3 Inertial Instability of Arbitrarily Meandering Crrents Governed by the Eccentrically Cyclogeostrophic Eqation HIDEO KAWAI* Shibagahara, Kse, Joyo, Kyoto Pref. 61-1, Japan (Received 11 Jne 1; in revised form 9 Jly ; accepted 5 Agst ) Using natral coordinates, we have derived a criterion for the inertial instability of arbitrarily meandering crrents. Sch crrents, governed by the eccentrically cyclogeostrophic eqation, are adopted as the basic crrent field for the parcel method. We assme that any virtal displacement which is given to a water parcel moving in the basic field has no inflence on this field. From the conservation of mechanical energy for a virtal displacement we derive an inertial instability freqency ω m = [( f + /r)z].5 for the eccentrically cyclogeostrophic crrent, where f is the Coriolis parameter, the velocity (always positive), r the radis of crvatre of a streamline (negative for an anticyclonic meander), and Z the vertical component of absolte vorticity. If ω m is negative, the eccentrically cyclogeostrophic crrent becomes nstable. Althogh the conventional, centrifgal instability criterion, derived from the conservation of anglar momentm in a circlarly symmetric crrent field, has a certain meaning for a monopolar vortex, it contains a radial shear vorticity that is difficlt to se in arbitrarily meandering crrents. The new criterion ω m contains a lateral shear vorticity that is applicable to arbitrarily meandering crrents. Examining instabilities of concentric rings with radii of 5 1 km, we consider reasons why the anticyclonic spersolid rotation has been very mch less freqently observed than the cyclonic spersolid rotation, despite a prediction of some common stability and a rapid change in radial velocity gradient for the former. Classifying eccentric streamlines into the large and small crvatre-gradient types, we point ot that the largegradient crvatre in anticyclonic rings is apt to be nstable. Keywords: Inertial instability, eccentrically cyclogeostrophic crrent, natral coordinates, parcel method, virtal deflection, instability freqency, ring, sbsolid rotation, spersolid rotation, zero potential vorticity. 1. Introdction While the inertial instability of circlarly symmetric flows has been treated since Rayleigh (1916), that of arbitrarily meandering crrents, having no symmetry with respect to any point or any line, has received little attention, as far as the athor is aware. To derive the instability freqency for sch arbitrarily meandering crrents, we apply the parcel method to the crrents that are governed by the eccentrically cyclogeostrophic eqation, expressed in terms of natral coordinates. In the parcel method we assme that any virtal displacement of a water parcel has no inflence on the basic crrent field arond the parcel. The basic crrent field means a steady crrent field that is to be examined to see whether it is stable or not with respect to the virtal displacement. Applying the parcel method to the basic crrent field * address: kawaihid@d.dion.ne.jp Copyright The Oceanographic Society of Japan. of meandering crrents, governed by the eccentrically cyclogeostrophic eqation, we find a relation between a virtal displacement of a water parcel across streamlines of the basic field and a virtal deflection-angle of the corse of the parcel, made to the streamline. From the conservation of mechanical energy for the virtal displacement we derive an inertial instability freqency for the eccentrically cyclogeostrophic crrent in terms of natral coordinates. When this instability freqency is imaginary, we can regard the basic crrent field as nstable. Using the inertial instability criterion we examine velocity strctre in anticyclonic and cyclonic rings with radii of 5 1 km.. Natral Coordinates We se two right-hand coordinate systems. One is the sal coordinate system (x, y, z): positive x and y are eastward and northward, respectively, on a horizontal plane; positive z is pward along a vertical line. The other is a natral coordinate system (ξ,, z): positive ξ is down- 163

2 stream along a streamline; along an orthogonal to streamlines; z the same as the above. ξ and are not crvilinear bt metric coordinates on a horizontal plane. The relation between nit vectors (i, j) along the (x, y)-axes and (i ξ, i ) along the (ξ, )-axes becomes i ξ icosθ + jsinθ, i isinθ + jcosθ, i ξ, (1) where θ is an angle of i ξ made to i, is the horizontal crrent velocity vector, and (the scalar of ) becomes positive in any case. The horizontal gradient operator H and the directional differentiation along the natral coordinate axes ( / ξ, / ) are given by H i / x + j / y = i ξ / ξ + i /, H + k / z () / ξ i ξ = cosθ / x + sinθ / y / i = sinθ / x + cosθ / y, (3) where k is the nit vector along the z-axis. Definition (1), () and (3) easily lead to the following relations: = i / x + j / y + k / z = i ξ / ξ + i / + k / z i ξ = (i / x + j / y + k / z) (icosθ + jsinθ) = i ξ θ/ z + k θ/ ξ i = i θ/ z + k θ/ i ξ = (i / x + j / y + k / z) (icosθ + jsinθ) = θ/ i = θ/ ξ = (i ξ ) = i ξ i ξ = i ξ θ/ z + i / z + k( θ/ ξ / ) (4) = (i ξ ) = i ξ + i ξ = / ξ + θ/. The radis of crvatre of a streamline is defined by 1/r θ/ ξ, (5) where r is measred from a point on the streamline toward the center of crvatre of the streamline (Fig. 1(b)). This measring direction is opposite to the measring direction of R in polar coordinates (Fig. 1(a)). Since r defined by (5) becomes positive (negative) for a cyclonically (anticyclonically) meandering streamline, the direction of r coincides with the positive (negative) -axis. Fig. 1. Difference in kinematic notations between horizontal polar coordinates (a) and natral coordinates (b). + signs in front of R, ξ and indicate the positive sense of each coordinate axis. Here and hereafter the cyclonic or anticyclonic deflection means that in the Northern Hemisphere. By the definition (3), the directional differentiation of sm, prodct and qotient of two or more fnctions, composite fnction, etc., are obtained nder the same rle as in the sal differentiation, even if θ is dependent on space coordinates and time. When considering the directional differentiation of the second order, however, the difference between the sal directional differentiation and the extended one becomes distinct (Kawai, 1957). 3. Eqations of Motion of an Unsteady State Althogh we regard the basic crrent field for the parcel method as steady, we need eqations that govern the motion of a virtally displaced parcel. This section is devoted to expressing eqations of motion of an nsteady state in terms of the natral coordinates. Assming hydrostatic eqilibrim of the basic crrent field, we can redce p/ρ at z in a top or pper layer with niform density to p/ρ = g (ς z), (6) where p is the presser, ρ the density, g the acceleration of gravity, and ς an elevation of the sea-srface of the basic field above a reference level. The eqation of motion for horizontal inviscid flow becomes t H Hp + + kz + =, ( 7) ρ where Z is the vertical component of absolte vorticity. From (4) and (5), 164 H. Kawai

3 Z = f + k = f +/r /, (8) where f is the Coriolis parameter. From (1), = ( iξ) t t = i ξ + i t θ. 9 t ( ) From (9), (8) and (6), scalar prodcts of (7) with i ξ and i become d 1 p d + g = + ς = ( 1 ρ ξ ξ ) dθ p d f + f g + 1 = + θ + ς ρ =, ( 11 ) respectively, where the material differentiation is defined by d/ / t + / ξ. (1) 4. Acceleration de to Virtal Displacement We regard the state of a field of the eccentrically cyclogeostrophic crrent as basic (Section 1). Neglecting / t and θ/ t contained in (1) and (11), respectively, we get eqations that govern the eccentrically cyclogeostrophic crrent, ξ b b f + r + gς b b = ( 13 ) g + ς =, ( 14) where the sbscript b indicates the basic state. Adopting the parcel method, we conventionally assme that any virtal displacement of a water parcel has no inflence on the basic crrent field arond the parcel. Using (1), (13), (11), (5) and (14), we derive the ξ- and -components of the eqation of motion for the parcel as follows: d b = g ς b = ξ ξ ( 15) θ b g f fb f = ς t r = + r. r ( 16 ) b Fig.. Geometry of virtal displacement. P is the initial position occpied by a water parcel. P d is the final position after a virtal displacement v = θ v ξ, where ξ = t means a displacement of the parcel along a basic streamline ξ for a minte time interval t. ξ and ξ d are basic streamlines passing throgh P and P d, respectively. θ v is a virtal deflection-angle made to ξ at P. C and C d are crvatre centers of streamlines, ξ at P and ξ d at P d, respectively. r and r d are radii of crvatre of streamlines, ξ at P and ξ d at P d, respectively. Althogh (16) does not contain any time derivative of v, from (19) below we shall find the term θ/ t in (16) eqivalent to d v /. Let a parcel, initially located at P, move with a crrent velocity along a basic streamline ξ that passes throgh P (Fig. ). Add to the moving parcel a virtal, minte, cross-stream displacement v for a minte time interval t. Denoting a virtal deflectionangle of the corse of the parcel, made to ξ, by θ v, we express v in terms of θ v as follows: v = θ v ξ = θ v t, (17) where ξ = t is a displacement of the parcel along the streamline ξ for t. From (17), d v / = θ v, (18) becase t is an independent variable, which can be sfficiently minte. From (18), 19 d v d θv d θ, = = ( ) becase θ v is an arbitrary variable independent of spatial coordinates. In the last term of (19), the sbscript has been removed from, in order that (19) may hold not only at P bt also at any point. Before sbstitting the Inertial Instability of Arbitrarily Meandering Crrents 165

4 last term of (19) by (16), sbscripts on the rightmost side of (16) shold be changed as below, becase the acceleration of (19) is indced by the virtal displacement, namely, by the vales of the qantities on the rightmost side of (16) at the displaced position P d. Sbstitte by c to denote a qantity conserved for a virtal displacement. Sbstitte the sbscript b by d to indicate qantities at P d. Sbstitte r by r d from the assmption of no inflence on the basic field by any virtal displacement of a parcel. Ths, sbstitting the last term of (19) by (16), which has been changed as mentioned above, we get The material differentiation of (5) becomes de = d g d d + ς g = + ς ξ =, ( 6) becase of relations, ς/ t = for the basic field and (1). This represents the conservation of mechanical energy of a parcel for the virtal displacement. From (6) with (5), ths c = + g (ς ς c ). (7) d v c θ d = f d c = ( )+ t r d. ( ) Expansion of ς c in powers of ξ and v abot ς, by sing (13) and (14), yield 5. Redcing Qantity Vales at Virtally Displaced Position This section is devoted to redcing vales of qantity at P d in () to those at P. 5.1 Qantities depending on basic crrent field at new position P d Vales of d, d and 1/r d at P d redce to their vales and gradients at P as follows: = + O + ξ ξ ξ, d v v + ( ) ( 1) d = + + v O v ξ ξ, + ( ) ( ) = + v O v 3 rd r r + ξ r + ( ξ, ). ( ) ξ Here O( ξ, v ) denotes a sm of negligible terms with the sqared ( ξ, v and v ξ) or frther higher (n ) powers of infinitesimal qantities ( ξ, v ) n, omitting some factors of finite qantities, for simplicity. 5. Qantities depending on property conserved by a parcel for virtal displacement Mechanical energy per nit mass of an incompressible water parcel is given by Sbstittion of (6) into (4) yields E = / + p/ρ + gz. (4) E = / + gς. (5) ς ς g( ς ςc)= g ξ g v O v ξ, = + ( ) ξ + f + v O ξ v 8,. r Sbstittion of (8) into (7) yields + ( ) ( ) c = + ξ + f + v O ξ, v, 9 r which leads to c ξ = 1+ + f + r v + ( ) ( ) 1 / = + ξ + f + v O ξ, v. r + O( ξ, ) + ( ) ( 3) 6. Inertial Instability Criterion Sbtracting (3) from (1), in consideration of (8), we get f ( d c ) = f Z v + O( ξ, v ). (31) Sbtracting (9) from (), in consideration of (8), we get d c = Z v + O( ξ, v ). (3) Mltiplying (3) by (3), we get ( d c )/r d = Z v /r + O( ξ, v ). (33) v 166 H. Kawai

5 Sbstitting (31) and (33) into (), and neglecting higher order terms, we redce () to d v = f + Z r v ωmv = ( 34) ω m = ( f + /r)z. (35) In (34) and (35), the sbscript has been removed from, r and Z, in order that they may hold not only at P bt also at any point. For concentric streamlines, we write r d = r v, which leads to v v = = + + O( v ). ( 36) r r r r r d 1 Taking the same procedres as mentioned before (34) except for sing (36) instead of (3), we obtain (34) and (35), too. Whether adopting the eccentrically cyclogeostrophic flow or the exactly cyclogeostrophic flow as the basic field, we apparently get the same form of ω m as was given by (35), bt expressions of Z are essentially different from each other, as will be shown between (41) and (4), below. If ω m >, (34) represents a simple harmonic motion with a freqency ω m. If all of the relative vorticities in (35) vanish, ω m redces to f, the freqency of the inertial oscillation. If ω m <, d v / increases with increasing v, which means the inertial instability. According to common meteorological textbooks, the inertial instability freqency for a zonal geostrophic wind ω z is given by ω z = f ( f / y) = f Z, (37) becase ( f / y) is eqivalent to Z for a zonal jet stream by (8). The sign of ω z given by (37) depends on the sign of Z only, bt the sign of ω m given by (35) depends on the sign of ( f + /r)z. A postscript to Rayleigh (1916) states, It may have been observed that according to what has been said above the stability of flid motion in cylindrical strata reqires only that the sqare of the circlation increase otwards. This, Rayleigh s stability criterion, is formlated as d(vr) /dr >, (38) according to Kloosterziel and van Heijst (1991, hereafter KH), where R is the radis, which is measred otward from the center of circlar streamlines to a point nder consideration, positive for both anticyclonic and cyclonic vortices, and V is the crrent velocity, which is positive (negative) for an anticyclonic (a cyclonic) vortex (Fig. 1(a)). Since R is always positive, division of (38) by R 3 leads to (V/R) (V/R + dv/dr) >. (39) Assming a circlarly symmetric basic state for the se of horizontal polar coordinates, several athors (e.g., KH; Carnevale et al., 1997; Potylitsin and Peltier, 1998; Orlandi and Carnevale, 1999) mentioned a criterion for the inertial instability in the f-plane, which sed to be called the centrifgal instability. Among them KH obtained a stability criterion d( fr / + VR) /dr >, (4) by applying a parcel method to the basic field of circlarly symmetric flows, governed by the exactly cyclogeostrophic eqation. Since R is positive, division of (4) by R 3 leads to ( f + V/R) ( f + V/R + dv/dr) >. (41) From (35) with (8), the stability criterion obtained in the present stdy becomes ω m = ( f + /r) ( f + /r / ) >. (4) This criterion (4) looks similar to that of (41) or (39), bt there is an essential difference between them, as will be mentioned in the next paragraph. In (39) and (41), R is always positive, directed otward in horizontal polar coordinates, bt V is positive (negative) for a cyclonic (an anticyclonic) vortex (Fig. 1(a)). By contrast, in (4), is always positive, bt r is positive (negative) for a cyclonically (an anticyclonically) meandering streamline, becase the center of crvatre of the streamline is on the left (right) of the streamline, facing downstream (Fig. 1(b)). Becase the difference between (4) and (41) is whether the basic crrent field is circlarly symmetric or not, this difference appears in the expression of shear vorticities sch as / in (4) and dv/dr in (41). Since is eqivalent to V, the essential difference in shear vorticity between (4) and (41) is the se of a variable in (4) instead of R in (41). Althogh the radial, shear vorticity dv/dr has a certain meaning for a monopolar vortex, it loses applicability to arbitrarily meandering crrents becase R is not an independent variable in the latter crrents. Adopting arbitrarily meandering crrents, governed by the eccentrically cyclogeostrophic eqation, as the basic field, and sing natral coordinates, we shall be able to apply this new criterion widely to varios crrent fields, obtained with advanced crrent measrement technology. Inertial Instability of Arbitrarily Meandering Crrents 167

6 7. Anglar-Velocity Gradient in Concentric Rings 7.1 KH stability criterion for concentric rings Becase r/ = 1 (Fig. 3) for both anticyclonic and cyclonic, concentric rings, (8) redces to Z = f + /r / = f + /r + / r. (43) Sbstittion of (43) into (35) yields ω m = ( f + /r)z = ( f + /r) ( f + /r + / r). (44) Since (44) is the same as the left side of (41) obtained by KH, it may be called the KH stability criterion for concentric rings. According to varios observations of anticyclonic and cyclonic rings by Tomosada (1975), Brown et al. (1983), Mied et al. (1983) and others, rings are more or less circlarly asymmetric, e.g., elliptical or D-shaped. Condcting laboratory experiments on barotropic vortices in a rotating flid, KH caght trianglar anticyclonic vortices on streakline photographs of small tracer particles. Despite the asymmetry of rings, Figs. 3 and 4, Tables 1 and and the analyses in Section 7 are based on schematically circlarized rings for simplicity. 7. Dependence of the sign of ω m on Z The soltion of (14), which is a qadratic eqation with respect to b, yields the azimthal crrent velocity fr 5. 4g ς = 1± 1+ f r r > for cyclones ; r < for anticyclones, 45 ( ) ( ) where the sbscript b in (14) has been omitted. In the same way as in (43), sbstittion of = r into the lateral gradient ς/ of (45) for both cyclonic and anticyclonic rings yields the radial gradient ς/ r. We sed to choose the + sign in front of the sqare root in (45), assming that the crrent decreases to zero when the horizontal pressre gradient vanishes. Althogh this assmption has not generally been accepted, it is valid for rings, becase any ring is characterized by a certain, vertical displacement of the main thermocline from that of srronding water, and hence by a remarkable radial pressre gradient. A vortex withot sch vertical displacement can no longer be classed as a ring. From (45), we can derive algebraically the traditional velocity limit for anticyclones, < cr, (46) where cr is the critical velocity defined by cr fr/. (47) Since vales of cr (Table 1) are mch higher than the vales of velocity measred in anticyclonic rings, (46) is also regarded as actally holding. From (46) with (47), f + /r > (48) for anticyclonic rings. For cyclonic rings r > and >, and hence (48) also holds. From Eqs. (48) and (44), the sign of ω m becomes the same as that of Z for both anticyclonic and cyclonic rings. Fig. 3. Radial profiles of velocity for the sbsolid (broken), solid (solid) and spersolid (chained) rotations in schematic, anticyclonic (a) and cyclonic (b) rings, which have positive and negative shear zones, divided by a maximm azimthal velocity m at r = ±r m (r < for anticyclonic rings). Fig. 4. An example of greater vales of / r of three types of rotations of cyclonic rings than those of anticyclonic rings. Typical vales of / r, which are indicated by the maximm slope of a tangent to velocity crves vs. r, are based on Fig. 5. They are measred by the nit of f = s 1, the Coriolis parameter at N. 168 H. Kawai

7 7.3 Stability of rotations in the main body of rings We call portions of a ring inside and otside of the velocity maximm ( m in Fig. 3) the main body and fringe, respectively. The main body ( / r > ) corresponds to the negative shear zone of anticyclonic rings and the positive shear zone of cyclonic rings. The fringe ( / r < ) corresponds to the positive shear zone of anticyclonic rings and the negative shear zone of cyclonic rings. In order to do the algebraic analyses leading to Fig. 5, we introdce three qantities W, X and Y, and express the radial gradient of anglar velocity by where From (48), (/r)/ r = ( / r /r)/r = W/r, (49) W X Y, X / r, Y /r. (5) Y >.5f. (51) Becase (49) with (5) characterizes radial profiles of azimthal velocity in the main body of rings, we classify rotations into three types as follows. (1) Sbsolid rotation is characterized by decreasing, absolte vales of anglar velocity with increasing r, namely, W/r <, which leads to / r <. Becase W > for anticyclonic rings (r < ), from (49) /r < / r <, namely, Y < X <. Hence, the lower limit of (43) becomes Z > f + /r = f + Y. Becase f + Y > from (51), Z becomes positive, which means that the sbsolid rotation is always stable. Becase W < for cyclonic rings (r > ), from (49) /r > / r >, namely, Y > X >. For both anticyclonic and cyclonic rings, /r > / r, namely, Y > X. () Solid rotation is characterized by a constant anglar velocity with increasing r, namely, W/r =, which leads to / r =. Becase W = for both anticyclonic and cyclonic rings, from (5) Y = X. Hence, (43) becomes Z = f + Y. Becase f + Y > from (51), Z becomes positive, which means that the solid rotation is always stable. (3) Spersolid rotation is characterized by increasing, absolte vales of anglar velocity with increasing r, namely, W/r >, which leads to / r >. Becase W < for anticyclonic rings (r < ), from (49) Table 1. Critical velocity cr = fr/ (cm s 1 ) for a radis of crvatre r (km) of anticyclonic rings. Table. Smmary of inertial instability of rings examined by ω m = ( f + /r) Z. Anticyclonic ring [r <, >, f + /r > ] Negative shear zone (Body) [ / r < ] sbsolid rotation [/r < / r < ] ω m > always stable solid rotation [/r = / r < ] ω m > always stable spersolid rotation [ / r < /r < ] If Z >, ω m > stable If Z =, ω m = netral If Z <, ω m < nstable Positive shear zone (Fringe) [ / r > ] ω m > always stable Cyclonic ring [r >, >, f + /r > f ] Positive shear zone (Body) [ / r > ] ω m > f always extremely stable Negative shear zone (Fringe) [ / r < ] If Z >, ω m > stable If Z =, ω m = netral If Z <, ω m < nstable ω m, KH stability freqency; Z = f + /r + / r, vertical component of absolte vorticity. Inertial Instability of Arbitrarily Meandering Crrents 169

8 Fig. 5. Geometrical analysis of ranges of parameters W (= / r /r) and Z (= f + /r + / r) for spersolid, solid and sbsolid rotations in the main body of anticyclonic and cyclonic rings, illstrated on the X (= / r) Y (= /r) coordinate plane. W indicates the extent to which the spersolid and sbsolid crves are warped (Fig. 4). Positive W- and Z-axes are directed to 45 clockwise from the positive X- and Y-axes, respectively. Each axis is gradated in.5f. / r < /r <, namely, X < Y <. The spersolid rotation becomes stable, netral or nstable, according as Z is positive, zero or negative. Becase W > for cyclonic rings (r > ), from (49) /r > / r >, namely, X > Y >. For both anticyclonic and cyclonic rings, / r > /r, namely, X > Y. 7.4 Stability of the fringe of rings (1) The positive shear zone of anticyclonic rings [ / r >, /r < ] always becomes stable, becase Z > as follows. From (43), Z = ( f + /r) + ( / r /r). Since ( f + /r) > from (48), and since ( / r /r) > from the ineqalities bracketed above, we get Z >. () The negative shear zone of cyclonic rings [ / r <, /r > ] becomes stable, netral or nstable, according as Z is positive, zero or negative. 7.5 Contradiction between observation and prediction by the stability criterion Figre 4 illstrates detailed radial profiles of velocity in the main body of rings. The spersolid crves in Fig. 4 appears to be warping otside the short arrows. Becase W defined by (5) is the difference in slope between / r (the slope of a tangent to the velocity crve) and /r; we can se W as an indicator of the extent to which the spersolid and sbsolid crves are warped. According to a geometrical analysis from Fig. 5, which was prepared from algebraic analyses among X, Y, Z and W (Sbsection 7.3), the pper limit of W for the stable (Z > ) spersolid rotation of anticyclonic rings is f. Althogh this vale is smaller than the pper limit of W for either spersolid or sbsolid rotation of cyclonic rings, it is greater than.5f for the sbsolid rotation of anticyclonic rings. Hence, common observations of anticyclonic rings in the spersolid rotation are predicted, as far as Z is positive. According to Kawai (199, sbsection 4.1), however, hardly any measrement has spported the continos existence of anticyclonic rings with a static core or in spersolid rotation. Most measrements have caght anticyclonic rings in solid or sbsolid rotation. Anticyclonic rings with a static core have scarcely been observed, nor have cyclonic rings withot a static core. These observational reslts of rings contradict the prediction of common observations of anticyclonic rings in spersolid rotation based on the above geometrical examination of W according to Fig. 5. This contradiction can be explained in terms of an observational limitation de to discrete measrements of velocity in connection with the difference in magnitde of W in velocity profiles. According to Fig. 5, the spersolid rotation of anticyclonic rings is stable when Z >, bt becomes nstable when Z <. By contrast, the spersolid rotation of cyclonic rings is extremely stable throghot the pper right qadrant, regardless of rotation types. Hence, W of the spersolid rotation of cyclonic rings occpies wider range than W of the stable, spersolid rotation of anticyclonic rings. Conseqently, even in case of missing velocity measrements arond the velocity maximm m, we have more opportnity to confirm the spersolid crve of cyclonic rings by obser- 17 H. Kawai

9 vations of remaining, remarkably warped crves. However, we have less opportnity to confirm the spersolid crve of anticyclonic rings by observations of remaining, less warped crves (Fig. 4). This is one way to interpret the contradiction between observation and prediction. 8. Crvatre Gradient in Eccentric Rings We can write / = ( / r) ( r/ ). (5) The term r/ at the end of (5) takes varios nmerical vales, which mst be obtained by analyzing detailed crrent measrements with modern technology. Sch vales are still navailable to the athor. First of all, we examine simple streamlines in Fig. 6. According to nmerical vales of r/, streamlines come in three types, as follows. (1) Normal-gradient crvatre type (concentric rings) is defined by r/ = 1. Althogh we have examined concentric rings in Section 7, to compare with eccentric types, this type is exemplified in Fig. 6(a) by the formlas r = r, dr = d, r/ = 1. (53) Here r is the radis of crvatre for a streamline, which passes throgh a reference point nder consideration (marked with in Fig. 6), and r is the radis of crvatre for a streamline, which passes throgh a point that is displaced by from the reference point. The +-axis is directed inward (otward) for cyclonic (anticyclonic) circlations for Fig. 6. Sbstittion of the last eqation of (53) into (5) yields / = / r, which redces (8) to (43). () Large-gradient crvatre type (eccentric rings) is defined by r/ < 1. An example is fond arond the bottom of eccentric circles in Fig. 6(b) with formlas r = r, dr = d < d, r/ =. (54) Sbstittion of the last eqation of (54) into (5) yields / = / r, which redces (8) to Z = f + /r + / r. (55) In the negative shear zone of anticyclonic concentric rings, two negative terms /r and / r in (43) redce vales of Z. Becase the radial shear term / r in (55) is twice / r in (43) for concentric rings, it works to redce Z frther. Hence, (55) shows that the large-gradient crvatre in anticyclonic rings is apt to have a negative vale of Z, and conseqently to bring an nstable state. Fig. 6. Examples of the relation between r and for concentric (a) and eccentric (b) rings with straight -axes along vertical diameters of circlar streamlines. +, r >, r > for cyclonic rings;, r <, r < for anticyclonic rings. Arond the bottom of eccentric arcs in Fig. 6(b) the interval between adjacent arcs is a half of that for concentric arcs in Fig. 6(a). This interval becomes longer toward both sides of the bottom center and then pward. When this interval for eccentric arcs is shorter than that for the concentric circles in Fig. 6(a), sch eccentric arcs belong to the large-gradient crvatre type. (3) Small-gradient crvatre type (eccentric rings) is defined by r/ > 1. An example is fond arond the top of eccentric circles in Fig. 6(b) with formlas r = r (/3), dr = (/3) d > d, r/ = /3. (56) Sbstittion of the last eqation of (56) into (5) yields / = (/3) / r, which redces (8) to Z = f + /r + (/3) / r. (57) Becase the radial shear term (/3) / r in (57) is smaller than / r in (43) for concentric rings, it works to relax a decrease in Z. Hence, (57) shows that the small-gradient crvatre in anticyclonic rings is apt to have a positive vale of Z, and conseqently to bring a stable state. Arond the top of eccentric circles in Fig. 6(b) the interval between adjacent arcs is 1.5 times that for concentric arcs in Fig. 6(a). This interval becomes shorter toward both sides of the top center and then downward. When this interval for eccentric arcs is longer than that for the concentric circles in Fig. 6(a), sch eccentric arcs belong to the small-gradient crvatre type. Inertial Instability of Arbitrarily Meandering Crrents 171

10 9. Remark What is called the zero-potential-vorticity assmption is eqivalent to the zero-absolte-vorticity assmption. Althogh sch an assmption has become an important method for analyzing nonlinear jets in rotating flids, the zero-absolte-vorticity assmption brings a netral state, according to inertial instability criteria for any type of crrents, sch as the zonal streams, concentric vortices and arbitrarily meandering crrents discssed in this paper. Becase a netral state is apt to lead to an nstable state owing to flctations, the zero-potentialvorticity jets are apt to become nstable. Acknowledgements The athor expresses his sincere thanks to anonymos referees who made critical bt constrctive comments on early drafts of this paper. References Brown, O. B., D. B. Olson, J. W. Brown and R. H. Evans (1983): Satellite infrared observation of the kinematics of a warmcore ring. Ast. J. Mar. Freshw. Res., 34, Carnevale, G. F., M. Briscolini, R. C. Kloosterziel and G. K. Vallis (1997): Three-dimensionally pertrbed vortex tbes in a rotating flow. J. Flid Mech., 341, Kawai, H. (1957): On the natral coordinate system and its applications to the Kroshio system. Bll. Tohok Regional Fisheries Res. Lab., 1, Kawai, H. (199): A model for the strctre of anticyclonic and cyclonic rings with niform dynamical energy density. Deep-Sea Res., 39, Sppl. 1, S1 S18. Kloosterziel, R. C. and G. J. F. van Heijst (1991): An experimental stdy of nstable barotropic vortices in a rotating flid. J. Flid Mech., 3, 1 4. Mied, R. P., G. J. Lindemann and J. M. Bergin (1983): Azimthal strctre of a cyclonic Glf Stream ring. J. Geophys. Res., 88, C4, Orlandi, P. and G. F. Carnevale (1999): Evoltion of isolated vortices in a rotating flid of finite depth. J. Flid Mech., 381, Potylitsin, P. G. and W. R. Peltier (1998): Stratification effects on the stability of colmnar vorticities on the f-plane. J. Flid Mech., 355, Rayleigh, Lord (1916): On the dynamics of revolving flids. Proc. R. Soc. Lond., A93, Tomosada, A. (1975): Observations of a warm eddy detached from the Kroshio east of Japan. Bll. Tokai Regional Fisheries Res. Lab., 81, (in Japanese with English abstract). 17 H. Kawai

Baroclinic Buoyancy-Inertia Joint Stability Parameter

Baroclinic Buoyancy-Inertia Joint Stability Parameter Jornal of Oceanography, Vol. 6, pp. 35 to 46, 5 Baroclinic Boyancy-Inertia Joint Stability Parameter HIDEO KAWAI* 3-8 Shibagahara, Kse, Joyo, Kyoto Pref. 6-, Japan (Receied 9 September 3; in reised form

More information

4 Primitive Equations

4 Primitive Equations 4 Primitive Eqations 4.1 Spherical coordinates 4.1.1 Usefl identities We now introdce the special case of spherical coordinates: (,, r) (longitde, latitde, radial distance from Earth s center), with 0

More information

1 Differential Equations for Solid Mechanics

1 Differential Equations for Solid Mechanics 1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to

More information

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1 y g j P(t) P(,y) r t0 i 4/4/006 Motion ( F.Robilliard) 1 Motion: We stdy in detail three cases of motion: 1. Motion in one dimension with constant acceleration niform linear motion.. Motion in two dimensions

More information

5. The Bernoulli Equation

5. The Bernoulli Equation 5. The Bernolli Eqation [This material relates predominantly to modles ELP034, ELP035] 5. Work and Energy 5. Bernolli s Eqation 5.3 An example of the se of Bernolli s eqation 5.4 Pressre head, velocity

More information

Setting The K Value And Polarization Mode Of The Delta Undulator

Setting The K Value And Polarization Mode Of The Delta Undulator LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions

More information

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University 9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)

More information

Prandl established a universal velocity profile for flow parallel to the bed given by

Prandl established a universal velocity profile for flow parallel to the bed given by EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient

More information

Efficiency Increase and Input Power Decrease of Converted Prototype Pump Performance

Efficiency Increase and Input Power Decrease of Converted Prototype Pump Performance International Jornal of Flid Machinery and Systems DOI: http://dx.doi.org/10.593/ijfms.016.9.3.05 Vol. 9, No. 3, Jly-September 016 ISSN (Online): 188-9554 Original Paper Efficiency Increase and Inpt Power

More information

Fluid Dynamics. Type of Flows Continuity Equation Bernoulli Equation Steady Flow Energy Equation Applications of Bernoulli Equation

Fluid Dynamics. Type of Flows Continuity Equation Bernoulli Equation Steady Flow Energy Equation Applications of Bernoulli Equation Tye of Flows Continity Eqation Bernolli Eqation Steady Flow Energy Eqation Alications of Bernolli Eqation Flid Dynamics Streamlines Lines having the direction of the flid velocity Flids cannot cross a

More information

FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL RESPONSE OF A CANTILEVER MODEL

FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL RESPONSE OF A CANTILEVER MODEL BBAA VI International Colloqim on: Blff Bodies Aerodynamics & Applications Milano, Italy, Jly, 2-24 28 FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL

More information

Gravitational Instability of a Nonrotating Galaxy *

Gravitational Instability of a Nonrotating Galaxy * SLAC-PUB-536 October 25 Gravitational Instability of a Nonrotating Galaxy * Alexander W. Chao ;) Stanford Linear Accelerator Center Abstract Gravitational instability of the distribtion of stars in a galaxy

More information

STATIC, STAGNATION, AND DYNAMIC PRESSURES

STATIC, STAGNATION, AND DYNAMIC PRESSURES STATIC, STAGNATION, AND DYNAMIC PRESSURES Bernolli eqation is g constant In this eqation is called static ressre, becase it is the ressre that wold be measred by an instrment that is static with resect

More information

OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1

OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1 OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIEL OF A POLYHERAL BOY WITH LINEARLY INCREASING ENSITY 1 V. POHÁNKA2 Abstract The formla for the comptation of the gravity field of a polyhedral body

More information

Experimental Study of an Impinging Round Jet

Experimental Study of an Impinging Round Jet Marie Crie ay Final Report : Experimental dy of an Impinging Rond Jet BOURDETTE Vincent Ph.D stdent at the Rovira i Virgili University (URV), Mechanical Engineering Department. Work carried ot dring a

More information

Section 7.4: Integration of Rational Functions by Partial Fractions

Section 7.4: Integration of Rational Functions by Partial Fractions Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to

More information

PLANETARY ORBITS. According to MATTER (Re-examined) Nainan K. Varghese,

PLANETARY ORBITS. According to MATTER (Re-examined) Nainan K. Varghese, PLANETARY ORBITS According to MATTER (Re-examined) Nainan K. Varghese, matterdoc@gmail.com http://www.matterdoc.info Abstract: It is an established fact that sn is a moving macro body in space. By simple

More information

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical

More information

1 JAXA Special Pblication JAXA-SP-1-E Small-scale trblence affects flow fields arond a blff body and therefore it governs characteristics of cross-sec

1 JAXA Special Pblication JAXA-SP-1-E Small-scale trblence affects flow fields arond a blff body and therefore it governs characteristics of cross-sec First International Symposim on Fltter and its Application, 1 11 IEXPERIMENTAL STUDY ON TURBULENCE PARTIAL SIMULATION FOR BLUFF BODY Hiroshi Katschi +1 and Hitoshi Yamada + +1 Yokohama National University,

More information

Elements of Coordinate System Transformations

Elements of Coordinate System Transformations B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and

More information

Experiment and mathematical model for the heat transfer in water around 4 C

Experiment and mathematical model for the heat transfer in water around 4 C Eropean Jornal of Physics PAPER Experiment and mathematical model for the heat transfer in water arond 4 C To cite this article: Naohisa Ogawa and Fmitoshi Kaneko 2017 Er. J. Phys. 38 025102 View the article

More information

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled. Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports

More information

MEASUREMENT OF TURBULENCE STATISTICS USING HOT WIRE ANEMOMETRY

MEASUREMENT OF TURBULENCE STATISTICS USING HOT WIRE ANEMOMETRY MEASUREMENT OF TURBULENCE STATISTICS USING HOT WIRE ANEMOMETRY Mrgan Thangadrai +, Atl Kmar Son *, Mritynjay Singh +, Sbhendra *, Vinoth Kmar ++, Ram Pyare Singh +, Pradip K Chatterjee + + Thermal Engineering,

More information

Momentum Equation. Necessary because body is not made up of a fixed assembly of particles Its volume is the same however Imaginary

Momentum Equation. Necessary because body is not made up of a fixed assembly of particles Its volume is the same however Imaginary Momentm Eqation Interest in the momentm eqation: Qantification of proplsion rates esign strctres for power generation esign of pipeline systems to withstand forces at bends and other places where the flow

More information

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA. PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 13 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA. PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 13 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS EDEXCEL NATIONAL CERTIFICATE/DIPLOMA PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 3 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS TUTORIAL - PIPE FLOW CONTENT Be able to determine the parameters of pipeline

More information

ρ u = u. (1) w z will become certain time, and at a certain point in space, the value of

ρ u = u. (1) w z will become certain time, and at a certain point in space, the value of THE CONDITIONS NECESSARY FOR DISCONTINUOUS MOTION IN GASES G I Taylor Proceedings of the Royal Society A vol LXXXIV (90) pp 37-377 The possibility of the propagation of a srface of discontinity in a gas

More information

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2 MATH 307 Vectors in Rn Dr. Neal, WKU Matrices of dimension 1 n can be thoght of as coordinates, or ectors, in n- dimensional space R n. We can perform special calclations on these ectors. In particlar,

More information

Second-Order Wave Equation

Second-Order Wave Equation Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order

More information

Spring Semester 2011 April 5, 2011

Spring Semester 2011 April 5, 2011 METR 130: Lectre 4 - Reynolds Averaged Conservation Eqations - Trblent Flxes (Definition and typical ABL profiles, CBL and SBL) - Trblence Closre Problem & Parameterization Spring Semester 011 April 5,

More information

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane Filomat 3:2 (27), 376 377 https://doi.org/.2298/fil7276a Pblished by Faclty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Conditions for Approaching

More information

called the potential flow, and function φ is called the velocity potential.

called the potential flow, and function φ is called the velocity potential. J. Szantr Lectre No. 3 Potential flows 1 If the flid flow is irrotational, i.e. everwhere or almost everwhere in the field of flow there is rot 0 it means that there eists a scalar fnction ϕ,, z), sch

More information

E ect Of Quadrant Bow On Delta Undulator Phase Errors

E ect Of Quadrant Bow On Delta Undulator Phase Errors LCLS-TN-15-1 E ect Of Qadrant Bow On Delta Undlator Phase Errors Zachary Wolf SLAC Febrary 18, 015 Abstract The Delta ndlator qadrants are tned individally and are then assembled to make the tned ndlator.

More information

Technical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty

Technical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte

More information

Curves - Foundation of Free-form Surfaces

Curves - Foundation of Free-form Surfaces Crves - Fondation of Free-form Srfaces Why Not Simply Use a Point Matrix to Represent a Crve? Storage isse and limited resoltion Comptation and transformation Difficlties in calclating the intersections

More information

Discontinuous Fluctuation Distribution for Time-Dependent Problems

Discontinuous Fluctuation Distribution for Time-Dependent Problems Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation

More information

Simplified Identification Scheme for Structures on a Flexible Base

Simplified Identification Scheme for Structures on a Flexible Base Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles

More information

08.06 Shooting Method for Ordinary Differential Equations

08.06 Shooting Method for Ordinary Differential Equations 8.6 Shooting Method for Ordinary Differential Eqations After reading this chapter, yo shold be able to 1. learn the shooting method algorithm to solve bondary vale problems, and. apply shooting method

More information

Subcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany

Subcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2

More information

FREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS

FREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS 7 TH INTERNATIONAL CONGRESS O THE AERONAUTICAL SCIENCES REQUENCY DOMAIN LUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS Yingsong G, Zhichn Yang Northwestern Polytechnical University, Xi an, P. R. China,

More information

The prediction of turbulence intensities in unsteady flow

The prediction of turbulence intensities in unsteady flow University of Wollongong Research Online Faclty of Engineering and Information Sciences - Papers: Part A Faclty of Engineering and Information Sciences 24 The prediction of trblence intensities in nsteady

More information

UNIT V BOUNDARY LAYER INTRODUCTION

UNIT V BOUNDARY LAYER INTRODUCTION UNIT V BOUNDARY LAYER INTRODUCTION The variation of velocity from zero to free-stream velocity in the direction normal to the bondary takes place in a narrow region in the vicinity of solid bondary. This

More information

Sareban: Evaluation of Three Common Algorithms for Structure Active Control

Sareban: Evaluation of Three Common Algorithms for Structure Active Control Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1638-1646 1638 Evalation of Three Common Algorithms for Strctre Active Control Mohammad Sareban Department of Civil Engineering Shahrood

More information

Applying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB

Applying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB IOS Jornal of Mathematics (IOS-JM) e-issn: 78-578, p-issn: 319-765X. Volme 13, Isse 6 Ver. II (Nov. - Dec. 17), PP 5-59 www.iosrjornals.org Applying Laminar and Trblent Flow and measring Velocity Profile

More information

The Dual of the Maximum Likelihood Method

The Dual of the Maximum Likelihood Method Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All

More information

Kragujevac J. Sci. 34 (2012) UDC 532.5: :537.63

Kragujevac J. Sci. 34 (2012) UDC 532.5: :537.63 5 Kragjevac J. Sci. 34 () 5-. UDC 53.5: 536.4:537.63 UNSTEADY MHD FLOW AND HEAT TRANSFER BETWEEN PARALLEL POROUS PLATES WITH EXPONENTIAL DECAYING PRESSURE GRADIENT Hazem A. Attia and Mostafa A. M. Abdeen

More information

Effects of modifications on the hydraulics of Denil fishways

Effects of modifications on the hydraulics of Denil fishways BOREAL ENVIRONMENT RESEARCH 5: 67 79 ISSN 1239-6095 Helsinki 28 March 2000 2000 Effects of modifications on the hydralics of Denil fishways Riitta Kamla 1) and Jan Bärthel 2) 1) Water Resorces and Environmental

More information

Obliqe Projection. A body is projected from a point with different angles of projections 0 0, 35 0, 45 0, 60 0 with the horizontal bt with same initial speed. Their respective horizontal ranges are R,

More information

Pulses on a Struck String

Pulses on a Struck String 8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a

More information

PIPELINE MECHANICAL DAMAGE CHARACTERIZATION BY MULTIPLE MAGNETIZATION LEVEL DECOUPLING

PIPELINE MECHANICAL DAMAGE CHARACTERIZATION BY MULTIPLE MAGNETIZATION LEVEL DECOUPLING PIPELINE MECHANICAL DAMAGE CHARACTERIZATION BY MULTIPLE MAGNETIZATION LEVEL DECOUPLING INTRODUCTION Richard 1. Davis & 1. Brce Nestleroth Battelle 505 King Ave Colmbs, OH 40201 Mechanical damage, cased

More information

Similarity Solution for MHD Flow of Non-Newtonian Fluids

Similarity Solution for MHD Flow of Non-Newtonian Fluids P P P P IJISET - International Jornal of Innovative Science, Engineering & Technology, Vol. Isse 6, Jne 06 ISSN (Online) 48 7968 Impact Factor (05) - 4. Similarity Soltion for MHD Flow of Non-Newtonian

More information

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method 17th World Conference on Nondestrctive Testing, 5-8 Oct 008, Shanghai, China Evalation of the Fiberglass-Reinforced Plastics Interfacial Behavior by sing Ultrasonic Wave Propagation Method Jnjie CHANG

More information

Chapter 1: Differential Form of Basic Equations

Chapter 1: Differential Form of Basic Equations MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)

More information

FEA Solution Procedure

FEA Solution Procedure EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations

More information

Krauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields.

Krauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields. Kraskopf, B, Lee,, & Osinga, H (28) odimension-one tangency bifrcations of global Poincaré maps of for-dimensional vector fields Early version, also known as pre-print Link to pblication record in Explore

More information

MODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION

MODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION MODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION Rssell Qadros, Krishnend Sinha Department of Aerospace Engineering Indian Institte of Technology Bombay Mmbai, India 476 Johan

More information

PhysicsAndMathsTutor.com

PhysicsAndMathsTutor.com C Integration - By sbstittion PhysicsAndMathsTtor.com. Using the sbstittion cos +, or otherwise, show that e cos + sin d e(e ) (Total marks). (a) Using the sbstittion cos, or otherwise, find the eact vale

More information

arxiv: v1 [physics.flu-dyn] 11 Mar 2011

arxiv: v1 [physics.flu-dyn] 11 Mar 2011 arxiv:1103.45v1 [physics.fl-dyn 11 Mar 011 Interaction of a magnetic dipole with a slowly moving electrically condcting plate Evgeny V. Votyakov Comptational Science Laboratory UCY-CompSci, Department

More information

4 Exact laminar boundary layer solutions

4 Exact laminar boundary layer solutions 4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved

More information

The spreading residue harmonic balance method for nonlinear vibration of an electrostatically actuated microbeam

The spreading residue harmonic balance method for nonlinear vibration of an electrostatically actuated microbeam J.L. Pan W.Y. Zh Nonlinear Sci. Lett. Vol.8 No. pp.- September The spreading reside harmonic balance method for nonlinear vibration of an electrostatically actated microbeam J. L. Pan W. Y. Zh * College

More information

The Heat Equation and the Li-Yau Harnack Inequality

The Heat Equation and the Li-Yau Harnack Inequality The Heat Eqation and the Li-Ya Harnack Ineqality Blake Hartley VIGRE Research Paper Abstract In this paper, we develop the necessary mathematics for nderstanding the Li-Ya Harnack ineqality. We begin with

More information

Study of the diffusion operator by the SPH method

Study of the diffusion operator by the SPH method IOSR Jornal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-684,p-ISSN: 2320-334X, Volme, Isse 5 Ver. I (Sep- Oct. 204), PP 96-0 Stdy of the diffsion operator by the SPH method Abdelabbar.Nait

More information

Transient Approach to Radiative Heat Transfer Free Convection Flow with Ramped Wall Temperature

Transient Approach to Radiative Heat Transfer Free Convection Flow with Ramped Wall Temperature Jornal of Applied Flid Mechanics, Vol. 5, No., pp. 9-1, 1. Available online at www.jafmonline.net, ISSN 175-57, EISSN 175-645. Transient Approach to Radiative Heat Transfer Free Convection Flow with Ramped

More information

PhysicsAndMathsTutor.com

PhysicsAndMathsTutor.com . Two smooth niform spheres S and T have eqal radii. The mass of S is 0. kg and the mass of T is 0.6 kg. The spheres are moving on a smooth horizontal plane and collide obliqely. Immediately before the

More information

The Linear Quadratic Regulator

The Linear Quadratic Regulator 10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.

More information

CONSIDER an array of N sensors (residing in threedimensional. Investigating Hyperhelical Array Manifold Curves Using the Complex Cartan Matrix

CONSIDER an array of N sensors (residing in threedimensional. Investigating Hyperhelical Array Manifold Curves Using the Complex Cartan Matrix IEEE JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING (SPECIAL ISSUE), VOL.?, NO.?,???? 1 Investigating yperhelical Array Manifold Crves Using the Complex Cartan Matrix Athanassios Manikas, Senior Member,

More information

The Bow Shock and the Magnetosheath

The Bow Shock and the Magnetosheath Chapter 6 The Bow Shock and the Magnetosheath The solar wind plasma travels sally at speeds which are faster than any flid plasma wave relative to the magnetosphere. Therefore a standing shock wave forms

More information

Material Transport with Air Jet

Material Transport with Air Jet Material Transport with Air Jet Dr. István Patkó Bdapest Tech Doberdó út 6, H-1034 Bdapest, Hngary patko@bmf.h Abstract: In the field of indstry, there are only a very few examples of material transport

More information

Frequency Estimation, Multiple Stationary Nonsinusoidal Resonances With Trend 1

Frequency Estimation, Multiple Stationary Nonsinusoidal Resonances With Trend 1 Freqency Estimation, Mltiple Stationary Nonsinsoidal Resonances With Trend 1 G. Larry Bretthorst Department of Chemistry, Washington University, St. Lois MO 6313 Abstract. In this paper, we address the

More information

Turbulence and boundary layers

Turbulence and boundary layers Trblence and bondary layers Weather and trblence Big whorls hae little whorls which feed on the elocity; and little whorls hae lesser whorls and so on to iscosity Lewis Fry Richardson Momentm eqations

More information

Numerical Simulation of Three Dimensional Flow in Water Tank of Marine Fish Larvae

Numerical Simulation of Three Dimensional Flow in Water Tank of Marine Fish Larvae Copyright c 27 ICCES ICCES, vol.4, no.1, pp.19-24, 27 Nmerical Simlation of Three Dimensional Flo in Water Tank of Marine Fish Larvae Shigeaki Shiotani 1, Atsshi Hagiara 2 and Yoshitaka Sakakra 3 Smmary

More information

IJAET International Journal of Application of Engineering and Technology ISSN: Vol.1 No.1

IJAET International Journal of Application of Engineering and Technology ISSN: Vol.1 No.1 IJAET International Jornal of Application of Engineering and Technology ISSN: 395-3594 Vol1 No1 ANALYSIS OF SUPERSONIC FLOWS IN THE E -LAVAL NOZZLE AT 1 INTO A SUENLY EXPANE UCT AT L/=WITH CAVITY ASPECT

More information

Restricted Three-Body Problem in Different Coordinate Systems

Restricted Three-Body Problem in Different Coordinate Systems Applied Mathematics 3 949-953 http://dx.doi.org/.436/am..394 Pblished Online September (http://www.scirp.org/jornal/am) Restricted Three-Body Problem in Different Coordinate Systems II-In Sidereal Spherical

More information

V. Hadron quantum numbers

V. Hadron quantum numbers V. Hadron qantm nmbers Characteristics of a hadron: 1) Mass 2) Qantm nmbers arising from space-time symmetries : total spin J, parity P, charge conjgation C. Common notation: 1 -- + 2 J P (e.g. for proton:

More information

Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Porous Rotating Disk

Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Porous Rotating Disk Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Poros otating Disk Hemant Poonia * * Asstt. Prof., Deptt. of Math, Stat & Physics, CCSHAU, Hisar-54.. C. Chadhary etd. Professor, Deptt.

More information

BLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students

BLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they

More information

FOUNTAIN codes [3], [4] provide an efficient solution

FOUNTAIN codes [3], [4] provide an efficient solution Inactivation Decoding of LT and Raptor Codes: Analysis and Code Design Francisco Lázaro, Stdent Member, IEEE, Gianligi Liva, Senior Member, IEEE, Gerhard Bach, Fellow, IEEE arxiv:176.5814v1 [cs.it 19 Jn

More information

Dynamics of a Holling-Tanner Model

Dynamics of a Holling-Tanner Model American Jornal of Engineering Research (AJER) 07 American Jornal of Engineering Research (AJER) e-issn: 30-0847 p-issn : 30-0936 Volme-6 Isse-4 pp-3-40 wwwajerorg Research Paper Open Access Dynamics of

More information

Chapter 9 Flow over Immersed Bodies

Chapter 9 Flow over Immersed Bodies 57:00 Mechanics o Flids and Transport Processes Chapter 9 Proessor Fred Stern Fall 01 1 Chapter 9 Flow over Immersed Bodies Flid lows are broadly categorized: 1. Internal lows sch as dcts/pipes, trbomachinery,

More information

Nonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator

Nonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator Proceedings of the 6 IEEE International Conference on Control Applications Mnich, Germany, October 4-6, 6 WeB16 Nonparametric Identification and Robst H Controller Synthesis for a Rotational/Translational

More information

Determining of temperature field in a L-shaped domain

Determining of temperature field in a L-shaped domain Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 0, (:-8 Determining of temperatre field in a L-shaped domain Oigo M. Zongo, Sié Kam, Kalifa Palm, and Alione Oedraogo

More information

The Theory of Virtual Particles as an Alternative to Special Relativity

The Theory of Virtual Particles as an Alternative to Special Relativity International Jornal of Physics, 017, Vol. 5, No. 4, 141-146 Available online at http://pbs.sciepb.com/ijp/5/4/6 Science and Edcation Pblishing DOI:10.1691/ijp-5-4-6 The Theory of Virtal Particles as an

More information

Lateral Load Capacity of Piles

Lateral Load Capacity of Piles Lateral Load Capacity of Piles M. T. DAVSSON, Department of Civil Engineering, University of llinois, Urbana Pile fondations sally find resistance to lateral loads from (a) passive soil resistance on the

More information

ECON3120/4120 Mathematics 2, spring 2009

ECON3120/4120 Mathematics 2, spring 2009 University of Oslo Department of Economics Arne Strøm ECON3/4 Mathematics, spring 9 Problem soltions for Seminar 4, 6 Febrary 9 (For practical reasons some of the soltions may inclde problem parts that

More information

A Model-Free Adaptive Control of Pulsed GTAW

A Model-Free Adaptive Control of Pulsed GTAW A Model-Free Adaptive Control of Plsed GTAW F.L. Lv 1, S.B. Chen 1, and S.W. Dai 1 Institte of Welding Technology, Shanghai Jiao Tong University, Shanghai 00030, P.R. China Department of Atomatic Control,

More information

1. INTRODUCTION. A solution for the dark matter mystery based on Euclidean relativity. Frédéric LASSIAILLE 2009 Page 1 14/05/2010. Frédéric LASSIAILLE

1. INTRODUCTION. A solution for the dark matter mystery based on Euclidean relativity. Frédéric LASSIAILLE 2009 Page 1 14/05/2010. Frédéric LASSIAILLE Frédéric LASSIAILLE 2009 Page 1 14/05/2010 Frédéric LASSIAILLE email: lmimi2003@hotmail.com http://lmi.chez-alice.fr/anglais A soltion for the dark matter mystery based on Eclidean relativity The stdy

More information

3.4-Miscellaneous Equations

3.4-Miscellaneous Equations .-Miscellaneos Eqations Factoring Higher Degree Polynomials: Many higher degree polynomials can be solved by factoring. Of particlar vale is the method of factoring by groping, however all types of factoring

More information

Regression Analysis of Octal Rings as Mechanical Force Transducers

Regression Analysis of Octal Rings as Mechanical Force Transducers Regression Analysis of Octal Rings as Mechanical Force Transdcers KHALED A. ABUHASEL* & ESSAM SOLIMAN** *Department of Mechanical Engineering College of Engineering, University of Bisha, Bisha, Kingdom

More information

The Real Stabilizability Radius of the Multi-Link Inverted Pendulum

The Real Stabilizability Radius of the Multi-Link Inverted Pendulum Proceedings of the 26 American Control Conference Minneapolis, Minnesota, USA, Jne 14-16, 26 WeC123 The Real Stabilizability Radis of the Mlti-Link Inerted Pendlm Simon Lam and Edward J Daison Abstract

More information

WEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY

WEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY APPLIED PHYSICS MEDICAL WEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY L. CÃPITANU, A. IAROVICI, J. ONIªORU Institte of Solid Mechanics, Romanian Academy, Constantin Mille 5, Bcharest Received

More information

Comments on Vertical Vorticity Advection

Comments on Vertical Vorticity Advection Comments on Vertical Vorticity Advection It shold be fairly intitive that ositive maima in vertical vorticity are associated with cyclones, and ths ositive cyclonic vorticity advection might be a sefl

More information

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018 Lectre 3 The dot prodct Dan Nichols nichols@math.mass.ed MATH 33, Spring 018 Uniersity of Massachsetts Janary 30, 018 () Last time: 3D space Right-hand rle, the three coordinate planes 3D coordinate system:

More information

Incompressible Viscoelastic Flow of a Generalised Oldroyed-B Fluid through Porous Medium between Two Infinite Parallel Plates in a Rotating System

Incompressible Viscoelastic Flow of a Generalised Oldroyed-B Fluid through Porous Medium between Two Infinite Parallel Plates in a Rotating System International Jornal of Compter Applications (97 8887) Volme 79 No., October Incompressible Viscoelastic Flow of a Generalised Oldroed-B Flid throgh Poros Medim between Two Infinite Parallel Plates in

More information

EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE

EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE 13 th World Conference on Earthqake Engineering Vancover, B.C., Canada Agst 1-6, 2004 Paper No. 3099 EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE Ellen M. RATHJE 1, Wen-Jong CHANG 2, Kenneth

More information

An Investigation into Estimating Type B Degrees of Freedom

An Investigation into Estimating Type B Degrees of Freedom An Investigation into Estimating Type B Degrees of H. Castrp President, Integrated Sciences Grop Jne, 00 Backgrond The degrees of freedom associated with an ncertainty estimate qantifies the amont of information

More information

Numerical verification of the existence of localization of the elastic energy for closely spaced rigid disks

Numerical verification of the existence of localization of the elastic energy for closely spaced rigid disks Nmerical verification of the existence of localization of the elastic energy for closely spaced rigid disks S. I. Rakin Siberian State University of transport Rssia, 6349, Novosibirsk, Dsy Kovalchk street,

More information

Cosmic Microwave Background Radiation. Carl W. Akerlof April 7, 2013

Cosmic Microwave Background Radiation. Carl W. Akerlof April 7, 2013 Cosmic Microwave Backgrond Radiation Carl W. Akerlof April 7, 013 Notes: Dry ice sblimation temperatre: Isopropyl alcohol freezing point: LNA operating voltage: 194.65 K 184.65 K 18.0 v he terrestrial

More information

Generalized Jinc functions and their application to focusing and diffraction of circular apertures

Generalized Jinc functions and their application to focusing and diffraction of circular apertures Qing Cao Vol. 20, No. 4/April 2003/J. Opt. Soc. Am. A 66 Generalized Jinc fnctions and their application to focsing and diffraction of circlar apertres Qing Cao Optische Nachrichtentechnik, FernUniversität

More information

Conceptual Questions. Problems. 852 CHAPTER 29 Magnetic Fields

Conceptual Questions. Problems. 852 CHAPTER 29 Magnetic Fields 852 CHAPTER 29 Magnetic Fields magnitde crrent, and the niform magnetic field points in the positive direction. Rank the loops by the magnitde of the torqe eerted on them by the field from largest to smallest.

More information

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers D.R. Espinoza-Trejo and D.U. Campos-Delgado Facltad de Ingeniería, CIEP, UASLP, espinoza trejo dr@aslp.mx Facltad de Ciencias,

More information