_{1}

The solution of problem convection-diffusion equation by way of control vol-ume method is considered. Approximate solution of problem is received. Three point scheme of high resolution is constructed.

Many physical problems involve the combination of convective and diffusive processes. Convection-diffusion problems arise frequently in many areas of applied sciences and engineering. They occur in fields where mathematical modeling is important such as physics, engineering and particularly in fluid dynamics and transport problems.

Mathematical models of physical, chemical, biological and environmental phenomena are governed by various forms of differential equations.

Convection-diffusion problems are governed by typical mathematical models, which are common in fluid and gas dynamics. Heat and mass transfer is conducted not only via diffusion, but appears due to motion of a medium, too.

When velocity is higher, that is, flow term is larger, a simple convection- diffusion problem is converted to convection-dominated diffusion problem because Peclet number is greater than two. A Peclet number is a dimensionless number relevant in the study of transport phenomena in fluid flows. It is defined to be the ratio of the rate of advection of a physical quantity by the flow to the rate of diffusion of the same quantity driven by an appropriate gradient. For diffusion of heat (thermal diffusion), Peclet number is defined as, Pе= F/D where F is the flow term and D is the diffusion term.

Finite volume methods are widely used in computational fluid dynamics. The elementary finite volume method uses a cell-centered mesh and finite-difference approximations of first order derivatives. This paper shows how the finite volume method is applied to a simple model of convective transport: the one- dimensional convection-diffusion equation.

There are two primary goals of this paper. The first is to apply the finite volume method to obtain approximately analytical solution. The second one is construct scheme which works for all mesh Peclet number. Readers interested in additional details, including application to the Navier-Stokes equations, should consult the classic text by Patankar [

Usually, solution of differential equations by numerical methods is obtained in the form of numbers. Here we will show a possibility of deriving of a solution of differential equations by control volume methods in the approximately-analytical form by the so-called moved node.

Such type research it is considered in work [

Let’s consider one-dimensional the convection-diffusion equation on a interval [ W , E ] :

d d x ( ρ u Φ ) = d d x ( Γ d Φ d x ) + S ( x ) (1)

with boundary conditions

d Φ ( W ) / d x = a 0 Φ ( W ) + a 1 , d Φ ( E ) / d x = b 0 Φ ( E ) + b 1 (2)

where u a stream velocity in a x direction, ρ a stream denseness, Γ ―a diffusivity, S ( x ) ―a given function (source), Φ ―unknown function. From an equation of continuity implies that F = ρ u = c o n s t .

Let’s consider the Equation (1) on segments [ W , E ] . For deriving of the approached analytical solution of a problem by the method of control volume we take an arbitrary point x ∈ [ W , E ] and control volume [ w , e ] . In

Let’s suppose that, the edge w is arranged in the middle between points W and x, and an edge e―in the middle between points x and E. Integrating the Equation (1) on control volume and substituting, derivatives the upwind scheme we will receive zero approach

( а Е + а W ) Φ 0 = а Е Φ E 0 + а W Φ W 0 + E − W 2 ⋅ S ( x ) (3)

here a E = Γ Е − х + max ( − F , 0 ) ; a W = Γ x − W + max ( F , 0 ) . As x ∈ [ W , E ] an arbitrary point from (3) we can define Φ 0 and receive the approached analytical solution of a problem (1) ( Φ W 0 , Φ E 0 ) are defined on a basis (2)). The source term is obtained by assuming that S has the uniform value of the control volume.

To improve approximate solution we take additional grids: x 1 = x + W 2 , x 2 = x + E 2 . Let’s write the upwind scheme of type (3) for a segment [ W , x ] , [ x 1 , x 2 ] and [ x , E ] . Let’s receive system from three equations. We exclude the received system, Φ 1 ( x 1 ) , Φ 1 ( x 2 ) and as a result we will receive the improved scheme:

[ β 1 + 1 + τ 1 + α 1 − 1 + γ 1 ] Φ 1 = β 1 + 1 + τ 1 Φ W 1 + α 1 − 1 + γ 1 Φ E 1 + E − W 4 ⋅ S ( x ) + 1 1 + τ 1 ⋅ x − W 2 ⋅ S ( W + x − W 2 ) + 1 1 + γ 1 ⋅ E − x 2 ⋅ S ( x + E − x 2 ) , (4)

where τ 1 = β 1 − / β 1 + , γ 1 = α 1 + / α 1 − , β 1 − = 2 D W + F − , β 1 + = 2 D W + F + , α 1 − = 2 D E + F − , α 1 + = 2 D E + F + , D E = Γ / ( E − x ) , D W = Γ / ( x − W ) , F − = max ( − F , 0 ) , F + = max ( F , 0 ) . In (4) Φ 1 improved value of unknown function in a point x ( Φ W 1 ≡ Φ W , Φ E 1 ≡ Φ E ) , ( Φ W 1 , Φ E 1 are defined on the basis of (2)).

Solving (4) rather Φ 1 again for solution improving we will arrive similarly: write the scheme (4) for a segment [ W , x ] , [ x 1 , x 2 ] and [ x , E ] , and receive the improved analytical solution. Also we will exclude unknowns in points x 1 and, x 2 etc. Continuing this process we will obtain

[ ( 1 − τ k ) β k + 1 − τ k 2 k + ( 1 − γ k ) α k − 1 − γ k 2 k ] Φ k = ( 1 − τ k ) β k + 1 − τ k 2 k Φ W k + ( 1 − γ k ) α k − 1 − γ k 2 k Φ E k + E − W 2 k + 1 ⋅ S ( x ) + 1 − τ k 1 − τ k 2 k ⋅ x − W 2 k ⋅ ∑ j = 1 2 k − 1 ∑ i = 1 j τ k i − 1 S ( W + j x − W 2 k ) + 1 − γ k 1 − γ k 2 k ⋅ E − x 2 k ⋅ ∑ j = 1 2 k − 1 ∑ i = 1 j γ k i − 1 S ( x + ( 2 k − j ) E − x 2 ) , (5)

where τ k = β k − / β k + , γ k = α k + / α k − , β k − = 2 k D W + F − , β k + = 2 k D W + F + , α k − = 2 k D E + F − , α k + = 2 k D E + F + .

In (5) Φ k improved value of unknown function in a point x ( Φ W k ≡ Φ W , Φ E k ≡ Φ E ) . Solving the Equation (5) rather we will receive the approached analytical solution of an initial problem.

On

On

The analytical scheme (5) allows not only receiving the approached analytical solution, but also gives the chance in creation of the qualitative scheme.

Comparisons of exact and difference solutions are on

Solution of problem 50 Φ ′ ( x ) = Φ ″ ( x ) + S ( x ) are on

In

Peclet number). Solid line is exact solutions of the problem. Rectangles under the scheme of Patankar, diamond at k = 2 and circles at k = 6 . From

Comparison of analytical solutions received finite difference [

The approached analytical solution to one-dimensional convection-diffusion problems is received. The algorithm of deriving of a solution is based on the control volume method. The analytical solution allows also constructing the compact scheme. It is the three point scheme of the high order of resolution. It is the scheme that is suitable for any Peclet numbers.

The author declares no conflicts of interest regarding the publication of this paper.

Umurdin, D. (2018) Computing Technology of a Method of Control Volume for Obtaining of the Approximate Analytical Solution to One- Dimensional Convection-Diffusion Problems. Open Access Library Journal, 5: e4962. https://doi.org/10.4236/oalib.1104962