Artificial boundary conditions for dispersive equations. Christophe Besse
|
|
- Sandra Lester
- 5 years ago
- Views:
Transcription
1 Artificial boundary conditions for dispersive equations by Christophe Besse Institut Mathématique de Toulouse, Université Toulouse 3, CNRS Groupe de travail MathOcéan Bordeaux INSTITUT de MATHEMATIQUES de TOULOUSE
2 Outline of the talk 1 Why Artifical Boundary Conditions? 2 Derivation of some TBC 3 Numerical approximation
3 Introduction Typical equations The Schrödinger Eq. in R d i tψ + ψ + V (x, t, ψ) ψ = 0, (x, t) R d [0; T ] (S) lim ψ(x, t) = 0, t [0; T ] x + ψ(x, 0) = ψ 0(x), x R d ψ(x, t): wave function, complex V can be real potential, V = V (x, t) C (R d R +, R) nonlinear functional, V = f(ψ) Ex: V = q ψ 2 a combination: V = V (x) + f(ψ) ψ 0 compact support in Ω
4 Introduction Typical equations The Korteweg-de-Vries Eq. in R tu 6u xu + xu 3 = 0, (x, t) R [0; T ] (K) lim u(x, t) = 0, t [0; T ] x + u(x, 0) = u 0(x), x R u(x, t): real function u 0 compact support in Ω
5 Domain trunction Problem : Mesh an unbounded domain (here in 1D) R d [0; T ]
6 Domain trunction Problem : Mesh an unbounded domain (here in 1D) R d [0; T ] Truncation R [0; T ] Ω T := ]x l, x r[ [0; T ] Introduction of a fictitious boundary Σ T := Σ [0; T ] with Σ := Ω = {x l, x r} BC on the boundaryσ T
7 Domain trunction Problem : Mesh an unbounded domain (here in 1D) R d [0; T ] Truncation R [0; T ] Ω T := ]x l, x r[ [0; T ] Introduction of a fictitious boundary Σ T := Σ [0; T ] with Σ := Ω = {x l, x r} BC on the boundaryσ T Expression of the boundary condtion with the help of the Dirichlet-to-Neumann map: nψ + Λ + ψ = 0, on Σ T.
8 Domain trunction Problem : Mesh an unbounded domain (here in 1D) R d [0; T ] Truncation R [0; T ] Ω T := ]x l, x r[ [0; T ] Introduction of a fictitious boundary Σ T := Σ [0; T ] with Σ := Ω = {x l, x r} BC on the boundaryσ T Expression of the boundary condtion with the help of the Dirichlet-to-Neumann map: nψ + Λ + ψ = 0, on Σ T.
9 What happen if we do not take BC into account Airy Equation tu + xu 3 = 0, on ( 15, 15) [0; T ], u(x, 0) = exp( x 2 ), x ( 15, 15), xu( 15, t) = 0, t > 0, u(15, t) = 0, xu(15, t) = 0, t > 0.
10 Outline of the talk 1 Why Artifical Boundary Conditions? 2 Derivation of some TBC 3 Numerical approximation
11 TBC for Wave Eq. Homogeneous Wave Equation in 1D t 2 ψ xψ 2 = 0, (x, t) R x [0; T ], (P) ψ(x, 0) = ψ 0(x), x R x, tψ(x, 0) = ψ 1(x), x R x. Hypothesis : supp(ψ 0,1) B(0, R). Solution Simply, with ψ(x, t) = 1 2 (ψ0(x + t) + ψ0(x t)) x+t x t ψ(x, t) = ϕ 1(x + t) + ϕ 2(x t), x 2ϕ 1(x) = ψ 0(x) + 2ϕ 2(x) = ψ 0(x) x ψ 1(y)dy, ψ 1(y)dy. ψ 1(y)dy.
12 TBC for Wave Eq. t x L R R L supp(ϕ 1(x + t)) supp(ϕ 2(x t)) t > 0 and x > R, the supports are disjoint. x = L x = L ψ(x, t) = ϕ 1(x + t) ψ(x, t) = ϕ 2(x t) xψ = tψ xψ = tψ Then, on x = ±L, the Neumann datum is expressed in function of the Dirichlet one nψ + tψ = 0 where n denotes the outwardly unit normal vector to Ω = ( L, L).
13 TBC for Wave Eq. The problem (P) is thus transformed in (P app) (P app) where Γ = { L, L}. t 2 ψ a xψ 2 a = 0, (x, t) Ω [0; T ], ψ a (x, 0) = ψ 0(x), x Ω, tψ a (x, 0) = ψ 1(x), x Ω, nψ a + tψ a = 0, (x, t) Γ [0; T ], ψ Ω (x, t) = ψ a (x, t), (x, t) Ω [0; T ]. Ths BCs do not perturb the solution: we call them Remark: Transparent Boundary Conditions (TBC) 2 t 2 x = ( t x)( t + x).
14 TBC for Wave Eq. Other way of derivation: Laplace transform L (u)(x, ω) = û(x, ω) = with frequency ω = σ + iτ, σ > 0 0 u(x, t)e ωt dt L ( tu)(x, ω) = ωû(x, ω) u(x, 0), L ( t 2 u)(x, ω) = ω 2 û(x, ω) ωu(x, 0) tu(x, 0). 1 Transmission problem Transmission conditions: continuity of traces (u and xu) Interior Problem ( t 2 2 x )v = 0, x Ω, t > 0, xv = xw, x Γ, t > 0, v(x, 0) = ψ 0 (x), x Ω, tv(x, 0) = ψ 1 (x), x Ω, Exterior problem ( t x)w 2 = 0, x Ω, t > 0, w(x, t) = v(x, t), x = ±L, t > 0, w(x, 0) = 0, x Ω, tw(x, 0) = 0, x Ω.
15 TBC for Wave Eq. 2 Laplace tranform w.r.t time for x > L : ω 2 ŵ 2 xŵ = 0. Solution: ŵ(x, ω) = A + (ω)e ω x + A (ω)e ω x 3 How to select the outgoing wave: solution of finite energy A + = 0. ŵ(x, ω) = e ω (x L) L (w(l, ))(ω). Taking derivative + continuity of trace xŵ(x, ω) x=l = ω ˆv(x, ω) x=l 4 Inverse Laplace tranform xw(x, t) x=l = tv(x, t) x=l. 5 Thus, we obtain the TBC for v (continuity of trace) nv + tv = 0
16 The Airy equation tu + xu 3 = 0, on R [0; T ]. u(x, 0) = u 0(x), x R. 1 Splitting between interior and exterior problems left exterior problem right exterior problem interior problem t=0
17 The Airy equation tu + xu 3 = 0, on R [0; T ]. u(x, 0) = u 0(x), x R. 1 Splitting between interior and exterior problems left exterior problem right exterior problem interior problem t=1 The B.Cs are different at left and right boundaries
18 The Airy equation 2 Transmission problem Interior problem ( t + x)u 3 = 0, x l < x < x r, t > 0, u(x, 0) = u 0 (x), x l < x < x r, x 2 u(x l, t) = x 2 v(x l, t), t > 0 xu(x r, t) = xw(x r, t), t > 0, xu(x 2 r, t) = xw(x 2 r, t), t > 0. Left exterior problem ( t + x 3 )v = 0, x < x l, t > 0, v(x, 0) = 0, x < x l, v(x l, t) = u(l, t), t > 0, xv(x l, t) = xu(l, t), t > 0, v 0, x. Right exterior problem ( t + x 3 )w = 0, x > xr, t > 0, w(x, 0) = 0, x > x r, w(x, x r) = u(x, x r), t > 0, w 0, x.
19 The Airy equation 3 Laplace transform w.r.t time tf + xf 3 = 0 ω ˆf + x 3 ˆf = 0 3 Solution: ˆf(x, ω) = c k (ω)e λk(ω)x with λ k (ω) = j k 1 3 ω, j = e 2iπ/3 k=1 Re λ 1(ω) < 0, Re λ 2(ω) > 0, Re λ 3(ω) > 0.
20 The Airy equation 4 Right exterior problem : c 2 = c 3 = 0 so ŵ(x, ω) = c 1e λ 1(ω)x. Using continuity, we have ŵ(x, ω) = L {u(x r, )}(ω)e λ 1(x x r). Taking first and second order derivatives and xŵ(x, ω) = λ 1L {u(x r, )}(ω)e λ 1(x x r) 2 xŵ(x, ω) = λ 2 1L {u(x r, )}(ω)e λ 1(x x r). Continuity of the first and second order derivatives xl {u(x r, )}(ω) = λ 1L {u(x r, )}(ω), 2 xl {u(x r, )}(ω) = λ 2 1L {u(x r, )}(ω).
21 The Airy equation 5 Inverse Laplace transform for right exterior problem u(x r, t) I 2/3 t u xx(x r, t) = 0, u x(x r, t) + I 1/3 t u xx(x r, t) = 0, where I α t with α is the fractional integral operator defined by I α t f(t) = 1 Γ(α) t 6 Left exterior problem : c 1 = 0 and thus, we have 0 (t τ) α 1 f(τ) dτ, t > 0. û(x l, ω) + 1 λ 1 û x(x l, ω) + 1 λ 2 1 û xx(x l, ω) = 0. 7 Inverse Laplace transform for left exterior problem u(x l, t) I 1/3 t u x(x l, t) + I 2/3 t u xx(x l, t) = 0.
22 The Airy equation The new IBVP is The Airy Eq. in Ω u t + u xxx = 0, (x, t) [x l, x r] (0, T ) u I 1/3 t u x + I 2/3 t u xx = 0, x = x l, t > 0, u I 2/3 t u xx = 0, x = x r, t > 0, u + I 1/3 t u xx = 0, x = x r, t > 0, u(x, 0) = u 0(x), x [x l, x r]
23 General Airy Eq. tu + U 1 xu + U 2 xu 3 = 0, on R [0; T ]. u(x, 0) = u 0(x), x R. We have to find the roots of the (depressed) cubic equation ω + U 1λ + U 2λ 3 = 0. We still have the separation property of the roots
24 General Airy Eq. The new IBVP is The General Airy Eq. in Ω u t + U 1u x + U 2u xxx = 0, (x, t) [x l, x r] (0, T ) ( ) ( u(t, x l ) U 2 L 1 u x(t, x l ) U 2 L 1 λ 1 (s) 2 s ( ) u(t, x r) L 1 1 λ 1 u (s) xx(t, x r) = 0, t > 0, 2 ( u x(t, x r) L 1 1 u(x, 0) = u 0(x), λ 1 (s) λ 1 (s) s ) u xx(t, x r) = 0, x = x r, t > 0, x [x l, x r]. ) u xx(t, x l ) = 0, t > 0, with λ k (ω) = j k 1 ζ(ω) 1 U U 2 j k 1 ζ(ω) and ζ(ω) = 1 ω ω /3 U 2 U 2 27 ( U1 U 2 ) 3 1/3.
25 The 1D Schrödinger Eq. TBC on {x l, x r} : Remark: nv + e iπ/4 1/2 t v = 0, on {x l, x r} [0; T ]. 2 x + i t = ( n + i i t)( n i i t) The Schrödinger Eq. in Ω i tψ + ψ = 0, (x, t) [x l, x r] [0, T ] (S) nψ + e iπ/4 1/2 t ψ = 0, on {x l, x r} [0, T ], ψ(x, 0) = ψ 0(x), x [x l, x r]
26 Outline of the talk 1 Why Artifical Boundary Conditions? 2 Derivation of some TBC 3 Numerical approximation
27 Standard scheme Crank Nicolson u n+1 u n + ( U 1 x + U 2 3 ) u n+1 + u n x = 0, on [x l, x r] [0; N] t 2 u 0 (x) = u 0(x), x [x l, x r], T.B.C., x {x l, x r}. If U 1 = 0, approximation of I α t by C. Zheng, X. Wen, and H. Han, (2008) β α,j = { ( t) α α(α + 1)Γ(α) I α t f(t m) m β α,m jf j, j=0 1, j = 0, (j + 1) α+1 + (j 1) α+1 2j α+1, j > 0. Followed path: Derivation of continuous TBC, then discretization.
28 Full discretization Collaboration with Matthias Ehrhardt (Wuppertal) and Ingrid Lacroix-Violet (Lille) Other way: Full discretization, then derivation of discrete TBC adapted to the numerical scheme. Discretization parameters (t n) 0 n N uniform discretization to [0, T ] s.t. t n = n t with t = T/N 0 = t 0 < t 1 < < t N 1 < t N = T. (x j) 0 j J uniform discretization to [a, b] s.t. x j = a + j x with x = (b a)/j u (n) j approximation to u(t n, x j) a = x 0 < x 1 < < x J 1 < x J = b.
29 Numerical scheme (R-CN) (R-CN): Case U 1 = 0 u (n+1) j u (n) j t Proprerties of the scheme (R-CN) inconditionally stable E R CN = O( x + t 2 ) 4 nodes scheme ( + U2 u (n+1) 2( x) 3 j+2 3u (n+1) j+1 + 3u (n+1) j ) u (n+1) j 1 ( + U2 u (n) 2( x) 3 j+2 3u(n) j+1 + 3u(n) j u (n) j 1 ) = 0.
30 TBCs for (R-CN) Definition of the Z transform û(z) = Z{(u n ) n}(z) = u k z k, z > R 1, k=0 Application of the Z transform to (R-CN) (û j = û j(z) = Z{(u (n) j ) n)}) ) û j+2 3û j+1 + (3 + 2( x)3 z 1 û j û j 1 = 0, 1 j J 2. U 2 t z + 1 Solution on the external domain û j = 3 k=1 where l = l(z) is solution to l 3 3l 2 + c k (z) l j k (z), j 1 ou j J 2, (3 + 2( x)3 U 2 t ) z 1 l 1 = 0. z + 1
31 TBCs for (R-CN) separation property of the discrete roots l 1(z) < 1, l 2(z) > 1, l 3(z) > 1, for all z, decay property = û j(z) = c 2(z) l j 2(z) + c 3(z) l j 3(z), j 1 û j(z) = c 1(z) l j 1(z), j J 2
32 TBCs for (R-CN) - Left B.C Scheme u (n+1) j û j satisfies u (n) j t = for j = 1 û j(z) = c 2(z) l j 2(z) + c 3(z) l j 3(z), j 1 ( + U2 u (n+1) 2( x) 3 j+2 3u (n+1) j+1 + 3u (n+1) j ) u (n+1) j 1 ( + U2 u (n) 2( x) 3 j+2 3u(n) j+1 + 3u(n) j u (n) j 1 = one need only one left BC û j+1(z) ( l 2(z) + l 3(z) ) û j(z) + l 2(z)l 3(z) û j 1(z) = 0. ) = 0. Z 1 {l 2(z)l 3(z)} d u (n) 0 Z 1 {l 2(z) + l 3(z))} d u (n) 1 + u (n) 2 = 0, n = 0, 1, 2..., where P d u (n) i = n k=0 P (k) u (n k) i.
33 TBCs for (R-CN) - right B.C scheme u (n+1) j u (n) j t û j verifies two relations = for j = J 2 û j(z) = c 1(z) l j 1(z), j J 2 ( + U2 u (n+1) 2( x) 3 j+2 3u (n+1) j+1 + 3u (n+1) j ) u (n+1) j 1 ( + U2 u (n) 2( x) 3 j+2 3u(n) j+1 + 3u(n) j u (n) j 1 = two BCs are needed û j+2(z) = l 1(z) 2 û j(z), et û j+1(z) = l 1(z) û j(z), u (n) J Z 1 {l 2 1(z)} d u (n) J 2 = 0, where P d u (n) i = n k=0 P (k) u (n k) i. ) = 0. u(n) J 1 Z 1 {l 1(z)} d u (n) J 2 = 0, n = 0, 1, 2...,
34 TBCs for (R-CN) - conclusion Numerical scheme for 1 j J 2 and n = 0, 1, 2... u (n+1) j u (n) j t ( + U2 u (n+1) 2( x) 3 j+2 3u (n+1) j+1 + 3u (n+1) j ) u (n+1) j 1 ( + U2 u (n) 2( x) 3 j+2 3u(n) j+1 + 3u(n) j u (n) j 1 ) = 0. Left BC Z 1 {l 2(z)l 3(z)} d u (n) 0 Z 1 {l 2(z) + l 3(z))} d u (n) 1 + u (n) 2 = 0, n = 0, 1, 2..., Right BC u (n) J Z 1 {l 2 1(z)} d u (n) J 2 = 0, u(n) J 1 Z 1 {l 1(z)} d u (n) J 2 = 0, n = 0, 1, 2...,
35 Scheme (C-CN) Numerical scheme (R-CN): Case U 1 0 u (n+1) j u (n) j t ( ) ( ) + U1 u (n+1) j+1 u (n+1) j 1 + U1 u (n) j+1 4 x 4 x u(n) j 1 ) ( + U2 u (n+1) 4( x) 3 j+2 2u (n+1) j+1 + 2u (n+1) j 1 + U2 4( x) 3 ( u (n) j+2 2u(n) j+1 + 2u(n) j 1 u(n) j 2 u (n+1) j 2 ) = 0. Properties inconditionally stable E C CN = O( x 2 + t 2 ) five nodes We still need two BCs on right boundary but also two on the left.
36 TBCs for (C-CN) Numerical scheme for 1 j J 2 and n = 0, 1, 2... u (n+1) j u (n) j t ( ) ( ) + U1 u (n+1) j+1 u (n+1) j 1 + U1 u (n) j+1 4 x 4 x u(n) j 1 ) ( + U2 u (n+1) 4( x) 3 j+2 2u (n+1) j+1 + 2u (n+1) j 1 + U2 4( x) 3 ( u (n) j+2 2u(n) j+1 + 2u(n) j 1 u(n) j 2 u (n+1) j 2 ) = 0. Left BC Right BC Z 1 {λ 1} d u (n) 0 Z 1 {λ 2} d u (n) 1 + u (n) 2 = 0 Z 1 {λ 2 1} d u (n) 0 Z 1 {λ 2 2} d u (n) 2 + 2Z 1 {λ 2} d u (n) 3 u (n) 4 = 0 Z 1 {ρ 1} d u (n) J 2 Z 1 {ρ 2} d u (n) J 1 + u(n) J = 0 Z 1 {ρ 2 1} d u (n) J 4 Z 1 {ρ 2 2} d u (n) J 2 + 2Z 1 {ρ 2} d u (n) J 1 u(n) J = 0 where λ 1 = l 3(z)l 4(z), λ 2 = l 3(z) + l 4(z), ρ 1 = l 1(z)l 2(z), ρ 2 = l 1(z) + l 2(z).
37 Numerical procedure for inverse Z. A. Zisowsky,, Discrete transparent boundary conditions for systems of evolution equations, PhD, Technische Universität Berlin, 2003 Z of (u n) n : U(z) = u k z k, z > R. Inverse Z: u n = 1 2iπ u n = ρn 2π k=0 2π 0 S ρ U(z)z n 1 dz, ρ > R. U(ρe iϕ )e inϕ dϕ ρn N 1 U N k=0 z = ρe iϕ Recall the rule of discrete Fourier transform of ( f n ) N 1 n=0 F k = F{f n}(k) = where ω N = e i 2π N. U k = U(ρω k N ) = U(z k ) N 1 n=0 ( ρe i 2π N k) e i 2π N kn N 1 f nω nk N, f n = F 1 {F k }(n) = u n u N n = ρ n F 1 {U k }(n), 0 n < N. k=0 F k ω nk N
38 Example 1 C. Zheng, X. Wen, H. Han,, Numerical Solution to a Linearized KdV Equation on Unbounded Domain, Numer. Meth. Part. Diff. Eqs., 2008 Equation u t + u xxx = 0, x R, u(0, x) = e x2, x R, u 0, x. Exact solution E(t, x) = 3 1 ( ) x Ai 3t 3, 3t u exact(t, x) = E(t, x) e x2, where denotes the convolution product. Discretization parameters T = 4, [a, b] = [ 6, 6], t = 4/2560, x = 12/5000, r = 1.001
39 Example 1
40 Example 1
41 Example 1 - (R-CN) rel.errl x rel.errtm 10 3 x 10 4 N=640 N=1280 N=2560 N=5120 N= N=640 N=1280 N=2560 N=5120 N= x x Error functions rel.errt m = ( ( max e (n)), rel.errl2 = t 0<n<N where e (n) is the relative error l 2 at time t = n t: e (n) = u (n) exact u num (n) / 2 u (n) exact 2 ) 1/2 N (e (n) ) 2 n=1
42 Example 1 - (C-CN) rel.errl rel.errtm x x 2 N= N=1280 N=2560 N=5120 N= x N= N=1280 N=2560 N=5120 N= x Error functions rel.errt m = ( ( max e (n)), rel.errl2 = t 0<n<N where e (n) is the relative error l 2 at time t = n t: e (n) = u (n) exact u num (n) / 2 u (n) exact 2 ) 1/2 N (e (n) ) 2 n=1
43 Example 1 - (R-CN) rel.errl rel.errtm t t t t Error functions rel.errt m = ( ( max e (n)), rel.errl2 = t 0<n<N where e (n) is the relative error l 2 at time t = n t: e (n) = u (n) exact u num (n) / 2 u (n) exact 2 ) 1/2 N (e (n) ) 2 n=1
44 Example 1 - (C-CN) rel.errl rel.errtm t t t t Error functions rel.errt m = ( ( max e (n)), rel.errl2 = t 0<n<N where e (n) is the relative error l 2 at time t = n t: e (n) = u (n) exact u num (n) / 2 u (n) exact 2 ) 1/2 N (e (n) ) 2 n=1
45 Example 2 W.L. Briggs, T. Sarie,, Finite difference solutions of dispersive partial differential equations, Math. Comput. Simul., 1983 Equation u t + u x + u xxx = 0, x R, u(0, x) = exp( 8(x 5) 2 ) sin(50π/4), x R, Exacte solution û exact(t, ξ) = û 0 exp( (iξ iξ 3 )t) Discretization parameters T = , [a, b] = [0, 10], t = T/2560, x = 10/5000
46 Example 2
47 Conclusion TBC for continuous and discrete linear Korteweg de Vries equation modified KdV : Zheng (06) u t ± 6u 2 u x + u xxx = 0 Use of inverse scattering to get exact TBC. Example (Zheng) : solitary waves generated by an initial Gaussian profile Numerical simulation u 0(x) = exp ( 1.5x 2 of a modified KdV equation 333 ). Future: Fig. extension 5 Evolution of toa Gaussian nonlinear wave packet problem: consider U1u x as a nonlinear term. of the wave field. We see that the interaction of solitons is very elastic. After collision, two solitons remains their shapes and velocities, but their phases are
Numerical aspects of the nonlinear Schro dinger equation. Christophe Besse
Numerical aspects of the nonlinear Schro dinger equation by Christophe Besse Laboratoire Paul Painleve, Universite Lille 1, CNRS, Inria Simpaf Team Lille Outline of the talk 1 Motivation 2 Numerical methods
More informationOpen boundary conditions and computational schemes for Schrödinger equations with general potentials and nonlinearities
1 Open boundary conditions and computational schemes for Schrödinger equations with general potentials and nonlinearities Xavier Antoine, Pauline Klein Institut Elie Cartan Nancy, Nancy-Université, CNRS
More informationAbsorbing boundary conditions for Schrödinger equations. Christophe Besse 1
Absorbing boundary conditions for Schrödinger equations by Christophe Besse 1 1 Laboratoire Paul Painlevé, Université Lille 1, Équipe INRIA Simpaf Workshop Stochastic Schrödinger Equation Paris Collaborators:
More informationNumerical methods for a fractional diffusion/anti-diffusion equation
Numerical methods for a fractional diffusion/anti-diffusion equation Afaf Bouharguane Institut de Mathématiques de Bordeaux (IMB), Université Bordeaux 1, France Berlin, November 2012 Afaf Bouharguane Numerical
More informationDiscrete Transparent Boundary Conditions for the Schrödinger Equation. Matthias Ehrhardt, Anton Arnold
Discrete Transparent Boundary Conditions for the Schrödinger Equation Matthias Ehrhardt, Anton Arnold Fakultät für Mathematik und Informatik, Universität des Saarlandes, Postfach 15 11 5, D-6641 Saarbrücken,
More informationNumerical solution of the Schrödinger equation on unbounded domains
Numerical solution of the Schrödinger equation on unbounded domains Maike Schulte Institute for Numerical and Applied Mathematics Westfälische Wilhelms-Universität Münster Cetraro, September 6 OUTLINE
More informationComputation of Nonlinear Schrödinger Equation on an Open Waveguide Terminated by a PML
Copyright 20 Tech Science Press CMES, vol.7, no.4, pp.347-362, 20 Computation of Nonlinear Schrödinger Equation on an Open Waveguide Terminated by a PML Jianxin Zhu, Zheqi Shen Abstract: It is known that
More informationContinuous limits and integrability for a semidiscrete system Zuo-nong Zhu Department of Mathematics, Shanghai Jiao Tong University, P R China
Continuous limits and integrability for a semidiscrete system Zuo-nong Zhu Department of Mathematics, Shanghai Jiao Tong University, P R China the 3th GCOE International Symposium, Tohoku University, 17-19
More informationPresenter: Noriyoshi Fukaya
Y. Martel, F. Merle, and T.-P. Tsai, Stability and Asymptotic Stability in the Energy Space of the Sum of N Solitons for Subcritical gkdv Equations, Comm. Math. Phys. 31 (00), 347-373. Presenter: Noriyoshi
More informationWave Equation With Homogeneous Boundary Conditions
Wave Equation With Homogeneous Boundary Conditions MATH 467 Partial Differential Equations J. Robert Buchanan Department of Mathematics Fall 018 Objectives In this lesson we will learn: how to solve the
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 informationA Friendly Review of Absorbing Boundary Conditions and Perfectly Matched Layers for Classical and Relativistic Quantum Waves Equations
A Friendly Review of Absorbing Boundary Conditions and Perfectly Matched Layers for Classical and Relativistic Quantum Waves Equations X. Antoine a, E. Lorin b,c and Q. Tang a a Institut Elie Cartan de
More informationDUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs
DUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs MATH 4354 Fall 2005 December 5, 2005 1 Duhamel s Principle
More informationSolutions of differential equations using transforms
Solutions of differential equations using transforms Process: Take transform of equation and boundary/initial conditions in one variable. Derivatives are turned into multiplication operators. Solve (hopefully
More informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More informationTowards accurate artificial boundary conditions for nonlinear PDEs through examples
Towards accurate artificial boundary conditions for nonlinear PDEs through examples Xavier Antoine, Christophe Besse, Jérémie Szeftel Abstract. The aim of this paper is to give a comprehensive review of
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
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 informationDiffusion on the half-line. The Dirichlet problem
Diffusion on the half-line The Dirichlet problem Consider the initial boundary value problem (IBVP) on the half line (, ): v t kv xx = v(x, ) = φ(x) v(, t) =. The solution will be obtained by the reflection
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationArtificial boundary conditions for the linearized Benjamin-Bona-Mahony equation
Artificial boundary conditions for the linearized Benamin-Bona-Mahony equation Christophe Besse, Benoît Mésognon-Gireau, Pascal Noble To cite this version: Christophe Besse, Benoît Mésognon-Gireau, Pascal
More informationA quantum heat equation 5th Spring School on Evolution Equations, TU Berlin
A quantum heat equation 5th Spring School on Evolution Equations, TU Berlin Mario Bukal A. Jüngel and D. Matthes ACROSS - Centre for Advanced Cooperative Systems Faculty of Electrical Engineering and Computing
More informationOptimized Schwarz method for the linearized KdV equation
Optimized Schwarz method for the linearized KdV equation Joao Guilherme Caldas Steinstraesser, Rodrigo Cienfuegos, José Daniel Galaz Mora, Antoine Rousseau To cite this version: Joao Guilherme Caldas Steinstraesser,
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationA Novel Nonlinear Evolution Equation Integrable by the Inverse Scattering Method
Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part, 384 39 A Novel Nonlinear Evolution Equation Integrable by the Inverse Scattering Method Vyacheslav VAKHNENKO and John PARKES
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More information6 Non-homogeneous Heat Problems
6 Non-homogeneous Heat Problems Up to this point all the problems we have considered for the heat or wave equation we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationBielefeld Course on Nonlinear Waves - June 29, Department of Mathematics University of North Carolina, Chapel Hill. Solitons on Manifolds
Joint work (on various projects) with Pierre Albin (UIUC), Hans Christianson (UNC), Jason Metcalfe (UNC), Michael Taylor (UNC), Laurent Thomann (Nantes) Department of Mathematics University of North Carolina,
More informationApplication of linear combination between cubic B-spline collocation methods with different basis for solving the KdV equation
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 3, 016, pp. 191-04 Application of linear combination between cubic B-spline collocation methods with different basis
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
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 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 informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationThe Orchestra of Partial Differential Equations. Adam Larios
The Orchestra of Partial Differential Equations Adam Larios 19 January 2017 Landscape Seminar Outline 1 Fourier Series 2 Some Easy Differential Equations 3 Some Not-So-Easy Differential Equations Outline
More informationTADAHIRO OH 0, 3 8 (T R), (1.5) The result in [2] is in fact stated for time-periodic functions: 0, 1 3 (T 2 ). (1.4)
PERIODIC L 4 -STRICHARTZ ESTIMATE FOR KDV TADAHIRO OH 1. Introduction In [], Bourgain proved global well-posedness of the periodic KdV in L T): u t + u xxx + uu x 0, x, t) T R. 1.1) The key ingredient
More informationFrom the N-body problem to the cubic NLS equation
From the N-body problem to the cubic NLS equation François Golse Université Paris 7 & Laboratoire J.-L. Lions golse@math.jussieu.fr Los Alamos CNLS, January 26th, 2005 Formal derivation by N.N. Bogolyubov
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationThe Helmholtz Equation
The Helmholtz Equation Seminar BEM on Wave Scattering Rene Rühr ETH Zürich October 28, 2010 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness
More informationRiesz bases of Floquet modes in semi-infinite periodic waveguides and implications
Riesz bases of Floquet modes in semi-infinite periodic waveguides and implications Thorsten Hohage joint work with Sofiane Soussi Institut für Numerische und Angewandte Mathematik Georg-August-Universität
More informationHylomorphic solitons and their dynamics
Hylomorphic solitons and their dynamics Vieri Benci Dipartimento di Matematica Applicata U. Dini Università di Pisa 18th May 2009 Vieri Benci (DMA-Pisa) Hylomorphic solitons 18th May 2009 1 / 50 Types
More informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationIntegral Representation Formula, Boundary Integral Operators and Calderón projection
Integral Representation Formula, Boundary Integral Operators and Calderón projection Seminar BEM on Wave Scattering Franziska Weber ETH Zürich October 22, 2010 Outline Integral Representation Formula Newton
More informationCollocation Method with Quintic B-Spline Method for Solving the Hirota equation
NTMSCI 1, No. 1, 1-1 (016) 1 Journal of Abstract and Computational Mathematics http://www.ntmsci.com/jacm Collocation Method with Quintic B-Spline Method for Solving the Hirota equation K. R. Raslan 1,
More informationCompacton-like solutions in some nonlocal hydrodynamic-type models
Compacton-like solutions in some nonlocal hydrodynamic-type models Vsevolod Vladimirov AGH University of Science and technology, Faculty of Applied Mathematics Protaras, October 26, 2008 WMS AGH Compactons
More informationMath 575-Lecture 26. KdV equation. Derivation of KdV
Math 575-Lecture 26 KdV equation We look at the KdV equations and the so-called integrable systems. The KdV equation can be written as u t + 3 2 uu x + 1 6 u xxx = 0. The constants 3/2 and 1/6 are not
More informationRecurrence in the KdV Equation? A. David Trubatch United States Military Academy/Montclair State University
Nonlinear Physics: Theory and Experiment IV Gallipoli, Lecce, Italy 27 June, 27 Recurrence in the KdV Equation? A. David Trubatch United States Military Academy/Montclair State University david.trubatch@usma.edu
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem and Its Difference Approximations 1-D Initial Boundary Value
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationStability of periodic waves in the defocusing cubic NLS equation
Stability of periodic waves in the defocusing cubic NLS equation Thierry Gallay 1 and Dmitry Pelinovsky 2 1 Institut Fourier, Université de Grenoble 1, 38402 Saint-Martin-d Hères, France 2 Department of
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationFORCED OSCILLATIONS OF A CLASS OF NONLINEAR DISPERSIVE WAVE EQUATIONS AND THEIR STABILITY
Jrl Syst Sci & Complexity (2007) 20: 284 292 FORCED OSCILLATIONS OF A CLASS OF NONLINEAR DISPERSIVE WAVE EQUATIONS AND THEIR STABILITY Muhammad USMAN Bingyu ZHANG Received: 14 January 2007 Abstract It
More informationInterior feedback stabilization of wave equations with dynamic boundary delay
Interior feedback stabilization of wave equations with dynamic boundary delay Stéphane Gerbi LAMA, Université Savoie Mont-Blanc, Chambéry, France Journée d EDP, 1 er Juin 2016 Equipe EDP-Contrôle, Université
More informationA NEW APPROACH FOR SOLITON SOLUTIONS OF RLW EQUATION AND (1+2)-DIMENSIONAL NONLINEAR SCHRÖDINGER S EQUATION
A NEW APPROACH FOR SOLITON SOLUTIONS OF RLW EQUATION AND (+2-DIMENSIONAL NONLINEAR SCHRÖDINGER S EQUATION ALI FILIZ ABDULLAH SONMEZOGLU MEHMET EKICI and DURGUN DURAN Communicated by Horia Cornean In this
More informationWeek 4 Lectures, Math 6451, Tanveer
1 Diffusion in n ecall that for scalar x, Week 4 Lectures, Math 6451, Tanveer S(x,t) = 1 exp [ x2 4πκt is a special solution to 1-D heat equation with properties S(x,t)dx = 1 for t >, and yet lim t +S(x,t)
More informationContinuous and Discrete Homotopy Operators with Applications in Integrability Testing. Willy Hereman
Continuous and Discrete Homotopy Operators with Applications in Integrability Testing Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado, USA http://www.mines.edu/fs
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationDISSIPATIVE INITIAL BOUNDARY VALUE PROBLEM FOR THE BBM-EQUATION
Electronic Journal of Differential Equations, Vol. 8(8), No. 149, pp. 1 1. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) DISSIPATIVE
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationOn the Whitham Equation
On the Whitham Equation Henrik Kalisch Department of Mathematics University of Bergen, Norway Joint work with: Handan Borluk, Denys Dutykh, Mats Ehrnström, Daulet Moldabayev, David Nicholls Research partially
More informationSpectrum and Exact Controllability of a Hybrid System of Elasticity.
Spectrum and Exact Controllability of a Hybrid System of Elasticity. D. Mercier, January 16, 28 Abstract We consider the exact controllability of a hybrid system consisting of an elastic beam, clamped
More informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationNumerical computation of the finite-genus solutions of the Korteweg-de Vries equation via Riemann Hilbert problems
Numerical computation of the finite-genus solutions of the Korteweg-de Vries equation via Riemann Hilbert problems Thomas Trogdon 1 and Bernard Deconinck Department of Applied Mathematics University of
More informationGroup analysis, nonlinear self-adjointness, conservation laws, and soliton solutions for the mkdv systems
ISSN 139-5113 Nonlinear Analysis: Modelling Control, 017, Vol., No. 3, 334 346 https://doi.org/10.15388/na.017.3.4 Group analysis, nonlinear self-adjointness, conservation laws, soliton solutions for the
More informationDynamics of Strongly Nonlinear Internal Solitary Waves in Shallow Water
Dynamics of Strongly Nonlinear Internal Solitary Waves in Shallow Water By Tae-Chang Jo and Wooyoung Choi We study the dynamics of large amplitude internal solitary waves in shallow water by using a strongly
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J NUMER ANAL Vol, No 3, pp 93 95 c 00 Society for Industrial and Applied Mathematics A RELAXATION SCHEME FOR THE NONLINEAR SCHRÖDINGER EQUATION CHRISTOPHE BESSE Abstract In this paper, we present
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 informationSome Aspects of Solutions of Partial Differential Equations
Some Aspects of Solutions of Partial Differential Equations K. Sakthivel Department of Mathematics Indian Institute of Space Science & Technology(IIST) Trivandrum - 695 547, Kerala Sakthivel@iist.ac.in
More informationTHE LAX PAIR FOR THE MKDV HIERARCHY. Peter A. Clarkson, Nalini Joshi & Marta Mazzocco
Séminaires & Congrès 14, 006, p. 53 64 THE LAX PAIR FOR THE MKDV HIERARCHY by Peter A. Clarkson, Nalini Joshi & Marta Mazzocco Abstract. In this paper we give an algorithmic method of deriving the Lax
More informationIntroduction and some preliminaries
1 Partial differential equations Introduction and some preliminaries A partial differential equation (PDE) is a relationship among partial derivatives of a function (or functions) of more than one variable.
More informationGLOBAL WELL-POSEDNESS OF NLS-KDV SYSTEMS FOR PERIODIC FUNCTIONS
Electronic Journal of Differential Equations, Vol. 66), o. 7, pp. 1. ISS: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu login: ftp) GLOBAL WELL-POSEDESS
More informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
More informationAdomian Polynomial and Elzaki Transform Method of Solving Third Order Korteweg-De Vries Equations
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 15, Number 3 (2019), pp 261 277 c Research India Publications http://www.ripublication.com Adomian Polynomial and Elzaki Transform
More informationComputer Physics Communications. Computational methods for the dynamics of the nonlinear Schrödinger/Gross-Pitaevskii equations
Computer Physics Communications (13) 1 3 Computer Physics Communications Computational methods for the dynamics of the nonlinear Schrödinger/Gross-Pitaevskii equations Xavier Antoine a,b, Weizhu Bao 1c,d,
More informationA converse to the Ambrosetti-Prodi theorem
A converse to the Ambrosetti-Prodi theorem Marta Calanchi, Università di Milano with Carlos Tomei and André Zaccur (PUC-Rio, Brazil) Varese, RISM, September 2015 The Ambrosetti-Prodi Theorem 1/14 The Ambrosetti-Prodi
More informationApplications of the Maximum Principle
Jim Lambers MAT 606 Spring Semester 2015-16 Lecture 26 Notes These notes correspond to Sections 7.4-7.6 in the text. Applications of the Maximum Principle The maximum principle for Laplace s equation is
More informationAn inverse source problem in optical molecular imaging
An inverse source problem in optical molecular imaging Plamen Stefanov 1 Gunther Uhlmann 2 1 2 University of Washington Formulation Direct Problem Singular Operators Inverse Problem Proof Conclusion Figure:
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 informationConvergence of operator splitting for the KdV equation
Convergence of operator splitting for the KdV equation H. Holden Norwegian University of Science and Technology Trondheim, Norway Joint work with K.H. Karlsen, N. H. Risebro, and T. Tao Institute for Mathematics
More informationApplications of the periodic unfolding method to multi-scale problems
Applications of the periodic unfolding method to multi-scale problems Doina Cioranescu Université Paris VI Santiago de Compostela December 14, 2009 Periodic unfolding method and multi-scale problems 1/56
More informationThermodynamic limit for a system of interacting fermions in a random medium. Pieces one-dimensional model
Thermodynamic limit for a system of interacting fermions in a random medium. Pieces one-dimensional model Nikolaj Veniaminov (in collaboration with Frédéric Klopp) CEREMADE, University of Paris IX Dauphine
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationHamiltonian partial differential equations and Painlevé transcendents
Winter School on PDEs St Etienne de Tinée February 2-6, 2015 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN SISSA, Trieste Lecture 2 Recall: the main goal is to compare
More informationFunctions of a Complex Variable (S1) Lecture 11. VII. Integral Transforms. Integral transforms from application of complex calculus
Functions of a Complex Variable (S1) Lecture 11 VII. Integral Transforms An introduction to Fourier and Laplace transformations Integral transforms from application of complex calculus Properties of Fourier
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More informationThere are five problems. Solve four of the five problems. Each problem is worth 25 points. A sheet of convenient formulae is provided.
Preliminary Examination (Solutions): Partial Differential Equations, 1 AM - 1 PM, Jan. 18, 16, oom Discovery Learning Center (DLC) Bechtel Collaboratory. Student ID: There are five problems. Solve four
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 informationAlmost Sure Well-Posedness of Cubic NLS on the Torus Below L 2
Almost Sure Well-Posedness of Cubic NLS on the Torus Below L 2 J. Colliander University of Toronto Sapporo, 23 November 2009 joint work with Tadahiro Oh (U. Toronto) 1 Introduction: Background, Motivation,
More informationExistence and exponential stability of the damped wave equation with a dynamic boundary condition and a delay term.
Existence and exponential stability of the damped wave equation with a dynamic boundary condition and a delay term. Stéphane Gerbi LAMA, Université de Savoie, Chambéry, France Fachbereich Mathematik und
More informationORBITAL STABILITY OF SOLITARY WAVES FOR A 2D-BOUSSINESQ SYSTEM
Electronic Journal of Differential Equations, Vol. 05 05, No. 76, pp. 7. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu ORBITAL STABILITY OF SOLITARY
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 informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationGroup Invariant Solutions of Complex Modified Korteweg-de Vries Equation
International Mathematical Forum, 4, 2009, no. 28, 1383-1388 Group Invariant Solutions of Complex Modified Korteweg-de Vries Equation Emanullah Hızel 1 Department of Mathematics, Faculty of Science and
More informationarxiv:patt-sol/ v1 25 Sep 1995
Reductive Perturbation Method, Multiple Time Solutions and the KdV Hierarchy R. A. Kraenkel 1, M. A. Manna 2, J. C. Montero 1, J. G. Pereira 1 1 Instituto de Física Teórica Universidade Estadual Paulista
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 information