Finite volume element method for analysis of unsteady reaction diffusion problems
|
|
- Bruce Lynch
- 6 years ago
- Views:
Transcription
1 Acta Mech Sin (2009) 25: DOI /s RESEARCH PAPER Finite volume element method for analysis of unsteady reaction diffusion problems Sutthisak Phongthanapanich Pramote Dechaumphai Received: 6 May 2008 / Revised: 6 October 2008 / Accepted: 5 December 2008 / Published online: 27 February 2009 The Chinese Society of Theoretical and Applied Mechanics and Springer-Verlag GmbH 2009 Abstract A finite volume element method is developed for analyzing unsteady scalar reaction diffusion problems in two dimensions. The method combines the concepts that are employed in the finite volume and the finite element method together. The finite volume method is used to discretize the unsteady reaction diffusion equation, while the finite element method is applied to estimate the gradient quantities at cell faces. Robustness and efficiency of the combined method have been evaluated on uniform rectangular grids by using available numerical solutions of the two-dimensional reaction diffusion problems. The numerical solutions demonstrate that the combined method is stable and can provide accurate solution without spurious oscillation along the highgradient boundary layers. Keywords Finite volume element method Explicit method Unsteady problem Singularly perturbed equation Reaction diffusion 1 Introduction Recent development of several reaction diffusion analysis schemes in biology, ecology, biochemistry, fluid mechanics S. Phongthanapanich Department of Mechanical Engineering Technology, College of Industrial Technology, King Mongkut s University of Technology North Bangkok, Bangkok 10800, Thailand sutthisakp@kmutnb.ac.th P. Dechaumphai (B) Mechanical Engineering Department, Chulalongkorn University, Bangkok 10330, Thailand fmepdc@eng.chula.ac.th and heat transfer has led to physically interesting and mathematically challenging problems related to the nontrivial numerical solution of the singularly perturbed initial/boundary value problem [1 9]. In general, solution of the singularly perturbed problem contains high-gradient boundary layers along the boundary of the domain [2]. Numerical approximation of such boundary layer flow behavior is not simple, and thus the development of accurate numerical method for analyzing this problem must be done carefully. Some examples of such difficulty are encountered in the cases where the reaction diffusion equation cannot be satisfied by a globally smooth function. In particular, if the unknown scalar variables vary rapidly within a thin layer, many standard numerical schemes may lead to inaccurate and unstable solutions [8 10]. In particular, a robust numerical method, in the sense that the error should not deteriorate as the singular perturbation parameter tends to zero, is required. During the past decade, many techniques of the finite element method for solving the singularly perturbed problems have been proposed and investigated [2,4,7,8,11,12]. The conventional Galerkin finite element method usually fails to analyze such unsteady reaction diffusion problems when either the diffusion coefficient becomes small or a small time step is employed in the time discretization. The stabilized finite element method overcame this trouble by employing the balanced residual perturbation term with some stabilization parameters. However, such method has great influence on the solution accuracy [13]. An improved finite element method so-called the Petrov Galerkin enriched method (PGEM), where the trial space is enriched by the multiscale functions, is later proposed by Franca et al. [12,14]. The method provides good solution accuracy as compared to the other finite element methods for solving the scalar reaction diffusion equation. Some other finite element methods, such as those use the Shishkin mesh, have also been
2 482 S. Phongthanapanich, P. Dechaumphai developed [2 4]. These methods apply robust algebraic error bounds, or employ the p/hp version of the finite element method on appropriate meshes [15,16]. At the same time, some researchers employed the finite difference algorithm to analyze the reaction diffusion problems [5,17]. Liao et al. proposed an efficient higher-order finite difference algorithm for solving the system of twodimensional reaction diffusion equations with nonlinear reaction term. The Padé approximation and Richardson extrapolation are used to achieve high-order solution accuracy in both the spatial and temporal domains. El Salam and Shehata [5] used the implicit Crank Nicolson technique to linearize the governing equations and then solved the system of algebraic equations by using the LU-decomposition method. Recently, many researchers proposed some mixed techniques of the finite element and finite volume methods for analyzing the parabolic equation [18 20]. Rui [18] presented a mixed covolume method for solving the parabolic equation from the general self-adjoint parabolic problem using a variable nondiagonal diffusion tensor and the lowest order Raviart Thomas mixed element. The basic idea of creating a covolume method is to seek for a good combination of the finite volume method and the MAC scheme. Ma et al. [19] proposed the semi-discrete and full-discrete symmetric finite volume schemes based on the linear finite elements for a class of parabolic problems. The symmetric finite volume schemes are analyzed on dual mesh (cell-vertex scheme) of triangular cells. Recently, Wang [20] proposed three schemes of the finite volume element method for solving two-dimensional parabolic equations. His first scheme is analogous to Douglas finite difference scheme with second-order splitting error; the second scheme has third-order splitting error; and the third scheme is an extended locally one-dimensional scheme. The finite volume schemes are also discretized on dual mesh and the nodal value at the center of control volume is approximated by linear finite element interpolation functions. It should be noted that the computational techniques for analyzing the hyperbolic equation are generally classified into the explicit and implicit (or semi-implicit) methods. The explicit method is a popular one because it is simple and requires less computer memory. However, the method is constrained by the CFL condition in order to stabilize the spatial error from growing without bound. Attempts to relief such constraint have been developed, such as one proposed by Wang and Liu [21]. On the other hand, the implicit method provides a more stable solution. However, large time step size may not be used because the solution accuracy degrades with time. The inversion of the coefficient matrix is another weakness of the method because it is a time consumable process. In addition, a large block of memory is required for the coefficient matrix formation. A finite volume element method, that combines the finite volume and the finite element methods, is proposed for analyzing the unsteady scalar reaction diffusion equation in this paper. The cell-centered finite volume method is used to analyze the problem on a uniform control volume. The concept of the finite element method is applied to estimate the gradient quantities herein by employing the method of weighted residuals to determine the nodal gradient quantities. This approach is different from those presented in Refs. [18 20], where the distribution of gradient quantities on element is simply approximated by using linear finite element interpolation functions. Finally, the trapezoidal rule is used to calculate the gradient quantities at cell faces. The robustness and accuracy of the combined method are examined by using reaction diffusion problems. The presentation of the paper starts from explanation of the theoretical formulation in Sect. 2. The proposed method is then evaluated in Sect. 3 by using three examples. The numerical results from these examples are used to evaluate the applicability of the combined method for analyzing both the steady and unsteady reaction diffusion problems. 2 Finite volume element formulation In this paper, the finite volume approach is used to discretize the following two-dimensional scalar reaction diffusion equation. φ ɛ 2 φ + κφ = f, in T, (1) t subject to the boundary conditions φ = g D, on D, (2) ɛ φ n = g N, on N, (3) with = D N with D N =. The initial condition is defined for x with R 2 by φ(x, 0) = φ 0 (x), where φ is the unknown scalar quantity, ɛ>0isthe diffusion coefficient, κ is the reaction coefficient, f = f (x) is the prescribed source term, and t (0, T ) for T <. The computational domain is first discretized into a collection of non-overlapping rectangular control volumes, i = 1,...,N, that completely cover the domain such that = N, =, and j =, if i = j. i=1 Equation (1) is then integrated over the arbitrary control volume to obtain ( ) φ ɛ 2 φ + κφ f dx = 0. (5) t (4)
3 Finite volume element method for analysis of unsteady reaction diffusion problems 483 The finite difference approximation and the divergence theorem are then applied to the temporal and spatial terms, respectively, of the above equation to yield φ(x, t n+1 )dx = φ(x, t n )dx + t ɛ κφ(x, t)dx + n i (ν) φ(ν,t)dν f (x)dx, (6) where n i (ν) is the unit normal vector of. In the conventional finite volume method, the approximation to the cell average of φ over control volume at time t n and t n+1 is represented by φi n+1 = 1 φ(x, t n+1 )dx, (7) φi n = 1 φ(x, t n )dx, where is the measure of. For an arbitrary control volume, the flux integral over appearing on the right-hand side of Eq. (6) could be approximated by the summation of fluxes passing through all the adjacent cell faces. Hence, by applying the quadrature integration formula on the spatial domain term, the flux integral over may be approximated by ɛ n i (ν) φ(ν,t)dν = ɛ 4 Ɣij nij φij n, (8) j=1 where Ɣ ij is the segment of boundary between the two adjacent cells and j, which is defined by = 4j=1 Ɣ ij and Ɣ ij = j. The integrations of the reaction and the source terms could be approximated by the mid-point quadrature rule to obtain κφ(x, t)dx = κφi n, (9) f (x)dx = f (x i ). (10) By substituting Eqs. (7) (10) into Eq. (6), a fully explicit finite volume scheme for solving Eq. (1) can be written in the form φ n+1 i = φ n i + ɛ t 4 Ɣij nij φij n t(κφn i f i ). j=1 (11) For unsteady computation, the second-order temporal accuracy is achieved by implementing the second-order accurate Runge Kutta time stepping method [22]. The treatment of boundary conditions in this paper is based on using the concept of ghost cell. The central differencing at the boundary is used for the Neumann boundary condition, while the Dirichlet boundary condition is prescribed at the center of the ghost cell. In this paper, the gradient term, φij n, is approximated by the weighted residuals method which is commonly used in the finite element technique [23]. Firstly, the φi n is assumed to distribute linearly over cell in the form φ n i = 4 N k (x) φk n, (12) k=1 where N k (x) denotes the linear interpolation functions for rectangular cell. By applying the weighted residuals method and the Gauss s theorem to Eq. (12), the gradient quantities at point J are obtained as φ n J,i = M 1 n i (ν)n J (ν)φi n dν N J (x) x φi n dx, (13) where M is the consistent mass matrix [23], and φ n J,i are the contributions of the gradient quantities in the cell to the gradient quantities at point J. As an example, Fig. 1 shows four rectangular cells surrounding point J. In order to compute the total gradient quantities at point J, Eq.(13) is applied to all cells surrounding the point J such that φ n J = I φ n J,i, (14) i=1 Fig. 1 Rectangular cells surrounding point J
4 484 S. Phongthanapanich, P. Dechaumphai where I is the number of the surrounding cells. Finally, the gradient quantities at Ɣ ij, φij n, are then computed by applying the trapezoidal rule along the common edge of cells i and j. Finally, for simplicity, the order of accuracy and stability of the explicit numerical scheme [24] given by Eq. (11) will be analyzed on uniform one-dimensional grid cell, = x. The numerical equation for ith cell, (x i 1/2, x i+1/2 ), may be written as following φ n+1 i = φ n i + ɛ t x ( φ x n i+1/2 φ x n i 1/2 t ( κφ n i f i ). (15) By using one-dimensional linear interpolation function, the gradient quantities at cell faces i 1/2 and i + 1/2 can be written by φ x φ x n i 1/2 n i+1/2 = 1 x (φn i φ n i 1 ), = 1 x (φn i+1 φn i ), ) (16) and by using these expressions on the right hand side of Eq. (15) togive φ n+1 i = φ n i + α(φ n i+1 2φn i + φ n i 1 ) t ( κφ n i f i ), (17) where α = ɛ t/( x) 2. Then, truncation error analysis by applying the Taylors series expansion on Eq. (17) evaluated at (x n, t n ) = (i, n) to provide that the order of accuracy of Eq. (17) with the second-order accurate Runge Kutta time stepping method is O( t 2, x 2 ). To analyze the stability of the numerical scheme, the discrete Fourier transform is applied to the homogeneous part of Eq. (17), term by term, to obtain φ n+1 =[1 2α(1 cos(ξ)) κ t] φ n [ ( ) ] ξ = 1 4α sin 2 κ t φ n. (18) 2 By choosing α<1/2, the amplification factor is given by ( ) ξ 1 4α sin2 κ t 2 eκ t. (19) Fig. 2 Numerical solutions on grid size at four different times of problem 3.1: a t = 0.25, b t = 0.50, c t = 0.75, d t = 1.0
5 Finite volume element method for analysis of unsteady reaction diffusion problems 485 Lastly, by applying the result of Eq. (18) n + 1 times and using the Parseval s Identity, to yield φ n+1 2, x e κ(n+1) t φ 0 2, x. (20) Hence by the definition of stability [24], Eq. (17) is stable with respect to the discrete L 2 norm. 3 Numerical examples To evaluate the robustness and accuracy of the proposed finite volume element method, three reaction diffusion problems subjected to homogeneous boundary conditions are examined. The first example is used to evaluate the capability of the proposed method for the singularly perturbed problem where the diffusion coefficient is very small as compared to the reaction coefficient. In the second example, the diffusion coefficient is given as one in order to examine the robustness of the method when a small time step is employed. The numerical results obtained from these two examples will be compared with the finite element solutions presented by Franca et al. [12]. The last example is similar to the first one but with non-zero homogeneous boundary conditions. The accuracy of the numerical solutions obtained from this example will be compared with the finite element solutions presented by Franca et al. [13]. 3.1 Reaction-dominated diffusion problem The first example is a reaction-dominated diffusion problem, with ɛ = 10 6, κ = 1, and f = 1. The initial condition is given by φ 0 (x) = 0. A unit square domain is divided into 400 uniform square grids (20 20). The problem is analyzed using the fixed time step of t = 0.1. Figure 2a d show the numerical solutions of the tested grids at four different times. The four computed profiles show that there is no spurious oscillation along the high-gradient boundary layers. For the purpose of comparison with the solutions presented in Ref. [12] (variant of the PGEM), the analysis is repeated but at different times of t = 0.1, and t = 7.5. At time t = 0.1, the computed profiles and the solutions at x = 0.5 are shown in Fig. 3a and b, respectively. Figure 4a and b show the similar behavior of the solution at the time of t = 7.5. Spurious oscillation does not occur along the high-gradient boundary layers of the solution profiles at these two times. The maximum values are comparable with the solutions from the method presented by Ref. [12]. Fig. 3 Numerical solution of problem 3.1: a At time t = 0.1, b along the line x = 0.5 Fig. 4 Numerical solution of problem 3.1: a At time t = 7.5, b along the line x = 0.5
6 486 S. Phongthanapanich, P. Dechaumphai Fig. 5 Numerical solutions on grid size at four different times of problem 3.2: a t = 0.01, b t = 0.05, c t = 0.1, d t = 1.0 Fig. 6 Numerical solution of problem 3.2: a At time t = 10 5, b along the line x = Reaction diffusion problem The second example is a reaction diffusion problem, with ɛ = 1, κ = 1, and f = 1. The initial condition is given by φ 0 (x) = 0, and the unit square domain is divided into 400 uniform square grids (20 20) in the same fashion as shown in the first example. The objective of this second example is to examine the robustness of the method using a small time step value of t = Figure 5a d shows the computed solutions of the tested grids at four different times of t = 0.01, 0.05, 0.1, and 1.0, respectively. The analysis of this problem is performed again at the two different times of t = 10 5, and t = 10 3, so that the solutions can be compared with the finite element solutions from Ref. [12]. At the time t = 10 5, the computed profiles and the solutions at x = 0.5 are shown in Fig. 6a and b, respectively. There is no spurious oscillation of the computed solution profiles along the high-gradient boundary layers at this earlier time. Figure 7a and b shows the same solution behavior but at the time t = 10 3.Again, the computed maximum values are comparable with the solutions from the method proposed by Ref. [12]. The present method provided less diffusive than the Variant of the PGEM.
7 Finite volume element method for analysis of unsteady reaction diffusion problems 487 Fig. 7 Numerical solution of problem 3.2: a At time t = 10 3, b along the line x = 0.5 Fig. 8 Numerical solutions on grid size at four different times of problem 3.3: a t = 0.25, b t = 0.50, c t = 0.75, d t = 1.0
8 488 S. Phongthanapanich, P. Dechaumphai Fig. 9 Numerical solution of problem 3.3: a At time t = 7.5, b along the line x = Reaction-dominated diffusion problem with non-zero boundary conditions The parameters for the last example are given by ɛ = 10 6, κ = 1, and f = 1. The initial condition is given by φ 0 (x) = 0. The boundary conditions along x = 0 and y = 1 are set as one, while the remaining boundaries are specified as zero. The computational domain and time step used in the first example are also used herein. Figure 8a d shows the numerical solutions of the tested grid at four different times of t = 0.25, 0.50, 0.75, and 1.0, respectively. The computed profiles show that there is no spurious oscillation along highgradient boundary layers. Again, for the purpose of comparison, the steady-state solution profile at time t = 7.5 is shown in Fig. 9a, and the corresponding solution at x = 0.5 are presented in Fig. 9b. These figures show that spurious oscillation does not occur along the high-gradient boundary layers at these final times. The computed maximum values are also comparable with those obtained from the method of Ref. [13]. 4 Conclusions This paper presents an explicit finite volume element method for analyzing the reaction diffusion equation on twodimensional domain. The theoretical formulation of the combined finite volume and finite element method was explained in details. The finite volume method was applied for discretizing the computational domain, while the finite element method was used to estimate the gradient quantities at the cell faces. The numerical test cases have shown that the combined method does not require any special numerical treatment to improve the solution stability. Three examples were used to evaluate the robustness and efficiency of the combined method. These examples show that the combined method can provide high-resolution profiles without spurious oscillation along the high-gradient boundary layers. Acknowledgments The authors are pleased to acknowledge the Thailand Research Fund (TRF) for supporting this research work. References 1. Apel, T., Lube, G.: Anisotropic mesh refinement for a singularly perturbed reaction diffusion model problem. Appl. Numer. Math. 26, (1998) 2. Li, J., Navon, I.M.: Uniformly convergent finite element meth-ods for singularly perturbed elliptic boundary value problems. I. Reaction diffusion type. Comput. Math. Appl. 35, (1998) 3. Castaings, W., Navon, I.M.: Mesh refinement strategies for solving singularly perturbed reaction diffusion problems. Comput. Math. Appl. 41, (2001) 4. Lin, B.T., Madden, N.: A finite element analysis of a coupled system of singularly perturbed reaction diffusion equations. Appl. Math. Comput. 148, (2004) 5. El-Salam, M.R.A., Shehata, M.H.: The numerical solution for reaction diffusion combustion with fuel consumption. Appl. Math. Comput. 160, (2005) 6. Xu, R.: A reaction diffusion predator prey model with stage structure and nonlocal delay. Appl. Math. Comput. 175, (2006) 7. Heineken, W., Warnecke, G.: Partitioning methods for reaction diffusion problems. Appl. Numer. Math. 56, (2006) 8. Parvazinia, M., Nassehi, V.: Multiscale finite element modeling of diffusion reaction equation using bubble functions with bilinear and triangular elements. Comput. Methods Appl. Mech. Eng. 196, (2007) 9. Moore, P.K.: Solving regularly and singularly perturbed reaction diffusion equations in three space dimensions. J. Comput. Phys. 224, (2007) 10. Roos, H.G., Stynes, M., Tobiska, L.: Numerical Methods for Singularly Perturbed Differential Equations. Springer-Verlag, Berlin (1996) 11. Xenophontos, C.: Optimal mesh design for the finite element approximation of reaction diffusion problems. Int. J. Numer. Meth. Eng. 53, (2002)
9 Finite volume element method for analysis of unsteady reaction diffusion problems Franca, L.P., Ramalho, J.V.A., Valentin, F.: Enriched finite element methods for unsteady reaction diffusion problems. Commun. Numer. Meth. Eng. 22, (2006) 13. Franca, L.P., Valentin, F.: On an improved unusual stabilized finite element method for the advective reactive diffusive equation. Comput. Methods Appl. Mech. Eng. 190, (2000) 14. Franca, L.P., Madureira, A.L., Valentin, F.: Towards multiscale functions: enriching finite element spaces with local but not bubble-like functions. Comput. Methods Appl. Mech. Eng. 194, (2005) 15. Schwab, C., Suri, M., Xenophontos, C.: The hp finite element method for problems in mechanics with boundary layers. Comput. Methods Appl. Mech. Eng. 157, (1996) 16. Schwab, C., Melenk, J.M.: hp FEM for reaction diffusion equations. I. Robust exponential convergence. SIAM J. Numer. Anal. 35, (1998) 17. Liao, W., Zhu, J., Khaliq, A.Q.M.: An effficient high-order algorithm for solving systems of reaction diffusion equations. Numer. Methods Partial Differ. Equ. 18, (2002) 18. Rui, H.: Symmetric mixed covolume methods for parabolic problems. Numer. Methods Partial Differ. Equ. 18, (2002) 19. Ma, X., Shu, S., Zhou, A.: Symmetric finite volume discretizations for parabolic problems. Comput. Methods Appl. Mech. Eng. 192, (2003) 20. Wang, T.: Alternating direction finite volume element methods for 2D parabolic partial differential equations. Numer. Methods Partial Differ. Equ. 24, (2008) 21. Wang, H., Liu, J.: Development of CFL-free, explicit schemes for multidimensional advection reaction equations. SIAM J. Sci. Comput. 23, (2001) 22. Shu, C.W., Osher, S.: Efficient implementation of essentially non-oscillatory shock-capturing schemes. J. Comput. Phys. 77, (1988) 23. Zienkiewicz, O.C., Taylor, R.L.: The Finite Element Method for Solid and Structural Mechanics 6th edn. Butterworth-Heinemann, Burlington (2005) 24. Richtmyer, R.D., Morton, K.W.: Difference Methods for Initialvalue Problems. Wiley-Interscience, New York (1967)
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
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 informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
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 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 informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 219 On finite element methods for 3D time dependent convection diffusion reaction equations
More informationOn finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion
On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion Volker John and Ellen Schmeyer FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationSUPERCONVERGENCE OF DG METHOD FOR ONE-DIMENSIONAL SINGULARLY PERTURBED PROBLEMS *1)
Journal of Computational Mathematics, Vol.5, No., 007, 185 00. SUPERCONVERGENCE OF DG METHOD FOR ONE-DIMENSIONAL SINGULARLY PERTURBED PROBLEMS *1) Ziqing Xie (College of Mathematics and Computer Science,
More informationSPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM. Mirjana Stojanović University of Novi Sad, Yugoslavia
GLASNIK MATEMATIČKI Vol. 37(57(2002, 393 403 SPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM Mirjana Stojanović University of Novi Sad, Yugoslavia Abstract. We introduce piecewise interpolating
More informationNumerical Methods for Partial Differential Equations CAAM 452. Spring 2005
Numerical Methods for Partial Differential Equations Instructor: Tim Warburton Class Location: Duncan Hall 1046 Class Time: 9:5am to 10:40am Office Hours: 10:45am to noon in DH 301 CAAM 45 Spring 005 Homeworks
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 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 informationOn the Comparison of the Finite Volume and Discontinuous Galerkin Methods
Diploma Thesis Institute for Numerical Simulation, TUHH On the Comparison of the Finite Volume and Discontinuous Galerkin Methods Corrected version Katja Baumbach August 17, 2006 Supervisor: Prof. Dr.
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
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 informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
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 informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationAdaptive Time Space Discretization for Combustion Problems
Konrad-Zuse-Zentrum fu r Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany JENS LANG 1,BODO ERDMANN,RAINER ROITZSCH Adaptive Time Space Discretization for Combustion Problems 1 Talk
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 Studies of Adaptive Finite Element Methods for Two Dimensional Convection-Dominated Problems
DOI 10.1007/s10915-009-9337-6 Numerical Studies of Adaptive Finite Element Methods for Two Dimensional Convection-Dominated Problems Pengtao Sun Long Chen Jinchao Xu Received: 17 April 2009 / Revised:
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 Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems
Finite Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems Zhimin Zhang Department of Mathematics Wayne State University Detroit, MI 48202 Received 19 July 2000; accepted
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 informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationA PARAMETER ROBUST NUMERICAL METHOD FOR A TWO DIMENSIONAL REACTION-DIFFUSION PROBLEM
A PARAMETER ROBUST NUMERICAL METHOD FOR A TWO DIMENSIONAL REACTION-DIFFUSION PROBLEM C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN Abstract. In this paper a singularly perturbed reaction diffusion partial
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 informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
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 high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
More informationNumerical Solutions of Laplacian Problems over L-Shaped Domains and Calculations of the Generalized Stress Intensity Factors
WCCM V Fifth World Congress on Computational Mechanics July 7-2, 2002, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Numerical Solutions of Laplacian Problems over L-Shaped Domains
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 informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationWeak Galerkin Finite Element Methods and Applications
Weak Galerkin Finite Element Methods and Applications Lin Mu mul1@ornl.gov Computational and Applied Mathematics Computationa Science and Mathematics Division Oak Ridge National Laboratory Georgia Institute
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 informationRECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS
RECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS Dean Wang April 30, 2015 24.505 Nuclear Reactor Physics Outline 2 Introduction and Background Coupled T-H/Neutronics Safety Analysis Numerical schemes
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 informationResearch Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation with Inconsistent Initial and Boundary Conditions
International Mathematics and Mathematical Sciences Volume 212, Article ID 82197, 13 pages doi:1.1155/212/82197 Research Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation
More informationHigh-Order Finite-Volume Methods! Phillip Colella! Computational Research Division! Lawrence Berkeley National Laboratory!
High-Order Finite-Volume Methods! Phillip Colella! Computational Research Division! Lawrence Berkeley National Laboratory! Why Higher Order?! Locally-refined grids, mapped-multiblock grids smooth except
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 informationA Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional Hyperbolic Telegraph Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(01) No.3,pp.59-66 A Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional
More informationChapter 5. Formulation of FEM for Unsteady Problems
Chapter 5 Formulation of FEM for Unsteady Problems Two alternatives for formulating time dependent problems are called coupled space-time formulation and semi-discrete formulation. The first one treats
More informationA Petrov-Galerkin Enriched Method: A Mass Conservative Finite Element Method For The Darcy Equation
University of Colorado at Denver and Health Sciences Center A Petrov-Galerkin Enriched Method: A Mass Conservative Finite Element Method For The Darcy Equation Gabriel R. Barrenechea, Leopoldo P. Franca,
More informationPartial Differential Equations
Next: Using Matlab Up: Numerical Analysis for Chemical Previous: Ordinary Differential Equations Subsections Finite Difference: Elliptic Equations The Laplace Equations Solution Techniques Boundary Conditions
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationq t = F q x. (1) is a flux of q due to diffusion. Although very complex parameterizations for F q
! Revised Tuesday, December 8, 015! 1 Chapter 7: Diffusion Copyright 015, David A. Randall 7.1! Introduction Diffusion is a macroscopic statistical description of microscopic advection. Here microscopic
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationComputation Fluid Dynamics
Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand
More informationQuintic B-Spline Galerkin Method for Numerical Solutions of the Burgers Equation
5 1 July 24, Antalya, Turkey Dynamical Systems and Applications, Proceedings, pp. 295 39 Quintic B-Spline Galerkin Method for Numerical Solutions of the Burgers Equation İdris Dağ 1,BülentSaka 2 and Ahmet
More informationSuperconvergence analysis of the SDFEM for elliptic problems with characteristic layers
Superconvergence analysis of the SDFEM for elliptic problems with characteristic layers S. Franz, T. Linß, H.-G. Roos Institut für Numerische Mathematik, Technische Universität Dresden, D-01062, Germany
More informationSelf-similar solutions for the diffraction of weak shocks
Self-similar solutions for the diffraction of weak shocks Allen M. Tesdall John K. Hunter Abstract. We numerically solve a problem for the unsteady transonic small disturbance equations that describes
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationFourier Type Error Analysis of the Direct Discontinuous Galerkin Method and Its Variations for Diffusion Equations
J Sci Comput (0) 5:68 655 DOI 0.007/s095-0-9564-5 Fourier Type Error Analysis of the Direct Discontinuous Galerkin Method and Its Variations for Diffusion Equations Mengping Zhang Jue Yan Received: 8 September
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 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 informationExtremum-Preserving Limiters for MUSCL and PPM
arxiv:0903.400v [physics.comp-ph] 7 Mar 009 Extremum-Preserving Limiters for MUSCL and PPM Michael Sekora Program in Applied and Computational Mathematics, Princeton University Princeton, NJ 08540, USA
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 informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationA numerical method for the wave equation subject to a non-local conservation condition
A numerical method for the wave equation subject to a non-local conservation condition Whye-Teong Ang Division of Engineering Mechanics School of Mechanical and Production Engineering Nanyang Technological
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationBeyond Wiener Askey Expansions: Handling Arbitrary PDFs
Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2005) DOI: 10.1007/s10915-005-9038-8 Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Xiaoliang Wan 1 and George Em Karniadakis 1
More informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationRening the submesh strategy in the two-level nite element method: application to the advection diusion equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 39:161 187 (DOI: 10.1002/d.219) Rening the submesh strategy in the two-level nite element method: application to
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationAdaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear Degenerate Reaction-Diffusion Problems
EPJ Web of Conferences 108, 0019 (016) DOI: 10.1051/ epjconf/ 0161080019 C Owned by the authors, published by EDP Sciences, 016 Adaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear
More informationSolving the One Dimensional Advection Diffusion Equation Using Mixed Discrete Least Squares Meshless Method
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 2015 Paper No. 292 Solving the One Dimensional Advection Diffusion Equation
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
More informationA parameter-free smoothness indicator for high-resolution finite element schemes
Cent. Eur. J. Math. 11(8) 2013 1478-1488 DOI: 10.2478/s11533-013-0254-4 Central European Journal of Mathematics A parameter-free smoothness indicator for high-resolution finite element schemes Research
More informationA PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM
A PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM Hyun Geun LEE 1, Jeong-Whan CHOI 1 and Junseok KIM 1 1) Department of Mathematics, Korea University, Seoul
More informationRobust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations
Robust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
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 informationNUMERICAL SOLUTION OF GENERAL SINGULAR PERTURBATION BOUNDARY VALUE PROBLEMS BASED ON ADAPTIVE CUBIC SPLINE
TWMS Jour. Pure Appl. Math., V.3, N.1, 2012, pp.11-21 NUMERICAL SOLUTION OF GENERAL SINGULAR PERTURBATION BOUNDARY VALUE PROBLEMS BASED ON ADAPTIVE CUBIC SPLINE R. MOHAMMADI 1 Abstract. We use adaptive
More informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
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 informationALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS
ALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS HAILIANG LIU AND MICHAEL POLLACK Abstract. In this work, we propose a high resolution Alternating Evolution Discontinuous
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationNumerical Methods for the Landau-Lifshitz-Gilbert Equation
Numerical Methods for the Landau-Lifshitz-Gilbert Equation L ubomír Baňas Department of Mathematical Analysis, Ghent University, 9000 Gent, Belgium lubo@cage.ugent.be http://cage.ugent.be/~lubo Abstract.
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationChapter 1. Introduction and Background. 1.1 Introduction
Chapter 1 Introduction and Background 1.1 Introduction Over the past several years the numerical approximation of partial differential equations (PDEs) has made important progress because of the rapid
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
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 informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
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 informationA Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws
A Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws Mehdi Dehghan, Rooholah Jazlanian Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University
More informationAccuracy-Preserving Source Term Quadrature for Third-Order Edge-Based Discretization
Preprint accepted in Journal of Computational Physics. https://doi.org/0.06/j.jcp.07.04.075 Accuracy-Preserving Source Term Quadrature for Third-Order Edge-Based Discretization Hiroaki Nishikawa and Yi
More informationIterated Defect Correction with B-Splines for a Class of Strongly Nonlinear Two-Point Boundary Value Problems
American Review of Mathematics and Statistics June 2016, Vol. 4, No. 1, pp. 31-44 ISSN: 2374-2348 (Print), 2374-2356 (Online) Copyright The Author(s). All Rights Reserved. Published by American Research
More informationNumerical Solution of partial differential equations
G. D. SMITH Brunei University Numerical Solution of partial differential equations FINITE DIFFERENCE METHODS THIRD EDITION CLARENDON PRESS OXFORD Contents NOTATION 1. INTRODUCTION AND FINITE-DIFFERENCE
More informationRESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS
Proceedings of ALGORITMY 2016 pp. 113 124 RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS VÍT DOLEJŠÍ AND FILIP ROSKOVEC Abstract.
More information