Domain decomposition schemes with high-order accuracy and unconditional stability
|
|
- Lewis Lamb
- 6 years ago
- Views:
Transcription
1 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 stability for solving parabolic equations are presented. The schemes are based on domain decomposition method, i.e., interface values between subdomains are computed by the explicit scheme; interior values are computed by the implicit scheme. The numerical stability and error are derived in the H norm in one dimensional case. Numerical results of both one and two dimensions examining the stability, accuracy, and parallelism of the procedure are also presented. Keywords: Domain decomposition, Finite difference, Parabolic equation, High-order accuracy, Unconditional stability. Introduction Domain decomposition is a powerful tool for devising parallel methods to solve time-dependent partial differential equations. There is rich literature on domain decomposition finite difference methods for solving parabolic equations on parallel computers. For the non-overlapping domain decomposition methods, the explicit nature of the calculation at the interface of sub-domain leads some domain decomposition schemes to be conditionally stable, which implies that they have to suffer from temporal step-size restrictions (see []-[5]). Schemes with unconditional stability as well as high-order accuracy being desired in the applications, many investigators have turned to improve the stability of the domain decomposition method. For example, the corrected explicit-implicit domain decomposition algorithms were presented in [6] and [7]. By adding the correction step to explicit-implicit domain decomposition methods, updating the interface solutions at each time level, the corrected methods were proved to Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, IN USA (whao@nd.edu). School of Mathematics Science and LPMC, Nanai University, Tianin, 30007, China (shhzhu@nanai.edu.cn).
2 have the unconditional stability. While the needless corrected domain decomposition schemes with unconditional stability were presented in [8] and [9]. All of these methods with unconditional stability reach the second order accuracy at most. The purpose of this paper is to present the domain decomposition finite difference procedure with third-order accuracy and unconditional stability. We first consider the following Dirichlet boundary problem U t U = 0, x (0, ), t (0, T], x (.) U(0, t) = U(, t) = 0, t (0, T], (.) U(x, 0) = U 0 (x), 0 x, (.3) where the initial function U 0 (x) satisfies the boundary condition, i.e., U 0 (0) = U 0 () = 0. Then we extend the method to the problem of two dimensional space. We will introduce two new finite difference schemes for solving (.)-(.3) in Section, and setch the domain decomposition procedure, the numerical stability and convergence in Section 3. In section 4, the proof of the unconditional stability and the error estimate will be given. Numerical examples and examination of the algorithm will be provided in Section 5. In Sections 6 and 7, we extend the method to the problem of two dimensional space and test some examples. Two new finite difference schemes Taing the usual h, τ mesh in x and t, and denoting the approximate value of U(x, t n ) U n by u n, where x = h and t n = nτ, we define the following operators + u n = un + un, u n = un un, τu n = τ (un un ). It is well nown that the following Taylor expansion resulting in the fully implicit finite difference scheme is valid. τ U + U ( ) τ h = + h U t (x, t ) + O(τ + h 4 ). (.4) Noticing that U t (x, t ) = U U n + Un τ + O(τ), (.5) and substituting U t (x, t ) into (.4), we obtain τ U + U ( ) τ h + + h U U n + Un τ = O(τ + h 4 + τh ). (.6)
3 3 Replacing by + and in (.5), we have and U+ = Un + Un + + τ U t (x +, t ) + O(τ 3 ) U = Un Un + τ U t (x, t ) + O(τ 3 ). Substituting U± into (.6), we can get τ U Un + Un + U + U n Un h + ( = τ U h t (x +, t ) + U t (x, t ) ( ) τ + h U ) + O(τ + h 4 + τh + τ3 U n + Un τ h ). (.7) By omitting the high order term, (.6) and (.7) yield two new finite difference schemes for (.) : τ u + u ( τ h + + h ( τ + h τ u + u h + ) u ) u u n + un τ = 0, (.8) u n + un τ r( τ u n + τ u + ) r( τu n τ u ) = 0, (.9) where r = τ. From the derivation of (.8) and (.9), we now that the truncation errors of (.8) and (.9) are O(τ +h 4 +τh h ( ) τ ) and O h +τ +h 4 +τh + τ3 h respectively. If r is any positive real number, the truncation errors become O(h 4 ) and O(h ). 3 Domain decomposition procedure and main results Suppose Jh =, Nτ = T. For simplicity, we will consider a domain decomposition which involves in decomposing (0, ) into only two subdomains, (0, x) and (x, ), where x = x for some integer ( < < J ). We use the explicit scheme (.9) to compute the solution value u and the implicit scheme (.8) to compute other solution values u ( ) respectively. The system can be written as L(u ) = 0, n N, 0 < < J, u = 0, = 0, J, (3.0)
4 4 where the linear operator L is L(u ) = τ u + u ( ) τ h + + h u u n + un τ,, τ u + u ( ) τ h + + h u u n + un τ r( τ u n + τ u + ) r( τu n τ u ), =. (3.). The resulting system of equations decouples into two disoint sets of equations corresponding to the subdomains. These systems can be solved in parallel. We are now in position to state two main theorems of this paper, which will be proved in the forthcoming sections. For the discrete function u n = {u n = 0,,..., J, u n 0 = un J = 0}, define u n = u n h, + u n = + u n h. = We have the following theorems. Theorem 3. (Stability) For any given r > 0, the finite difference solutions of the parallel scheme (3.0)-(3.) satisfy max +u n n =0 r + r + ( +u + u u 0 ). Theorem 3. (Convergence) Let e n = U(x, t n ) u n. For any given r > 0, the finite difference solutions of the parallel scheme (3.0)-(3.) satisfy max n +e n C( + e + e e 0 + h 4 ), where C is a positive constant independent of h and τ. Since the finite difference scheme (3.0) has three time levels, besides taing u 0 = U 0(h), thus e 0 = 0, we need to find other methods to solve u,. In order to match high order accuracy, we can use either fourth order explicit schemes such as the impact scheme or the high-order parallel iterative method [0] to compute u, i.e., let e O(h4 ). Then according to Theorem 3., scheme (3.0) will reach third-order accuracy. 4 Proof of stability and convergence We first state three auxiliary lemmas. The stability and convergence results are then derived.
5 5 Lemma 4. (Discrete Poincare Inequality) For the discrete function u n = {u n = 0,,..., J, un 0 = u n J = 0}, there exists u n + u n h. Lemma 4. (Discrete Green Theorem) If u and v are discrete functions on the set {x = 0,,, J}, then we have J u + v = v u u 0 v + u J v J. = = Lemma 4.3 For any given r > 0, f L (0, T; L (0, )), the finite difference solutions of system L(u ) = f with the Dirichlet boundary condition satisfy ( max +u n + u + n r + ) u r + u 0 + Th max r n fn. Lemma 4. and Lemma 4. are proved in []. Before proving the stability and convergence, we will give the proof of Lemma Proof of Lemma 4.3 Denoting w = u boundary condition as w r + u u n, we can rewrite L(u ) = f with the Dirichlet ( + w r + u + + ) ( r + r (w ) (w w n ) r(w+ n w + ) r(wn w ) = f u = 0, = 0, J. w n ) = f τ,, τ, =, Multiplying the above equations by w h, =,,, J and summing them up respect to, we have = (w rw ) h r = ( w + u h + + ) r (w+ n w + )h rw (w n w From Lemma 4., we have = )h = = (w w n )w h f w τh. (4.)
6 6 J w + u = u w = = = ( + u ) + ( + u n ) ( + u + u n ). =0 Then (4.) becomes w + r [ =0 ( + + )[ r ( + u = (w +rhw (w + + w ) = rhw =0 (w+ n + wn ) + ) h =0 =0 ] ( + u n ) h + r ) h (w n ) h + = w f = ( (w ) + (w+ n rh ) + (w τh ) + (w n ) = ) + = =0 ( + u (w ] w n ) h (w ) + (f + u n ) h τ) rh(w ) + rh (wn + ) + (w n ) + w + τ f. (4.3) Noting that = r r ( + u + u n ) h + rhw (w + + w =0 =0,, ( + w we can simplify (4.3) as follows w + r ( r [ (w + +rh ) + (w ) [ (w ) + h + rh ) + (w [ + u + u n ] + r ) ) ( w w n + w w n ) ) =0,, ] + rh(w ), (wn + ) + (w n ] ) τ f ( + w ) h h
7 7 Since w 0, we have =0,, ] ( r [ + u + u n + [ (w + +rh ) + (w ) Summing up respect to n, we get Thus ( r +u r +u + ( 4 + 4r + u + u + Lemma 4.3 is proved. 4. Proof of stability ( + w ) h 0, and w w n 0, then 4 + 4r (wn + ) + (w n ) ) ( w w n ) ] τ f ) w + rh (w + ) + (w ) 4r ) w + rh (w + ) + (w ) + Tτ max n fn. ( r + ) u r + u 0 + Th max r n fn. Let f = 0 in Lemma 4.3, thus Theorem 3. is proved. 4.3 Proof of convergence According to the parallel scheme (3.0)-(3.), the errors e ( n < N) follows the below equations: τ e + e ( τ h + + h τ e + e ( τ h + + h ) e ) e e n + en τ = G,, e n + en τ r( τ e n + τe + ) r( τe n τe ) = Φ h + G, =, e = 0, = 0, J, where Φ = r U t (x, t ), G O(h 4 ). In order to get the error estimates for e, we assume that e = p + q, where p and q are the solutions of the following problems respectively.
8 8 problem I τ p + p ( τ h + + h τ p + p ( τ h + + h ) p ) p p n + pn τ = G,, p n + pn τ r( τ p n + τp + ) r( τp n τp ) = G, =, (4.4) p = 0, = 0, J, p 0 = e0, p = e,. problem II τ q + q ( τ h + + h τ q + q ( τ h + + h ) q ) q q n + qn τ = 0,, q n + qn τ r( τ q n + τ q + ) r( τq n τ q ) = Φ h, =, q = 0, = 0, J, q 0 = q = 0,. From lemma 4.3, we can obtain the estimate for p, i.e., ( + p n + p + r + ) p r + p 0 + Th max r n Gn ( r + ) r + ( + e + e e 0 ) + O(h 0 ). (4.5) For estimating q, we first consider q ( n N )which satisfies the following equations Then the formula of q q = + q h = 0,, + q h = Φ h, q 0 = q J = 0. is 0, = 0, J, J Φ h 4,, J + J J Φ h 4, < < J. i=
9 9 Hence + q O(h 4 ) and q q n O(τh 3 ). Define q = q q, then we have τ q + q ( τ h + + h τ q + q ( τ h + + h ) q ) q q n + qn τ = R,, q n + qn τ r( τ q n + τ q + ) r( τ q n τ q ) = R, =, q = 0, = 0, J, q 0 = q 0, q = q,, where R = ( τ τ q + + h ( τ τ q + + h r( τ q n τ q ) q ) q ), =, q n + qn τ,, q n + q n τ r( τ q n + τ q + ) and R O(h 3 ). From Lemma 4.3, we can obtain the estimate of + q. ( + q + q + r + ) q r + q 0 + Th max r n Rn O(h 8 ). Thus, Combining with (4.5), we get + q + q + + q O(h 4 ). + e + p + + q C( + e + e e 0 + h 4 ), where C is a positive constant independent of h and τ. Theorem 3. is proved. 5 Numerical experiments In this section, some numerical results are presented to show the stability, accuracy, and parallelism of the scheme described above, and the computational costs are also presented. All the experiments are run on a cluster consisting of a manager that uses one core of a Xeon 540 processor and up to 3 computing nodes, each containing two Xeon 540 processors running 64-bit Linux, i.e., each node consists of 8 processing cores.
10 0 We consider the problem defined in equations (.)-(.3) with U 0 (x) = sin(πx). Obviously the exact solution of the equations is U(x, t) = e πt sin(πx). First, we verify the stability of the scheme by taing the step size h = 0 3, r =, 0, 00, 000 with two subdomains. FIG. clearly shows that the norm of u n doesn t occur blowing up even if r is large enough. This explains the unconditional stability of the scheme. r= r=0 u n 0.5 u n time t r= time t r=000 u n 0.5 u n time t time t Figure : The infinity norm of u n v.s. t Second, we examine the numerical errors in the solutions. Table shows that the errors for each case are roughly of the same order of magnitude, and the errors appear to be O(h 3 ) in each case. Third, we test the speed-up for the scheme. Here we tae h = 0 5, τ = 0 6, i.e., r = 0 4, T = 0.5 and list the time for computing and the speed-up in Table 5, which shows that the scheme has a coemptive parallelism. 6 Extension to two dimensional case In this section, U(x, y, t) will be a solution of the following Dirichlet boundary problem on Ω = (0, ) (0, ),
11 Table : Numerical errors for different grid points r = 0, T = processors 4 processors 8 processors J u n U n /h 3 J u n U n /h 3 J u n U n /h r = 00, T = processors 4 processors 8 processors J u n U n /h 3 J u n U n /h 3 J u n U n /h e e e e e e e e e e e e+0 r = 000, T = processors 4 processors 8 processors J u n U n /h 3 J u n U n /h 3 J u n U n /h e e e e e e e e e e e e+04 U U = 0, (x, y) Ω, t (0, T], t (6.6) U(x, y, t) = 0, (x, y) Ω, t (0, T], (6.7) U(x, y, 0) = U 0 (x, y), (x, y) Ω, (6.8) where the initial function U 0 (x, y) satisfies the boundary condition, and U = U x + U. We extend the domain decomposition method stated in Section 3 y to the above problem. Tae Ω = {(x, y) Ω : x < x}, Ω = {(x, y) Ω : x > x}. Let x i = ih same as in Section 3, and let y = h. Suppose that there exists an integer such that x = x (see Fig. ). In analogy with Section 3 we call points (x i, y, t n ) as boundary points if (x i, y ) Ω, and interface points if i =. Otherwise, we call them interior points. The values u n i, will approximate U(x i, y, t n ) Ui, n. We denote + u i, = u i+, u i, + u i,, δ + δ u i, = u i,+ u i, + u i,,
12 Table : Comparison of time for different processors Processors n Time t n (seconds) Speed-up t /t n Figure : domain decomposition for two dimensional case and Then Thus we have τ u n i, = un i, un i,. τ + U i, h = U x (x i, y, t) + h 4 U x 4 (x i, y, t) + O(h 4 ), δ + δ U i, h = U y (x i, y, t) + h 4 U y 4 (x i, y, t) + O(h 4 ), δ + δ + U i, h 4 = 4 U x y (x i, y, t) + O(h ). + U i, + δ + δ U i, h + 6h δ +δ + U i, = U(x i, y, t) + h ( 4 U x 4 (x i, y, t) + 4 U x y (x i, y, t) + 4 U y 4 (x i, y, t) ) + O(h 4 ) = U(x i, y, t) + h U t (x i, y, t) + O(h 4 ). (6.9)
13 3 Moreover, τ U i, = U t (x i, y, t ) τ U t (x i, y, t ) + O(τ ). (6.0) Subtract (6.0) with (6.9), we obtain τ Ui, + Ui, + δ + δ Ui, h 6h δ +δ + Ui, ( ) τ = + h U t (x i, y, t ) + O(τ + h 4 ). Therefore, the fourth order finite difference scheme for (6.6) is τ u i, + u i, + δ + δ u i, h 6h δ +δ + u i, ( ) τ + + h u i, u n i, + un i, τ = 0, (6.) which is a nine-point finite difference scheme with three levels (see Fig. 3). Figure 3: nine point scheme Denote u i, = u i, + u i,+ + u i+, + u i+,+ 4u i, = + u i, + δ + δ u i, + δ +δ + u i,. Then (6.) is equivalent to τ u i, u i, + + u i, + δ + δ u i, 3h + ( τ + h ) u i, u n i, + un i, τ = 0,
14 4 which can be used for computing u i, when i. When i =, u i±, in (6.) are replaced by u n i±, un i±,. In this case, we need to solve a tridiagonal matrix using Thomas algorithm for computing u,. It is straightforward to extend the two-subdomain results to many subdomains on x direction by cutting the whole domain into vertical strips. Moreover, we can also extend the domain decomposition scheme to three dimensional case by dividing x axis into many subdomains. 7 Two dimensional numerical experiments We consider the problem defined in equations (6.6)-(6.8) with U 0 (x, y) = sin(πx)sin(πy). Obviously the exact solution of the equations is U(x, y, t) = e πt sin(πx)sin(πy). The computing time and the speed-up of the scheme ( ) are shown in Table 3 with h = 0 3, τ = 0 and T =. Table 3: Comparison of time for two dimensional case Processors n Time t n (seconds) Speed-up t /t n 0h3m9s - 6h0m9s.7 4 h5m4s h3m54s m47s m3s m58s Lota-Volterra system The competitive Lota-Volterra equations are a simple model of the population dynamics of species competing for some common resource. Given two populations, u and v, with logistic dynamics, the Lota-Volterra formulation adds an additional term to account for the species interactions. Thus the competitive Lota-Volterra equations are: u t = u + u( v), (x, y) Ω v t = v v( u), (x, y) Ω v = u = 0, (x, y) Ω. We choose random values for u(x, y) and v(x, y) as initial conditions shown in FIG 4. We tae h = 0.00, τ = 0.0, T = and show the computing time in Table 4. FIG 5 shows the numerical solution for T =.
15 5 u(x,y,t) time t= y x v(x,y,t) time t= y x Figure 4: Initial condition for the Lota-Volterra system Table 4: Comparison of time for Lota-Volterra system Processors n Time t n (seconds) Speed-up t /t n h5m3s - h45m8s hm36s h58m57s h33m37s m8s ms m33s m8s 75.6 Conclusion We presented domain decomposition finite difference schemes with unconditional stability and third-order accuracy for the parabolic system. Error estimate and stability of the numerical solutions have been derived for one dimensional case. The scheme is easy to implement the parallelism and is extended in
16 6 y u(x,y,t) time t= x x v(x,y,t) time t= x y x Figure 5: Numerical solutions for T = two and three dimensional case. The numerical results demonstrate the good performance of the parallel scheme, namely, unconditional stability, the third order accuracy and high degree of parallelism. References [] C.N. Dawson, Q. Du, and T. F. Dupont, A finite difference domain decomposition algorithm for numerical solution of the Heat equation,math. Comp., Vol. 57, 63 7 (99). [] G. Yuan, S. Zhu and L. Shen, Domain decomposition algorithm based on the group explicit formula for the heat equation, Int. J. Comput. Math, Vol. 8, (005). [3] S. Zhu, Z. Yu and J. Zhao, A high-order parallel finite difference algorithm, Applied Math. and Comp., Vol. 83, (006). [4] C.N. Dawson and T.F. Dupont, Explicit/implicit conservative domain decomposition procedures for parabolic problems based on bloc-centered finite differences, SIAM J. Numer. Anal., Vol. 3, (994).
17 7 [5] G. Yuan and L. Shen, Stability and convergence of the explicit-implicit conservative domain decomposition procedure for parabolic problems, Computers and Math. with Applications, Vol. 47, (004). [6] H. Shi and H. Liao, Unconditional stability of corrected explicit-implicit domain decomposition algorithms for parallel approximation of heat equations, SIAM J. Numer. Anal., Vol. 44, (006). [7] H. Liao, H. Shi and Z. Sun, Corrected explicit-implicit domain decomposition algorithms for two-dimensional semilinear parabolic equations, Science in China Series A, Vol. 5, (009). [8] Z. Sheng, G. Yuan and X. Hang, Unconditional stability of parallel difference schemes with second order accuracy for parabolic equation, Applied Math. and Comp., Vol. 84, (007). [9] S. Zhu, Conservative domain decomposition procedure with unconditional stability and second-order accuracy, Applied Math. and Comp., Vol.6, (00). [0] W. Hao and S. Zhu, Parallel iterative methods for parabolic equations, Int. J. Comput. Math, Vol. 86, (009). [] Y. Zhou, Applications of discrete functional analysis to the finite difference method, International Academic publishers (99).
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 and Its Difference Approximations 1-D Initial Boundary Value
More informationA FINITE DIFFERENCE DOMAIN DECOMPOSITION ALGORITHM FOR NUMERICAL SOLUTION OF THE HEAT EQUATION
mathematics of computation volume 57, number 195 july 1991, pages 63-71 A FINITE DIFFERENCE DOMAIN DECOMPOSITION ALGORITHM FOR NUMERICAL SOLUTION OF THE HEAT EQUATION CLINT N. DAWSON, QIANG DU, AND TODD
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 informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
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 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 informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
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 informationA stencil of the finite-difference method for the 2D convection diffusion equation and its new iterative scheme
International Journal of Computer Mathematics Vol. 87, No. 11, September 2010, 2588 2600 A stencil of the finite-difference method for the 2D convection diffusion equation and its new iterative scheme
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 informationSECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 3 (06), pp. 9 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or
More informationParallel Galerkin Domain Decomposition Procedures for Parabolic Equation on General Domain
Parallel Galerkin Domain Decomposition Procedures for Parabolic Equation on General Domain Keying Ma, 1 Tongjun Sun, 1 Danping Yang 1 School of Mathematics, Shandong University, Jinan 50100, People s Republic
More informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
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 informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
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 informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More informationECE539 - Advanced Theory of Semiconductors and Semiconductor Devices. Numerical Methods and Simulation / Umberto Ravaioli
ECE539 - Advanced Theory of Semiconductors and Semiconductor Devices 1 General concepts Numerical Methods and Simulation / Umberto Ravaioli Introduction to the Numerical Solution of Partial Differential
More informationBlock-tridiagonal matrices
Block-tridiagonal matrices. p.1/31 Block-tridiagonal matrices - where do these arise? - as a result of a particular mesh-point ordering - as a part of a factorization procedure, for example when we compute
More informationA CCD-ADI method for unsteady convection-diffusion equations
A CCD-ADI method for unsteady convection-diffusion equations Hai-Wei Sun, Leonard Z. Li Department of Mathematics, University of Macau, Macao Abstract With a combined compact difference scheme for the
More informationConvergence Behavior of a Two-Level Optimized Schwarz Preconditioner
Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Olivier Dubois 1 and Martin J. Gander 2 1 IMA, University of Minnesota, 207 Church St. SE, Minneapolis, MN 55455 dubois@ima.umn.edu
More informationAn Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations.
An Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations by Tong Chen A thesis submitted in conformity with the requirements
More informationAn additive average Schwarz method for the plate bending problem
J. Numer. Math., Vol. 10, No. 2, pp. 109 125 (2002) c VSP 2002 Prepared using jnm.sty [Version: 02.02.2002 v1.2] An additive average Schwarz method for the plate bending problem X. Feng and T. Rahman Abstract
More informationCode: 101MAT4 101MT4B. Today s topics Finite-difference method in 2D. Poisson equation Wave equation
Code: MAT MTB Today s topics Finite-difference method in D Poisson equation Wave equation Finite-difference method for elliptic PDEs in D Recall that u is the short version of u x + u y Dirichlet BVP:
More informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
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 in a 2D Box Region Let us consider a model problem of parabolic
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationON APPROXIMATION OF LAPLACIAN EIGENPROBLEM OVER A REGULAR HEXAGON WITH ZERO BOUNDARY CONDITIONS 1) 1. Introduction
Journal of Computational Mathematics, Vol., No., 4, 75 86. ON APPROXIMATION OF LAPLACIAN EIGENPROBLEM OVER A REGULAR HEXAGON WITH ZERO BOUNDARY CONDITIONS ) Jia-chang Sun (Parallel Computing Division,
More informationA fourth-order finite difference scheme for the numerical solution of 1D linear hyperbolic equation
203 (203) - Available online at www.ispacs.com/cna Volume 203, Year 203 Article ID cna-0048, Pages doi:0.5899/203/cna-0048 Research Article A fourth-order finite difference scheme for the numerical solution
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 informationEfficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems
Efficient smoothers for all-at-once multigrid methods for Poisson and Stoes control problems Stefan Taacs stefan.taacs@numa.uni-linz.ac.at, WWW home page: http://www.numa.uni-linz.ac.at/~stefant/j3362/
More informationDomain decomposition for the Jacobi-Davidson method: practical strategies
Chapter 4 Domain decomposition for the Jacobi-Davidson method: practical strategies Abstract The Jacobi-Davidson method is an iterative method for the computation of solutions of large eigenvalue problems.
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationSynchronization and Bifurcation Analysis in Coupled Networks of Discrete-Time Systems
Commun. Theor. Phys. (Beijing, China) 48 (2007) pp. 871 876 c International Academic Publishers Vol. 48, No. 5, November 15, 2007 Synchronization and Bifurcation Analysis in Coupled Networks of Discrete-Time
More information= ( 1 P + S V P) 1. Speedup T S /T V p
Numerical Simulation - Homework Solutions 2013 1. Amdahl s law. (28%) Consider parallel processing with common memory. The number of parallel processor is 254. Explain why Amdahl s law for this situation
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 informationA SEMI-LAGRANGIAN RUNGE-KUTTA METHOD FOR TIME-DEPENDENT PARTIAL DIFFERENTIAL EQUATIONS
Journal of pplied nalysis and Computation Website:http://jaac-online.com/ Volume 3, Number 3, ugust 23 pp. 25 263 SEMI-LGRNGIN RUNGE-KUTT METHOD FOR TIME-DEPENDENT PRTIL DIFFERENTIL EQUTIONS Daniel X.
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 informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationRandom Graphs. EECS 126 (UC Berkeley) Spring 2019
Random Graphs EECS 126 (UC Bereley) Spring 2019 1 Introduction In this note, we will briefly introduce the subject of random graphs, also nown as Erdös-Rényi random graphs. Given a positive integer n and
More informationMulti-Factor Finite Differences
February 17, 2017 Aims and outline Finite differences for more than one direction The θ-method, explicit, implicit, Crank-Nicolson Iterative solution of discretised equations Alternating directions implicit
More informationFinite Difference Method of Fractional Parabolic Partial Differential Equations with Variable Coefficients
International Journal of Contemporary Mathematical Sciences Vol. 9, 014, no. 16, 767-776 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.1988/ijcms.014.411118 Finite Difference Method of Fractional Parabolic
More informationPARAREAL TIME DISCRETIZATION FOR PARABOLIC CONTROL PROBLEM
PARAREAL TIME DISCRETIZATION FOR PARABOLIC CONTROL PROBLEM Daoud S Daoud Dept. of Mathematics Eastern Mediterranean University Famagusta-North Cyprus via Mersin 10-Turkey. e mail: daoud.daoud@emu.edu.tr
More informationNatural Boundary Integral Method and Its Applications
Natural Boundary Integral Method and Its Applications By De-hao Yu State Key Laboratory of Scientific and Engineering Computing Institute of Computational Mathematics and Scientific/Engineering Computing
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
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 informationarxiv: v1 [math.na] 11 Jul 2018
An Additive Overlapping Domain Decomposition Method for the Helmholtz Equation Wei Leng Lili Ju arxiv:1807.04180v1 [math.na] 11 Jul 2018 July 12, 2018 Abstract In this paper, we propose and analyze an
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationA FOURTH-ORDER ACCURATE FINITE DIFFERENCE SCHEME FOR THE EXTENDED-FISHER-KOLMOGOROV EQUATION. Tlili Kadri and Khaled Omrani
Bull. Korean Math. Soc. 0 (0), No. 0, pp. 0 https://doi.org/0.434/bkms.b603 pissn: 05-8634 / eissn: 34-306 A FOURTH-ORDER ACCURATE FINITE DIFFERENCE SCHEME FOR THE EXTENDED-FISHER-KOLMOGOROV EQUATION Tlili
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
More informationKrylov Implicit Integration Factor Methods for Semilinear Fourth-Order Equations
mathematics Article Krylov Implicit Integration Factor Methods for Semilinear Fourth-Order Equations Michael Machen and Yong-Tao Zhang * Department of Applied and Computational Mathematics and Statistics,
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 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 informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationA local Crank-Nicolson method for solving the heat equation. Abdurishit ABUDUWALI, Michio SAKAKIHARA and Hiroshi NIKI
HIROSHIMA MATH. J. 24 (1994), 1-13 A local Crank-Nicolson method for solving the heat equation Abdurishit ABUDUWALI, Michio SAKAKIHARA and Hiroshi NIKI (Received February 6, 1992) 1. Introduction du A
More informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
More informationA Fast High Order Algorithm for 3D Helmholtz Equation with Dirichlet Boundary
Applied and Computational Mathematics 2018; 7(4): 180-187 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20180704.11 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Fast High Order
More informationNumerical Algorithms for Visual Computing II 2010/11 Example Solutions for Assignment 6
Numerical Algorithms for Visual Computing II 00/ Example Solutions for Assignment 6 Problem (Matrix Stability Infusion). The matrix A of the arising matrix notation U n+ = AU n takes the following form,
More informationNumerical methods for weak solution of wave equation with van der Pol type boundary conditions
Numerical methods for weak solution of wave equation with van der Pol type boundary conditions Jun Liu a, Yu Huang b, Haiwei Sun c, Mingqing Xiao a, a Department of Mathematics, Southern Illinois University,
More informationCONVERGENCE OF GAUGE METHOD FOR INCOMPRESSIBLE FLOW CHENG WANG AND JIAN-GUO LIU
MATHEMATICS OF COMPUTATION Volume 69, Number 232, Pages 135{1407 S 0025-571(00)0124-5 Article electronically published on March 24, 2000 CONVERGENCE OF GAUGE METHOD FOR INCOMPRESSIBLE FLOW CHENG WANG AND
More informationA parallel algorithm for the heat equation with derivative boundary conditions
International Mathematical Forum, 2, 2007, no. 12, 565-574 A parallel algorithm for the heat equation with derivative boundary conditions M. Akram University College of Information Technology University
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationPreconditioning for Nonsymmetry and Time-dependence
Preconditioning for Nonsymmetry and Time-dependence Andy Wathen Oxford University, UK joint work with Jen Pestana and Elle McDonald Jeju, Korea, 2015 p.1/24 Iterative methods For self-adjoint problems/symmetric
More informationAdditive Schwarz Methods for Hyperbolic Equations
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03044-2 Additive Schwarz Methods for Hyperbolic Equations Yunhai Wu, Xiao-Chuan Cai, and David E. Keyes 1. Introduction In recent years, there
More informationThe Milne error estimator for stiff problems
13 R. Tshelametse / SAJPAM. Volume 4 (2009) 13-28 The Milne error estimator for stiff problems Ronald Tshelametse Department of Mathematics University of Botswana Private Bag 0022 Gaborone, Botswana. E-mail
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 informationCompact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems
Compact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems Daniel Ritter 1, Ulrich Rüde 1, Björn Gmeiner 1, Rochus Schmid 2 Copper Mountain, March 18th, 2013
More informationExponentially Fitted Error Correction Methods for Solving Initial Value Problems
KYUNGPOOK Math. J. 52(2012), 167-177 http://dx.doi.org/10.5666/kmj.2012.52.2.167 Exponentially Fitted Error Correction Methods for Solving Initial Value Problems Sangdong Kim and Philsu Kim Department
More informationFinite Difference Method for PDE. Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36
Finite Difference Method for PDE Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36 1 Classification of the Partial Differential Equations Consider a scalar second order partial
More informationAn interior-point trust-region polynomial algorithm for convex programming
An interior-point trust-region polynomial algorithm for convex programming Ye LU and Ya-xiang YUAN Abstract. An interior-point trust-region algorithm is proposed for minimization of a convex quadratic
More informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More informationUltraconvergence of ZZ Patch Recovery at Mesh Symmetry Points
Ultraconvergence of ZZ Patch Recovery at Mesh Symmetry Points Zhimin Zhang and Runchang Lin Department of Mathematics, Wayne State University Abstract. The ultraconvergence property of the Zienkiewicz-Zhu
More informationBoundary Value Problems and Approximate Solutions
Boundary Value Problems and Approximate Solutions Gebreslassie Tesfy, Venketeswara Rao J*, Ataklti Araya and Daniel Tesfay Department of Mathematics, College of Natural and Computational Scineces, Mekelle
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
More informationSolving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions.
Journal of Mathematical Modeling Vol 1, No 1, 213, pp 28-4 JMM Solving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions Farshid
More informationPerformance Analysis of Parallel Alternating Directions Algorithm for Time Dependent Problems
Performance Analysis of Parallel Alternating Directions Algorithm for Time Dependent Problems Ivan Lirkov 1, Marcin Paprzycki 2, and Maria Ganzha 2 1 Institute of Information and Communication Technologies,
More informationError formulas for divided difference expansions and numerical differentiation
Error formulas for divided difference expansions and numerical differentiation Michael S. Floater Abstract: We derive an expression for the remainder in divided difference expansions and use it to give
More informationOn the efficiency of the Peaceman-Rachford ADI-dG method for wave-type methods
On the efficiency of the Peaceman-Rachford ADI-dG method for wave-type methods Marlis Hochbruck, Jonas Köhler CRC Preprint 2017/34 (revised), March 2018 KARLSRUHE INSTITUTE OF TECHNOLOGY KIT The Research
More informationAnalysis of Two-Grid Methods for Nonlinear Parabolic Equations by Expanded Mixed Finite Element Methods
Advances in Applied athematics and echanics Adv. Appl. ath. ech., Vol. 1, No. 6, pp. 830-844 DOI: 10.408/aamm.09-m09S09 December 009 Analysis of Two-Grid ethods for Nonlinear Parabolic Equations by Expanded
More informationCRANK-NICOLSON FINITE DIFFERENCE METHOD FOR SOLVING TIME-FRACTIONAL DIFFUSION EQUATION
Journal of Fractional Calculus and Applications, Vol. 2. Jan. 2012, No. 2, pp. 1-9. ISSN: 2090-5858. http://www.fcaj.webs.com/ CRANK-NICOLSON FINITE DIFFERENCE METHOD FOR SOLVING TIME-FRACTIONAL DIFFUSION
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More informationDraft TAYLOR SERIES METHOD FOR SYSTEM OF PARTICLES INTERACTING VIA LENNARD-JONES POTENTIAL. Nikolai Shegunov, Ivan Hristov
TAYLOR SERIES METHOD FOR SYSTEM OF PARTICLES INTERACTING VIA LENNARD-JONES POTENTIAL Nikolai Shegunov, Ivan Hristov March 9, Seminar in IMI,BAS,Soa 1/25 Mathematical model We consider the Hamiltonian :
More informationarxiv: v1 [math.na] 15 Nov 2017
Noname manuscript No. (will be inserted by the editor) An HDG method with orthogonal projections in facet integrals Issei Oikawa arxiv:1711.05396v1 [math.na] 15 Nov 2017 Received: date / Accepted: date
More informationMATH 425, HOMEWORK 5, SOLUTIONS
MATH 425, HOMEWORK 5, SOLUTIONS Exercise (Uniqueness for the heat equation on R) Suppose that the functions u, u 2 : R x R t R solve: t u k 2 xu = 0, x R, t > 0 u (x, 0) = φ(x), x R and t u 2 k 2 xu 2
More 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 informationDYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS. Wei Feng
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS SUPPLEMENT 7 pp. 36 37 DYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS Wei Feng Mathematics and Statistics Department
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationDOMAIN DECOMPOSITION APPROACHES FOR MESH GENERATION VIA THE EQUIDISTRIBUTION PRINCIPLE
DOMAIN DECOMPOSITION APPROACHES FOR MESH GENERATION VIA THE EQUIDISTRIBUTION PRINCIPLE MARTIN J. GANDER AND RONALD D. HAYNES Abstract. Moving mesh methods based on the equidistribution principle are powerful
More information