arxiv: v1 [math-ph] 9 Jan 2012
|
|
- Emma Wiggins
- 6 years ago
- Views:
Transcription
1 Splitting methods for Levitron Problems Jürgen Geiser Karl Felix Lüskow arxiv: v1 [math-ph] 9 Jan 01 December 1, 01 Abstract In this paper we describe splitting methods for solving Levitron, which is motivated to simulate magnetostatic traps of neutral atoms or ion traps. The idea is to levitate a magnetic spinning top in the air repelled by a base magnet. The main problem is the stability of the reduced Hamiltonian, while it is not defined at the relative equilibrium. Here it is important to derive stable numerical schemes with high accuracy. For the numerical studies, we propose novel splitting schemes and analyze their behavior. We deal with a Verlet integrator and improve its accuracy with iterative and extrapolation ideas. Such a Hamiltonian splitting method, can be seen as geometric integrator and saves computational time while decoupling the full equation system. Experiments based on the Levitron model are discussed. Keywords splitting method, Verlet integrator, iterative and extrapolation methods, Levitron problem. AMS subject classifications. 65M1, 65L06, 65P10. 1 Introduction We are motivated to simulate a Levitron, which is a magnetic spinning top and can levitate in a magnetic field. The main problem of such a nonlinear problem is to achieve a stability for the calculation of the critical splint rate. While the stability of Levitrons are discussed in the work of [] and their dynamics in [?], we concentrate on improving the standard time-integrator schemes for the reduced Hamiltonian systems. It is important to derive stable numerical schemes with high accuracy to compute the non-dissipative equation of motions. For the numerical studies, we propose novel splitting schemes and analyze their behavior. We deal with a standard Verlet integrator and improve its accuracy with University of Greifswald, Institute of Physics, Felix-Hausdorff-Str. 6, D Greifswald, Germany, juergen.geiser@uni-greifswald.de University of Greifswald, Institute of Physics, Felix-Hausdorff-Str. 6, D Greifswald, Germany, tt-karl@freenet.de 1
2 1 INTRODUCTION iterative and extrapolation ideas. Such a Hamiltonian splitting method, can be seen as geometric integrator and saves computational time while decoupling the full equation system, see the splitting ideas in the overview article [1]. In the following we describe the reduced model of Gans [6] and an extension based on a novel idea of magentic field of Dullin [] for a disk. 1.1 Hamiltonian of Gans In the paper, we deal with the following problem reduced Hamiltonian: H = 1 p 1 + p + p + p 4 a + p 5 P 6 cos q 4 a sin + p 6 q 4 c ] Ψ Ψ Ψ M [sin q 4 cos q 5 sin q 5 + cos q 4 + q 1 1 The evolution of the dynamical variable uq, p including q and p themselves is given by the Poisson bracket, t uq, p = u H p u p H = A + Buq, p. For the non-separable Hamiltonian of 1, we have: q = H p p, q = The same is given for: p 1, p, p, p4 a, p5 p6 cos q4 a sin q 4, p6cos q 4+a/c sin 4 p5 cos q4 a sin q 4 ṗ = H p, q Ψ Ψ = M sin q 4 cos q cos q 4 1 M sin Ψ Ψ q 4 sin q 5 + cos q 5 1 Ψ Ψ Ψ M cos q 4 sin q 5 + cos q 5 sin q 4 q 1 M sin q 4 cos q 5 Ψ sin q 5 Ψ 1 A and B are Lie operators, or vector fields Ψ, M sin q 4 cos q 5 + cos q 4 Ψ 1, Ψ + cos q 4, p 6p 5 p 6 cos q 4 a sin q 4 cos q 4p 5 p 6 cos q 4 a sin q 4,, 0 4 A = H p B = H p 5 The transfer to the operators are given in the following description.
3 1 INTRODUCTION The exponential operators e ha and e hb are then just shift operators, with T h is a symmetric second order splitting method: T,V V h t = e t/b e ta e t/b. 6 and corresponds to the velocity form of the Verlet algorithm VV. Further the splitting scheme: T,P V h t = e t/a e tb e t/a. 7 and corresponds to the position-form of the Verlet algorithm PV. See also the derivation of the Verlet algorithm in Appendix 7. T,V V h t = S AB h, the symplectic Verlet or leap-frog algorithm is given as: We start with q 0, p 0 t = qt n, pt n t : q 1, p 1 t = e h/b q 0, p 0 t = I 1 h i H p i, q i q 0, p 0 t, 8 p i q, v t = e ha q 1, v 1 t = I + h i H p p i, q i q 1, p 1 t, 9 i q, v t = e h/b q, p t = I 1 h i H p i, q i q, p t. 10 p i And the substitution is given the algorithm for one time-step n n + 1 and we obtain: qt n+1, vt n+1 t = q, v t. 1. Higher order Expansion of Verlet-algorithm In the following we extend the Verlet algorithm with respect to higher order terms. Such terms are important in the application with iterative schemes to achieve higher order schemes. We start with q 0, p 0 t = qt n, pt n t : q 1, p 1 t = e h/b q 0, p 0 t = I + N j=1 1 1 j! h i H p i, q i p i jq0, p 0 t, 11
4 ITERATIVE SCHEMES FOR COUPLED PROBLEMS 4 q, v t = e ha q 1, v 1 t = I + N j=1 1 H h j! p p i, q i jq0, p 0 t, 1 i i q, v t = e h/b q, p t = I + N j=1 1 1 j! h i H p i, q i jq0, p 0 t. 1 p i And the substitution is given the algorithm for one time-step n n + 1 and we obtain: qt n+1, vt n+1 t = q, v t. Iterative Schemes for coupled problems Based on the nonlinear equations, we have to deal with linearization or nonlinear averaging techniques. In the following, we discuss the fixed point iteration and Newton s method. We solve the nonlinear problem: where F : IR n IR n..1 Fixed-point iteration F x = 0, 14 The nonlinear equations can be formulated as fixed-point problems: x = Kx, 15 where K is the fixed-point map and is nonlinear, e.g. Kx = x F x. A solution of 16 is called fix-point of the map K. The fix-point iteration is given as: x i+1 = Kx i, 16 and is called nonlinear Richardson iteration, Picard iteration, or the method of successive substitution. Definition.1 Let Ω IR n and let G : Ω IR m. G is Lipschitz continuous on Ω with Lipschitz constant γ if for all x, y Ω. Gx Gy γ x y, 17
5 ITERATIVE SCHEMES FOR COUPLED PROBLEMS 5 For the convergence we have to assume that K be a contraction map on Ω with Lipschitz constant γ < 1. Algorithm.1 We apply the fix-point iterative scheme to decouple the nonseparable Hamiltonian problem and 4. q i = H p p i 1, q i 1, t [t n, t n+1 ] 18 ṗ i = H p i 1, q i 1, t [t n, t n+1 ] 19 pt 0 = p 0, qt 0 = q 0, 0 the starting solutions for the i-th iterative steps are given as : p i 1 t, q i 1 t are the solutions of the i 1 th iterative step and we have the initial condition for the fix-point iteration: p 0 t, q 0 t t = pt n, qt n t We assume that we have convergent results after i = 1..., m iterative steps or with the stopping criterion: max p i+1 p i, q i+1 q i err, while is the Euclidean norm or a simple vector-norm, e.g. L. Iterative Verlet applied to the Hamiltonian and 4 : We start with q 0, p 0 t = qt n, pt n t : The iterative scheme is given as: q i t = qt n + h H p ptn h H p i 1t, q i 1 t, qt n, 1 p i t = pt n h H p i 1t, q i 1 t h H pt n h H p i 1t, q i 1 t, qt n + h H pt n h H p p i 1t, q i 1 t, qt n, where I = or 4. for t [t n, t n+1 ], h = t n+1 t n, n = 0, 1,..., N, i = 1,,..., I, 4 For the fix-point iteration, we have the problem of the initialization, means the start of the iterative scheme. We can improve the starting solution with a preprocessing method, which derives a first improved initial solution. Improve Initialization Process
6 ITERATIVE SCHEMES FOR COUPLED PROBLEMS 6 Algorithm. To improve the initial solution we can start with: 1. We initialize with a result of the explicit Euler-method : q 0, p 0 t = qt n+1 Euler1st, pt n+1 Euler1st t :. We initialize with a result of the explicit RK-method : q 0, p 0 t = qt n+1 RK4th, pt n+1 RK4th t :. Newton s method We solve the nonlinear operator equation 14. While F : D X Y with the Banach spaces X, Y is given with the norms X and Y. Let F be at least once continuous differentiable, further we assume x 0 is a starting solution of the unknown solution x. Then the successive linearization lead to the general Newton s method: F x i x i = F x i, 5 where x i = x i+1 x i and i = 0, 1,,.... The method derive the solution of a nonlinear problem by solving the following algorithm. Algorithm. By considering the sequential splitting method we obtain the following algorithm. We apply the equations q H p p, q = 0, t [tn, t n+1 ] 6 ṗ + H p, q = 0, t [tn, t n+1 ] 7 p n = pt n, q n = qt n, 8 where p = p 1,..., p 6 t and q = q 1,..., q 6 t into the Newtons-formula we have: and we can compute F p, q = ẋ + H x 9 x x k+1 = x k DF x k 1 F x k, 0 where DF x is the Jacobian matrix and k = 0, 1,.... We stop the iterations when we obtain : x k+1 x k err, where err is an error bound, e.g. err = 10 4.
7 SPLITTING METHODS 7 The solution vector F is given as: F x = where x = q 1,..., q 6, p 1,..., p 6 t and F x = Fq x F p x F q,1 x F q, x. F q,6 x F p,1 x F p, x. F p,6 x = q H p ṗ H 1. The Jacobian matrix for the equation system is given as : DF x = F 1 x 1 F 1 x... F x 1 F x.... F 1 x 1 F 1 x... F 1 x 1 F x 1 F 1 x 1 where x = x 1,..., x 1 t = q 1,..., q 6, p 1,..., p 6 t. Splitting Methods In the following, we discuss the different splitting schemes. The simplest such symmetric product is If one naively assumes that T h = S AB h or T h = S BA h. 4 T h = e ta+b + Ch + Dh 4 +, 5 then a Richardson extrapolation would only give a third-order algorithm. 1 [ k k T k h/k T h ] = e ta+b + Oh 4, 6 1
8 4 ITERATIVE MPE METHOD 8 Thus for a given set of n whole numbers {k i } one can have a nth-order approximation e ta+b = n i=1 c i T ki t For orders four to ten, one has explicitly: k i T 4 t = 1 T t + 4 T T 6 t = 1 4 T t T T 8 t = 1 60 T t T t + Oh n+1. 7 t T T 10 t = T t T T 4 t 4 4 Iterative MPE method t, 8 40 T t t T 5 t, T T t 5 t, 40 4 t. 41 Based on the nonlinear problem, we extend the MPE method to an iterative scheme. Algorithm 4.1 Iterative Verlet applied to the Hamiltonian and 4 : We start with q 0, p 0 t = qt n, pt n t : The iterative scheme for computing T p i, q i, h = q i h, p i h t is given as: T p i, q i, h = e h/b e ha e h/b q 0, p 0 t = I + I + I + N 1 1 j! h H p j i 1, q i 1 4 pt n j=1 N 1 H h j! p p j i 1, q i 1 t n j=1 N 1 1 j! h H p jqt i 1, q i 1 n pt n, pt n t j=1 We have the higher order schemes given as:
9 5 NUMERICAL EXAMPLES 9 T 4 p i, q i, h = 1 T p i, q i, h + 4 T p i, q i, h T 6 p i, q i, h = 1 4 T p i, q i, h T p i, q i, h T p i, q i, h T 8 p i, q i, h = 1 60 T p i, q i, h T T p i, q i, h T 4 p i, q i, h p i, q i, h T 10 p i, q i, h = T p i, q i, h T T 4 p i, q i, h 4 p i, q i, h T T p i, q i, h 5 p i, q i, h. 46 where for t [t n, t n+1 ], h = t n+1 t n, n = 0, 1,..., N, i = 1,,..., I, and we have a stopping criterion or a fixed number of iterative steps, e.g. I = or 4. 5 Numerical Examples The Levitron is described on the base of rigid body theory. With the convention of Goldstein [9] for the Euler angles the angular velocity ω φ is along the z-axis of the system, ω θ along the line of nodes and ω ψ along the z -axis. Transforming them into body coordinates one gets φ sin θ sin ψ + θ cos ψ ω = φ sin θ cos ψ + θ sin ψ 47 φ cos θ + ψ Finally the kinetic energy can be written as T = 1 [ mẋ + ẏ + ż + A θ + φ sin θ + C ψ + φ cos θ ] 48 The potential energy U is given by the sum of the gravitational energy and the interaction potential of the Levitron in the magnetic field of the base plate: U = mgz µsin ψ sin θ Φ x + cos ψ sin θ Φ y + cos θ Φ z 49 with mu as the magnetic moment of the top and Φ the magneto-static potential. Following Gans [6] we uses the potential of a ring dipole as approximation for a magnetized plane with a centered unmagnetized hole. Furthermore we
10 5 NUMERICAL EXAMPLES 10 introduced a nondimensionalization for the variables and the magneto-static potential: Z Ψ = 1 + Z / X + Y Z Z Z 7/ Lengths were scaled by the radius R of the base plane, mass were measured in units of m and energy in units of mgh. Therefore the one time unit is R/g. So the dimensionless Hamiltonian with q = X, Y, Z, θ, ψ, φ is given by H = 1 p 1 + p + p + p 4 a + p 5 P 6 cos q 4 51 c M [sin q 4 cos q 5 Ψ 1 sin q 5 Ψ a sin + p 6 q 4 Ψ + cos q 4 ] + q 5 with a and c as the nondimensionalized inertial parameter and M as the ratio of gravitational and magnetic energy. Solving the equations of motion and 4 numerically with the methods described above, one is able to plot the movement of the center of mass like it is shown in figure 1. In the plot the first 5 seconds of a stable trajectory were plotted. With a longer calculation we tested, that the top would levitate for more than one minute z y x y x Figure 1: stable trajectory of the center of mass in D and D The axis in the plot show the nondimensional variables X,Y and Z. The trajectory starts at the equilibrium point q 1, q, q = 0, 0, 1.7. This trajectory was calculated with the fourth-order RungeKutta method with a small timestep of 10 5 units of time, where one time unit is about 59, 5ms, because of the nondimensionalization. For further considerations this trajectory will be used as a reference solution. The errors of the different time-steps with Runge-Kutta are given in Figure. The same equations were solved with the iterative Verlet algorithm described before. Due to the long computation time needed, we simulated only 1000 timesteps and compare the trajectory with the reference solution from the
11 5 NUMERICAL EXAMPLES x error error x t t Figure : Errors of the numerical scheme: Runge-Kutta method explicit 4th order. Runge-Kutta algorithm. In figure is shown how the trajectory of the same initial conditions looks like with the Verlet algorithm y z y x x Figure : trajectory calculated with Verlet algorithm This were done for one, two and four iterations per timestep, to see whether how many iterations are reasonable. The results are shown in Figure 4. In a first comparison, we deal with the second order Verlet algorithm and improve the scheme with iterative steps. In a first initialisation process, one can see that the errors are very similar to 1 or iterative steps. The reduction of the error is possible with the improvement to higher order initialisation scheme, e.g. start with a first approximate solution with a RK scheme, or apply extrapolation schemes. The following tables 1 and should give an impression of the timescales of the problem and the errors. Remark 5.1 Obviously the iterations does not improve the algorithm, when only using a lower order initialisation. By the way, it is sufficient to apply one iterative step in in comparison with the Runge-Kutta algorithm.
12 5 NUMERICAL EXAMPLES 1 Figure 4: Errors of the numerical scheme: Iterative Verlet method. timestep number of steps computing time stability Runge-Kutta min sec ok ok sec ok Table 1: Stability and Computational Time with 4th order explicit RungeKutta. Also we have a benefit in reducing the computational time instead of applying only Runge-Kutta schemes. We tried to improve the solution with a extrapolation scheme in fourth order. We have a view at the errors this algorithm produces in comparison with the Runge-Kutta Solution with small time-steps 10 5 time units per step. In Figure 5, we presented the results of the 4th MPE method with different time-steps and compared it with the Runge-Kutta solution. Figure 5: Errors of the numerical scheme: 4th oder Extrapolation Scheme with Verlet method a Kernel h = Also we tested the 6th order MPE method with different time-steps and
13 6 CONCLUSIONS AND DISCUSSIONS 1 Verlet timestep iterations per step 1 4 stability ok ok ok computing time 67min 10min 19min mean error maximal error Table : Stability and Computational Time with nd order Verlet Scheme. compared it with the Runge-Kutta solution, see Figure 6. Figure 6: Errors of the numerical scheme: 6th order Extrapolation Scheme with Verlet method a Kernel h = Like for Runge-Kutta we want to give an impression of the time scales for this extrapolation schemes, see Table. Extrapolation 4th order Extrapolation 6th order timestep number of steps computing time 14min 14min 9min 7min mean error maximal error Table : Errors and Computational Time with 4th order MPE scheme unsing Verlet Scheme as Kernel. 6 Conclusions and Discussions In the paper, we have presented a model to simulate a Levitron. Based on the given Hamiltonian system, which is nonlinear, we present novel and simpler schemes based on splitting ideas to solve the equation systems. In future, we concentrate on the numerical analysis and embedding higher order splitting kernels to the extraopoation schemes.
14 7 APPENDIX 14 7 Appendix For example, the evolution of any dynamical variable uq, p including q and p themselves is given by the Poisson bracket, t uq, p = For a separable Hamiltonian, A and B are Lie operators, or vector fields u H p u p H = A + Buq, p. 5 Hp, q = p + V q, 54 m A = v where we have abbreviated H p B = aq v 55 H p, q = v = p/m and p, q = aq = V q/m. The exponential operators e ha and e hb are then just shift operators. Sh = e h/b e ha e h/b That is also given as a Verlet-algorithm in the following scheme. We start with q 0, v 0 t = qt n, vt n t : q 1, v 1 t = e h/b q 0, v 0 t = I + 1 h i aq v i q 0, v 0 t 56 = q 0, v haq 0 t, 57 q, v t = e ha q 1, v 1 t = I + h v i q 1, v 1 t 58 i i = q 1 + hv 1, v 1 t, 59 q, v t = e h/b q, v t = I + 1 h i aq v i q, v t 60 = q, v + 1 haq 1 t. 61 And the substitution is given the algorithm for one time-step n n + 1: q, v t = q 0 + hv 0 + h aq 0, v 0 + h aq 0 + h aq 0 + hv 0 + h aq 0 t, 6
15 REFERENCES 15 while qt n+1, vt n+1 t = q, v t. Iterative Verlet Algorithm In the abstract version of H p p i 1, q i 1 q, H p i 1, q i 1. Algorithm 7.1 We have the iterative Verlet Algorithm: 1. We start with the initialisation : p 0 t n+1, q 0 t n+1 t = pt n, qt n t and i = 0. The iterative step is given as: i = i + 1 and we have: q n+1 i = q n + h H p pn 1 h H pn+1 i 1, qn+1 i 1, qn+1 i 1, 6 p i t n+1 = p n h H pn+1 i 1, qn+1 i, 64 p i t n+1 = p n h H pn+1 i 1, qn + h H p pn 1 h H pn+1 i 1, qn+1 i 1, qn+1 i 1, where h = t n+1 t n is the local time-step. We compute the stopping criterion: max pt n+1 i p n+1 i 1, qtn+1 i q n+1 i 1 err or we stop after i= I, while I is the maximal iterative step.. The result is given as: pt n+1, qt n+1 t = p i t n+1, q i t n+1 t and n = n + 1 if n > N, while N is the maximal time-step, we stop else we go to step 1. References [1] S. Chin and J. Geiser. Multi-product operator splitting as a general method of solving autonomous and non-autonomous equations. IMA J. Numer. Anal., first published online January 1, 011. [] B. Davis. Integral Transform and Their Applications. Applied Mathematical Sciences, 5, Springer Verlag, New York, Heidelberg, Berlin, [] H.R. Dullin and R. Easton. Stability of Levitron. Physica D: Nonlinear Phenomena, vol. 16, no. 1-, 1-17, [4] K.-J. Engel, R. Nagel, One-Parameter Semigroups for Linear Evolution Equations. Springer-Verlag, Heidelberg, New York, 000. [5] I. Farago and J. Geiser. Iterative Operator-Splitting Methods for Linear Problems. Preprint No. 104 of the Weierstrass Institute for Applied Analysis and Stochastics, International Journal of Computational Science and Engineering, accepted September 007.
16 REFERENCES 16 [6] R.F. Gans, T.B. Jones, and M. Washizu. Dynamics of the Levitron. J. Phys. D., 1, , [7] J. Geiser. Higher order splitting methods for differential equations: Theory and applications of a fourth order method. Numerical Mathematics: Theory, Methods and Applications. Global Science Press, Hong Kong, China, accepted, April 008. [8] J. Geiser and L. Noack. Iterative operator-splitting methods for nonlinear differential equations and applications of deposition processes Preprint 008-4, Humboldt University of Berlin, Department of Mathematics, Germany, 008. [9] H. Goldstein, Ch.P. Poole, and J. Safko. Classical mechanics. Addison Wesley, San Francisco, USA, 00. [10] M. Hieber, A. Holderrieth and F. Neubrander. Regularized semigroups and systems of linear partial differential equations. Annali della Scuola Normale Superiore di Pisa - Classe di Scienze, Ser.4, 19:, 6-79, 199. [11] M. Hochbruck and A. Ostermann. Explicit Exponential Runge-Kutta Methods for Semilinear Parabolic Problems. SIAM Journal on Numerical Analysis, 4:, , 005. [1] E. Hansen and A. Ostermann. Exponential splitting for unbounded operators. Mathematics of Computation, accepted, 008. [1] Hildebrand, F.B Introduction to Numerical Analysis. Second Edition, Dover Edition. [14] T. Jahnke and C. Lubich. Error bounds for exponential operator splittings. BIT Numerical Mathematics, 40:4, , 000. [15] J. Kanney, C. Miller and C. Kelley. Convergence of iterative split-operator approaches for approximating nonlinear reactive transport problems. Advances in Water Resources, 6:47 61, 00. [16] C.T. Kelly. Iterative Methods for Linear and Nonlinear Equations. Frontiers in Applied Mathematics, SIAM, Philadelphia, USA, [17] G.I Marchuk. Some applications of splitting-up methods to the solution of problems in mathematical physics. Aplikace Matematiky, 1, 10-1, [18] G. Strang. On the construction and comparison of difference schemes. SIAM J. Numer. Anal., 5, , [19] I. Najfeld and T.F. Havel. Derivatives of the matrix exponential and their computation. Adv. Appl. Math, ftp://ftp.das.harvard.edu/pub/cheatham/tr--94.ps.gz, 1995.
17 REFERENCES 17 [0] Strang, G On the construction and comparison of difference schemes. SIAM J. Numer. Anal., 5, [1] M. Suzuki. General Decomposition Theory of Ordered Exponentials. Proc. Japan Acad., 69, Ser. B, 161, 199. [] H. Goldstein. Klassische Mechanik. Akademische Verlagsgesellschaft, Wiesbaden, 1981.
THE Hamilton equations of motion constitute a system
Proceedings of the World Congress on Engineering 0 Vol I WCE 0, July 4-6, 0, London, U.K. Systematic Improvement of Splitting Methods for the Hamilton Equations Asif Mushtaq, Anne Kværnø, and Kåre Olaussen
More informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
More informationHigh-order actions and their applications to honor our friend and
High-order actions and their applications to honor our friend and collaborator Siu A. Chin Universitat Politecnica de Catalunya, Barcelona, Spain Lecture: Magnus Expansion and Suzuki s Method Outline of
More informationMulti-product expansion, Suzuki s method and the Magnus integrator for solving time-dependent problems
Multi-product expansion, Suzuki s method and the Magnus integrator for solving time-dependent problems Jürgen Geiser and Siu A. Chin geiser@mathematik.hu-berlin.de Abstract. In this paper we discuss the
More informationPhysical Dynamics (SPA5304) Lecture Plan 2018
Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle
More informationAppell-Hamel dynamical system: a nonlinear test of the Chetaev and the vakonomic model
ZAMM Z. Angew. Math. Mech. 87, No. 10, 692 697 (2007) / DOI 10.1002/zamm.200710344 Appell-Hamel dynamical system: a nonlinear test of the Chetaev and the vakonomic model Shan-Shan Xu 1, Shu-Min Li 1,2,
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 informationSANGRADO PAGINA 17CMX24CM. PhD Thesis. Splitting methods for autonomous and non-autonomous perturbed equations LOMO A AJUSTAR (AHORA 4CM)
SANGRADO A3 PAGINA 17CMX24CM LOMO A AJUSTAR (AHORA 4CM) We have considered the numerical integration of non-autonomous separable parabolic equations using high order splitting methods with complex coefficients.
More informationSemi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations
Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang
More informationHigher-Order Difference and Higher-Order Splitting Methods for 2D Parabolic Problems with Mixed Derivatives
International Mathematical Forum, 2, 2007, no. 67, 3339-3350 Higher-Order Difference and Higher-Order Splitting Methods for 2D Parabolic Problems with Mixed Derivatives Jürgen Geiser Department of Mathematics
More informationTHEWAY in which time evolution acts on observables may be rephrased in
Appendix F Implementing evolution F.1 Koopmania THEWAY in which time evolution acts on observables may be rephrased in the language of functional analysis, by introducing the Koopman operator, whose action
More informationIntroduction to symplectic integrator methods for solving diverse dynamical equations
Introduction to symplectic integrator methods for solving diverse dynamical equations Siu A. Chin Department of Physics and Astronomy Texas A&M University College Station, TX, 77840, USA Third Lecture
More informationEfficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators
MATEMATIKA, 8, Volume 3, Number, c Penerbit UTM Press. All rights reserved Efficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators Annie Gorgey and Nor Azian Aini Mat Department of Mathematics,
More information1 Ordinary differential equations
Numerical Analysis Seminar Frühjahrssemester 08 Lecturers: Prof. M. Torrilhon, Prof. D. Kressner The Störmer-Verlet method F. Crivelli (flcrivel@student.ethz.ch May 8, 2008 Introduction During this talk
More informationPhysics 115/242 Comparison of methods for integrating the simple harmonic oscillator.
Physics 115/4 Comparison of methods for integrating the simple harmonic oscillator. Peter Young I. THE SIMPLE HARMONIC OSCILLATOR The energy (sometimes called the Hamiltonian ) of the simple harmonic oscillator
More informationPHY 5246: Theoretical Dynamics, Fall November 16 th, 2015 Assignment # 11, Solutions. p θ = L θ = mr2 θ, p φ = L θ = mr2 sin 2 θ φ.
PHY 5246: Theoretical Dynamics, Fall 215 November 16 th, 215 Assignment # 11, Solutions 1 Graded problems Problem 1 1.a) The Lagrangian is L = 1 2 m(ṙ2 +r 2 θ2 +r 2 sin 2 θ φ 2 ) V(r), (1) and the conjugate
More informationSYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS
BIT 0006-3835/00/4004-0726 $15.00 2000, Vol. 40, No. 4, pp. 726 734 c Swets & Zeitlinger SYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS E. HAIRER Section de mathématiques, Université
More informationReducing round-off errors in symmetric multistep methods
Reducing round-off errors in symmetric multistep methods Paola Console a, Ernst Hairer a a Section de Mathématiques, Université de Genève, 2-4 rue du Lièvre, CH-1211 Genève 4, Switzerland. (Paola.Console@unige.ch,
More informationIntroduction to numerical methods for Ordinary Differential Equations
Introduction to numerical methods for Ordinary Differential Equations Charles-Edouard Bréhier Abstract The aims of these lecture notes are the following. We introduce Euler numerical schemes, and prove
More informationSMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract
SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm Abstract Dedicated to Prof. Michele Caputo on the occasion of his 8th birthday We consider ordinary fractional
More informationTheoretical physics. Deterministic chaos in classical physics. Martin Scholtz
Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton
More informationRunga-Kutta Schemes. Exact evolution over a small time step: Expand both sides in a small time increment: d(δt) F (x(t + δt),t+ δt) Ft + FF x ( t)
Runga-Kutta Schemes Exact evolution over a small time step: x(t + t) =x(t)+ t 0 d(δt) F (x(t + δt),t+ δt) Expand both sides in a small time increment: x(t + t) =x(t)+x (t) t + 1 2 x (t)( t) 2 + 1 6 x (t)+
More informationAccuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations
COMMUICATIOS I COMPUTATIOAL PHYSICS Vol. 5, o. -4, pp. 779-79 Commun. Comput. Phys. February 009 Accuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations Tao
More informationClassical Mechanics Comprehensive Exam Solution
Classical Mechanics Comprehensive Exam Solution January 31, 011, 1:00 pm 5:pm Solve the following six problems. In the following problems, e x, e y, and e z are unit vectors in the x, y, and z directions,
More informationAssignment 8. [η j, η k ] = J jk
Assignment 8 Goldstein 9.8 Prove directly that the transformation is canonical and find a generating function. Q 1 = q 1, P 1 = p 1 p Q = p, P = q 1 q We can establish that the transformation is canonical
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationNotes on symplectic geometry and group actions
Notes on symplectic geometry and group actions Peter Hochs November 5, 2013 Contents 1 Example of Hamiltonian mechanics 1 2 Proper actions and Hausdorff quotients 4 3 N particles in R 3 7 4 Commuting actions
More informationPhysics 351 Wednesday, February 14, 2018
Physics 351 Wednesday, February 14, 2018 HW4 due Friday. For HW help, Bill is in DRL 3N6 Wed 4 7pm. Grace is in DRL 2C2 Thu 5:30 8:30pm. Respond at pollev.com/phys351 or text PHYS351 to 37607 once to join,
More information10/24/07. Optics Rotation: Levitron. Nicki Bornhauser
10/24/07 Optics Rotation: Levitron Nicki Bornhauser Levitron - Overview What is the Levitron? Theory Basic principles An simple model Experiment Magnetic field Levitation Results of measurement Wobbling
More informationFourth-Order Splitting Methods for Time-Dependant Differential Equations
NUMERICAL MATHEMATICS: Theory, Methods and Applications Numer. Math. Theor. Meth. Appl., Vol. 1, No. 3, pp. 321-339 (2008) Fourth-Order Splitting Methods for Time-Dependant Differential Equations Jürgen
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
More informationExponential Integrators
Exponential Integrators John C. Bowman and Malcolm Roberts (University of Alberta) June 11, 2009 www.math.ualberta.ca/ bowman/talks 1 Outline Exponential Integrators Exponential Euler History Generalizations
More informationExponential Integrators
Exponential Integrators John C. Bowman (University of Alberta) May 22, 2007 www.math.ualberta.ca/ bowman/talks 1 Exponential Integrators Outline Exponential Euler History Generalizations Stationary Green
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationOn nonexpansive and accretive operators in Banach spaces
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 10 (2017), 3437 3446 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa On nonexpansive and accretive
More informationAn Exactly Solvable 3 Body Problem
An Exactly Solvable 3 Body Problem The most famous n-body problem is one where particles interact by an inverse square-law force. However, there is a class of exactly solvable n-body problems in which
More informationGeneralized mid-point and trapezoidal rules for Lipschitzian RHS
Generalized mid-point and trapezoidal rules for Lipschitzian RHS Richard Hasenfelder, Andreas Griewank, Tom Streubel Humboldt Universität zu Berlin 7th International Conference on Algorithmic Differentiation
More informationTHE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS
THE MAGNUS METHOD FOR SOLVING OSCILLATORY LIE-TYPE ORDINARY DIFFERENTIAL EQUATIONS MARIANNA KHANAMIRYAN Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road,
More informationNumerical Methods for the Landau-Lifshitz-Gilbert Equation
Numerical Methods for the Landau-Lifshitz-Gilbert Equation L ubomír Baňas Department of Mathematical Analysis, Ghent University, 9000 Gent, Belgium lubo@cage.ugent.be http://cage.ugent.be/~lubo Abstract.
More informationSpurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Spurious Chaotic Solutions of Dierential Equations Sigitas Keras DAMTP 994/NA6 September 994 Department of Applied Mathematics and Theoretical Physics
More informationLeap Frog Solar System
Leap Frog Solar System David-Alexander Robinson Sch. 08332461 20th November 2011 Contents 1 Introduction & Theory 2 1.1 The Leap Frog Integrator......................... 2 1.2 Class.....................................
More informationNONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis
LECTURE 3: FLOWS NONLINEAR DYNAMICS AND CHAOS Patrick E McSharr Sstems Analsis, Modelling & Prediction Group www.eng.o.ac.uk/samp patrick@mcsharr.net Tel: +44 83 74 Numerical integration Stabilit analsis
More informationMCE 366 System Dynamics, Spring Problem Set 2. Solutions to Set 2
MCE 366 System Dynamics, Spring 2012 Problem Set 2 Reading: Chapter 2, Sections 2.3 and 2.4, Chapter 3, Sections 3.1 and 3.2 Problems: 2.22, 2.24, 2.26, 2.31, 3.4(a, b, d), 3.5 Solutions to Set 2 2.22
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationM3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011
M3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011 Text for the course: Professor Darryl D Holm 25 October 2011 Imperial College London d.holm@ic.ac.uk http://www.ma.ic.ac.uk/~dholm/ Geometric Mechanics
More informationECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability
ECE 4/5 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability Spring 014 Instructor: Kai Sun 1 Transient Stability The ability of the power system to maintain synchronism
More informationReducing round-off errors in rigid body dynamics
Reducing round-off errors in rigid body dynamics Gilles Vilmart INRIA Rennes, ENS Cachan Bretagne, avenue Robert Schuman, 357 Bruz, France Université de Genève, Section de mathématiques, -4 rue du Lièvre,
More informationHamiltonian flow in phase space and Liouville s theorem (Lecture 5)
Hamiltonian flow in phase space and Liouville s theorem (Lecture 5) January 26, 2016 90/441 Lecture outline We will discuss the Hamiltonian flow in the phase space. This flow represents a time dependent
More informationBackward error analysis
Backward error analysis Brynjulf Owren July 28, 2015 Introduction. The main source for these notes is the monograph by Hairer, Lubich and Wanner [2]. Consider ODEs in R d of the form ẏ = f(y), y(0) = y
More informationCanonical Forms for BiHamiltonian Systems
Canonical Forms for BiHamiltonian Systems Peter J. Olver Dedicated to the Memory of Jean-Louis Verdier BiHamiltonian systems were first defined in the fundamental paper of Magri, [5], which deduced the
More informationDIFFERENTIABLE SEMIFLOWS FOR DIFFERENTIAL EQUATIONS WITH STATE-DEPENDENT DELAYS
UNIVERSITATIS IAGELLONICAE ACTA MATHEMATICA, FASCICULUS XLI 23 DIFFERENTIABLE SEMIFLOWS FOR DIFFERENTIAL EQUATIONS WITH STATE-DEPENDENT DELAYS by Hans-Otto Walther Careful modelling of systems governed
More informationModified Equations for Variational Integrators
Modified Equations for Variational Integrators Mats Vermeeren Technische Universität Berlin Groningen December 18, 2018 Mats Vermeeren (TU Berlin) Modified equations for variational integrators December
More informationThe Koopmanoperator action on a state space function a(x) is to replace it by
Appendix A9 Implementing evolution A9. Koopmania The Koopmanoperator action on a state space function a(x) is to replace it by its downstream value time t later, a(x) a(x(t)), evaluated at the trajectory
More informationThe Milne error estimator for stiff problems
13 R. Tshelametse / SAJPAM. Volume 4 (2009) 13-28 The Milne error estimator for stiff problems Ronald Tshelametse Department of Mathematics University of Botswana Private Bag 0022 Gaborone, Botswana. E-mail
More informationDIMENSION SPLITTING FOR TIME DEPENDENT OPERATORS. Eskil Hansen and Alexander Ostermann
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume X, Number 0X, XX 00X pp. 1 11 DIMENSION SPLITTING FOR TIME DEPENDENT OPERATORS Eskil Hansen and Alexander Ostermann Institut
More informationC Spoke k. Spoke k+1
Numerical Accuracy Case Studies: D Rimless Spoked Wheel and a 3D Passive-Dynamic Model of Walking Michael J. Coleman Andy Ruina Anindya Chatterjee 1 Department of Theoretical and Applied Mechanics Cornell
More informationNonlocal problems for the generalized Bagley-Torvik fractional differential equation
Nonlocal problems for the generalized Bagley-Torvik fractional differential equation Svatoslav Staněk Workshop on differential equations Malá Morávka, 28. 5. 212 () s 1 / 32 Overview 1) Introduction 2)
More informationChapter 11 ORDINARY DIFFERENTIAL EQUATIONS
Chapter 11 ORDINARY DIFFERENTIAL EQUATIONS The general form of a first order differential equations is = f(x, y) with initial condition y(a) = y a We seek the solution y = y(x) for x > a This is shown
More informationCOMPLETE ALL ROUGH WORKINGS IN THE ANSWER BOOK AND CROSS THROUGH ANY WORK WHICH IS NOT TO BE ASSESSED.
BSc/MSci EXAMINATION PHY-304 Time Allowed: Physical Dynamics 2 hours 30 minutes Date: 28 th May 2009 Time: 10:00 Instructions: Answer ALL questions in section A. Answer ONLY TWO questions from section
More informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More informationNumerical Integration of Equations of Motion
GraSMech course 2009-2010 Computer-aided analysis of rigid and flexible multibody systems Numerical Integration of Equations of Motion Prof. Olivier Verlinden (FPMs) Olivier.Verlinden@fpms.ac.be Prof.
More informationSolving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions.
Journal of Mathematical Modeling Vol 1, No 1, 213, pp 28-4 JMM Solving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions Farshid
More informationIntegration of Vlasov-type equations
Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Plasma the fourth state of matter 99% of the visible matter in the universe is made of plasma
More informationStudy the Numerical Methods for Solving System of Equation
Study the Numerical Methods for Solving System of Equation Ravi Kumar 1, Mr. Raj Kumar Duhan 2 1 M. Tech. (M.E), 4 th Semester, UIET MDU Rohtak 2 Assistant Professor, Dept. of Mechanical Engg., UIET MDU
More informationNumerical Simulations of N Point Vortices on Bounded Domains
Cornell University Mathematics Department Senior Thesis Numerical Simulations of N Point Vortices on Bounded Domains Author: Sy Yeu Chou Bachelor of Arts May 2014, Cornell University Thesis Advisor: Robert
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 2: Basic tools and concepts Nonlinear Single-Particle Dynamics in High Energy Accelerators This course consists of eight lectures: 1.
More informationAsymptotic stability of an evolutionary nonlinear Boltzmann-type equation
Acta Polytechnica Hungarica Vol. 14, No. 5, 217 Asymptotic stability of an evolutionary nonlinear Boltzmann-type equation Roksana Brodnicka, Henryk Gacki Institute of Mathematics, University of Silesia
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationExponential integrators
Exponential integrators Marlis Hochbruck Heinrich-Heine University Düsseldorf, Germany Alexander Ostermann University Innsbruck, Austria Helsinki, May 25 p.1 Outline Time dependent Schrödinger equations
More informationEstimating Accuracy in Classical Molecular Simulation
Estimating Accuracy in Classical Molecular Simulation University of Illinois Urbana-Champaign Department of Computer Science Institute for Mathematics and its Applications July 2007 Acknowledgements Ben
More informationHamiltonian Dynamics
Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationPH 610/710-2A: Advanced Classical Mechanics I. Fall Semester 2007
PH 610/710-2A: Advanced Classical Mechanics I Fall Semester 2007 Time and location: Tuesdays & Thursdays 8:00am 9:15am (EB 144) Instructor and office hours: Dr. Renato Camata, camata@uab.edu CH 306, (205)
More informationHopf Bifurcation and Limit Cycle Analysis of the Rikitake System
ISSN 749-3889 (print), 749-3897 (online) International Journal of Nonlinear Science Vol.4(0) No.,pp.-5 Hopf Bifurcation and Limit Cycle Analysis of the Rikitake System Xuedi Wang, Tianyu Yang, Wei Xu Nonlinear
More informationGEOMETRIC INTEGRATION OF ORDINARY DIFFERENTIAL EQUATIONS ON MANIFOLDS
BIT 0006-3835/01/4105-0996 $16.00 2001, Vol. 41, No. 5, pp. 996 1007 c Swets & Zeitlinger GEOMETRIC INTEGRATION OF ORDINARY DIFFERENTIAL EQUATIONS ON MANIFOLDS E. HAIRER Section de mathématiques, Université
More informationDunkl operators and Clifford algebras II
Translation operator for the Clifford Research Group Department of Mathematical Analysis Ghent University Hong Kong, March, 2011 Translation operator for the Hermite polynomials Translation operator for
More informationRigid body simulation. Once we consider an object with spatial extent, particle system simulation is no longer sufficient
Rigid body dynamics Rigid body simulation Once we consider an object with spatial extent, particle system simulation is no longer sufficient Rigid body simulation Unconstrained system no contact Constrained
More informationMULTI-VARIABLE POLYNOMIAL SOLUTIONS TO PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS
MULTI-VARIABLE POLYNOMIAL SOLUTIONS TO PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS J. MC LAUGHLIN Abstract. Solving Pell s equation is of relevance in finding fundamental units in real
More informationAN EXISTENCE-UNIQUENESS THEOREM FOR A CLASS OF BOUNDARY VALUE PROBLEMS
Fixed Point Theory, 13(2012), No. 2, 583-592 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.html AN EXISTENCE-UNIQUENESS THEOREM FOR A CLASS OF BOUNDARY VALUE PROBLEMS M.R. MOKHTARZADEH, M.R. POURNAKI,1 AND
More informationNUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS
TWMS J Pure Appl Math V5, N2, 24, pp22-228 NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS S ASADI, AH BORZABADI Abstract In this paper, Haar wavelet benefits are applied to the delay
More informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 14, 2009) I. INTRODUCTION One frequently obtains detailed dynamical
More informationu = (A + F )u, u(0) = η, (1.1)
A CONVERGENCE ANALYSIS OF THE PEACEMAN RACHFORD SCHEME FOR SEMILINEAR EVOLUTION EQUATIONS ESKIL HANSEN AND ERIK HENNINGSSON Abstract. The Peaceman Rachford scheme is a commonly used splitting method for
More informationMagnetostatics III Magnetic Vector Potential (Griffiths Chapter 5: Section 4)
Dr. Alain Brizard Electromagnetic Theory I PY ) Magnetostatics III Magnetic Vector Potential Griffiths Chapter 5: Section ) Vector Potential The magnetic field B was written previously as Br) = Ar), 1)
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationSplitting Approach to Coupled Navier-Stokes. and Molecular Dynamics Equations
Journal of Computations & Modelling, vol., no.,, -33 ISSN: 79-765 (print), 79-885 (online) Scienpress Ltd, Splitting Approach to Coupled Navier-Stokes and Molecular Dynamics Equations Jürgen Geiser Abstract
More informationNumerical Analysis II. Problem Sheet 8
P. Grohs S. Hosseini Ž. Kereta Spring Term 15 Numerical Analysis II ETH Zürich D-MATH Problem Sheet 8 Problem 8.1 The dierential equation Autonomisation and Strang-Splitting ẏ = A(t)y, y(t ) = y, A(t)
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationStability of Equilibrium Positions of Mechanical Systems with Switched Force Fields
SCIETIFIC PUBLICATIOS OF THE STATE UIVERSITY OF OVI PAZAR SER. A: APPL. MATH. IFORM. AD MECH. vol. 4, 2 2012, 35-39 Stability of Equilibrium Positions of Mechanical Systems with Switched Force Fields A.
More informationMechatronics 1: ME 392Q-6 & 348C 31-Aug-07 M.D. Bryant. Analogous Systems. e(t) Se: e. ef = p/i. q = p /I, p = " q C " R p I + e(t)
V + - K R + - - k b V R V L L J + V C M B Analogous Systems i = q. + ω = θ. C -. λ/l = q v = x F T. Se: e e(t) e = p/i R: R 1 I: I e C = q/c C = dq/dt e I = dp/dt Identical dierential equations & bond
More informationKinematics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Kinematics Semester 1, / 15
Kinematics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-17 B. Bona (DAUIN) Kinematics Semester 1, 2016-17 1 / 15 Introduction The kinematic quantities used to represent a body frame are: position
More informationDynamical Low-Rank Approximation to the Solution of Wave Equations
Dynamical Low-Rank Approximation to the Solution of Wave Equations Julia Schweitzer joint work with Marlis Hochbruck INSTITUT FÜR ANGEWANDTE UND NUMERISCHE MATHEMATIK 1 KIT Universität des Landes Baden-Württemberg
More informationHigh-order time-splitting methods for Schrödinger equations
High-order time-splitting methods for Schrödinger equations and M. Thalhammer Department of Mathematics, University of Innsbruck Conference on Scientific Computing 2009 Geneva, Switzerland Nonlinear Schrödinger
More informationDISCRETE METHODS AND EXPONENTIAL DICHOTOMY OF SEMIGROUPS. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXXIII, 2(2004), pp. 97 205 97 DISCRETE METHODS AND EXPONENTIAL DICHOTOMY OF SEMIGROUPS A. L. SASU Abstract. The aim of this paper is to characterize the uniform exponential
More informationEnergy and Equations of Motion
Energy and Equations of Motion V. Tanrıverdi tanriverdivedat@googlemail.com Physics Department, Middle East Technical University, Ankara / TURKEY Abstract. From the total time derivative of energy an equation
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationq 1 F m d p q 2 Figure 1: An automated crane with the relevant kinematic and dynamic definitions.
Robotics II March 7, 018 Exercise 1 An automated crane can be seen as a mechanical system with two degrees of freedom that moves along a horizontal rail subject to the actuation force F, and that transports
More informationHamiltonian aspects of fluid dynamics
Hamiltonian aspects of fluid dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS 01/29/08, 01/31/08 Joris Vankerschaver (CDS) Hamiltonian aspects of fluid dynamics 01/29/08, 01/31/08 1 / 34 Outline
More informationX i t react. ~min i max i. R ij smallest. X j. Physical processes by characteristic timescale. largest. t diff ~ L2 D. t sound. ~ L a. t flow.
Physical processes by characteristic timescale Diffusive timescale t diff ~ L2 D largest Sound crossing timescale t sound ~ L a Flow timescale t flow ~ L u Free fall timescale Cooling timescale Reaction
More informationTechnische Universität Berlin
Technische Universität Berlin Institut für Mathematik M7 - A Skateboard(v1.) Andreas Steinbrecher Preprint 216/8 Preprint-Reihe des Instituts für Mathematik Technische Universität Berlin http://www.math.tu-berlin.de/preprints
More information