Abstract

A newly computational method based on the Coiflet wavelet and homotopy analysis method is developed, which inherits the great nonlinear treatment of the homotopy analysis technique and the local high-precision capability of the wavelet approach, to give solutions to the classic problem of channel flow with moving walls. The basic principle of this suggested technique and the specific solving process are presented in detail. Its validity and efficiency are then checked via rigid comparisons with other computational approaches. It is found that the homotopy-based convergence-control parameter and the wavelet-based resolution level of Coiflet are two effective ways to improve on accuracies of solutions.

1. Introduction

Liquid flow in porous channels has received extensive attention owing to its wide industrial and engineering applications. The earliest study on channel flow with porous walls was performed by Berman [1], who obtained a perturbation solution for small Reynolds number. Morduchow [2] formulated a closed-form solution of Berman’s problem [1] over the full injection range by means of the method of averages. Berman’s solution was then extended by Proudman [3] to the large Reynolds number case and further generalized by Wah [4] for both small and large Reynolds number cases. Other extended studies have been done by various researchers on different aspects such as the theoretical solution at high suction Reynolds number [5], experimental observation for injection and suction effects [6], exact solution at high suction or injection velocity [7], complete solution for large injection with the viscous layer effect [8], and high-precision series solution for the uniformly porous tube [9]. A significant extension was given by Brady and Acrivos [10] for inclusion of the accelerating walls into the classic channel flow problem. Watson et al. [11] further investigated the case that the channel walls undergo both acceleration and porosity. Dauenhauer and Majdalani [12] examined the unsteady channel flow caused by uniformly time-dependent wall expansion ratio in the presence of porous and moving walls. Their work was further extended by Xu et al. [13] to the nonuniform time-dependent wall expansion ratio case.

In fact, most of fluid flow problems are nonlinear problems, which bring great difficulties to give their solutions. Various analytical or numerical methods have been proposed to solve these nonlinear flow problems such as the Runge–Kutta method [14, 15], the Akbari–Ganji method (AGM) [16, 17], the homotopy perturbation method (HPM), the collocation method (CM), the finite element method (FEM) [18], the Galerkin finite element method (GFEM) [19], and other numerical methods [20, 21]. Homotopy method is widely used because of its superiority in solving nonlinear problems. At the same time, many homotopy-based methods have also been proposed. Liu studied a parameter identification problem of two-phase porous media fractional flow equations [22] and then solved a nonlinear convection-diffusion equation [23] by a wavelet multiscale-homotopy method. In addition, based on the homotopy method, a multigrid-homotopy method was used to solve a nonlinear inverse problem governed by a convection-diffusion equation [24], and the numerical identification of diffusion parameters in a nonlinear convection-diffusion equation was investigated by an adaptive homotopy perturbation method [25].

Recently, the homotopy analysis method (HAM) [26] was developed for handling strong nonlinear problems, which possesses several advantages such as freedom to choose solution expressions, independence of small or large parameters, and adjustment of solution convergence. It has been used to solve many kinds of nonlinear problems in different fields [27, 28]. However, solutions by this technique are purely based on analytical computation; as a result, it is very time consuming when nonlinear problems are very complicated. On the contrary, the Coiflet wavelet technique [29] received applausive attention owing to its excellent features including feasible wavelet transformation, high computational efficiency, and superior expression of frequency and time information. This technique was initially used in signal analysis [30] and recently extended by Zhou et al. [31, 32] and Liu et al. [33, 34] to some classic nonlinear homogeneous problems in applied mathematics and mechanics. Very recently, the homotopy-based Coiflet wavelet method was developed by Yang and Liao [35, 36] for capturing solutions with high precision to nonlinear equations. One-dimensional and two-dimensional Bratus equations subjected to homogeneous boundary conditions were considered in their work as illustrative examples. Their proposed technique was then extended by Yu and his collaborators to solve the nonlinear plate bending problem [37] and cavity flow problem [38] with nonhomogeneous boundary conditions.

Though several successful applications of the homotopy-Coiflet wavelet technique have been achieved, many aspects of this approach are not yet clear. Further efforts are still needed to improve on its computational capability and efficiency. The aim of this paper is to extend the homotopy-Coiflet wavelet method to solve the classic porous channel flow problem governed by a fourth-order strong nonlinear system subjected to nonhomogeneous boundary conditions. The basic idea of this approach will be exhibited in detail. A range of wall expansion ratios for both injection and suction cases, especially for large injection and suction, will be taken into account. The accuracy and efficiency of our results will be then checked via a rigid comparison with previous studies.

2. Mathematical Description

As shown in Figure 1, liquid flow is driven through a two-dimensional rectangular channel bounded by a porous wall at . The channel walls undergo uniform expansion or contraction. The constant ratio of fluid injection and suction through the two opposing channel walls is applied. The fixed Cartesian coordinates in the space are used, with the -axis extending along the length of the channel and the -axis in the wall-normal direction.

The fourth-order nonlinear system describing fluid flow in this porous channel with moving walls is written as (refer to Dauenhauer and Majdalani [12])subject to the boundary conditionswhere prime denotes the derivation with , is the dimensionless stream function and the function of , is the wall expansion and extraction ratio, and is the Reynolds number with positive representing the injection case and vice versa. Note that this equation is derived directly from the Navier–Stokes equation via scaling and vorticity transformations without reduction. As a result, its solution can be deemed as the exact one of this famous equation.

2.1. Linearization

Here, the homotopy analysis method is employed to linearize nonlinear equation (1) and its boundary conditions (2). In doing so, we construct the following zeroth-order deformation equation:subject to the boundary conditionswhere is an embedding parameter evolving from 0 to 1, is the auxiliary linear operator, is the initial guess that satisfies all given boundary conditions, is the convergence-control parameter, is the mapping function of , and is the nonlinear operator, which is defined by

Since the computational domain is , the highest-order derivative of equation (1) is fourth-order, and the quantity of boundary conditions is four; we choose as the auxiliary linear operator for this problem.

From equation (3), it is readily known and . To get the exact solution , we expand the solution of equation (3) bywhere is the high-order subfunction. Clearly, is acquired at .

Note that the initial guess should satisfy the given boundary conditions. According to equation (4), the initial guess is chosen as

To determine , we construct the high-order deformation equation, by differentiating equation (3) times with respect to , then dividing by , and finally setting , aswherein which denotes a group of subfunctions .

The corresponding high-order boundary conditions then take the form

2.2. Coiflet Wavelet Expression

Different from the pure HAM technique, we select the Coiflet wavelet function rather than the power series given by Xu et al. [13] as the solution expression of .

Following Wang [39] and Liu [40], an arbitrary function is approximated by Coiflet wavelet aswhere is the vanishing moment, is the first-order vanishing moment, and is the wavelet basis.

In practical applications of wavelet approximation in boundary value problems, boundary modification is needed for handling boundary conditions. The modified boundary technique [32, 34, 39] for Coiflet wavelet is constructed asin which

Note that and are coefficients determined by the following coefficient matrices and P1via the relation and with . Taking the boundary conditions (4) into consideration, it is known that satisfies the inhomogeneous Dirichlet and homogeneous Neumann types of boundary conditions. Hence, following Yu et al. [37, 38, 4143], by modifying the coefficient matrices, the Coiflet wavelet basis is selected as

Furthermore, and can be expressed by Coiflet wavelet aswhere

Note that the boundary conditions on both sides of the Coiflet wavelet expansion of are rectified due to the boundary condition constraints at and .

Generally, from Coiflet wavelet approximation (16), the -order derivatives with respect to are obtained by the following formula as

Therefore, th-order deformation equation (8) is expressed as

Following Yang and Liao [35, 36], the Galerkin method is employed to solve -order deformation equation (19). Multiplying both sides of equation (19) by and then integrating on the domain , we obtain the th-order wavelet-Galerkin deformation equation aswhereand the connection coefficient matrices are defined asin which and are the generalized connection coefficients defined by

Note that the generalized connection coefficients change as the auxiliary linear operator and boundary conditions vary.

By solving the th-order wavelet-Galerlin equations from the order , the wavelet homotopy solution can be completely determined.

3. Validation

3.1. Validation in the PDE Case

In this section, we will present a 2D linear problem to verify the accuracy and efficiency of our proposed method. The governing equation considered here is as follows:where

The corresponding homogeneous Neumann boundary conditions are

Note that equation (24) has the analytical solution

In our approach, by modifying the coefficient matrices, the Coiflet wavelet basis is expressed by

Furthermore, can be expressed by Coiflet wavelet as

Therefore, equation (24) is expressed as

As described above, multiplying both sides of equation (30) by and then applying the multiple integrals on the domain , we obtain the wavelet-Galerkin equation aswhere vec() is a straighten vector operator which pulls a matrix into a vector row by row and is a matrix defined bywhere the coefficient matrices are obtained through

Here, denotes the operator of the Kronecker tensor product, are matrices, and the generalized connection coefficients and are defined as equation (23).

To illustrate the error level of our solution, we define the average square error based on the exact solution (27) as follows:

The solution of equation (24) is shown in Figure 2, whose residual error is listed in Table 1. It is found that the average square error decreases rapidly with the resolution level increasing. For , especially, the error reaches to within 1 minute CPU time. This completely demonstrates the accuracy and great efficiency of our proposed method.

3.2. Convergence and Error Analysis in the Nonlinear Case

To further verify the accuracy of our proposed method on nonlinear problems, we will calculate a special case of the fourth-order nonlinear system describing fluid flow in this porous channel and compare it with other numerical results.

Firstly, to test the convergence of our solution, we define the following residual error function:where and denote the -order and the -order wavelet solutions, respectively.

Then, we used the numerical solution by means of the software package BVP4c in software Matlab as a reference solution since the original problem did not have an exact analytical solution. Therefore, we define the relative error as follows:

Based on the definitions of Dauenhauer and Majdalani [12], the axial velocity and the normal velocity are given by

The residual error functions of the axial and normal velocities are therefore expressed as

We start our verification for , corresponding to the case of fluid flow through the porous channel with fixed walls. The injection Reynolds number is chosen as , and the convergence-control parameters for and are designated as and , respectively. As shown in Figure 3, the residual errors for both and reduce very rapidly as the computational order evolves. The least residual error reaches about at the computational order when the resolution level . This figure clearly shows that error reduction can be achieved either by the increment of the resolution level or the increase of the computational order.

We then check results of the axial velocity and the normal velocity for different prescribed points at certain computational orders. As shown in Table 2 and Table 3, the computational results at each discrete point approach to the fixed values very quickly. Furthermore, we compare these results with the numerical ones, and the very excellent agreement is found, as shown in Figure 4. The corresponding computational time and the relative error for different resolution levels are listed in Table 4. It is seen that the computational efficiency is high enough to give highly accurate solutions in one second, and our wavelet solution is very close to the numerical result, and the relative error can reach to . Note that the wavelet-homotopy computational parameters used here are for and for . The computational order is .

4. Results and Discussion

As discussed above, investigations of the special case of the fourth-order nonlinear system show excellent computational efficiency of our proposed method. Therefore, we apply it to the common cases of the fourth-order nonlinear system for a range of wall expansion ratios with different injection or suction Reynolds number . Without loss of generality, the resolution level is used for all the following calculations, but the convergence-control parameter needs to be reasonably selected based on the nonlinear strength of the problem at different and .

It has been known that the traditional homotopy analysis method is very powerful to give highly accurate solutions of nonlinear problems. Here, we compare the computational capability between the traditional HAM technique and our proposed method, firstly. As moderate physical parameters are prescribed, both approaches give very excellent results. However, our proposed technique is superior to the traditional HAM technique for large physical parameters such as and , as shown in Figures 5 and 6. In these special cases, the homotopy Padé technique [13] has to be used to optimize the HAM results for high-precision solutions.

Tables 5 and 6 give the comparison of CPU times for the maximum values of the axial velocity by the wavelet-homotopy method and HAM technique with different for and , respectively. Note that the residual errors used here are just an indication to illustrate solution validity. It is seen that our proposed technique has superior computational efficiency over the traditional HAM technique; especially, nonlinearity becomes very strong. For example, in the case of and , the HAM technique requires 28,450.3 seconds to complete the calculation, but the wavelet-homotopy technique only takes 0.6066 second. Note that the resolution level is employed in our computation with the convergent-control parameters being listed in Tables 7 and 8.

Based on the above analysis, we further obtain the solutions for various physical parameters readily by means of our proposed technique, as shown in Figures 7 and 8. It can be seen that the values of the axial velocity reach the highest at the middle section of the porous channel which is , and as the wall expansion ratios, , increase, the maximum value of the axial velocity becomes gradually larger; however, this situation is reversed near the walls. In addition, increasing wall expansion causes more rapid change in the value of the axial velocity from the walls to the middle section of the porous channel, while increasing wall contraction results in a smoother axial velocity curve with a wider range of plateaus. What is noteworthy is that flow reversal happens near the walls at , and 2.5, and , cases, which is the result of the rapid expansion of the wall creating a suction area near the wall, so the fluid is forced to move in the direction of the wall to occupy the space created by the wall expansion process. This is also consistent with the findings by Dauenhauer and Majdalani [12].

5. Conclusions

In this paper, the wavelet-homotopy method has been developed to obtain highly accurate solutions of fourth-order nonlinear boundary value problems describing fluid flows in the porous channel with moving walls. By far, now, homotopy-wavelet technique has been investigated for this kind of problems. The full solution procedure has been given in detail. We have performed wavelet approximation on the homogeneous Newman boundary condition by modifying the Coiflet wavelet basis function. The efficiency and accuracy of the Coiflet wavelet method have been verified by an example of a linear two-dimensional problem. Furthermore, the application of our homotopy-based Coiflet wavelet method to nonlinear problems has been carried out. By comparing with the results of the traditional homotopy analysis method, the reliability and efficiency of our Coiflet wavelet-homotopy method in solving strong nonlinear problems are further explained. The solutions of classic problems with liquid flow in porous channels for various physical parameters have been obtained readily by our proposed method. It is expected that this technique is expected to give accurate solutions for more complicated problems with strong nonlinearity.

Nomenclature

:Wall expansion and extraction ratio
:Generalized connection coefficient
:Connection coefficient matrix in the wavelet-Galerkin deformation equation
:Auxiliary linear operator
:Nonlinear operator
:Mapping function of in HAM
:Coiflet wavelet basis
:Convergence-control parameter
:Relative error
:Average square error
:Dimensionless stream function
:Initial guess
:Coiflet wavelet basis under Neumann types of boundary conditions
:First-order vanishing moment of the Coiflet wavelet
:Vanishing moment of the Coiflet wavelet
:Coefficient in modified Coiflet wavelet basis
:Coefficient matrix in modified Coiflet wavelet basis
:Embedding parameter in HAM
:Reynolds number
:Residual error
:Axial velocity and normal velocity.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was partially supported by the National Natural Science Foundation of China (Grant no. 11872241).