Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
|
|
- Arline Murphy
- 5 years ago
- Views:
Transcription
1 Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada Physics seminar November 10, 2011 A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
2 Outline 1 Conservation Laws Definition and Examples Applications 2 Variational Principles 3 Symmetries and Noether s Theorem Symmetries and Variational Symmetries Noether s Theorem, Examples Limitations of Noether s Theorem 4 Direct Construction Method for Conservation Laws The Idea Completeness A Detailed Example 5 Examples of Construction of Conservation Laws Symbolic Software (Maple) KdV Surfactant Dynamics Equations 6 Conclusions A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
3 Outline 1 Conservation Laws Definition and Examples Applications 2 Variational Principles 3 Symmetries and Noether s Theorem Symmetries and Variational Symmetries Noether s Theorem, Examples Limitations of Noether s Theorem 4 Direct Construction Method for Conservation Laws The Idea Completeness A Detailed Example 5 Examples of Construction of Conservation Laws Symbolic Software (Maple) KdV Surfactant Dynamics Equations 6 Conclusions A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
4 Conservation Laws: Some Definitions Notation: The total derivative by x: The partial derivative by x assuming all chain rules carried out as needed. Dependent variables should be taken care of. Example: u = u(x, y), g = g(u), then D x[xg] [xg(u(x, y))] x = x [xg] + g [xg] u [g(u)] u(x, y). x A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
5 Conservation Laws: Some Definitions Conservation law: Given: some model. Independent variables: x = (t, x, y,...); dependent variables: u = (u, v,...). A conservation law: a divergence expression equal to zero D tθ(x, u,...) + D xψ 1 (x, u,...) + D y Ψ 2 (x, u,...) + = 0. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
6 Conservation Laws: Some Definitions Conservation law: Given: some model. Independent variables: x = (t, x, y,...); dependent variables: u = (u, v,...). A conservation law: a divergence expression equal to zero D tθ(x, u,...) + D xψ 1 (x, u,...) + D y Ψ 2 (x, u,...) + = 0. Time-independent: D i Ψ i div x,y,...ψ = 0. Time-dependent: A conserved quantity: D tθ + div x,y,...ψ = 0. D t Θ dv = 0. V A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
7 ODE Models: Conserved Quantities An ODE: Dependent variable: u = u(t); A conservation law D tf (t, u, u,...) = d dt DtF (t, u, u,...) = 0. yields a conserved quantity (constant of motion): F (t, u, u,...) = C = const. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
8 ODE Models: Conserved Quantities An ODE: Dependent variable: u = u(t); A conservation law D tf (t, u, u,...) = d dt DtF (t, u, u,...) = 0. yields a conserved quantity (constant of motion): F (t, u, u,...) = C = const. Example: Harmonic oscillator, spring-mass system Independent variable: t, dependent: x(t). ODE: ẍ(t) + ω 2 x(t) = 0; ω 2 = k/m = const. ( d mẋ 2 (t) Conservation law: + kx ) 2 (t) = 0. dt 2 2 Conserved quantity: energy. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
9 ODE Models: Conserved Quantities Example: ODE integration An ODE: Three independent conserved quantities: K (x) = 2 (K (x)) 2 K(x) (K (x)) 2 K (x). K(x)K (x) KK (K ) 2 = C1, KK ln K (K ) 2 ln K = C 2, yield complete ODE integration. xkk + KK (K ) 2 x = C 3. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
10 PDE Models Example 1: small string oscillations, 1D wave equation Independent variables: x, t; dependent: u(x, t). u tt = c 2 u xx, c 2 = T /ρ. Conservation laws: Momentum: D t(ρu t) D x(tu x) = 0; Conserved quantity: total momentum M = ρu tdx = const. ) D x(tu tu x) = 0; ( ρu 2 Energy: D t t 2 + Tu2 x 2 Conserved quantity: total energy E = ( ρut Tu2 x 2 ) dx = const. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
11 PDE Models Example 2: Adiabatic motion of an ideal gas in 3D Independent variables: t; x = (x 1, x 2, x 3 ) D R 3. Dependent: ρ(x, t), v 1 (x, t), v 2 (x, t), v 3 (x, t), p(x, t). Equations: D tρ + D j (ρv j ) = 0, ρ(d t + v j D j )v i + D i p = 0, i = 1, 2, 3, ρ(d t + v j D j )p + γρpd j v j = 0. Conservation laws: Mass: D tρ + D j (ρv j ) = 0, Momentum: D t(ρv i ) + D j (ρv i v j + pδ ij ) = 0, i = 1, 2, 3, Energy: D t(e) + D j ( v j (E + p) ) = 0, E = 1 2 ρ v 2 + p γ 1. Angular momentum, etc. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
12 Applications of Conservation Laws ODEs Constants of motion. Integration. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
13 Applications of Conservation Laws ODEs Constants of motion. Integration. PDEs Direct physical meaning. Constants of motion. Differential constraints (div B = 0, etc.). Analysis: existence, uniqueness, stability. Nonlocally related PDE systems, exact solutions. Potentials, stream functions, etc.: V = (u, v), div V = u x + v y = 0, { u = Φy, v = Φ x. Modern numerical methods often require conserved forms. An infinite number of conservation laws can indicate integrability / linearization. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
14 Outline 1 Conservation Laws Definition and Examples Applications 2 Variational Principles 3 Symmetries and Noether s Theorem Symmetries and Variational Symmetries Noether s Theorem, Examples Limitations of Noether s Theorem 4 Direct Construction Method for Conservation Laws The Idea Completeness A Detailed Example 5 Examples of Construction of Conservation Laws Symbolic Software (Maple) KdV Surfactant Dynamics Equations 6 Conclusions A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
15 Variational Principles Action integral J[U] = L(x, U, U,..., k U)dx. Ω Principle of extremal action Variation of U: U(x) U(x) + εv(x). Require: variation of action δj J[U + εv] J[U] = Ω (δl) dx = O(ε2 ). A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
16 Variational Principles Action integral J[U] = L(x, U, U,..., k U)dx. Ω Principle of extremal action Variation of U: U(x) U(x) + εv(x). Require: variation of action δj J[U + εv] J[U] = Ω (δl) dx = O(ε2 ). Variation of Lagrangian δl = L(x, U + εv, U + ε v,..., k U + ε k v) L(x, U, U,..., k U) = ( ) L[U] ε U v σ + L[U] v σ σ Uj σ j + + L[U] v Uj σ j σ 1 j k + O(ε 2 ) 1 j k = ε(v σ E U σ (L[U]) +...) + O(ε 2 ), where E U σ are the Euler operators. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
17 Variational Principles Action integral J[U] = L(x, U, U,..., k U)dx. Ω Principle of extremal action Variation of U: U(x) U(x) + εv(x). Require: variation of action δj J[U + εv] J[U] = Ω (δl) dx = O(ε2 ). Euler-Lagrange equations: E u σ (L[u]) = L[u] u σ + + ( 1)k D j1 D jk L[u] u σ j 1 j k = 0, σ = 1,..., m. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
18 Variational Principles Example 1: Harmonic oscillator, x = x(t) L = 1 2 mẋ kx 2, E xl = m(ẍ + ω 2 x) = 0. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
19 Variational Principles Example 1: Harmonic oscillator, x = x(t) L = 1 2 mẋ kx 2, E xl = m(ẍ + ω 2 x) = 0. Example 2: Wave equation for u(x, t) L = 1 2 ρut T ux 2, E ul = ρ(u tt c 2 u xx) = 0. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
20 Variational Principles Example 1: Harmonic oscillator, x = x(t) L = 1 2 mẋ kx 2, E xl = m(ẍ + ω 2 x) = 0. Example 2: Wave equation for u(x, t) L = 1 2 ρut T ux 2, E ul = ρ(u tt c 2 u xx) = 0. Example 3: Klein-Gordon nonlinear equations for u(x, t) L = 1 2 utux + h(u), EuL = utx + h (u) = 0. Many other non-dissipative systems have variational formulations. A practical way to tell whether a DE system has a variational formulation: its linearization must be self-adjoint (symmetric). Relatively few models are! A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
21 Outline 1 Conservation Laws Definition and Examples Applications 2 Variational Principles 3 Symmetries and Noether s Theorem Symmetries and Variational Symmetries Noether s Theorem, Examples Limitations of Noether s Theorem 4 Direct Construction Method for Conservation Laws The Idea Completeness A Detailed Example 5 Examples of Construction of Conservation Laws Symbolic Software (Maple) KdV Surfactant Dynamics Equations 6 Conclusions A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
22 Symmetries of Differential Equations Consider a general DE system R σ [u] = R σ (x, u, u,..., k u) = 0, σ = 1,..., N with variables x = (x 1,..., x n ), u = u(x) = (u 1,..., u m ). Definition A transformation x = f (x, u; a) = x + aξ(x, u) + O(a 2 ), u = g(x, u; a) = u + aη(x, u) + O(a 2 ). depending on a parameter a is a point symmetry of R σ [u] if the equation is the same in new variables x, u. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
23 Symmetries of Differential Equations Consider a general DE system R σ [u] = R σ (x, u, u,..., k u) = 0, σ = 1,..., N with variables x = (x 1,..., x n ), u = u(x) = (u 1,..., u m ). Definition A transformation x = f (x, u; a) = x + aξ(x, u) + O(a 2 ), u = g(x, u; a) = u + aη(x, u) + O(a 2 ). depending on a parameter a is a point symmetry of R σ [u] if the equation is the same in new variables x, u. Example 1: translations The translation leaves KdV invariant: x = x + C, t = t, u = u u t + uu x + u xxx = 0 = u t + u u x + u x x x. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
24 Symmetries of Differential Equations Consider a general DE system R σ [u] = R σ (x, u, u,..., k u) = 0, σ = 1,..., N with variables x = (x 1,..., x n ), u = u(x) = (u 1,..., u m ). Definition A transformation x = f (x, u; a) = x + aξ(x, u) + O(a 2 ), u = g(x, u; a) = u + aη(x, u) + O(a 2 ). depending on a parameter a is a point symmetry of R σ [u] if the equation is the same in new variables x, u. Example 2: scaling Same for the scaling: One has x = αx, t = α 3 t, u = αu. u t + uu x + u xxx = 0 = u t + u u x + u x x x. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
25 Variational Symmetries Consider a general DE system R σ [u] = R σ (x, u, u,..., k u) = 0, σ = 1,..., N that follows from a variational principle with J[U] = Ω L(x, U, U,..., k U)dx. Definition A symmetry of R σ [u] given by x = f (x, u; a) = x + aξ(x, u) + O(a 2 ), u = g(x, u; a) = u + aη(x, u) + O(a 2 ). is a variational symmetry of R σ [u] if it preserves the action J[U]. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
26 Variational Symmetries Consider a general DE system R σ [u] = R σ (x, u, u,..., k u) = 0, σ = 1,..., N that follows from a variational principle with J[U] = Ω L(x, U, U,..., k U)dx. Definition A symmetry of R σ [u] given by x = f (x, u; a) = x + aξ(x, u) + O(a 2 ), u = g(x, u; a) = u + aη(x, u) + O(a 2 ). is a variational symmetry of R σ [u] if it preserves the action J[U]. Example 1: translations for the wave equation u tt = (T /ρ) 2 u xx, L = 1 2 ρut T ux 2. The translation x = x + C, t = t, u = u is evidently a variational symmetry. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
27 Variational Symmetries Consider a general DE system R σ [u] = R σ (x, u, u,..., k u) = 0, σ = 1,..., N that follows from a variational principle with J[U] = Ω L(x, U, U,..., k U)dx. Definition A symmetry of R σ [u] given by x = f (x, u; a) = x + aξ(x, u) + O(a 2 ), u = g(x, u; a) = u + aη(x, u) + O(a 2 ). is a variational symmetry of R σ [u] if it preserves the action J[U]. Example 2: scaling for the wave equation u tt = (T /ρ) 2 u xx, L = 1 2 ρut T ux 2. The scaling x = x, t = t, u = u/α is not a variational symmetry: L = α 2 L, J = α 2 J. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
28 Noether s Theorem (restricted to point symmetries) Theorem Given: a PDE system R σ [u] = R σ (x, u, u,..., k u) = 0, σ = 1,..., N, following from a variational principle; a variational symmetry (x i ) = f i (x, u; a) = x i + aξ i (x, u) + O(a 2 ), (u σ ) = g σ (x, u; a) = u σ + aη σ (x, u) + O(a 2 ). Then the system R σ [u] has a conservation law D i Φ i [u] = 0. In particular, D i Φ i [u] Λ σ[u]r σ [u] = 0, where the multipliers are given by Λ σ = η σ (x, u) uσ x i ξ i (x, u). A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
29 Noether s Theorem: Examples Example: translation symmetry for the harmonic oscillator Equation: ẍ(t) + ω 2 x(t) = 0, ω 2 = k/m. Symmetry: t = t + a, ξ = 1; x = x, η = 0, Multiplier (integrating factor): Λ = η ẋ(t)ξ = ẋ; Conservation law: ΛR = ẋ(ẍ(t) + ω 2 x(t)) = 1 m ( d mẋ 2 (t) + kx ) 2 (t) = 0. dt 2 2 A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
30 Limitations of Noether s Theorem The given DE system, as written, must be variational. Numbers of PDEs and dependent variables must coincide. Dissipative systems are not variational. If single PDE, must be of even order. If a PDE system is not variational, artifices sometimes can make it variational! A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
31 Limitations of Noether s Theorem The given DE system, as written, must be variational. Numbers of PDEs and dependent variables must coincide. Dissipative systems are not variational. If single PDE, must be of even order. If a PDE system is not variational, artifices sometimes can make it variational! Example 1: The use of multipliers. The PDE u tt + H (u x)u xx + H(u x) = 0, as written, does not admit a variational principle. However, the equivalent PDE e x [u tt + H (u x)u xx + H(u x)] = 0, does follow from a variational principle! A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
32 Limitations of Noether s Theorem The given DE system, as written, must be variational. Numbers of PDEs and dependent variables must coincide. Dissipative systems are not variational. If single PDE, must be of even order. If a PDE system is not variational, artifices sometimes can make it variational! Example 2: The use of a transformation. The PDE e x u tt e 3x (u + u x) 2 (u + 2u x + u xx) = 0, as written, does not admit a variational principle. But the point transformation x = x, t = t, u (x, t ) = y(x, t) = e x u(x, t), maps it into the self-adjoint PDE y tt (y x) 2 y xx = 0, which is the Euler Lagrange equation for the Lagrangian L[Y ] = y 2 t /2 y 4 x /12. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
33 Limitations of Noether s Theorem The given DE system, as written, must be variational. Numbers of PDEs and dependent variables must coincide. Dissipative systems are not variational. If single PDE, must be of even order. If a PDE system is not variational, artifices sometimes can make it variational! Example 3: The use of a differential substitution. The KdV equation u t + uu x + u xxx = 0, as written, does not admit a variational principle. However, a substitution u = v x yields a variational PDE v xt + v xv xx + v xxxx = 0. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
34 Limitations of Noether s Theorem (ctd.)... The use of an artificial additional equation. Any PDE system can be made variational, by appending an adjoint equation! Example: The diffusion equation u t u xx = 0 is dissipative, hence not self-adjoint. Its adjoint equation: w t + w xx = 0. But the PDE system u t u xx = 0, is self-adjoint! Lagrangian: ũ t + ũ xx = 0 L = u tw w tu + 2u xw x. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
35 Outline 1 Conservation Laws Definition and Examples Applications 2 Variational Principles 3 Symmetries and Noether s Theorem Symmetries and Variational Symmetries Noether s Theorem, Examples Limitations of Noether s Theorem 4 Direct Construction Method for Conservation Laws The Idea Completeness A Detailed Example 5 Examples of Construction of Conservation Laws Symbolic Software (Maple) KdV Surfactant Dynamics Equations 6 Conclusions A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
36 The Idea of the Direct Construction Method Definition The Euler operator with respect to U j : E U j = D U j i U j i + + ( 1) s D i1... D is U j i 1...i s +, j = 1,..., m. Consider a general DE system R σ [u] = R σ (x, u, u,..., k u) = 0, variables x = (x 1,..., x n ), u = u(x) = (u 1,..., u m ). σ = 1,..., N with Theorem The equations E U j F (x, U, U,..., s U) 0, hold for arbitrary U(x) if and only if j = 1,..., m for some functions Ψ i (x, U,...). F (x, U, U,..., s U) D i Ψ i A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
37 The Idea of the Direct Construction Method Consider a general DE system R σ [u] = R σ (x, u, u,..., k u) = 0, σ = 1,..., N with variables x = (x 1,..., x n ), u = u(x) = (u 1,..., u m ). Direct Construction Method Specify dependence of multipliers: Λ σ = Λ σ(x, U,...), σ = 1,..., N. Solve the set of determining equations E U j (Λ σr σ ) 0, j = 1,..., m, for arbitrary U(x) (off of solution set!) to find all such sets of multipliers. Find the corresponding fluxes Φ i (x, U,...) satisfying the identity Λ σr σ D i Φ i. Each set of fluxes and multipliers yields a local conservation law D i Φ i (x, u,...) = 0, holding for all solutions u(x) of the given PDE system. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
38 Completeness Completeness of the Direct Construction Method For the majority of physical DE systems (in particular, all systems in solved form), all conservation laws follow from linear combinations of equations! Λ σr σ D i Φ i. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
39 A Detailed Example Consider a nonlinear telegraph system for u 1 = u(x, t), u 2 = v(x, t): R 1 [u, v] = v t (u 2 + 1)u x u = 0, R 2 [u, v] = u t v x = 0. Multiplier ansatz: Λ 1 = ξ(x, t, U, V ), Λ 2 = φ(x, t, U, V ). A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
40 A Detailed Example Consider a nonlinear telegraph system for u 1 = u(x, t), u 2 = v(x, t): R 1 [u, v] = v t (u 2 + 1)u x u = 0, R 2 [u, v] = u t v x = 0. Multiplier ansatz: Λ 1 = ξ(x, t, U, V ), Λ 2 = φ(x, t, U, V ). Determining equations: [ E U ξ(x, t, U, V )(Vt (U 2 + 1)U x U) + φ(x, t, U, V )(U t V ] x) 0, [ E V ξ(x, t, U, V )(Vt (U 2 + 1)U x U) + φ(x, t, U, V )(U t V ] x) 0. Euler operators: E U = U Dx D t, U x U t E V = V Dx D t. V x V t A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
41 A Detailed Example Consider a nonlinear telegraph system for u 1 = u(x, t), u 2 = v(x, t): R 1 [u, v] = v t (u 2 + 1)u x u = 0, R 2 [u, v] = u t v x = 0. Multiplier ansatz: Λ 1 = ξ(x, t, U, V ), Λ 2 = φ(x, t, U, V ). Determining equations: [ E U ξ(x, t, U, V )(Vt (U 2 + 1)U x U) + φ(x, t, U, V )(U t V ] x) 0, [ E V ξ(x, t, U, V )(Vt (U 2 + 1)U x U) + φ(x, t, U, V )(U t V ] x) 0. Split determining equations: φ V ξ U = 0, φ U (U 2 + 1)ξ V = 0, φ x ξ t Uξ V = 0, (U 2 + 1)ξ x φ t Uξ U ξ = 0. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
42 A Detailed Example Consider a nonlinear telegraph system for u 1 = u(x, t), u 2 = v(x, t): R 1 [u, v] = v t (u 2 + 1)u x u = 0, R 2 [u, v] = u t v x = 0. Multiplier ansatz: Λ 1 = ξ(x, t, U, V ), Λ 2 = φ(x, t, U, V ). Solution: five sets of multipliers (ξ, φ) = 0 1 t x 1 2 t2 1 t e x+ 1 2 U2 +V e x+ 1 2 U2 V Ue x+ 1 2 U2 +V Ue x+ 1 2 U2 V A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
43 A Detailed Example Consider a nonlinear telegraph system for u 1 = u(x, t), u 2 = v(x, t): R 1 [u, v] = v t (u 2 + 1)u x u = 0, R 2 [u, v] = u t v x = 0. Multiplier ansatz: Λ 1 = ξ(x, t, U, V ), Λ 2 = φ(x, t, U, V ). Resulting five conservation laws: D tu D xv = 0, D t[(x 1 2 t2 )u + tv] + D x[( 1 2 t2 x)v t( 1 3 u3 + u)] = 0, D t[v tu] + D x[tv ( 1 3 u3 + u)] = 0, D t[e x+ 2 1 u2 +v ] + D x[ ue x+ 1 2 u2 +v ] = 0, D t[e x+ 1 2 u2 v ] + D x[ue x+ 1 2 u2 v ] = 0. To obtain further conservation laws, extend the multiplier ansatz... A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
44 Outline 1 Conservation Laws Definition and Examples Applications 2 Variational Principles 3 Symmetries and Noether s Theorem Symmetries and Variational Symmetries Noether s Theorem, Examples Limitations of Noether s Theorem 4 Direct Construction Method for Conservation Laws The Idea Completeness A Detailed Example 5 Examples of Construction of Conservation Laws Symbolic Software (Maple) KdV Surfactant Dynamics Equations 6 Conclusions A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
45 Symbolic Software for Computation of Conservation Laws Example of use of the GeM package for Maple for the KdV. Use the module: with(gem): Declare variables: gem_decl_vars(indeps=[x,t], deps=[u(x,t)]); Declare the equation: gem_decl_eqs([diff(u(x,t),t)=u(x,t)*diff(u(x,t),x) +diff(u(x,t),x,x,x)], solve_for=[diff(u(x,t),t)]); Generate determining equations: det_eqs:=gem_conslaw_det_eqs([x,t, U(x,t), diff(u(x,t),x), diff(u(x,t),x,x), diff(u(x,t),x,x,x)]): Reduce the overdetermined system: CL_multipliers:=gem_conslaw_multipliers(); simplified_eqs:=detools[rifsimp](det_eqs, CL_multipliers, mindim=1); A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
46 Symbolic Software for Computation of Conservation Laws Example of use of the GeM package for Maple for the KdV. Solve determining equations: multipliers_sol:=pdsolve(simplified_eqs[solved]); Obtain corresponding conservation law fluxes/densities: gem_get_cl_fluxes(multipliers_sol, method=*****); A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
47 Examples of Symbolic Computation of Conservation Laws Example 1 KdV equation: u t + uu x + u xxx = 0. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
48 Example 2: Conserved Form of Surfactant Dynamics Equations Ref.: C. Kallendorf, A.S., M. Oberlack, Y.Wang, 2011 Surfactants = surface active agents. Two-phase interface described by a level set function Φ; S = { x R 3 : Φ(x) = 0 }. Unit normal: n = Φ. Concentration: c. Φ A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
49 Example 2: Conserved Form of Surfactant Dynamics Equations Ref.: C. Kallendorf, A.S., M. Oberlack, Y.Wang, 2011 Surfactants = surface active agents. Two-phase interface described by a level set function Φ; S = { x R 3 : Φ(x) = 0 }. Unit normal: n = Φ. Concentration: c. Φ Surfactant Dynamics Equations: Incompressibility of the flow: Interface transport by the flow: u = 0. Surfactant transport: Φ t + u (Φ) = 0. c t + u i c x i cn i n j u i x j α(δ ij n i n j ) x j ( (δ ik n i n k ) c ) = 0. x k A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
50 Example 2: Conserved Form of Surfactant Dynamics Equations Result The surfactant dynamics system can be written in a fully conserved form: where u = 0. Φ t + u (Φ) = 0. (c F(Φ) Φ ) + ( ) A i F(Φ) Φ t x i A i = cu i α and F is an arbitrary sufficiently smooth function. = 0, ( (δ ik n i n k ) c ), i = 1, 2, 3, x k A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
51 Conclusions and Open Problems Conclusions Divergence-type conservation laws are useful in analysis and numerics. Conservation laws can be obtained systematically through the Direct Construction Method, which employs multipliers and Euler operators. The method is implemented in a symbolic package GeM for Maple. For variational DE systems, conservation laws correspond to variational symmetries. Noether s theorem is usually not a preferred way to derive unknown conservation laws. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
52 Conclusions and Open Problems Conclusions Divergence-type conservation laws are useful in analysis and numerics. Conservation laws can be obtained systematically through the Direct Construction Method, which employs multipliers and Euler operators. The method is implemented in a symbolic package GeM for Maple. For variational DE systems, conservation laws correspond to variational symmetries. Noether s theorem is usually not a preferred way to derive unknown conservation laws. Open problems Computations become hard for conservation laws of high orders. Can we tell anything in advance about highest order of the conservation law for a given DE system? For complicated DEs and multipliers, computation of fluxes/densities can become challenging. A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
53 Some references Bluman, G.W., Cheviakov, A.F., and Anco, S.C. (2010). Applications of Symmetry Methods to Partial Differential Equations. Springer: Applied Mathematical Sciences, Vol Anco, S.C. and Bluman, G.W. (1997). Direct construction of conservation laws from field equations. Phys. Rev. Lett. 78, Anco, S.C. and Bluman, G.W. (2002). Direct construction method for conservation laws of partial differential equations. Part I: Examples of conservation law classifications. Eur. J. Appl. Math. 13, Cheviakov, A.F. (2007). GeM software package for computation of symmetries and conservation laws of differential equations. Comput. Phys. Commun. 176, Kallendorf, C., Cheviakov, A.F., Oberlack, M., and Wang, Y. Conservation Laws of Surfactant Transport Equations. Submitted, A. Shevyakov (Math&Stat) Conservation Laws 10 Nov / 33
Alexei 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 informationConservation Laws of Surfactant Transport Equations
Conservation Laws of Surfactant Transport Equations Alexei Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada Winter 2011 CMS Meeting Dec. 10, 2011 A. Cheviakov
More informationConservation Laws for Nonlinear Equations: Theory, Computation, and Examples
Conservation Laws for Nonlinear Equations: Theory, Computation, and Examples Alexei Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada Seminar February 2015
More informationConstruction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem
Proceedings of 4th Workshop Group Analysis of Differential Equations & Integrability 2009, 1 23 Construction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem George W. BLUMAN,
More informationA Recursion Formula for the Construction of Local Conservation Laws of Differential Equations
A Recursion Formula for the Construction of Local Conservation Laws of Differential Equations Alexei Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada December
More informationNatural States and Symmetry Properties of. Two-Dimensional Ciarlet-Mooney-Rivlin. Nonlinear Constitutive Models
Natural States and Symmetry Properties of Two-Dimensional Ciarlet-Mooney-Rivlin Nonlinear Constitutive Models Alexei Cheviakov, Department of Mathematics and Statistics, Univ. Saskatchewan, Canada Jean-François
More informationSymmetry Methods for Differential Equations and Conservation Laws. Peter J. Olver University of Minnesota
Symmetry Methods for Differential Equations and Conservation Laws Peter J. Olver University of Minnesota http://www.math.umn.edu/ olver Santiago, November, 2010 Symmetry Groups of Differential Equations
More informationSymbolic Computation of Nonlocal Symmetries and Nonlocal Conservation Laws of Partial Differential Equations Using the GeM Package for Maple
Symbolic Computation of Nonlocal Symmetries and Nonlocal Conservation Laws of Partial Differential Equations Using the GeM Package for Maple Alexei F. Cheviakov Abstract The use of the symbolic software
More informationA symmetry-based method for constructing nonlocally related partial differential equation systems
A symmetry-based method for constructing nonlocally related partial differential equation systems George W. Bluman and Zhengzheng Yang Citation: Journal of Mathematical Physics 54, 093504 (2013); doi:
More informationA comparison of conservation law construction approaches for the two-dimensional incompressible Mooney-Rivlin hyperelasticity model
JOURNAL OF MATHEMATICAL PHYSICS 56, 505 (05) A comparison of conservation law construction approaches for the two-dimensional incompressible Mooney-Rivlin hyperelasticity model A. F. Cheviakov and S. St.
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 informationConservation Laws of Fluid Dynamics Models
Conservation Laws of Fluid Dynamics Models Prof. Alexei Cheviakov (Alt. English spelling: Alexey Shevyakov) Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada June
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 informationNew methods of reduction for ordinary differential equations
IMA Journal of Applied Mathematics (2001) 66, 111 125 New methods of reduction for ordinary differential equations C. MURIEL AND J. L. ROMERO Departamento de Matemáticas, Universidad de Cádiz, PO Box 40,
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 informationApplications of Symmetries and Conservation Laws to the Study of Nonlinear Elasticity Equations
Applications of Symmetries and Conservation Laws to the Study of Nonlinear Elasticity Equations A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment of the Requirements
More informationB.7 Lie Groups and Differential Equations
96 B.7. LIE GROUPS AND DIFFERENTIAL EQUATIONS B.7 Lie Groups and Differential Equations Peter J. Olver in Minneapolis, MN (U.S.A.) mailto:olver@ima.umn.edu The applications of Lie groups to solve differential
More informationOn the Linearization of Second-Order Dif ferential and Dif ference Equations
Symmetry, Integrability and Geometry: Methods and Applications Vol. (006), Paper 065, 15 pages On the Linearization of Second-Order Dif ferential and Dif ference Equations Vladimir DORODNITSYN Keldysh
More informationGroup Actions and Cohomology in the Calculus of Variations
Group Actions and Cohomology in the Calculus of Variations JUHA POHJANPELTO Oregon State and Aalto Universities Focused Research Workshop on Exterior Differential Systems and Lie Theory Fields Institute,
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
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 informationConservation laws for the geodesic equations of the canonical connection on Lie groups in dimensions two and three
Appl Math Inf Sci 7 No 1 311-318 (013) 311 Applied Mathematics & Information Sciences An International Journal Conservation laws for the geodesic equations of the canonical connection on Lie groups in
More informationDirect construction method for conservation laws of partial differential equations Part I: Examples of conservation law classifications
Euro. Jnl of Applied Mathematics (2002), vol. 13, pp. 545 566. c 2002 Cambridge University Press DOI: 10.1017/S095679250100465X Printed in the United Kingdom 545 Direct construction method for conservation
More informationSymmetry Classification of KdV-Type Nonlinear Evolution Equations
Proceedings of Institute of Mathematics of NAS of Ukraine 2004, Vol. 50, Part 1, 125 130 Symmetry Classification of KdV-Type Nonlinear Evolution Equations Faruk GÜNGÖR, Victor LAHNO and Renat ZHDANOV Department
More informationSymbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multiple Space Dimensions
Symbolic Computation of Conservation Laws of Nonlinear Partial Differential Equations in Multiple Space Dimensions Willy Hereman and Douglas Poole Department of Mathematical and Computer Sciences Colorado
More informationOn Reduction and Q-conditional (Nonclassical) Symmetry
Symmetry in Nonlinear Mathematical Physics 1997, V.2, 437 443. On Reduction and Q-conditional (Nonclassical) Symmetry Roman POPOVYCH Institute of Mathematics of the National Academy of Sciences of Ukraine,
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More 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 informationDirect and inverse problems on conservation laws
Direct and inverse problems on conservation laws Roman O. Popovych Wolfgang Pauli Institute, Vienna, Austria & Institute of Mathematics of NAS of Ukraine, Kyiv, Ukraine joint work with Alexander Bihlo,
More informationA Note on Nonclassical Symmetries of a Class of Nonlinear Partial Differential Equations and Compatibility
Commun. Theor. Phys. (Beijing, China) 52 (2009) pp. 398 402 c Chinese Physical Society and IOP Publishing Ltd Vol. 52, No. 3, September 15, 2009 A Note on Nonclassical Symmetries of a Class of Nonlinear
More informationSymmetries and reduction techniques for dissipative models
Symmetries and reduction techniques for dissipative models M. Ruggieri and A. Valenti Dipartimento di Matematica e Informatica Università di Catania viale A. Doria 6, 95125 Catania, Italy Fourth Workshop
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 informationA NONLINEAR GENERALIZATION OF THE CAMASSA-HOLM EQUATION WITH PEAKON SOLUTIONS
A NONLINEAR GENERALIZATION OF THE CAMASSA-HOLM EQUATION WITH PEAKON SOLUTIONS STEPHEN C. ANCO 1, ELENA RECIO 1,2, MARÍA L. GANDARIAS2, MARÍA S. BRUZÓN2 1 department of mathematics and statistics brock
More informationResearch Article Equivalent Lagrangians: Generalization, Transformation Maps, and Applications
Journal of Applied Mathematics Volume 01, Article ID 86048, 19 pages doi:10.1155/01/86048 Research Article Equivalent Lagrangians: Generalization, Transformation Maps, and Applications N. Wilson and A.
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 information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
More informationNonlinear Equations for Finite-Amplitude Wave Propagation in Fiber-Reinforced Hyperelastic Media
Nonlinear Equations for Finite-Amplitude Wave Propagation in Fiber-Reinforced Hyperelastic Media Alexei F. Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada
More informationNonlocally related PDE systems for one-dimensional nonlinear elastodynamics
J Eng Math (2008) 62:203 22 DOI 0.007/s0665-008-922-7 Nonlocally related PDE systems for one-dimensional nonlinear elastodynamics G. Bluman A. F. Cheviakov J.-F. Ganghoffer Received: 5 July 2006 / Accepted:
More informationGroup classification of nonlinear wave equations
JOURNAL OF MATHEMATICAL PHYSICS 46, 053301 2005 Group classification of nonlinear wave equations V. Lahno a State Pedagogical University, 36000 Poltava, Ukraine R. Zhdanov b Institute of Mathematics of
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 informationNonlocal Symmetry and Generating Solutions for the Inhomogeneous Burgers Equation
Proceedings of Institute of Mathematics of NAS of Ukraine 004, Vol. 50, Part, 77 8 Nonlocal Symmetry and Generating Solutions for the Inhomogeneous Burgers Equation Valentyn TYCHYNIN and Olexandr RASIN
More informationDifferential characteristic set algorithm for PDEs symmetry computation, classification, decision and extension
Symposium in honor of George Bluman UBC May 12 16, 2014 Differential characteristic set algorithm for PDEs symmetry computation, classification, decision and extension Chaolu Shanghai Maritime University
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationSymbolic Computation and New Soliton-Like Solutions of the 1+2D Calogero-Bogoyavlenskii-Schif Equation
MM Research Preprints, 85 93 MMRC, AMSS, Academia Sinica, Beijing No., December 003 85 Symbolic Computation and New Soliton-Like Solutions of the 1+D Calogero-Bogoyavlenskii-Schif Equation Zhenya Yan Key
More informationPartial Differential Equations, Winter 2015
Partial Differential Equations, Winter 215 Homework #2 Due: Thursday, February 12th, 215 1. (Chapter 2.1) Solve u xx + u xt 2u tt =, u(x, ) = φ(x), u t (x, ) = ψ(x). Hint: Factor the operator as we did
More informationBenjamin Bona Mahony Equation with Variable Coefficients: Conservation Laws
Symmetry 2014, 6, 1026-1036; doi:10.3390/sym6041026 OPEN ACCESS symmetry ISSN 2073-8994 www.mdpi.com/journal/symmetry Article Benjamin Bona Mahony Equation with Variable Coefficients: Conservation Laws
More informationx ct x + t , and the characteristics for the associated transport equation would be given by the solution of the ode dx dt = 1 4. ξ = x + t 4.
. The solution is ( 2 e x+ct + e x ct) + 2c x+ct x ct sin(s)dx ( e x+ct + e x ct) + ( cos(x + ct) + cos(x ct)) 2 2c 2. To solve the PDE u xx 3u xt 4u tt =, you can first fact the differential operat to
More informationSYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS. Willy Hereman
. SYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS Willy Hereman Dept. of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado
More informationSuperintegrability? Hidden linearity? Classical quantization? Symmetries and more symmetries!
Superintegrability? Hidden linearity? Classical quantization? Symmetries and more symmetries! Maria Clara Nucci University of Perugia & INFN-Perugia, Italy Conference on Nonlinear Mathematical Physics:
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationNonlinear Wave Propagation in 1D Random Media
Nonlinear Wave Propagation in 1D Random Media J. B. Thoo, Yuba College UCI Computational and Applied Math Seminar Winter 22 Typeset by FoilTEX Background O Doherty and Anstey (Geophysical Prospecting,
More informationTraffic Flow Problems
Traffic Flow Problems Nicodemus Banagaaya Supervisor : Dr. J.H.M. ten Thije Boonkkamp October 15, 2009 Outline Introduction Mathematical model derivation Godunov Scheme for the Greenberg Traffic model.
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More informationOn second order sufficient optimality conditions for quasilinear elliptic boundary control problems
On second order sufficient optimality conditions for quasilinear elliptic boundary control problems Vili Dhamo Technische Universität Berlin Joint work with Eduardo Casas Workshop on PDE Constrained Optimization
More informationGeneralized evolutionary equations and their invariant solutions
Generalized evolutionary equations and their invariant solutions Rodica Cimpoiasu, Radu Constantinescu University of Craiova, 13 A.I.Cuza Street, 200585 Craiova, Dolj, Romania Abstract The paper appplies
More informationÉlie Cartan s Theory of Moving Frames
Élie Cartan s Theory of Moving Frames Orn Arnaldsson Department of Mathematics University of Minnesota, Minneapolis Special Topics Seminar, Spring 2014 The Story of Symmetry Felix Klein and Sophus Lie
More information1.5 First Order PDEs and Method of Characteristics
1.5. FIRST ORDER PDES AND METHOD OF CHARACTERISTICS 35 1.5 First Order PDEs and Method of Characteristics We finish this introductory chapter by discussing the solutions of some first order PDEs, more
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
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 informationPotential Symmetries and Differential Forms. for Wave Dissipation Equation
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 42, 2061-2066 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2013.36163 Potential Symmetries and Differential Forms for Wave Dissipation
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 informationPARTIAL DIFFERENTIAL EQUATIONS. Lecturer: D.M.A. Stuart MT 2007
PARTIAL DIFFERENTIAL EQUATIONS Lecturer: D.M.A. Stuart MT 2007 In addition to the sets of lecture notes written by previous lecturers ([1, 2]) the books [4, 7] are very good for the PDE topics in the course.
More informationA symmetry analysis of Richard s equation describing flow in porous media. Ron Wiltshire
BULLETIN OF THE GREEK MATHEMATICAL SOCIETY Volume 51, 005 93 103) A symmetry analysis of Richard s equation describing flow in porous media Ron Wiltshire Abstract A summary of the classical Lie method
More informationexamples of equations: what and why intrinsic view, physical origin, probability, geometry
Lecture 1 Introduction examples of equations: what and why intrinsic view, physical origin, probability, geometry Intrinsic/abstract F ( x, Du, D u, D 3 u, = 0 Recall algebraic equations such as linear
More informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
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 informationSymmetry and Exact Solutions of (2+1)-Dimensional Generalized Sasa Satsuma Equation via a Modified Direct Method
Commun. Theor. Phys. Beijing, China 51 2009 pp. 97 978 c Chinese Physical Society and IOP Publishing Ltd Vol. 51, No., June 15, 2009 Symmetry and Exact Solutions of 2+1-Dimensional Generalized Sasa Satsuma
More informationTravelling wave solutions for a CBS equation in dimensions
AMERICAN CONFERENCE ON APPLIED MATHEMATICS (MATH '8), Harvard, Massachusetts, USA, March -6, 8 Travelling wave solutions for a CBS equation in + dimensions MARIA LUZ GANDARIAS University of Cádiz Department
More informationConnecting Euler and Lagrange systems as nonlocally related systems of dynamical nonlinear elasticity
Arch. Mech., 63, 4, pp. 1 20, Warszawa 2011 Connecting Euler and Lagrange systems as nonlocally related systems of dynamical nonlinear elasticity G. BLUMAN 1), J. F. GANGHOFFER 2) 1) Mathematics Department
More informationPartial Differential Equations Separation of Variables. 1 Partial Differential Equations and Operators
PDE-SEP-HEAT-1 Partial Differential Equations Separation of Variables 1 Partial Differential Equations and Operators et C = C(R 2 ) be the collection of infinitely differentiable functions from the plane
More informationMultidimensional partial differential equation systems: Nonlocal symmetries, nonlocal conservation laws, exact solutions
JOURNAL OF MATHEMATICAL PHYSICS 51, 103522 2010 Multidimensional partial differential equation systems: Nonlocal symmetries, nonlocal conservation laws, exact solutions Alexei F. Cheviakov 1,a and George
More informationThe first three (of infinitely many) conservation laws for (1) are (3) (4) D t (u) =D x (3u 2 + u 2x ); D t (u 2 )=D x (4u 3 u 2 x +2uu 2x ); D t (u 3
Invariants and Symmetries for Partial Differential Equations and Lattices Λ Ünal Göktaοs y Willy Hereman y Abstract Methods for the computation of invariants and symmetries of nonlinear evolution, wave,
More informationANALYSIS OF A NONLINEAR SURFACE WIND WAVES MODEL VIA LIE GROUP METHOD
Electronic Journal of Differential Equations, Vol. 206 (206), No. 228, pp. 8. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ANALYSIS OF A NONLINEAR SURFACE WIND WAVES MODEL
More informationAlgorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations
dmitry.lyakhov@kaust.edu.sa Symbolic Computations May 4, 2018 1 / 25 Algorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations Dmitry A. Lyakhov 1 1 Computational Sciences
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 informationPotential symmetry and invariant solutions of Fokker-Planck equation in. cylindrical coordinates related to magnetic field diffusion
Potential symmetry and invariant solutions of Fokker-Planck equation in cylindrical coordinates related to magnetic field diffusion in magnetohydrodynamics including the Hall current A. H. Khater a,1,
More informationarxiv: v1 [math-ph] 25 Jul Preliminaries
arxiv:1107.4877v1 [math-ph] 25 Jul 2011 Nonlinear self-adjointness and conservation laws Nail H. Ibragimov Department of Mathematics and Science, Blekinge Institute of Technology, 371 79 Karlskrona, Sweden
More informationConditional Symmetry Reduction and Invariant Solutions of Nonlinear Wave Equations
Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part 1, 229 233 Conditional Symmetry Reduction and Invariant Solutions of Nonlinear Wave Equations Ivan M. TSYFRA Institute of Geophysics
More informationSymmetry Properties of Two-Dimensional Ciarlet-Mooney-Rivlin Constitutive Models in Nonlinear Elastodynamics
Symmetry Properties of Two-Dimensional Ciarlet-Mooney-Rivlin Constitutive Models in Nonlinear Elastodynamics A. F. Cheviakov a, J.-F. Ganghoffer b a Department of Mathematics and Statistics, University
More informationApplications of Lie Group Analysis to the Equations of Motion of Inclined Unsagged Cables
Applied Mathematical Sciences, Vol., 008, no. 46, 59-69 Applications of Lie Group Analysis to the Equations of Motion of Inclined Unsagged Cables Waraporn Chatanin Department of Mathematics King Mongkut
More informationFree energy concept Free energy approach LBM implementation Parameters
BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH
More informationConditional symmetries of the equations of mathematical physics
W.I. Fushchych, Scientific Works 2003, Vol. 5, 9 16. Conditional symmetries of the equations of mathematical physics W.I. FUSHCHYCH We briefly present the results of research in conditional symmetries
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 informationComplete Description of Turbulence in Terms of Hopf Functional and LMN Hierarchy: New Symmetries and Invariant Solutions
Complete Description of Turbulence in Terms of Hopf Functional and LMN Hierarchy: New Symmetries and Invariant Solutions Marta Wacławczyk Abstract This paper deals with two methods for the full statistical
More informationLie and Non-Lie Symmetries of Nonlinear Diffusion Equations with Convection Term
Symmetry in Nonlinear Mathematical Physics 1997, V.2, 444 449. Lie and Non-Lie Symmetries of Nonlinear Diffusion Equations with Convection Term Roman CHERNIHA and Mykola SEROV Institute of Mathematics
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 information2.3 Calculus of variations
2.3 Calculus of variations 2.3.1 Euler-Lagrange equation The action functional S[x(t)] = L(x, ẋ, t) dt (2.3.1) which maps a curve x(t) to a number, can be expanded in a Taylor series { S[x(t) + δx(t)]
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 informationMATH1231 CALCULUS. Session II Dr John Roberts (based on notes of A./Prof. Bruce Henry) Red Center Room 3065
MATH1231 CALCULUS Session II 2007. Dr John Roberts (based on notes of A./Prof. Bruce Henry) Red Center Room 3065 Jag.Roberts@unsw.edu.au MATH1231 CALCULUS p.1/66 Overview Systematic Integration Techniques
More informationMATH 819 FALL We considered solutions of this equation on the domain Ū, where
MATH 89 FALL. The D linear wave equation weak solutions We have considered the initial value problem for the wave equation in one space dimension: (a) (b) (c) u tt u xx = f(x, t) u(x, ) = g(x), u t (x,
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 informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationMath 5440 Problem Set 3 Solutions
Mth 544 Mth 544 Problem Set 3 Solutions Aron Fogelson Fll, 25 1: Logn, 1.5 # 2) Repet the derivtion for the eqution of motion of vibrting string when, in ddition, the verticl motion is retrded by dmping
More informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
More informationSymmetry classification of KdV-type nonlinear evolution equations
arxiv:nlin/0201063v2 [nlin.si] 16 Sep 2003 Symmetry classification of KdV-type nonlinear evolution equations F. Güngör Department of Mathematics, Faculty of Science and Letters, Istanbul Technical University,
More informationMA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation
MA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation Lecture 12 Lecture 12 MA 201, PDE (2016) 1 / 24 Formal Solution of
More informationPoint-transformation structures on classes of differential equations
Point-transformation structures on classes of differential equations Roman O. Popovych Wolfgang Pauli Institute, Vienna, Austria & Institute of Mathematics of NAS of Ukraine, Kyiv, Ukraine joint work with
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More information