NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE
|
|
- Kevin Hopkins
- 5 years ago
- Views:
Transcription
1 Global and Stocastic Analysis Vol. 4 No. 1, January 2017, 1-10 NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE K. PHANEENDRA AND E. SIVA PRASAD Abstract. Tis paper presents a class of accurate non standard finite difference metod based on cubic spline function on uniform mes for te numerical solution of a second order singularly perturbed two point boundary value problems associate wit boundary layer teory. A fitting factor known to be artificial viscosity was introduced in discretization equation of non standard finite difference sceme and its value was obtained from te teory of singular perturbations. We ave solved te discretization equation using discrete invariant imbedding. Te proposed sceme is analyzed for Convergence. Te numerical results for several test examples demonstrate te efficiency of te proposed metod. 1. Introduction Singularly perturbed two point boundary value problems occurs frequently in many areas of applied matematics and matematical pysics suc as fluid mecanics, elasticity, optimal control, cemical-reactor teory, aerodynamics, reaction diffusion process, geopysics and many oter related areas. Solutions of tis type of equations exibit boundary layer penomena; tat is, te solution of tese problems varies rapidly in some parts and varies slowly in some oter parts. Te numerical treatment of singularly perturbed differential equations is far from trivial and also gives major computational difficulties due to te presence of boundary and/or interior layers. Standard numerical metods fail to approximate teir solutions because of te boundary layer beavior. So to capture te layers, two classical approaces are used. Te first approac is to design a special mes known as Siskin mes [15] or graded meses [3] and te second approac is a fitting factor is introduced in te existing finite difference metods known as fitted operator metod [9] to control te beavior in te boundary layer. Tis paper uses te second approac to solve te singularly perturbed two point boundary value problems. To describe tis metod, we consider singularly perturbed two point boundary value problems of te form y x + pxy x + qxyx = rx, 0 x Matematics Subject Classification. Primary 65L11; Secondary 65L12. Key words and prases. Singularly perturbed boundary value problem, Fitting factor, Cubic spline, Tridiagonal system, Invariant imbedding algoritm. 1 1
2 2 K. PHANEENDRA AND E. SIVA PRASAD wit boundary conditions y0 = α, y1 = β 1.2 were 0 < << 1, px, qx and rx are bounded continuous functions in 0, 1, and α, β are finite constants. Tis type of problem was solved asymptotically by Bellman [1], Bender and Orszag [2], Kevorkian and Cole [6], Nayfe [10], O Malley [11] and numerically by Kreiss [7], Miller [9], Kadalbajoo and Kumar [5], Rasidinia etc. [14], Lin and Vancouver [8], Reddy [4, 13]. In tis paper, a fitted non standard finite difference sceme using cubic spline functions on a uniform mes was presented for solving linearly singularly perturbed two-point boundary value problems wit boundary layer at one end point. In section 2, cubic spline was discussed briefly. In section 3, numerical metod using cubic spline was presented. In section 4, convergence analysis of te metod was discussed. To justify te proposed metod, test examples are illustrated in section 5. Finally, in section 6, summary and conclusions of te metod are given. 2. Description of te metod Discretize te interval [0,1] into N equal subintervals of mes size = 1 N, so tat x i = i, i =0, 1, 2,..., N wit 0 = x 0, 1 = x N. Let yx be te exact solution and y i be an approximation to yx i obtained by te polynomial cubic splines i x passing troug te points x i, y i andx i+1, y i+1. From te teory of te Splines, S i x satisfies not only te interpolatory conditions at x i andx i+1 but also te continuity of first derivative at te common nodes x i, y i are fulfilled. For eac i t segment, te cubic polynomial spline function S i xas te form S i x = a i + b i x x i + c i x x i 2 + d i x x i 3, i = 0, 1, 2,..., N were a i, b i, c i and d i are constants. A cubic spline function Sx of class C 2 [a, b] interpolates yx at te grid pointsx i for i = 0,1,2...,N. To develop expressions for te four coefficients of Eq. 2.1 in terms ofy i,y i+1,m i and M i+1. We first define S i x i = y i,s i x i+1 = y i+1,s i x i = M i,s i x i+1 = M i+1. From algebraic manipulation, we obtain te following expression: a i = y i, b i = y i+1 y i M i+1 M i, c i = M i 6 2, d i = M i+1 M i 6 Were i = 0, 1, 2...N-1. Using te continuity of cubic spline Sx and its first derivative atx i, y i, tat iss i 1 x i = S i x i, we obtain te following relations for i =1,2,...,N-1. M i+1 + 4M i + M i 1 = 6 2 y i+1 2y i + y i Numerical Sceme 3.1. Left end boundary layer. We assume tat qx 0, px M > 0 trougout te interval [0, 1], were M is positive constant. Under tese assumptions, Eq. 1.1 as unique solution and display boundary layer at x = 0 for small values of. 2
3 NON STANDARD FDM FOR SPP USING CUBIC SPLINE 3 From te teory of singular perturbation it is known tat te solution of Eq. 1.1 and Eq. 1.2 is of te form O Malley [11] yx = y 0 x + p0 px α y 00 e x 0 px pxdx qx + O 3.1 were y 0 x is te solution of reduced problem pxy 0x + qxy 0 x = rx,wit y 0 1 = β. Expanding px and qx about te point 0 up to first term using Taylor s series, Eq. 3.1 gives From Eq. 3.2 we ave i.e., terefore yx = y 0 x + α y 0 0 e yx i = y 0 x i + α y 0 0 e yi = y 0 i + α y 0 0 e lim yi = y 00 + α y 0 0 e 0 p0 p0 p0 q0 p0x + O 3.2 q0 p0x i + O, q0 p0i + O, p 2 0 q0 iρ p0 + O, wereρ = 3.3 At te grid pointx i, te proposed differential equation Eq. 1.1 may be discretized by M i = rx i px i y ix qx i yx i Te non standard finite differences of y i 1, y i+1 and y i are y i 1 y i+1 + 4y i 3y i 1 y i+1 3y i+1 4y i + y i 1 y i y i+1 y i 1. Substituting te values of M i, M i 1 and M i+1 along wit above differences in Eq. 2.2, we get 6 y 2 i+1 2y i + y i 1 = pi+1 + 2pi q i 1 + 3pi 1 y i 1 + 2pi+1 4q i 2pi 1 y i + 3pi+1 2pi q i 1 + pi 1 + r i+1 + 4r i + r i 1 Introducing fitting factor σρ in Eq. 3.4, we get 6σρ 2 y i+1 2y i + y i 1 = + 3pi+1 2pi q i 1 + pi 1 pi+1 + 2pi q i 1 + 3pi 1 + r i+1 + 4r i + r i 1 y i 1 + 2pi q i 2pi 1 y i 3.5 3
4 4 K. PHANEENDRA AND E. SIVA PRASAD Multiplying Eq. 3.5 by and taking limit as 0, we get σ p0 lim yi + 1 2yi + yi 1 = lim yi 1 yi ρ By substituting Eq. 3.3 in to Eq. 3.6, we get σ = ρ 4 p0 cot p0 2 q0 p0 ρ is te required fitting factor in te left end boundary layer. From Eq.3.5 we get te following tridiagonal system were E i y i 1 + F i y i + G i y i+1 = H i for i = 1, 2,..., N E i = σ+ 4 p i p i 12 p i q i 1, F i = 2σ 3 p i p i q i G i = σ + 12 p i 1 3 p i 4 p i q i+1, H i = [r i r i + r i+1 ] px i = p i,qx i = q i, rx i = r i, for i = 0, 1,...,N. We solve te tridiagonal system 3.8, using metod of invariant imbedding algoritm Rigt-end boundary layer. Furter, if we assume tat qx 0, px M < 0 trougout te interval [0, 1] were M is negative constant. Under tese assumptions, 1.1 as unique solution and display boundary layer at x =1 for small values of. From te teory of singular perturbation it is known tat te solution of Eq. 1.1 and Eq. 1.2 is of te form yx = y 0 x + p1 px α y 01 e x 0 px pxdx qx + O 3.9 were y 0 x is te solution of pxy 0x + qxy 0 x = rx,wit y 0 0 = α. Expanding px and qx about te point 1 up to first term using Taylor s series, Eq. 3.9 gives yx = y 0 x + α y 0 1 e yx i = y 0 x i + α y 0 1 e i.e., yi = y 0 i + α y 0 1 e p1 p1 p1 q1 p11 x + O 3.10 q1 p11 x i + O, q1 p11 i + O, terefore lim 00 + α y 0 1 e 0 p 2 1 q1 p1 1 iρ + O 3.11 were ρ = Now, we consider cubic spline finite difference sceme and introduce te fitting 4
5 NON STANDARD FDM FOR SPP USING CUBIC SPLINE 5 factor σρ as: 6σρ 2 y i+1 2y i + y i 1 = + 3pi+1 2pi q i 1 + pi 1 pi+1 + 2pi q i 1 + 3pi 1 y i 1 + 2pi+1 + r i+1 + 4r i + r i 1 Multiplying 3.12 by and taking limit as 0, we get 4q i 2pi 1 y i 3.12 σ p0 lim yi + 1 2yi + yi 1 = lim yi 1 yi ρ By substituting Eq in to Eq we get σ = ρ 4 p0 cot p1 2 q1 p1 ρ 2 is te fitting factor in te rigt end boundary layer. From Eq.3.12 we get te following tridiagonal system were 3.14 E i y i 1 + F i y i + G i y i+1 = H i for i = 1, 2,..., N E i = σ+ 4 p i p i 12 p i q i 1, F i = 2σ 3 p i p i q i G i = σ + 12 p i 1 3 p i 4 p i q i+1, H i = [r i r i + r i+1 ] px i = p i,qx i = q i, rx i = r i, for i = 0, 1,...,N To solve te above system 3.15, we used metod of invariant imbedding algoritm. 4. Convergence Analysis Incorporating te given boundary conditions we obtain te matrix equation as D + P Y + Q + T = σ σ σ 2σ σ were D = [ σ, 2σ, σ] = 0 σ 2σ σ σ 2σ and v 1 w z 2 v 2 w P = [z i, v i, w i ] = 0 z 3 v 3 w z N 1 v N 1 5
6 6 K. PHANEENDRA AND E. SIVA PRASAD were z i = 4 p i p i 12 p i q i 1, v i = 3 p i p i q w i = 12 p i 1 3 p i 4 p i q i+1fori = 1, 2, 3, 4...N 1 and Q = [q 1 + σ + z 1 α, q 2, q 3,..., q N 2, q N 1 + σ + w N 1 β] T were q i = 2 6 r i 1 + 4r i + r i+1 i = 1, 2, 3, 4,..., N 1 T = 0 4 and Y = [Y 1, Y 2, Y 3,..., Y N 1 ] T, T = [T 1, T 2,..., T N 1 ] T, O = [0, 0,..., 0] T are associated vectors of equation 4.1. Let y = [y 1, y 2,..., y N 1 ] T = Y wic satisfies te equation D + P y + Q = Lete i = y i Y i, i = 11N 1be te discretization error so tat E = [e 1, e 2,..., e N 1 ] T = y Y. Subtracting 4.1 from 4.2, we obtain te error equation D + P E = T 4.3 Let px C 1 and qx C 2 were C 1, C 2 are positive constants. If P i,j be te i, j t element of te matrix D + P, ten P i,i+1 σ + 3 C C 2, i = 1, 2,..., N 2 P i,i 1 σ + 3 C C 2i = 2, 3,..., N 1 Tus for sufficiently small i i.e. as 0 we ave P i,i+1 < σ, i = 1, 2,...,.N 2 P i,i 1 < σ, i = 2, 3,..., N Hence D + P is irreducible see Ref. [16]. Let S i be te sum of te elements of te it row of te matrix D + P, ten we ave S i = σ 4 p i 1 3 p i + 12 p i q i + q i+1 for i = 1 S i = 2 6 qi-1 + 4q i + q i+1 for i = 2, 3,..., N 2 S i = σ 12 p i p i + 4 p i q i + q i 1 for i = N 1 Let C 1 = min px and 1 i N C 1 = max px, C 2 1 i N = min qx and 1 i N C 2 = max qx. 1 i N Since 0 < << 1 and O,it is verify tat for sufficiently small, D + P is monotone [16, 17]. Hence D + P 1 exists and D + P 1 0. Tus from Eq.4.3, we ave E D + P 1 T 4.5 6
7 NON STANDARD FDM FOR SPP USING CUBIC SPLINE 7 For sufficiently small, we ave Let D + P 1 i,k D + P 1 = Since D + P 1 i,k Hence Furtermore, N 1 k=1 S i > 2 C 2 for i = 1 S i > 2 C 2 for i = n 1 and S i > 2 C 2 for i = 2, 3,..., n 2. be te i, kt element of D + P 1 and we define N 1 max 1 i N 1 k=1 N 1 0 and k=1 D + P 1 i,k and T = max T, i N 1 D + P 1 i,k.s k = 1 for i = 1, 2, 3,..., N 1. D + P 1 i,k 1 S i < 1 2 C 2, i = D + P 1 i,k 1 S i < 1 2 C 2, i = N D + P 1 i,k 1 min 2 i N 2 S i By te elp of Eqs , from 4.5, we obtain Hence te metod is second order convergence. < 1 2 C 2 for i = N E O Numerical examples Tis section presents test examples to demonstrate te efficiency of te metod computationally. We consider two numerical examples wit left-end boundary layer and one problem wit rigt-end boundary layer of te underlying interval were considered. Tese problems were cosen because tey ave been widely discussed in te literature and exact solutions were available for comparison. Te maximum absolute errors wit and witout fitting factor were presented to support te given metod. Example 1. [2], Consider te boundary value problem from Bender and Orszag y x + y x yx = 0; x [0, 1] wit y0=1 and y1=1. Clearly tis problem as a boundary layer at x = 1. Te exact solution is given by yx = [em2 1 e m1x + 1 e m1 e m2x ] e m2 e m1 were m 1 = /2 and m 2 = /2. Table 1 gives te maximum absolute errors for different values of and wit and witout fitting factor. 7
8 8 K. PHANEENDRA AND E. SIVA PRASAD Example 2. Consider te following non-omogeneous problem from fluid dynamics for fluid of small viscosity y x + y x = 1 + 2x; x [0, 1] wit y0 = 0 and y1 = 1. Te exact solution is given by yx = xx e x/ 1 e 1/. Table 2 gives te maximum absolute errors for different values of and wit and witout fitting factor. Example 3. Consider te boundary problem y x y x 1 + yx = 0; x [0, 1] wit y0 = 1 + exp 1 + / and y1 = 1 + 1/e Clearly, tis problem as a boundary layer at x = 1. Te exact solution is given by yx = e 1+x 1/ + e x Table 3 gives te maximum absolute errors for different values of and wit and witout fitting factor. 6. Summary and conclusions Non standard fitted finite difference metod for singular perturbation problems using cubic splines was developed. We ave introduced a fitting factor wic is called artificial viscosity in te spline difference sceme to control te rapid canges in te boundary layer. Convergence analysis of te metod is discussed and it sows tat our metod is second order convergent. We ave presented maximum absolute errors for te standard test examples cosen from literature wit te proposed metod wit and witout fitting factor. From te numerical results, it sows te importance of te fitting factor introduced in te spline difference sceme. References 1. Bellman,R.: Perturbation Tecniques in Matematics, Pysics and Engineering, Holt, Rineart & Winston, New York, Bender, C.M., Orszag, S.A.: Advanced Matematical Metods for Scientists and Engineers, McGraw-Hill, New York, Gartland, Jr., E.C.: Graded-mes difference sceme for singular perturbed two point boundary value problems, Mat. Comp Kadalbajoo, M.K., Reddy Y.N.: A non asymptotic metod for general singular perturbation problems, Journal of Optimization Teory and Applications, Kadalbajoo, M.K., Devendra Kumar: Initial value tecnique for singularly perturbed two point boundary value problems using an exponentially fitted finite difference sceme, Computers and Matematics wit Applications, Kevorkian, J., Cole, J.D.: Perturbation Metods in Applied Matematics, Springer-Verlag, New York, Kreiss, B. Kreiss, H.O.: Numerical metods for singular perturbation problems, SIAM Journal on Numerical Analysis
9 NON STANDARD FDM FOR SPP USING CUBIC SPLINE 9 Table 1. Maximum absolute errors in te solution of Example 1 δ\ Wit fitting factor Witout fitting factor Table 2. Maximum absolute errors in te solution of Example 2 δ\ Wit fitting factor Witout fitting factor Table 3. Maximum absolute errors in te solution of Example 3 δ\ Wit fitting factor Witout fitting factor
10 10 K. PHANEENDRA AND E. SIVA PRASAD 8. Lin, P., Vancouver: A Numerical solution of quasilinear singularly perturbed ordinary differential equation witout turning points, Applied Matematics and Mecanics, Miller,J.J.H., O Riordan, E. and Siskin, G.I.: Fitted Numerical Metods for Singular Perturbation Problems, World Scientific, River Edge, NJ Nayfe A.H.: Introduction to Perturbation Tecniques, Wiley, New York, O Malley, R.E.: Introduction to Singular Perturbations, Academic Press, New York, Reinardt, H.J.: Singular Perturbations of difference metods for linear ordinary differential equations, Applicable Anal., Reddy Y.N., Awoke A.: An Exponentially Fitted Special second order finite difference metod for singular perturbation problems, Applied Matematics and Computation, Rasidinia, J., Moammadi, R. and Gasemi, M.: Cubic spline solution of singularly perturbed boundary value problems wit significant first derivatives, Applied Matematics and Computation Siskins, G.I.: A difference sceme for a singular perturbed equation of parabolic type wit a discontinuous initial condition, Sov. Mat. Dokl Varga, R.S.: Matrix Iterative Analysis, Prentice-Hall, Englewood Cliffs, New Jersey, Young, D.M.: Iterative Solutions of Large Linear Systems, Academic press, New York, K. Paneendra: Department of Matematics, University College of Science, Saifabad, Osmania University, Hyderabad, address: kollojupaneendra@yaoo.co.in URL: ttp:// jondoe optional E.Siva Prasad: Department of Matematics, Kavikulguru Institute of Tecnology & Science, Ramtek, Maarastra, address: emineni@yaoo.co.in 10
Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationAn Exponentially Fitted Non Symmetric Finite Difference Method for Singular Perturbation Problems
An Exponentially Fitted Non Symmetric Finite Difference Method for Singular Perturbation Problems GBSL SOUJANYA National Institute of Technology Department of Mathematics Warangal INDIA gbslsoujanya@gmail.com
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More information[21] B. J. McCartin, Theory of exponential splines, Journal of Approximation Theory 66 (1991) 1 23.
15 T. Ak, S. B. G. Karakoc, H. Triki, Numerical simulation for treatment of dispersive sallow water waves wit Rosenau- KdV equation, Te European Pysical Journal Plus 131 (1 (16 356-37. 16 M. Abbas, A.
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationParametric Spline Method for Solving Bratu s Problem
ISSN 749-3889 print, 749-3897 online International Journal of Nonlinear Science Vol4202 No,pp3-0 Parametric Spline Metod for Solving Bratu s Problem M Zarebnia, Z Sarvari 2,2 Department of Matematics,
More informationResearch Article Cubic Spline Iterative Method for Poisson s Equation in Cylindrical Polar Coordinates
International Scolarly Researc Network ISRN Matematical Pysics Volume 202, Article ID 2456, pages doi:0.5402/202/2456 Researc Article Cubic Spline Iterative Metod for Poisson s Equation in Cylindrical
More informationDIRECTLY SOLVING SECOND ORDER LINEAR BOUNDARY VALUE PROBLEMS OF ORDINARY DIFFERENTIAL EQUATIONS. Ra ft Abdelrahim 1, Z. Omar 2
International Journal of Pure and Applied Matematics Volume 6 No. 6, -9 ISSN: - (printed version); ISSN: -95 (on-line version) url: ttp://www.ijpam.eu doi:.7/ijpam.v6i. PAijpam.eu DIRECTLY SOLVING SECOND
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationNumerical Solution of Convection Diffusion Problem Using Non- Standard Finite Difference Method and Comparison With Standard Finite Difference Methods
IOSR Journal of Matematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 12, Issue 3 Ver. IV (May. - Jun. 2016), PP 94-109 www.iosrjournals.org Numerical Solution of Convection Diffusion Problem
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationTwo Step Hybrid Block Method with Two Generalized Off-step Points for Solving Second Ordinary Order Differential Equations Directly
Global Journal of Pure and Applied Matematics. ISSN 0973-768 Volume 2, Number 2 (206), pp. 59-535 Researc India Publications ttp://www.ripublication.com/gjpam.tm Two Step Hybrid Block Metod wit Two Generalized
More informationMore on generalized inverses of partitioned matrices with Banachiewicz-Schur forms
More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
More informationCONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL SINGULARLY PERTURBED DIFFERENTIAL EQUATIONS ON EQUIDISTRIBUTING MESHES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Number, Pages 57 74 c 5 Institute for Scientific Computing and Information CONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationChapter 1D - Rational Expressions
- Capter 1D Capter 1D - Rational Expressions Definition of a Rational Expression A rational expression is te quotient of two polynomials. (Recall: A function px is a polynomial in x of degree n, if tere
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
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 informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationDepartment of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801
RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused
More informationCalculus I - Spring 2014
NAME: Calculus I - Spring 04 Midterm Exam I, Marc 5, 04 In all non-multiple coice problems you are required to sow all your work and provide te necessary explanations everywere to get full credit. In all
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
More informationarxiv: v1 [math.dg] 4 Feb 2015
CENTROID OF TRIANGLES ASSOCIATED WITH A CURVE arxiv:1502.01205v1 [mat.dg] 4 Feb 2015 Dong-Soo Kim and Dong Seo Kim Abstract. Arcimedes sowed tat te area between a parabola and any cord AB on te parabola
More informationThe Sinc-Collocation Method for Solving the Telegraph Equation
Te Sinc-Collocation Metod for Solving te Telegrap Equation E. Hesameddini *1, E. Asadolaifard Department of Matematics, Faculty of Basic Sciences, Siraz University of Tecnology, Siraz, Iran *1 esameddini@sutec.ac.ir;
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationarxiv: v1 [math.pr] 28 Dec 2018
Approximating Sepp s constants for te Slepian process Jack Noonan a, Anatoly Zigljavsky a, a Scool of Matematics, Cardiff University, Cardiff, CF4 4AG, UK arxiv:8.0v [mat.pr] 8 Dec 08 Abstract Slepian
More informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More information1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow
1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation
More informationChapter 2 Limits and Continuity
4 Section. Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 6) Quick Review.. f () ( ) () 4 0. f () 4( ) 4. f () sin sin 0 4. f (). 4 4 4 6. c c c 7. 8. c d d c d d c d c 9. 8 ( )(
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 informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationMAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points
MAT 15 Test #2 Name Solution Guide Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points Use te grap of a function sown ere as you respond to questions 1 to 8. 1. lim f (x) 0 2. lim
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More informationALGEBRA AND TRIGONOMETRY REVIEW by Dr TEBOU, FIU. A. Fundamental identities Throughout this section, a and b denotes arbitrary real numbers.
ALGEBRA AND TRIGONOMETRY REVIEW by Dr TEBOU, FIU A. Fundamental identities Trougout tis section, a and b denotes arbitrary real numbers. i) Square of a sum: (a+b) =a +ab+b ii) Square of a difference: (a-b)
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More informationGRID CONVERGENCE ERROR ANALYSIS FOR MIXED-ORDER NUMERICAL SCHEMES
GRID CONVERGENCE ERROR ANALYSIS FOR MIXED-ORDER NUMERICAL SCHEMES Cristoper J. Roy Sandia National Laboratories* P. O. Box 5800, MS 085 Albuquerque, NM 8785-085 AIAA Paper 00-606 Abstract New developments
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationNumerical Solution of Second Order Singularly Perturbed Differential Difference Equation with Negative Shift. 1 Introduction
ISSN 749-3889 print), 749-3897 online) International Journal of Nonlinear Science Vol.8204) No.3,pp.200-209 Numerical Solution of Second Order Singularly Perturbed Differential Difference Equation with
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationSection 5.2 Series Solution Near Ordinary Point
DE Section 5.2 Series Solution Near Ordinary Point Page 1 of 5 Section 5.2 Series Solution Near Ordinary Point We are interested in second order homogeneous linear differential equations with variable
More informationApplication of Quintic B-splines Collocation Method on Some Rosenau Type Nonlinear Higher Order Evolution Equations.
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(2012) No.2,pp.142-152 Application of Quintic B-splines Collocation Metod on Some Rosenau Type Nonlinear Higer
More informationNumerical Solution to Parabolic PDE Using Implicit Finite Difference Approach
Numerical Solution to arabolic DE Using Implicit Finite Difference Approac Jon Amoa-Mensa, Francis Oene Boateng, Kwame Bonsu Department of Matematics and Statistics, Sunyani Tecnical University, Sunyani,
More informationNumerical study of the Benjamin-Bona-Mahony-Burgers equation
Bol. Soc. Paran. Mat. (3s.) v. 35 1 (017): 17 138. c SPM ISSN-175-1188 on line ISSN-0037871 in press SPM: www.spm.uem.br/bspm doi:10.569/bspm.v35i1.8804 Numerical study of te Benjamin-Bona-Maony-Burgers
More informationInvestigation of Tangent Polynomials with a Computer Algebra System The AMATYC Review, Vol. 14, No. 1, Fall 1992, pp
Investigation of Tangent Polynomials wit a Computer Algebra System Te AMATYC Review, Vol. 14, No. 1, Fall 199, pp. -7. Jon H. Matews California State University Fullerton Fullerton, CA 9834 By Russell
More informationOn convexity of polynomial paths and generalized majorizations
On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal
More informationAn Approximation to the Solution of the Brusselator System by Adomian Decomposition Method and Comparing the Results with Runge-Kutta Method
Int. J. Contemp. Mat. Sciences, Vol. 2, 27, no. 2, 983-989 An Approximation to te Solution of te Brusselator System by Adomian Decomposition Metod and Comparing te Results wit Runge-Kutta Metod J. Biazar
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More information5.1 introduction problem : Given a function f(x), find a polynomial approximation p n (x).
capter 5 : polynomial approximation and interpolation 5 introduction problem : Given a function f(x), find a polynomial approximation p n (x) Z b Z application : f(x)dx b p n(x)dx, a a one solution : Te
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationOSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix
Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More information1 1. Rationalize the denominator and fully simplify the radical expression 3 3. Solution: = 1 = 3 3 = 2
MTH - Spring 04 Exam Review (Solutions) Exam : February 5t 6:00-7:0 Tis exam review contains questions similar to tose you sould expect to see on Exam. Te questions included in tis review, owever, are
More informationMaterial for Difference Quotient
Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient
More informationNUMERICAL DIFFERENTIATION
NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationDerivatives and Rates of Change
Section.1 Derivatives and Rates of Cange 2016 Kiryl Tsiscanka Derivatives and Rates of Cange Measuring te Rate of Increase of Blood Alcool Concentration Biomedical scientists ave studied te cemical and
More informationDe-Coupler Design for an Interacting Tanks System
IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 2278-1676,p-ISSN: 2320-3331, Volume 7, Issue 3 (Sep. - Oct. 2013), PP 77-81 De-Coupler Design for an Interacting Tanks System
More informationDerivatives. By: OpenStaxCollege
By: OpenStaxCollege Te average teen in te United States opens a refrigerator door an estimated 25 times per day. Supposedly, tis average is up from 10 years ago wen te average teenager opened a refrigerator
More informationTime (hours) Morphine sulfate (mg)
Mat Xa Fall 2002 Review Notes Limits and Definition of Derivative Important Information: 1 According to te most recent information from te Registrar, te Xa final exam will be eld from 9:15 am to 12:15
More informationQuantum Theory of the Atomic Nucleus
G. Gamow, ZP, 51, 204 1928 Quantum Teory of te tomic Nucleus G. Gamow (Received 1928) It as often been suggested tat non Coulomb attractive forces play a very important role inside atomic nuclei. We can
More informationNumerical Solution of One Dimensional Nonlinear Longitudinal Oscillations in a Class of Generalized Functions
Proc. of te 8t WSEAS Int. Conf. on Matematical Metods and Computational Tecniques in Electrical Engineering, Bucarest, October 16-17, 2006 219 Numerical Solution of One Dimensional Nonlinear Longitudinal
More informationA SHORT INTRODUCTION TO BANACH LATTICES AND
CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
More informationComputer Derivations of Numerical Differentiation Formulae. Int. J. of Math. Education in Sci. and Tech., V 34, No 2 (March-April 2003), pp
Computer Derivations o Numerical Dierentiation Formulae By Jon H. Matews Department o Matematics Caliornia State University Fullerton USA Int. J. o Mat. Education in Sci. and Tec. V No (Marc-April ) pp.8-87.
More informationSECTION 1.10: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES
(Section.0: Difference Quotients).0. SECTION.0: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES Define average rate of cange (and average velocity) algebraically and grapically. Be able to identify, construct,
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
More information1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist
Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter
More informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
More informationSolutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014
Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
More informationSome Applications of Fractional Step Runge-Kutta Methods
Some Applications of Fractional Step Runge-Kutta Metods JORGE, J.C., PORTERO, L. Dpto. de Matematica e Informatica Universidad Publica de Navarra Campus Arrosadia, s/n 3006 Pamplona Navarra SPAIN Abstract:
More informationBlueprint End-of-Course Algebra II Test
Blueprint End-of-Course Algebra II Test for te 2001 Matematics Standards of Learning Revised July 2005 Tis revised blueprint will be effective wit te fall 2005 administration of te Standards of Learning
More informationA New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions
Numerisce Matematik manuscript No. will be inserted by te editor A New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions Wang Ming 1, Zong-ci Si 2, Jincao Xu1,3 1 LMAM, Scool of Matematical
More information