Convergence and Error Bound Analysis for the Space-Time CESE Method
|
|
- Roland Doyle
- 6 years ago
- Views:
Transcription
1 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 Wayne State University Detroit, MI Department of Mathematics and Statistics University of Michigan Dearborn Dearborn, MI 4818 Received 17 September 1999; accepted 31 August 000 In this work, we study the convergence behavior of a recently developed space-time conservation element and solution element method for solving conservation laws. In particular, we apply the method to a onedimensional time-dependent convection-diffusion equation possibly with high Peclet number. We prove that the scheme converges and we obtain an error bound. This method performs well even for strong convection dominance over diffusion with good long-time accuracy. Numerical simulations are performed to verify the results. c 001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 64 78, 001 Keywords: Space-time CESE method; a µ scheme; convection diffusion equation I. INTRODUCTION In this article, we study the convergent behavior of a recently developed numerical method: the space-time conservation element and solution element method the space-time CESE. This new method is proposed in [1] for conservation laws and has substantial differences in both concept and methodology from well-known methods for solving conservation laws, such as finite difference, finite volume, and finite element methods. Mathematically, conservation laws are represented by a set of integral equations. The differential form of these laws is obtained from the integral form with the assumption that the physical Correspondence to: Jennifer Zhao, Dept. of Mathematics and Statistics, Univ. of Michigan-Dearborn, 4901 Evergreen Rd, Dearborn, MI xich@umich.edu, xich c 001 John Wiley & Sons, Inc.
2 CONSERVATION ELEMENT/SOLUTION ELEMENT METHOD 65 solution is smooth. But, for a physical solution in a region of rapid change, this smoothness assumption may be difficult to realize. Thus, a method designed to obtain numerical solutions to the differential form without enforcing flux conservation may be at a fundamental disadvantage in modeling many physical phenomena, especially when obtaining accurate solutions for local and long-time behaviors. By contrast, the space-time CESE method was developed to enforce flux-conservation locally and globally in space and time. It always retains the basic physical reality of flux conservation, even in a region involving discontinuities. In summary, the spacetime CESE method has the following major advantages. First, it enforces both local and global flux conservation in space and time. Second, space and time are unified and are treated on the same footing. Third, dependent variables and their derivatives are solved simultaneously with comparable accuracy. The goal of this article is to conduct convergence and error bound analysis for this method. To our best knowledge, this work is the first approach in this direction. In order to provide a clear description of the method and conduct our study without involving too many technical difficulties, we carry out our study by applying the space-time CESE method to the following onedimensional time-dependent convection-diffusion equation with possible convection dominance over diffusion: u t + a u x u µ x As shown in the next section, the a µ scheme thus obtained is explicit and has two independent marching variables. It has the unusual property that its stability condition is independent of the diffusion coefficient µ, and there are no more mesh requirements other than the stability condition for convergence. We prove that this scheme is convergent, and that the global error for the solution in L norm of space is of order x. Most related works can be found in [1] and []. In [1], the a µ scheme is derived and its stability is discussed. In both [1] and [], the space-time CESE method is extended to other equations, such as the 1D time-dependent Navier Stokes equations of a perfect gas, and the 1D convection-diffusion equation of nonconstant velocity and viscosity coefficients. There is not much mathematical study to the method in [1] or []; that is what this article is for. The order of the article is organized as following. In Section II, we describe the a µ scheme. In Section III, we prove its convergence and derive an error bound. Numerical simulations are performed in Section IV. II. A µ SCHEME In this section, we apply the space-time CESE method proposed in [1] to the following Cauchy problem of one-dimensional convection-diffusion equation: u t + a u x µ u x 0, 0 <t<t,ux, 0 u 0x,.1 a is the convection velocity, µ>0 is the diffusion coefficient, both are assumed to be constants, and T is a positive constant representing the simulation ending time. Let x 1 x, x t be considered as the coordinates in the half-plane of a two-dimensional Euclidean space R + {x, t : <x<, 0 t< }. Using the Gauss divergence theorem in R +, Eq..1 can be written as the differential form of the following conservation
3 66 YANG, YU, AND ZHAO law: SV h ds 0.. Here SV is the boundary of an arbitrary space-time region V in R+, h au µ u x,u is a current density vector in R+, ds dσ n with dσ and n being, respectively, the area and the outward unit normal of a surface element on SV. We use x to represent the space grid size, t the time grid size, and x j,t n a particular grid point, x j j x, <j<+,t n n t, n 0,.3 see Fig. 1. We denote by SEj, n the interior of a space-time region bounded by a dashed curve and call it a solution element at grid point x j,t n ; see Fig. A. For our discussion, the exact size of this region does not matter. For any x, t SEj, n, weuseux, t; j, n and Hx, t; j, n to approximate ux, t and hx, t: Ux, t; j, n U n j +U x n j x x j +U t n j t t n, x, t SEj, n,.4 U n j, U x n j, and U t n j are independent constants in SEj, n, which approximate ux j,t n, u x x j,t n, and u t x j,t n, respectively. We require Ux, t; j, n to satisfy.1 within SEj, n, this implies that Therefore, We also define U t n j au x n j..5 Ux, t; j, n U n j +U x n j x x j at t n, x, t SEj, n..6 Hx, t; j, n aux, t; j, n µu x x, t; j, n,ux, t; j, n..7 FIG. 1. Staggered grids at full and half steps.
4 CONSERVATION ELEMENT/SOLUTION ELEMENT METHOD 67 FIG.. Solution elements and conservation elements. Let R+ be divided into nonoverlapping rectangular regions referred to as conservation elements CE. The CE with its top-right vertex being the mesh point j, n Ω is denoted by CE j, n; see Fig. B, and the CE whose top-left vertex is the mesh point j, n is denoted by CE + j, n; see Fig. C. Similarly CE + j 1/,n+1/ is illustrated in Fig. D and CE j +1/,n+ 1/ is illustrated in Fig. E. We approximate. by F + j, n H ds 0,F j, n H ds 0..8 SCE +j,n SCE j,n Because each SCE ± j, n is a simple closed curve in R+, the surface integral in.8 can be converted to a line integral. Let G U, au µu x, dr dx, dt..9 Then Therefore, F + j, n SCE +j,n H ds ± G dr..10 G dr, F j, n SCE j,n G dr,.11 with the understanding that the integral is evaluated in the counterclockwise direction. By evaluating.11, we have that 4 x F ±j, n ± 1 [1 τ + δu x n j +1 τ δu x n 1/ j±1/ ] If we define + 1 ± τ x Uj n U n 1/ j±1/,.1 τ a t 4µ t,δ x x..13 qj, n U n j x 4 U x n j,.14 then solving the two equations: 4 x F 4 +j, n 0, x F j, n 0,.15
5 68 YANG, YU, AND ZHAO leads to the a µ scheme: qj, n Q + qj 1,n 1+Q + Q + Q Q + qj, n 1 + Q qj +1,n 1,.16 Q τ 1 τ δ 1 τ 1 τ +δ 1 τ1 τ δ 1 τ +δ,q 1 1 τ 1 τ δ 1 τ 1 τ +δ 1+τ1 τ δ 1 τ +δ..17 This scheme is explicit and has two independent marching variables Uj n the following initial data for.16: x and 4 U x n j. We take Uj 0 u 0 x j, U x 0 j u 0x j..18 x III. CONVERGENCE AND ERROR BOUND ANALYSIS In this section, we show that the a µ scheme.16 converges and we derive its exact rate of convergence error bound. We first give the following stability result. Lemma 3.1. The a µ scheme.16 is stable if and only if 1 τ + δ 0,τ < Proof. The proof follows from a standard Von Neumann analysis. For details, please see [1]. To provide the convergence and error bound results, we start by introducing some notations. If W, W 1, W 0, W 1, W j +, 3. each W j is a vector, we define W W j x, 3.3 W j represents the standard L -norm for a vector. We define the local truncation error of.16 as j+ ɛ n ɛj, n j+ ɛ1 j, n,n 0, 3.4 ɛ j, n [ ɛ 1 j, n ux j,t n Q uxj 1,t n + x 4 u xx j 1,t n + Q uxj+1,t n 1 x 4 u xx j+1,t n ] uxj,t +Q + Q + Q Q + n x 4 u xx j,t n, ɛ j, n x [ 4 u xx j,t n Q uxj 1,t n + x 4 u xx j 1,t n + Q uxj+1,t n x 4 u xx j+1,t n
6 CONSERVATION ELEMENT/SOLUTION ELEMENT METHOD 69 ] uxj,t +Q + Q + Q Q + x 4 u xx j,t n, 3.6 [ ] [ ] ξ ξ ξ, ζ, 3.7 ζ ζ 1 and ux, t is the exact solution of.1. We now have the following results. Lemma 3.. Let ux, t be the exact solution of.1, and we assume that ux, t is sufficiently regular. Then for any n and j, we have { 1 ɛ 1 j, n t u xxx j,t n a u tt x j,t n t + a [ u xxx x j,t n x u xxxxx j,t n 1 τ δ x a t ] 1 τ t + δ } + C 1 t 3 + x 4, 3.8 ɛ j, n 4δ1 τ { 1 τ 1 τ + δ x + δ 16µ1 τ x u xt x j,t n a1 τ δ 16µ1 τ x u xxx x j,t n + a1 τ x u xx x j,t n + 1 τ + δ 8µ 4δ1 τ t + C 1 1+ τ δ 1 τ 1+ 1 τ δ x 1+τ + C 1 τ δ 1 τ 1+ 1 τ } δ x, τ C 1,C are constants involving higher-order derivatives of ux, t evaluated at some point near x j,t n, but independent of the grid size and coefficients. Proof. 3.8 and 3.9 can be obtained by the standard Taylor expansion under the assumption that ux, t is sufficiently smooth. The proof is quite tedious and technical, and thus is omitted here. Next, we derive an error bound. Define the error vector between the exact solutions and approximate solutions at time level t n as E n, E 1,n, E0,n, E1,n, Ej, n j+, 3.10 Ej, n ux j,t n Uj n x 4 [u xx j,t n U x n j ]. 3.11
7 70 YANG, YU, AND ZHAO By combining.16, 3.5, and 3.6, we get Ej, n Q + Ej 1,n 1 + Q Ej +1,n 1 +Q + Q + Q Q + Ej, n 1 + ɛj, n. 3.1 The following theorem gives the error bound result. Theorem 3.1. Let T be a fixed time, and ux, t be the exact solution of.1, which is assumed to be sufficiently regular. Then for any t n T,the a µ scheme.16 is convergent, if the stability condition 3.1 is satisfied, and 1/ E n Ex j,t n x [ ux j,t n Uj n + ] x 4 u xx j,t n U x n j x C T maxτ, aτ, x, 3.13 C gx, 0 L 1 R 1 + x g x x, 0 L 1 R 1 + x g xt x, t L 1 R 1 0,T 1/ + g t x, t L1 R 1 0,T, 3.14 and gx, t K k1 0 α k K k ux, t C αk x α, 3.15 k t k α k C αk represents bounded positive constant independent of τ,δ, and K is a fixed positive integer. Since ux, t is assumed to be sufficiently regular, C is a bounded, positive constant. Proof. We first define a Fourier transform Ên θ of E n.forθ π, π], we define i 1. Then Ê n θ 1 π + Ên 1 θ Ê n θ π e ijθ Ej, n e ijθ ux j,t n Uj n 1 + x π e ijθ 4 u xx j,t n U x n j, 3.16 Ê n θ 1 e ijθ Q + Ej 1,n 1 π + Q Ej +1,n 1+Q + Q + Q Q + Ej, n 1 + ɛj, n
8 CONSERVATION ELEMENT/SOLUTION ELEMENT METHOD 71 + Qθ Ê n 1 θ+ 1 e ijθ ɛj, n π Qθ Ê n 1 θ+ˆɛ n θ, 3.17 and From 3.17, we have that Clearly, Thus, we have that Qθ e iθ/ Q + + e iθ/ Q cos θ iτ sin θ i1 τ δ sin θ i 1 τ sin θ 1 τ +θ 1 τ +δ 1 τ +δ cos θ + iτ sin θ ˆɛ n θ 1 π + ˆɛ n 1 θ ˆɛ n θ e ijθ ɛj, n 1 π e ijθ ɛ 1 j, n π e ijθ ɛ j, n Ê n θ Qθ Ê n 1 θ+ˆɛ n θ Qθ 4 Ê n θ+qθ ˆɛ n 1 θ+ˆɛ n θ n Qθ n Ê 0 θ+ Qθ n lˆɛ l θ. 3.0 Ê n θ l1 Ê 0 θ n Qθ n lˆɛ l θ. 3. It is proved in [1] that, as long as the stability condition 3.1 is satisfied, then l1 Qθ 1. Qθ represents the matrix- norm of Qθ. By using the fact that E n Ên θ L π,+π see [3], we have that E n n Ên θ L π,π Qθ n lˆɛ l θ l1 n Qθ n l ˆɛ l θ L π,π l1 n ɛ l x 1/ l1 n l1 L π,π Qθ n l ˆɛ l θ L π,π
9 7 YANG, YU, AND ZHAO n l1 [by using the results of 3.8 and 3.9] ɛ 1 j, l + ɛ j, l x a t + a x t + a t + δ µ x3 + aδ µ x3 + δ µ x3 + x t K n l1 k1 0 α k k In Appendix A, we prove that K k1 0 α k k C αk u k C αk x α k t k α k 1/ u k. x α k t k α k x j,t n x x j,t n x 1/, 3.3 gx j,t n x. 3.4 gx j,t n x C. 3.5 So, we get from 3.3 that E n C n a t + a x t + a t + δ µ x3 + x t C T 1 a t + a x t + a t + δ t µ x3 + x t. 3.6 If we use the definition τ a t 4µ t,δ x x, and plug them into 3.6, then 3.13 follows. Remark 3.1. As described in Lemma 3.1, the a µ scheme is not sensitive to µ. Our numerical results confirm this. Remark 3.. Theorem 3. shows several unusual features regarding the convergent behavior of a µ scheme: I Other than the stability condition for x, t, there are no further requirements for x or t regarding the convergence of a µ scheme. II The error bound or convergence rate does not depend on µ. This is confirmed by our numerical experiments. III The error bound depends on T. But unlike standard schemes for time-dependent cases, which involve a factor of e T, the dependence is linear. IV As required by the stability result, τ satisfies 0 <τ<1, so the factor maxτ,aτ in 3.13 is really bounded by maxa, 1, which is a fixed constant. Remark 3.3. Theorem 3. indicates that the global convergence rate for the solution in the L space norm is of order x, but we cannot conclude any global convergence rate for its derivative.
10 CONSERVATION ELEMENT/SOLUTION ELEMENT METHOD 73 TABLE I. Relative errors between exact solution and numerical solution for Example 4.1 with different diffusion coefficients. The temporal grid size t 0.9 x. Simulation ends at T 1.5. µ 1 µ 10 µ e e e-05.6e e-05.58e e e e e e e e-06.03e-06.57e e-05.68e e e e e e e e-06 IV. NUMERICAL EXPERIMENTS In this section, numerical experiments are conducted to support our results in previous sections. In particular, we want to show numerically that the method performs well for small and large Peclet numbers and for short and long simulation time as mentioned in Remark 3.. In the following experiments, we show relative errors between the exact solution and numerical solution in the L 0,T; L [b, d] norm: d 1/ max 0 tn T b ut n,x U n x dx solution error d 1/, 4.1 max 0 tn T b ut n,x dx derivative error max 0 tn T d b max 0 tn T utn,x x d b Ux n x dx utn,x x 1/ 1/, 4. dx the integration in space is evaluated by the composite trapezoidal rule, and U n,u n x are, respectively, the approximate solution and its derivative at t n. Example 4.1. Find u satisfying u t + a u x µ u fx, t, x, t b, d 0,T], 4.3 x ub, t g 0, u x d, t g 1,t [0,T], 4.4 TABLE II. Relative errors between exact solution and numerical solution for Example 4.1 with different diffusion coefficients. The temporal grid size t 0.9 x. Simulation ends at T 1.5. µ 10 6 µ µ e e e e e e e e e e e e e e-06.68e e-06.68e e e e e e e e-06
11 74 YANG, YU, AND ZHAO TABLE III. Relative errors between exact solution and numerical solution for Example 4.1 with different diffusion coefficients. The temporal grid size t 0.9 x. Simulation ends at T 1.5. µ 10 0 µ µ e e e e e e e e e e e e e e-06.68e e-06.68e e e e e e e e-06 TABLE IV. Relative errors between exact solution and numerical solution for Example 4.1 with different diffusion coefficients. The temporal grid size t 0.9 x. Simulation ends at T 6.0. µ 1 µ 10 6 µ e e e e e e e e e-05.70e e-05.68e e-06.50e e e e e e e e e e e-06 ux, 0 u 0 x,x [b, d], 4.5 [b, d] [0, 1],a 1, and the source function fx, t and boundary values g 0 and g 1 are chosen such that the exact solution is ux, t sinx+t. Many values for the ending simulation time T and for the diffusion coefficient µ are given in the numerical results. Tables I III and IV VI show the errors in the solution and its derivative for various simulation ending times, T 1.5,T 6.0,T 1.0, and T 4.0 for Example 4.1. Note that the orders of convergence are consistent with our theorem. We do not observe any indication from the results that the accuracy deteriorates with the Peclet number and time; this is also consistent with our results. When the convergence condition is violated, for example, when a t/ x 1, the convergence can be unpredictable, see Tables VII and VIII. The conclusion is, when the convergence condition is violated, the accuracy can still be good for convection-equal-to-diffusion problems, but can be arbitrarily bad for convection-dominated-diffusion problems. TABLE V. Relative errors between exact solution and numerical solution for Example 4.1 with different diffusion coefficients. The temporal grid size t 0.9 x. Simulation ends at T 1.0. µ 1 µ 10 6 µ e e e e e e e e e-05.70e e-05.68e e-06.50e e e e e e e e e e e-06
12 CONSERVATION ELEMENT/SOLUTION ELEMENT METHOD 75 TABLE VI. Relative errors between exact solution and numerical solution for Example 4.1 with different diffusion coefficients. The temporal grid size t 0.9 x. Simulation ends at T 4.0. µ 1 µ 10 6 µ e e e e e e e e e-05.70e e-05.68e e-06.50e e e e e e e e e e e-06 TABLE VII. Relative errors between exact solution and numerical solution for Example 4.1 with different diffusion coefficients. The temporal grid size t x, which violates the convergence condition and leads to unpredictable results for the derivative. Simulation ends at T 1.5. µ 1 µ 10 4 µ e e e e e e e e e e e e e-06.44e-06.9e e e e e e e e e e+1 Example 4.. Find u satisfying u t + a u x µ u fx, t, x, t b, d 0,T], 4.6 x ub, t g 0, u x d, t g 1,t [0,T], 4.7 ux, 0 u 0 x,x [b, d], 4.8 [b, d] [0, 1],a 1, and the source function fx, t and boundary values g 0 and g 1 are chosen such that the exact solution is ux, t e x+t. Many values for the ending simulation time T and for the diffusion coefficient µ are given in the numerical results. Numerical results for Example 4. are shown in Tables IX XII for simulation ending times at T 1,T 10,T 100, and T 300, and Peclet numbers 1, 10 6, and , respectively. Clearly the accuracy does not show any deterioration with respect to time. It seems that there TABLE VIII. Relative errors between exact solution and numerical solution for Example 4.1 with different diffusion coefficients. The temporal grid size t 1.5 x, which violates the convergence condition and leads to unpredictable results. Simulation ends at T 5.0. µ 1 µ 10 4 µ e e e e e e e e e e e e e e-05 Infinity Infinity Infinity Infinity e e-06 Infinity Infinity Infinity Infinity
13 76 YANG, YU, AND ZHAO TABLE IX. Relative errors between exact solution and numerical solution for Example 4. with different diffusion coefficients. The temporal grid size t 0.9 x. Simulation ends at T 1.0. µ 1 µ 10 6 µ e e e e e e e e e e e e e e-07.69e e-06.69e-06.78e e e e e e e-07 TABLE X. Relative errors between exact solution and numerical solution for Example 4. with different diffusion coefficients. The temporal grid size t 0.9 x. Simulation ends at T µ 1 µ 10 6 µ e e e-05.6e e e e-06 3.e e e e e e e-07.69e e-06.69e e e-07.01e e e e-07.48e-07 is a slight deterioration in accuracy when the diffusion coefficient is decreased from 1 to However, when the diffusion coefficient is further decreased from 10 6 to , the accuracy remains the same. Besides, the solution and its derivative are approximated at the same time with about the same accuracy. As for Example 4.1, the accuracy remains the same when the simulation ending time T is increased. This implies good long-time behavior of the numerical solution obtained by the space-time CESE scheme. TABLE XI. Relative errors between exact solution and numerical solution for Example 4. with different diffusion coefficients. The temporal grid size t 0.9 x. Simulation ends at T µ 1 µ 10 6 µ e e e-05.3e e e e-06 3.e e e e e e e-07.69e e-06.69e e e-07.01e e e e-07.47e-07 TABLE XII. Relative errors between exact solution and numerical solution for Example 4. with different diffusion coefficients. The temporal grid size t 0.9 x. Simulation ends at T µ 1 µ 10 6 µ e e e-05.e e e e-06 3.e e e e e e e-07.69e e-06.69e e e-07.01e e e e-07.47e-07
14 CONSERVATION ELEMENT/SOLUTION ELEMENT METHOD 77 APPENDIX A In this appendix, we prove Eq. 3.5, which is gx j,t n x C, A1 x j is some point in x j 1,x j ],t n is some point in t n 1,t n ]. C is a positive constant defined in Notice that gx, t is a sufficiently regular and positive function. Proof. We start with + gx j,t n x gx, 0 dx I + II, A Then I So t n I II g t x j,t dt x 0 t n x j 0 t n x j I x We now estimate II: 0 gx j,t n x gx j, 0 x + gx j, 0 x; gx, 0 dx A3 x j gx j, 0 gx, 0 dx. A4 g t x j,t g t x, t dt dx + x j x x j g xt w, t dw dx dt + T 0 g xt x, t dx dt + x j t n x j 0 t n x j 0 T 0 g t x, t dxdt g t x, t dx dt. g t x, t dx dt A5 x g xt L1 R 1 [0,T ] + g t L1 R 1 [0,T ]. A6 II x x j x j x g x w, 0 dw dx x j g x w, 0 dw x g x x, 0 L1 R. A7
15 78 YANG, YU, AND ZHAO Equation 3.5 follows from A.6 and A.7. References 1. Sin-Chung Chang, The method of space-time conservation element and solution element A new approach for solving the Navier Stokes and Euler equations, J Comp Ph , S. C. Chang and W. M. To, NASA TM , August 1991 unpublished. 3. B. Noble and J. W. Daniel, Applied linear algebra, nd Ed., Prentice Hall, Englewood Cliffs, NJ, J. W. Thomas, Numerical partial differential equations, Springer Verlag, New York, 1995.
Numerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
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 informationMathematics of Physics and Engineering II: Homework problems
Mathematics of Physics and Engineering II: Homework problems Homework. Problem. Consider four points in R 3 : P (,, ), Q(,, 2), R(,, ), S( + a,, 2a), where a is a real number. () Compute the coordinates
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 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 informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More 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 informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationA positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations
A positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations Sashank Srinivasan a, Jonathan Poggie a, Xiangxiong Zhang b, a School of Aeronautics and Astronautics,
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationMATH 6B Spring 2017 Vector Calculus II Study Guide Final Exam Chapters 8, 9, and Sections 11.1, 11.2, 11.7, 12.2, 12.3.
MATH 6B pring 2017 Vector Calculus II tudy Guide Final Exam Chapters 8, 9, and ections 11.1, 11.2, 11.7, 12.2, 12.3. Before starting with the summary of the main concepts covered in the quarter, here there
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
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 informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
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 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 informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationTHE WAVE EQUATION. d = 1: D Alembert s formula We begin with the initial value problem in 1 space dimension { u = utt u xx = 0, in (0, ) R, (2)
THE WAVE EQUATION () The free wave equation takes the form u := ( t x )u = 0, u : R t R d x R In the literature, the operator := t x is called the D Alembertian on R +d. Later we shall also consider the
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
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 informationMathematics (Course B) Lent Term 2005 Examples Sheet 2
N12d Natural Sciences, Part IA Dr M. G. Worster Mathematics (Course B) Lent Term 2005 Examples Sheet 2 Please communicate any errors in this sheet to Dr Worster at M.G.Worster@damtp.cam.ac.uk. Note that
More informationMath 234 Final Exam (with answers) Spring 2017
Math 234 Final Exam (with answers) pring 217 1. onsider the points A = (1, 2, 3), B = (1, 2, 2), and = (2, 1, 4). (a) [6 points] Find the area of the triangle formed by A, B, and. olution: One way to solve
More informationA proof for the full Fourier series on [ π, π] is given here.
niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More 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 informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
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 informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
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 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 information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationThe Wiener Itô Chaos Expansion
1 The Wiener Itô Chaos Expansion The celebrated Wiener Itô chaos expansion is fundamental in stochastic analysis. In particular, it plays a crucial role in the Malliavin calculus as it is presented in
More informationQualifying Examination
Summer 24 Day. Monday, September 5, 24 You have three hours to complete this exam. Work all problems. Start each problem on a All problems are 2 points. Please email any electronic files in support of
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 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 information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationMA 441 Advanced Engineering Mathematics I Assignments - Spring 2014
MA 441 Advanced Engineering Mathematics I Assignments - Spring 2014 Dr. E. Jacobs The main texts for this course are Calculus by James Stewart and Fundamentals of Differential Equations by Nagle, Saff
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
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 informationAn Exponential High-Order Compact ADI Method for 3D Unsteady Convection Diffusion Problems
An Exponential High-Order Compact ADI Method for 3D Unsteady Convection Diffusion Problems Yongbin Ge, 1 Zhen F. Tian, 2 Jun Zhang 3 1 Institute of Applied Mathematics and Mechanics, Ningxia University,
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More informationIMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS
Approved for public release; distribution is unlimited IMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS Sean L. Gibbons Captain, United
More informationBasics of Discretization Methods
Basics of Discretization Methods In the finite difference approach, the continuous problem domain is discretized, so that the dependent variables are considered to exist only at discrete points. Derivatives
More informationMath 575-Lecture Fundamental solutions for point force: Stokeslet, stresslet
Math 575-Lecture 2 Fundamental solutions for point force: tokeslet, stresslet Consider putting a point singular force with strength F in a tokes flow of infinite extent. We then need to solve P µ U = δ(x)f,
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationConvergence of Summation-by-Parts Finite Difference Methods for the Wave Equation
J Sci Comput (2) 71:219 245 DOI 1./s1915-16-297-3 Convergence of Summation-by-Parts Finite Difference Methods for the Wave Equation Siyang Wang 1 Gunilla Kreiss 1 Received: 15 January 216 / Revised: 16
More informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
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 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 informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
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 informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationFourth Order Convergence of Compact Finite Difference Solver for 2D Incompressible Flow
Fourth Order Convergence of Compact Finite Difference Solver for D Incompressible Flow Cheng Wang 1 Institute for Scientific Computing and Applied Mathematics and Department of Mathematics, Indiana University,
More informationModule Two: Differential Calculus(continued) synopsis of results and problems (student copy)
Module Two: Differential Calculus(continued) synopsis of results and problems (student copy) Srikanth K S 1 Syllabus Taylor s and Maclaurin s theorems for function of one variable(statement only)- problems.
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 informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
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 informationMathematics 324 Riemann Zeta Function August 5, 2005
Mathematics 324 Riemann Zeta Function August 5, 25 In this note we give an introduction to the Riemann zeta function, which connects the ideas of real analysis with the arithmetic of the integers. Define
More information1 + lim. n n+1. f(x) = x + 1, x 1. and we check that f is increasing, instead. Using the quotient rule, we easily find that. 1 (x + 1) 1 x (x + 1) 2 =
Chapter 5 Sequences and series 5. Sequences Definition 5. (Sequence). A sequence is a function which is defined on the set N of natural numbers. Since such a function is uniquely determined by its values
More information3.4. Monotonicity of Advection Schemes
3.4. Monotonicity of Advection Schemes 3.4.1. Concept of Monotonicity When numerical schemes are used to advect a monotonic function, e.g., a monotonically decreasing function of x, the numerical solutions
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
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 informationNumerical Analysis of Differential Equations Numerical Solution of Elliptic Boundary Value
Numerical Analysis of Differential Equations 188 5 Numerical Solution of Elliptic Boundary Value Problems 5 Numerical Solution of Elliptic Boundary Value Problems TU Bergakademie Freiberg, SS 2012 Numerical
More informationSolutions for Euler and Navier-Stokes Equations in Finite and Infinite Series of Time
Solutions for Euler and Navier-Stokes Equations in Finite and Infinite Series of Time Valdir Monteiro dos Santos Godoi valdir.msgodoi@gmail.com Abstract We present solutions for the Euler and Navier-Stokes
More informationChapter 2 Boundary and Initial Data
Chapter 2 Boundary and Initial Data Abstract This chapter introduces the notions of boundary and initial value problems. Some operator notation is developed in order to represent boundary and initial value
More informationMAT 211 Final Exam. Spring Jennings. Show your work!
MAT 211 Final Exam. pring 215. Jennings. how your work! Hessian D = f xx f yy (f xy ) 2 (for optimization). Polar coordinates x = r cos(θ), y = r sin(θ), da = r dr dθ. ylindrical coordinates x = r cos(θ),
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 informationFundamentals Physics
Fundamentals Physics And Differential Equations 1 Dynamics Dynamics of a material point Ideal case, but often sufficient Dynamics of a solid Including rotation, torques 2 Position, Velocity, Acceleration
More informationBoundary conditions and estimates for the linearized Navier-Stokes equations on staggered grids
Boundary conditions and estimates for the linearized Navier-Stokes equations on staggered grids Wendy Kress and Jonas Nilsson Department of Scientific Computing Uppsala University Box S-75 4 Uppsala Sweden
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationThe Divergence Theorem Stokes Theorem Applications of Vector Calculus. Calculus. Vector Calculus (III)
Calculus Vector Calculus (III) Outline 1 The Divergence Theorem 2 Stokes Theorem 3 Applications of Vector Calculus The Divergence Theorem (I) Recall that at the end of section 12.5, we had rewritten Green
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More information1. (30 points) In the x-y plane, find and classify all local maxima, local minima, and saddle points of the function. f(x, y) = 3y 2 2y 3 3x 2 + 6xy.
APPM 35 FINAL EXAM FALL 13 INSTUTIONS: Electronic devices, books, and crib sheets are not permitted. Write your name and your instructor s name on the front of your bluebook. Work all problems. Show your
More informationA.1 THE SAMPLED TIME DOMAIN AND THE Z TRANSFORM. 0 δ(t)dt = 1, (A.1) δ(t)dt =
APPENDIX A THE Z TRANSFORM One of the most useful techniques in engineering or scientific analysis is transforming a problem from the time domain to the frequency domain ( 3). Using a Fourier or Laplace
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
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 informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
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 information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
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 informationDomain decomposition schemes with high-order accuracy and unconditional stability
Domain decomposition schemes with high-order accuracy and unconditional stability Wenrui Hao Shaohong Zhu March 7, 0 Abstract Parallel finite difference schemes with high-order accuracy and unconditional
More informationChapter 5: Bases in Hilbert Spaces
Orthogonal bases, general theory The Fourier basis in L 2 (T) Applications of Fourier series Chapter 5: Bases in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University
More information