norges teknisk-naturvitenskapelige universitet Preserving energy resp. dissipation in numerical PDEs, using the Average Vector Field method
|
|
- Jemimah Pope
- 5 years ago
- Views:
Transcription
1 norges teknisk-naturvitenskapelige universitet Preserving energy resp dissipation in numerical PDEs, using the Average Vector Field method by E Celledoni,V Grimm, RI McLachlan, DI McLaren, DRJ O Neale, B Owren, and GRW Quispel preprint numerics no 7/9 norwegian university of science and technology trondheim, norway This report has URL Address: Department of Mathematical Sciences, Norwegian University of Science and Technology, N-7491 Trondheim, Norway
2
3 Preserving energy resp dissipation in numerical PDEs, using the Average Vector Field method E Celledoni,V Grimm, RI McLachlan, DI McLaren, DRJ O Neale, B Owren, and GRW Quispel September 4, 9 We give a systematic method for discretising Hamiltonian partial differential equations (PDEs) with constant symplectic structure, while preserving their total energy exactly The same method, applied to PDEs with dissipative constant structure, also preserves the correct monotonic decrease of energy 1
4 BIT Numerical Mathematics (6)46:- c Springer 6 DOI:117/s x PRESERVING ENERGY RESP DISSIPATION IN NUMERICAL PDEs, USING THE AVERAGE VECTOR FIELD METHOD E CELLEDONI 1, V GRIMM, RI MCLACHLAN 3, DI MCLAREN 4, DRJ O NEALE 5, B OWREN 6, GRW QUISPEL 7 1 Department of Mathematical Sciences, NTNU 7491 Trondheim, Norway elenac@mathntnuno Department of Mathematics, Friedrich-Alexander-Universität Erlangen-Nürnberg, Martensstr 3, D 9158 Erlangen, Germany grimm@amuni-erlangende Abstract 3 Institute of Fundamental Sciences, Massey University, Private Bag 11- Palmerston North, New Zealand rmclachlan@masseyacnz 4 Department of Mathematics and Statistics, La Trobe University, Victoria 386, Australia DMcLaren@latrobeeduau 5 Department of Mathematics and Statistics, La Trobe University, Victoria 386, Australia DO Neale@latrobeeduau 6 Department of Mathematical Sciences, NTNU 7491 Trondheim, Norway BrynjulfOwren@mathntnuno 7 Department of Mathematics and Statistics, La Trobe University, Victoria 386, Australia RQuispel@latrobeeduau We give a systematic method for discretising Hamiltonian partial differential equations (PDEs) with constant symplectic structure, while preserving their energy exactly The same method, applied to PDEs with constant dissipative structure, also preserves the correct monotonic decrease of energy The method is illustrated by many examples In the Hamiltonian case these include: the sine-gordon, Korteweg-de Vries, nonlinear Schrödinger, (linear) -dependent Schrödinger, and Maxwell equations In the dissipative case the examples are: the Allen-Cahn, Cahn-Hilliard, Ginzburg-Landau, and Heat equations AMS subject classification (): 65L1, 65M6, 65N, 65P1 Key words: Average vector field method, Hamiltonian PDEs, dissipative PDEs, integration 1 Introduction The opening line of Anna Karenina, All happy families resemble one another, but each unhappy family is unhappy in its own way, is a useful metaphor for the relationship between computational ordinary differential equations (ODEs) and computational partial differential equations (PDEs) ODEs are a happy family perhaps they do not resemble each other, but, at the very least, we can treat Received ** Revised ** Communicated by **
5 E CELLEDONI ET AL them by a relatively small compendium of computational techniques PDEs are a huge and motley collection of problems, each unhappy in its own way (Quote from A Iserles book [8]) Whereas there is much truth in the above quote, in this paper we set out to convince the reader that, as far as conservation or dissipation of energy is concerned, many PDEs form part of one big happy family (cf also [9]) that, after a very straightforward and uniform semi-discretisation, may actually be solved by a single unique geometric integration method the so-called average vector field method while preserving the correct conservation, respectively, dissipation of energy The concept of energy has far-reaching importance throughout the physical sciences [4] Therefore a single procedure, as presented here, that correctly conserves, resp dissipates, energy for linear as well as nonlinear, low-order as well as high-order, PDEs would seem to be worth while We consider evolutionary PDEs with independent variables (x, t) R d R, functions u belonging to a Banach space B with values 1 u(x, t) R m, and PDEs of the form (11) u = D δh δu, where D is a constant linear differential operator, the dot denotes t, and (1) H[u] = H(x; u (n) ) dx where Ω is a subset of R d δh R, and dx = dx 1 dx dx d δu derivative of H in the sense that d (13) dɛ H[u + ɛv] δh ɛ= = δu v dx, for all u, v B (cf [16]) For example, if d = m = 1, (14) H[u] = H(x; u, u x,u xx, ) dx, then (15) Ω Ω δh δu = H ( ) ( ) H H u x + x, u x u xx when the boundary terms are zero Similarly, for general d and m, we obtain (16) δh = H δu l u l d k=1 x k ( H u l,k Ω ) +, l =1,, m is the variational 1 Although it is generally real-valued, the function u may also be complex-valued, for example, the nonlinear Schrödinger equation
6 ENERGY CONSERVING PDES 3 We consider Hamiltonian systems of the form (11), where D is a constant skew symmetric operator (cf [16]) and H the energy (Hamiltonian) In this case, we prefer to designate the differential operator in (11) with S instead of D The PDE preserves the energy because S is skew-adjoint with respect to the L inner product, ie (17) usu dx =, u B Ω The system (11) has I : B R as an integral if İ = δi Ω δu S δh δu dx = Integrals C with D δc δu = are called Casimirs Besides PDEs of type (11) where D is skew-adjoint, we also consider PDEs of type (11) where D is a constant negative (semi)definite operator with respect to the L inner product, ie (18) udu dx, u B Ω In this case, we prefer to designate the differential operator D with N and the function H is a Lyapunov function, since then the system (11), ie (19) u = N δh δu, δh δu N δh has H as a Lyapunov function, ie Ḣ = Ω δu dx We will refer to systems (11) with a skew-adjoint S and an energy H as conservative and to systems (11) with a negative (semi)definite operator N and a Lyapunov function H as dissipative Conservative PDEs (11) can be semi-discretised in skew-gradient form (11) u = S H(u), S T = S, when D = S is skew-adjoint u R k, and here, and in the following, we will always denote the discretisations with bars H is chosen in such a way that H x is an approximation to H Lemma 11 Let (111) H[u] = Ω H(x; u (n) )dx, and let H x be any consistent (finite difference) approximation to H (where x := x 1 x x d ) Then the discrete analogue of the variational derivative δh is given by H The proof is given in the appendix δu It is worth noting that the above lemma also applies directly when the approximation to H is obtained by a spectral discretization, since such an approximation can be viewed as a finite difference approximation where the finite difference
7 4 E CELLEDONI ET AL stencil has the same number of entries as the number of grid points on which it is defined The operator is the standard gradient, which replaces the variational derivative because we are now working in a finite (although large) number of dimensions (cf eg (16)) When dealing with (semi-)discrete systems we use the notation u j,n where the index j corresponds to increments in space and n to increments in That is, the point u j,n is the discrete equivalent of u(a + j x, t + n t) where x [a, b] and where t is the initial In most of the equations we present, one of the indices is held constant, in which case, for simplicity, we drop it from the notation For example, we use u j to refer to the values of u at different points in space and at a fixed level Theorem 1 Let S (resp N ) be any consistent constant skew (resp negative-definite) matrix approximation to S (resp N ) Let H x be any consistent (finite difference) approximation to H Finally, let (11) f(u) := S H(u) (resp f(u) := N H(u)), and let u n be the solution of the average vector field (AVF) method (113) u n+1 u n t = 1 f((1 ξ)u n + ξu n+1 ) dξ, applied to equation (11) Then the semidiscrete energy H is preserved exactly (resp dissipated monotonically): H(u n+1 )=H(u n ) (resp H(u n+1 ) H(u n )) H is preserved since (114) H = ( H ) T S H = Discretisations of this type can be given for pseudospectral, finite-element, Galerkin and finite-difference methods (cf [1, 13]); for simplicity s sake, we will concentrate on finite-difference methods, though we include one example of a pseudospectral method for good measure The AVF method was recently [17] shown to preserve the energy H exactly for any vector field f of the form f(u) =S H(u), where H is an arbitrary function, and S is any constant skew matrix The AVF method is related to discrete gradient methods (cf [11]) If D is a constant negative-definite operator, then the dissipative PDE (11) can be discretized in the form (115) u = N H(u), The relationship of (113) to Runge-Kutta methods was explored in [3]
8 ENERGY CONSERVING PDES 5 where N is a negative (semi)definite matrix and H is a discretisation as above That is, H is a Lyapunov-function for the semi-discretized system, since (116) H = ( H ) T N H The AVF method (113) again preserves this structure, ie we have (117) H(u n+1 ) H(u n ), and H is a Lyapunov function for the discrete system Taking the scalar product of (113) with 1 H((1 ξ)u n + ξu n+1 ) dξ on both sides of the equation yields (118) ie (119) and therefore (1) 1 t 1 (u n+1 u n ) H((1 ξ)u n + ξu n+1 ) dξ, 1 1 d t dξ H((1 ξ)u n + ξu n+1 ) dξ, 1 t (H(u n+1) H(u n )) Our purpose is to show that the procedure described above, namely 1 Discretize the energy functional H using any (consistent) approximation H x Discretize D by a constant skew-symmetric (resp negative (semi)definite) matrix 3 Apply the AVF method can be generally applied and leads, in a systematic way, to energy-preserving methods for conservative PDEs and energy-dissipating methods for dissipative PDEs We shall demonstrate the procedure by going through several well-known nonlinear and linear PDEs step by step In particular we give examples of how to discretise nonlinear conservative PDEs (in subsection 1), linear conservative PDEs (in subsection ), nonlinear dissipative PDEs (in subsection 31), and linear dissipative PDEs (in subsection 3) Conservative PDEs 1 Nonlinear conservative PDEs Example 1 Sine-Gordon equation: Continuous: (1) ϕ t = ϕ α sin ϕ x
9 6 E CELLEDONI ET AL The Sine-Gordon equation is of type (11) with () H = where u := ( ϕ π ) and [ 1 π + 1 (3) S = ( ) ϕ + α (1 cos ϕ)] dx, x ( 1 1 (Note that it follows that π = ϕ t ) Boundary conditions: periodic, u(,t) = u(,t) Semi-discrete: finite differences 3 (4) H fd = [ ] 1 1 π j + ( x) (ϕ j+1 ϕ j ) + α (1 cos ϕ j ) j ) (5) S = ( id id ) The resulting system of ordinary differential equations is [ ] [ ] π (6) = S H ϕ π fd = 1, x Lϕ α sin ϕ where L is the circulant matrix 1 1 L = We have used the bold variables ϕ and π for the finite dimensional vectors [ϕ 1,ϕ,, ϕ N ], et cetera, which replace the functions π and ϕ in the (semi-) discrete case Where necessary, we will write ϕ n, et cetera to denote the vector ϕ at t + n t The integral in the AVF method can be calculated exactly to give 4 (7) 1 t [ ] ϕn+1 ϕ n = π n+1 π n [ (π n+1 + π n )/ L(ϕ n+1 + ϕ n )/ α(cos ϕ n+1 cos ϕ n )/(ϕ n+1 ϕ n ) 3 Summations of the form P j mean P N 1 j= unless stated otherwise 4 For numerical computations, care must be taken to avoid problems when the difference ϕ n+1 ϕ n in the denominator of (7) becomes small We used the sum-to-product identity cos a cos b = sin((a + b)/) sin((a b)/) to give a more numerically amenable expression ]
10 ENERGY CONSERVING PDES 7 Semi-discrete: spectral discretization Instead of using finite differences for the discretization of the spatial derivative in (), one may use a spectral discretization This can be thought of as replacing ϕ with its Fourier series, truncated after N terms, where N is the number of spatial intervals, and differentiating the Fourier series This can be calculated, using the discrete Fourier transform 5 (DFT), as F 1 N d NF N ϕ where F N is the matrix of DFT coefficients with entries given by [F N ] n,k = ωn nk,ω N = e iπ/n Additionally, [F 1 N ] n,k = 1 N ω nk N and d N is a diagonal matrix whose (non-zero) entries are the scaled wave-numbers 6 [d N ] k,k = iπk/l, k = N +1,, N, (for N even), where l = b a is the extent of the spatial domain; that is l/n = x (For more details on properties of the DFT and its application to spectral methods see [] and [18]) (8) H sp = [ 1 π j + 1 [ F 1 N d NF N ϕ ] ] j + α(1 cos ϕ j), j (9) S = ( id id ) The resulting system of ODEs is then given by (1) [ ϕ π ] [ = S H sp = (F 1 N π D N F N ) (F 1 d NF N ϕ) α sin ϕ where [D N ] n,k = θ k Again, the integral in the AVF method can be calculated exactly to give N ], (11) (1) ϕ n+1 ϕ n t π n+1 π n t =(π n+1 + π n )/, = (F 1 N D N F N ) (F 1 N d NF N )(ϕ n+1 + ϕ n )/ α(cos ϕ n+1 cos ϕ n )/(ϕ n+1 ϕ n ) Initial conditions and numerical data for both discretizations: Spatial domain, number N of spatial intervals, and -step size t used were 7 x [, ], N =, t =1, parameter: α =1 5 In practice, one uses the fast Fourier transform algorithm to calculate the DFTs in O(N log N) operations 6 Care must be taken with the ordering of the wave numbers since different computer packages use different effective orderings of the DFT/IDFT matrices in their algorithms Additionally, one must ensure that all modes of the Fourier spectrum are treated symmetrically for N even, this requires replacing the k = N entry with zero For the FFT/IFFT algorithms in Matlab the vector of wave numbers is [,, N 1,, N +1,, 1] 7 Here and below, if x [a, b], then x = b a N, and x j = a + j x, j =, 1,, N
11 8 E CELLEDONI ET AL 1 x 1!4 SG equation, Nspace= 1, t = 1 avf Midpoint 1!1 1! SG equation, Nspace= 1, t = 1 avf Midpoint Energy Error!1!!3!4 Global Error 1!3 1!4 1!5! ! Figure 1: Sine-Gordon equation with finite differences semi-discretization: Energy error (left) and global error (right) vs, for AVF and implicit midpoint integrators Energy Error # '(#!!&!!#!$!% SG equation, Nspace= 1, t= 1 ( )*+$,-/-1 Global Error #!!$ #!!% #!!' SG equation, Nspace= 1, t= 1 )*+$,-/-1 ( #!!"!&!" (! " #! #" $! $" %! #!!& (! " #! #" $! $" %! Figure : Sine-Gordon equation with spectral semi-discretization: Energy error (left) and global error (right) vs, for AVF and implicit midpoint integrators Initial conditions: (13) ϕ(x, ) =, π(x, ) = 8 cosh(x) Right-moving kink and left-moving anti-kink solution Numerical comparisons of the AVF method with the well known (symplectic) implicit midpoint integrator 8 are given in figure 1 for the finite differences discretization, and in figure for the spectral discretization Example Korteweg-de Vries equation: Continuous: (14) u t = 6u u x 3 u x 3, 8 Recall that the implicit midpoint integrator is given by u n+1 u n t = f un+un+1
12 ENERGY CONSERVING PDES 9 x 1!3 KdV equation, Nspace= 41, t = 1 18 KdV equation, Nspace= 41, t = 1 16 Energy Error!!4!6!8 avf Midpoint Global Error avf Midpoint!1 4!1! Figure 3: Korteweg-de Vries equation: Energy error (left) and global error (right) vs, for AVF and implicit midpoint integrators (15) H = [ ] 1 (u x) u 3 dx, (16) S = x Boundary conditions: periodic, u(,t) = u(,t) Semi-discrete: (17) H = [ ] 1 ( x) (u j+1 u j ) u 3 j, j (18) S = 1 x Initial conditions and numerical data: x [, ], N = 4, t =1 Initial condition: u(x, ) = 6 sech (x) (for two solitons) Example 3 Nonlinear Schrödinger equation: Continuous: (19) ( ) ( u i t u = i )( δh δu δh δu ),
13 1 E CELLEDONI ET AL 16 x 1!6 14 NLS equation, Nspace= 1, t = 5 1!1 NLS equation, Nspace= 1, t = 5 avf Midpoint Energy Error avf Midpoint Global Error 1! 1!3 1!4! ! Figure 4: Nonlinear Schrödinger equation: Energy error (left) and global error (right) vs, for AVF and implicit midpoint integrators where u denotes the complex conjugate of u [ () H = u x (1) S = ( i i + γ u 4 ] ) dx, Boundary conditions: periodic, u(,t) = u(,t) Semi-discrete: () H = [ 1 ( x) u j+1 u j + γ ] u j 4, j ( id (3) S = i id Initial conditions and numerical data: x [, ], N =, t =5, parameter: γ =1 Initial conditions: ) (4) { Ru(x, ) = exp ( (x 1) / ), Iu(x, ) = exp ( x / ) Example 4 Nonlinear Wave Equation: Continuous:
14 ENERGY CONSERVING PDES x 1!6 NLS equation, Nspace= 1, t = 5 Probability Error avf Midpt! Figure 5: Nonlinear Schrödinger equation: Total probability error vs, for AVF and implicit midpoint integrators The D wave equation (5) ϕ t V (ϕ) = ϕ, ϕ = ϕ(x, y, t), (x, y) [ 1, 1] [ 1, 1], t, ϕ is a Hamiltonian PDE with Hamiltonian function 1 1 [ ] 1 (6) H = (π + ϕ x + ϕ y)+v (ϕ) dx dy, 1 1 and the operator S is the canonical symplectic matrix Boundary conditions: periodic Semi-discrete: We discretise the Hamiltonian in space with a tensor product Lagrange quadrature formula based on p + 1 Gauss-Lobatto-Legendre (GLL) quadrature nodes in each space direction We obtain (7) ( H = 1 p p p ) ( p ) w j1 w j π j 1,j + d j1,kϕ k,j + d j,mϕ j1,m + 1 ϕ4 j 1,j, j 1= j = where d j1,k = dl k(x) dx k= m=, and l k (x) is the k-th Lagrange basis function based x=xj1 on the GLL quadrature nodes x,, x p, and with w,, w p the corresponding quadrature weights The numerical approximation is p p (8) ϕ p (x, y, t) = ϕ k,m (t)l k (x)l m (y), k= m= and has the property ϕ p (x j1,y j,t)=ϕ j1,j (t), so that the data can be stored in the (p + 1) (p + 1) matrix with entries ϕ j1,j Initial conditions and numerical data: (x, y) [ 1, 1], V(ϕ) = ϕ4 4
15 1 E CELLEDONI ET AL T = T =315 %&!"!' &'!"!(!$ %"!$!"!" $ $ " "!$! "#$ "!"#$!!!!!"#$ " "#$! T =6875 T = 1!$! "#$ "!"#$!!!!!"#$ " "#$! )*!"!+ &'!"!( ( %" '!$ &!" % $ " "!%!!$! "#$ "!"#$!!!!!"#$ " "#$! "#$ "!"#$!!!!!"#$ " "#$! Figure 6: Snapshots of the solution of the D wave equation at different s AVF method with step-size t = 65 Space discretization with 6 Gauss Lobatto nodes in each space direction Numerical solution interpolated on a equidistant grid of 1 nodes in each space direction Initial condition: ϕ(x, y, ) = sech(1x)sech(1y), π(x, y, ) = In figure 6 we show some snapshots of the solution shown in figure 7 Linear conservative PDEs Example 5 (Linear) Time-dependent Schrödinger Equation: Continuous: The energy error is (9) u t = u i iv (x)u x This equation is bi-hamiltonian, ie it has independent symplectic structures The first Hamiltonian formulation has [ π ] (3) H 1 = u x V (x) u dx and π (31) S 1 = ( i i )
16 ENERGY CONSERVING PDES 13!() *+&!!&& "6+789/+/:;8,-51+89<+81=+5=/&)>!!!() 1/34+/5!&!&()!"!"()!'! " # $ % &!,-/ Figure 7: The D wave equation (5) MATLAB routine ode15s with absolute and relative tolerance 1 14 (dashed line), and AVF method with step size t = 1/( 5 ) (solid line) Energy error versus Time interval [, 1] Space discretization with 6 Gauss Lobatto nodes in each space direction The second Hamiltonian formulation has (3) H = and ( (33) S = π π u dx x V (x) x + V (x) ) Boundary conditions: periodic, u( π, t) = u(π, t) Semi-discrete: (34) H 1 = [ 1 ] ( x) u j+1 u j V (x j ) u j, j ( id (35) S 1 = i id The second semi-discretization is ) (36) H = j u j, ( A (37) S = i A ),
17 14 E CELLEDONI ET AL 5 x 1 15 Energy error Figure 8: Linear Schrödinger equation: Error in energy H 1 x vs, AVF method where (38) A = V V V Both discretizations result in the same semi-discrete system and the AVF method (which in the linear case coincides with the midpoint rule) therefore preserves both H 1 and H, as well as the two symplectic structures Initial conditions and numerical data: x [ π, π], N = 5, t =1, V(x) =1 cos(x) Initial conditions: (39) Ru(x, ) = e ( x ), Iu(x, ) = Example 6 Maxwell s Equations (1d): Maxwell equation Continuous: We first look at the one-dimensional (4) [ E t B ] = [ c x c x ][ δh δe δh δb ], (41) H(E, B) = 1 c 1 ( E + B ) dx,
18 ENERGY CONSERVING PDES 15 and (4) S = [ x x S is skew-adjoint on {(E, B) C 1 : E() = E(1) = } (and therefore on the Sobolev space H 1 ) Boundary conditions: { E(,t)=E(1,t)=, (43) ] B B x (,t)= x (1,t)= Semi-discrete: We now obtain H by discretizing H in a simple way by applying the trapezoidal rule to the integral (41) at the points x j = 1 N j and dividing by x, that is (44) N 1 ( H(E 1, E N 1,B,,B N )= c 1 ) E j + c 1 N 1 ( 4 B + c 1 ) B j + c 1 4 B N, j=1 j=1 where we have already used that E(x,t) = E(x N,t) = operator S is discretized with central differences yielding [ ] N 1,N 1 G (45) S = G T N+1,N+1, The differential where the (N 1) (N + 1) matrix G is given by (46) G = 1 x Initial conditions and numerical data: Note that the Neumann boundary conditions are satisfied at least to order 1 in space, despite the fact that we only intended to somehow approximate the energy H The numerical experiments confirm that the discrete energy H x is preserved to machine precision Figure 9 shows the error of the AVF method for the Maxwell equation with N = 1, t =1, c = 1, and initial value E(x, ) = e 1(x 1 ), B(x, ) = e 1(x 1 ) Example 7 Maxwell s Equations (3D): Continuous: We consider Maxwell s equations in CGS units for the electromagnetic field in a vacuum [ ] [ ][ ] B c B (47) = t E c E
19 16 E CELLEDONI ET AL 6 x Energy Error Figure 9: One-dimensional Maxwell equation: energy error vs, AVF integrator with the operator (48) := z y z x y x This equation has two Hamiltonian formulations of type (11) The first Hamiltonian formulation has the helicity Hamiltonian ( (49) H 1 = c 1 BT ( B)+c 1 ) ET ( E) dxdydz Q (cf [1]) and the operator (5) S 1 = [ I3 I 3 ], where I 3 designates the 3 3 unit matrix The second Hamiltonian formulation has the Hamiltonian ( (51) H = c 1 BT B + c 1 ) ET E dxdydz Q (cf [7]) and the operator (5) S = [ ] Boundary condition: periodic on the unit cube Q
20 ENERGY CONSERVING PDES 17 Semi-discrete: On a regular grid with lexicographical ordering in every component of E (resp B) and concatenating the discretized components gives one discrete vector E h (resp B h ), the operator is represented by a matrix A The discretisation in the first case is given by (53) H 1 = c 1 ET h AE h + c 1 BT h AB h and the obvious discretisation of S 1 For the quadratic Hamiltonian, (54) H = c 1 ET h E h + c 1 BT h B h and (55) S = [ A A ] Both discretisations result in the same semi-discrete system [ ] [ ][ ] Ḃh A Bh (56) = A E h Ėh and the average vector field method preserves both H 1 and H Initial conditions and numerical data: Preservation of H 1 and H is numerically confirmed by an experiment with random initial data on a regular grid with 3 points in every direction and constant c = 1 The result of the AVF method with t =1 can be seen in Figure 1 3 Dissipative PDEs 31 Nonlinear dissipative PDEs Example 31 Allen - Cahn equation: Continuous: (31) u t = du xx + u u 3, (3) H = [ 1 d (u x) 1 u + 1 ] 4 u4 dx, (33) N = 1 Boundary conditions: Neumann, u x (,t)=,u x (1,t) =
21 18 E CELLEDONI ET AL 5 15 Energies Figure 1: Three-dimensional Maxwell equation: plots of energies H 1 x (dash-dot) and H x (dash) vs 1!7 AC equation, Nspace= 11, t = 1 avf!15 AC equation, Nspace= 11, t = 1 avf!16 1!8!17 Global Error 1!9 1!1 Energy!18!19! 1!11!1! 1! ! Figure 31: Allen-Cahn equation: Global error (left) and energy (right) vs, AVF integrator Semi-discrete: (34) H = N [ j= d ( x) (u j+1 u j ) 1 u j u4 j ], (35) N = id Initial conditions and numerical data: x [, 1], N = 1, t =1, parameter: d =1 Initial condition: u(x, ) = cos(πx)
22 ENERGY CONSERVING PDES 19 1!6 CH equation, Nspace= 51, t = avf CH equation, Nspace= 51, t = avf! Global Error 1!7 1!8 Energy!4!6!8!1!1 1! ! Figure 3: Cahn-Hilliard equation: Global error (left) and energy (right) vs, AVF integrator Example 3 Cahn-Hilliard equation: Continuous: (36) (37) H = u t = x ( pu + ru 3 + qu xx ), [ 1 pu ru4 1 ] q (u x) dx, (38) N = x Boundary condition: periodic, u(,t) = u(1,t) Semi-discrete: (39) H = j [ 1 pu j ru4 j 1 ] q ( x) (u j+1 u j ), (31) N = 1 ( x) Initial conditions and numerical data: x [, 1], N = 5, t =1/1, parameters: p = 1, q= 1, r=1 Initial condition: u(x, ) = 1 sin(πx) + 1 cos(4πx) + 6 sin(4πx) + cos(1πx)
23 E CELLEDONI ET AL Example 33 Ginzburg-Landau equation: Continuous: A Ginzburg-Landau equation arising in a model of traffic flow is given by (311) u t = ( x ɛ x )[ 6u + x u u 3] and is a slight modification of the model considered in [14] and [15] The equation can be written as (31) with u t = N δh δu, (313) N = x ɛ x and (314) H = [ 3u 1 ( ) ] u 1 x 4 u4 dx Boundary condition: u(±5,t) = and u xx (±5,t) = Semi-discrete: (315) H = and N 1 j=1 [ 3u j 1 ( uj+1 u j x (316) N = A ɛb, where (317) A = 1 x ) ] 1 4 u4 j, u N = is a discretisation of x, and (318) B = 1 ( x)
24 ENERGY CONSERVING PDES 1 Lyapunov function Lyapunov function Figure 33: Ginzburg-Landau equation: Plots of the energy function H x computed with the AVF method (left) and Matlab s ode45 (right) Note that ode45 does not exhibit the correct monotonic decrease in energy is a discretisation of x The average vector field method preserves the decay of function H in contrast to some standard integrators Initial conditions and numerical data: x [ 5, 5], N = 1, t =1, parameter: ɛ = 1 6 Initial condition: u(x, ) = e 1(x 1 ) In Figure 33, we compare the AVF method with Matlab s ode45, the latter not preserving the monotonic decay of H 3 Linear dissipative PDEs Example 34 Heat equation: Continuous: The heat equation (319) u t = u xx, is a dissipative PDE and can be written in the form (11), ie (3) u t = N δh 1 1 δu, u t = N δh δu, with the Lyapunov functions H 1 (u) = 1 1 u x dx and H (u) = 1 1 u dx and the operators N 1 = 1 and N = x, respectively Boundary conditions: u(,t) = u(1,t) = Semi-discrete: N 1 1 (31) H 1 = u ( x) 1 + (u j u j 1 ) + u N 1 j=
25 E CELLEDONI ET AL First Lyapunov function Second Lyapunov function Figure 34: Heat equation: plots of Lyapunov functions H 1 x (left) and H x (right) vs, AVF integrator and (3) H = N 1 j=1 1 (u j), as well as (33) N = 1 ( x) and the obvious discretisation of N 1 With these choices, both discretisations yield identical semi-discrete equations of motion and therefore H 1 and H are simultaneously Lyapunov functions of the semi-discrete system and therefore, the AVF integrator preserves both Lyapunov functions Initial conditions and numerical data: (34) x [, 1], N = 5, t =5 Initial condition: u(x, ) = x(1 x) This system is numerically illustrated in Figure 34, where the monotonic decrease of the Lyapunov functions for the heat equation in (31) and (3) is shown 4 Concluding Remarks The concept of energy, ie its preservation or dissipation, has far reaching consequences in the physical sciences Therefore many methods to preserve energy, and several to preserve the correct dissipation of energy (eg [6, 11]), have
26 ENERGY CONSERVING PDES 3 been proposed for ordinary differential equations Surprisingly, when partial differential equations are considered, a unified way to discuss the preservation or correct dissipation of energy is missing and similar ideas are often developed from scratch (eg [5, 1]) In this paper, we have presented a systematic and unified way to discretise partial differential equations and to preserve their correct energy preservation, or dissipation, by the average vector-field method For the equations treated in this paper, one can replace the average vector-field method by any energy-preserving B-series method, while retaining the advantageous properties of energy preservation or dissipation More generally, geometric integrators for Hamiltonian or non-hamiltonian PDEs with non-constant matrix D can be constructed using discrete gradient methods, cf [11] Acknowledgements This research is supported by the Australian Research Council We are grateful to Will Wright for many useful discussions REFERENCES 1 N Anderson and AM Arthurs, Helicity and variational principles for Maxwell s equations, Int J Electronics, 54 (1983), pp WL Briggs and VE Henson, The DFT: An Owner s Manual for the Discrete Fourier Transform, SIAM, E Celledoni, RI McLachlan, DI McLaren, B Owren, GRW Quispel, and WM Wright, Energy-preserving Runge-Kutta methods, Preprint, submitted, 8 4 RP Feynman Conservation of energy Chapter 4, Volume 1, The Feynman Lectures on Physics, Addison-Wesley Pub Co, D Furihata, Finite Difference Schemes for u t = ` x α δg δu that Inherit Energy Conservation or Dissipation Property, J Comput Phys, 156 (1999), pp V Grimm and GRW Quispel Geometric integration methods that preserve Lyapunov functions, BIT, 45 (5), pp T Hirono, W Lui, S Seki and Y Yoshikuni, A three-dimensional fourth-order Finite-Difference Time-Domain scheme using a symplectic integrator propagator, IEEE Trans Microwave Theory and Tech, 49 (1), pp A Iserles, A First Course in the Numerical Analysis of Differential Equations, Cambridge University Press, S Klainerman, Partial differential equations in The Princeton Companion to Mathematics, T Gowers, Ed Princeton University Press, 8, pp T Matsuo New conservative schemes with discrete variational derivatives for nonlinear wave equations, J Comput Appl Math, 3 (7), pp RI McLachlan, GRW Quispel, and N Robidoux Geometric integration using discrete gradients, Phil Trans Roy Soc A, 357 (1999), pp RI McLachlan and N Robidoux Antisymmetry, pseudospectral methods, weighted residual discretizations, and energy conserving partial differential equations, Preprint, 13 RI McLachlan and N Robidoux Antisymmetry, pseudospectral methods, and conservative PDEs, in International Conference on Differential Equations, Vol 1, (Berlin, 1999), World Sci Publ, River Edge, NJ,, pp
27 4 E CELLEDONI ET AL 14 Eq(7) in: T Nagatani Thermodynamic theory for the jamming transition in traffic flow, Phys Rev E, 58 (1998), pp T Nagatani Time-dependent Ginzburg Landau equation for the jamming transition in traffic flow, Physica A, 58 (1998), pp P Olver, Applications of Lie Groups to Differential Equations, Second Edition, Springer-Verlag, New York, GRW Quispel and DI McLaren A new class of energy-preserving numerical integration methods, J Phys A: Math Theor, 41 (8), 456 (7pp) 18 LN Trefethen, Spectral Methods in MATLAB, SIAM, 5 APPENDIX PROOF OF LEMMA 11 We denote the consistent discretisation of (51) H[u] = H(x; u (n) ) dx by Ω (5) H [u] = Ω a n H(x n ; u (n+k) ) x Here H = H x, and u (n+k) denotes the set (53) { u1 (n + k 1,1 ),, u 1 (n + k 1,p ); u (n + k,1 ),, u (n + k,p ); ; u m (n + k m,1 ),, u m (n + k m,p ) }, where n and the k i,j are discrete vectors in Z d, and Ω is some discrete approximation to Ω The discrete analogue of the variational derivative δh is then defined by δu d (δh ) (54) v n x dε H (u + εv) ε= =: Ω = (v δh δu δu ) n x
28 ENERGY CONSERVING PDES 5 We now calculate the lhs of eq(54): (55) d dε H (u + εv) ε= = d a n H(x n ; u 1 (n + k 1,1 )+εv 1 (n + k 1,1 ), dε Ω, u m (n + k m,p )+εv m (n + k m,p )) ε= = Ω [ an k1,1h (x n k1,1; u 1 (n),, u m (n + k m,p k 1,1 )) +a n km,p H mp+1 (x n km,p ; u 1 (n + k 1,1 k m,p ),, u m (n)) ] v n x = [ S kq,r H(x n, u 1 (n + k 1,1 ),, u m (n + k m,p )) ] v n x u n Ω q,r where S a denotes the shift over the multi-index a It follows that (56) d dε H (u + εv) ε= =(v H ) Since eqs (54) and (56) must hold for all vectors v, we have (57) and hence (58) δh δu x H, δh δu H While the proof above is for a scalar function u, the case where u is vectorvalued follows in a similar fashion
Energy Preserving Numerical Integration Methods
Energy Preserving Numerical Integration Methods Department of Mathematics and Statistics, La Trobe University Supported by Australian Research Council Professorial Fellowship Collaborators: Celledoni (Trondheim)
More informationarxiv: v1 [math.na] 21 Feb 2012
Preserving energy resp. dissipation in numerical PDEs using the Average Vector Field method E. Celledoni a, V. Grimm b, R. I. McLachlan c, D. I. McLaren d, D. O Neale d,e, B. Owren a, G. R. W. Quispel
More informationENERGY PRESERVATION (DISSIPATION) IN NUMERICAL ODEs and PDEs
ENERGY PRESERVATION (DISSIPATION) IN NUMERICAL ODEs and PDEs E. CELLEDONI Department of Mathematical Sciences, NTNU 749 Trondheim, Norway. email: elenac@math.ntnu.no Abstract This note considers energy-preservation
More informationOn the construction of discrete gradients
On the construction of discrete gradients Elizabeth L. Mansfield School of Mathematics, Statistics and Actuarial Science, University of Kent, Canterbury, CT2 7NF, U.K. G. Reinout W. Quispel Department
More informationENERGY PRESERVING INTEGRATION OF BI-HAMILTONIAN PARTIAL DIFFERENTIAL EQUATIONS
TWMS J. App. Eng. Math. V., N.,, pp. 75-86 ENERGY PRESERVING INTEGRATION OF BI-HAMILTONIAN PARTIAL DIFFERENTIAL EQUATIONS B. KARASÖZEN, G. ŞİMŞEK Abstract. The energy preserving average vector field (AVF
More informationENERGY PRESERVING METHODS FOR VOLTERRA LATTICE EQUATION
TWMS J. App. Eng. Math. V., N.,, pp. 9- ENERGY PRESERVING METHODS FOR VOLTERRA LATTICE EQUATION B. KARASÖZEN, Ö. ERDEM, Abstract. We investigate linear energy preserving methods for the Volterra lattice
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationSOLVING ODE s NUMERICALLY WHILE PRESERVING ALL FIRST INTEGRALS
SOLVING ODE s NUMERICALLY WHILE PRESERVING ALL FIRST INTEGRALS G.R.W. QUISPEL 1,2 and H.W. CAPEL 3 Abstract. Using Discrete Gradient Methods (Quispel & Turner, J. Phys. A29 (1996) L341-L349) we construct
More informationSpectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions
Applied Mathematical Sciences, Vol. 1, 2007, no. 5, 211-218 Spectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions M. Javidi a and A. Golbabai b a Department
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationWRT in 2D: Poisson Example
WRT in 2D: Poisson Example Consider 2 u f on [, L x [, L y with u. WRT: For all v X N, find u X N a(v, u) such that v u dv v f dv. Follows from strong form plus integration by parts: ( ) 2 u v + 2 u dx
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationGeneralised Summation-by-Parts Operators and Variable Coefficients
Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation
More informationDiscrete Variational Derivative Method
.... Discrete Variational Derivative Method one of structure preserving methods for PDEs Daisue Furihata Osaa Univ. 211.7.12 furihata@cmc.osaa-u.ac.jp (Osaa Univ.)Discrete Variational Derivative Method
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 informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationGEOMETRIC INTEGRATION METHODS THAT PRESERVE LYAPUNOV FUNCTIONS
BIT Numerical Mathematics 25, Vol. 45, No., pp. 79 723 c Springer GEOMETRIC INTEGRATION METHODS THAT PRESERVE LYAPUNOV FUNCTIONS V. GRIMM and G.R.W. QUISPEL Department of Mathematics, and Centre for Mathematics
More informationEnergy-preserving Pseudo-spectral Methods for Klein-Gordon-Schrödinger Equations
.. Energy-preserving Pseudo-spectral Methods for Klein-Gordon-Schrödinger Equations Jingjing Zhang Henan Polytechnic University December 9, 211 Outline.1 Background.2 Aim, Methodology And Motivation.3
More informationSymmetry Reductions of (2+1) dimensional Equal Width. Wave Equation
Authors: Symmetry Reductions of (2+1) dimensional Equal Width 1. Dr. S. Padmasekaran Wave Equation Asst. Professor, Department of Mathematics Periyar University, Salem 2. M.G. RANI Periyar University,
More informationSpatial discretization of partial differential equations with integrals
Spatial discretization of partial differential equations with integrals Robert I. McLachlan IFS, Massey University, Palmerston North, New Zealand E-mail: r.mclachlan@massey.ac.nz We consider the problem
More informationOn universality of critical behaviour in Hamiltonian PDEs
Riemann - Hilbert Problems, Integrability and Asymptotics Trieste, September 23, 2005 On universality of critical behaviour in Hamiltonian PDEs Boris DUBROVIN SISSA (Trieste) 1 Main subject: Hamiltonian
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
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 informationDiagonalization of the Coupled-Mode System.
Diagonalization of the Coupled-Mode System. Marina Chugunova joint work with Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaborators: Mason A. Porter, California Institute
More informationSufficient conditions for functions to form Riesz bases in L 2 and applications to nonlinear boundary-value problems
Electronic Journal of Differential Equations, Vol. 200(200), No. 74, pp. 0. ISSN: 072-669. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Sufficient conditions
More informationEnergy-Preserving Runge-Kutta methods
Energy-Preserving Runge-Kutta methods Fasma Diele, Brigida Pace Istituto per le Applicazioni del Calcolo M. Picone, CNR, Via Amendola 122, 70126 Bari, Italy f.diele@ba.iac.cnr.it b.pace@ba.iac.cnr.it SDS2010,
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 informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationGalerkin method for the numerical solution of the RLW equation using quintic B-splines
Journal of Computational and Applied Mathematics 19 (26) 532 547 www.elsevier.com/locate/cam Galerkin method for the numerical solution of the RLW equation using quintic B-splines İdris Dağ a,, Bülent
More informationDispersion relations, stability and linearization
Dispersion relations, stability and linearization 1 Dispersion relations Suppose that u(x, t) is a function with domain { < x 0}, and it satisfies a linear, constant coefficient partial differential
More informationNumerical simulation for a nonlinear partial differential equation with variable coefficients by means of the discrete variational derivative method
Journal of Computational and Applied Mathematics 94 (6) 45 459 www.elsevier.com/locate/cam Numerical simulation for a nonlinear partial differential equation with variable coefficients by means of the
More informationOrthogonal Symmetric Toeplitz Matrices
Orthogonal Symmetric Toeplitz Matrices Albrecht Böttcher In Memory of Georgii Litvinchuk (1931-2006 Abstract We show that the number of orthogonal and symmetric Toeplitz matrices of a given order is finite
More informationNumerical Methods I Orthogonal Polynomials
Numerical Methods I Orthogonal Polynomials Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Nov. 4th and 11th, 2010 A. Donev (Courant Institute)
More informationTHE 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 informationThe L p -dissipativity of first order partial differential operators
The L p -dissipativity of first order partial differential operators A. Cialdea V. Maz ya n Memory of Vladimir. Smirnov Abstract. We find necessary and sufficient conditions for the L p -dissipativity
More informationCONSTRUCTING OF CONSTRAINT PRESERVING SCHEME FOR EINSTEIN EQUATIONS
CONSTRUCTING OF CONSTRAINT PRESERVING SCHEME FOR EINSTEIN EQUATIONS TAKUYA TSUCHIYA AND GEN YONEDA arxiv:1610.04543v3 [gr-qc] 14 Jul 2017 Abstract. We propose a new numerical scheme of evolution for the
More informationNotes on the Inverse Scattering Transform and Solitons. November 28, 2005 (check for updates/corrections!)
Notes on the Inverse Scattering Transform and Solitons Math 418 November 28, 2005 (check for updates/corrections!) Among the nonlinear wave equations are very special ones called integrable equations.
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationarxiv:nlin/ v1 [nlin.si] 25 Sep 2006
Remarks on the conserved densities of the Camassa-Holm equation Amitava Choudhuri 1, B. Talukdar 1a and S. Ghosh 1 Department of Physics, Visva-Bharati University, Santiniketan 73135, India Patha Bhavana,
More informationNumerical Analysis Preliminary Exam 10.00am 1.00pm, January 19, 2018
Numerical Analysis Preliminary Exam 0.00am.00pm, January 9, 208 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationThe elliptic sinh-gordon equation in the half plane
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 8 25), 63 73 Research Article The elliptic sinh-gordon equation in the half plane Guenbo Hwang Department of Mathematics, Daegu University, Gyeongsan
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationMultisoliton Interaction of Perturbed Manakov System: Effects of External Potentials
Multisoliton Interaction of Perturbed Manakov System: Effects of External Potentials Michail D. Todorov Faculty of Applied Mathematics and Computer Science Technical University of Sofia, Bulgaria (Work
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationINTRODUCTION TO COMPUTER METHODS FOR O.D.E.
INTRODUCTION TO COMPUTER METHODS FOR O.D.E. 0. Introduction. The goal of this handout is to introduce some of the ideas behind the basic computer algorithms to approximate solutions to differential equations.
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationSOME PROPERTIES OF SYMPLECTIC RUNGE-KUTTA METHODS
SOME PROPERTIES OF SYMPLECTIC RUNGE-KUTTA METHODS ERNST HAIRER AND PIERRE LEONE Abstract. We prove that to every rational function R(z) satisfying R( z)r(z) = 1, there exists a symplectic Runge-Kutta method
More informationNumerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method
Numerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method Y. Azari Keywords: Local RBF-based finite difference (LRBF-FD), Global RBF collocation, sine-gordon
More informationImplicit-explicit exponential integrators
Implicit-explicit exponential integrators Bawfeh Kingsley Kometa joint work with Elena Celledoni MaGIC 2011 Finse, March 1-4 1 Introduction Motivation 2 semi-lagrangian Runge-Kutta exponential integrators
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationA SIXTH ORDER AVERAGED VECTOR FIELD METHOD *
Journal of Computational Mathematics Vol.34, No.5, 6, 479 498. http://www.global-sci.org/jcm doi:.48/jcm.6-m5-65 A SIXTH ORDER AVERAGED VECTOR FIELD METHOD * Haochen Li Jiangsu Key Laboratory for NSLSCS,
More informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationExact solutions through symmetry reductions for a new integrable equation
Exact solutions through symmetry reductions for a new integrable equation MARIA LUZ GANDARIAS University of Cádiz Department of Mathematics PO.BOX, 1151 Puerto Real, Cádiz SPAIN marialuz.gandarias@uca.es
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More information1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationEQUIVALENCE OF STURM-LIOUVILLE PROBLEM WITH FINITELY MANY ABDULLAH KABLAN, MEHMET AKİF ÇETİN
International Journal of Analysis and Applications Volume 16 Number 1 (2018) 25-37 URL: https://doi.org/10.28924/2291-8639 DOI: 10.28924/2291-8639-16-2018-25 EQUIVALENCE OF STURM-LIOUVILLE PROBLEM WITH
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 informationHamiltonian partial differential equations and Painlevé transcendents
The 6th TIMS-OCAMI-WASEDA Joint International Workshop on Integrable Systems and Mathematical Physics March 22-26, 2014 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN
More informationON THE HÖLDER CONTINUITY OF MATRIX FUNCTIONS FOR NORMAL MATRICES
Volume 10 (2009), Issue 4, Article 91, 5 pp. ON THE HÖLDER CONTINUITY O MATRIX UNCTIONS OR NORMAL MATRICES THOMAS P. WIHLER MATHEMATICS INSTITUTE UNIVERSITY O BERN SIDLERSTRASSE 5, CH-3012 BERN SWITZERLAND.
More informationA fast and well-conditioned spectral method: The US method
A fast and well-conditioned spectral method: The US method Alex Townsend University of Oxford Leslie Fox Prize, 24th of June 23 Partially based on: S. Olver & T., A fast and well-conditioned spectral method,
More informationA Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional Hyperbolic Telegraph Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(01) No.3,pp.59-66 A Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
More informationITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS
ITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS LONG CHEN In this chapter we discuss iterative methods for solving the finite element discretization of semi-linear elliptic equations of the form: find
More informationProlongation structure for nonlinear integrable couplings of a KdV soliton hierarchy
Prolongation structure for nonlinear integrable couplings of a KdV soliton hierarchy Yu Fa-Jun School of Mathematics and Systematic Sciences, Shenyang Normal University, Shenyang 110034, China Received
More informationSymplectic integration with Runge-Kutta methods, AARMS summer school 2015
Symplectic integration with Runge-Kutta methods, AARMS summer school 2015 Elena Celledoni July 13, 2015 1 Hamiltonian systems and their properties We consider a Hamiltonian system in the form ẏ = J H(y)
More informationExponential integrators for semilinear parabolic problems
Exponential integrators for semilinear parabolic problems Marlis Hochbruck Heinrich-Heine University Düsseldorf Germany Innsbruck, October 2004 p. Outline Exponential integrators general class of methods
More informationSome recent results on controllability of coupled parabolic systems: Towards a Kalman condition
Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition F. Ammar Khodja Clermont-Ferrand, June 2011 GOAL: 1 Show the important differences between scalar and non
More informationSpecial classical solutions: Solitons
Special classical solutions: Solitons by Suresh Govindarajan, Department of Physics, IIT Madras September 25, 2014 The Lagrangian density for a single scalar field is given by L = 1 2 µφ µ φ Uφ), 1) where
More informationEnergy Stable Model Order Reduction for the Allen-Cahn Equation
Energy Stable Model Order Reduction for the Allen-Cahn Equation Bülent Karasözen joint work with Murat Uzunca & Tugba Küçükseyhan Institute of Applied Mathematics & Department of Mathematics, METU, Ankara
More informationLinearized methods for ordinary di erential equations
Applied Mathematics and Computation 104 (1999) 109±19 www.elsevier.nl/locate/amc Linearized methods for ordinary di erential equations J.I. Ramos 1 Departamento de Lenguajes y Ciencias de la Computacion,
More informationFree Vibration Analysis of Kirchoff Plates with Damaged Boundaries by the Chebyshev Collocation Method. Eric A. Butcher and Ma en Sari
Free Vibration Analysis of Kirchoff Plates with Damaged Boundaries by the Chebyshev Collocation Method Eric A. Butcher and Ma en Sari Department of Mechanical and Aerospace Engineering, New Mexico State
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
More informationExact Solutions of The Regularized Long-Wave Equation: The Hirota Direct Method Approach to Partially Integrable Equations
Thai Journal of Mathematics Volume 5(2007) Number 2 : 273 279 www.math.science.cmu.ac.th/thaijournal Exact Solutions of The Regularized Long-Wave Equation: The Hirota Direct Method Approach to Partially
More informationThe Spectral-Element Method: Introduction
The Spectral-Element Method: Introduction Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Computational Seismology 1 / 59 Outline 1 Introduction 2 Lagrange
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
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 informationA Spectral Time-Domain Method for Computational Electrodynamics
A Spectral Time-Domain Method for Computational Electrodynamics James V. Lambers Abstract Block Krylov subspace spectral (KSS) methods have previously been applied to the variable-coefficient heat equation
More informationKrylov Subspace Methods for the Evaluation of Matrix Functions. Applications and Algorithms
Krylov Subspace Methods for the Evaluation of Matrix Functions. Applications and Algorithms 2. First Results and Algorithms Michael Eiermann Institut für Numerische Mathematik und Optimierung Technische
More informationDmitry Pelinovsky Department of Mathematics, McMaster University, Canada
Spectrum of the linearized NLS problem Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaboration: Marina Chugunova (McMaster, Canada) Scipio Cuccagna (Modena and Reggio Emilia,
More informationLinear and Multilinear Algebra. Linear maps preserving rank of tensor products of matrices
Linear maps preserving rank of tensor products of matrices Journal: Manuscript ID: GLMA-0-0 Manuscript Type: Original Article Date Submitted by the Author: -Aug-0 Complete List of Authors: Zheng, Baodong;
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationBreaking soliton equations and negative-order breaking soliton equations of typical and higher orders
Pramana J. Phys. (2016) 87: 68 DOI 10.1007/s12043-016-1273-z c Indian Academy of Sciences Breaking soliton equations and negative-order breaking soliton equations of typical and higher orders ABDUL-MAJID
More informationAbout Integrable Non-Abelian Laurent ODEs
About Integrable Non-Abelian Laurent ODEs T. Wolf, Brock University September 12, 2013 Outline Non-commutative ODEs First Integrals and Lax Pairs Symmetries Pre-Hamiltonian Operators Recursion Operators
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationNew Diagonal-Norm Summation-by-Parts Operators for the First Derivative with Increased Order of Accuracy
AIAA Aviation -6 June 5, Dallas, TX nd AIAA Computational Fluid Dynamics Conference AIAA 5-94 New Diagonal-Norm Summation-by-Parts Operators for the First Derivative with Increased Order of Accuracy David
More informationON MATRIX VALUED SQUARE INTEGRABLE POSITIVE DEFINITE FUNCTIONS
1 2 3 ON MATRIX VALUED SQUARE INTERABLE POSITIVE DEFINITE FUNCTIONS HONYU HE Abstract. In this paper, we study matrix valued positive definite functions on a unimodular group. We generalize two important
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationCOMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE
BIT 6-85//41-11 $16., Vol. 4, No. 1, pp. 11 118 c Swets & Zeitlinger COMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE WALTER GAUTSCHI Department of Computer Sciences,
More informationNumerical Simulations on Two Nonlinear Biharmonic Evolution Equations
Numerical Simulations on Two Nonlinear Biharmonic Evolution Equations Ming-Jun Lai, Chun Liu, and Paul Wenston Abstract We numerically simulate the following two nonlinear evolution equations with a fourth
More informationMATH 425, HOMEWORK 5, SOLUTIONS
MATH 425, HOMEWORK 5, SOLUTIONS Exercise (Uniqueness for the heat equation on R) Suppose that the functions u, u 2 : R x R t R solve: t u k 2 xu = 0, x R, t > 0 u (x, 0) = φ(x), x R and t u 2 k 2 xu 2
More informationA MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE
Pacific Journal of Applied Mathematics Volume 1, Number 2, pp. 69 75 ISSN PJAM c 2008 Nova Science Publishers, Inc. A MULTI-COMPONENT LAX INTEGRABLE HIERARCHY WITH HAMILTONIAN STRUCTURE Wen-Xiu Ma Department
More information