Mathematical Problems in Engineering
Volume 2019 (2019), Article ID 3467276, 8 pages
https://doi.org/10.1155/2019/3467276

Research Article

Similarity Solution and Heat Transfer Characteristics for a Class of Nonlinear Convection-Diffusion Equation with Initial Value Conditions

Yunbin Xu

School of Mathematics and Statistics, Yulin University, Shaanxi, Yulin 719000, China

Correspondence should be addressed to Yunbin Xu; xuyunbin886@sina.com

Received 25 February 2019; Revised 7 May 2019; Accepted 16 May 2019; Published 23 May 2019

Academic Editor: Vincenzo Vespri

Copyright © 2019 Yunbin Xu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A class of nonlinear convection-diffusion equation is studied in this paper. The partial differential equation is converted into nonlinear ordinary differential equation by introducing a similarity transformation. The asymptotic analytical solutions are obtained by using double-parameter transformation perturbation expansion method (DPTPEM). The influences of convection functional coefficient and power law index on the heat transport characteristics are discussed and shown graphically. The comparison with the numerical results is presented and it is found to be in excellent agreement. The method and technique used in this paper have the significance in studying other engineering problems.

1. Introduction

Convection-diffusion equations arose from various fields of applied sciences such as heat transfer problems in a draining film [1] or a nanofluid filled enclosure [2], radial transport in a porous medium [3], and water transport in soils [4] and have received extensive attentions during the past several decades. For convection-diffusion equations, only few cases with special initial or boundary value conditions have analytical solutions. Therefore, most of the main concerns were the study of the qualitative properties of the solutions [5–13] and numerical study [14–23] for convection-diffusion equations. However, the very important approximate solution of convection-diffusion equation has not been well solved. In this paper we present similarity solutions for a class of convection-diffusion equation, which are then solved by DPTPEM. The method was first proposed by Yan Zhang et al. [24, 25]. In this paper, we can find that the DPTPEM can be successfully used to solve the studied problem, and the approximate solutions obtained by DPTPEM agree very well the numerical solutions obtained by bvp4c with Matlab [26].

2. Mathematical Formulation

In this paper, we consider a class of nonlinear convection-diffusion equation as follows [20]:

subject to the following initial value conditions:

where , is positive power law exponent. Convection functional coefficient is assumed to be a real-valued continuous differential function defined on , , for . is the heat density per unit area.

3. Converting into Nonlinear Ordinary Differential Equation Boundary Value Problems

Introduce the following similarity variables [20]:

The transformed convection-diffusion equation together with the initial value conditions given by (1)-(2) can be written as

Conversely, if is a solution of (4)-(5), then must be a similarity solution to (1)-(2).

Let be a solution to (4)-(5). If is strictly increasing in , then the function inverse to exists. Here holds on . Then (4) can be changed into the following nonlinear ordinary differential equation:

and integration from to 1 yields

Set

Combine (7) and (8) to get

The corresponding boundary conditions are

where is heat diffusion flux. The derivation process indicates that only the positive solutions of (9)-(10) are physically significant.

4. A Brief Introduction to DPTPEM

The method was first proposed by Yan Zhang and Liancun Zheng [24, 25].

We consider a class of ordinary differential equation initial or boundary value problem without small parameter as follows:

where is order differential operator without small parameter , is a constant, and is the region containing the origin.

Introduce an embedded parameter transformation as follows:

where is an artificial small parameter, is independent variable parameter, is dependent variable parameter, and and are undetermined constants. Substituting this transformation (12) into (11), we can obtain the following nonlinear initial value problem:

where is a function of , is the region containing the origin. The solution of is expressed in form of power series by applying perturbation expansion method

Combine (12) and (14) to get the solution of (11) in terms of power series. The undetermined parameters are determined according to the boundary conditions of the original problem.

The basic idea of the DPTPEM consists of three steps. Firstly, introducing an artificial small parameter , the independent variable and the unknown function are transformed simultaneously, and the problem is transformed into a new one related to small parameter . Secondly, the transformed new differential equation is expanded in the form of power series of and decomposed into the sum of several solution components. Then, by using known initial or boundary value conditions, we try to find out the solution components of each order separately and combine them to get the solution of the new problem. Thirdly, by substituting the solution of the new equation into the transformation equation (12), the small parameter is eliminated and the solution of the original problem is obtained.

5. Approximate Analytical Results of (9)-(10) and Discussion

The initial condition is assumed as follows:

In order to solve (9), we transform the dependent variable and independent variable as follows according to formula (12):

where is an artificial small parameter.

Let and substituting this transformation (16) into (9)-(10), we can obtain the following nonlinear initial value problems:

Expanding in a power series development, we can get

The solution of (19) can be obtained by expanding in a power series development near as follows:

Substituting (20) into (19) and equating the coefficients of , we can get the following expressions:

Then, we can obtain in the following.

Combine (16), (20), and the above ,… to get solution of (9)-(10).

Similarly, let ; we obtain approximate analytical solution of (9)-(10) as follows:

Let ; approximate analytical solution of (9)-(10) is presented as follows:

It is obvious that we can promptly obtain the value of by applying in (23)-(25) for each fixed . Based on (23)-(25) and the corresponding we can easily obtain the graph of the heat diffusion flux distribution for different and . The results are presented in Tables 1–5 and Figures 1–5.

Table 1: Values of for , , , , , , and .

0.2

1.013

0.6

0.7815

1.0

0.6792

2.5

0.5701

5.0

0.5342

10.0

0.5169

Table 2: Values of for , and .

0.2

0.9122

0.6

0.6532

1.0

0.5352

2.5

0.4090

5.0

0.3691

10.0

0.3506

Table 3: Values of for , and .

0.2

0.8643

0.6

0.5927

1.0

0.4670

2.5

0.3300

5.0

0.2870

10.0

0.2676

Table 4: Values of for and .

2.0

0.5885

2.0

0.4300

2.0

0.3529

Table 5: Values of for and .

8.0

0.5212

8.0

0.3551

8.0

0.2723

Figure 1: Flux distribution for .

Figure 2: Flux distribution for .

Figure 3: Flux distribution for .

Figure 4: Flux distribution for

Figure 5: Flux distribution for

It can be seen from Figures 1–3 that the heat diffusion flux decreases with increase of for specific , the physical meaning is that heat diffusion flux is a decreasing function of , which means that the profiles exhibited by a big power law index possess a smaller diffusion. Figures 4 and 5 indicate the heat diffusion flux decrease with the decrease of .

In order to verify the efficiency and reliability of approximate analytical solution obtained by using DPTPEM, a comparison of approximate analytical solution and numerical one obtained by bvp4c with Matlab is presented in Tables 6 and 7 and Figures 6 and 7. It is obvious that excellent agreement exists for approximate analytical solution and numerical one.

Table 6: Comparison of for and .

0.0
0.5169
0.5169

0.1
0.5119
0.5119

0.2
0.4968
0.4968

0.3
0.4715
0.4715

0.4
0.4359
0.4359

0.5
0.3899
0.3899

0.6
0.3334
0.3334

0.7
0.2663
0.2663

0.8
0.1884
0.1884

0.9
0.0997
0.0997

1.0
0
0

Table 7: Comparison of for and .

0.0
0.3529
0.3551

0.1
0.3528
0.3550

0.2
0.3518
0.3540

0.3
0.3483
0.3506

0.4
0.3405
0.3427

0.5
0.3255
0.3278

0.6
0.3000
0.3023

0.7
0.2601
0.2624

0.8
0.2008
0.2031

0.9
0.1164
0.1186

1.0
0
0

Figure 6: The comparison of the results of approximate analytical solution and numerical solution for

Figure 7: The comparison of the results of approximate analytical solution and numerical solution for

6. Conclusions

In the paper, a class of nonlinear convection-diffusion was studied. The partial differential equation and corresponding initial value conditions were transformed into a class of singular nonlinear boundary value problems of ordinary differential equation when similarity variables were introduced. An efficient approximate analytical method named DPTPEM was applied to solve these nonlinear problems. The reliability and effectiveness of the DPTPEM were verified by comparing approximate results with the numerical solutions. The effects of convection functional coefficient and power law index on transfer behavior were presented. The results show that heat diffusion flux is an increasing function of convection functional coefficient and a decreasing function of power law index .

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The author declares that there are no conflicts of interest.

Acknowledgments

The research was supported by a grant from the National Natural Science Foundation of China (no. 11501496) and the Natural Science Basic Research Plan in Shaanxi Province of China (no. 2014JQ2-1003).

References

J. Isenberg and C. Gutfinger, “Heat transfer to a draining film,” International Journal of Heat and Mass Transfer, vol. 16, no. 2, pp. 505–512, 1973.

M. Sheikholeslami, R. Ellahi, M. Hassan, and S. Soleimani, “A study of natural convection heat transfer in a nanofluid filled enclosure with elliptic inner cylinder,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 24, no. 8, pp. 1906–1927, 2014.

E. J. Veling, “Radial transport in a porous medium with dirichlet, neumann and robin-type inhomogeneous boundary values and general initial data: analytical solution and evaluation,” Journal of Engineering Mathematics, vol. 75, pp. 173–189, 2012.

J.-Y. Parlange, “Water transport in soils.,” Annual review of fluid mechanics, Vol. 12, pp. 77–102, 1980.

X. Yu and Q. S. Zhang, “Asymptotically self-similar solution for the convection-diffusion equation,” Applied Mechanics and Materials, vol. 623, pp. 97–100, 2014.

M. o. Krupski, “Convection-diffusion equations with random initial conditions,” Journal of Mathematical Analysis and Applications, vol. 470, no. 2, pp. 1194–1221, 2019.

B. I. Olajuwon and A. O. Popoola, “On self-similar solution of non-linear diffusion equation with convection term,” Journal of Mathematics and Statistics, vol. 2, no. 4, pp. 453–456, 2006.

G. Lu and L. Jiang, “A remark on removable singularity for nonlinear convection-diffusion equation,” Nonlinear Analysis: Real World Applications, vol. 27, pp. 1–25, 2016.

L. I. Ignat and D. Stan, “Asymptotic behavior of solutions to fractional diffusion-convection equations,” Journal Of The London Mathematical Society-Second Series, vol. 97, no. 2, pp. 258–281, 2018.

A. Gmira and L. Veron, “Large time behaviour of the solutions of a semilinear parabolic equation in ,” Journal of Differential Equations, vol. 53, no. 2, pp. 258–276, 1984.

M. Escobedo and E. Zuazua, “Large time behavior for convection-diffusion equations in RN,” Journal of Functional Analysis, vol. 100, pp. 119–161, 1991.

K. Kalli and K. N. Soltanov, “Existence and behavior of solutions for convection-diffusion equations with third type boundary condition,” TWMS Journal of Pure and Applied Mathematics, vol. 8, no. 2, pp. 209–222, 2017.

L. Greco, G. Moscariello, and G. Zecca, “Very weak solutions to elliptic equations with singular convection term,” Journal of Mathematical Analysis and Applications, vol. 457, no. 2, pp. 1376–1387, 2018.

S. V. Siryk, “A note on the application of the Guermond-Pasquetti mass lumping correction technique for convection-diffusion problems,” Journal of Computational Physics, vol. 376, pp. 1273–1291, 2019.

R. Bürger and I. Kröker, “Computational uncertainty quantification for some strongly degenerate parabolic convection-diffusion equations,” Journal of Computational and Applied Mathematics, vol. 348, pp. 490–508, 2019.

R. Costa, J. M. Nóbrega, S. Clain, G. J. Machado, and R. Loubère, “Very high-order accurate finite volume scheme for the convection-diffusion equation with general boundary conditions on arbitrary curved boundaries,” International Journal for Numerical Methods in Engineering, vol. 117, no. 2, pp. 188–220, 2019.

Z. Chen, S.-H. Chou, and D. Y. Kwak, “Characteristic-mixed covolume methods for advection-dominated diffusion problems,” Numerical Linear Algebra with Applications, vol. 13, no. 9, pp. 677–697, 2006.

S. May, “Spacetime discontinuous Galerkin methods for convection-diffusion equations,” Bulletin of the Brazilian Mathematical Society, vol. 47, no. 2, pp. 561–573, 2016.

A. V. Gladky, “Stability of difference splitting schemes for convection diffusion equation,” Cybernetics and Systems Analysis, vol. 53, no. 2, pp. 193–203, 2017.

L. C. Zheng, X. X. Zhang, and J. C. He, “Free boundary value problems for a class of generalized diffusion equation,” Journal of University of Science and Technology Beijing, vol. 9, no. 6, pp. 422–425, 2002.

M. R. Sidi Ammi and I. Jamiai, “Finite difference and Legendre spectral method for a time-fractional diffusion-convection equation for image restoration,” Discrete and Continuous Dynamical Systems - Series S, vol. 11, no. 1, pp. 103–117, 2018.

R. C. Lin, X. Ye, S. Zhang, and P. Zhu, “A weak Galerkin finite element method for singularly perturbed convection-diffusion-reaction problems,” SIAM Journal on Numerical Analysis, vol. 56, no. 3, pp. 1482–1497, 2018.

M. Saqib, S. Hasnain, and D. S. Mashat, “Computational solutions of two dimensional convection diffusion equation using crank-nicolson and time efficient ADI,” American Journal of Computational Mathematics, vol. 7, no. 3, pp. 208–227, 2017.

Y. Zhang, L. C. Zheng, and X. X. Zhang, “Analytical approximate solution for Marangoni convection in a liquid layer with coupled boundary conditions,” Acta Physica Sinica, vol. 58, no. 8, pp. 5501–5506, 2009.

L. C. Zheng and X. X. Zhang, Modern Analysis Methods in Nonlinear Partial Differential Equations, Science Press, Beijing, China, 2011.

X. Si, H. Li, Y. Shen, and L. Zheng, “Effects of nonlinear velocity slip and temperature jump on pseudo-plastic power-law fluid over moving permeable surface in presence of magnetic field,” Applied Mathematics and Mechanics-English Edition, vol. 38, no. 3, pp. 333–342, 2017.