Advances in Mathematical Physics
Volume 2019 (2019), Article ID 7351042, 7 pages
https://doi.org/10.1155/2019/7351042

Research Article

Flux Transport Characteristics of Free Boundary Value Problems for a Class of Generalized Convection-Diffusion Equation

Yunbin Xu and Meihua Wei

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

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

Received 17 December 2018; Accepted 15 January 2019; Published 3 February 2019

Academic Editor: Yannis Dimakopoulos

Copyright © 2019 Yunbin Xu and Meihua Wei. 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.

The similarity transformation is introduced for studying free boundary value problems for a class of generalized convection-diffusion equation. A class of singular nonlinear boundary value problems are obtained and solved by using Adomian decomposition method (ADM). The approximate solution can be expressed in terms of a rapid convergent power series with easily computable terms. The efficiency and reliability of the approximate solution are verified by numerical ones. The effects of the variable thermal conduction , convection functional coefficient , power law exponent , and parameter on the flux transport characteristics are presented graphically and analyzed in detail.

1. Introduction

Convection-diffusion equation is a class of very important equations; it often appears in many research fields such as hydrodynamics, transport, electronics, energy, and environmental science [1–4]. In the process of dealing with practical problems, for many mathematical models, especially partial differential equations, it is difficult to obtain their analytical solutions in general. On the one hand, some scholars consider the existence, uniqueness, or nonuniqueness of solutions for the convection-diffusion equations (for example, see [5–13]). On the other hand, others focus on the numerical solution of the convection-diffusion equation by all kinds of methods, for instance, the spectral element method [14], the finite element method [15–17], the finite difference method [18, 19], and the Runge-Kutta method [20].

However, all of the above-mentioned problems paid attention only to the qualitative properties of the solutions, such as the existence and uniqueness, or of the numerical ones. The very important approximate solution [21, 22] of convection-diffusion equation has not been well solved. In this paper we present similarity solutions for generalized convection-diffusion equation with free boundary conditions, which are then solved using ADM. ADM [23–25] has been shown as a useful way of obtaining accurate and computable solutions to operator equations involving nonlinear terms. The characteristic of ADM is to decompose the nonlinear terms in the equations into a peculiar series of polynomials which are the so-called Adomian polynomials. The solution of the equations is then considered as the sum of a series rapidly converging to an accurate solution. The convergence analysis and proof of ADM for nonlinear problems have been studied by some scholars (for example, see [26–30]). In this paper, we can find that the approximate solution agrees very well with the numerical solution, which shows the reliability and validity of the present work. The effects of the convection functional coefficient , variable thermal conduction , and power law index on the flux transport characteristics are discussed by graph in detail.

2. Mathematical Formulation

In the paper, we consider a class of generalized convection-diffusion equation with free boundary conditions as follows:

Here , is an unknown function and must be determined as a part of the solution; is given positive power law exponent; and is given free parameter. is assumed to be a real-valued continuous differential function defined on , , for . and are supposed to be real-valued nonnegative continuous differential functions defined on is convection functional coefficient, and is variable thermal conduction. is the heat density per unit area.

3. Transformation of (1)–(4)

Introducing the similarity transformation [31–34],

Substituting (5) into (1)–(4) and choosing , , , the nonlinear ordinary differential equation boundary value problem with respect to is derived as

Conversely, if is a solution of (6)–(8), then must be a similarity solution to (1)–(4).

Let be a solution to (6)–(8) and it is strictly decreasing in ; then its inverse function exists. Setting , then (6)–(8) can be converted into the following singular nonlinear ordinary differential equation of two-point boundary value problem

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

4. Adomian Approximate Results and Discussion

Integrate equation (9) twice from 0 to for and to obtain

where and is a twofold integral operator. Letting

is presented as follows according to ADM:

The iterative formula of is presented as follows:

where . We can obtain Adomian polynomial for and according to the iterative formula of as follows:

Substituting the derived Adomian polynomial (15) into (13), we obtain

We express the solution of (10) and (13) by a four-ordered approximate solution and the higher order can be determined in a parallel manner.

It is obvious that we can promptly obtain the value of by applying in (17) for each fixed and each fixed . Based on (17) and the corresponding we can easily obtain the graph of the diffusion flux distribution for , , and different values of and . The results are given in Figures 1–3.

Figure 1: Flux distribution for

Figure 2: Flux distribution for

Figure 3: Flux distribution for

It can be seen from Figure 1 that the diffusion flux decreases with increase of for , , and specific . Figure 2 indicates that for , , and specific , the diffusion flux appears with an utterly different characteristic with the changing .There exists a subinterval such that diffusion flux increases with the increase of for but decreases with the increase of for . Figure 3 describes the behavior of with the changes of parameter ; from Figure 3 we note that the diffusion flux increases sharply with the increasing of .

Similarly, (9)-(10) can be solved for , , . The corresponding approximate solutions are, respectively, given in (18)–(20).

It is very easy to obtain the value of by applying in (18)–(20). Based on (18)–(20) and the corresponding we can easily describe the graph of the diffusion flux distribution for and . The results are illustrated in Figure 4.

Figure 4: Flux distribution for , , , .

Considering , , we have . Figure 4 indicates that the diffusion flux decreases with the decrease of .

Similarly, it is easy to solve (9)-(10) for , , , . The corresponding approximate solutions are, respectively, given in (18) and (21)-(22).

We can obtain easily the value of by employing in (21)-(22). Based on (18), (21)-(22), and the corresponding , it is easy to get the graph of the diffusion flux distribution for and , . The results are illustrated in Figure 5.

Figure 5: Flux distribution for , , , .

Considering for and , it is seen that from Figure 5 the diffusion flux increases with the decrease of in the domain near zero but increases with the increase of in the domain near one.

In order to verify the efficiency and reliability of approximate solutions, approximate solutions obtained by the ADM are compared with the numerical ones obtained by the finite difference method in Figures 6–11. It is seen that the approximate solution is highly in agreement with the numerical solution.

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

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

Figure 8: The comparison of the results of approximate solution and numerical solution for , , ,

Figure 9: The comparison of the results of approximate solution and numerical solution for , , , .

Figure 10: The comparison of the results of approximate solution and numerical solution for , , , .

Figure 11: The comparison of the results of approximate solution and numerical solution for , , , .

5. Conclusions

This paper presented a similarity analysis for a class of generalized convection-diffusion with free boundary conditions. The partial differential equation together with the free boundary conditions was changed into a singular nonlinear ordinary differential equation of two-point boundary value problem by using similarity transformation. An efficient approximate technique of the problem was presented. The approximate solution of the problem was obtained and the corresponding transfer behavior is discussed. The numerical results show that the approximate solution proposed in this paper can be successfully used to reveal the physical nature of the studied problem.

Data Availability

All the data in this paper are correct and reliable. They can be obtained by using MATLAB.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Both authors read and approved the final paper.

Acknowledgments

The research was supported by a grant from the National Natural Science Foundation of China (No. 11501496), the Natural Science Basic Research Plan in Shaanxi Province of China (No. 2014JQ2-1003), and the Doctor Start-up Research Fund of Yulin University (No. 13GK04).

References

P. von Böckh and T. Wetzel, Heat Transfer: Basics and Practice, Springer, Berlin, Germany, 2012.

Y. Q. Ren, A. J. Cheng, and H. Wang, “A uniformly optimal-order estimate for finite volume method for advection-diffusion equations,” Numerical Methods for Partial Differential Equations, vol. 30, no. 1, pp. 17–43, 2014.

K. Yapici and S. Obut, “Benchmark results for natural and mixed convection heat transfer in a cavity,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 25, no. 5, pp. 998–1029, 2015.

A. Muhammad, A. J. Chamkha, S. Iqbal, and M. Ahmad, “Effects of temperature-dependent viscosity and thermal conductivity on mixed convection flow along a magnetized vertical surface,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 26, no. 5, pp. 1580–1592, 2016.

I. Pankratova and A. Piatnitski, “On the behaviour at infinity of solutions to stationary convection-diffusion equations in a cylinder,” Discrete and Continuous Dynamical Systems - Series B, vol. 11, no. 4, pp. 935–970, 2009.

I. Pankratova and A. Piatnitski, “Homogenization of convection-diffusion equation in infinite cylinder,” Networks and Heterogeneous Media, vol. 6, no. 1, pp. 111–126, 2011.

I. Pettersson and A. Piatnitski, “Stationary convection-diffusion equation in an infinite cylinder,” Journal of Differential Equations, vol. 264, no. 7, pp. 4456–4487, 2018.

M. D. Surnachev, “On the uniqueness of solutions to stationary convection-diffusion equations with generalized divergence-free drift,” Complex Variables and Elliptic Equations. An International Journal, vol. 63, no. 7-8, pp. 1168–1184, 2018.

K. N. Soltanov, “Existence and nonexistence of the global solutions some nonlinear elliptic parabolic equations,” Differential Equations, vol. 29, no. 4, pp. 550–563, 1993.

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.

M. H. Guan, L. C. Zheng, and X. X. Zhang, “The similarity solution to a generalized diffusion equation with convection,” Advances in Dynamical Systems and Applications, vol. 1, no. 2, pp. 183–189, 2006.

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.

G. T. Mekuria and J. A. Rao, “Adaptive finite element method for steady convection-diffusion equation,” American Journal of Computational Mathematics, vol. 6, no. 3, pp. 275–285, 2016.

W. Qiu and K. Shi, “An HDG method for convection diffusion equation,” Journal of Scientific Computing, vol. 66, no. 1, pp. 346–357, 2016.

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.

X. F. He and K. Wang, “New finite difference methods for singularly perturbed convection-diffusion equations,” Taiwannese Journal of Mathematics, vol. 22, no. 4, pp. 949–978, 2018.

X. Piao, H.-J. Choi, S. D. Kim et al., “A fast singly diagonally implicit Runge-Kutta method for solving 1D unsteady convection-diffusion equations,” Numerical Methods for Partial Differential Equations, vol. 30, no. 7, pp. 788–812, 2013.

S. Ali Waleeda, “Approximation solution for non-linear volterra integral equation by using Ado-mian decomposition method,” Al-Mustansiriyah Journal of Science, vol. 21, no. 2, pp. 153–159, 2010.

S. S. Ray and R. K. Bera, “An approximate solution of a nonlinear fractional differential equation by adomian decomposition method,” Applied Mathematics and Computation, vol. 167, no. 1, pp. 561–571, 2005.

G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, Kluwer Academic Publishers, Boston, Mass, USA, 1994.

G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Aanlysis and Applications, vol. 135, no. 2, pp. 501–544, 1988.

G. Adomian, “On the solution of partial differential equations with specified boundary conditions,” Mathematical and Computer Modelling, vol. 140, no. 2, pp. 569–581, 1994.

J. Sunday, “Convergence analysis and implementation of adomian decomposition method on second-order oscillatory problems,” Asian Research Journal of Mathematics, vol. 2, no. 5, pp. 1–12, 2017.

A. Abdelrazec and D. Pelinovsky, Convergence of ADM for Initial Value Problems, Wiley Periodicals, 2009.

M. M. Hosseini and H. Nasabzadeh, “On the convergence of Adomian decomposition method,” Applied Mathematics and Computation, vol. 182, no. 1, pp. 536–543, 2006.

Y. Cherruault, G. Saccomandi, and B. Some, “New results for convergence of Adomian's method applied to integral equations,” Mathematical and Computer Modelling, vol. 16, no. 2, pp. 85–93, 1992.

G. Adomian, Nonlinear Stochastic Systems Theory and Applications to Physics, Academic Press, Dordrecht, Netherlands, 1989.

C. Q. Dai, Y. Y. Wang, and J. F. Zhang, “Analytical spatiotemporal localizations for the generalized (3+1)-dimensional nonlinear Schrödinger equation,” Optics Letters, vol. 35, no. 9, pp. 1437–1439, 2010.

C.-Q. Dai and J.-F. Zhang, “Controllable dynamical behaviors for spatiotemporal bright solitons on continuous wave background,” Nonlinear Dynamics, vol. 73, no. 3, pp. 2049–2057, 2013.

Y.-Y. Wang, C.-Q. Dai, G.-Q. Zhou, Y. Fan, and L. Chen, “Rogue wave and combined breather with repeatedly excited behaviors in the dispersion/diffraction decreasing medium,” Nonlinear Dynamics, vol. 87, no. 1, pp. 67–73, 2017.

C.-Q. Dai, J. Liu, Y. Fan, and D.-G. Yu, “Two-dimensional localized Peregrine solution and breather excited in a variable-coefficient nonlinear Schrodinger equation with partial nonlocality,” Nonlinear Dynamics, vol. 88, no. 2, pp. 1373–1383, 2017.