Lecture 5: Travelling Waves
|
|
- Silas Hancock
- 5 years ago
- Views:
Transcription
1 Computational Biology Group (CoBI), D-BSSE, ETHZ Lecture 5: Travelling Waves Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015
2 26. Oktober / 68 Contents 1 Introduction to Travelling Waves Definitions Propagating Wave Solutions Initial conditions Axisymmetric form of Fisher equation 2 Modeling Examples The calcium-stimulated-calcium release mechanism Calcium Waves on Amphibian Eggs Models for Somitogenesis 3 Wave Pinning
3 Introduction to Travelling Waves Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
4 26. Oktober / 68 Travelling Wave-like Phenomena in Biology Calcium Waves on Amphibian Eggs
5 26. Oktober / 68 Travelling Wave-like Phenomena in Biology Somitogenesis Movie by Andy Oates
6 26. Oktober / 68 Travelling Wave Reference Frame Travelling Wave travels without change of shape. if u(x, t) represents a travelling wave, the shape of the solution will be the same for all time speed of propagation is a constant, which we denote c. If we look at this wave in a travelling wave frame moving at speed c then this wave will appear stationary.
7 26. Oktober / 68 Travelling Wave in one spatial dimension In one spatial dimension, x, diffusion of a molecule with concentration c(x, t) can be described as du dt = D d 2 u = D u. (1) dx 2 In mathematical terms, u(x, t) is a travelling wave that moves at constant speed c in the positive x-direction, if u(x, t) = u(x ct) = u(z), z = x ct. (2) z is referred to as wave variable.
8 26. Oktober / 68 PDE Set of ODEs can be rewritten as PDE Set of ODEs u(x, t) = u(x ct) = u(z), u t u x = c du dz z = x ct = du dz. (3) Partial differential equations in x and t become sets of ordinary differential equations.
9 26. Oktober / 68 Biological Constraints on PDEs To be physically realistic u(z) has to be bounded for all z and non-negative.
10 26. Oktober / 68 Linear Parabolic PDEs It can be shown that there are no physically realistic travelling wave solutions for reaction-diffusion equations of the form if f (u) is linear. u t = f (u) + Du xx (4) In biological systems f (u) is typically non-linear.
11 Propagating Wave Solutions Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
12 Example: Logistic Growth Consider an ODE that describes logistic growth of the variable u, du dt = u(1 u). (5) The general solution with integration constant c is given by u(x, t) = c exp (t) 1 + c exp (t). (6) If we chose as initial conditions u(x, 0) = exp (x) (7) then we obtain u(x, 0) = c 1 + c = 1 c = exp ( x). (8) 1 + exp (x) Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
13 26. Oktober / 68 u(x, t) = can be written as c exp (t) 1 + c exp (t) u(x, t) = Using the wave variable z = t x, we then have u(z) = with c = exp ( x) (9) exp (t x) 1 + exp (t x). (10) exp (z) 1 + exp (z). (11)
14 26. Oktober / 68 Travelling Wave Solution u(z) = exp (z) 1 + exp (z) If z = x t = const then the shape does not change, i.e. if dx dt = 1. The shape of the wave thus does not change if one travels with the wave at speed 1. This wave depends on the initial conditions and is highly unstable.
15 26. Oktober / 68 Example: Fisher Equation Consider the Fisher equation u t = ku(1 u) + Du xx. (12) This equation can be non-dimensionalized by using 1/k as timescale and D/k as length scale. The non-dimensionalized Fisher equation reads u t = u(1 u) + u xx. (13)
16 26. Oktober / 68 Transformation to travelling wave coordinates Consider the non-dimensionalized Fisher equation u t = u(1 u) + u xx. (14) We write u(x, t) = u(x ct) = U(z), z = x ct, c 0. (15) and upon substitution into Eqn.(47) we obtain U + cu + U(1 U) = 0 (16)
17 26. Oktober / 68 PDE Set of ODEs U + cu + U(1 U) = 0 can be written as set of ODEs, U = V V = cv U(1 U). (17)
18 26. Oktober / 68 Phase Plane Analysis Determine Nullclines Determine Steady States Determine Stability of Steady States Determine Trajectories and Phase Vectors
19 26. Oktober / 68 Nullclines and Steady States U = V = 0 V = cv U(1 U) = 0. (18) Nullclines: The U-nullcline is given by V = 0 The V-nullcline is given by V = U(1 U) c Steady states: (U,V) = (0,0) and (U, V) = (1,0).
20 26. Oktober / 68 Stability of Steady States To determine the stability of the steady states we linearize the set of equations at these steady states (U, V ) and determine the eigenvalues of the Jacobian, J = ( 0 1 (1 2U ) c ). (19) as λ ± = 1 2 ( ) tr(j) ± tr(j) 2 4det(J) (20)
21 26. Oktober / 68 Stability of Steady States λ ± = 1 2 ( ) tr(j) ± tr(j) 2 4det(J) Steady States: 1 (U, V ) = (1, 0): saddle node (1 positive and 1 negative eigenvalue) 2 (U, V ) = (0, 0): stable node (2 real, negative eigenvalues) if c 2 spiral (2 conjugate complex eigenvalues with negative real part) if c < 2.
22 26. Oktober / 68 Phase Plane Fisher Equation: U = V V = cv U(1 U). Steady States: 1. (U, V ) = (1, 0): saddle node 2. (U, V ) = (0, 0): stable stst Phase Plane Trajectory: dv du = cv U(1 U) V
23 26. Oktober / 68 Phase Plane Since U can assume negative values if there is a stable spiral at (U, V ) = (0, 0) we require c 2 for physically realistic travelling wave solutions.
24 26. Oktober / 68 For travelling waves to exist we require one stable node and one saddle node. A key question at this stage is: What kind of initial conditions u(x, 0) will evolve to a travelling wave solution? And, if such solution exists, what is its wave speed c?
25 26. Oktober / 68 Initial conditions & Wave Speed Kolmogoroff et al (1937) proved that if u(x, 0) has compact support, u(x, 0) = u 0 (x) 0, u 0 (x) = { 1 if x x1 0 if x x 2 (21) where x 1 < x 2 and u 0 (x) is continuous in x 1 < x < x 2, then the solution u(x, t) of the Fisher Equation evolves to a travelling wave solutions with c = c min = 2.
26 More general Initial Conditions For other initial data the solution critically depends on the behaviour of u(x, 0) as x ±. To see this, consider the leading edge of the evolving wave where since u is small u 2 << u such that Consider now as initial condition We then have if u t u + u xx. (22) u(x, 0) A exp ( αx); α > 0 as x. (23) u(x, t) = A exp ( α(x ct)) (24) αcu = u + α 2 u c = 1 α + α (25) Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
27 26. Oktober / 68 Wave Speed depends on Initial Conditions For initial condition u(x, 0) A exp ( αx) as x the wave speed c = 1 α + α depends on the initial conditions α > 0.
28 26. Oktober / 68 General Equation u t = f (u) + u xx u(x, 0) A exp ( αx), α > 0 (26) can be linearized as u t = f (u)u + u xx u(x, 0) A exp ( αx) (27) with solution if i.e. u(x, t) = A exp ( α(x ct)) (28) αcu = f (u)u + α 2 u (29) c = f (u) α + α (30)
29 26. Oktober / 68 The minimal travelling wave speed The minimal wave speed of can be determined as by differentiating c with respect to α c = f (u) α + α (31) c min = 2 f (u) (32) dc dα = 1 α 2 f (u) + 1 = 0 α = ± f (u) (33)
30 26. Oktober / 68 The travelling wave solution u t = f (u) + u xx u(x, 0) A exp ( αx), α > 0 (34) with f (u) having two zeros u 1, u 2 > u 1 If f (u 1 ) > 0 and f (u 2 ) < 0 wavefront solutions monotonically evolve with u going monotonically from u 1 to u 2 with wave speed c c min = 2 f (u) (35)
31 The flatter the wave, the faster it moves Recall that the V-nullcline is given by V = f (U) c. Moreover, du dz = V = f (U) c. Thus the flatter the wave, the faster it moves. Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
32 26. Oktober / 68 Stability of travelling wave Solutions Introduce perturbation around travelling wave solution u(z, t) = u c (z) + ωv(z, t) 0 < w 1 (36) The travelling wave solution u c (z) is stable if If lim v = 0. (37) t lim v = du c t dz. (38) then there are small translations along the x-axis.
33 26. Oktober / 68 Stability of travelling wave Solutions u(z, t) = u c (z) + ωv(z, t) 0 < ω 1 (39) Ansatz: v(z, t) = g(z) exp ( λt) (40) Determine eigenvalues to evaluate long-term behaviour of v(z, t).
34 26. Oktober / 68 Waves in 3D Ansatz: transform to spherical coordinate system x = r cos(ψ) sin(θ) y = r sin(ψ) sin(θ) z = r cos(θ) ψ [0, 2π) the azimuthal angle, and θ [0, π] the polar angle. Spherical coordinate system u(t, x, y, z) = u(t, r cos(ψ) sin(θ), r sin(ψ) sin(θ), r cos(θ))
35 26. Oktober / 68 The diffusion equation then reads where s is the spherical Laplace operator, s = 1 r 2 r (r 2 r ) + 1 r 2 sin(θ) du dt = f (u) + D su. (41) 2 θ (sin(θ) θ ) + 1 r 2 sin 2 (θ) ψ 2. (42)
36 26. Oktober / 68 Axisymmetry θ = 0 The spherical Laplace operator then simplifies to ψ = 0 (43) a = 2 r r r (44)
37 26. Oktober / 68 Transformation to travelling wave coordinates We write u(r, t) = u(r ct) = U(z), z = r ct, c 0. (45) u t u r but 1 r cannot be transformed. = c du dz = du dz (46)
38 26. Oktober / 68 Axisymmetric form of Fisher equation u t = u(1 u) + 1 r u r + u rr (47) does not become an ordinary differential equation in the variable z = r ct. does not possess travelling wavefront solutions in which a wave spreads out with constant speed, because of the u r /r term. wavespeed c(r) is a function of r: increases monotonically with r and reaches c(r) 2 for r large.
39 The calcium-stimulated-calcium release mechanism Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
40 Elevated levels of Ca 2+ in the cytoplasm stimulate further release of Ca 2+ from the endoplasmatic reticulum (ER) and sarcoplasmatic reticulum (SR) where available. Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
41 26. Oktober / 68 The kinetics for Ca 2+ can therefore be described by the following phenomenological model u = A(u) r(u) + L = f (u) (48) A(u) = k 1u 2 k 2 +u 2 : autocatalytic Ca 2+ accumulation r(u): linear resequestration process L: constant leakage flux
42 26. Oktober / 68 Graphical Analysis u = k 1u 2 k 2 + u 2 k 3u + L = f (u)
43 26. Oktober / 68 Importantly f (u) has three zeros, and in the following we will consider a simplified version, i.e. f (u) = A(u u 1 )(u 2 u)(u u 3 ), u 1 < u 2 < u 3, A > 0.
44 26. Oktober / 68 Since Ca 2+ can diffuse we have to consider a reaction-diffusion equation to adequately describe the Ca 2+ kinetics and we thus have As before we write u = f (u) + D d 2 u dx 2. (49) u(x, t) = u(x ct) = U(z), z = x ct, c 0. (50) and upon substitution into Eqn.(49) DU + cu + f (U) = 0 (51)
45 26. Oktober / 68 DU + cu + f (U) = 0 We need to have U = 0 at U = u 1 and U = u 3. Ansatz: U = α(u u 1 )(U u 3 ) and obtain upon substitution (U u 1 )(U u 3 )[Dα 2 (2U u 1 u 3 ) + cα A(U u 2 )] = 0. We thus require 2Da 2 A = 0 Da 2 (u 1 + u 3 ) cα Au 2 = 0 A α = 2D, c = A 2D (u 1 2u 2 + u 3 ).
46 26. Oktober / 68 As solution we thus obtain U(z) = u 3 + Ku 1 exp (a(u 3 u 1 )z) 1 + K exp (a(u 3 u 1 )z) (52) where K is an integration constant that determines the origin in the z-plane. The sign of c depends on the reaction kinetics, i.e. on u 1, u 2, u 3.
47 Calcium Waves on Amphibian Eggs Models by Cheer et al (1987) and Lane et al (1987) Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
48 A number of wave-like events can be observed upon fertilization. There are for example both chemical and mechanical waves which propagate on the surface of many vertebrate eggs. These waves arise from a combination of local reactions and long-range diffusion. Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
49 26. Oktober / 68 Waves are on the surface of a sphere, i.e. Ca 2+ wavefront is a ring that is propagating over the surface. We therefore consider diffusion in 3 dimensions and use the spherical coordinate system, i.e. u(t, x, y, z) = u(t, r cos(ψ) sin(θ), r sin(ψ) sin(θ), r cos(θ)) ψ [0, 2π) the azimuthal angle, and θ [0, π] the polar angle.
50 26. Oktober / 68 The diffusion equation then reads where s is the spherical Laplace operator, s = 1 r 2 r (r 2 r ) + 1 r 2 sin(θ) du dt = f (u) + D su. (53) θ (sin(θ) θ ) + 1 r 2 sin 2 (θ) ψ 2. (54) Given the inherent symmetries the spherical Laplace operator simplifies to 1 s = r 2 sin(θ) θ (sin(θ) θ ) = 1 ( 2 r 2 θ 2 + cot(θ) ). (55) θ 2
51 26. Oktober / 68 Cheer et al (1987) use the phenomenological description f (u) = A(u u 1 )(u u 2 )(u u 3 ) where A is a positive parameter to describe the excitable kinetics. At each fixed θ we obtain a wavefront solution of the form u(θ, t) = U(z), z = Rθ ct. (56)
52 26. Oktober / 68 We therefore have DU + (c + D R cot (θ))u + f (U) = 0 (57) This equation can be solved in the same way as Eqn.(51) as long as we set θ = const. We then obtain for the wave speed in analogy to Eqn.(52) AD c = 2 (u 1 2u 2 + u 3 ) D cot θ. (58) R
53 26. Oktober / 68 c = AD 2 (u 1 2u 2 + u 3 ) D cot θ. R The wave speed thus increases as the wave moves from the animal pole (θ = 0) to the vegetal pole (θ = π). In reality calcium waves slow down as they move towards the vegetal pole. There must therefore be important cortical properties that have been neglected by the model and which prevent the speeding up tendencies for propagating waves on the surface of a sphere.
54 Models for Somitogenesis Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
55 Somitogenesis Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
56 26. Oktober / 68 Clock and Wavefront Model Cooke, J., & Zeeman, E. C. (1976). A clock and wavefront model for control of the number of repeated structures during animal morphogenesis. Journal of theoretical biology, 58(2),
57 26. Oktober / 68 Somitogenesis Goldbeter, A., Gonze, D.,& Pourquie, O. (2007). Sharp developmental thresholds defined through bistability by antagonistic gradients of retinoic acid and FGF signaling. Developmental dynamics, 236(6), 1495Ð1508
58 26. Oktober / 68 Bistability model for Somitogenesis Goldbeter, A., Gonze, D.,& Pourquie, O. (2007). Sharp developmental thresholds defined through bistability by antagonistic gradients of retinoic acid and FGF signaling. Developmental dynamics, 236(6), 1495Ð1508
59 Wave Pinning Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
60 26. Oktober / 68 Wave Pinning Wave Pinning If a travelling wave eventually halts, the phenomenon is referred to as wave pinning.
61 Example: Cell Polarization Membrane translocation cycle of a typical Rho family GTPase. The inactive cytosolic form (B) diffuses much faster than the active membrane-bound form (A) and is approximately uniformly distributed. Mori et al. Biophys J 94, , (2008). Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
62 26. Oktober / 68 Conserved System a t = D 2 a a + f (a, b) x 2 b t = D 2 a b f (a, b) x 2 with Db D a and f (a, b) leading to bistability (at least two stable steady states (a, b ), (a +, b + )), e.g. f (a, b) = b ( k 0 + ) γa2 K 2 + a 2 δa For details, see Mori, Y., Jilkine, A. & Edelstein-Keshet, L. Wave-pinning and cell polarity from a bistable reaction-diffusion system. Biophys J 94, , doi: /biophysj (2008).
63 Single Variable Case (fixed b) a t = D 2 a a + f (a, b) x 2 with f (a, b) leading to bistability (at least two stable steady states (a, a + ). Suppose a develops a propagating front on 0 < x < L, sufficiently far from the boundaries. Then we can approximate a region of validity by < z <. Let A(z) = A(x ct) = a(x, t). A + ca + f (a, b) = 0 Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
64 Wave Speed A + ca + f (a, b) = 0 Multiply by A and integrate from to A (A + ca + f (a, b))dz = 0 Since A (± ) = 0 (flat part of wave) and A( ) = a and A( ) = a +, this integrates to give c [A ] 2 dz = f (a, b)a dz = a+ a f (a, b)da a+ a f (a, b)da c = ( ) 2 A z dz Computational Biology Group (CoBI), D-BSSE, ETHZ Prof Dagmar Iber, PhD DPhil Lecture 5 MSc Oktober / 68
65 26. Oktober / 68 Wave Speed Wave Speed c(b) = a+ a f (a, b)da ( ) 2 A z dz The denominator is always positive, while the numerator changes sign dependent on b.
66 26. Oktober / 68 Conditions for Wave Pinning At least two stable steady states for b min < b < b max a+ a f (a, b)da changes sign at b = b c c(b) = a+ a f (a, b)da ( ) 2 A z dz Mori, Y., Jilkine, A. & Edelstein-Keshet, L. Biophys J 94, , (2008). Pattern formation much faster than by diffusion (Turing Pattern)!!
67 26. Oktober / 68 Example for Wave Pinning Dashed line: initial conditions Mori, Y., Jilkine, A. & Edelstein-Keshet, L. Biophys J 94, , (2008).
68 26. Oktober / 68 Thanks!! Thanks for your attention! Slides for this talk will be available at:
Introduction to Travelling Waves Modeling Examples. Travelling Waves. Dagmar Iber. Computational Biology (CoBI), D-BSSE, ETHZ
Introduction to Travelling Waves Modeling Examples Travelling Waves Dagmar Iber Computational Biology (CoBI), D-BSSE, ETHZ 1 2 Introduction to Travelling Waves Modeling Examples Outline 1 Introduction
More informationFigure 1: Ca2+ wave in a Xenopus oocyte following fertilization. Time goes from top left to bottom right. From Fall et al., 2002.
1 Traveling Fronts Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 2 When mature Xenopus oocytes (frog
More information1 Existence of Travelling Wave Fronts for a Reaction-Diffusion Equation with Quadratic- Type Kinetics
1 Existence of Travelling Wave Fronts for a Reaction-Diffusion Equation with Quadratic- Type Kinetics Theorem. Consider the equation u t = Du xx + f(u) with f(0) = f(1) = 0, f(u) > 0 on 0 < u < 1, f (0)
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for COMPUTATIONAL BIOLOGY A. COURSE CODES: FFR 110, FIM740GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for COMPUTATIONAL BIOLOGY A COURSE CODES: FFR 110, FIM740GU, PhD Time: Place: Teachers: Allowed material: Not allowed: June 8, 2018, at 08 30 12 30 Johanneberg Kristian
More informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More informationDiffusion - The Heat Equation
Chapter 6 Diffusion - The Heat Equation 6.1 Goal Understand how to model a simple diffusion process and apply it to derive the heat equation in one dimension. We begin with the fundamental conservation
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationCell polarity models
Cell polarity models Universal features of polarizing cells 1. Ability to sense both steep and shallow external gradients (as small as 1% 2%) in vast range of concentrations. Polarization leads to an amplification
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More informationASYMPTOTIC AND BIFURCATION ANALYSIS OF WAVE-PINNING IN A REACTION-DIFFUSION MODEL FOR CELL POLARIZATION
ASYMPTOTIC AND BIFURCATION ANALYSIS OF WAVE-PINNING IN A REACTION-DIFFUSION MODEL FOR CELL POLARIZATION YOICHIRO MORI, ALEXANDRA JILKINE, AND LEAH EDELSTEIN-KESHET Abstract. We describe and analyze a bistable
More informationLecture 18: Bistable Fronts PHYS 221A, Spring 2017
Lecture 18: Bistable Fronts PHYS 221A, Spring 2017 Lectures: P. H. Diamond Notes: Xiang Fan June 15, 2017 1 Introduction In the previous lectures, we learned about Turing Patterns. Turing Instability is
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationSPATIO-TEMPORAL MODELLING IN BIOLOGY
SPATIO-TEMPORAL MODELLING IN BIOLOGY Prof Dagmar Iber, PhD DPhil ((Vorname Nachname)) 04/10/16 1 Challenge: Integration across scales Butcher et al (2004) Nat Biotech, 22, 1253-1259 INTERDISCIPLINARY WORK
More informationHopf bifurcations, and Some variations of diffusive logistic equation JUNPING SHIddd
Hopf bifurcations, and Some variations of diffusive logistic equation JUNPING SHIddd College of William and Mary Williamsburg, Virginia 23187 Mathematical Applications in Ecology and Evolution Workshop
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 27 May, 2004 1.30 to 3.30 PAPER 64 ASTROPHYSICAL FLUID DYNAMICS Attempt THREE questions. There are four questions in total. The questions carry equal weight. Candidates
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationMagnetic waves in a two-component model of galactic dynamo: metastability and stochastic generation
Center for Turbulence Research Annual Research Briefs 006 363 Magnetic waves in a two-component model of galactic dynamo: metastability and stochastic generation By S. Fedotov AND S. Abarzhi 1. Motivation
More informationDispersion relations, stability and linearization
Dispersion relations, stability and linearization 1 Dispersion relations Suppose that u(x, t) is a function with domain { < x 0}, and it satisfies a linear, constant coefficient partial differential
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More information(Refer Slide Time: 1:20) (Refer Slide Time: 1:24 min)
Engineering Chemistry - 1 Prof. K. Mangala Sunder Department of Chemistry Indian Institute of Technology, Madras Lecture - 5 Module 1: Atoms and Molecules Harmonic Oscillator (Continued) (Refer Slide Time:
More informationFundamentals of Transport Processes Prof. Kumaran Indian Institute of Science, Bangalore Chemical Engineering
Fundamentals of Transport Processes Prof. Kumaran Indian Institute of Science, Bangalore Chemical Engineering Module No # 05 Lecture No # 25 Mass and Energy Conservation Cartesian Co-ordinates Welcome
More informationModelling and Mathematical Methods in Process and Chemical Engineering
Modelling and Mathematical Methods in Process and Chemical Engineering Solution Series 3 1. Population dynamics: Gendercide The system admits two steady states The Jacobi matrix is ẋ = (1 p)xy k 1 x ẏ
More informationEmergence of polarity in early C. elegans Embryo. Abstract
Emergence of polarity in early C. elegans Embryo Pin-Yi Li Department of Physics, University of Illinois Urbana-Champaign (Dated: May 13, 2018) Abstract The establishment of cell polarity during embryogenesis
More informationENGI Duffing s Equation Page 4.65
ENGI 940 4. - Duffing s Equation Page 4.65 4. Duffing s Equation Among the simplest models of damped non-linear forced oscillations of a mechanical or electrical system with a cubic stiffness term is Duffing
More informationMath 46, Applied Math (Spring 2008): Final
Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the
More informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationF11AE1 1. C = ρν r r. r u z r
F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions
More informationODE, part 2. Dynamical systems, differential equations
ODE, part 2 Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Dynamical systems, differential equations Consider a system of n first order equations du dt = f(u, t),
More informationF1.9AB2 1. r 2 θ2 + sin 2 α. and. p θ = mr 2 θ. p2 θ. (d) In light of the information in part (c) above, we can express the Hamiltonian in the form
F1.9AB2 1 Question 1 (20 Marks) A cone of semi-angle α has its axis vertical and vertex downwards, as in Figure 1 (overleaf). A point mass m slides without friction on the inside of the cone under the
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
More informationMath 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More informationModels Involving Interactions between Predator and Prey Populations
Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predator-prey models are used to show the intricate
More informationInstabilities In A Reaction Diffusion Model: Spatially Homogeneous And Distributed Systems
Applied Mathematics E-Notes, 10(010), 136-146 c ISSN 1607-510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Instabilities In A Reaction Diffusion Model: Spatially Homogeneous And
More informationLIMIT CYCLE OSCILLATORS
MCB 137 EXCITABLE & OSCILLATORY SYSTEMS WINTER 2008 LIMIT CYCLE OSCILLATORS The Fitzhugh-Nagumo Equations The best example of an excitable phenomenon is the firing of a nerve: according to the Hodgkin
More informationDispersion relations, linearization and linearized dynamics in PDE models
Dispersion relations, linearization and linearized dynamics in PDE models 1 Dispersion relations Suppose that u(x, t) is a function with domain { < x 0}, and it satisfies a linear, constant coefficient
More informationIs chaos possible in 1d? - yes - no - I don t know. What is the long term behavior for the following system if x(0) = π/2?
Is chaos possible in 1d? - yes - no - I don t know What is the long term behavior for the following system if x(0) = π/2? In the insect outbreak problem, what kind of bifurcation occurs at fixed value
More informationMATH H53 : Final exam
MATH H53 : Final exam 11 May, 18 Name: You have 18 minutes to answer the questions. Use of calculators or any electronic items is not permitted. Answer the questions in the space provided. If you run out
More informationFront Propagation. Chapter Reaction-Diffusion Systems Single component systems
Chapter 8 Front Propagation 8.1 Reaction-Diffusion Systems We ve studied simple N = 1 dynamical systems of the form du dt = R(u). (8.1) Recall that the dynamics evolves u(t) monotonically toward the first
More informationBifurcation Analysis of Non-linear Differential Equations
Bifurcation Analysis of Non-linear Differential Equations Caitlin McCann 0064570 Supervisor: Dr. Vasiev September 01 - May 013 Contents 1 Introduction 3 Definitions 4 3 Ordinary Differential Equations
More informationIntroduction to PDEs and Numerical Methods: Exam 1
Prof Dr Thomas Sonar, Institute of Analysis Winter Semester 2003/4 17122003 Introduction to PDEs and Numerical Methods: Exam 1 To obtain full points explain your solutions thoroughly and self-consistently
More informationQ ( q(m, t 0 ) n) S t.
THE HEAT EQUATION The main equations that we will be dealing with are the heat equation, the wave equation, and the potential equation. We use simple physical principles to show how these equations are
More informationLecture 15: Biological Waves
Lecture 15: Biological Waves Jonathan A. Sherratt Contents 1 Wave Fronts I: Modelling Epidermal Wound Healing 2 1.1 Epidermal Wound Healing....................... 2 1.2 A Mathematical Model.........................
More informationAMATH 353 Lecture 9. Weston Barger. How to classify PDEs as linear/nonlinear, order, homogeneous or non-homogeneous.
AMATH 353 ecture 9 Weston Barger 1 Exam What you need to know: How to classify PDEs as linear/nonlinear, order, homogeneous or non-homogeneous. The definitions for traveling wave, standing wave, wave train
More informationExamples of Excitable Media. Excitable Media. Characteristics of Excitable Media. Behavior of Excitable Media. Part 2: Cellular Automata 9/7/04
Examples of Excitable Media Excitable Media Slime mold amoebas Cardiac tissue (& other muscle tissue) Cortical tissue Certain chemical systems (e.g., BZ reaction) Hodgepodge machine 9/7/04 1 9/7/04 2 Characteristics
More informationMath 337, Summer 2010 Assignment 5
Math 337, Summer Assignment 5 Dr. T Hillen, University of Alberta Exercise.. Consider Laplace s equation r r r u + u r r θ = in a semi-circular disk of radius a centered at the origin with boundary conditions
More informationRobot Control Basics CS 685
Robot Control Basics CS 685 Control basics Use some concepts from control theory to understand and learn how to control robots Control Theory general field studies control and understanding of behavior
More informationLecture 10: The quality of the Estimate
Computational Biology Group (CoBI), D-BSSE, ETHZ Lecture 10: The quality of the Estimate Prof Dagmar Iber, PhD DPhil BSc Biotechnology 2015/16 May 12, 2016 2 / 102 Contents 1 Review of Previous Lectures
More informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationMath 251 December 14, 2005 Answer Key to Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Name Section Math 51 December 14, 5 Answer Key to Final Exam There are 1 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning
More informationSolutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the
More informationPEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity
PEAT8002 - SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity Nick Rawlinson Research School of Earth Sciences Australian National University Anisotropy Introduction Most of the theoretical
More informationIntroduction to Physiology V - Coupling and Propagation
Introduction to Physiology V - Coupling and Propagation J. P. Keener Mathematics Department Coupling and Propagation p./33 Spatially Extended Excitable Media Neurons and axons Coupling and Propagation
More informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
More informationPattern formation and Turing instability
Pattern formation and Turing instability. Gurarie Topics: - Pattern formation through symmetry breaing and loss of stability - Activator-inhibitor systems with diffusion Turing proposed a mechanism for
More informationMath 2930 Worksheet Final Exam Review
Math 293 Worksheet Final Exam Review Week 14 November 3th, 217 Question 1. (* Solve the initial value problem y y = 2xe x, y( = 1 Question 2. (* Consider the differential equation: y = y y 3. (a Find the
More informationMATH 126 FINAL EXAM. Name:
MATH 126 FINAL EXAM Name: Exam policies: Closed book, closed notes, no external resources, individual work. Please write your name on the exam and on each page you detach. Unless stated otherwise, you
More informationarxiv: v1 [q-bio.cb] 12 Sep 2018
A COUPLED BULK-SURFACE MODEL FOR CELL POLARISATION Davide Cusseddu 1, Leah Edelstein-Keshet 2, John A. Mackenzie 3, Stéphanie Portet 4, Anotida Madzvamuse 1 arxiv:1809.04535v1 [q-bio.cb] 12 Sep 2018 Abstract:
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat
More informationLecture 37: Principal Axes, Translations, and Eulerian Angles
Lecture 37: Principal Axes, Translations, and Eulerian Angles When Can We Find Principal Axes? We can always write down the cubic equation that one must solve to determine the principal moments But if
More informationcarries the circle w 1 onto the circle z R and sends w = 0 to z = a. The function u(s(w)) is harmonic in the unit circle w 1 and we obtain
4. Poisson formula In fact we can write down a formula for the values of u in the interior using only the values on the boundary, in the case when E is a closed disk. First note that (3.5) determines the
More informationIntroduction to Biomathematics. Problem sheet
Introction to Biomathematics Problem sheet 1. A model for population growth is given in non-dimensional units in the form = u1 u ) u0) > 0. a) Sketch the graph of the function fu) = u1 u ) against u for
More informationPropagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R
Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract We consider semilinear
More informationIntroduction. Dagmar Iber Jörg Stelling. CSB Deterministic, SS 2015, 1.
Introduction Dagmar Iber Jörg Stelling joerg.stelling@bsse.ethz.ch CSB Deterministic, SS 2015, 1 Origins of Systems Biology On this assumption of the passage of blood, made as a basis for argument, and
More informationLow Emittance Machines
Advanced Accelerator Physics Course RHUL, Egham, UK September 2017 Low Emittance Machines Part 1: Beam Dynamics with Synchrotron Radiation Andy Wolski The Cockcroft Institute, and the University of Liverpool,
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationClassical Mechanics. Luis Anchordoqui
1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body
More informationLimit Cycles II. Prof. Ned Wingreen MOL 410/510. How to prove a closed orbit exists?
Limit Cycles II Prof. Ned Wingreen MOL 410/510 How to prove a closed orbit eists? numerically Poincaré-Bendison Theorem: If 1. R is a closed, bounded subset of the plane. = f( ) is a continuously differentiable
More informationFinite difference method for solving Advection-Diffusion Problem in 1D
Finite difference method for solving Advection-Diffusion Problem in 1D Author : Osei K. Tweneboah MATH 5370: Final Project Outline 1 Advection-Diffusion Problem Stationary Advection-Diffusion Problem in
More informationMATH3203 Lecture 1 Mathematical Modelling and ODEs
MATH3203 Lecture 1 Mathematical Modelling and ODEs Dion Weatherley Earth Systems Science Computational Centre, University of Queensland February 27, 2006 Abstract Contents 1 Mathematical Modelling 2 1.1
More information0 = p. 2 x + 2 w. z +ν w
Solution (Elliptical pipe flow (a Using the Navier Stokes equations in three dimensional cartesian coordinates, given that u =, v = and w = w(x,y only, and assuming no body force, we are left with = p
More informationDiffusion, Reaction, and Biological pattern formation
Diffusion, Reaction, and Biological pattern formation Morphogenesis and positional information How do cells know what to do? Fundamental questions How do proteins in a cell segregate to front or back?
More informationMathematical Models of Biological Systems
Mathematical Models of Biological Systems CH924 November 25, 2012 Course outline Syllabus: Week 6: Stability of first-order autonomous differential equation, genetic switch/clock, introduction to enzyme
More informationLecture Notes for PHY 405 Classical Mechanics
Lecture Notes for PHY 405 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics September 1, 2005 Chapter
More informationx + ye z2 + ze y2, y + xe z2 + ze x2, z and where T is the
1.(8pts) Find F ds where F = x + ye z + ze y, y + xe z + ze x, z and where T is the T surface in the pictures. (The two pictures are two views of the same surface.) The boundary of T is the unit circle
More information18.325: Vortex Dynamics
8.35: Vortex Dynamics Problem Sheet. Fluid is barotropic which means p = p(. The Euler equation, in presence of a conservative body force, is Du Dt = p χ. This can be written, on use of a vector identity,
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationDynamical Systems and Chaos Part II: Biology Applications. Lecture 6: Population dynamics. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part II: Biology Applications Lecture 6: Population dynamics Ilya Potapov Mathematics Department, TUT Room TD325 Living things are dynamical systems Dynamical systems theory
More informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More informationVectors, matrices, eigenvalues and eigenvectors
Vectors, matrices, eigenvalues and eigenvectors 1 ( ) ( ) ( ) Scaling a vector: 0.5V 2 0.5 2 1 = 0.5 = = 1 0.5 1 0.5 ( ) ( ) ( ) ( ) Adding two vectors: V + W 2 1 2 + 1 3 = + = = 1 3 1 + 3 4 ( ) ( ) a
More information20 The Hydrogen Atom. Ze2 r R (20.1) H( r, R) = h2 2m 2 r h2 2M 2 R
20 The Hydrogen Atom 1. We want to solve the time independent Schrödinger Equation for the hydrogen atom. 2. There are two particles in the system, an electron and a nucleus, and so we can write the Hamiltonian
More informationMoving Finite Elements; From Shallow Water equations to Aggregation of microglia in Alzheimer s disease
Moving Finite Elements; From Shallow Water equations to Aggregation of microglia in Alzheimer s disease Abigail Wacher Past collaborators on this topic: Ian Sobey (Oxford Computing Laboratory) Keith Miller
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationGradient, Divergence and Curl in Curvilinear Coordinates
Gradient, Divergence and Curl in Curvilinear Coordinates Although cartesian orthogonal coordinates are very intuitive and easy to use, it is often found more convenient to work with other coordinate systems.
More informationWave Pinning, Actin Waves, and LPA
Wave Pinning, Actin Waves, and LPA MCB 2012 William R. Holmes Intercellular Waves Weiner et al., 2007, PLoS Biology Dynamic Hem1 waves in neutrophils Questions? How can such waves / pulses form? What molecular
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More information+ i. cos(t) + 2 sin(t) + c 2.
MATH HOMEWORK #7 PART A SOLUTIONS Problem 7.6.. Consider the system x = 5 x. a Express the general solution of the given system of equations in terms of realvalued functions. b Draw a direction field,
More informationBreakdown of Pattern Formation in Activator-Inhibitor Systems and Unfolding of a Singular Equilibrium
Breakdown of Pattern Formation in Activator-Inhibitor Systems and Unfolding of a Singular Equilibrium Izumi Takagi (Mathematical Institute, Tohoku University) joint work with Kanako Suzuki (Institute for
More informationLecture 2c: Satellite Orbits
Lecture 2c: Satellite Orbits Outline 1. Newton s Laws of Mo3on 2. Newton s Law of Universal Gravita3on 3. Kepler s Laws 4. Pu>ng Newton and Kepler s Laws together and applying them to the Earth-satellite
More information7a3 2. (c) πa 3 (d) πa 3 (e) πa3
1.(6pts) Find the integral x, y, z d S where H is the part of the upper hemisphere of H x 2 + y 2 + z 2 = a 2 above the plane z = a and the normal points up. ( 2 π ) Useful Facts: cos = 1 and ds = ±a sin
More informationMATH20411 PDEs and Vector Calculus B
MATH2411 PDEs and Vector Calculus B Dr Stefan Güttel Acknowledgement The lecture notes and other course materials are based on notes provided by Dr Catherine Powell. SECTION 1: Introctory Material MATH2411
More information