ORE Open Research Exeter
|
|
- Harriet Carter
- 5 years ago
- Views:
Transcription
1 ORE Open Reearch Exeter TITLE The onet of convection in rotating circular cylinder with experimental boundary condition AUTHORS Zhang, Keke; Liao, X. JOURNAL Journal of Fluid Mechanic DEPOSITED IN ORE 4 March 213 Thi verion available at COPYRIGHT AND REUSE Open Reearch Exeter make thi work available in accordance with publiher policie. A NOTE ON VERSIONS The verion preented here may differ from the publihed verion. If citing, you are advied to conult the publihed verion for pagination, volume/iue and date of publication
2 J. Fluid Mech. (29), vol. 622, pp c 29 Cambridge Univerity Pre doi:1.117/s x Printed in the United Kingdom 63 The onet of convection in rotating circular cylinder with experimental boundary condition KEKE ZHANG 1 AND XINHAO LIAO 2 1 Department of Mathematical Science, Univerity of Exeter, Exeter EX4 4QE, UK 2 Shanghai Atronomical Obervatory, Chinee Academy of Science, Shanghai 23, China (Received 1 Augut 28 and in revied form 27 October 28) Convective intabilitie in a fluid-filled circular cylinder heated from below and rotating about it vertical axi are invetigated both analytically and numerically under experimental boundary condition. It i found that there exit two different form of convective intabilitie: convection-driven inertial wave for mall and moderate Prandtl number and wall-localized travelling wave for large Prandtl number. Aymptotic olution for both form of convection are derived and numerical imulation for the ame problem are alo performed, howing a atifactory quantitative agreement between the aymptotic and numerical analye. 1. Introduction Motivated by the wih to undertand the fundamental dynamic taking place in planetary fluid interior and atmophere, convection in rotating cylindrical geometry ha been extenively tudied both experimentally and theoretically (ee, for example, Daniel 198; Zhong, Ecke & Steinberg 1991; Herrmann & Bue 1993; Kuo & Cro 1993; Aurnou et al. 23; Plaut 23; Young & Read 28). In particular, careful numerical tudie of convective intabilitie in a fluid contained in a circular cylinder uniformly heated from below and rotating about it vertical axi (Goldtein et al. 1993, 1994) reveal an exceedingly complicated and eemingly incomprehenible behaviour at mall and moderate Prandtl number. We have recently undertaken the tudy of convective intabilitie in a rotating circular cylinder with tre-free end (Zhang, Liao & Bue 27). The aumption of the tre-free boundary condition permit eparable olution and, conequently, lead to a much implified mathematical problem governed by ordinary differential equation. It reveal that inertial convection i preferred for fluid with mall and moderate Prandtl number if the tre-free condition i aumed. Since inertial convection i energetically driven by thermal buoyancy againt vicou diipation that occur largely in the vicou boundary layer (ee, for example, Zhang 1995; Bue & Simitev 24), the type of velocity boundary condition play, both phyically and mathematically, an eential role in determining the main propertie of convection. The Ekman boundary layer reulting from the non-lip top and bottom not only introduce trong ma fluxe into the interior but alo bring dominant term in the aymptotic addre for correpondence: kzhang@ex.ac.uk
3 64 K. Zhang and X. Liao olution of inertial convection in rapidly rotating cylinder. It i therefore deirable to gain a better undertanding of convection under experimental boundary condition, the problem then being governed by partial differential equation, providing a ueful theoretical framework and guidance for experimental tudie of the problem. Thu the primary objective of thi tudy i to undertand, through both numerical and aymptotic analyi, convective intabilitie in a rotating cylinder under experimental boundary condition. Our focu will be mainly on intabilitie at mall and moderate Prandtl number. However, the aymptotic olution for large Prandtl number i alo briefly dicued for the purpoe of completene. It i fortunate that on recognizing and utilizing the phyical nature of convective intabilitie, we are able to derive everal relatively imple aymptotic formula that correctly decribe the key feature of convection for all value of the Prandtl number and the apect ratio and which agree atifactorily with the reult of numerical analyi. In what follow we hall begin by preenting the mathematical equation in 2. The numerical analyi i briefly preented in 3, while the correponding aymptotic analyi i dicued in 4, with a ummary and everal remark given in Mathematical formulation Conideration i given to convective intabilitie in a Bouineq fluid confined in a circular cylinder of radiu o and height d with contant thermal diffuivity κ, thermal expanion coefficient α and kinematic vicoity ν. The circular cylinder rotate uniformly about it axi with a contant vertical angular velocity Ω in the preence of a contant gravitational field g = gẑ, and i heated uniformly from below to produce an untable vertical temperature gradient, T = β ẑ, where β i a poitive contant. Cylindrical coordinate (, φ, z) with the correponding unit vector (ŝ, ˆφ, ẑ, whereẑ i parallel to the axi of rotation) are ued. Making ue of the depth d a the length cale, Ω 1 a the unit of time and βd 3 Ω/κ a the unit of temperature fluctuation, the problem of linear convective intability i governed by the equation u t +2ẑ u = p + RΘẑ + E 2 u, (2.1) u =, (2.2) Pr Θ = E(ẑ u + 2 Θ), (2.3) t where u i the three-dimenional velocity field, and the temperature deviation from the purely conductive tate T (z) i repreented by Θ. In cylindrical coordinate, we write olution of the equation in the form (u,θ)(, φ, z, t)= (u,θ) (, z) exp [i(2σt + mφ), where m i the azimuthal wavenumber aumed to be poitive and σ i the half-frequency which i ued becaue of it property σ < 1 in the limit E =. The three non-dimenional phyical parameter the Rayleigh number R, the Prandtl number Pr andtheekmannumbere are defined a R = αβgd2 Ωκ, Pr = ν κ, E = ν Ωd. 2 The geometric parameter i given by the apect ratio, defined a = /d. Appropriately for experiment of the problem (ee, for example, Goldtein et al. 1993), we hall aume no-lip velocity and perfectly conducting condition on the
4 Convection in rotating circular cylinder 65 top and bottom of the cylinder, u = Θ = at z =, 1. (2.4) On the idewall, we aume no-lip and perfectly inulating condition which require u = Θ = at =. (2.5) Equation (2.1) (2.3) ubject to the boundary condition in (2.4) (2.5) form the convective intability problem which will be olved by both aymptotic and numerical analyi at, a in practical computation and phyical experiment, a mall but fixed E for different value of Pr. 3. Numerical analyi Compared to the numerical analyi for tre-free end, which allow olution in eparable variable, the numerical analyi for the no-lip condition require olution of partial differential equation and hence i much more complicated. In order to determine the mot untable mode of convective intabilitie, we mut conider diturbance of all poible azimuthal wavenumber. Axiymmetric convection may be phyically preferred in rotating pherical ytem (ee, for example, Net, Carcia &Sánchez 28). However, our extenive numerical and aymptotic tudie how that axiymmetric (m =) intabilitie are not preferred for any value of Pr when E 1and = O(1). In conequence, we hall preent the numerical and aymptotic analyi only for non-axiymmetric (m ) olution. For non-axiymmetric intabilitie with m, the velocity vector u atifying (2.2) can be expreed in term of two calar potential Ψ and Φ (Marqué 199): u = 1 ( Ψ Φ φ ŝ + z Ψ ) ˆφ 1 Φ ẑ. (3.1) φ In term of Ψ and Φ, the non-lip condition on the idewall and the end become Ψ = Ψ = Φ = at = ; Ψ = Φ z = Φ = at z =, 1. (3.2) Making ue of the expreion (3.1) and applying ẑ and ŝ onto (2.1), we can derive the three independent non-dimenional calar equation: ( )[ ( 1 t E 2 Φ ) ) ( 2 2 Ψ Φ =, (3.3) z z 2 z φ [ ( t E Ψ z φ 2E + 1 ) [ 2 Ψ 2 z [ ( ) 1 2 Φ z 2 [ E 2 Pr t ( 2 1 ( 2 2 Θ E Φ φ z 2 ) Ψ z ) Φ R Θ φ =, (3.4) =. (3.5) The above equation are then olved numerically by uing the Galerkin-type method in which, for example, Ψ and Φ are expanded in term of the Chebyhev function
5 66 K. Zhang and X. Liao R c (1. 1 6, 3,1) (1. 1 4, 3,1) (1. 1 3, 5,1) ( , 8,1) (5. 1 3, 7,2) (1. 1 2, 5,4) ( , 4,6) ( , 3,8) Pr Figure 1. The critical Rayleigh number R c i plotted a a function of Pr in a rotating cylinder with =2andE =1 4. Solid line repreent the numerical olution while dahed line reult from the aymptotic olution. Of the three number in the bracket, the firt i Pr, the econd i m and the third i the number of node in the direction. T j (x): [ J Ψ = m (1 š) 2 (1 ž 2 ) j= L Ψ mj l T j (ž) T l (š) e i(mφ+2σt), (3.6) l= [ J L Φ = m+1 (1 š)(1 ž 2 ) 2 Φ mj l T j (ž) T l (š) e i(mφ+2σt), (3.7) j= l= where ž =2z 1, š =2/ 1andJ = L = O(1) for achieving an accuracy within 1 %, the factor m and m+1 are impoed for regularity of the olution at = and the complex coefficient like Ψ mj l are obtained, together with the Rayleigh number R and the half-frequency σ, by a tandard numerical procedure. The primary aim of our numerical analyi i to provide a valuable comparion with the reult of the aymptotic analyi valid only for mall E. Some typical reult of the numerical analyi, along with the aymptotic reult, are preented in figure 1, howing the critical Rayleigh number R c for the mot untable mode of convective intabilitie a a function of Pr in a rotating cylinder with =2 and E =1 4. The convection pattern for three typical Prandtl number, Pr=.25,.25 and 1., are depicted in figure 2. We hall dicu the relevant detail when comparing them with the aymptotic olution. 4. Aymptotic analyi A key aumption in the aymptotic analyi i that the velocity of convection at leading order for mall or moderate Pr can be repreented by a ingle inertial-wave mode while buoyancy force appear only at the next order to drive the inertial wave againt the effect of vicou damping, leading to an aymptotic expanion in the form } u = u + u b + u 1, p= p + p b + p 1, (4.1) σ = σ + σ 1, Θ= Θ +, R = R 1 +, where u 1, p 1 and σ 1 repreent mall interior vicou correction to the leading-order invicid interior olution u, p and σ, while u b and p b denote olution of the vicou boundary layer which are non-zero only in the vicinity of the bounding urface of
6 Convection in rotating circular cylinder 67 (a) (c) (e) (b) z (d) ( f ) Figure 2. Contour of u φ in a horizontal plane at z =.75 (upper panel) and in a vertical plane (lower panel): (a,b) forpr =.25, (c,d) forpr =.25 and (e,f )forpr =1. with =2 and E =1 4. the cylinder. More preciely, u = O(1), u b = O(1) and u 1 1, with the flux from the ocillatory boundary layer u b providing the matching condition for the econdary interior flow u 1.ItiofimportancetonotethatE 1/2 cannot be ued a an expanion parameter in thi problem becaue the thickne of an ocillatory vicou boundary layer i not generally of order E 1/2, the flux from the ocillatory boundary layer i not generally of order E 1/2 and the internal vicou contribution i not generally of the order E 1/2 maller than that of the ocillatory boundary layer (Zhang & Liao 28). In other word, the effect of patial and temporal non-uniformitie obcure the form of the aymptotic expanion even for the firt vicou corrective term. After ubtitution of the expanion into (2.1) (2.3), the leading-order problem, ubject to the uual invicid boundary condition, i given by u +2ẑ u + p =, (4.2) t u =, (4.3) Pr Θ = E(ẑ u + 2 Θ ), (4.4) t which decribe non-diipative thermal inertial wave with olution in the form [ i ŝ u = 2 ( σ ξ ) 1 σ 2 [ 1 ξ ˆφ u = 2 ( ) 1 σ 2 J m 1( )+ m(1 σ ) J m ( ) J m 1( ) m(1 σ ) J m ( ) co(πz) e i(mφ+2σt), (4.5) co(πz) e i(mφ+2σ t), (4.6) ẑ u = iπ J m ( ) in(πz) e i(mφ+2σt), 2σ (4.7) N 2im ( ξ 2 βmn) 2 Qmn J m (ξ) in(πz) J m (β mn / ) Θ = e i(mφ+2σt), π 2 [π 2 +(β mn / ) 2 +2PrE 1 σ i J m (β mn ) (4.8) n=1 p = J m ( )co(πz) e i(mφ+2σ t), (4.9)
7 68 K. Zhang and X. Liao where = ξ/ (the implet vertical tructure that i relevant to convective intabilitie i taken), J m (x) denote the tandard Beel function, β mn are the root of J m(β mn ) = with <β m1 <β m2 <β m3, Q mn =(πβ mn ) 2 /[2σ 2(ξ 2 βmn) 2 2 (βmn 2 m 2 ), N O(1) for achieving an accuracy within 1 % and ξ i a olution of [ ξj m 1 (ξ)+ σ ( ) 2 1/2 ξ 1+ 1 σ π mj m(ξ) =, (4.1) σ = ±[1 + (ξ/π) 2 1/2. (4.11) Different value of ξ or σ correpond to different inertial mode. At the next order, the equation governing the perturbation of the interior flow are i2σ u 1 +2ẑ u 1 + p 1 = R 1 ẑθ + E 2 u i2σ 1 u, (4.12) u 1 =. (4.13) It i traightforward to derive the olvability condition required for the inhomogeneou equation (4.12): [ E 1/2 p ˆn ( ˆn u b )dη ds = S V [ R1 u ẑθ 2iσ 1 u 2 E u 2 dv, (4.14) where u and p denote the complex conjugate of u and p, V volume integral over the cylinder, repreent a urface integral over it bounding S urface with the normal vector ˆn and η i the tretched boundary-layer coordinate. In deriving (4.14), the ma flux from the Ekman boundary layer i ued to provide the matching condition for u 1 at the outer edge of the layer. Each integral in (4.14) can be evaluated. Upon uing (4.5) (4.8) together with the propertie of the Beel function, we can how that N [ 2 u πq mn mjm (ξ) ẑθ dv =, π 2 +(β mn / ) 2 +2iPrE 1 σ σ (4.15) V V n=1 u 2 dv = σ 2 π 2 V u 2 dv = π[(π )2 + m(m σ )J ( ) m(ξ) 2. (4.16) 4σ 2 1 σ 2 In order to evaluate the urface integral on the left-hand ide of (4.14), we mut derive the boundary-layer flow u b on the bounding urface of a cylinder. Conider firt the Ekman boundary layer on the bottom urface at z = decribed by ) ) (i2σ 2 ẑ u η 2 b =2u b, (i2σ 2 u η 2 b = 2ẑ u b, (4.17) where η = E 1/2 z and u b = u at z = o that the no-lip condition i obeyed. Note that η = i at the bottom urface while η = i at the outer edge of the boundary layer but till located at the bottom in term of the coordinate z. The boundary-layer olution to (4.17) on the bottom urface can be expreed in the form {[ 1 u b = 4 ( ξ(σ 1) ) J 1 σ 2 m 1 ( )+ 2m(1 σ ) J m ( ) (iŝ + ˆφ)e (1+i)η 1+σ [ } ξ(σ +1) + J m 1 ( ) (iŝ ˆφ)e (1 i)η 1 σ e i(mφ+2σt), (4.18)
8 Convection in rotating circular cylinder 69 by which, together with (4.9), we can obtain the contribution from the vicou boundary layer at z = to the integral on the left-hand ide of (4.14): 2π [ [ (p) πj 2 z= ˆn ( ˆn u b )dη d dφ = ( m(ξ) ) 1 σ 2 8σ 2 {( 1+σ + 1 σ )C 1 2σ ( 1+σ 1 σ )C 2 +i[( 1+σ 1 σ )C 1 2σ ( 1+σ + 1 σ )C 2 }, (4.19) where C 1 =(1+σ 2)(m2 + π 2 2 ) 2mσ and C 2 = π m(m σ ). Becaue of the ymmetry between the top and the bottom of the cylinder, the contribution from the vicou boundary layer at z = 1 i exactly the ame a that from the boundary layer at z =. Similarly, we can derive the contribution from the idewall boundary layer to the integral on the left-hand ide of (4.14): 2π 1 [ (p) = ˆn ( ˆn u b )dη dz dφ = π 4σ σ ( i+ σ ) (m π 2 )J 2 σ m(ξ). (4.2) Subtitution of (4.15), (4.16), (4.19) and (4.2) into the olvability condition (4.14) yield a complex equation whoe real part determine R 1 while the imaginary part give rie to the vicou correction σ 1 : [ E 1/2 N ( π 2 + β R 1 = mn/ 2 2) 1 Q mn 4m 2 ( 1 σ 2 n=1 π2 + βmn/ 2 2) 2 +(2σ Pr/E) 2 { E 1/2 π 2 [(π ) 2 + m(m σ ) 2σ m[ 1+σ + 1 σ 1 σ 2 σ 2 [ } +(m 2 + π 2 2 σ (1 σ 2 ) ) +(1+σ ) 3/2 +(1 σ ) 3/2, (4.21) ( ) 1 σ 2 1/2E {[ 1/2 N 2σ =2σ Q mn π 2 2 ( + m(m σ ) π2 + βmn/ 2 2) 2 +(2σ Pr/E) 2 n=1 [8m 2 R PrE 3/2 σ ( 1 σ 2 ) 1/2 +2σ m[ 1+σ 1 σ +(m 2 + π 2 2 ) [ σ 1 σ 2 σ } (1 + σ ) 3/2 +(1 σ ) 3/2. (4.22) In addition to the expreion for the Rayleigh number R 1 and the frequency σ,the leading-order velocity of convection can be explicitly written a
9 7 K. Zhang and X. Liao { [ i ŝ u = 4 ( σ ξ ) 2 1 σ 2 J m 1 ( ) + m(1 σ ) J m ( ) co πz [ ξ(σ 1) + J m 1 ( ) + 2m(1 σ ) [e χ J m ( ) +z e χ + (1 z) + ξ(σ +1) J m 1 ( ) [ e χ z e } χ (1 z) e i(2σt+mφ), (4.23) { [ 1 ˆφ u = 4 ( ξjm 1 ( ) ) 2 m(1 σ ) J 1 σ 2 m ( ) co πz + mj m ( ) ) co πze χ( 2σ [ ξ(σ 1) + J m 1 ( ) + 2m(1 σ ) [e χ J m ( ) +z e χ + (1 z) ξ(σ +1) J m 1 ( ) [ e χ z e } χ (1 z) e i(2σt+mφ), (4.24) ẑ u = iπ 2σ [ Jm ( ) J m (ξ) e χ( ) in (πz) e i(2σt+mφ), (4.25) where χ + =(1 + i) E 1 (1 + σ ), χ =(1 i) E 1 (1 σ ) and χ = σ E 1 (1 + iσ / σ ). A complete analytical olution of the problem i decribed by (4.21) (4.25). A indicated by (4.21), the Rayleigh number R 1 required to initiate inertial convection increae rapidly with increaing Pr and the idewall-localized convection would become phyically preferred at large Pr. It i worth mentioning that the aymptotic analyi (Herrmann & Bue 1993) wa concerned with the leading-order olution for which the effect of the non-lip condition wa not conidered. Since cylindrical geometry play an inignificant role in the wall-localized convective flow, we can imply re-cale the aymptotic olution in a channel obtained by Liao, Zhang & Chang (26), giving rie to R =2π E 1/3, (4.26) σ =[π E/Pr 29.6E 4/3 Pr 1, (4.27) m = [ π E 1/3, (4.28) which ha taken into account the non-lip boundary condition while it leading order i the ame a that given by Herrmann & Bue (1993). For any given value of, Pr and E, the following three tep are uually required to find the mot untable mode of convective intabilitie: (a) ue (4.21) to calculate the mallet R 1 required for the non-axiymmetric inertial convection over different value of σ,i.e. over different inertial mode; (b) ue (4.26) to calculate the econd Rayleigh number required for the wall-localized convection and (c) compare the two Rayleigh number obtained at exactly the ame, Pr and E. Of phyical ignificance i the minimum value of the Rayleigh number, which i denoted by R c and will be referred to a the critical Rayleigh number, at which convective intabilitie firt et in when the Rayleigh number R gradually increae from zero. The R c a a function of Pr, calculated from both the aymptotic and the numerical analyi, i hown in figure 1 for E =1 4 and = 2: all the mot untable mode of convection are non-axiymmetric. The azimuthal wavenumber
10 Convection in rotating circular cylinder 71 m c and the correponding node in the radial direction for everal typical Prandtl number are alo hown in figure 1. For example, the aymptotic formula (4.21) (4.22) at Pr=.25 yield R c =3.41,m c =8,σ c =.474 with σ =.4773, while the numerical analyi give R c =3.58,m c =8,σ c =.473, with the patial tructure of it flow diplayed in figure 2(a,b). Note that there are no noticeable difference between the numerical olution and the correponding analytical olution (4.23) (4.25) when they are plotted. When Pr increae, the critical wavenumber m c may decreae at the expene of increaing the radial complexity of the flow. At Pr=.25, the aymptotic expreion (4.21) (4.22) give rie to R c =2.4,m c =4,σ c =.248 with σ =.2634 and ix node in the radial direction, whoe patial tructure i hown in figure 2(c,d). The correponding numerical olution i characterized by R c =21.6,m c =4,σ c =.247, howing nearly the identical patial tructure. When Pr increae to about.7 at E =1 4, a hown in figure 1, the aymptotic olution given by (4.21), which i till mathematically valid and give a atifactory approximation to the problem, ceae to be phyically relevant becaue the walllocalized convection given by (4.26) (4.28) become preferred. A typical wall-localized olution at Pr=1. fore =1 4 and = 2 i hown in figure 2(e,f ). In thi cae, the aymptotic analyi give R c =63.6,m c =9.56 and σ c =.2, while the numerical analyi yield R c =62.6,m c =1 and σ c =.24. Moreover, the agreement between the aymptotic and numerical analyi i expectedly better for maller value of E. For example, at E =1 5 with =2 and Pr=.1, the aymptotic expreion (4.21) (4.22) yield R c =15.,m c =2,σ c =.15 while the numerical analyi give R c =15.4,m c =2,σ c =.15. Our analyi reveal that there exit three different regime of convection in a rotating cylinder at a given mall E. The firt i the regime of inertial convection at very mall Pr, i.e. Pr/E O(1). A diplayed in the left portion of figure 1 for Pr 1 4, a well a indicated by the aymptotic expreion (4.21), the critical Rayleigh number R c i nearly independent of Pr with R c = O(E 1/2 ). In thi regime, the temperature Θ i marked by approximately the ame phae a that of the vertical flow u z in the convective heat tranport. But thi ituation i changed when Pr/E>O(1), which i hown in the middle portion of figure 1 for 1 4 <Pr<1 1. Although the inertial effect i till predominant, the temperature Θ i marked by a large phae hift comparing to the vertical flow u z and, conequently, the R c increae ignificantly with increaing Pr. In thi moderate-pr regime, the aymptotic expreion (4.21) indicate that there exit no imple aymptotic relation between R c and E or between R c and Pr. It i alo in thi regime that convection pattern i not only highly complex but alo trongly dependent on Pr. Thi complicated behaviour tem from the exitence of a large manifold of two-dimenional inertial mode that may be excited by convective intabilitie in thi regime, indicating how difficult any numerical attempt to pinpoint it coherent tructure will be. For larger value of Pr, which i hown in the right portion of figure 1 for Pr>1 1, the trong vicou effect reulting from the walllocalized boundary-layer flow i dynamically dominant, leading to the relatively imple aymptotic dependence. Our calculation for different value of E ugget that thi picture i qualitatively correct for = O(1) at any given mall Ekman number E. 5. Summary and remark The problem of convection in a rotating circular cylinder heated from below and rotating about a vertical axi ha been tudied a a paradigm for undertanding the general dynamic of rotating convection. We have undertaken both numerical and
11 72 K. Zhang and X. Liao analytical tudie of convective intabilitie under experimental boundary condition, howing that the intabilitie are either of inertial-wave type or of wall-localized type. In the aymptotic tudy, the primary aumption i that the convective flow at mall or moderate Pr can be repreented by an inertial-wave mode. In the numerical tudy, we have olved the partial differential equation governing the convection problem. Agreement between the analytical and numerical olution ha been hown to be atifactory. There are phyical and mathematical conequence of our aymptotic analyi for convection in rotating cylinder. Phyically peaking, when a cylinder rotate rapidly (E 1), the fundamental dynamic of fluid motion can be intuitively illutrated by the Proudman Taylor theorem tating that infiniteimal teady motion in rotating invicid fluid are two-dimenional with repect to the direction of ẑ. It follow that u z / z = and that, upon applying the boundary condition on the top and bottom of a cylinder, we mut conclude that u z = within the cylinder. In other word, convection cannot take place: the effect of rotation trongly contrain and tabilize the ytem. The contraint mut be broken o that convection can occur. Our aymptotic analyi how that it can be broken either by inertial effect for fluid with mall Prandtl number or by vicou effect for fluid with large Prandtl number in connection with the boundary-layer-type convection. Mathematically peaking, the baic aumption that the convection velocity at leading order for mall or moderate Pr can be repreented by an inertial-wave mode implie that the preent tability analyi i quite different from the uual one uch a that for the Rayleigh Bénard problem. In our tability analyi, it i not a quetion of imply finding a wavenumber m that minimize the Rayleigh number and then determine the mot untable convection mode. It i a quetion of finding the two-dimenional tructure of a leading-order flow that give rie to the minimum Rayleigh number. Thi explain why it i o difficult to elucidate the parametric dependence of the problem numerically. Moreover, ince the aymptotic formula can be readily evaluated, they would provide a helpful guidance for the experimental tudie of convection in rotating circular cylinder. Thi tudy repreent the firt aymptotic analyi for convective intabilitie in a rotating cylinder with experimental boundary condition for mall and moderate Prandtl number. It i ignificant to note that the two et of aymptotic formula given in 4 cover the whole parameter regime of the problem: all value of the apect ratio and the Prandtl number Pr at any fixed mall Ekman number E. Epecially, it i omewhat urpriing that the eemingly incomprehenible and highly complicated convection revealed numerically, for example, by Goldtein et al. (1993, 1994) can be atifactorily decribed by everal relatively imple analytical formula. Keke Zhang i upported by UK PPARC/NERC grant and Xinhao Liao i upported by NSFC/16333, STCSM/8XD1452 and CAS grant. The numerical computation i upported by SSC. We would like to thank Profeor F. H. Bue for hi helpful dicuion. REFERENCES Aurnou,J.M.,Andreadi,S.,Zhu,L.&OlonP.L.23. Experiment on convection in Earth core tangent cylinder. Earth Planet. Sci. Lett. 212, Bue,F.H.&SimitevR.24. Inertial convection in rotating fluid phere. J. Fluid Mech. 498, Daniel, P. G The effect of centrifugal acceleration on axiymmetric convection in a hallow rotating cylinder or annulu. J. Fluid Mech. 99,
12 Convection in rotating circular cylinder 73 Goldtein, H. F., Knobloch, E., Mercader, I. & Net, M Convection in a rotating cylinder. Part 1. Linear theory for moderate Prandtl number. J. Fluid Mech. 248, Goldtein, H. F., Knobloch, E., Mercader, I. & Net, M Convection in a rotating cylinder. Part 2. Linear theory for mall Prandtl number. J. Fluid Mech. 248, Herrmann, J. & Bue, F. H Aymptotic theory of wall-localized convection in a rotating fluid layer. J. Fluid Mech. 255, Kuo,E.Y.&Cro,M.C Travelling-wave wall tate in rotating Rayleigh Benard convection. Phy. Rev. E 47, Liao, X., Zhang, K. & Chang, Y. 26. On boundary-layer convection in a rotating fluid layer. J. Fluid Mech. 549, Marqué, F On boundary condition for velocity potential in confined flow: application to Couette flow. Phy. Fluid A 2, Net, M., Carcia, F. Sanchez, J. On the onet of low-prandtl-number convection in rotating pherical hell: non-lip boundary condition. J. Fluid Mech. 61, Plaut, E. 23. Nonlinear dynamic of travelling Rayleigh Bénard convection: effect of the boundary condition and of the topology. Phy. Rev. E 67, 4633, Young, R. M. B. & Read, P. L. 28. Breeding and predictability in the baroclinic rotating annulu uing a perfect model. Nonlinear Proce. Geophy. 15, Zhang, K On coupling between the Poincaré equation and the heat equation: non-lip boundary condition. J. Fluid Mech. 284, Zhang, K. & Liao, X. 28. On the initial value problem in a rotating circular cylinder. J. Fluid Mech. 61, Zhang, K., Liao, X. & Bue, F. H. 27. Aymptotic theory of inertial convection in a rotating cylinder. J. Fluid Mech. 575, Zhong, F., Ecke, R. E. & Steinberg, V Aymmetric mode and the tranition to vortex tructure in rotating Rayleigh Benard convection. Phy. Rev. Lett. 67,
PRESSURE WORK EFFECTS IN UNSTEADY CONVECTIVELY DRIVEN FLOW ALONG A VERTICAL PLATE
Proceeding of 3ICCHMT 3 rd International Conference on Computational Heat and Ma Tranfer May 6 3, 3, Banff, CANADA Paper Number 87 PRESSURE WORK EFFECTS IN UNSTEADY CONVECTIVELY DRIVEN FLOW ALONG A VERTICAL
More informationStudy of a Freely Falling Ellipse with a Variety of Aspect Ratios and Initial Angles
Study of a Freely Falling Ellipe with a Variety of Apect Ratio and Initial Angle Dedy Zulhidayat Noor*, Ming-Jyh Chern*, Tzyy-Leng Horng** *Department of Mechanical Engineering, National Taiwan Univerity
More informationTo appear in International Journal of Numerical Methods in Fluids in Stability analysis of numerical interface conditions in uid-structure therm
To appear in International Journal of Numerical Method in Fluid in 997. Stability analyi of numerical interface condition in uid-tructure thermal analyi M. B. Gile Oxford Univerity Computing Laboratory
More informationBogoliubov Transformation in Classical Mechanics
Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How
More informationLecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)
Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.
More informationFluid-structure coupling analysis and simulation of viscosity effect. on Coriolis mass flowmeter
APCOM & ISCM 11-14 th December, 2013, Singapore luid-tructure coupling analyi and imulation of vicoity effect on Corioli ma flowmeter *Luo Rongmo, and Wu Jian National Metrology Centre, A*STAR, 1 Science
More informationJump condition at the boundary between a porous catalyst and a homogeneous fluid
From the SelectedWork of Francico J. Valde-Parada 2005 Jump condition at the boundary between a porou catalyt and a homogeneou fluid Francico J. Valde-Parada J. Alberto Ochoa-Tapia Available at: http://work.bepre.com/francico_j_valde_parada/12/
More informationNumerical Simulations of Coriolis Flow Meters for Low Reynolds Number Flows
MAPAN - Journal Numerical of Metrology Simulation Society of of Corioli India, Vol. Flow 26, Meter No. 3, 2011; for Low pp. Reynold 225-235 Number Flow ORIGINAL ARTICLE Numerical Simulation of Corioli
More information696 Fu Jing-Li et al Vol. 12 form in generalized coordinate Q ffiq dt = 0 ( = 1; ;n): (3) For nonholonomic ytem, ffiq are not independent of
Vol 12 No 7, July 2003 cfl 2003 Chin. Phy. Soc. 1009-1963/2003/12(07)/0695-05 Chinee Phyic and IOP Publihing Ltd Lie ymmetrie and conerved quantitie of controllable nonholonomic dynamical ytem Fu Jing-Li(ΛΠ±)
More informationMarch 18, 2014 Academic Year 2013/14
POLITONG - SHANGHAI BASIC AUTOMATIC CONTROL Exam grade March 8, 4 Academic Year 3/4 NAME (Pinyin/Italian)... STUDENT ID Ue only thee page (including the back) for anwer. Do not ue additional heet. Ue of
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction
More informationSingular perturbation theory
Singular perturbation theory Marc R. Rouel June 21, 2004 1 Introduction When we apply the teady-tate approximation (SSA) in chemical kinetic, we typically argue that ome of the intermediate are highly
More informationAMS 212B Perturbation Methods Lecture 20 Part 1 Copyright by Hongyun Wang, UCSC. is the kinematic viscosity and ˆp = p ρ 0
Lecture Part 1 Copyright by Hongyun Wang, UCSC Prandtl boundary layer Navier-Stoke equation: Conervation of ma: ρ t + ( ρ u) = Balance of momentum: u ρ t + u = p+ µδ u + ( λ + µ ) u where µ i the firt
More informationFinal Comprehensive Exam Physical Mechanics Friday December 15, Total 100 Points Time to complete the test: 120 minutes
Final Comprehenive Exam Phyical Mechanic Friday December 15, 000 Total 100 Point Time to complete the tet: 10 minute Pleae Read the Quetion Carefully and Be Sure to Anwer All Part! In cae that you have
More informationON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang
Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang
More informationOnline supplementary information
Electronic Supplementary Material (ESI) for Soft Matter. Thi journal i The Royal Society of Chemitry 15 Online upplementary information Governing Equation For the vicou flow, we aume that the liquid thickne
More informationSupplementary Figures
Supplementary Figure Supplementary Figure S1: Extraction of the SOF. The tandard deviation of meaured V xy at aturated tate (between 2.4 ka/m and 12 ka/m), V 2 d Vxy( H, j, hm ) Vxy( H, j, hm ) 2. The
More informationMAE 101A. Homework 3 Solutions 2/5/2018
MAE 101A Homework 3 Solution /5/018 Munon 3.6: What preure gradient along the treamline, /d, i required to accelerate water upward in a vertical pipe at a rate of 30 ft/? What i the anwer if the flow i
More informationOMNÈS REPRESENTATIONS WITH INELASTIC EFFECTS FOR HADRONIC FORM FACTORS
FIELD THEORY, GRAVITATION AND PARTICLE PHYSICS OMNÈS REPRESENTATIONS WITH INELASTIC EFFECTS FOR HADRONIC FORM FACTORS IRINEL CAPRINI National Intitute of Phyic and Nuclear Engineering POB MG 6, Bucharet,
More informationV = 4 3 πr3. d dt V = d ( 4 dv dt. = 4 3 π d dt r3 dv π 3r2 dv. dt = 4πr 2 dr
0.1 Related Rate In many phyical ituation we have a relationhip between multiple quantitie, and we know the rate at which one of the quantitie i changing. Oftentime we can ue thi relationhip a a convenient
More informationMath 273 Solutions to Review Problems for Exam 1
Math 7 Solution to Review Problem for Exam True or Fale? Circle ONE anwer for each Hint: For effective tudy, explain why if true and give a counterexample if fale (a) T or F : If a b and b c, then a c
More informationA Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems
A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationControl Systems Analysis and Design by the Root-Locus Method
6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If
More informationin a circular cylindrical cavity K. Kakazu Department of Physics, University of the Ryukyus, Okinawa , Japan Y. S. Kim
Quantization of electromagnetic eld in a circular cylindrical cavity K. Kakazu Department of Phyic, Univerity of the Ryukyu, Okinawa 903-0, Japan Y. S. Kim Department of Phyic, Univerity of Maryland, College
More informationTHE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER
Proceeding of IMAC XXXI Conference & Expoition on Structural Dynamic February -4 Garden Grove CA USA THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER Yung-Sheng Hu Neil S Ferguon
More informationCHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS
Copyright 22 IFAC 5th Triennial World Congre, Barcelona, Spain CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Tritan Pérez Graham C. Goodwin Maria M. Serón Department of Electrical
More informationConvergence criteria and optimization techniques for beam moments
Pure Appl. Opt. 7 (1998) 1221 1230. Printed in the UK PII: S0963-9659(98)90684-5 Convergence criteria and optimization technique for beam moment G Gbur and P S Carney Department of Phyic and Atronomy and
More informationRecent progress in fire-structure analysis
EJSE Special Iue: Selected Key Note paper from MDCMS 1 1t International Conference on Modern Deign, Contruction and Maintenance of Structure - Hanoi, Vietnam, December 2007 Recent progre in fire-tructure
More informationTHE IDENTIFICATION OF THE OPERATING REGIMES OF THE CONTROLLERS BY THE HELP OF THE PHASE TRAJECTORY
Mariu M. B LA Aurel Vlaicu Univerity of Arad, Engineering Faculty Bd. Revolu iei nr. 77, 3030, Arad, Romania, E-mail: mariu.bala@ieee.org THE IDENTIFICATION OF THE OPERATING REGIMES OF THE CONTROLLERS
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world leading publiher of Open Acce book Built by cientit, for cientit 3,5 8,.7 M Open acce book available International author and editor Download Our author are among the 5 Countrie
More informationMolecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions
Original Paper orma, 5, 9 7, Molecular Dynamic Simulation of Nonequilibrium Effect ociated with Thermally ctivated Exothermic Reaction Jerzy GORECKI and Joanna Natalia GORECK Intitute of Phyical Chemitry,
More informationOn the Isentropic Forchheimer s Sound Waves Propagation in a Cylindrical Tube Filled with a Porous Media
5th WSEAS Int. Conf. on FLUID MECHANICS (FLUIDS') Acapulco, Mexico, January 5-7, On the Ientropic Forchheimer Sound Wave Propagation in a Cylindrical Tube Filled with a Porou Media H. M. Dwairi Civil Engineering
More informationSIMPLIFIED MODEL FOR EPICYCLIC GEAR INERTIAL CHARACTERISTICS
UNIVERSITY OF PITESTI SCIENTIFIC BULLETIN FACULTY OF ECHANICS AND TECHNOLOGY AUTOOTIVE erie, year XVII, no. ( 3 ) SIPLIFIED ODEL FOR EPICYCLIC GEAR INERTIAL CHARACTERISTICS Ciobotaru, Ticuşor *, Feraru,
More informationThe continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.
1 I. CONTINUOUS TIME RANDOM WALK The continuou time random walk (CTRW) wa introduced by Montroll and Wei 1. Unlike dicrete time random walk treated o far, in the CTRW the number of jump n made by the walker
More informationCHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3
More informationChapter 1 Basic Description of Laser Diode Dynamics by Spatially Averaged Rate Equations: Conditions of Validity
Chapter 1 Baic Decription of Laer Diode Dynamic by Spatially Averaged Rate Equation: Condition of Validity A laer diode i a device in which an electric current input i converted to an output of photon.
More informationAdvanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment
Journal of Multidiciplinary Engineering Science and Technology (JMEST) ISSN: 59- Vol. Iue, January - 5 Advanced D-Partitioning Analyi and it Comparion with the haritonov Theorem Aement amen M. Yanev Profeor,
More informationIII.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES
III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SBSTANCES. Work purpoe The analyi of the behaviour of a ferroelectric ubtance placed in an eternal electric field; the dependence of the electrical polariation
More informationOnset of Stern Type Thermohaline Convection with Variable Gravity Using Positive Operator
International Journal of Pure and Applied Mathematical Science. ISSN 97-988 Volume 9 Number (6 pp. 7-5 Reearch India Publication http://www.ripublication.com Onet of Stern ype hermohaline Convection with
More informationMODERN CONTROL SYSTEMS
MODERN CONTROL SYSTEMS Lecture 1 Root Locu Emam Fathy Department of Electrical and Control Engineering email: emfmz@aat.edu http://www.aat.edu/cv.php?dip_unit=346&er=68525 1 Introduction What i root locu?
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. VIII Decoupling Control - M. Fikar
DECOUPLING CONTROL M. Fikar Department of Proce Control, Faculty of Chemical and Food Technology, Slovak Univerity of Technology in Bratilava, Radlinkého 9, SK-812 37 Bratilava, Slovakia Keyword: Decoupling:
More informationSimilarity of the Temperature Profile formed by Fluid Flow along a Wall
Similarity of the emperature Profile formed by Fluid Flow along a Wall David Weyburne 1 AFRL/RYDH 2241 Avionic Circle WPAFB OH 45433 ABSRAC A new approach to the tudy of imilarity of temperature profile
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationBUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL
J. Autral. Math. Soc. Ser. B 4(999), 332 349 BUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL J. LEE and J.-M. VANDEN-BROECK 2 (Received 22 April 995; revied 23 April 996)
More informationSpot-on: Safe Fuel/Air Compression
Spot-on: Safe Fuel/Air Compreion Problem preented by Robert Hart and Kevin Hughe Veeder-Root Participant: Jeffrey Bank Joeph Fehribach Alitair Fitt John Ockendon Colin Pleae Don Schwendeman Burt Tilley
More informationPHYS 110B - HW #6 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased
PHYS B - HW #6 Spring 4, Solution by David Pace Any referenced equation are from Griffith Problem tatement are paraphraed. Problem. from Griffith Show that the following, A µo ɛ o A V + A ρ ɛ o Eq..4 A
More informationEFFECT OF VOLUMETRIC HEAT LOSS ON TRIPLE-FLAME PROPAGATION
Proceeding of the Combution Intitute, Volume 29, 2002/pp. 1559 1564 EFFECT OF VOLUMETRIC HEAT LOSS ON TRIPLE-FLAME PROPAGATION R. DAOU, J. DAOU and J. DOLD Department of Mathematic, UMIST, Mancheter M60
More informationonline learning Unit Workbook 4 RLC Transients
online learning Pearon BTC Higher National in lectrical and lectronic ngineering (QCF) Unit 5: lectrical & lectronic Principle Unit Workbook 4 in a erie of 4 for thi unit Learning Outcome: RLC Tranient
More informationMath Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK
ppendix 5 Scientific Notation It i difficult to work with very large or very mall number when they are written in common decimal notation. Uually it i poible to accommodate uch number by changing the SI
More informationANALYTICAL BEARING MODEL FOR ANALYSIS OF INNER LOAD DISTRIBUTION AND ESTIMATION OF OPERATIONAL LUBRICATION REGIME
58 th ICMD 017 6-8 September 017, Prague, Czech Republic ANALYTICAL BEARING MODEL FOR ANALYSIS OF INNER LOAD DISTRIBUTION AND ESTIMATION OF OPERATIONAL LUBRICATION REGIME Jakub CHMELAŘ 1, Votěch DYNYBYL
More informationThe Hassenpflug Matrix Tensor Notation
The Haenpflug Matrix Tenor Notation D.N.J. El Dept of Mech Mechatron Eng Univ of Stellenboch, South Africa e-mail: dnjel@un.ac.za 2009/09/01 Abtract Thi i a ample document to illutrate the typeetting of
More informationA Simplified Methodology for the Synthesis of Adaptive Flight Control Systems
A Simplified Methodology for the Synthei of Adaptive Flight Control Sytem J.ROUSHANIAN, F.NADJAFI Department of Mechanical Engineering KNT Univerity of Technology 3Mirdamad St. Tehran IRAN Abtract- A implified
More informationBernoulli s equation may be developed as a special form of the momentum or energy equation.
BERNOULLI S EQUATION Bernoulli equation may be developed a a pecial form of the momentum or energy equation. Here, we will develop it a pecial cae of momentum equation. Conider a teady incompreible flow
More informationSMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD
SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD S.P. Teeuwen, I. Erlich U. Bachmann Univerity of Duiburg, Germany Department of Electrical Power Sytem
More informationLecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004
18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem
More informationIsentropic Sound Waves Propagation in a Tube Filled with a Porous Media
INTERNATIONAL JOURNAL OF ECHANICS Ientropic Sound Wave Propagation in a Tube Filled with a Porou edia H.. Duwairi Abtract A rigid frame, cylindrical capillary theory of ound propagation in porou media
More informationPhysics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2014
Phyic 7 Graduate Quantum Mechanic Solution to inal Eam all 0 Each quetion i worth 5 point with point for each part marked eparately Some poibly ueful formula appear at the end of the tet In four dimenion
More informationCritical Height of Slopes in Homogeneous Soil: the Variational Solution
Critical Height of Slope in Homogeneou Soil: the Variational Solution Chen, Rong State Key Laboratory of Coatal and Offhore Engineering & Intitute of Geotechnical Engineering, Dalian Univerity of Technology,
More informationS_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS
S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS by Michelle Gretzinger, Daniel Zyngier and Thoma Marlin INTRODUCTION One of the challenge to the engineer learning proce control i relating theoretical
More informationEFFECT ON PERSISTENCE OF INTRA-SPECIFIC COMPETITION IN COMPETITION MODELS
Electronic Journal of Differential Equation, Vol. 2007(2007, No. 25, pp. 0. ISSN: 072-669. URL: http://ejde.math.txtate.edu or http://ejde.math.unt.edu ftp ejde.math.txtate.edu (login: ftp EFFECT ON PERSISTENCE
More informationFiltered Leapfrog Time Integration with Enhanced Stability Properties
Journal of Applied Mathematic and Phyic, 6, 4, 4-7 Publihed Online July 6 in SciRe. http://www.cirp.org/ournal/amp http://dx.doi.org/.46/amp.6.474 Filtered Leapfrog Time Integration with Enhanced Stability
More informationA novel protocol for linearization of the Poisson-Boltzmann equation
Ann. Univ. Sofia, Fac. Chem. Pharm. 16 (14) 59-64 [arxiv 141.118] A novel protocol for linearization of the Poion-Boltzmann equation Roumen Tekov Department of Phyical Chemitry, Univerity of Sofia, 1164
More informationAt the end of this lesson, the students should be able to understand:
Intructional Objective: At the end of thi leon, the tudent hould be able to undertand: Baic failure mechanim of riveted joint. Concept of deign of a riveted joint. 1. Strength of riveted joint: Strength
More informationIntroduction to Laplace Transform Techniques in Circuit Analysis
Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found
More informationUSING NONLINEAR CONTROL ALGORITHMS TO IMPROVE THE QUALITY OF SHAKING TABLE TESTS
October 12-17, 28, Beijing, China USING NONLINEAR CONTR ALGORITHMS TO IMPROVE THE QUALITY OF SHAKING TABLE TESTS T.Y. Yang 1 and A. Schellenberg 2 1 Pot Doctoral Scholar, Dept. of Civil and Env. Eng.,
More informationModeling of seed magnetic island formation
EUROFUSION WPMST1-PR(16) 15656 IG Miron et al. Modeling of eed magnetic iland formation Preprint of Paper to be ubmitted for publication in 43rd European Phyical Society Conference on Plama Phyic (EPS)
More informationPHYS 110B - HW #2 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased
PHYS 11B - HW # Spring 4, Solution by David Pace Any referenced equation are from Griffith Problem tatement are paraphraed [1.] Problem 7. from Griffith A capacitor capacitance, C i charged to potential
More informationNCAAPMT Calculus Challenge Challenge #3 Due: October 26, 2011
NCAAPMT Calculu Challenge 011 01 Challenge #3 Due: October 6, 011 A Model of Traffic Flow Everyone ha at ome time been on a multi-lane highway and encountered road contruction that required the traffic
More informationSERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lectures 41-48)
Chapter 5 SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lecture 41-48) 5.1 Introduction Power ytem hould enure good quality of electric power upply, which mean voltage and current waveform hould
More informationCHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL
98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i
More informationTHE THERMOELASTIC SQUARE
HE HERMOELASIC SQUARE A mnemonic for remembering thermodynamic identitie he tate of a material i the collection of variable uch a tre, train, temperature, entropy. A variable i a tate variable if it integral
More informationAnnex-A: RTTOV9 Cloud validation
RTTOV-91 Science and Validation Plan Annex-A: RTTOV9 Cloud validation Author O Embury C J Merchant The Univerity of Edinburgh Intitute for Atmo. & Environ. Science Crew Building King Building Edinburgh
More informationarxiv: v2 [nucl-th] 3 May 2018
DAMTP-207-44 An Alpha Particle Model for Carbon-2 J. I. Rawlinon arxiv:72.05658v2 [nucl-th] 3 May 208 Department of Applied Mathematic and Theoretical Phyic, Univerity of Cambridge, Wilberforce Road, Cambridge
More information84 ZHANG Jing-Shang Vol. 39 of which would emit 5 He rather than 3 He. 5 He i untable and eparated into n + pontaneouly, which can alo be treated a if
Commun. Theor. Phy. (Beijing, China) 39 (003) pp. 83{88 c International Academic Publiher Vol. 39, No. 1, January 15, 003 Theoretical Analyi of Neutron Double-Dierential Cro Section of n+ 11 B at 14. MeV
More information/University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2009
Lecture 0 /6/09 /Univerity of Wahington Department of Chemitry Chemitry 453 Winter Quarter 009. Wave Function and Molecule Can quantum mechanic explain the tructure of molecule by determining wave function
More informationStochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions
Stochatic Optimization with Inequality Contraint Uing Simultaneou Perturbation and Penalty Function I-Jeng Wang* and Jame C. Spall** The John Hopkin Univerity Applied Phyic Laboratory 11100 John Hopkin
More informationDEVELOPMENT OF A STRUCTURED THERMOCLINE THERMAL ENERGY STORAGE SYSTEM
DEVELOPMENT OF A STRUCTURED THERMOCLINE THERMAL ENERGY STORAGE SYSTEM Brad M. Brown Matt N. Straer R. Paneer Selvam Univerity of Arkana Department of Civil Engineering 4190 Bell Engineering Center Fayetteville,
More informationNearshore Sediment Transport Modeling: Collaborative Studies with the U. S. Naval Research Laboratory
Nearhore Sediment Tranport Modeling: Collaborative Studie with the U. S. Naval Reearch Laboratory Donald N. Slinn Department of Civil and Coatal Engineering, Univerity of Florida Gaineville, FL 32611-6590,
More informationModeling of Transport and Reaction in a Catalytic Bed Using a Catalyst Particle Model.
Excerpt from the Proceeding of the COMSOL Conference 2010 Boton Modeling of Tranport and Reaction in a Catalytic Bed Uing a Catalyt Particle Model. F. Allain *,1, A.G. Dixon 1 1 Worceter Polytechnic Intitute
More informationELECTROMAGNETIC FIELD IN THE PRESENCE OF A THREE-LAYERED SPHERICAL REGION
Progre In Electromagnetic Reearch, PIER 45, 103 11, 004 ELECTROMAGNETIC FIELD IN THE PRESENCE OF A THREE-LAYERED SPHERICAL REGION K. Li and S.-O. Park School of Engineering Information and Communication
More informationAutomatic Control Systems. Part III: Root Locus Technique
www.pdhcenter.com PDH Coure E40 www.pdhonline.org Automatic Control Sytem Part III: Root Locu Technique By Shih-Min Hu, Ph.D., P.E. Page of 30 www.pdhcenter.com PDH Coure E40 www.pdhonline.org VI. Root
More informationOn the Curv ature of Space ²
Gener al Relativity and Gravitation, Vol. 31, No. 1, 1999 On the Curv ature of Space ² By A. Friedman in Peterburg * With one gure. Received on 9. June 19 1. 1. In their well-known work on general comologic
More informationEmittance limitations due to collective effects for the TOTEM beams
LHC Project ote 45 June 0, 004 Elia.Metral@cern.ch Andre.Verdier@cern.ch Emittance limitation due to collective effect for the TOTEM beam E. Métral and A. Verdier, AB-ABP, CER Keyword: TOTEM, collective
More informationConduction Heat transfer: Unsteady state
Conduction Heat tranfer: Unteady tate Chapter Objective For olving the ituation that Where temperature do not change with poition. In a imple lab geometry where temperature vary alo with poition. Near
More informationGain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays
Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,
More informationNULL HELIX AND k-type NULL SLANT HELICES IN E 4 1
REVISTA DE LA UNIÓN MATEMÁTICA ARGENTINA Vol. 57, No. 1, 2016, Page 71 83 Publihed online: March 3, 2016 NULL HELIX AND k-type NULL SLANT HELICES IN E 4 1 JINHUA QIAN AND YOUNG HO KIM Abtract. We tudy
More informationControl of Delayed Integrating Processes Using Two Feedback Controllers R MS Approach
Proceeding of the 7th WSEAS International Conference on SYSTEM SCIENCE and SIMULATION in ENGINEERING (ICOSSSE '8) Control of Delayed Integrating Procee Uing Two Feedback Controller R MS Approach LIBOR
More informationDepartment of Aerospace Engineering and Engineering Mechanics, University of Texas, Austin,
Buckling mode of elatic thin film on elatic ubtrate Haixia Mei and Rui Huang * Department of Aeropace Engineering and Engineering Mechanic, Univerity of Texa, Autin, TX 78712 Jun Young Chung and Chritopher
More informationOn the Normal-Mode Frequency Spectrum of Kinetic Magnetohydrodynamics. J.J. Ramos. December, 2014
PSFC/JA-14-33 On the Normal-Mode Frequency Spectrum of Kinetic Magnetohydrodynamic J.J. Ramo December, 14 Plama Science and Fuion Center Maachuett Intitute of Technology Cambridge MA 139, U.S.A. Thi work
More informationCalculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient
Ŕ periodica polytechnica Mechanical Engineering 54/1 21 15 2 doi: 1.3311/pp.me.21-1.3 web: http:// www.pp.bme.hu/ me c Periodica Polytechnica 21 RESERCH RTICLE Calculation of the temperature of boundary
More informationStandard Guide for Conducting Ruggedness Tests 1
Deignation: E 69 89 (Reapproved 996) Standard Guide for Conducting Ruggedne Tet AMERICA SOCIETY FOR TESTIG AD MATERIALS 00 Barr Harbor Dr., Wet Conhohocken, PA 948 Reprinted from the Annual Book of ASTM
More informationPI control system design for Electromagnetic Molding Machine based on Linear Programing
PI control ytem deign for Electromagnetic Molding Machine baed on Linear Programing Takayuki Ihizaki, Kenji Kahima, Jun-ichi Imura*, Atuhi Katoh and Hirohi Morita** Abtract In thi paper, we deign a PI
More informationA SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES. Sanghyun Cho
A SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES Sanghyun Cho Abtract. We prove a implified verion of the Nah-Moer implicit function theorem in weighted Banach pace. We relax the
More informationTHE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY STABILIZING CONTROLLERS. Tatsuya Hoshikawa, Kou Yamada and Yuko Tatsumi
International Journal of Innovative Computing, Information Control ICIC International c 206 ISSN 349-498 Volume 2, Number 2, April 206 pp. 357 370 THE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY
More information1. Basic introduction to electromagnetic field. wave properties and particulate properties.
Lecture Baic Radiometric Quantitie. The Beer-Bouguer-Lambert law. Concept of extinction cattering plu aborption and emiion. Schwarzchild equation. Objective:. Baic introduction to electromagnetic field:
More informationA BATCH-ARRIVAL QUEUE WITH MULTIPLE SERVERS AND FUZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH
Mathematical and Computational Application Vol. 11 No. pp. 181-191 006. Aociation for Scientific Reearch A BATCH-ARRIVA QEE WITH MTIPE SERVERS AND FZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH Jau-Chuan
More informationRelationship between surface velocity divergence and gas transfer in open-channel flows with submerged simulated vegetation
IOP Conference Serie: Earth and Environmental Science PAPER OPEN ACCESS Relationhip between urface velocity divergence and ga tranfer in open-channel flow with ubmerged imulated vegetation To cite thi
More informationCake ltration analysis the eect of the relationship between the pore liquid pressure and the cake compressive stress
Chemical Engineering Science 56 (21) 5361 5369 www.elevier.com/locate/ce Cake ltration analyi the eect of the relationhip between the pore liquid preure and the cake compreive tre C. Tien, S. K. Teoh,
More informationEuler-Bernoulli Beams
Euler-Bernoulli Beam The Euler-Bernoulli beam theory wa etablihed around 750 with contribution from Leonard Euler and Daniel Bernoulli. Bernoulli provided an expreion for the train energy in beam bending,
More information