Implicit-explicit variational integration of highly oscillatory problems
|
|
- Gwendolyn Jacobs
- 5 years ago
- Views:
Transcription
1 Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na]. Department of Matematics, UCSD, astern@mat.ucsd.edu
2 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Te Problem: Multiple Time Scales vs. Numerical Integrators Many systems in Lagrangian mecanics ave components acting on different time scales, for example:. Elasticity: Several spatial elements of varying stiffness, resulting from irregular meses and/or inomogeneous materials.. Planetary Dynamics: N-body problem wit nonlinear gravitational forces, arising from pairwise inverse-square potentials. Multiple time scales result from te different distances between te bodies, mass ratios, etc. 3. Higly Oscillatory Problems: Potential energy can be split into a fast linear oscillatory component and a slow nonlinear component. Tese problems are widely encountered in modeling molecular dynamics, but ave also been used to model oter diverse applications, for example, in computer animation. Existing metods tend to waste too muc computational effort evaluating te non-stiff forces e.g., Störmer/Verlet, implicit midpoint or suffer from resonance instability problems e.g., Verlet-I/r-RESPA, AVI, many trigonometric integrators. Can we develop geometric numerical integrators tat are stable wit respect to te stiff forces, but still efficient wit respect to te non-stiff forces?
3 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Review: Te Discrete Lagrangian and Action Sum Suppose we ave a mecanical system on a configuration manifold Q, specified by a Lagrangian L: T Q R. Given a set of discrete time points t < < t N wit uniform step size, we wis to compute a numerical approximation q n q t n, n =,..., N, to te continuous trajectory qt. To construct a variational integrator for tis problem, we define a discrete Lagrangian L : Q Q R, replacing tangent vectors by pairs of consecutive configuration points, so tat in some sense we ave te approximation tn+ L qn, q n+ L q, q dt. Ten te action integral over [t, t N ] is approximated by te discrete action sum tn S [q] = N n= L qn, q n+ tn t L q, q dt.
4 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Review: Te Discrete Euler Lagrange Equations and Legendre Transform If we apply Hamilton s principle to tis action sum, so tat δs [q] = for fixedendpoint variations, ten tis yields te discrete Euler Lagrange equations D L qn, q n+ + D L qn, q n =, n =,..., N, were D and D denote partial differentiation in te first and second arguments, respectively. Tis defines a two-step numerical metod on Q Q, mapping q n, q n qn, q n+. Te equivalent one-step metod on te cotangent bundle T Q, mapping q n, p n qn+, p n+, is defined by te discrete Legendre transform p n = D L qn, q n+, pn+ = D L qn, q n+, were te first equation updates q, and te second updates p. Here, L is a generating function for te symplectic map q n, p n qn+, p n+. 3
5 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Review: Te Trapezoidal and Midpoint Discrete Lagrangians Consider a Lagrangian of te form L q, q = qt M q V q, were Q = R d, M is a constant d d mass matrix, and V : Q R is a potential. If we use linear interpolation of q wit trapezoidal quadrature to approximate te contribution of V to te action integral, we get L trap qn, q n+ = qn+ q T n qn+ q n M V q n + V qn+, wic we call te trapezoidal discrete Lagrangian. It is straigtforward to see tat te discrete Euler Lagrange equations for L trap correspond to te explicit Störmer/Verlet metod. Alternatively, if we use midpoint quadrature to approximate te integral of te potential, tis yields te midpoint discrete Lagrangian, L mid qn, q n+ = qn+ q T n qn+ q n qn + q n+ M V, for wic te resulting integrator is te implicit midpoint metod. 4
6 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Te IMEX Discrete Lagrangian Suppose tat we ave a Lagrangian of te form L q, q = qt M q Uq W q, were U is a slow potential and W is a fast potential, for te configuration space Q = R d. Define te IMEX discrete Lagrangian L IMEX qn+ q T n qn+ q n qn, q n+ = M U q n + U qn+ qn + q n+ W, using explicit trapezoidal approximation for te slow potential and implicit midpoint approximation for te fast potential. Te discrete Euler Lagrange equations give a two-step variational integrator on Q Q, [ q n+ q n + q n = M U q n + W qn + q n + ] W qn + q n+. 5
7 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Te IMEX Algoritm on T Q Te discrete Legendre transform for te IMEX discrete Lagrangian is qn+ q n p n = M qn+ q n p n+ = M If we introduce te intermediate stages, + U q n + W U q n+ W qn + q n+, qn + q n+. p + n = p n U q n, p n+ = p n+ + U q n+, ten tis gives us te following impulse-type integration algoritm: Step : p + n = p n U q n, explicit kick Step : q n+ pn+ p = q n + M + n + pn+, implicit midpoint = p n + qn + q n+ W, Step 3: p n+ = p n+ U q n+. explicit kick 6
8 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Application to Higly Oscillatory Problems For igly oscillatory problems on Q = R d, we start by taking a quadratic fast potential W q = qt Ω q, Ω R d d symmetric and positive semidefinite. A prototypical Ω is given by te block-diagonal matrix Ω =, ω, ωi were some of te degrees of freedom are subjected to an oscillatory force wit constant fast frequency ω. We also denote te slow force gq = Uq and assume, witout loss of generality, tat te constant mass matrix is given by M = I. Terefore, te nonlinear system we wis to approximate numerically is q + Ω q = gq. 7
9 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 IMEX as Störmer/Verlet wit a Modified Mass Matrix Applying te IMEX metod to tis example, we get te discrete Lagrangian L IMEX qn, q n+ = qn+ q T n qn+ q n U q n + U qn+ qn + q T n+ Ω qn + q n+, and so te two-step IMEX sceme is given by te discrete Euler Lagrange equations q n+ q n + q n + 4 Ω q n+ + q n + q n = g q n. Rearranging terms, we can rewrite tis as [ ] I + qn+ 4 Ω q n + q n + Ω q n = g q n, wic is equivalent to Störmer/Verlet wit a modified mass matrix I + Ω/. 8
10 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Te Discrete Lagrangian and Modified Mass Matrix Proposition. Suppose we ave a Lagrangian L q, q = qt M q qt Ω q and its corresponding midpoint discrete Lagrangian L mid. Next, define te modified Lagrangian L q, q = qt M q qt Ω q, aving te same quadratic potential but a different mass matrix M, and take its trapezoidal discrete Lagrangian L trap. Ten Lmid L trap wen M = M + Ω/. Proof. Te midpoint discrete Lagrangian is given by L mid qn, q n+ = qn+ q T n qn+ q n M qn + q T n+ Ω qn + q n+. Now, notice tat we can rearrange te terms qn + q T n+ Ω qn + q n+ qn+ q T n = Ω qn+ q n qt nω q n qt n+ Ω q n+ qn+ q T n Ω qn+ q n = qt nω q n qt n+ Ω q n+. 9
11 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Terefore te discrete Lagrangian can be written in te trapezoidal form [ L mid qn+ q T ] qn+ n Ω q n qn, q n+ = M + qt nω q n + qt n+ Ω q n+, wic is precisely L trap qn, q n+ wen M = M + Ω/. Corollary. Consider a igly oscillatory system wit an arbitrary slow potential U, quadratic fast potential W q = qt Ω q, and constant mass matrix M = I, so tat te Lagrangian L and IMEX discrete Lagrangian L IMEX are defined as above. Next, take te modified Lagrangian L wit te same potentials but different mass matrix M. Ten L IMEX L trap wen M = I + Ω/.
12 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Analysis of Linear Resonance Stability For a armonic oscillator wit unit mass and frequency ν, te Störmer/Verlet metod is linearly stable if and only if ν ; for a system wit constant mass m and spring constant ν, tis condition generalizes to ν 4m. Linear model problem: Let Uq = qt q, and W q = qt Ω q, were Ω = ωi for some ω, and again let M = I. Teorem. Te IMEX metod is linearly stable, for te system described above, if and only if i.e., if and only if is a stable time step size for te slow oscillator alone. Proof. As proved in te previous section, te IMEX metod for tis system is equivalent to Störmer/Verlet wit te modified mass matrix I + Ω/. Now, tis modified oscillatory system as constant mass m = +ω/ and spring constant ν = +ω. Terefore, te necessary and sufficient condition for linear stability is + ω ω, and since te ω terms cancel on bot sides, tis is equivalent to. 4, or Tis sows tat, in contrast to multiple-time-stepping metods, te IMEX metod does not exibit linear resonance instability.
13 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Numerical Experiment: Coupled Linear Oscillators r!respa Variational IMEX Maximum absolute energy error!!!3! ! /" Maximum energy error of r-respa and variational IMEX, integrated over te time interval [, ] for a range of parameters ω. Te r-respa metod exibits resonance instability near integer values of ω/π, wile te variational IMEX metod remains stable.
14 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Numerical Experiment: Te Fermi Pasta Ulam Problem q q q m q m stiff armonic soft nonlinear Source: Hairer, Lubic, and Wanner 6. Cain wit alternating soft nonlinear stiff linear springs I + I + I I I I Slow energy excange among te stiff springs. 5 3
15 Stern, IMEX variational integration of igly oscillatory problems I + I + I I I + II + I 3 I 3 I I I + I + I3 I3 I I.6 I3 I.8.4 I I + I + I3.6 I I + I + I3 I3 I I3 I Reference solution: Störmer/Verlet wit time step size = IMEX wit = IMEX wit = IMEX wit = Störmer/Verlet wit = IMEX wit =.3..5 Störmer/Verlet wit = I3 I I + I + I Structured Integrators Worksop, May 7, IMEX wit = IMEX wit =.3. 4
16 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, IMEX wit =. 3 4 IMEX wit =.3 Numerical simulation of te FPU problem for T = 4, wic sows te beavior of te IMEX metod on te ω scale. For =., we already ave ω = 5, yet te oscillatory beavior and adiabatic invariant are qualitatively correct. By contrast, for =.3, te metod as begun to blow up; oscillatory coupling is a drawback of implicit midpoint metods for large time steps. 5
17 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Analysis of Slow Energy Excange First, let us rewrite te fast oscillatory system q + Ω q = as te first-order system Ω q Ω Ωq =, ṗ Ω p so it follows tat te exact solution satisfies Ωqt + cos Ω sin Ω Ωqt =. pt + sin Ω cos Ω pt Now, in tese coordinates, te implicit midpoint metod as te expression I Ω/ Ωqn+ I Ω/ Ωqn =. Ω/ I Ω/ I p n+ Terefore, if we take te skew matrix A = Ωqn+ p n+ p n Ω, it follows tat Ω = I A/ Ωqn I + A/. Notice tat te expression I A/ I + A/ = cay A is te Cayley transform. p n 6
18 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Implicit Midpoint and Modified Frequency Because cay and exp are bot maps from sod SOd, te midpoint metod corresponds to a modified rotation matrix, were Ω is replaced by Ω suc tat Ω/ = tan Ω/ Tis gives anoter interpretation of te resonance stability wen Ω = ωi and Ω = ωi. We always ave ω < π, since te Cayley transform maps to a rotation by π only in te limit as ω. Terefore, we never encounter resonance points for finite ω. We now write te variational IMEX metod as te following modified impulse sceme: Step : p + n = p n U q n, explicit kick Step : Ωqn+ p n+ cos Ω sin Ω Ωqn =, implicit midpoint sin Ω cos Ω Step 3: p n+ = p n+ U q n+. explicit kick In oter words, te implicit midpoint step doesn t just smear out te oscillations in fact, it resolves te oscillations of some modified problem. Because te propagation matrix is still special ortogonal, tis does not affect te fast energy component. p + n 7
19 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Consistency of Slow Energy Excange Hairer et al. ave analyzed te slow energy excange beavior of impulse metods using modulated Fourier expansion. By casting IMEX as a modified impulse metod, we were able to use tese same tecniques to prove te following teorem: Teorem. Let te variational IMEX metod be applied to te igly oscillatory problem above, and suppose te numerical solution remains bounded. Ten te ordinary differential equation [...] describing te slow energy excange in te numerical solution, is consistent wit tat for [...] te exact solution; tis olds up to order O ω 3. In fact, tis is not true for eiter Störmer/Verlet or implicit midpoint so IMEX is not merely ceaper, but also better for tese applications. Te only trigonometric/exponential metod saring tis property is Deuflard/impulse, wic also suffers from resonance instability problems. Te most comparable integrators in tis respect appear to be multi-force trigonometric metods, but tese require at least twice as many force evaluations as IMEX. 8
20 Stern, IMEX variational integration of igly oscillatory problems Structured Integrators Worksop, May 7, 9 Summary Te variational IMEX metod is developed by splitting te discrete Lagrangian into slow explicit and fast implicit components, and applying separate quadrature rules trapezoidal and midpoint, respectively. For igly oscillatory problems, tis is equivalent to Störmer/Verlet wit a modified mass matrix. Tis leads to unconditional linear stability in te fast modes, and in particular, te absence of linear resonance instability. Te Fermi Pasta Ulam example demonstrates tat te variational IMEX metod does not attain its stability merely by smooting out te fast frequencies, in a way tat migt destroy te structure of any fast-slow nonlinear coupling. Rater, despite te fact tat it does not resolve te fast frequencies, te metod is still capable of capturing te complex multiscale interactions seen in te FPU problem. 9
arxiv: v3 [math.na] 22 Apr 2009
IMPLICIT-EXPLICIT VARIATIONAL INTEGRATION OF HIGHLY OSCILLATORY PROBLEMS ARI STERN AND EITAN GRINSPUN arxiv:88.39v3 [mat.na] Apr 9 ABSTRACT. In tis paper we derive a variational integrator for certain
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationLecture 2: Symplectic integrators
Geometric Numerical Integration TU Müncen Ernst Hairer January February 010 Lecture : Symplectic integrators Table of contents 1 Basic symplectic integration scemes 1 Symplectic Runge Kutta metods 4 3
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationMODIFIED DIFFERENTIAL EQUATIONS. Dedicated to Prof. Michel Crouzeix
ESAIM: PROCEEDINGS, September 2007, Vol.21, 16-20 Gabriel Caloz & Monique Dauge, Editors MODIFIED DIFFERENTIAL EQUATIONS Pilippe Cartier 1, Ernst Hairer 2 and Gilles Vilmart 1,2 Dedicated to Prof. Micel
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationAMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.
Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error
More informationRecent Progress in the Integration of Poisson Systems via the Mid Point Rule and Runge Kutta Algorithm
Recent Progress in te Integration of Poisson Systems via te Mid Point Rule and Runge Kutta Algoritm Klaus Bucner, Mircea Craioveanu and Mircea Puta Abstract Some recent progress in te integration of Poisson
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationApplication of numerical integration methods to continuously variable transmission dynamics models a
ttps://doi.org/10.1051/ssconf/2018440005 Application of numerical integration metods to continuously variable transmission dynamics models a Orlov Stepan 1 1 Peter te Great St. Petersburg Polytecnic University
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Matematics National Institute of Tecnology Durgapur Durgapur-7109 email: anita.buie@gmail.com 1 . Capter 8 Numerical Solution of Ordinary
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationMAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016
MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationUniversity Mathematics 2
University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Matematics and Computer Science. ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI3097 TU) Tuesday January 9 008, 9:00-:00
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationMath 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0
3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,
More informationBlanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS
Opuscula Matematica Vol. 26 No. 3 26 Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationNotes on Multigrid Methods
Notes on Multigrid Metods Qingai Zang April, 17 Motivation of multigrids. Te convergence rates of classical iterative metod depend on te grid spacing, or problem size. In contrast, convergence rates of
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationQuasiperiodic phenomena in the Van der Pol - Mathieu equation
Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van
More informationFractional Derivatives as Binomial Limits
Fractional Derivatives as Binomial Limits Researc Question: Can te limit form of te iger-order derivative be extended to fractional orders? (atematics) Word Count: 669 words Contents - IRODUCIO... Error!
More informationRecall from our discussion of continuity in lecture a function is continuous at a point x = a if and only if
Computational Aspects of its. Keeping te simple simple. Recall by elementary functions we mean :Polynomials (including linear and quadratic equations) Eponentials Logaritms Trig Functions Rational Functions
More informationChapter Seven The Quantum Mechanical Simple Harmonic Oscillator
Capter Seven Te Quantum Mecanical Simple Harmonic Oscillator Introduction Te potential energy function for a classical, simple armonic oscillator is given by ZÐBÑ œ 5B were 5 is te spring constant. Suc
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationOn convexity of polynomial paths and generalized majorizations
On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More informationMATH1131/1141 Calculus Test S1 v8a
MATH/ Calculus Test 8 S v8a October, 7 Tese solutions were written by Joann Blanco, typed by Brendan Trin and edited by Mattew Yan and Henderson Ko Please be etical wit tis resource It is for te use of
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare ttp://ocw.mit.edu 5.74 Introductory Quantum Mecanics II Spring 9 For information about citing tese materials or our Terms of Use, visit: ttp://ocw.mit.edu/terms. Andrei Tokmakoff, MIT
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationIntegral Calculus, dealing with areas and volumes, and approximate areas under and between curves.
Calculus can be divided into two ke areas: Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and minima problems Integral
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationMath Module Preliminary Test Solutions
SSEA Summer 207 Mat Module Preliminar Test Solutions. [3 points] Find all values of tat satisf =. Solution: = ( ) = ( ) = ( ) =. Tis means ( ) is positive. Tat is, 0, wic implies. 2. [6 points] Find all
More informationarxiv: v1 [math.na] 8 Mar 2011
Variational integrators for electric circuits Sina Ober-Blöbaum a, Molei Tao b, Mulin Ceng c, Houman Owadi b,c, Jerrold E. Marsden b,c, a Computational Dynamics and Optimal Control, University of Paderborn,
More informationSolution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.
December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More informationMath Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim
Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationCHAPTER 3: Derivatives
CHAPTER 3: Derivatives 3.1: Derivatives, Tangent Lines, and Rates of Cange 3.2: Derivative Functions and Differentiability 3.3: Tecniques of Differentiation 3.4: Derivatives of Trigonometric Functions
More informationGradient Descent etc.
1 Gradient Descent etc EE 13: Networked estimation and control Prof Kan) I DERIVATIVE Consider f : R R x fx) Te derivative is defined as d fx) = lim dx fx + ) fx) Te cain rule states tat if d d f gx) )
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationMath 242: Principles of Analysis Fall 2016 Homework 7 Part B Solutions
Mat 22: Principles of Analysis Fall 206 Homework 7 Part B Solutions. Sow tat f(x) = x 2 is not uniformly continuous on R. Solution. Te equation is equivalent to f(x) = 0 were f(x) = x 2 sin(x) 3. Since
More informationName: Answer Key No calculators. Show your work! 1. (21 points) All answers should either be,, a (finite) real number, or DNE ( does not exist ).
Mat - Final Exam August 3 rd, Name: Answer Key No calculators. Sow your work!. points) All answers sould eiter be,, a finite) real number, or DNE does not exist ). a) Use te grap of te function to evaluate
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationTangent Lines-1. Tangent Lines
Tangent Lines- Tangent Lines In geometry, te tangent line to a circle wit centre O at a point A on te circle is defined to be te perpendicular line at A to te line OA. Te tangent lines ave te special property
More informationFINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS
J. KSIAM Vol.22, No.2, 101 113, 2018 ttp://dx.doi.org/10.12941/jksiam.2018.22.101 FINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS DEOK-KYU JANG AND JAE-HONG PYO DEPARTMENT
More informationA Reconsideration of Matter Waves
A Reconsideration of Matter Waves by Roger Ellman Abstract Matter waves were discovered in te early 20t century from teir wavelengt, predicted by DeBroglie, Planck's constant divided by te particle's momentum,
More information1 (10) 2 (10) 3 (10) 4 (10) 5 (10) 6 (10) Total (60)
First Name: OSU Number: Last Name: Signature: OKLAHOMA STATE UNIVERSITY Department of Matematics MATH 2144 (Calculus I) Instructor: Dr. Matias Sculze MIDTERM 1 September 17, 2008 Duration: 50 minutes No
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationChapter 1D - Rational Expressions
- Capter 1D Capter 1D - Rational Expressions Definition of a Rational Expression A rational expression is te quotient of two polynomials. (Recall: A function px is a polynomial in x of degree n, if tere
More information1. For d=3,2 from ε<< ε F it follows that ετ >> e-e h, i.e.,
Quasiparticle decay rate at T = 0 in a clean Fermi Liquid. ω +ω Fermi Sea τ e e ( ) F ( ) log( ) Conclusions:. For d=3, from > e-e, i.e., tat te qusiparticles are well determined
More informationExam 1 Solutions. x(x 2) (x + 1)(x 2) = x
Eam Solutions Question (0%) Consider f() = 2 2 2 2. (a) By calculating relevant its, determine te equations of all vertical asymptotes of te grap of f(). If tere are none, say so. f() = ( 2) ( + )( 2)
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationMaterial for Difference Quotient
Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient
More information1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).
. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use
More informationSome Review Problems for First Midterm Mathematics 1300, Calculus 1
Some Review Problems for First Midterm Matematics 00, Calculus. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd,
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationA First-Order System Approach for Diffusion Equation. I. Second-Order Residual-Distribution Schemes
A First-Order System Approac for Diffusion Equation. I. Second-Order Residual-Distribution Scemes Hiroaki Nisikawa W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department of Aerospace
More informationMATH 1A Midterm Practice September 29, 2014
MATH A Midterm Practice September 9, 04 Name: Problem. (True/False) If a function f : R R is injective, ten f as an inverse. Solution: True. If f is injective, ten it as an inverse since tere does not
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationLecture 10: Carnot theorem
ecture 0: Carnot teorem Feb 7, 005 Equivalence of Kelvin and Clausius formulations ast time we learned tat te Second aw can be formulated in two ways. e Kelvin formulation: No process is possible wose
More informationAnalytic Functions. Differentiable Functions of a Complex Variable
Analytic Functions Differentiable Functions of a Complex Variable In tis capter, we sall generalize te ideas for polynomials power series of a complex variable we developed in te previous capter to general
More informationarxiv: v1 [math.na] 7 Mar 2019
Local Fourier analysis for mixed finite-element metods for te Stokes equations Yunui He a,, Scott P. MacLaclan a a Department of Matematics and Statistics, Memorial University of Newfoundland, St. Jon
More informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
More informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
More informationA CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume XLVIII, Number 3, September 2003 A CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION GH. MICULA, E. SANTI, AND M. G. CIMORONI Dedicated to
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationCONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION
CONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION SÖREN BARTELS AND ANDREAS PROHL Abstract. Te Landau-Lifsitz-Gilbert equation describes dynamics of ferromagnetism,
More informationQuaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 01 Aug 08, 2016.
Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals Gary D. Simpson gsim1887@aol.com rev 1 Aug 8, 216 Summary Definitions are presented for "quaternion functions" of a quaternion. Polynomial
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
More information