A novel method for a class of parameterized singularly perturbed boundary value problems

In this paper a novel approach is presented for solving parameterized singularly perturbed two-point boundary value problems with a boundary layer. By the boundary layer correction technique, the original problem is converted into two non-singularly perturbed problems which can be solved using traditional numerical methods, such as Runge–Kutta methods. Several non-linear problems are solved to demonstrate the applicability of the method. Numerical experiments indicate the high accuracy and the efficiency of the new method.     

Analytical solutions of the boundary layer flow of power-law fluid over a power-law stretching surface

This article discusses analytical solutions for a nonlinear problem arising in the boundary layer flow of power-law fluid over a power-law stretching surface. Using perturbation method analytical solution is presented for linear stretching surface. This solution covers large range of shear thinning and shear thickening fluids and matches excellently with the numerical solution. Furthermore, some new exact solutions are found for particular combination of m (power-law stretching index) and n (power-law fluid index). This leads to generalize the case of linear stretching to nonlinear stretching surface. The effects of fluid index n on the boundary layer thickness and the skin friction for nonlinear stretching surface is analyzed and discussed. It is observed that the boundary layer thickness and the skin friction coefficient increase as non-linear parameter n decreases. This study gives a new dimension to obtain analytical solutions asymptotically for highly nonlinear problems which to the best of our knowledge has not been examined so far.     

An efficient method for solving non-linear singularly perturbed two points boundary-value problems of fractional order

In this paper, we discuss a numerical solution of a class of non-linear fractional singularly perturbed two points boundary-value problem. The method of solution consists of solving reduced problem and boundary layer correction problem. A series method is used to solve the boundary layer correction problem, and then the series solutions is approximated by the Pade’ approximant of order [m, m]. Some theoretical results are established and proved. Two numerical examples are discussed to illustrate the efficiency of the present scheme.     

Heat transfer in a Walter’s liquid B fluid over an impermeable stretching sheet with non-uniform heat source/sink and elastic deformation

In this present article an analysis is carried out to study the boundary layer flow behavior and heat transfer characteristics in Walter’s liquid B fluid flow. The stretching sheet is assumed to be impermeable, the effects of viscous dissipation, non-uniform heat source/sink in the presence and in the absence of elastic deformation (which was escaped from attention of researchers while formulating the viscoelastic boundary layer flow problems)on heat transfer are addressed. The basic boundary layer equations for momentum and heat transfer, which are non-linear partial differential equations, are converted into non-linear ordinary differential equations by means of similarity transformation. Analytical solutions are obtained for the resulting boundary value problems. The effects of viscous dissipation, Prandtl number, Eckert number and non-uniform heat source/sink on heat transfer (in the presence and in the absence of elastic deformation) are shown in several plots and discussed. Analytical expressions for the wall frictional drag coefficient, non-dimensional wall temperature gradient and non-dimensional wall temperature are obtained and are tabulated for various values of the governing parameters. The present study reveals that, the presence of work done by deformation in the energy equation yields an augment in the fluid’s temperature.     

Application of orthogonal collocation on finite elements for solving non-linear boundary value problems

A simple, convenient and easy approach to solve non-linear boundary value problems (BVP) using orthogonal collocation on finite elements (OCFE) is presented. The algorithm is the conjunction of finite element method (FEM) and orthogonal collocation method (OCM). The stability of the numerical results is checked by a novel algorithm which not only justifies the stability of the results but also checks the convergence of the method. The method is applied to the non-symmetric boundary value problems having Dirichlet’s and mixed Robbin’s boundary conditions.     

Heat and mass transfer in MHD non-Darcian flow of a micropolar fluid over a stretching sheet embedded in a porous media with non-uniform heat source and thermal radiation

A mathematical analysis has been carried out to study magnetohydrodynamic boundary layer flow, heat and mass transfer characteristic on steady two-dimensional flow of a micropolar fluid over a stretching sheet embedded in a non-Darcian porous medium with uniform magnetic field. Momentum boundary layer equation takes into account of transverse magnetic field whereas energy equation takes into account of Ohmic dissipation due to transverse magnetic field, thermal radiation and non-uniform source effects. An analysis has been performed for heating process namely the prescribed wall heat flux (PHF case). The governing system of partial differential equations is first transformed into a system of non-linear ordinary differential equations using similarity transformation. The transformed equations are non-linear coupled differential equations which are then linearized by quasi-linearization method and solved very efficiently by finite-difference method. Favorable comparisons with previously published work on various special cases of the problem are obtained. The effects of various physical parameters on velocity, temperature, concentration distributions are presented graphically and in tabular form.     

Numerical treatment of singularly perturbed two point boundary value problems exhibiting boundary layers

A boundary value method for solving a class of nonlinear singularly perturbed two point boundary value problems with a boundary layer at one end is proposed. Using singular perturbation analysis the method consists of solving two problems; namely, a reduced problem and a boundary layer correction problem. We use Pade’ approximation to obtain the solution of the latter problem and to satisfy the condition at infinity. Numerical examples will be given to illustrate the method.     

A note on the use of modified Adomian decomposition method for solving singular boundary value problems of higher-order ordinary differential equations

This paper extend the work [Yahya Qaid Hasan, Liu Ming Zhu. Solving singular boundary value problems of higher-order ordinary differential equations by modified Adomian decomposition method. Commun Nonlinear Sci Numer Simul. doi :10.1016/j.cnsns.2008.09.027] to high order of singular boundary value problems. Solution of these problems is considered by proposed modification of Adomian decomposition method. The proposed method can be applied to linear and nonlinear problems. Some examples are presented to show the ability of the method for linear and non-linear ordinary differential equation.     

B-spline solution of non-linear singular boundary value problems arising in physiology

We use B-spline functions to develop a numerical method for computing approximations to the solution of non-linear singular boundary value problems associated with physiology science. The original differential equation is modified at singular point then the boundary value problem is treated by using B-spline approximation. The numerical method is tested for its efficiency by considering three model problems from physiology.     

On the Asymptotic and Numerical Analyses of Exponentially III-Conditioned Singularly Perturbed Boundary Value Problems

Asymptotic and numerical methods are used to study several classes of singularly perturbed boundary value problems for which the underlying homogeneous operators have exponentially small eigenvalues. Examples considered include the familiar boundary layer resonance problems and some extensions and certain linearized equations associated with metastable internal layer motion. For the boundary layer resonance problems, a systematic projection method, motivated by the work of De Groen [1], is used to analytically calculate high-order asymptotic solutions. This method justifies and extends some previous results obtained from the variational method of Grasman and Matkowsky [2]. A numerical approach, based on an integral equation formulation, is used to accurately compute boundary layer resonance solutions and their associated exponentially small eigenvalues. For various examples, the numerical results are shown to compare very favorably with two-term asymptotic results. Finally, some Sturm-Liouville operators with exponentially small spectral gap widths are studied. One such problem is applied to analyzing metastable internal layer motion for a certain forced Burgers equation.     

A singular perturbation approach to solve Burgers–Huxley equation via monotone finite difference scheme on layer-adaptive mesh

This paper deals with a numerical method for solving one-dimensional unsteady Burgers–Huxley equation with the viscosity coefficient ε. The parameter ε takes any values from the half open interval (0, 1]. At small values of the parameter ε, an outflow boundary layer is produced in the neighborhood of right part of the lateral surface of the domain and the problem can be considered as a non-linear singularly perturbed problem with a singular perturbation parameter ε. Using singular perturbation analysis, asymptotic bounds for the derivatives of the solution are established by decomposing the solution into smooth and singular components. We construct a numerical scheme that comprises of implicit-Euler method to discretize in temporal direction on uniform mesh and a monotone hybrid finite difference operator to discretize the spatial variable with piecewise uniform Shishkin mesh. To obtain better accuracy, we use central finite difference scheme in the boundary layer region. Shishkin meshes are refined in the boundary layer region, therefore stability constraint is satisfied by proposed scheme. Quasilinearization process is used to tackle the non-linearity and it is shown that quasilinearization process converges quadratically. The method has been shown to be first order uniformly accurate in the temporal variable, and in the spatial direction it is first order parameter uniform convergent in the outside region of boundary layer, and almost second order parameter uniform convergent in the boundary layer region. Accuracy and uniform convergence of the proposed method is demonstrated by numerical examples and comparison of numerical results made with the other existing methods.     

Smooth boundary integration method in a linear two-dimensional problem of incompressible fluid and solid

This paper presents a new application of a theoretical and computational method of smooth boundary integration which belongs to the methods of boundary integral equations. Smooth integration is not a method of approximation. In its final analytical form, a smooth-kernel integral equation is computerized easily and accurately.

Smooth integration is associated with a “pressure-vorticity” formulation which covers linear problems in elasticity and fluid mechanics. The solution presented herein is essentially the same as that reported in an earlier paper for regular elasticity. The constraint of incompressibility does not cause difficulties in the pressure-vorticity formulation.

The linear fluid mechanics problem formulated and solved in this paper covers Stokes' problem of a slow viscous flow, and has a wider interpretation. The translational inertia forces are incorporated in the linear problem, as in Euler's dynamic theory of inviscid flow. The centrifugal inertia forces are left for the non-linear problem. The linear problem is perceived as a step in solution of the non-linear problems.     

Solutions of MHD boundary layer equations with transverse magnetic field and arbitrary pressure gradient

Certain solutions of Magnetohydrodynamic boundary layer equations for a flat plate with a transverse magnetic field fixed relative to the fluid have been obtained using a power-series method given by Hassan. This power series solution has all the desirable qualities of Görtler series and, in addition, its zeroth order term which is governed by a non-linear total differential equation can be given in closed form. The first order term governed by a linear total differential equation has been integrated numerically. The results are tabulated for various values of S, the interaction parameter. The results show a decrease in the boundary layer thickness with a consequent increase in skin friction as the strength of magnetic field is increased.     

Heat transfer in a liquid film over an unsteady stretching surface with viscous dissipation in presence of external magnetic field

This paper presents a mathematical analysis of MHD flow and heat transfer to a laminar liquid film from a horizontal stretching surface. The flow of a thin fluid film and subsequent heat transfer from the stretching surface is investigated with the aid of similarity transformation. The transformation enables to reduce the unsteady boundary layer equations to a system of non-linear ordinary differential equations. Numerical solution of resulting non-linear differential equations is found by using efficient shooting technique. Boundary layer thickness is explored numerically for some typical values of the unsteadiness parameter S and Prandtl number Pr, Eckert number Ec and Magnetic parameter Mn. Present analysis shows that the combined effect of magnetic field and viscous dissipation is to enhance the thermal boundary layer thickness.     

Spectral methods for some singularly perturbed third order ordinary differential equations

Spectral methods with interface point are presented to deal with some singularly perturbed third order boundary value problems of reaction-diffusion and convection-diffusion types. First, linear equations are considered and then non-linear equations. To solve non-linear equations, Newton’s method of quasi-linearization is applied. The problem is reduced to two systems of ordinary differential equations. And, then, each system is solved using spectral collocation methods. Our numerical experiments show that the proposed methods are produce highly accurate solutions in little computer time when compared with the other methods available in the literature.      

Boundary Integral Equation Solution of Non-linear Plane Potential Problems

This paper presents the boundary integral equation (BIE) formulation,and numerical solution procedure for two-dimensional problemsgoverned by Laplac's equation and subject to non-linear boundaryconditions. The introduction of non-linear terms constitutesa fundamental extension of the BIE method, as previous applicationshave been restricted entirely to linear problems. Furthermore,non-linearities necessitate the use of iterative solution techniqueswhich present the conceptual disadvantage that a solution isnot guaranteed. However, such difficulties were not encounteredwith the Newton—Raphson method employed in this study.The various features of the BIE technique are illustrated bythe application to a physical problem which is of significancein heat exchanger design.     

Method of successive approximations for three-dimensional laminar boundary layer problems (locally self-similar case) : PMM vol. 37, n6, 1973, pp. 974–983

We present an analytical method for the computation of problems of incompressible boundary layer theory based on an application of the method of successive approximations. The system of equations is reduced to a form suitable for integration. Parameters characterizing the external flow and the body geometry are contained only in the coefficients of the system and do not enter into the boundary conditions. The transformed momentum equations are integrated across the boundary layer from a current value to infinity with the boundary conditions taken into account. If the integration is made from zero to infinity, then the equations pass over into the Kármán relations. Integrating the system of equations a second time, using the boundary conditions at the wall, we obtain a system of nonlinear integro-differential equations. To solve this system of equations we apply the method of successive approximations. To satisfy the boundary Conditions at infinity we introduce, at each step of the iterations, unknown “governing” functions. From the conditions at the outer side of the boundary layer we obtain additional equations for their determination. With the iterational algorithm formulated in this way, the boundary conditions, both on the body and at the outer side of the boundary layer; are satisfied automatically.We consider a locally self-similar approximation. In this case, relative to the “governing” functions, we obtain an algebraic system of equations. We write out the solution in the first approximation. The results obtained in the first approximation are compared with the results of finite-difference computations for a wide range of problems. The results obtained in this paper are compared with those obtained in [1] for the flow in the neighborhood of a stagnation point. An indication is given of the nonuniqueness of the solutions of the three-dimensional boundary layer equations.     

合成展开法应用于球壳对称弯曲的边界层问题

本文推广钱伟长在[5]中提出的合成展开法分析双参数边界层问题. 对于受均布荷载作用的球壳对称变形问题,其非线性平衡方程可以写成(2.3a),(2.3b):式中ε与δ是待定参数.当δ=1,ε是小参数时,这是第一边界层问题:当δ与ε都县小参数时.这是第二边界层问题. 对于上述问题,我们假定ε,δ和p满足ε3pδ=1-ε在这个条件下,应用推广的合成展开法,求出上述问题具有固定边界条件情况的渐近解.     

Lattice Boltzmann simulation for high-speed compressible viscous flows with a boundary layer

In this study, the lattice Boltzmann method is employed for simulating high-speed compressible viscous flows with a boundary layer. The coupled double-distribution-function lattice Boltzmann method proposed by Li et al. (2007) is employed because of its good numerical stability and non-free-parameter feature. The non-uniform mesh construction near the wall boundary in fine grids is combined with an appropriate wall boundary treatment for the finite difference method in order to obtain accurate spatial resolution in the boundary layer problem. Three typical problems in high-speed viscous flows are solved in the lattice Boltzmann simulation, i.e., the compressible boundary layer problem, shock wave problem, and shock boundary layer interaction problem. In addition, in-depth comparisons are made with the non-oscillatory and non-free-parameter dissipation (NND) scheme and second order upwind scheme in the present lattice Boltzmann model. Our simulation results indicate the great potential of the lattice Boltzmann method for simulating high-speed compressible viscous flows with a boundary layer. Further research is needed (e.g., better numerical models and appropriate finite difference schemes) because the lattice Boltzmann method is still immature for high-speed compressible viscous flow applications.     

On the global existence and uniqueness of solutions to the nonstationary boundary layer system

In this paper, we study the problem of boundary layer for nonstationary flows of viscous incompressible fluids. There are some open problems in the field of boundary layer. The method used here is mainly based on a transformation which reduces the boundary layer system to an initial-boundary value problem for a single quasilinear parabolic equation. We prove the existence of weak solutions to the modified nonstationary boundary layer system. Moreover, the stability and uniqueness of weak solutions are discussed.