﻿
Reach Us  +44-7482-875032
The Finite Difference Methods for Parabolic Partial Differential Equations

# Journal of Applied & Computational MathematicsOpen Access

All submissions of the EM system will be redirected to Online Manuscript Submission System. Authors are requested to submit articles directly to Online Manuscript Submission System of respective journal.

# The Finite Difference Methods for Parabolic Partial Differential Equations

Nigatie Y*
Department of Mathematics, College of Natural and Computational Sciences, University of Gondar, Gondar, Ethiopia
*Corresponding Author: Nigatie Y, Department of Mathematics, College of Natural and Computational Sciences, University of Gondar, Gondar, Ethiopia, Tel: +251913080434, Email: [email protected]

Received Date: Aug 18, 2018 / Accepted Date: Sep 12, 2018 / Published Date: Sep 27, 2018

Keywords: Parabolic; Heat equation; Finite difference; Bender- Schmidt; Crank-Nicolson

#### Introduction

Parabolic partial differential equations

The well-known parabolic partial differential equation is the one dimensional heat conduction equation . The solution of this equation is a function u(x,t) which is defined for values of x from 0 to l and for values of t from 0 to ∞ [2-4]. The solution is not defined in a closed domain but advances in an open ended region from initial values, satisfying the prescribed boundary conditions (Figure 1). In general the study of pressure waves in a fluid, propagation of heat and unsteady state problems lead to parabolic type of equations .

Solution of one-dimensional Heat conduction equation

Consider the heat conduction which is of the form 

cut=uxx (1)

where C is a constant Let the plane (x,t) be divided in to smaller rectangles by means of the sets of lines 

x=ih,i=0,1,2,.

t=jk,j=0,1,2,.

Using the finite difference approximations [8-11] (2)

and (3)

Hence the finite difference analogue in eqn. (1) is (4)

Which can be rewritten as

ui,j+1=ui,j+λ(ui−1,j −2ui,j+ui+1,j)

where this formula expresses the unknown function value at the (i,j+1)th interior point in terms of the known function values and hence it is called the explicit formula. It can be shown that this formula is valid only  for For becomes in eqn. (4) (5)

The above formula in eqn. (5) is called Bender-Schmidt recurrence relation. In formula in eqn. (4), we have used the function values along the jth row only in the approximation of uxx is replaced by the  average of its finite difference approximations on the jth and (i,j+1)th rows. Thus and hence in eqn. (1) is replaced by This implies (6)

where , on the left side in eqn. (6) we have three unknowns and on the right side all the three quantities are known. The above formula in eqn. (6) which is an implicit scheme is called Crank-Nicolson formula and it is convergent for all finite values of λ. If there are N internal mesh points on each row, then the formula in eqn. (6) gives N simultaneous equations for the N unknowns in terms of the given boundary values. Similarly the internal mesh points on all rows can be calculated .

Example 1: Use the Bender-Schmidt recurrence relation to solve the equation uxx=2ut. With the conditions u(x,0)=4x−x2,u(0,t)=u(4,t)=0 Solution Taking h=1 we obtain Also u(0,0)=0,u(1,0)=3,u(2,0)=4,u(3,0)=3 and u(4,0)=0. For this the first time step, k=1.

For this the first time step, k=1.

Using Bender-Schmidt recurrence relation, we obtain For k=2, we have For k=3 we have For k=4 we have Similarly with k=5 we obtain The computation can be continued to any number of time steps.

Example 2: Solve the equation ut=uxx, subject to the conditions,

u(x,0)=sinπx,0 ≤ x ≤ 1,u(0,t)=u(1,t)=0.

Using

(a) Bender-Schmidt Method

(b) Crank-Nicolson Method

Carry out the computations for two levels, taking Solution

Here and so that Also and all boundary values are zero in Figure 2.

(a) Using Bender-Schmidt formula

Using in eqn. (4) of this section

we have ui,j+1=λui−1,j+(1−2λ)ui,jui+1,j

This becomes For i=1,2; j=0 (b) By using Crank-Nicolson method

From eqn. (6) above we have For i=1,2; j=0 Solving these equations, we find

u1,1=u2,1=0.67

For i=1,2; j=1

−u0,2+10u1,2−u2,2=u0,1+ 6u1,1+ u2,1

⇒ 10u1,2−u2,2=4.69

−u1,2+10u2,2−u3,2=u1,1+6u2,1+ u3,1

⇒−u1,2+10u2,2=4.69

Again solving the above two equation we can obtain

u1,2=u2,2=0.52.

Solution of two dimensional heat equations

ut =c2(uxx+uyy).

The method employed for the solution of one dimensional heat equation can be readily extended to the solution two dimensional heat equations in eqn. (7).

Consider a square region 0 ≤ x ≤ y ≤ a and assume that u is known at all points within and on the boundary of this square.

If h be the step size then a mesh point

(x,y,t) =(ih,jh,nl) may be denoted as (i,j,n).

Replacing the derivatives in eqn. (7) by their finite difference approximations, we get (8)

Where . This equation needs the five points available on the nth plane (Figure 3).

The computation process consists point by point evaluation in the (n+1)th plane using the points on the nth plane. This method is known as ADE (Alternating Direct Explicit) method.

Example: Solve the equation ut=uxx+uyy.

Subject to the initial conditions u(x,y,0) =sin2πxsin2πy,0 ≤ x,y ≤ 1 and the conditions u(x,y,t)=0, t ≥ 0, on the boundaries using ADE method and (Calculate the results for one time level).

Solution: In eqn. (8) becomes The mesh points and the computation model is given in Figure 4 below

At the zeroth level (n=0), the initial and boundary conditions are And

ui,0,0=u0,j,0+u3,j,0=ui,3,0=0;i,j=0,1,2,3

Now we can calculate the mesh values at the first level.

For n=0, from (∗) we obtain (i) Putting i=j=1 in (∗∗) above (ii) Putting i=2,j=1 in(∗∗) above (iii) Putting i=1,j=2 in (∗∗) above (iv) Putting i=j=2 in (∗∗) above Similarly the mesh values at the second and higher levels can be calculated.

#### References

Citation: Nigatie Y (2018) The Finite Difference Methods for Parabolic Partial Differential Equations. J Appl Computat Math 7: 418. DOI: 10.4172/21689679.1000418

Copyright: © 2018 Nigatie Y. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Select your language of interest to view the total content in your interested language

##### Recommended Journals
Viewmore
###### Article Usage
• Total views: 1125
• [From(publication date): 0-0 - Aug 23, 2019]
• Breakdown by view type
• HTML page views: 1020
• PDF downloads: 105

## Post your comment  Can't read the image? click here to refresh
###### Peer Reviewed Journals

Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals

International Conferences 2019-20

Meet Inspiring Speakers and Experts at our 3000+ Global Annual Meetings

Top