Numerical Solutions to Partial Differential Equations
|
|
- Allyson Quinn
- 6 years ago
- Views:
Transcription
1 Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University
2 Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear Hyperbolic Partial Differential Equation n 1 Scalar case (u R 1 ), u t + a i u xi + b u = ψ 0, i=1 where a i, b, ψ 0 are real functions of x = (x 1,..., x n ) and t. 2 Vector case (u = (u 1,, u p ) T R p ), n u t + A i u xi + B u = ψ 0, i=1 where A i, B R p p, ψ 0 R p are real functions of t and x = (x 1,..., x n ), and α R n, A(x, t) = n i=1 α ia i (x, t) is real diagonalizable, i.e. A(x, t) has p linearly independent eigenvectors corresponding to real eigenvalues. 3 If A(x, t) = n i=1 α ia i (x, t) has p mutually different real eigenvalues, the system is called strictly hyperbolic. 2 / 40
3 Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 2nd Order Linear Hyperbolic Partial Differential Equation A general 2nd order scalar equation (u R 1 ), n n n u tt + 2 a i u xi t + b 0 u t a i u xi x + b i u xi + cu = ψ 0, i=1 i,=1 where a i, a i, b i, c and ψ 0 are real functions x = (x 1,..., x n ) and t, (a i ) is real symmetric positive definite. Define v = u, v 0 = u t, v i = u xi, then the above 2nd order scalar equation transforms into a first order linear system of partial differential equations for v = (v, v 0, v 1,, v n ) T n Av t + A i v xi + Bv = ψ 0, i=1 i=1 3 / 40
4 Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 2nd Order Scalar Transforms to n-d 1st Order System (p = n + 2) a i a 1i a ni A = 0 0 a 11 a 1n, A i = 0 a 1i 0 0, a n1 a nn 0 a ni c b 0 b 1 b n ψ 0 B = , ψ 0 = / 40
5 Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 2nd Order Scalar Transforms to n-d 1st Order System (p = n + 2) 1 Let R T AR = I (A is real symmetric positive definite); 2 Introduce a new variable w = R 1 v; 3 Denote Âi = R T A i R, ˆB = R T BR, ˆψ 0 = R T ψ; 4 Â(x, t) = n i=1 α iâi(x, t) is a real symmetric matrix and thus is real diagonizable for all real α i, i = 1,..., n; 5 The 2nd order scalar equation now transforms into a 1st order linear hyperbolic system of partial differential equations for w R (n+2) : n w t + Â i w xi + ˆBw = ˆψ 0. i=1 5 / 40
6 Introduction to Hyperbolic Equations The Hyperbolic Equations Standard Form of n-d 1st Order Linear Hyperbolic Equations 1 The standard form of 1st order linear hyperbolic equation: n u t + a i u xi = ψ, (ψ = ψ 0 b u). i=1 2 The standard form of 1st order linear hyperbolic system: n u t + A i u xi = ψ, (ψ = ψ 0 B u); i=1 3 The equation (system) is said to be homogeneous, if ψ = 0; 6 / 40
7 Introduction to Hyperbolic Equations The Hyperbolic Equations Standard Form of n-d 1st Order Linear Hyperbolic Equations 4 In general, a higher order linear hyperbolic equation (system of equations) can always be transformed into a first order linear hyperbolic system of equations. 5 An equation (system) is said to be nonlinear, if at least one of the coefficients depends on the unknown or the right hand side term is a nonlinear function of the unknown. 7 / 40
8 Introduction to Hyperbolic Equations Example on Hyperbolic Equations and Balance Laws An Example of a Balance Law Substance balance in flowing fluid in a 1D pipe. 1 u(x, t): substance density, measured by mass per unit length. 2 f (x, t, u): mass flux, measured by mass per unit time. 3 ψ(x, t, u): the density of the mass source (or sink), measured by mass per unit length per unit time. 4 Balance law (integral form), for given x l < x r and t b < t a, xr x l u(x, t a ) dx = ta xr t b f (x r, t, u(x r, t)) dt + x l u(x, t b ) dx + ta xr t b x l ta t b f (x l, t, u(x l, t)) dt ψ(x, t, u(x, t)) dx dt. 8 / 40
9 Introduction to Hyperbolic Equations Example on Hyperbolic Equations and Balance Laws An Example of a Balance Law 5 Suppose u and f are sufficiently smooth, we are led to ta xr t b x l [u(x, t) t + f (x, t, u(x, t)) x ψ(x, t, u(x, t))] dx dt = 0, for all given x l < x r and t b < t a ; 6 or equivalently: u(x, t) t + f (x, t, u(x, t)) x = ψ(x, t, u(x, t)), this is the balance law in differential form. 9 / 40
10 Introduction to Hyperbolic Equations Example on Hyperbolic Equations and Balance Laws Convection Equation and Convection-diffusion Equation The simplest example is the 1D convection equation: u t + a u x = 0, in which the flowing velocity of the fluid is a constant a, the flux is simply f (x, t, u) = f (u) = a u, and the source term is zero. Another simple example is the 1D convection-diffusion equation: u t + a u x = σu xx, in which diffusion as well as convection is considered, and the flux is of the form f (x, t, u, u x ) = a u σ u x, where σ is the diffusion parameter. Note: convection equation is also called advection equation, convection-diffusion equation is also called advection-diffusion equation. 10 / 40
11 Introduction to Hyperbolic Equations Characteristics and Riemann Invariants Characteristics and Riemann Invariants In hyperbolic equations (systems), the information propagates at finite characteristic speeds. Consider a 1D constant-coefficient hyperbolic system u t + A u x = 0. 1 AR = RΛ, Λ = diag(λ 1,, λ p ) with λ 1 λ 2 λ p. 2 Define w = R 1 u, the system is transformed to w t + Λ w x = 0. 3 p families of characteristics given by p characteristic equations dx i dt = λ i, i = 1,, p. 11 / 40
12 Introduction to Hyperbolic Equations Characteristics and Riemann Invariants Riemann Invariants and General Solutions of Initial Value Problems 4 w i is a constant on each characteristics of the ith family: dw i (x i (t), t) dt = [w i,t + λ i w i,x ] (x i (t), t) = 0. 5 w i (x, t) = w i (x λ i t, 0), i = 1,, p. 6 w i, i = 1,, p, are called the Riemann invariants of the system with respect to the characteristics of (λ i, ξ i ), with Aξ i = λ i ξ i, i = 1,, p. 12 / 40
13 Introduction to Hyperbolic Equations Characteristics and Riemann Invariants Riemann Invariants and General Solutions of Initial Value Problems 7 Denote ξ i the ith column of R, and η i the ith row of R 1. 8 By definition η i u = w i, in particular w i (x, 0) = η i u 0 (x). 9 By definition u(x, t) = Rw(x, t), thus, by 5 and 8, u(x, t) = p ξ i w i (x, t) = i=1 p ξ i η i u 0 (x λ i t). i=1 13 / 40
14 Introduction to Hyperbolic Equations Characteristics and Riemann Invariants Domains of Dependance, Influence and Determination Since u(x, t) = p i=1 ξi η i u 0 (x λ i t), we define 1 domain of dependance of a set Ω T R R + : D(Ω T ) = {y R : y = x λ i t, i = 1,, p, (x, t) Ω T }. 2 domain of influence of a set Ω R: I (Ω) = {(x, t) R R + : 1 i p, s.t. x λ i t Ω}, 3 domain of determination of a set Ω R: K(Ω) = {(x, t) R R + : x i = x λ i t Ω, i = 1,, p}. 4 We can also define D(Ω T, t 0 ), I (Ω, t 0 ) and K(Ω, t 0 ) by replacing λ i t by λ i (t t 0 ) in the above definitions. 14 / 40
15 Introduction to Hyperbolic Equations Characteristics and Riemann Invariants Boundary Conditions in Initial-Boundary Value Problems For an initial-boundary value problem defined on (0, 1) R +, suppose λ 1 <... < λ l < 0 < λ r < λ p (l = r 1 or l = r 2). By the characteristics of the system, boundary conditions should be imposed in the following way: 1 p r + 1 linearly independent boundary conditions on the left boundary 0, since {w i } p i=r picks up information from left. 2 l linearly independent boundary conditions on the right boundary 1, since {w i } l i=1 picks up information from right. 15 / 40
16 Introduction to Hyperbolic Equations Characteristics and Riemann Invariants Boundary Conditions in Initial-Boundary Value Problems 3 The p r + 1 boundary conditions on the left boundary 0 must contain sufficient information for {w i } p i=r. 4 The l boundary conditions on the right boundary 1 must contain sufficient information for {w i } l i=1. 5 The domains of dependance, influence and determination can also be defined for initial-boundary value problems. 16 / 40
17 Introduction to Hyperbolic Equations Characteristics and Riemann Invariants The Characteristic Method for the Convection Equation 1 Convection equation: u t (x, t) + a(x, t) u x (x, t) = 0, x I = (x l, x r ) R, t > 0; 2 Characteristics: x (t) = a(x, t), x I, t > 0; 3 Solution is a constant on a characteristic curve: du(x(t), t) dt = u t + u dx x dt = u t + a u x = 0. 4 Suppose a(x l, t) > 0 and a(x r, t) > 0 for all t / 40
18 Introduction to Hyperbolic Equations Characteristics and Riemann Invariants The Characteristic Method for the Convection Equation 5 Initial-boundary conditions u(x, 0) = u 0 (x), u(x l, t) = u l (t). 6 Take x l x 1 < < x N < x r, and 0 t 1 < < t M. 7 Solving the characteristic equation with initial conditions: x(0) = x i, i = 1,, N, and x(t m ) = x l, m = 1,, M, to obtain characteristic curves; 8 On these (approximate) characteristic curves, set U i (x, t) := u 0 (x i ), i = 1,, N, and U m (x, t) := u(x l, t m ) = u l (t m ), m = 1,, M respectively. 18 / 40
19 Introduction to Hyperbolic Equations Characteristics and Riemann Invariants The Characteristic Method for the Convection Equation In particular, if a(x, t) a is a constant, then, we have 9 for a > 0, u 0 ( (x at), if x at x l, u(x, t) = u l t x x ) l, if x at < x l. a 10 if a < 0, u 0 (x ( at), if x at x r, u(x, t) = u r t x x ) r, if x at > x r. a 19 / 40
20 Introduction to Hyperbolic Equations Characteristics and Riemann Invariants Inhomogeneous Convection Equation and Characteristic Method 1 Convection equation: u t (x, t) + a(x, t) u x (x, t) = ψ(x, t), x I = (x l, x r ) R, t > 0; 2 Characteristics: x (t) = a(x, t), x I, t > 0; 3 On a characteristic curve, the solution satisifes: du(x(t), t) dt = u t + u dx x dt = u t + a u x = ψ(x(t), t). 4 Suppose a(x l, t) > 0 and a(x r, t) > 0 for all t / 40
21 Introduction to Hyperbolic Equations Characteristics and Riemann Invariants Inhomogeneous Convection Equation and Characteristic Method 5 Initial-boundary conditions u(x, 0) = u 0 (x), u(x l, t) = u l (t). 6 Take x l x 1 < < x N < x r, and 0 t 1 < < t M. 7 Solving the characteristic equation with initial conditions: x i (0) = x i, i = 1,, N, and x m (t m ) = x l, m = 1,, M, to obtain characteristic curves; 8 Along these (approximate) characteristic curves, solving ODEs systems U i (t) = ψ(x i (t), t) and U m (t) = ψ(x m (t), t) with initial conditions U i (0) := u 0 (x i ), i = 1,, N, and U m (t m ) := u(x l, t m ) = u l (t m ), m = 1,, M respectively. 21 / 40
22 FDM for 1D 1st Order Linear Hyperbolic Equation Characteristic Lines and CFL Condition The Domain of Dependance of a Difference Scheme 1 Convection equation: u t (x, t) + a u x (x, t) = 0, (a = constant). 2 The simplest scheme: U m+1 U m + a Um U 1 m = 0; τ h 3 or equivalently: U m+1 = (1 ν)u m + νu 1 m, (ν = aτ/h). 4 For P = (x, t m+1 ), D(P) = Q = x a t m+1, the domain of dependance of the scheme D h (P) = {x m 1,, x 1, x }. t m+1 P m Q > m 1 Q Q x 22 / 40
23 FDM for 1D 1st Order Linear Hyperbolic Equation Characteristic Lines and CFL Condition The CFL Condition of a Difference Scheme Simple observation: if Q = Q > < x m 1, or Q = Q > x, then U m+1 can not properly approximate u(x, t m+1 ). CFL condition: a necessary condition for a numerical scheme to converge is that the domain of dependence of the the real solution is contained, at least in the sense of limit, in the domains of dependence of the numerical scheme. t m+1 P m Q > m 1 Q Q x 23 / 40
24 FDM for 1D 1st Order Linear Hyperbolic Equation Characteristic Lines and CFL Condition The CFL Condition is a Necessary Condition for Stability The CFL condition is often used as a necessary condition for the stability of finite difference schemes for hyperbolic differential equations. In general, the CFL condition is not a sufficient condition for the stability of a numerical method. 24 / 40
25 FDM for 1D 1st Order Linear Hyperbolic Equation Characteristic Lines and CFL Condition CFL Condition is not a Sufficient Condition for Stability Consider the finite difference scheme of the convection equation U m+1 U m + a Um +1 Um 1 = 0, a R 1 \ {0}. τ 2h The CFL condition for the scheme is ν = a τ/h 1. The Fourier modes solution U m = λ m k with amplification factor λ k = 1 i ν sin kπ N. Since λ k > 1 for all ν 0, the scheme is always unstable, whether the CFL condition is satisfied or not. kπ ei N Note: For initial value problems or initial-boundary value problems with periodic boundary conditions of a constant coefficient convection equation, a necessary and sufficient condition for L 2 -stability is the so called von Neumann condition: λ k 1 + Kτ, for all N + 1 k N. 25 / 40
26 FDM for 1D 1st Order Linear Hyperbolic Equation The Upwind Scheme The Upwind Scheme for the Convection Equation For the convection equation u t (x, t) + a(x, t) u x (x, t) = 0, by the CFL condition, the simplest difference scheme is { U m+1 U m ν m + U m, if a m 0, = U m ν m U m, if a m 0, where ν m = a m τ/h satisfies ν m 1, or equivalently { U m+1 (1 + ν m )U m ν m U+1 m =, if am 0, (1 ν m )U m + ν m U 1 m, if am 0, which may be viewed as obtained by using +t / t to approximate / t; +x / x (a < 0) or x / x (a > 0) to approximate / x. 26 / 40
27 FDM for 1D 1st Order Linear Hyperbolic Equation The Upwind Scheme The Upwind Scheme for the Convection Equation Characteristic method + Linear interpolation Upwind scheme. Let a > 0 be a constant, assume ν = aτ/h 1. 1 By the characteristic method: u(x, t m+1 ) = u(x a τ, t m ); 2 By linear interpolation: u(x a τ, t m ) x (x aτ) h 3 Hence u(x, t m+1 ) νu m 1 + (1 ν)um ; u 1 m + (x aτ) x 1 u m ; h 4 This leads to the scheme: U m+1 = (1 ν)u m + νu m 1, which is exactly the upwind scheme. 27 / 40
28 FDM for 1D 1st Order Linear Hyperbolic Equation The Upwind Scheme The Truncation Error of the Upwind Scheme O(τ + h) By Taylor series expansion, we have T m := um+1 u m + a um u m [ 1 u τ h t + a u ] m x = 1 2 [τu tt a h u xx ] m + 1 [ τ 2 u ttt + a h 2 ] m u xxx + 6 = [ 1 2 a h (1 ν)u xx a h2 (1 ν 2 )u xxx + ] m, if a > 0. Note, here a is a constant, so u tt = a 2 u xx and ν = aτ/h, and τu tt = ahνu xx. Similarly τ k t k+1 u = ( 1) k+1 ah k ν k x k+1 u. 28 / 40
29 FDM for 1D 1st Order Linear Hyperbolic Equation The Upwind Scheme The Truncation Error of the Upwind Scheme O(τ + h) Similarly, we have T m := um+1 u m + a um +1 [ um u τ h t + a u ] m x = 1 2 [τu tt + a h u xx ] m + 1 [ τ 2 u ttt + a h 2 ] m u xxx + 6 = [ 1 2 a h (1 + ν)u xx a h2 (1 ν 2 )u xxx + ] m, if a < 0. Note, here a is a constant, so u tt = a 2 u xx and ν = aτ/h, and τu tt = ahνu xx. Similarly τ k t k+1 u = ( 1) k+1 ah k ν k x k+1 u. 29 / 40
30 FDM for 1D 1st Order Linear Hyperbolic Equation The Upwind Scheme The Stability and Convergence of the Upwind Scheme The error e m = U m u m satisfies the following error equation e m+1 = { (1 ν )e m + ν e+1 m τt m, if a < 0, (1 ν )e m + ν e 1 m τt m, if a > 0. If the CFL condition is satisfied, i.e. ν 1, then, we have E m+1 := max e m+1 E m + τ max T m E 0 + t max max T m, m, (m + 1)τ t max. 30 / 40
31 FDM for 1D 1st Order Linear Hyperbolic Equation The Upwind Scheme The Solutions to Nonlinear Hyperbolic Equations are Generally not Smooth Therefore, if the CFL condition is satisfied, the upwind scheme is stable in L norm, and its convergence rate is O(τ + h), if the solution u is sufficiently smooth. However, in hyperbolic problems, discontinuities are commonly seen in the physical solutions. In fact, the discontinuity creation and propagation is a focus of investigation on nonlinear problems. 31 / 40
32 Nonlinear Convection Equation and Blow Up of Its Classical Solution 1 Nonlinear convection equation: u t (x, t) + a(u(x, t)) u x (x, t) = 0; 2 Characteristic equation: x (t) = a(u(x, t)), t > 0; 3 Solution is a constant on a characteristic curve: du(x(t), t) dt = u t + u dx x dt = u t + a u x = 0. 4 On a characteristic curve, a(u(x, t)) is a constant. 5 Characteristic curves are straight lines (with different slopes). 6 Suppose for ξ 1 < ξ 2, a 1 := a(u 0 (ξ 1 )) > a 2 := a(u 0 (ξ 2 )), then, at time t = (ξ 2 ξ 1 )/(a 1 a 2 ) > 0, the two corresponding characteristic lines intersect: (ξ 1 + a 1 t, t) = (ξ 2 + a 2 t, t). 7 In such a case, the classical solution blows up.
33 Weak Solutions for Nonlinear Conservation Laws Since discontinuities are inevitable, we should take it into consideration by introducing the concept of weak solutions. Definition u L 1 loc (R R +) is called a weak solution to the following 1D initial value problem of the conservation law u t (x, t) + f (u(x, t)) x = 0, u(x, 0) = u 0 (x), if u satisfies the initial condition, and for all fixed 0 t b < t a and < x l < x r <, the following equation holds xr x l u(x, t a )dx = xr x l u(x, t b )dx + ta t b f (u(x l, t))dt ta f (u(x r, t))dt. t b
34 Shock Speed and Rankine-Hugoniot Jump Condition 1 Suppose there is an isolated discontinuity (i.e. a shock) x(t) propagating at speed s(t) (i.e. x (t) = s(t)); 2 By the integral form conservation law, for x 1 = x(t 1 ) and x 1 + x = x(t 1 + t), we have x1 + x x 1 = u(x, t 1 + t)dx t1 + t t 1 f (u(x 1, t))dt x1 + x x 1 t1 + t t 1 u(x, t 1 )dx f (u(x 1 + x, t))dt. 3 Denote u l = u(x 1 0, t 1 ), u r = u(x 1 + 0, t 1 ), then, we have (u l u r ) x = (f (u l ) f (u r )) t + O( t 2 ). 4 Since x (t) = s(t), we are lead to the Rankine-Hugoniot ump condition: s[u] = [f ].
35 FDM for 1D 1st Order Linear Hyperbolic Equation Nonlinear Hyperbolic Equations and Weak Solutions Weak Solutions of Non-viscous Burger s Equation Consider the initial value problem of the Burgers equation { u t + 1 u 2 1, if x < 0; = 0; u(x, 0) = 2 x 0, if x 0. By the Rankine-Hugoniot ump condition, s = f (u r ) f (u l ) u r u l = 1 2 u2 r 1 2 u2 l u r u l = = 1 2. Hence the weak solution of the problem is given by 1, if x < t/2; u(x, t) = 0, if x t/2. 35 / 40
36 FDM for 1D 1st Order Linear Hyperbolic Equation Nonlinear Hyperbolic Equations and Weak Solutions Weak Solutions of Non-viscous Burger s Equation Consider the initial value problem of the equation u 2 t + 2 u 3 3 x = 0; u(x, 0) = { 1, if x < 0; 0, if x 0. Define w = u 2, by the Rankine-Hugoniot ump condition, s = f (w r ) f (w l ) w r w l = 2 3 w r 3/2 2 3 w 3/2 l = w r w l 0 1 = 2 3. Hence the weak solution of the problem is given by 1, if x < 2t/3; u(x, t) = 0, if x 2t/3. 36 / 40
37 FDM for 1D 1st Order Linear Hyperbolic Equation Nonlinear Hyperbolic Equations and Weak Solutions Weak Solutions of Non-viscous Burger s Equation For smooth u, both the Burger s equation u t + 1 u 2 2 x = 0 and the equation u 2 t + 2 u 3 3 x = 0 are equivalent to the first order nonlinear hyperbolic equation u t + uu x = / 40
38 FDM for 1D 1st Order Linear Hyperbolic Equation Nonlinear Hyperbolic Equations and Weak Solutions Key Ingredients of a Conservation Law The key ingredients of a conservation law are the conservative quantities and the corresponding fluxes; the integral form equations based on physical conservation laws; Entropy condition to distinguish the physical solution from the others. Remark: For a nonlinear problem, the weak solution is generally not unique even for a conservation law (see Exercise 3.4). 38 / 40
39 FDM for 1D 1st Order Linear Hyperbolic Equation Nonlinear Hyperbolic Equations and Weak Solutions L 2 Is a Better Space to Investigate Hyperbolic Equations The upwind scheme of the convection equation satisfies the maximum principle, if CFL condition is satisfied; However, for general hyperbolic equations (systems), the solution do not satisfy the maximum principle. Hyperbolic equations (systems) are often used to characterize the propagation and evolution of waves. Fourier modes are waves of various frequencies propagating at their own characteristic speeds. It is natural to apply the Fourier method to analyze the L 2 stability and accuracy of finite difference schemes for hyperbolic equations (systems). 39 / 40
40 SK 3µ1, 3 Thank You!
Numerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem in a 2D Box Region Let us consider a model problem of parabolic
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Modified Equation of a Difference Scheme What is a Modified Equation of a Difference
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 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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More information0.3.4 Burgers Equation and Nonlinear Wave
16 CONTENTS Solution to step (discontinuity) initial condition u(x, 0) = ul if X < 0 u r if X > 0, (80) u(x, t) = u L + (u L u R ) ( 1 1 π X 4νt e Y 2 dy ) (81) 0.3.4 Burgers Equation and Nonlinear Wave
More informationMath Partial Differential Equations 1
Math 9 - Partial Differential Equations Homework 5 and Answers. The one-dimensional shallow water equations are h t + (hv) x, v t + ( v + h) x, or equivalently for classical solutions, h t + (hv) x, (hv)
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationAMath 574 February 11, 2011
AMath 574 February 11, 2011 Today: Entropy conditions and functions Lax-Wendroff theorem Wednesday February 23: Nonlinear systems Reading: Chapter 13 R.J. LeVeque, University of Washington AMath 574, February
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationModeling Rarefaction and Shock waves
April 30, 2013 Inroduction Inroduction Rarefaction and shock waves are combinations of two wave fronts created from the initial disturbance of the medium. Inroduction Rarefaction and shock waves are combinations
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems Professor Dr E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8, 2014 1 / 56 We study some basic
More informationNotes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
Outline This lecture Diffusion and advection-diffusion Riemann problem for advection Diagonalization of hyperbolic system, reduction to advection equations Characteristics and Riemann problem for acoustics
More informationNon-linear Scalar Equations
Non-linear Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro August 24, 2014 1 / 44 Overview Here
More information1 Upwind scheme for advection equation with variable. 2 Modified equations: numerical dissipation and dispersion
1 Upwind sceme for advection equation wit variable coefficient Consider te equation u t + a(x)u x Applying te upwind sceme, we ave u n 1 = a (un u n 1), a 0 u n 1 = a (un +1 u n ) a < 0. CFL condition
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, x < 1, 1 + x, 1 < x < 0, ψ(x) = 1 x, 0 < x < 1, 0, x > 1, so that it verifies ψ 0, ψ(x) = 0 if x 1 and ψ(x)dx = 1. Consider (ψ j )
More informationNumerical Methods for Hyperbolic Conservation Laws Lecture 4
Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Wen Shen Department of Mathematics, Penn State University Email: wxs7@psu.edu Oxford, Spring, 018 Lecture Notes online: http://personal.psu.edu/wxs7/notesnumcons/
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2017 Randall J. LeVeque Applied Mathematics University of Washington January 4, 2017 http://faculty.washington.edu/rjl/classes/am574w2017
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationQuasi-linear first order equations. Consider the nonlinear transport equation
Quasi-linear first order equations Consider the nonlinear transport equation u t + c(u)u x = 0, u(x, 0) = f (x) < x < Quasi-linear first order equations Consider the nonlinear transport equation u t +
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationHyperbolic Systems of Conservation Laws. I - Basic Concepts
Hyperbolic Systems of Conservation Laws I - Basic Concepts Alberto Bressan Mathematics Department, Penn State University Alberto Bressan (Penn State) Hyperbolic Systems of Conservation Laws 1 / 27 The
More informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationLecture Notes on Hyperbolic Conservation Laws
Lecture Notes on Hyperbolic Conservation Laws Alberto Bressan Department of Mathematics, Penn State University, University Park, Pa. 16802, USA. bressan@math.psu.edu May 21, 2009 Abstract These notes provide
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationScalar Conservation Laws and First Order Equations Introduction. Consider equations of the form. (1) u t + q(u) x =0, x R, t > 0.
Scalar Conservation Laws and First Order Equations Introduction. Consider equations of the form (1) u t + q(u) x =, x R, t >. In general, u = u(x, t) represents the density or the concentration of a physical
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 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 informationLecture 5.7 Compressible Euler Equations
Lecture 5.7 Compressible Euler Equations Nomenclature Density u, v, w Velocity components p E t H u, v, w e S=c v ln p - c M Pressure Total energy/unit volume Total enthalpy Conserved variables Internal
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2011 Randall J. LeVeque Applied Mathematics University of Washington January 3, 2011 R.J. LeVeque, University of Washington AMath 574,
More informationHyperbolic Systems of Conservation Laws
Hyperbolic Systems of Conservation Laws III - Uniqueness and continuous dependence and viscous approximations Alberto Bressan Mathematics Department, Penn State University http://www.math.psu.edu/bressan/
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationMath 126 Final Exam Solutions
Math 126 Final Exam Solutions 1. (a) Give an example of a linear homogeneous PE, a linear inhomogeneous PE, and a nonlinear PE. [3 points] Solution. Poisson s equation u = f is linear homogeneous when
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan
Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global
More informationNumerische Mathematik
Numer. Math. (2007 106:369 425 DOI 10.1007/s00211-007-0069-y Numerische Mathematik Hyperbolic balance laws: Riemann invariants and the generalized Riemann problem Matania Ben-Artzi Jiequan Li Received:
More informationThe Hopf equation. The Hopf equation A toy model of fluid mechanics
The Hopf equation A toy model of fluid mechanics 1. Main physical features Mathematical description of a continuous medium At the microscopic level, a fluid is a collection of interacting particles (Van
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationConvergence and Error Bound Analysis for the Space-Time CESE Method
Convergence and Error Bound Analysis for the Space-Time CESE Method Daoqi Yang, 1 Shengtao Yu, Jennifer Zhao 3 1 Department of Mathematics Wayne State University Detroit, MI 480 Department of Mechanics
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 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 informationPropagation of discontinuities in solutions of First Order Partial Differential Equations
Propagation of discontinuities in solutions of First Order Partial Differential Equations Phoolan Prasad Department of Mathematics Indian Institute of Science, Bangalore 560 012 E-mail: prasad@math.iisc.ernet.in
More informationThe Discontinuous Galerkin Method for Hyperbolic Problems
Chapter 2 The Discontinuous Galerkin Method for Hyperbolic Problems In this chapter we shall specify the types of problems we consider, introduce most of our notation, and recall some theory on the DG
More informationCLASSROOM NOTES PART II: SPECIAL TOPICS. APM526, Spring 2018 Last update: Apr 11
CLASSROOM NOTES PART II: SPECIAL TOPICS APM526, Spring 2018 Last update: Apr 11 1 Function Space Methods General Setting: Projection into finite dimensional subspaces t u = F (u), u(t = 0) = u I, F : B
More informationChp 4: Non-linear Conservation Laws; the Scalar Case. By Prof. Dinshaw S. Balsara
Chp 4: Non-linear Conservation Laws; the Scalar Case By Prof. Dinshaw S. Balsara 1 4.1) Introduction We have seen that monotonicity preserving reconstruction and iemann solvers are essential building blocks
More informationConservation and dissipation principles for PDEs
Conservation and dissipation principles for PDEs Modeling through conservation laws The notion of conservation - of number, energy, mass, momentum - is a fundamental principle that can be used to derive
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 informationAPM526 CLASSROOM NOTES PART IB. APM526, Spring 2018 Last update: MAR 12
APM526 CLASSROOM NOTES PART IB APM526, Spring 2018 Last update: MAR 12 1 PARTICLE METHODS 1. Particle based methods: Large systems of ODEs for individual particles. 2. Partial differential equations: Equations
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 informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
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 informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Discretization of Boundary Conditions Discretization of Boundary Conditions On
More informationHyperbolic Partial Differential Equations
Chapter 1 Hyperbolic Partial Differential Equations We begin our study of finite difference methods for partial differential equations by considering the important class of partial differential equations
More informationA stochastic particle system for the Burgers equation.
A stochastic particle system for the Burgers equation. Alexei Novikov Department of Mathematics Penn State University with Gautam Iyer (Carnegie Mellon) supported by NSF Burgers equation t u t + u x u
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationResearch Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation with Inconsistent Initial and Boundary Conditions
International Mathematics and Mathematical Sciences Volume 212, Article ID 82197, 13 pages doi:1.1155/212/82197 Research Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation
More informationCoupling conditions for transport problems on networks governed by conservation laws
Coupling conditions for transport problems on networks governed by conservation laws Michael Herty IPAM, LA, April 2009 (RWTH 2009) Transport Eq s on Networks 1 / 41 Outline of the Talk Scope: Boundary
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationWaves and characteristics: Overview 5-1
Waves and characteristics: Overview 5-1 Chapter 5: Waves and characteristics Overview Physics and accounting: use example of sound waves to illustrate method of linearization and counting of variables
More informationHIGH ORDER FINITE VOLUME SCHEMES BASED ON RECONSTRUCTION OF STATES FOR SOLVING HYPERBOLIC SYSTEMS WITH NONCONSERVATIVE PRODUCTS
MATHEMATICS OF COMPUTATION Volume 75, Number 255, July 26, Pages 113 1134 S 25-5718(6)1851-5 Article electronically published on March 21, 26 HIGH ORDER FINITE VOLUME SCHEMES BASED ON RECONSTRUCTION OF
More informationLast time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler:
Lecture 7 18.086 Last time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler: U j,n+1 t U j,n = U j+1,n 2U j,n + U j 1,n x 2 Expected accuracy: O(Δt) in time,
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 informationStability of Patterns
Stability of Patterns Steve Schecter Mathematics Department North Carolina State University 1 2 Overview Motivation Systems of viscous conservation laws and their variants are partial differential equations
More informationPart 1. The diffusion equation
Differential Equations FMNN10 Graded Project #3 c G Söderlind 2016 2017 Published 2017-11-27. Instruction in computer lab 2017-11-30/2017-12-06/07. Project due date: Monday 2017-12-11 at 12:00:00. Goals.
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
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 informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, < 1, 1 +, 1 < < 0, ψ() = 1, 0 < < 1, 0, > 1, so that it verifies ψ 0, ψ() = 0 if 1 and ψ()d = 1. Consider (ψ j ) j 1 constructed as
More informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
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 informationThe Riemann problem. The Riemann problem Rarefaction waves and shock waves
The Riemann problem Rarefaction waves and shock waves 1. An illuminating example A Heaviside function as initial datum Solving the Riemann problem for the Hopf equation consists in describing the solutions
More informationSistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas
Sistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas Pedro González-Casanova Henríquez Unidad de Investigación en Cómputo Aplicado DGSCA, UNAM Elvira Torondel y Ricard Garrido
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 informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More information12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis
ATH 337, by T. Lakoba, University of Vermont 113 12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis 12.1 Formulation of the IBVP and the minimax property of its
More informationNumerical Integration of Linear and Nonlinear Wave Equations
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Dissertations, Theses, and Student Research Papers in Mathematics Mathematics, Department of 12-2004 Numerical Integration
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More informationSTRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN PROBLEM FOR A NON-STRICTLY HYPERBOLIC SYSTEM WITH FLUX APPROXIMATION
Electronic Journal of Differential Equations, Vol. 216 (216, No. 126, pp. 1 16. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu STRUCTURAL STABILITY OF SOLUTIONS TO THE RIEMANN
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 informationLINEAR SECOND-ORDER EQUATIONS
LINEAR SECON-ORER EQUATIONS Classification In two independent variables x and y, the general form is Au xx + 2Bu xy + Cu yy + u x + Eu y + Fu + G = 0. The coefficients are continuous functions of (x, y)
More informationOn the Cauchy Problems for Polymer Flooding with Gravitation
On the Cauchy Problems for Polymer Flooding with Gravitation Wen Shen Mathematics Department, Penn State University Abstract We study two systems of conservation laws for polymer flooding in secondary
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 informationHomework For each of the following matrices, find the minimal polynomial and determine whether the matrix is diagonalizable.
Math 5327 Fall 2018 Homework 7 1. For each of the following matrices, find the minimal polynomial and determine whether the matrix is diagonalizable. 3 1 0 (a) A = 1 2 0 1 1 0 x 3 1 0 Solution: 1 x 2 0
More information