Flow Induced Vibration
|
|
- Madlyn Long
- 6 years ago
- Views:
Transcription
1 Flow Induced Vbraton Project Progress Report Date: 16 th November, 2005 Submtted by Subhrajt Bhattacharya Roll no.: 02ME101 Done under the gudance of Prof. Anrvan Dasgupta Department of Mechancal Engneerng, T Kharagpur Department of Mechancal Engneerng, Indan Insttute of Technology, Kharagpur
2 Certfcate Ths s to certfy that the report enttled Flow Induced Vbraton submtted by Subhrajt Bhattacharya to the Department of Mechancal Engneerng, T Kharagpur, s a bona fde record of work carred out under my supervson and gudance. Prof. Anrvan Dasgupta, Dept. of Mechancal Engneerng, T Kharagpur. Date:
3 1. Introducton: The problem of flud-structure nteracton s encountered n varous Engneerng applcatons. The present problem deals wth vbratons nduced nto structures due to flow takng place on ts surface, and hence analyzng the stablty of the flow. Though the present approach to the problem s grossly smplfed and no substantal smulaton or numercal soluton have yet been deduced, the works done tll now have great potental n applcatons lke flow over underwater vehcle, etc. 2. Problem defnton and am: Our present am s to analyze the stablty of flud flows over a flexble flat surface. By the term stablty we mean that we try to determne the crtcal Reynolds number for a gven flow profle over the structure. We wll consder flexblty of the plate and frame the correspondng equatons. However for the purpose of determnng the numercal soluton we have presently nvestgated only the case of flow over a rgd surface to check f we obtan the standard results for stablty of parallel flow over rgd surfaces avalable n standard text. 3. Orgn of Turbulence and ways to analyze stablty of a flow: Flow nduced vbraton s caused n structures by forcng due to tme varant pressure actng on the surface of the structure. If the flow s magned to be a lnear superposton of a steady state lamnar flow and a perturbaton flow feld, the lamnar flow won t cause the forcng on the structure as t acts as a tme nvarant pressure on the structure s surface. It s the perturbaton flow feld that vares wth tme and hence produces forcng on the structure s surface. Ths same perturbaton flow feld s the sole cause of Turbulence n the flow. The orgn of the perturbaton flow feld may be somethng lke a very small dsturbance caused n the lamnar flow feld. A flow s sad to be stable f for any small ntal dsturbance (.e. perturbaton) added to the lamnar flow feld, the perturbaton flow feld gradually des down wth tme. It wll be termed as a Turbulent flow f the perturbaton flow feld gets magnfed wth tme. Hence our prmary approach wll be to nvestgate the nature of varaton of a perturbaton flow feld wth tme.
4 . The Orr-Somerfeld equaton: The two-dmensonal ncompressble flow over the plate s governed by the Naver Stoke s equatons, 2 2 u u u 1 1 p u v F u u + + = X + ν t x y ρ ρ x x y (1a) 2 2 v v v 1 1 p v v + u + v = FY + ν t x y ρ ρ y x y (1b) and the contnuty equaton, u v + = 0 (2) x y Let us consder a steady-state lamnar flow over an nfntely long plate. Y fg - 1 X For ths flow, u = U(y) v = 0 (3) p = P(x) Now let the perturbaton feld be descrbed by the perturbaton components denoted by a prme upon the steady-state lamnar varables. Hence the fnal flow feld wll be descrbed by, u = U(y) + u (x, y, t) v = v (x, y, t) () p = P(x) + p (x, y, t) Now, both the steady-state lamnar feld (3) and the superposed feld () satsfes equaton (1) and (2). Hence by substtutng them n (1) and (2) and performng some smplfcaton, we obtan the equatons governng the perturbaton flow feld, u' u' U 1 p' 2 + U + v' + = ν u ' (5a) t x y ρ x v' v' 1 p' 2 + U + = ν v' (5b) t x ρ y u' v' + = 0 (5c) x y From these equatons, wth a known velocty profle U(y), we may obtan solutons for u', v', and p'. Our present am wll be to determne for a gven velocty profle U(y), the coeffcent of vscosty ν and the boundary & ntal condtons of the perturbaton felds, whether or not the perturbaton components de down wth tme.
5 In order to satsfy eqn. (5c), we defne a stream functon ψ ( xyt,, ) such that ψ ψ u ' = and v ' = (6) y x It s now assumed that the dsturbances, and hence ψ s superposton of several perodc dsturbances (perodc n x) propagatng along the drecton of flow. Hence, we substtute, ( α x βt) ψ ( xyt,, ) = φ( ye ) (7) 2π ψ s sad to be perodc on x wth frequency α and wavelength λx =. α Though α can be assumed to be real, β beng the tme frequency should be assumed to be complex n order to keep the possblty of non-perodc magnfcaton or decay of ψ wth tme. The fnal ψ wll be superposton of all the solutons of ψ. We defne the complex velocty of propagaton of the dsturbance as β c= = cr + c (8) α [ ct ( x crt ψ φ( y) e α + )] = It may be noted that for c < 0, the soluton of ψ des down wth tme, and hence so does the perturbaton components. Thus the flow s stable for c < 0 and tends to become turbulent for c > 0. We put u' and v' n terms of φ, α and β n (5a) and (5b) and elmnate p' to obtan a sngle equaton. We non-dmensonalze the equatons by redefnng y as y δ, U as U U, c m c as U, where U m s the free-stream velocty of the flow and δ s the boundary layer m thckness. The non-dmensonalzed equaton hence obtaned s called the Orr- Sommerfeld equaton: where, R U δ ν 2 ( U c)( φ αφ) U φ= ( φ 2αφ 2 + αφ ) (9) α R m = s the Reynolds Number. A trval soluton to eqn.(9) s φ = 0. For non-trval soluton, for a gven α and R, we obtan an egenvalue of c and the correspondng egenfuncton φ. Hence our mmedate target s to fnd an egenvalue soluton of (9).
6 5. Boundary condtons for the Orr-Sommerfeld equaton: A. Boundary condtons for flow over a rgd, statc plate: For ths case, At y = 0, v' = 0 and u' = 0 φ (0) = 0 and φ (0) = 0 At y, v' = 0 and u' = 0 φ( ) = 0 and φ ( ) = 0 B. Boundary condtons for flow over a flexble plate modeled as a beam: p'(x,0,t) We model the plate as nfnte flexble beam. Hence, the dfferental equaton governng the moton of the beam s gven by, 2 fg - 2 λ w w = p ( x,0, t) EI (10) 2 t x Where, w denotes the dsplacement of the beam n Y drecton, λ s mass per unt length of the beam. The forcng on the beam s due to the tme-varant perturbaton pressure. It s to be noted n ths case that equatons (5), (9) and (10) gets coupled. One way of solvng the equaton wll be as follows. ( α x βt) As p' s of the form of e, we may wrte for eqn. (10), ( α x βt) wxt (, ) = ke (11) Substtutng ths w n (10) and hence substtutng the hence obtaned p'(x,0,t) nto (5a), and Substtutng v' n terms of φ n the same (5a) equaton, all at y = 0, we obtan, ρ αu φ α ν β φ νφ k + = α λβ Hence the boundary condtons become, At y = 0, 2 (0) (0) ( ) (0) (0) 2 ( EIα ) v' = w and u' = 0 β φ(0) = k and φ (0) = 0, where k s the non-dmensonalzed k. α At y, v' = 0 and u' = 0 φ( ) = 0 and φ ( ) = 0 Hence, as t can be seen, one of the boundary condtons s a bt more complex nvolvng φ, φ, φ and c. (12)
7 The numercal technque for solvng the Orr-Sommerfeld equaton for the egenvalues and egenfunctons n ether of the cases wll reman smlar, except that the boundary condtons are modfed. However tll date we have only nvestgated the case A,.e. the case of flow over a rgd, statc plate. 6. Numercal methods attempted for solvng the Orr-Sommerfeld equaton: A. Galerkn s Method: Ths method, though can handle the case of flow over a rgd, statc plate satsfactorly, ts applcaton n solvng the case of flexble plate s dffcult. The method for the case of flow over a rgd, statc plate s descrbed below n bref. We denote eqn.(9) by Γ ( φ) = 0, where Γ denotes the operator. We wrte φ as a lnear superposton of several functons φ that satsfy the boundary y condtons gven n 5.A. Such functons were chosen to be of the form, ye η φ = µ. Hence we wrte, We now defne an error, n φ( y) = ξφ( y) = 1 ey ( ) ( ) n =Γ ξφ y = 1 We need to choose the values of ξ n such a way that ths error s mnmzed. We do that by solvng the set of n equatons, e, φ = 0, = 1 to n (1) where, χ, χ denotes the nner product gven by, j e, φ e = φ dy The n equatons n (1) have ξ as the unknowns and can be represented n the matrx form as, ξ1 [ M ] = 0 (15) ξ n Where the matrx [M] contans the unknown c. For non-trval soluton of ths equaton, we must have, [ M ] = 0 (16) 0 (13)
8 In general we ll obtan n solutons for c. We should chose the one for whch the ξ s are n such that Γ ξφ ( y) s mnmum. = 1 Ths method was mplemented n Mathemata 5.1 and a rather unsatsfactory result was obtaned probably due to the followng reasons:. Due to lmtaton of computatonal power, n had to be lmted to n order to get a satsfactory and accurate ntegraton value of the nner product and soluton of c from (16).. The choce of η and µ n choosng the functonal forms were done arbtrarly. We determned the egenvalues c for dfferent α and R and plotted the contour n the α-r plane for whch c = 0. Ths contour wll mark the margn between the stable and unstable zones of α and R. However the contour c = 0 could not be found satsfactorly n the fst quadrant of α-r plane. But on plottng a contour plot of c, the followng was obtaned. The plot demonstrates that the basc shape of the standard results s beng approached by the soluton: fg - 3 The horzontal axs denotes R and the vertcal axs denotes α.
9 B. Automated search of egenvalues Integraton by Runge-Kutta: Ths method prmarly proposed by Betchov & Crmnale [2], deals the regons above and below the boundary layer separately. We frst nvestgate the Orr-Sommerfeld equaton (9) for y > 1. In ths regon, the non-dmensonalsed velocty s U = 1. Hence the equaton s modfed as, 2 2 (1 c)( φ αφ) = ( φ 2αφ + αφ) for y > 1 (17) α R Ths equaton beng th order lnear n φ wth constant coeffcent has smple analytcal soluton gven by, where, p j are the roots of, Hence, p1 j φ( y) = Aje j= 1 p y (18) (1 c) ( p α ) = ( p α ) (19) α R = α, p1 R p3 = α 1 + (1 c) α 1 2, = α R p3 = α 1 + (1 c) α From boundary condton at, as y, φ = 0 and φ = 0. Hence we have, A 1 = A 3 = 0 for y > 1. p2y py φ( y) = A e + A e. 2 Now we argue that, snce for y > 1, the soluton s a lnear superposton of two modes p2 y p y e and e, the soluton for y < 1 wll also be lnear superposton of these two same modes. Hence now our task s to fnd the soluton of these two modes n the regon y < 1. We attan ths by performng two ntegraton passes from y = 1 to y = 0 ndependently. We used Runge-Kutta method for ths numercal ntegraton. Integraton Pass I: We start from y = 1 wth A 2 = 0 and A = 1(or some other value). p p 2 p Therefore we take the ntal values φ (1) = e, φ (1) = pe, φ (1) = pe & φ (1) = pe 3 p and move on ntegratng towards y = 0. Let the soluton obtaned n ths process be called φ ( y). Integraton Pass : Smlarly wth A 2 = 1(or some other value) and A = 0 we obtan the second pass ntegraton φ ( y). Hence the fnal soluton s of the form φ( y) = a φ ( y) + a φ ( y). I I I 1 2
10 From the boundary condtons at y = 0, φ(0) = aiφi (0) + aφ (0) = 0 and φ (0) = aiφ I (0) + aφ (0) = 0, for non-trval soluton of a I and a, we must have, φi(0) φ(0) = 0 φ (0) φ (0) I (20) It s to be noted that the only unknown n (20) for a gven α and R s c. Hence from 20 we obtan the egenvalue c. The correspondng egenvector gves a I and a, and hence the fnal soluton of φ ( y). Search for Egenvalue n c-plane: However t was not possble to solve c explctly from (20). Hence we assumed some c and determned the solutons from the two ntegraton passes, φ I ( y) and φ ( y). For those partcular solutons we defned, φi(0) φ(0) f() c = φ (0) φ (0) I As we know f(c) should converge to 0, we used an teraton scheme [2] as follows. We start wth an arbtrary value of c and some small value of δ c and go on updatng t usng the followng teraton, f( c+ δc) c= λ 1 δc f() c c c+ c δ c µ c where, λ = 1, 0.8 or 0.5 dependng on whether the convergence s fast, moderate or slow. 1 and, µ =. The teraton contnues tll c reaches a substantally small value. 1
11 Lke before, we once agan searched for the contour of c = 0 n the α-r plane. The results obtaned, though not satsfactory, s descrbed n the followng plot of 20 ponts: fg - The prmary reason for the devaton of ponts from a sngle smooth contour s the hgh oscllaton of the soluton φ ( y) as the second numercal ntegraton approaches y = 0. Ths fact s well demonstrated n the followng plot of the ampltude of φ ( y) vs. y for a partcular α and R: fg - 5
12 The small ampltude of φ ( y) s due to the ntal choce of a small A 2 for the second Integraton pass. C. Modfed second Integraton Pass usng Statons n between: Ths method, as explaned by Betchov & Crmnale [2], was mplemented n order to reduce the oscllaton of the second soluton. The man prncple of ths method s based on the choce of some statons n between y = 1 and y = 0. At these statons the second ntegraton s paused and t s updated by lnearly combnng wth the frst ntegraton φ I ( y) to make A = 0. We are presently workng on ths technque and hope to obtan some satsfactory result very soon. Once we obtan a soluton for the case of flow over a rgd, statc plate, we ll attempt the soluton for the case of flow over a flexble plate (modeled as an nfnte beam) on the smlar lne. References: [1] Dr. Hermann Schlchtng, Boundary Layer Theory, McGraw Hll, [2] Robert Betchov, Wllam O. Crmnale, Jr., Stablty of Parallel Flows, Academc Press, 1967.
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More information12. The Hamilton-Jacobi Equation Michael Fowler
1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationΔ x. u(x,t) Fig. Schematic view of elastic bar undergoing axial motions
ME67 - Handout 4 Vbratons of Contnuous Systems Axal vbratons of elastc bars The fgure shows a unform elastc bar of length and cross secton A. The bar materal propertes are ts densty ρ and elastc modulus
More informationΔ x. u(x,t) Fig. Schematic view of elastic bar undergoing axial motions
ME67 - Handout 4 Vbratons of Contnuous Systems Axal vbratons of elastc bars The fgure shows a unform elastc bar of length and cross secton A. The bar materal propertes are ts densty ρ and elastc modulus
More informationTHEOREMS OF QUANTUM MECHANICS
THEOREMS OF QUANTUM MECHANICS In order to develop methods to treat many-electron systems (atoms & molecules), many of the theorems of quantum mechancs are useful. Useful Notaton The matrx element A mn
More informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationElectrical double layer: revisit based on boundary conditions
Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationProf. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model
EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More informationA particle in a state of uniform motion remain in that state of motion unless acted upon by external force.
The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,
More informationPARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY
POZNAN UNIVE RSITY OF TE CHNOLOGY ACADE MIC JOURNALS No 86 Electrcal Engneerng 6 Volodymyr KONOVAL* Roman PRYTULA** PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY Ths paper provdes a
More informationSolution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method
Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another well-known teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222
More informationThree views of mechanics
Three vews of mechancs John Hubbard, n L. Gross s course February 1, 211 1 Introducton A mechancal system s manfold wth a Remannan metrc K : T M R called knetc energy and a functon V : M R called potental
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationConsistency & Convergence
/9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency
More informationPrinciples of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer
Prncples of Food and Boprocess Engneerng (FS 31) Solutons to Example Problems on Heat Transfer 1. We start wth Fourer s law of heat conducton: Q = k A ( T/ x) Rearrangng, we get: Q/A = k ( T/ x) Here,
More informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
More informationSnce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t
8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More informationChapter 8. Potential Energy and Conservation of Energy
Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal
More information6.3.4 Modified Euler s method of integration
6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationPHYS 705: Classical Mechanics. Canonical Transformation II
1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m
More informationLimited Dependent Variables
Lmted Dependent Varables. What f the left-hand sde varable s not a contnuous thng spread from mnus nfnty to plus nfnty? That s, gven a model = f (, β, ε, where a. s bounded below at zero, such as wages
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationIn this section is given an overview of the common elasticity models.
Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp
More informationGrid Generation around a Cylinder by Complex Potential Functions
Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around
More information14 The Postulates of Quantum mechanics
14 The Postulates of Quantum mechancs Postulate 1: The state of a system s descrbed completely n terms of a state vector Ψ(r, t), whch s quadratcally ntegrable. Postulate 2: To every physcally observable
More information10. Canonical Transformations Michael Fowler
10. Canoncal Transformatons Mchael Fowler Pont Transformatons It s clear that Lagrange s equatons are correct for any reasonable choce of parameters labelng the system confguraton. Let s call our frst
More informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationSTATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION
STATIC ANALYSIS OF TWO-LERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION Ákos József Lengyel István Ecsed Assstant Lecturer Emertus Professor Insttute of Appled Mechancs Unversty of Mskolc Mskolc-Egyetemváros
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationLagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013
Lagrange Multplers Monday, 5 September 013 Sometmes t s convenent to use redundant coordnates, and to effect the varaton of the acton consstent wth the constrants va the method of Lagrange undetermned
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More information5.04, Principles of Inorganic Chemistry II MIT Department of Chemistry Lecture 32: Vibrational Spectroscopy and the IR
5.0, Prncples of Inorganc Chemstry II MIT Department of Chemstry Lecture 3: Vbratonal Spectroscopy and the IR Vbratonal spectroscopy s confned to the 00-5000 cm - spectral regon. The absorpton of a photon
More informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationLECTURE 9 CANONICAL CORRELATION ANALYSIS
LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of
More informationMechanics Physics 151
Mechancs Physcs 5 Lecture 0 Canoncal Transformatons (Chapter 9) What We Dd Last Tme Hamlton s Prncple n the Hamltonan formalsm Dervaton was smple δi δ p H(, p, t) = 0 Adonal end-pont constrants δ t ( )
More information2.3 Nilpotent endomorphisms
s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms
More informationQuantum Mechanics for Scientists and Engineers. David Miller
Quantum Mechancs for Scentsts and Engneers Davd Mller Types of lnear operators Types of lnear operators Blnear expanson of operators Blnear expanson of lnear operators We know that we can expand functons
More informationρ some λ THE INVERSE POWER METHOD (or INVERSE ITERATION) , for , or (more usually) to
THE INVERSE POWER METHOD (or INVERSE ITERATION) -- applcaton of the Power method to A some fxed constant ρ (whch s called a shft), x λ ρ If the egenpars of A are { ( λ, x ) } ( ), or (more usually) to,
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationInexact Newton Methods for Inverse Eigenvalue Problems
Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More information5 The Rational Canonical Form
5 The Ratonal Canoncal Form Here p s a monc rreducble factor of the mnmum polynomal m T and s not necessarly of degree one Let F p denote the feld constructed earler n the course, consstng of all matrces
More information), it produces a response (output function g (x)
Lnear Systems Revew Notes adapted from notes by Mchael Braun Typcally n electrcal engneerng, one s concerned wth functons of tme, such as a voltage waveform System descrpton s therefore defned n the domans
More informationWaveguides and resonant cavities
Wavegudes and resonant cavtes February 8, 014 Essentally, a wavegude s a conductng tube of unform cross-secton and a cavty s a wavegude wth end caps. The dmensons of the gude or cavty are chosen to transmt,
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationPES 1120 Spring 2014, Spendier Lecture 6/Page 1
PES 110 Sprng 014, Spender Lecture 6/Page 1 Lecture today: Chapter 1) Electrc feld due to charge dstrbutons -> charged rod -> charged rng We ntroduced the electrc feld, E. I defned t as an nvsble aura
More informationwhere the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt
Physcs 543 Quantum Mechancs II Fall 998 Hartree-Fock and the Self-consstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes.
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationChapter 3 Differentiation and Integration
MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton
More informationFundamental loop-current method using virtual voltage sources technique for special cases
Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,
More informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationPHYS 705: Classical Mechanics. Hamilton-Jacobi Equation
1 PHYS 705: Classcal Mechancs Hamlton-Jacob Equaton Hamlton-Jacob Equaton There s also a very elegant relaton between the Hamltonan Formulaton of Mechancs and Quantum Mechancs. To do that, we need to derve
More informationThe Second Anti-Mathima on Game Theory
The Second Ant-Mathma on Game Theory Ath. Kehagas December 1 2006 1 Introducton In ths note we wll examne the noton of game equlbrum for three types of games 1. 2-player 2-acton zero-sum games 2. 2-player
More informationNovember 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project
SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω
More informationHow Differential Equations Arise. Newton s Second Law of Motion
page 1 CHAPTER 1 Frst-Order Dfferental Equatons Among all of the mathematcal dscplnes the theory of dfferental equatons s the most mportant. It furnshes the explanaton of all those elementary manfestatons
More informationPHYS 705: Classical Mechanics. Newtonian Mechanics
1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]
More informationOptimal Control of Temperature in Fluid Flow
Kawahara Lab. 5 March. 27 Optmal Control of Temperature n Flud Flow Dasuke YAMAZAKI Department of Cvl Engneerng, Chuo Unversty Kasuga -3-27, Bunkyou-ku, Tokyo 2-855, Japan E-mal : d33422@educ.kc.chuo-u.ac.jp
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationk t+1 + c t A t k t, t=0
Macro II (UC3M, MA/PhD Econ) Professor: Matthas Kredler Fnal Exam 6 May 208 You have 50 mnutes to complete the exam There are 80 ponts n total The exam has 4 pages If somethng n the queston s unclear,
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationCausal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim
Causal Damonds M. Aghl, L. Bombell, B. Plgrm Introducton The correcton to volume of a causal nterval due to curvature of spacetme has been done by Myrhem [] and recently by Gbbons & Solodukhn [] and later
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationEigenvalues of Random Graphs
Spectral Graph Theory Lecture 2 Egenvalues of Random Graphs Danel A. Spelman November 4, 202 2. Introducton In ths lecture, we consder a random graph on n vertces n whch each edge s chosen to be n the
More informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationTurbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH
Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant
More information