We establish the conditions for the compute of the stability restriction and local accuracy on the time step and we prove the consistency and local truncation error by using
θ-scheme and 3-level scheme for Heat Equation with smooth initial conditions and for some parameter
θ∈[0,1].

Global Truncation Local Accuracy Stability Restriction1. Introduction

In this paper we have considered the heat equation with. Using -scheme and 3-level scheme in space we compute the order of local accuracy in space and time and stability restriction as a function of on the time step. Much attention has been paid to the development, analysis and implementation of accurate methods for the numerical solution of this problem in the literature. Many problems are modeled by smooth initial conditions and Dirichlet boundary conditions. A number of procedures have been suggested (see, for instance [1] - [3] ). We can say that three classes of solution techniques have emerged for solution of PDE: the finite difference techniques, the finite element methods and the spectral techniques (see [4] and [5] ). The last one has the advantage of high accuracy attained by the resulting discretization for a given number of nodes [6] - [8] .

We consider Scheme (1) for the 1D heat equation for some parameter. We compute the order of local accuracy in space and time as a function of and its the stability restriction. Until, we compute the solution with some fixed error with the smallest amount of CPU time, and finally we can see this findings producing the relevant convergence and efficiency plot. For the 3-level scheme we consider (11) for the 1D heat equation and we compute the local truncation error. For different values of and we find the stability criterion of the scheme and its accuracy.

be the -scheme applied to the one-dimensional heat equation

Now for the order of local accuracy in space and time as a function of we write the local truncation error. In time we have

where represents the exact solution of the heat equation. Now we perform Taylor expansion of at.

We can write the LHS of (1) as

here represents the derivative with respect to time, of. On the RHS, we have a centered difference approximating second derivative of

As we are solving the heat equation, the previous expression is

Now, at time we have

therefore, applying Taylor expansion with respect to we can write

So (8) becomes

Here RHS of (1) becomes

After the elimination of some terms we have

Now simplifying we obtain

Cancelling and moving all terms to the right side, we get

Scheme (10) is first order in time, second order in space. If for example, it becomes second order, this is due to cancellation of the.

Stability Restriction as a Function of

Here we will apply Von Neumann stability. Let and. Then Equation (1) can be written as

Now dividing by we have

Therefore by using and we can rewrite the expression as

By using the identity we have

We can say this scheme is stable only for. Now, let. The inequality is

thus

Now multiplying by the denominator we have

The expression in the absolute value becomes

Therefore by the Von Neumann stability condition, the scheme is stable if.

In this case we can say the following about the best combination for and. In order to have both local accuracy and stability, the optimal value of is and therefore this scheme represents the Crank-Ni- cholson scheme. Here and appear in the form.

In Figure 1 the convergence plot equation (varying the radio r) is

with matrix A described in the heat equation. We can say the scheme is unconditionally stable. We can see in Figure 1 that we have a linear convergence with respect to r.

3. Three-Level Scheme

We start by computing the stability restriction one has to impose on. We apply Von Neumannstability analysis to the scheme.

Let

where

E vs. r for 1D-heat equation, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7402655x86.png" xlink:type="simple"/></inline-formula>with initial temperature <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7402655x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7402655x87.png" xlink:type="simple"/></inline-formula> on<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7402655x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7402655x87.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7402655x88.png" xlink:type="simple"/></inline-formula>

and

By using (12) and (13) we can rewrite (11) as

or as

The local truncation error for this scheme is as follow.

where represents the exact solution of the heat equation. Therefore we have

Now expanding operator on the left side, we can isolate the forward difference in time at, then

however,

Expanded this expression becomes

Finally we have

4. Stability Criterion for the Three-Level Scheme and Its Accuracy When <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7402655x103.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7402655x103.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7402655x104.png" xlink:type="simple"/></inline-formula>

By using Equation (14) we have

Now applying Von Neumann stability again, the aim is to use and , therefore

Multiplying both sides by and write we obtained

Using the cosine identity that we have

We have a quadratic equation in, where, therefore

After some cancellations, we can write

Here, if all we need

Acknowledgements

We would like to thank the referee for his valuable suggestions that improved the presentation of this paper.

ReferencesStrikwerda, J.C. (1989) Finite Difference Schemes and Partial Differential Equations. Wadsworth & Brooks, Siam, Madison, Wisconsin.McRea, G.J. and Godin, W.R. (1967) Numerical Solution of Atmospheric Diffusion for Chemically Reacting Flows. Journal of Computational Physics, 77, 1-42.Cárdenas Alzate, P.P. (2014) A Survey of the Implementation of Numerical Schemes for Linear Advection Equation. Advances in Pure Mathematics, 4, 467-479.
http://dx.doi.org/10.4236/apm.2014.48052Cárdenas Alzate, P.P. (2014) A Survey of the Implementation of Numerical Schemes for the Heat Equation Using forward Euler in Time. Journal of Applied Mathematics and Physics, 2, 1153-1158.
http://dx.doi.org/10.4236/jamp.2014.213135Hundsdorfer, W. and Koren, B. (1995) A Positive Finite-Difference Advection Scheme Applied on Locally Refined Grids. Journal of Computational Physics, 117, 35-36.
http://dx.doi.org/10.1006/jcph.1995.1042Canuto, C. and Hussaini, M. (1988) Spectral Methods in Fluids Dynamics. Springer Series in Computational Physics, Springer-Verlag, Berlin.
http://dx.doi.org/10.1007/978-3-642-84108-8Dehghan, M. (2007) The One-Dimensional Heat Equation Subject to a Boundary Integral Specification. Chaos, Solitons & Fractals, 32, 661-675.
http://dx.doi.org/10.1155/MPE.2005.61Lu, X., et al. (2005) A New Analytical Method to Solve the Heat Equation for a Multi-Dimensional Composite Slab. Journal of Physics, 38, 2873.