UMEÅ UNIVERSITY September 1, 2016 Computational Science and Engineering Modeling and Simulation. Dynamical systems. Peter Olsson
|
|
- Lionel Garrison
- 5 years ago
- Views:
Transcription
1 UMEÅ UNIVERSITY Septembe, 26 Computational Science and Engineeing Modeling and Simulation Dynamical systems Pete Olsson
2 Continuous population models fo a single species. Continuous gowth models The simplest possible: dn dt = biths deaths+migation. dn dt = bn dn N(t) = N e (b d)t, can t be tue foeve. Thee should be a self-limiting pocess when the population gets too lage. Intoduce a caying capacity, K: This is called logistic gowth in a population. dn ( dt = N N ). () K.. Steady states and lineaization Thee ae two steady states (equilibium states), N = and N = K. The fist one is unstable since lineaization aound N = gives dn dt = N and so N gows exponentially. The second steady state is stable. To see this, wite N = K + n, and lineaize: d(k +n) dt ( = (K +n) K +n ) K dn ( dt K n ) = n, K which gives an exponential decay, n(t) = n()e t. The caying capacity detemines the size of the steady state population wheeas / is the chaacteistic time scale of the esponse of the model to a change in the population.
3 ..2 Geneal fomulation In the geneal fomulation we have dn = f(n), (2) dt whee f(n) is a non-linea function of N. Then steady state solutions, N, ae solutions to f(n) =, they ae linealy stable if f (N ) < and linealy unstable if f (N ) >. dn/dt N To show this, conside the deviations fom N : n = N N and dn/dt = dn/dt: dn dt = f(n) f(n )+nf (N ) = nf (N ), with the solution n(t) = n()e f (N )t, which inceases if f (N ) >..2 Insect outbeak model: Spuce Budwom We ae now consideing the Spuce Budwom model which is a pactical model with two linealy stable steady state populations. The dynamics is given by dn dt = BN ( N ) p(n), KB whee p epesents pedation by bids. We now specialize to the following fom fo pedation: Pedation N p(n) = BN2 A 2 +N 2 which has a change fom low to high pedation at an appoximate theshold value N c = A. The dynamics now becomes dn dt = BN ( N KB ) BN2 A 2 +N 2. 2
4 .2. Nondimensionalisation This model has fou paametes, B, K B, B, and A and it is then essentially impossible to analyze how the behavio changes with the paamete values. To simplify things we intoduce the dimensionless quantities which leads to the equation u = N A, = A B B, q = K B A, τ = Bt A, ( du = u u ) u2 q +u = f(u;,q). 2 This now depends on only two paametes and q which ae pue numbes. Note that the time scale is also changed. Also note that thee ae othe possible ways to get a dimensionless population, the most obvious one being ũ = N/K B. Diffeent fomulations may give diffeent insights..2.2 Gaphical solution The steady states ae solutions of f(u;,q) = u ( u q ) = u2 +u 2, and since u = clealy is a solution this may be ewitten ( u q ) = u +u u +u 2 ( u q. 2 u ) The ight hand side of the equation is a cuve with a non-tivial shape which is independent of and q. The left hand side on the othe hand depends on the paametes but is simply the staight line (/q)u fom (,) to (q,). The solution to the equation is given by the intesections of these cuves as shown above. Fom diffeent and q one gets diffeent staight lines, and it is clea that the equation can have eithe one o thee solutions. Fo the case with thee solutions we denote them by u, u 2, and u 3. 3
5 To investigate which of these steady states solutions that ae stable, we conside du/ f(u;,q) and the sign of f/ u. Since f >, u u2 u = u 2 is linealy unstable. The othe two, u = u and u = u 3 ae stable as they have negative deivatives. f(u).5. =.4 q = 3 f(u) = u ( u q ) u 2 +u u.2.3 Hysteesis It is inteesting to examine the effect of a gadual change of the paamete as it tuns out that it leads to hysteesis. Stating with = the only solution is u = u =. When inceases the u - equilibium inceases monotonically until u themaximumofthecuve iseached +u 2 whee = h. When is inceased futhe this steady state disappeas and the equilibium value jumps to u 3. If is deceased again, the solution stays at u 3 even when the solutions u and u 2 eappea. It is only when = l that u 3 disappeas that the system jumps to the u -equilibium. h l u u 3 u l h.3 Delay Models One of the deficiencies of Eq. (2) is that the bith ate is assumed to act instantaneously wheeas thee is typically a time delay fom one geneation to the next because of the time to each matuity and the time fom conception to bith. One way to genealize Eq. () is to wite [ dn dt = N(t) N(t T) ], (3) K 4
6 whee the egulatoy effect depends on the population at an ealie time t T athe than the population at t. A moe easonable model would be to say that the effect should eally be an aveage ove past populations, e.g. dn [ dt = N(t) t ] w(t t )N(t )dt. K The weighting facto w will tend to zeo fo lage negative and positive t and will have a maximum at some time T. Fo w appoaching a delta function at time T, we will ecove Eq. (3). It tuns out that this model can show stable limit cycle peiodic solutions fo cetain values of the poduct T. This model has two steady states, N(t) = and N(t) = K. It is again convenient to use dimensionless units, u = N K, τ = t, T = T, which gives du = u(τ)[ u(τ T) ]..3. Linea stability analysis Using u = +n this becomes dn = [+n(τ)][ n(τ T)], and lineaizing, we get dn = n(τ T). Look fo solutions of the fom n(τ) = ce λτ λ = e λ T. The condition fo the solution to be unstable is Reλ >. That is clealy not possible if λ is eal-valued. We instead take λ = µ+iω and conside the eal and imaginay pats of λ = e λ T: µ = e µ T cosω T, ω = e µ T sinω T. (4) 5
7 We now equie µ >. Also note that if ω is a solution to Eq. (4) then so is ω, so we can take ω > without loss of geneality. Fom the fist equation wegetcosω T < which leadstoω T > π/2. Multiplying thesecond equation by T the condition ω T > π/2 implies Te µ T sinω T > π/2. Since sinω T and e µ T < when µ >, this leads to T T > π/2. When the constant u(t) = (i.e. n(t) = ) is unstable, the altenative becomes an oscillatoy behavio a stable limit cycle and that happens to be the case in this model..3.2 Voluntay execise Ty a simulation of Eq. (3) with T = 2, K = and =.5,., and.5. Use t =., stat with the initial condition N(t) = +.(t T/2), t T. Confim that the petubation fom N = K descibed above dies out when the bith ate is taken to be =.5 (i.e. T = ) but gows to stable oscillations fo the lage. Also detemine the peiod of the oscillations which is expected to be close to 4T fo T not too fa above π/2. These equations may also be studied with values close to π/2, e.g. T =.99π/2 and T =.π/2, but the simulations then have to be un substantially longe to see the espective t behavios. 6
8 2 Continuous Models fo Inteacting Populations Thegenealpictueisthattheeisawholewebofinteactingspecies. Wewill hee go fo the simplest possible situation and conside only two inteacting species. We ae mostly inteested in the pedato-pey situation whee one species feeds on the othe. 2. The Lotka-Voltea model To be concete we will conside abbits and foxes. The assumptions of the model ae: In the absence of foxes the abbits beed with dr/dt = ar, i.e. an exponential gowth, without bounds. The pesence of foxes means that the gowth ate is educed by a tem bfr. The elative gowth ate of the fox population is popotional to the population of abbits, df/dt = cfr. The foxes have a cetain death ate, df. Taken togethe this becomes dr dt df dt = ar bfr, = cfr df. This is the Lotka-Voltea system of equations. It was used by Voltea in 926 as a simple model to explain the oscillatoy levels of cetain fish catches. The same equations had aleady in 92 been deived by Lotka to descibe a hypotetical chemical eaction. 2.2 Rewite in dimensionless units To simplify the analysis we ewite the equations with f = bf a, = cr d, τ = at, α = d a, 7
9 and they become d df It is easy to see that this gives two steady states: 2.3 Phase space plot = ( f), (5) = αf( ). (6) (,f) = (,), (,f) = (,). The figue to the ight shows one possible solution fo the Lotka-Voltea equations with the paamete α =. This is fo the initial values () = 2.5 and f() =.. This shows oscillatoy behavio as a consequence of the coupling between the two species. A diffeent kind of figue a phase spaceplot isshown totheight. That is obtained by plotting the two vaiables hee the non-dimensional densities of foxes and abbits against one anothe. That kind of plot gives mostly the same infomation (though the time scale is eliminated) but as we will see it is vey convenient fo analyzing these kinds of systems. f f. τ Even though the Lotka-Voltea equations give the desied oscillations it has a seious dawbacks the behavio is vey sensitive to petubations it theefoe mainly woks as a stating point fo developing moe ealistic models. 8
10 The sensitivity to petubations is clealy seen in egions whee the cuves ae close togethe. E.g. aound (,f) = (,.2) in this figue, a small decease in f takes us to a diffeent cuve with athe diffeent popeties. We now tun to two diffeent methods that help us detemine popeties of systems of coupled diffeential equations f without finding the exact solution. These methods ae () phase space analyses and (2) stability analyses of the steady states. The latte becomes somewhat moe involved fo the case of two coupled vaiables than what we saw in the pevious chapte. 2.4 Gaphical phase space analysis Fom the Lotka-Voltea equations we see that the ate of change of and f is only a function of these same vaiables and this suggests plotting ( d, df ) as a vecto field. (We have hee suppessed the infomation of the magnitude of these vectos, only keeping the diection of the field at each point.) The figues below show these vecto fields fo both the odinay Lotka- Voltea equations to the left and some modified equations to the ight. By following the vectos it is possible to constuct the possible tajectoies. 2 2 f f Besides the vectos the ight figue also shows the nullclines which ae the cuves in the (,f)-plane whee d/ = o df/ =. The points whee two nullclines meet ae the steady states, and the analysis of thei popeties is impotant fo detemining the behavio of the system. 9
11 2.5 Linea analysis of the steady state Fo the case of only one vaiable we lineaize by witing f(x) f(x )+(x x )f (x ), andthis is especially useful when f(x ) = when the behavio is govened by the deivative f (x ). Fo two vaiables we wite f = (f,f 2 ) and x = (x,x 2 ) and the coesponding expession become f(x) f(x )+(x x ) f, (7) x x whee the ightmost symbol is a way to wite the Jacobian fo f = x x ( f x f 2 x 2 f 2 f x x 2 ) x = A. (8) To apply this to ou pesent situation we use the compact notation d df dx = f(x), = f (,f), = f 2 (,f). (Hee the symbol f unfotunately has two diffeent meanings, as it is both the function and the density of foxes. The meaning should howeve hopefully be clea fom the context.) Let x denote a steady state, i.e. a state with f(x ) =. It then follows that the dynamics may be witten dx f(x )+A(x x ) = A x, whee x is the deviation fom the steady state x.
12 2.6 Stable o unstable In the one-dimensional case we found that the stability of the steady state was govened by f (x ). In the pesent case it is instead the eigenvalues of the matix that detemine the behavio. Fo the following, we change the notation and let x denote the deviation fom the fixed point. Thee ae then eigenvectos x ν and eigenvalues λ ν such that Ax ν = λ ν x ν. An abitay initial state may now be witten in tems of the eigenvectos as 2 x() = c ν ()x ν, ν= and the time dependence then becomes dx = Ax = ν c ν Ax ν = ν c ν λ ν x ν. The solution becomes x(τ) = ν c ν (τ)x ν, whee c ν (τ) = c ν ()e λντ. 2.7 Eigenvalues The eigenvalues ae detemined fom the equation With this becomes which gives det(a λi) =. A = ( a a 2 a 2 a 22 (a λ)(a 22 λ) a 2 a 2 =, a a 22 a 2 a 2 λ(a +a 22 )+λ 2 =, and, with TA = a +a 22, we find λ ± = ( ) TA± (TA) 2 4 deta. (9) 2 ),
13 This shows that the behavio is detemined by TA and deta. Thee ae then seveal diffeent possible behavios. A complete list is given in Appendix A of Muay, Mathematical Biology. We hee just focus on the cases whee the eigenvalues ae complex, which is what we get if 4 deta > (TA) 2. (The equation in Fig. A.2 appeas to be wong.) In the phase space plot the solutions ae ciculating aound the fixed point in one of the following ways: A stable spial towads to fixed point, if TA <. An ellipse centeed at the fixed point, if TA =. An unstable spial out fom the fixed point, if TA >. 2.8 Stability analysis of the Lotka-Voltea equations We now wanty to apply this stability analysis to the Lotka-Voltea equations. Ou equations ae d = f (,f) = f, df = f 2 (,f) = α(f f), and the Jacobian becomes ( ) ( ) f / f A = / f f =. f 2 / f 2 / f αf α( ) This should now be evaluated at the espective fixed points The tivial fixed point (,f) = (,) At (,f) = (,) we get A = ( α TA = α and deta = α. The eigenvalues ae then λ ± = ( ) α± ( α) 2 +4α = 2 2 ( α±(+α)), which gives ) λ + =, λ = α, and since one of these eigenvalues has the eal component >, x(τ) gows exponentially. 2,
14 2.8.2 The non-tivial fixed point (,f) = (,) At (,f) = (,), which is the othe fixed point, we get A = ( α and TA = and deta = α. The eigenvalues ae and the solution is λ ± = ±i α, x(τ) = c x e i ατ +c 2 x 2 e i ατ, which shows that the solution close to the neighbohood of (,) is peiodic with peiod 2π/ α. Since Reλ = fo both eigenvalues the steady state is neutally stable. 2.9 Realistic Pedato-Pey Models One of the unealistic assumptions in the Lotka-Voltea model is that gowth of abbits is unbounded in the absence of pedation. The othe is that thee is no limit to the pey consumption. The task in the compute lab is to modify the equations in diffeent ways to see how the behavio changes. ), 3
15 2. Realistic Pedato-Pey Models One of the unealistic assumptions in the Lotka-Voltea model is that gowth of abbits is unbounded in the absence of pedation. The othe is that thee is no limit to the pey consumption. Define modified equations: d df = f (,f) = ( g) f +s, = f 2 (,f) = α f +s αf. To detemine the non-tivial fixed point which we denote by (,f ) we use df/ = in the second equation and get /( + s ) = which gives = /( s). In the fist equation d/ = gives f = ( g )(+s ). To detemine the eigenvalues we fist get the following esult d d f +s = f +s s f (+s) = 2 f (+s) 2, which we use in Eq. (8) to get A = 2g f (+s) 2 α f (+s) 2 +s α +s α, () 2.. How does the behaviou depend on g and s? The next step to analyze this model is to get out the eigenvalues and detemine fo what paamete values diffeent kinds of behavious should be expected. To detemine the eigenvalues we would need to put in the expessions fo and f in Eq. () and then use the diffeent elements in A in Eq. (9). This quickly becomes athe messy and we theefoe tun to a gaphical method to detemine the egion in paamete space hee the new paametes g and s whee a cetain behavio would be expected. The method consists of calculating the popeties fo many diffeent paametesatthesametime. Inmatlab(osome othe plotting pogam) one defines two vectos with values of g and s between cetainlimits. Onethencalculates, f, g = /Rgass det A > (ta) 2 T A > s = /R sat
16 TA, and deta fo each of these paamete paisandfinally plotsthepoints(g,s) which fulfill cetain conditions. The figue shows two such egions. The solid dots ae the points fo which 4 deta > (TA) 2, and thus give complex eigenvalues. The geen squaes ae the points fo which TA > and whee we expect the fixed point to be unstable. In the egion whee both these conditions ae fulfilled thee is a possibility of a stable limit cycle. The following figues ae the esults obtained with g =.25 and s =.5. One of the cuves stats athe close to the(unstable) fixed point and spials outwad until it stabilizes as it appoaches the stable limit cycle. The othe cuve stats outside the limit cycle and spials inwads. The figues below show the time dependence fo both and f, stating fom the diffeent stating points, but eventually eaching the same behavio f τ f f τ 5
17 3 Discete Population Models fo a Single Species 3. Intoduction: simple models Many species have no ovelap between successive geneations, population gowth take place though discete steps: N t+ = N t F(N t ) = f(n t ). The simplest case: suppose that F(N t ) = > : N t+ = N t N t = t N. () This is usually not vey ealistic, though it could wok fo the ealy stages of gowth of cetain bacteia. 3.2 Moe geneal One geneally expects the function f(n t ) to have a maximum at N m and decease fo N t > N m. One such model is ( N t+ = N t N ) t, K but has the dawback that N t+ < if N t > K, which is of couse not ealistic. Anothe bette model is: [ ( N t+ = N t exp N t K )], >, K >, which is a modification of Eq. () with a motality facto 3.3 Gaphical solution exp( N t /K). 6
18 Giventhefunctionf(N t )andthestating value N the sequence N, N 2,..., may be detemined gaphically: The steady state solutions N ae intesections of the cuve N t+ = f(n t ) and the staight line N t+ = N t. If N is a stable o unstable solution is detemined by the eigenvalue of the system f (N t ), since that detemine the effect of a small petubation fom the steady state: f(n +δ) = N +δf (N ). A small petubation will lead to oscillations if f (N ) <. A small petubation will die out if f (N ) <. 3.4 Discete Logistic Model and Chaos Rescale with u t = N t /K andexamine the behavio of u t+ = u t ( u t ), >. We ae inteested in solutions with u t >, only. The steady states and the coesponding eigenvalues λ = f (u ) ae u =, λ =, u 2 =, λ 2 = 2. ut The cuves hee ae fo =,.5, The solid dots show the espective steady state solutions, u Stable steady state The detemining facto fo the stability is the eigenvalue, λ = f (u 2) = 2 ; the system is stable fo λ <. The figues hee show the behavio fo u t 7
19 = 2.8 and the initial value u =.72. The successive iteations take u t towads the stable solution u = ( )/.643. ut+ ut u t 2 t Unstable steady state The value = 3.2 on the othe hand gives λ =.2 and we expect the solution u to be unstable. Stating at u =.72, which is close to u 2 = ( )/.688, we see, just as expected, that successive iteations take us away fom the steady state. As the figues below show this gives a new behavio with an oscillation with peiod=2. ut+ ut u t 2 t Peiod doubling To examine this behavio in moe detail we conside the map fom u t to u t+2 defined by C u t+ = u t ( u t ), u t+2 = u t+ ( u t+ ). ut+2, ut+ A B The figue to the ight is fo = 3.2 and we see that thee ae then thee solutions to u t+2 = u t, which we denote by u A, u t 8
20 u B, and u C. Of these u B is an unstable solution since it has f (u B ) >. (Note that u B = u 2 ). On the othe hand, u A and u C ae stable. Note what this means: u t = u A will give u t+2 = u A, wheeas u t = u C similaly gives u t+2 = u C. Togethe with the ealie figue we may conclude that that system oscillates, u A, u C, u A, u C,... This phenomenon which is hee seen as a change fom a simple steady state to an oscillation between two diffeent values is called a bifucation. The discussion above may now be taken one moe step in the same diection as shown in this figue, which is C fo = 3.5: when inceases the solutions u A and u C move somewhat, and the B A steady states eventually become unstable as f (u C ) < and f (u A ) <. When that happens we again get a peiod doubling and the system epeats itself afte u t 4 units of time. It is then possible to epeat the discussion in this section fo u t+4. We would then find fou solutions to u t+4 = u t which ae stable fo = 3.5 but tun unstable at a slightly lage. This instability gives an oscillation of peiod eight, and it now seems that this peiod doubling may be continued without limit. It tuns out that the distance in between successive bifucations deceases apidly, and as we appoach c the oscillations have peiod 2 k with k. Fo > c the behavio is apeiodic and we thee ente the egion of chaos Chaos It tuns out that the behavio changes in unexpected ways as a function of : Fo < 3 thee is a unique solution ( )/. Fo 3 < + 6( 3.45) the system has peiodic fluctuates between two values. Fo + 6 < < 3.54 (appoximately) the system has peiodic oscillations between fou values. 9 ut+2, ut+
21 Fo 3.54 < < 3.57 the system oscillates between 8, 6, 32, values, etc. At 3.57 is the onset of chaos. We can no longe see any oscillations of finite peiod and slight vaiations in the initial value yields damatically diffeent esults ove time. ut
MATH 415, WEEK 3: Parameter-Dependence and Bifurcations
MATH 415, WEEK 3: Paamete-Dependence and Bifucations 1 A Note on Paamete Dependence We should pause to make a bief note about the ole played in the study of dynamical systems by the system s paametes.
More informationF-IF Logistic Growth Model, Abstract Version
F-IF Logistic Gowth Model, Abstact Vesion Alignments to Content Standads: F-IFB4 Task An impotant example of a model often used in biology o ecology to model population gowth is called the logistic gowth
More informationMagnetic Field. Conference 6. Physics 102 General Physics II
Physics 102 Confeence 6 Magnetic Field Confeence 6 Physics 102 Geneal Physics II Monday, Mach 3d, 2014 6.1 Quiz Poblem 6.1 Think about the magnetic field associated with an infinite, cuent caying wie.
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More informationHandout: IS/LM Model
Econ 32 - IS/L odel Notes Handout: IS/L odel IS Cuve Deivation Figue 4-4 in the textbook explains one deivation of the IS cuve. This deivation uses the Induced Savings Function fom Chapte 3. Hee, I descibe
More informationBifurcation Analysis for the Delay Logistic Equation with Two Delays
IOSR Jounal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue 5 Ve. IV (Sep. - Oct. 05), PP 53-58 www.iosjounals.og Bifucation Analysis fo the Delay Logistic Equation with Two Delays
More informationCHAPTER 3. Section 1. Modeling Population Growth
CHAPTER 3 Section 1. Modeling Population Gowth 1.1. The equation of the Malthusian model is Pt) = Ce t. Apply the initial condition P) = 1. Then 1 = Ce,oC = 1. Next apply the condition P1) = 3. Then 3
More informationRotational Motion. Lecture 6. Chapter 4. Physics I. Course website:
Lectue 6 Chapte 4 Physics I Rotational Motion Couse website: http://faculty.uml.edu/andiy_danylov/teaching/physicsi Today we ae going to discuss: Chapte 4: Unifom Cicula Motion: Section 4.4 Nonunifom Cicula
More information3.1 Random variables
3 Chapte III Random Vaiables 3 Random vaiables A sample space S may be difficult to descibe if the elements of S ae not numbes discuss how we can use a ule by which an element s of S may be associated
More informationA Crash Course in (2 2) Matrices
A Cash Couse in ( ) Matices Seveal weeks woth of matix algeba in an hou (Relax, we will only stuy the simplest case, that of matices) Review topics: What is a matix (pl matices)? A matix is a ectangula
More informationPhysics 2B Chapter 22 Notes - Magnetic Field Spring 2018
Physics B Chapte Notes - Magnetic Field Sping 018 Magnetic Field fom a Long Staight Cuent-Caying Wie In Chapte 11 we looked at Isaac Newton s Law of Gavitation, which established that a gavitational field
More informationEM Boundary Value Problems
EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do
More information( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.
9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can
More informationAbsorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere
Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in
More information6 PROBABILITY GENERATING FUNCTIONS
6 PROBABILITY GENERATING FUNCTIONS Cetain deivations pesented in this couse have been somewhat heavy on algeba. Fo example, detemining the expectation of the Binomial distibution (page 5.1 tuned out to
More informationChapter 9 Dynamic stability analysis III Lateral motion (Lectures 33 and 34)
Pof. E.G. Tulapukaa Stability and contol Chapte 9 Dynamic stability analysis Lateal motion (Lectues 33 and 34) Keywods : Lateal dynamic stability - state vaiable fom of equations, chaacteistic equation
More informationScattering in Three Dimensions
Scatteing in Thee Dimensions Scatteing expeiments ae an impotant souce of infomation about quantum systems, anging in enegy fom vey low enegy chemical eactions to the highest possible enegies at the LHC.
More informationNon-Linear Dynamics Homework Solutions Week 2
Non-Linea Dynamics Homewok Solutions Week Chis Small Mach, 7 Please email me at smach9@evegeen.edu with any questions o concens eguading these solutions. Fo the ececises fom section., we sketch all qualitatively
More informationASTR415: Problem Set #6
ASTR45: Poblem Set #6 Cuan D. Muhlbege Univesity of Mayland (Dated: May 7, 27) Using existing implementations of the leapfog and Runge-Kutta methods fo solving coupled odinay diffeential equations, seveal
More informationworking pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50
woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,
More informationSurveillance Points in High Dimensional Spaces
Société de Calcul Mathématique SA Tools fo decision help since 995 Suveillance Points in High Dimensional Spaces by Benad Beauzamy Januay 06 Abstact Let us conside any compute softwae, elying upon a lage
More informationSAMPLING DELAY AND BACKLASH IN BALANCING SYSTEMS
PERIODICA POLYTECHNICA SER. MECH. ENG. VOL. 44, NO., PP. 77 84 () SAMPLING DELAY AND BACKLASH IN BALANCING SYSTEMS László E. KOLLÁR, Gáo STÉPÁN and S. John HOGAN Depatment of Applied Mechanics Technical
More informationNumerical Integration
MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.
More informationPROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.
POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and
More informationChapter 3 Optical Systems with Annular Pupils
Chapte 3 Optical Systems with Annula Pupils 3 INTRODUCTION In this chapte, we discuss the imaging popeties of a system with an annula pupil in a manne simila to those fo a system with a cicula pupil The
More informationThe Substring Search Problem
The Substing Seach Poblem One algoithm which is used in a vaiety of applications is the family of substing seach algoithms. These algoithms allow a use to detemine if, given two chaacte stings, one is
More informationChapter 7-8 Rotational Motion
Chapte 7-8 Rotational Motion What is a Rigid Body? Rotational Kinematics Angula Velocity ω and Acceleation α Unifom Rotational Motion: Kinematics Unifom Cicula Motion: Kinematics and Dynamics The Toque,
More informationEFFECTS OF FRINGING FIELDS ON SINGLE PARTICLE DYNAMICS. M. Bassetti and C. Biscari INFN-LNF, CP 13, Frascati (RM), Italy
Fascati Physics Seies Vol. X (998), pp. 47-54 4 th Advanced ICFA Beam Dynamics Wokshop, Fascati, Oct. -5, 997 EFFECTS OF FRININ FIELDS ON SINLE PARTICLE DYNAMICS M. Bassetti and C. Biscai INFN-LNF, CP
More informationComputational Methods of Solid Mechanics. Project report
Computational Methods of Solid Mechanics Poject epot Due on Dec. 6, 25 Pof. Allan F. Bowe Weilin Deng Simulation of adhesive contact with molecula potential Poject desciption In the poject, we will investigate
More informationWhen two numbers are written as the product of their prime factors, they are in factored form.
10 1 Study Guide Pages 420 425 Factos Because 3 4 12, we say that 3 and 4 ae factos of 12. In othe wods, factos ae the numbes you multiply to get a poduct. Since 2 6 12, 2 and 6 ae also factos of 12. The
More informationPhysics 2A Chapter 10 - Moment of Inertia Fall 2018
Physics Chapte 0 - oment of netia Fall 08 The moment of inetia of a otating object is a measue of its otational inetia in the same way that the mass of an object is a measue of its inetia fo linea motion.
More informationInternet Appendix for A Bayesian Approach to Real Options: The Case of Distinguishing Between Temporary and Permanent Shocks
Intenet Appendix fo A Bayesian Appoach to Real Options: The Case of Distinguishing Between Tempoay and Pemanent Shocks Steven R. Genadie Gaduate School of Business, Stanfod Univesity Andey Malenko Gaduate
More informationOn the Comparison of Stability Analysis with Phase Portrait for a Discrete Prey-Predator System
Intenational Jounal of Applied Engineeing Reseach ISSN 0973-4562 Volume 12, Nume 24 (2017) pp. 15273-15277 Reseach India Pulications. http://www.ipulication.com On the Compaison of Staility Analysis with
More informationExplosive Contagion in Networks (Supplementary Information)
Eplosive Contagion in Netwoks (Supplementay Infomation) Jesús Gómez-Gadeñes,, Laua Loteo, Segei N. Taaskin, and Fancisco J. Péez-Reche Institute fo Biocomputation and Physics of Comple Systems (BIFI),
More informationTo Feel a Force Chapter 7 Static equilibrium - torque and friction
To eel a oce Chapte 7 Chapte 7: Static fiction, toque and static equilibium A. Review of foce vectos Between the eath and a small mass, gavitational foces of equal magnitude and opposite diection act on
More informationOSCILLATIONS AND GRAVITATION
1. SIMPLE HARMONIC MOTION Simple hamonic motion is any motion that is equivalent to a single component of unifom cicula motion. In this situation the velocity is always geatest in the middle of the motion,
More informationON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0},
ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION E. J. IONASCU and A. A. STANCU Abstact. We ae inteested in constucting concete independent events in puely atomic pobability
More information7.2. Coulomb s Law. The Electric Force
Coulomb s aw Recall that chaged objects attact some objects and epel othes at a distance, without making any contact with those objects Electic foce,, o the foce acting between two chaged objects, is somewhat
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II September 15, 2012 Prof. Alan Guth PROBLEM SET 2
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Depatment Physics 8.07: Electomagnetism II Septembe 5, 202 Pof. Alan Guth PROBLEM SET 2 DUE DATE: Monday, Septembe 24, 202. Eithe hand it in at the lectue,
More informationNuclear and Particle Physics - Lecture 20 The shell model
1 Intoduction Nuclea and Paticle Physics - Lectue 0 The shell model It is appaent that the semi-empiical mass fomula does a good job of descibing tends but not the non-smooth behaviou of the binding enegy.
More informationPulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects
Pulse Neuton Neuton (PNN) tool logging fo poosity Some theoetical aspects Intoduction Pehaps the most citicism of Pulse Neuton Neuon (PNN) logging methods has been chage that PNN is to sensitive to the
More informationHomework 7 Solutions
Homewok 7 olutions Phys 4 Octobe 3, 208. Let s talk about a space monkey. As the space monkey is oiginally obiting in a cicula obit and is massive, its tajectoy satisfies m mon 2 G m mon + L 2 2m mon 2
More informationequilibrium in the money market
Sahoko KAJI --- Open Economy Macoeconomics ectue Notes I I A Review of Closed Economy Macoeconomics We begin by eviewing some of the basics of closed economy macoeconomics that ae indispensable in undestanding
More information-Δ u = λ u. u(x,y) = u 1. (x) u 2. (y) u(r,θ) = R(r) Θ(θ) Δu = 2 u + 2 u. r = x 2 + y 2. tan(θ) = y/x. r cos(θ) = cos(θ) r.
The Laplace opeato in pola coodinates We now conside the Laplace opeato with Diichlet bounday conditions on a cicula egion Ω {(x,y) x + y A }. Ou goal is to compute eigenvalues and eigenfunctions of the
More informationPES 3950/PHYS 6950: Homework Assignment 6
PES 3950/PHYS 6950: Homewok Assignment 6 Handed out: Monday Apil 7 Due in: Wednesday May 6, at the stat of class at 3:05 pm shap Show all woking and easoning to eceive full points. Question 1 [5 points]
More informationME 3600 Control Systems Frequency Domain Analysis
ME 3600 Contol Systems Fequency Domain Analysis The fequency esponse of a system is defined as the steady-state esponse of the system to a sinusoidal (hamonic) input. Fo linea systems, the esulting steady-state
More informationAppendix 2. Equilibria, trophic indices, and stability of tri-trophic model with dynamic stoichiometry of plants.
OIKO O15875 all. R. huin J. B. Diehl. and NisbetR. M. 2007. Food quality nutient limitation of seconday poduction and the stength of tophic cascades. Oikos 000: 000 000. Appendix 2. Equilibia tophic indices
More informationGeometry of the homogeneous and isotropic spaces
Geomety of the homogeneous and isotopic spaces H. Sonoda Septembe 2000; last evised Octobe 2009 Abstact We summaize the aspects of the geomety of the homogeneous and isotopic spaces which ae most elevant
More informationPhysics 221 Lecture 41 Nonlinear Absorption and Refraction
Physics 221 Lectue 41 Nonlinea Absoption and Refaction Refeences Meye-Aendt, pp. 97-98. Boyd, Nonlinea Optics, 1.4 Yaiv, Optical Waves in Cystals, p. 22 (Table of cystal symmeties) 1. Intoductoy Remaks.
More informationSection 8.2 Polar Coordinates
Section 8. Pola Coodinates 467 Section 8. Pola Coodinates The coodinate system we ae most familia with is called the Catesian coodinate system, a ectangula plane divided into fou quadants by the hoizontal
More informationA New Approach to General Relativity
Apeion, Vol. 14, No. 3, July 7 7 A New Appoach to Geneal Relativity Ali Rıza Şahin Gaziosmanpaşa, Istanbul Tukey E-mail: aizasahin@gmail.com Hee we pesent a new point of view fo geneal elativity and/o
More information763620SS STATISTICAL PHYSICS Solutions 2 Autumn 2012
763620SS STATISTICAL PHYSICS Solutions 2 Autumn 2012 1. Continuous Random Walk Conside a continuous one-dimensional andom walk. Let w(s i ds i be the pobability that the length of the i th displacement
More informationPHYS 705: Classical Mechanics. Small Oscillations
PHYS 705: Classical Mechanics Small Oscillations Fomulation of the Poblem Assumptions: V q - A consevative system with depending on position only - The tansfomation equation defining does not dep on time
More informationPhysics 161 Fall 2011 Extra Credit 2 Investigating Black Holes - Solutions The Following is Worth 50 Points!!!
Physics 161 Fall 011 Exta Cedit Investigating Black Holes - olutions The Following is Woth 50 Points!!! This exta cedit assignment will investigate vaious popeties of black holes that we didn t have time
More informationMath 124B February 02, 2012
Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial
More informationMAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS
The 8 th Intenational Confeence of the Slovenian Society fo Non-Destuctive Testing»pplication of Contempoay Non-Destuctive Testing in Engineeing«Septembe 1-3, 5, Potoož, Slovenia, pp. 17-1 MGNETIC FIELD
More informationAnalysis of simple branching trees with TI-92
Analysis of simple banching tees with TI-9 Dušan Pagon, Univesity of Maibo, Slovenia Abstact. In the complex plane we stat at the cente of the coodinate system with a vetical segment of the length one
More informationPartition Functions. Chris Clark July 18, 2006
Patition Functions Chis Clak July 18, 2006 1 Intoduction Patition functions ae useful because it is easy to deive expectation values of paametes of the system fom them. Below is a list of the mao examples.
More informationPhysics 235 Chapter 5. Chapter 5 Gravitation
Chapte 5 Gavitation In this Chapte we will eview the popeties of the gavitational foce. The gavitational foce has been discussed in geat detail in you intoductoy physics couses, and we will pimaily focus
More informationLecture 8 - Gauss s Law
Lectue 8 - Gauss s Law A Puzzle... Example Calculate the potential enegy, pe ion, fo an infinite 1D ionic cystal with sepaation a; that is, a ow of equally spaced chages of magnitude e and altenating sign.
More informationOn the integration of the equations of hydrodynamics
Uebe die Integation de hydodynamischen Gleichungen J f eine u angew Math 56 (859) -0 On the integation of the equations of hydodynamics (By A Clebsch at Calsuhe) Tanslated by D H Delphenich In a pevious
More informationMath 301: The Erdős-Stone-Simonovitz Theorem and Extremal Numbers for Bipartite Graphs
Math 30: The Edős-Stone-Simonovitz Theoem and Extemal Numbes fo Bipatite Gaphs May Radcliffe The Edős-Stone-Simonovitz Theoem Recall, in class we poved Tuán s Gaph Theoem, namely Theoem Tuán s Theoem Let
More informationMotion in One Dimension
Motion in One Dimension Intoduction: In this lab, you will investigate the motion of a olling cat as it tavels in a staight line. Although this setup may seem ovesimplified, you will soon see that a detailed
More informationChem 453/544 Fall /08/03. Exam #1 Solutions
Chem 453/544 Fall 3 /8/3 Exam # Solutions. ( points) Use the genealized compessibility diagam povided on the last page to estimate ove what ange of pessues A at oom tempeatue confoms to the ideal gas law
More informationThe Poisson bracket and magnetic monopoles
FYST420 Advanced electodynamics Olli Aleksante Koskivaaa Final poject ollikoskivaaa@gmail.com The Poisson backet and magnetic monopoles Abstact: In this wok magnetic monopoles ae studied using the Poisson
More informationAs is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.
Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.
More informationMATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form
MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE ANDRAS VASY We conside second ode constant coefficient scala linea PDEs on R n. These have the fom Lu = f L = a ij xi xj + b i xi + c i whee a ij b i and
More informationVector d is a linear vector function of vector d when the following relationships hold:
Appendix 4 Dyadic Analysis DEFINITION ecto d is a linea vecto function of vecto d when the following elationships hold: d x = a xxd x + a xy d y + a xz d z d y = a yxd x + a yy d y + a yz d z d z = a zxd
More information2 Governing Equations
2 Govening Equations This chapte develops the govening equations of motion fo a homogeneous isotopic elastic solid, using the linea thee-dimensional theoy of elasticity in cylindical coodinates. At fist,
More informationChapter 2: Introduction to Implicit Equations
Habeman MTH 11 Section V: Paametic and Implicit Equations Chapte : Intoduction to Implicit Equations When we descibe cuves on the coodinate plane with algebaic equations, we can define the elationship
More informationEncapsulation theory: radial encapsulation. Edmund Kirwan *
Encapsulation theoy: adial encapsulation. Edmund Kiwan * www.edmundkiwan.com Abstact This pape intoduces the concept of adial encapsulation, wheeby dependencies ae constained to act fom subsets towads
More informationDo Managers Do Good With Other People s Money? Online Appendix
Do Manages Do Good With Othe People s Money? Online Appendix Ing-Haw Cheng Haison Hong Kelly Shue Abstact This is the Online Appendix fo Cheng, Hong and Shue 2013) containing details of the model. Datmouth
More informationA Bijective Approach to the Permutational Power of a Priority Queue
A Bijective Appoach to the Pemutational Powe of a Pioity Queue Ia M. Gessel Kuang-Yeh Wang Depatment of Mathematics Bandeis Univesity Waltham, MA 02254-9110 Abstact A pioity queue tansfoms an input pemutation
More informationMath Notes on Kepler s first law 1. r(t) kp(t)
Math 7 - Notes on Keple s fist law Planetay motion and Keple s Laws We conside the motion of a single planet about the sun; fo simplicity, we assign coodinates in R 3 so that the position of the sun is
More informationFall 2016 Semester METR 3113 Atmospheric Dynamics I: Introduction to Atmospheric Kinematics and Dynamics
Fall 06 Semeste METR 33 Atmospheic Dynamics I: Intoduction to Atmospheic Kinematics Dynamics Lectue 7 Octobe 3 06 Topics: Scale analysis of the equations of hoizontal motion Geostophic appoximation eostophic
More informationBoundary Layers and Singular Perturbation Lectures 16 and 17 Boundary Layers and Singular Perturbation. x% 0 Ω0æ + Kx% 1 Ω0æ + ` : 0. (9.
Lectues 16 and 17 Bounday Layes and Singula Petubation A Regula Petubation In some physical poblems, the solution is dependent on a paamete K. When the paamete K is vey small, it is natual to expect that
More informationExploration of the three-person duel
Exploation of the thee-peson duel Andy Paish 15 August 2006 1 The duel Pictue a duel: two shootes facing one anothe, taking tuns fiing at one anothe, each with a fixed pobability of hitting his opponent.
More informationHopefully Helpful Hints for Gauss s Law
Hopefully Helpful Hints fo Gauss s Law As befoe, thee ae things you need to know about Gauss s Law. In no paticula ode, they ae: a.) In the context of Gauss s Law, at a diffeential level, the electic flux
More information4/18/2005. Statistical Learning Theory
Statistical Leaning Theoy Statistical Leaning Theoy A model of supevised leaning consists of: a Envionment - Supplying a vecto x with a fixed but unknown pdf F x (x b Teache. It povides a desied esponse
More informationCompactly Supported Radial Basis Functions
Chapte 4 Compactly Suppoted Radial Basis Functions As we saw ealie, compactly suppoted functions Φ that ae tuly stictly conditionally positive definite of ode m > do not exist The compact suppot automatically
More information15 Solving the Laplace equation by Fourier method
5 Solving the Laplace equation by Fouie method I aleady intoduced two o thee dimensional heat equation, when I deived it, ecall that it taes the fom u t = α 2 u + F, (5.) whee u: [0, ) D R, D R is the
More informationAST 121S: The origin and evolution of the Universe. Introduction to Mathematical Handout 1
Please ead this fist... AST S: The oigin and evolution of the Univese Intoduction to Mathematical Handout This is an unusually long hand-out and one which uses in places mathematics that you may not be
More informationOn a quantity that is analogous to potential and a theorem that relates to it
Su une quantité analogue au potential et su un théoème y elatif C R Acad Sci 7 (87) 34-39 On a quantity that is analogous to potential and a theoem that elates to it By R CLAUSIUS Tanslated by D H Delphenich
More informationLecture 28: Convergence of Random Variables and Related Theorems
EE50: Pobability Foundations fo Electical Enginees July-Novembe 205 Lectue 28: Convegence of Random Vaiables and Related Theoems Lectue:. Kishna Jagannathan Scibe: Gopal, Sudhasan, Ajay, Swamy, Kolla An
More informationRigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018
Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining
More informationClassical Worm algorithms (WA)
Classical Wom algoithms (WA) WA was oiginally intoduced fo quantum statistical models by Pokof ev, Svistunov and Tupitsyn (997), and late genealized to classical models by Pokof ev and Svistunov (200).
More informationGraphs of Sine and Cosine Functions
Gaphs of Sine and Cosine Functions In pevious sections, we defined the tigonometic o cicula functions in tems of the movement of a point aound the cicumfeence of a unit cicle, o the angle fomed by the
More informationChapter 13 Gravitation
Chapte 13 Gavitation In this chapte we will exploe the following topics: -Newton s law of gavitation, which descibes the attactive foce between two point masses and its application to extended objects
More informationChapter 5 Linear Equations: Basic Theory and Practice
Chapte 5 inea Equations: Basic Theoy and actice In this chapte and the next, we ae inteested in the linea algebaic equation AX = b, (5-1) whee A is an m n matix, X is an n 1 vecto to be solved fo, and
More informationStress Intensity Factor
S 47 Factue Mechanics http://imechanicaog/node/7448 Zhigang Suo Stess Intensity Facto We have modeled a body by using the linea elastic theoy We have modeled a cack in the body by a flat plane, and the
More information! E da = 4πkQ enc, has E under the integral sign, so it is not ordinarily an
Physics 142 Electostatics 2 Page 1 Electostatics 2 Electicity is just oganized lightning. Geoge Calin A tick that sometimes woks: calculating E fom Gauss s law Gauss s law,! E da = 4πkQ enc, has E unde
More informationThe Strain Compatibility Equations in Polar Coordinates RAWB, Last Update 27/12/07
The Stain Compatibility Equations in Pola Coodinates RAWB Last Update 7//7 In D thee is just one compatibility equation. In D polas it is (Equ.) whee denotes the enineein shea (twice the tensoial shea)
More informationTheWaveandHelmholtzEquations
TheWaveandHelmholtzEquations Ramani Duaiswami The Univesity of Mayland, College Pak Febuay 3, 2006 Abstact CMSC828D notes (adapted fom mateial witten with Nail Gumeov). Wok in pogess 1 Acoustic Waves 1.1
More information1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More informationME 210 Applied Mathematics for Mechanical Engineers
Tangent and Ac Length of a Cuve The tangent to a cuve C at a point A on it is defined as the limiting position of the staight line L though A and B, as B appoaches A along the cuve as illustated in the
More informationRight-handed screw dislocation in an isotropic solid
Dislocation Mechanics Elastic Popeties of Isolated Dislocations Ou study of dislocations to this point has focused on thei geomety and thei ole in accommodating plastic defomation though thei motion. We
More informationLecture 2 Date:
Lectue 2 Date: 5.1.217 Definition of Some TL Paametes Examples of Tansmission Lines Tansmission Lines (contd.) Fo a lossless tansmission line the second ode diffeential equation fo phasos ae: LC 2 d I
More informationSMT 2013 Team Test Solutions February 2, 2013
1 Let f 1 (n) be the numbe of divisos that n has, and define f k (n) = f 1 (f k 1 (n)) Compute the smallest intege k such that f k (013 013 ) = Answe: 4 Solution: We know that 013 013 = 3 013 11 013 61
More informationAdvanced Problems of Lateral- Directional Dynamics!
Advanced Poblems of Lateal- Diectional Dynamics! Robet Stengel, Aicaft Flight Dynamics! MAE 331, 216 Leaning Objectives 4 th -ode dynamics! Steady-state esponse to contol! Tansfe functions! Fequency esponse!
More information2 x 8 2 x 2 SKILLS Determine whether the given value is a solution of the. equation. (a) x 2 (b) x 4. (a) x 2 (b) x 4 (a) x 4 (b) x 8
5 CHAPTER Fundamentals When solving equations that involve absolute values, we usually take cases. EXAMPLE An Absolute Value Equation Solve the equation 0 x 5 0 3. SOLUTION By the definition of absolute
More information