TOPICAL PROBLEMS OF FLUID MECHANICS 245 ON NUMERICAL APPROXIMATION OF FLUID-STRUCTURE INTERACTIONS OF AIR FLOW WITH A MODEL OF VOCAL FOLDS
|
|
- Abigayle Parrish
- 5 years ago
- Views:
Transcription
1 TOPICAL PROBLEMS OF FLUID MECHANICS 245 ON NUMERICAL APPROXIMATION OF FLUID-STRUCTURE INTERACTIONS OF AIR FLOW WITH A MODEL OF VOCAL FOLDS J. Valášek 1, J. Horáček 2, P. Sváček 1 1 Department of Technical Mathematic, Faculty of Mechanical Engineering, Czech Technical Univerity in Prague, Karlovo nám. 13, Praha 2, , Czech Republic 2 Intitute of Thermomechanic, Czech Academy of Science, Dolejškova 5, Praha 8 Abtract Thi paper deal with flow driven vibration of an elatic body. Our goal i to develop and mathematically decribe a implified model of the human vocal fold. The developed numerical cheme for vicou incompreible fluid flow in ALE formulation and the elatic body are implemented by two olver, pecific for each domain. The tudied problem i coupled by Dirichlet-Neumann boundary condition. Both olver are baed on the finite element method. Particularly, for the fluid model the crogrid element are ued. Numerical reult focu on the verification of the developed program. Keyword: finite element method, 2D Navier-Stoke equation, vocal fold, aeroelaticity 1 Introduction The flow driven vibration of elatic bodie i a problem invetigated and olved in many technical application. Thi paper deal biomechanic of human vocal fold, ee e.g. [11]. Human voice i one of baic human being characteritic and it play an important role in the quality of a human life. The air flow from lung excite the vibration of the human fold which etting i influenced by the human mucle and caue the human voice production. Therefore the attention i further devoted for better undertanding of thi complex proce. Thi problem wa tudied in many paper, for example ee [6], [10]. Neverthele, a uitable benchmark for fluid-tructure interaction (FSI) problem i till miing. Here, we focu on the development of numerical method, careful verification of of all it component, and on verification of the olution of the coupled problem. The verification i baed on a benchmark problem introduced here, where only low Reynold number i involved and thu the convergence of the numerical method can be eaily tudied. The numerical method i baed on finite element method (FEM). For the fluid olver the cro-grid P1 element and for the elatic part P1 element are implemented. The ALE method i involved for handling with the time-dependent domain. For the coupled problem the emi-implicit cheme i ued. Figure 1: Schema of vocal fold model with boundarie marked before and after deformation.
2 246 Prague, February 11-13, Mathematical model For the ake of implicity our tudy i retricted to a 2D model problem. A cheme of the ued model i hown in Figure 1, where Ω ref denote the reference repreentation of the tructure, Ωf ref i the domain occupied by the fluid at the time intant t = 0 (alo reference) and Γ Wref = Γ W0 i the common interface. At intant time t the fluid domain Ω f ref turn to Ωf t. The deformated domain Ω = Ω ref at time t i handled in reference coordinate. Elatic tructure. The vocal fold model i conidered a the 2D compliant tructure with the aumption of mall diplacement. Uing the generalized Hook law, the Cauchy tre tenor τij i given by, ee for example [1], τij = C ijkl e kl, (1) where C ijkl denote the fourth order tenor of elatic coefficient and e jk i the train tenor given by e jk = 1 ( uj + u ) k. (2) 2 x k x j Under the aumption of the iotropic body equation (1) can be written in the form τ ij = λ div u + 2µ e ij, (3) where u denote the diplacement, λ, µ are Lame coefficient dependent on the Young modulu of elaticity E and the Poion ratio σ a λ = E σ (1 + σ )(1 2σ ), E µ = 2(1 + σ ). The deformation of the tructure i given by the equation of motion ρ 2 u t 2 τ ij (u) x j = f in Ω (0, T), (4) where the vector f i the volume denity of an acting force and ρ i the tructure denity. Equation (4) i upplied with the following initial and the boundary condition a) u(x, 0) = u 0 (x) for x Ω, u b) t (x, 0) = u 1(x) for x Ω, (5) c) u(x, t) = u Dir (x, t) for x Γ Dir, t (0, T), d) τ ij(x, t) n j(x) = q i (x, t), for x Γ W t, t (0, T), where the boundarie Γ Wt, Γ Dir are hown in Figure 1. ALE method. During the interaction the vocal fold change it hape, thu the fluid domain i alo deformed. Therefore we ue o called Arbitrary-Lagrangian method (ALE) to enable the flow imulation in the time-dependent domain. ALE method i the generalization of the Eulerian and the Lagrangian method. The baic aumption i that there exit difeomorphim A t : Ω f ref Ω f t, e.g. X Ω f ref x = A t (X) Ω f t. In addition mapping A t mut meet the condition A t C(Ω f ref ), A t ( Ω f ref t ) = Ωf t, t (0, T). (6) Uing he ALE mapping the domain velocity of the deformation i defined by w D (x, t) = ŵ D (A 1 t (x), t), t (0, T), x Ω f t, (7) where ŵ D i quantity defined on Ω f ref a ŵ D(X, t) = t A t(x), derivative atifie t (0, T), X Ω f ref. The ALE D A f f(x, t) = Dt t (A t(x), t) = f t (x, t) + w D(x, t) f(x, t), (8) where w D (x, t) i given by (7). Futher information i poible to find for example in the article [7].
3 TOPICAL PROBLEMS OF FLUID MECHANICS 247 Figure 2: The computational domain Ω f t at the time intant t and it boundarie. Flow model. The domain Ω f t at arbitrary time t i hown in Figure 2. The Γ f In i the inlet part of the boundary, Γ f Out i the outlet part of the boundary, Γf Dir repreent olid wall of the glottal channel urface and the only time dependent boundary Γ Wt repreent the interface on the vibrating vocal fold. The motion of the vicou incompreible fluid in the time dependent domain Ω f t i modelled by the Navier-Stoke equation in the ALE form D A v Dt + ((v w D ) )v ν f v + p = g f in Ω f t, div v = 0 in Ω f t, where v denote the fluid velocity and p i the kinematic preure, ν f i the kinematic vicoity of the fluid and g f i vector of acting force. The ytem of equation (9) i equipped with the boundary condition a) v(x, t) = 0 for x Γ f Dir, t (0, T), b) v(x, t) = v Dir (x, t) for x Γ f In, t (0, T), (10) c) (p(x, t) p ref ) n ν f v (x, t) = 0, n for x Γf Out, t (0, T), d) v(x, t) = w D (x, t) for x Γ Wt. The boundary condition on the interface Γ Wt i decribed later. Coupling condition. Deformation and location of interface Γ Wt in time t i dependent on (previou) acting of aerodynamical and elatic force. On the other hand thee force are dependent on poition and velocity of thi interface. Therefore we peak about a coupled problem. The hape of the interface Γ Wt in time t i the reult of force equilibrium between fluid and tructure. It poition i given by the deformation u, o Γ Wt = { x R 2 x = X + u(x, t), X Γ Wref }. (11) On the other hand aerodynamic force work on tructure interface via preure and hear force and the equilibrium i given by the equation 2 τij(x) n j(x) = j=1 2 σ f ij (x) n j(x), i = 1, 2, x Γ Wt, X Γ Wref, (12) j=1 where σ f ij i tre tenor of fluid n j (X) denote the component of the outer normal of the tructure domain to the interface Γ Wt=0 and x i computed from the known deformation u a x = X+u(X, t). The equation (12) can be rewritten a the Neumann boundary condition (5 d) q i (X, t) = 2 σ f ij (x) n j(x), i = 1, 2, x Γ Wt, t (0, T). (13) j=1 (9)
4 248 Prague, February 11-13, Numerical approximation The preented mathematical model repreented by the partial differential equation (4), (9) i dicretized in pace by the finite element method and in time by the Newmark and the backward Euler method, repectively. The time interval (0, T) i divided into the equiditant partition with a contant time tep t. Elatic tructure. The olution u i ought in the pace V = V V, where V = {f W 1,2 (Ω ) f = 0 on Γ Dir }, and W k,p (Ω) denote the Sobolev pace. With the ue of Hook law (1) and the ymmetry of tre tenor equation (4) i reformulated in the weak ene ρ 2u t 2 Φ dx + C ijkl e kl(u)e ij(φ) dx = f Φ dx + q Φ ds. (14) Ω Ω Ω Futhermore, the pace v i approximated by the finite dimenional pace V h V with the dimenion 2N h. Thu the dicrete olution can be written a a linear combination of bai function Φ j V h, i.e. u h (x, t) = 2N h j=1 α j(t)φ j (x). Uing thi and α = (α i ) the equation (3) now ha the form Γ Neu M T α + K T α = b(t), (15) where the element of the matrice M = (b ij ), K = (k ij ) and b = (b i ) are given by m ij = ρ Φ j Φ i dx, k ij = C klmn e mn(φ j )e kl(φ i ) dx, Ω Ω (b(t)) i = b i (t) = f Φ i dx + q Φ i ds. (16) Ω In the practical computation the Lagrange fnite element of the firt order are ued baed on a triangulation, which give the firt order of accuracy in pace. Newmark method. For the time dicretization the Newmark method i ued for the olution of ordinary differential equation of the econd order in a prototype form y (t) = f(t, y(t), y (t)) for t (0, T), y(0) = y 0, y (0) = y 0. (17) Γ Neu It application lead to the numerical cheme y n+1 = y n + ty n + t 2 (βf n+1 + ( 1 2 β)f n ) (18) y n+1 = y n + t (γf n+1 + (1 γ)f n ), (19) where f n = f(t n, y n, y n), f n+1 = f(t n+1, y n+1, y n+1) and the cheme i of econd order accuracy in time for the choice of the parameter β = 1 4, γ = 1 2, ee for example [2]. In the cae of equation (15) the value of y n+1, f n+1 i replaced by y n+1 obtained from the ytem ( M + β( t) 2 K ) y n+1 = gn+1. (20) }{{} =A Sytem of equation (20) i eay to be olved for example by conjugate gradient method, becaue the matrix A i compoed of two poitive defined matrice K, M. ALE method. The ALE mapping i treated uing the elatic analogy, becaue ALE method in fact map Ω f ref to deformated Ωf t, for dicrete contruction of ALE mapping i ued afore mentioned tationary part of the numerical cheme for elatic body, which i implified to K T d = 0, (21)
5 TOPICAL PROBLEMS OF FLUID MECHANICS 249 where d tay for the diplacement of the vertice of the triangulation τ h at the time intant t n+1. A boundary condition it i precribed Dirichlet boundary condition, which equal deformation of common interface on Γ Wtn+1 and otherwie equal zero d i = u(a 1 t n+1 (X i ), t n+1 ), X i Γ Wtn+1. (22) The Lame coefficient, which are ued during aembling of the matrix K T artificial and baed on article [5], according to thi formula in (21), i choen λ + µ = 0, λ = C avgdiam τ h diam K, (23) where C > 0 i a uitable contant, diam K i the diameter of the triangle K τ h and avgdiam τ h i the average element diameter of the triangulation τ h. So Lame contant are recaled on every triangle with repect to it area. Thi choice enure better behaviour at the interface, where a refined meh i ued. Flow model. For the flow model we ue the invere procedure - firtly the equation (9) i dicretized in time and then in pace. For time dicretization we ued backward Euler method uitable for ordinary differential equation of firt order. So the ALE derivative i approximated by D A v Dt (x n+1, t n+1 ) vn+1 (x n+1 ) v n (x n ). (24) t Denoting v n (x n+1 ) = v n (x n ), where x n = A tn (A 1 t n+1 (x n+1 )). the cheme read v n+1 v n t + ((v n+1 w n+1 D ) )vn+1 ν f v n+1 + p n+1 = g f,n+1, (25) v n+1 = 0. Then we continue with the weak formulation, where the velocity olution v in time t n+1 i ought in the functional pace X = W 1,2 (Ω f t n+1 ) and q M = L 2 (Ω f t n+1 ). Further, the pace of the tet function i defined by X = {f W 1,2 (Ω f t n+1 ) f = 0 on Γ f Dir Γf In Γf W tn+1 } W 1,2 (Ω f t n+1 ). For better arrangement it i written v n+1 a v and Ω f = Ω f t n+1. The weak formulation in pace i acquired by multiplication of the firt equation (25) by Φ X, integration over the whole domain Ω f t n+1 and by uing the Green theorem, that reult in the following equation ( v v n ), ϕ + ((v w D ) )v, ϕ) t Ω f + ν f ( v, ϕ) Ω f (p, div ϕ) Ω f = Ω f (q, div v) Ω f = 0, = (g f, ϕ) Ω f (p ref, ϕ n) L 2 (Γ f Out ), (26) where (, ) Ω tay for the calar product of L 2 (Ω f ) or [L 2 (Ω f )] 2 pace. The FEM then approximate pace X and M by the finite dimenion pace X h and M h, o the olution v v h can be expreed a v h (x) = 2Nh vel j=1 β j ϕ j (x), p h (x) = 2N vel h +N p h j=2n vel h +1 γ j q j (x). (27) Now by uing relation (27) in equation (26) we get the ytem ( ( ) ( A(v h ) B β g B 0) T =, (28) γ 0)
6 250 Prague, February 11-13, 2015 Figure 3: The crogrid finite element for incompreible problem (P1-crogrid/P1). vertice repreent preure variable and arrow are component of velocity. Dot at where A(v h ) = 1 t M + C(v h ) + D. The element of the matrice M = (m ij), C = (c ij ), D = (d ij ) are given by m ij = (ϕ j, ϕ i ) Ω f, c ij = ((v h w D ) )ϕ j, ϕ i ) Ω f, d ij = ν f ( ϕ j, ϕ i ) Ω f, (29) b ij = ( q j, div ϕ i ) Ω f, g i = (g f, ϕ i ) Ω f (p ref, ϕ i n) L2 (Γ f Out ) + ( un t, ϕ i) Ω f. The ytem of equation (28) i non-linear. For it olution the linearization v h = vn i ued. For the olution of the ytem (28) the mathematical library UMFPACK i employed, ee [3]. One of the apect of the FEM i that pace X h, M h cannot be choen randomly, but they mut atify the well-known Babuka-Brezzi condition, ee [4]. In thi article P1-crogrid/P1 element (ee Figure 3) are ued, which according to [8] atify it. Coupled problem. Lat part of our algoritm i the evaluation of the aerodynamical force. For that the velocity and preure value at adjacent triangle are ued. It mean, that we have according to (13) ( qi (X, t) = σ f ij (x) n j(x) = p + 2ν f 1 ( vi + v )) j n 2 x j x j, (30) i where n = (n 1, n 2) i the outer normal to boundary of tructure Γ Wt. Then it i evaluated in the middle of the common ide and with weight of one half added to the value in the vertice on thi ide. The value of aerodynamical force computed in vertice on the interface are then ued a dicrete verion of Neumann boundary condition of elatic olver. The algoritm olving coupled problem i implemented in emi-implicit form: We tart with the initial value v 0, p 0, u 0, q 0, A t0 and Ω f t 0. Then for n = 0, 1,... we proceed in computation in the following tep: 1. Baed on the preented cheme (20) it i acquired new olution u n+1 on n+1-th time layer, where q n+1 i extrapolated from q n. 2. From the known deformation u n+1 the ALE mapping A tn+1 i contructed by (21) and Ω f t n+1 i determined. Afterward we et w n+1 D (x) At (X) Atn (X) n+1 t, where x = A tn+1 (X). 3. We olve (29) and get v n+1, p n+1 defined on Ω f t n The value of the aerodynamical force on the interface are determined by (30) at time t n+1 from the known value v n+1, p n We et n := n + 1 and continue with the firt tep.
7 TOPICAL PROBLEMS OF FLUID MECHANICS Numerical reult The numerical reult for fluid flow interacting with the vocal fold model M5 uggeted by paper [9] are preented. The model M5 together with the triangulation i hown in Figure 4. Here, only one half of the channel wa ued a the computational domain with ymmetric boundary condition precribed at the top of the fluid domain (y = m). Figure 4: The triangulation of the computational domain Ωf0 and of the vocal fold model M5 (dimenion in [m]). Modal analyi. The eigenfrequencie of the vocal fold model wa determined by modal analyi. The olution of the ytem of ordinary differential equation Mu + Ku = 0, where the matrice M and K are given by (16), it i ought in the form u = ei ωj t uj, which lead to a generalized eigenvalue problem (K ωj2 M)uj = 0 = det(k ωj2 M) = 0. (31) The reult are hown in Table 1, where the Young modulu of elaticity wa E = 12 kpa, σ = 0, 4 inide Ωref and E = 100 kpa, σ = 0, 4 in a thin layer along the interface ΓWref. Frequency f1 f2 f3 f4 f5 [Hz] Frequency f1 f2 f3 x-direction y-direction Table 1: Left: The lowet eigenfrequencie obtained from the modal analyi of model M5. Right: The ignificant eigenfrequencie obtained from the pectral analyi of the time ignal at the point A. Dynamic of elatic tructure, energy conervation tet. Futhermore, the olution of the ytem (15) with zero body force and non-zero initial condition wa approximated by FEM, and the time ignal of a chooen point A [x = m, y = m] wa analyzed by Fourier tran 0 formation, ee Table 2. The initial condition were u1 = 0, u0 (x, y) = (y ) The frequencie obtained by the Fourier tranformation are hown in Table 2, where the motion of the point A wa analyzed. The frequencie agree well with the reult ummarized in Table 1. Another poibility, how to verify olver of elatic body, i to control energy conervation of vibrating tructure. The total energy E conit of um of kinetic and potential energy E = Ekin + Epot. Thee energie Ekin, Epot are given and then approximated in dicrete form a Z Z 1 T Ekin = ρ u dx u Mu, Epot = Cijkl eij (u)ekl (u) dx ut Ku. (32) Ω Ω
8 252 Prague, February 11-13, 2015 Figure 5: Graph of energy conervation, potential energy Epot i plotted with red colour and croe, kinetic energy Ekin i pictured by blue colour and behaviour of total energy E i hown with black colour. The reult hown in Figure 5 confirm energy conervation of implemented procedure (no dumping wa conidered). Flow in the channel with the fixed model of vocal fold. The flow olver wa teted for flow in the domain Ωfref without any interaction. Figure 6 how the approximation of the development of x-component velocity between time 0.05 and 0.3. The kinematic vicoity wa et ν f = 1, m/2 and the inlet boundary condition wa a parabolic profile: (Yi )(0.017 Yi ) vdir = 0.025, [Xi, Yi ] ΓfIn. (33) 0 Figure 6: The flow x-velocity pattern around the fixed model at ix different time intant t = to + j 0.05, (j = 0,..., 5). Figure are ordered from the left to the right.
9 TOPICAL PROBLEMS OF FLUID MECHANICS 253 The flow acceleration can be een in Figure 6 in the narrowet part of the channel. Reynold number for thi etting wa approximately 10. FSI tet. For teting interaction between fluid and elatic body, we adopted previou tet and to the elatic body we precribed all initial and boundary condition equal zero. The interaction wa enabled after time t = 0.1, when the velocity profile behind the vocal fold wa developed. The flow excited a periodic vibration of vocal fold with mall amplitude. Figure 7: Time evolution of x-diplacement of point A for the FSI tet (left) and it Fourier tranformation with dominant frequency f = 51 Hz (right). Figure 8: Time evolution of y-diplacement of point A for the FSI tet (left) and it Fourier tranformation with dominant frequency f = 51 Hz (right). Figure 9: Comparion of the x-diplacement of the point A computed with different time tep, t = (blue), t = (red with croe), for the FSI tet.
10 254 Prague, February 11-13, 2015 Figure 7 and 8 how time ignal of x-diplacement and y-diplacement of the point A together with their Fourier tranformation, repectively. The comparion of the time ignal computed by uing two different time tep in numerical procedure are in hown in Figure 9 demotrating a very good agreetment between thee olution. 5 Concluion The paper decribed the formulation of the FSI problem demontrated by the example of human fold vibration in airflow. The formulation ue FEM and the ALE method to imulate behaviour of the coupled ytem. The developed numerical cheme were implemented in an own program, which i able to handle complex geometrie. In the end the each part of the olver wa teted by imply tet, which verifie it baic functionality in the interaction problem. The acquired reult proved convergence of the computational cheme for low Reynold number and mall diplacement. Acknowledgment The financial upport for the preent project wa partly provided by the Czech Science Foundation under the Grant No. 101/11/0207 and project S GS 13/174/OHK2/3T/12. Reference [1] Brdička, M., Samek, L., Sopko, B.:Mechanika kontinua. Academia, Praha, [2] Curnier, A., Computational Method in Solid Mechanic, Kluwer Academic Publihing Group, Dodrecht, [3] Davi, T.: Umfpack. [mathematical library] Verion Univerity of Florida, USA. Available at: [4] Girault, V., Raviart, P., A.:Finite Element Method for Navier-Stoke Equation. Springer, Berlin, [5] Hadrava, M. & Feitauer M. & Horáček & Koík, A.: Dicontinuou Galerkin Method for the Problem of Linear Elaticity with Application to the Fluid-Structure Interaction. AIP Conference Proceeding, vol. 1558, pp ; [6] Koík, A. & Feitauer, M. & Horáček, J. & Sváček, P.: Numerical Simulation of Interaction of Human Vocal Fold and Fluid Flow Springer Proceeding in Phyic, vol. 139, pp , [7] Nomura, T. & Hughe, T. J. R.: An arbitrary Lagrangian-Eulerian finite element method for interaction of fluid and a rigid body. Computer Method in Applied Mechanic and Engineering, vol. 95, pp , [8] Quarteroni, A. & Valli, A.: Numerical Approximation of Partial Differential Equation. Springer, Berlin, [9] Scherer et al.: Intraglottal preure profile for a ymmetric and oblique glotti with a divergence angle of 10 degree Journal Acoutic Society of America, 4, 109, April [10] Sváček, P. & Horáček, J.: Numerical Simulation of Glottal Flow in Interaction with Self Ocillating Vocal Fold: Comparion of Finite Element Approximation with a Simplified Model, Communication in computational phyic, vol. 12, Iue: 3, pp , [11] Titze, I. R.: The Myoelatic Aerodynamic Theory of Phonation. National Center for Voice and Speech, USA, 2006.
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION REP- RESENTED BY HUMAN VOCAL FOLD
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION REP- RESENTED BY HUMAN VOCAL FOLD J. Valášek 1, P. Sváček 1 1 Department of Technical Mathematics, Faculty of Mechanical Engineering, Czech Technical
More informationThe Multilayer Impedance Pump Model
12 Chapter 2 The Multilayer Impedance Pump Model 2.1 Phyical model The MIP wa a luid-illed elatic tube with an excitation zone located aymmetrically with repect to the length o the pump. The pump had an
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 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 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 informationTWO DIMENSIONAL SIMULATION OF FLUID-STRUCTURE INTERACTION USING DGFEM
Proceedings of ALGORITMY 212 pp. 85 94 TWO DIMENSIONAL SIMULATION OF FLUID-STRUCTURE INTERACTION USING DGFEM JAROSLAVA HASNEDLOVÁ-PROKOPOVÁ,, MILOSLAV FEISTAUER, JAROMÍR HORÁČEK, ADAM KOSÍK, AND VÁCLAV
More informationFinite Element Truss Problem
6. rue Uing FEA Finite Element ru Problem We tarted thi erie of lecture looking at tru problem. We limited the dicuion to tatically determinate tructure and olved for the force in element and reaction
More informationOne Class of Splitting Iterative Schemes
One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi
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 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 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 informationDomain Optimization Analysis in Linear Elastic Problems * (Approach Using Traction Method)
Domain Optimization Analyi in Linear Elatic Problem * (Approach Uing Traction Method) Hideyuki AZEGAMI * and Zhi Chang WU *2 We preent a numerical analyi and reult uing the traction method for optimizing
More informationFinite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization
Finite Element Analyi of a Fiber Bragg Grating Accelerometer for Performance Optimization N. Baumallick*, P. Biwa, K. Dagupta and S. Bandyopadhyay Fiber Optic Laboratory, Central Gla and Ceramic Reearch
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 informationApproximate Analytical Solution for Quadratic Riccati Differential Equation
Iranian J. of Numerical Analyi and Optimization Vol 3, No. 2, 2013), pp 21-31 Approximate Analytical Solution for Quadratic Riccati Differential Equation H. Aminikhah Abtract In thi paper, we introduce
More informationChapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog
Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
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 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 informationA FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: CORRESPONDENCE: ABSTRACT
A FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: Zenon Medina-Cetina International Centre for Geohazard / Norwegian Geotechnical Intitute Roger
More informationLateral vibration of footbridges under crowd-loading: Continuous crowd modeling approach
ateral vibration of footbridge under crowd-loading: Continuou crowd modeling approach Joanna Bodgi, a, Silvano Erlicher,b and Pierre Argoul,c Intitut NAVIER, ENPC, 6 et 8 av. B. Pacal, Cité Decarte, Champ
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 informationBasics of a Quartz Crystal Microbalance
Baic of a Quartz Crytal Microbalance Introduction Thi document provide an introduction to the quartz crytal microbalance (QCM) which i an intrument that allow a uer to monitor mall ma change on an electrode.
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 informationLinear Motion, Speed & Velocity
Add Important Linear Motion, Speed & Velocity Page: 136 Linear Motion, Speed & Velocity NGSS Standard: N/A MA Curriculum Framework (006): 1.1, 1. AP Phyic 1 Learning Objective: 3.A.1.1, 3.A.1.3 Knowledge/Undertanding
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 informationNUMERICAL SIMULATION OF DESICCATION CRACKING PROCESS BY WEAK COUPLING OF DESICCATION AND FRACTURE
Geotec., Cont. Mat. & Env., ISSN:86-990, Japan, DOI: http://dx.doi.org/0.660/07.33.535 NUMERICAL SIMULATION OF DESICCATION CRACKING PROCESS BY WEAK COUPLING OF DESICCATION AND FRACTURE *Sayako Hirobe and
More informationOn the regularity to the solutions of the Navier Stokes equations via one velocity component
On the regularity to the olution of the Navier Stoke equation via one velocity component Milan Pokorný and Yong Zhou. Mathematical Intitute of Charle Univerity, Sokolovká 83, 86 75 Praha 8, Czech Republic
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 informationReformulation of Block Implicit Linear Multistep Method into Runge Kutta Type Method for Initial Value Problem
International Journal of Science and Technology Volume 4 No. 4, April, 05 Reformulation of Block Implicit Linear Multitep Method into Runge Kutta Type Method for Initial Value Problem Muhammad R., Y. A
More informationCRACK TIP STRESS FIELDS FOR ANISOTROPIC MATERIALS WITH CUBIC SYMMETRY
CRACK TIP TRE FIELD FOR ANIOTROPIC MATERIAL WITH CUBIC YMMETRY D.E. Lempidaki, N.P. O Dowd, E.P. Buo Department of Mechanical Engineering, Imperial College London, outh Kenington Campu, London, W7 AZ United
More informationResonance investigation of pump-turbine during startup process
IOP Conference Serie: Earth and Environmental Science OPEN ACCESS Reonance invetigation of pump-turbine during tartup proce To cite thi article: L Y He et al 2014 IOP Conf. Ser.: Earth Environ. Sci. 22
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 informationChapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem
Chapter 5 Conitency, Zero Stability, and the Dahlquit Equivalence Theorem In Chapter 2 we dicued convergence of numerical method and gave an experimental method for finding the rate of convergence (aka,
More informationModeling the scalar wave equation with Nyström methods
GEOPHYSICS, VOL. 71, NO. 5 SEPTEMBER-OCTOBER 200 ; P. T151 T158, 7 FIGS. 10.1190/1.2335505 Modeling the calar wave equation with Nytröm method Jing-Bo Chen 1 ABSTRACT High-accuracy numerical cheme for
More informationInteractive Processes Underlying the Production of Voice
APSIPA ASC 11 Xi an Interactive Procee Underlying the Production of Voice Tokihiko Kaburagi Kyuhu Univerity, Fukuoka, Japan E-mail: kabu@deign.kyuhu-u.ac.jp Tel/Fax: +81-9-553-457/45 Abtract Thi paper
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 informationMODELLING OF FRICTIONAL SOIL DAMPING IN FINITE ELEMENT ANALYSIS
MODELLING OF FRICTIONAL SOIL DAMPING IN FINITE ELEMENT ANALYSIS S. VAN BAARS Department of Science, Technology and Communication, Univerity of Luxembourg, Luxembourg ABSTRACT: In oil dynamic, the oil i
More informationResearch Article A New Kind of Weak Solution of Non-Newtonian Fluid Equation
Hindawi Function Space Volume 2017, Article ID 7916730, 8 page http://doi.org/10.1155/2017/7916730 Reearch Article A New Kind of Weak Solution of Non-Newtonian Fluid Equation Huahui Zhan 1 and Bifen Xu
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 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 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 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 4. The Laplace Transform Method
Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination
More informationCorrection for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002
Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in
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 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 informationinto a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get
Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}
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 informationDigital Control System
Digital Control Sytem - A D D A Micro ADC DAC Proceor Correction Element Proce Clock Meaurement A: Analog D: Digital Continuou Controller and Digital Control Rt - c Plant yt Continuou Controller Digital
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 informationMath 201 Lecture 17: Discontinuous and Periodic Functions
Math 2 Lecture 7: Dicontinuou and Periodic Function Feb. 5, 22 Many example here are taken from the textbook. he firt number in () refer to the problem number in the UA Cutom edition, the econd number
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 informationOverflow from last lecture: Ewald construction and Brillouin zones Structure factor
Lecture 5: Overflow from lat lecture: Ewald contruction and Brillouin zone Structure factor Review Conider direct lattice defined by vector R = u 1 a 1 + u 2 a 2 + u 3 a 3 where u 1, u 2, u 3 are integer
More informationThe Influence of the Load Condition upon the Radial Distribution of Electromagnetic Vibration and Noise in a Three-Phase Squirrel-Cage Induction Motor
The Influence of the Load Condition upon the Radial Ditribution of Electromagnetic Vibration and Noie in a Three-Phae Squirrel-Cage Induction Motor Yuta Sato 1, Iao Hirotuka 1, Kazuo Tuboi 1, Maanori Nakamura
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 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 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 informationEvolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis
Proceeding of 01 4th International Conference on Machine Learning and Computing IPCSIT vol. 5 (01) (01) IACSIT Pre, Singapore Evolutionary Algorithm Baed Fixed Order Robut Controller Deign and Robutne
More informationStresses near a plate vertex due to a shear force on one of the edges
Stree near a plate vertex due to a hear force on one of the edge P.C.J. Hoogenboom Delft Univerity of Technology, Faculty of Civil Engineering and Geocience, Delft, the Netherland A cloed form olution
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 informationDetermination of the local contrast of interference fringe patterns using continuous wavelet transform
Determination of the local contrat of interference fringe pattern uing continuou wavelet tranform Jong Kwang Hyok, Kim Chol Su Intitute of Optic, Department of Phyic, Kim Il Sung Univerity, Pyongyang,
More informationThe combined Laplace-homotopy analysis method for partial differential equations
Available online at wwwir-publicationcom/jmc J Math Computer Sci 6 (26), 88 2 Reearch Article The combined Laplace-homotopy analyi method for partial differential equation Javad Vahidi Department of Mathematic,
More informationDimension Effect on Dynamic Stress Equilibrium in SHPB Tests
International Journal of Material Phyic. ISSN 97-39X Volume 5, Numer 1 (1), pp. 15- International Reearch Pulication Houe http://www.irphoue.com Dimenion Effect on Dynamic Stre Equilirium in SHPB Tet Department
More informationEffects of soil structure interaction on behavior of reinforced concrete structures
Journal of Structural Engineering & Applied Mechanic 18 Volume 1 Iue 1 Page 8-33 http://doi.org/1.3146/jeam.18.1833 www.goldenlightpublih.com RESEARCH ARTICLE Effect of oil tructure interaction on behavior
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 informationSimple Observer Based Synchronization of Lorenz System with Parametric Uncertainty
IOSR Journal of Electrical and Electronic Engineering (IOSR-JEEE) ISSN: 78-676Volume, Iue 6 (Nov. - Dec. 0), PP 4-0 Simple Oberver Baed Synchronization of Lorenz Sytem with Parametric Uncertainty Manih
More informationCompact finite-difference approximations for anisotropic image smoothing and painting
CWP-593 Compact finite-difference approximation for aniotropic image moothing and painting Dave Hale Center for Wave Phenomena, Colorado School of Mine, Golden CO 80401, USA ABSTRACT Finite-difference
More informationModule 4: Time Response of discrete time systems Lecture Note 1
Digital Control Module 4 Lecture Module 4: ime Repone of dicrete time ytem Lecture Note ime Repone of dicrete time ytem Abolute tability i a baic requirement of all control ytem. Apart from that, good
More informationOnline Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat
Online Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat Thi Online Appendix contain the proof of our reult for the undicounted limit dicued in Section 2 of the paper,
More informationNORMAL MODE ANALYSIS OF SECOND-ORDER PROJECTION METHODS FOR INCOMPRESSIBLE FLOWS. Jae-Hong Pyo and Jie Shen. (Communicated by Z.
DISCRETE AND CONTINUOUS Webite: http://aimcience.org DYNAMICAL SYSTEMS SERIES B Volume 5 Number 3 Augut 005 pp. 817 840 NORMAL MODE ANALYSIS OF SECOND-ORDER PROJECTION METHODS FOR INCOMPRESSIBLE FLOWS
More informationMucus Transport in the Human Lung Airways: Effect of Porosity Parameter and Air Velocity
Mucu Tranport in the Human Lung Airway: Effect of Poroity Parameter and Air Velocity V.S. Verma 1, Vikah Rana Department of Mathematic and Statitic, DDU Gorakhpur Univerity, Gorakhpur-73009, U.P. (India)
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 informationHyperbolic Partial Differential Equations
Hyperbolic Partial Differential Equation Evolution equation aociated with irreverible phyical procee like diffuion heat conduction lead to parabolic partial differential equation. When the equation i a
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 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 informationAccuracy of Symmetric Partitioned Runge-Kutta Methods for Differential Equations on Lie-Groups
Accuracy of Symmetric Partitioned Runge-Kutta Method for Differential Equation on Lie-Group WUB / 11-23 BUW-IMACM 11-19, Michael Günther, Franceco Knechtli and Michèle Wandelt Bergiche Univerität Wuppertal,
More informationUNITS FOR THERMOMECHANICS
UNITS FOR THERMOMECHANICS 1. Conitent Unit. Every calculation require a conitent et of unit. Hitorically, one et of unit wa ued for mechanic and an apparently unrelated et of unit wa ued for heat. For
More informationFinding the location of switched capacitor banks in distribution systems based on wavelet transform
UPEC00 3t Aug - 3rd Sept 00 Finding the location of witched capacitor bank in ditribution ytem baed on wavelet tranform Bahram nohad Shahid Chamran Univerity in Ahvaz bahramnohad@yahoo.com Mehrdad keramatzadeh
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 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 informationIEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation
IEOR 316: Fall 213, Profeor Whitt Topic for Dicuion: Tueday, November 19 Alternating Renewal Procee and The Renewal Equation 1 Alternating Renewal Procee An alternating renewal proce alternate between
More informationA Coupled Momentum Method to Model Blood Flow in Deformable Arteries
COMPUTATIONAL MECHANICS WCCM VI in conjunction with APCOM 4, Sept. 5-, 24, Beijing, China 24 Tinghua Univerity Pre & Springer-Verlag A Coupled Momentum Method to Model Blood Flow in Deformle Arterie Alberto
More informationEXTENDED STABILITY MARGINS ON CONTROLLER DESIGN FOR NONLINEAR INPUT DELAY SYSTEMS. Otto J. Roesch, Hubert Roth, Asif Iqbal
EXTENDED STABILITY MARGINS ON CONTROLLER DESIGN FOR NONLINEAR INPUT DELAY SYSTEMS Otto J. Roech, Hubert Roth, Aif Iqbal Intitute of Automatic Control Engineering Univerity Siegen, Germany {otto.roech,
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 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 informationINVESTIGATION OF THE TURBULENT MASS TRANSPORT DURING THE MIXING OF A STABLE STRATIFICATION WITH A FREE JET USING CFD-METHODS
INVESTIGATION OF THE TURBULENT MASS TRANSPORT DURING THE MIXING OF A STABLE STRATIFICATION WITH A FREE JET USING CFD-METHODS I1. Armin Zirkel, I2. Eckart Laurien Intitute of Nuclear Technology and Energy
More informationMulti-dimensional Fuzzy Euler Approximation
Mathematica Aeterna, Vol 7, 2017, no 2, 163-176 Multi-dimenional Fuzzy Euler Approximation Yangyang Hao College of Mathematic and Information Science Hebei Univerity, Baoding 071002, China hdhyywa@163com
More informationA Buckling Problem for Graphene Sheets. J. Galagher 1, Y. Milman 2, S. Ryan 3, D. Golovaty 3, P. Wilber 3, and A. Buldum 4
A Buckling Problem for Graphene Sheet J. Galagher 1, Y. Milman 2, S. Ryan 3, D. Golovaty 3, P. Wilber 3, and A. Buldum 4 1 Department of Phyic, Rocheter Intitute of Technology, Rocheter, NY 14623, USA
More informationConvective Heat Transfer
Convective Heat Tranfer Example 1. Melt Spinning of Polymer fiber 2. Heat tranfer in a Condener 3. Temperature control of a Re-entry vehicle Fiber pinning The fiber pinning proce preent a unique engineering
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 informationME 375 EXAM #1 Tuesday February 21, 2006
ME 375 EXAM #1 Tueday February 1, 006 Diviion Adam 11:30 / Savran :30 (circle one) Name Intruction (1) Thi i a cloed book examination, but you are allowed one 8.5x11 crib heet. () You have one hour to
More informationUnified Design Method for Flexure and Debonding in FRP Retrofitted RC Beams
Unified Deign Method for Flexure and Debonding in FRP Retrofitted RC Beam G.X. Guan, Ph.D. 1 ; and C.J. Burgoyne 2 Abtract Flexural retrofitting of reinforced concrete (RC) beam uing fibre reinforced polymer
More informationAdaptive Control of Level in Water Tank: Simulation Study
Adaptive Control of Level in Water Tank: Simulation Study Jiri Vojteek and Petr Dotal Abtract An adaptive control i a popular, o called modern control method which could be ued for variou type of ytem
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 informationSymmetry Lecture 9. 1 Gellmann-Nishijima relation
Symmetry Lecture 9 1 Gellmann-Nihijima relation In the lat lecture we found that the Gell-mann and Nihijima relation related Baryon number, charge, and the third component of iopin. Q = [(1/2)B + T 3 ]
More information1 Parity. 2 Time reversal. Even. Odd. Symmetry Lecture 9
Even Odd Symmetry Lecture 9 1 Parity The normal mode of a tring have either even or odd ymmetry. Thi alo occur for tationary tate in Quantum Mechanic. The tranformation i called partiy. We previouly found
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 informationA Simplified Dynamics Block Diagram for a Four-Axis Stabilized Platform
A Simplified Dynamic Block Diagram for a FourAxi Stabilized Platform Author: Hendrik Daniël Mouton a Univerity of Cape Town, Rondeboch, Cape Town, South Africa, 770 Abtract: t i relatively traightforward
More informationMidterm Test Nov 10, 2010 Student Number:
Mathematic 265 Section: 03 Verion A Full Name: Midterm Tet Nov 0, 200 Student Number: Intruction: There are 6 page in thi tet (including thi cover page).. Caution: There may (or may not) be more than one
More information