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.     

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.     

Einige laminare Grenzschichtströmungen,berechnet mittels einer neuen Reihenmethode

Summary A new and formally exact series solution of the general laminar boundary layer problem for arbitrary outer pressure distribution [1] is discussed briefly. The main part of the paper is concerned with the application of this solution to the numerical evaluation of instructive examples of boundary layer flows. Application of this new and rigorous series method requires only a short and simple numerical calculation, involving the application of tables of certain universal functions.      

Numerical solution of singularly perturbed convection-diffusion problem using parameter uniform B-spline collocation method

This paper is concerned with a numerical scheme to solve a singularly perturbed convection-diffusion problem. The solution of this problem exhibits the boundary layer on the right-hand side of the domain due to the presence of singular perturbation parameter ε. The scheme involves B-spline collocation method and appropriate piecewise-uniform Shishkin mesh. Bounds are established for the derivative of the analytical solution. Moreover, the present method is boundary layer resolving as well as second-order uniformly convergent in the maximum norm. A comprehensive analysis has been given to prove the uniform convergence with respect to singular perturbation parameter. Several numerical examples are also given to demonstrate the efficiency of B-spline collocation method and to validate the theoretical aspects.     

Falkner-Skan equation for flow past a moving wedge with suction or injection

The characteristics of steady two-dimensional laminar boundary layer flow of a viscous and incompressible fluid past a moving wedge with suction or injection are theoretically investigated. The transformed boundary layer equations are solved numerically using an implicit finite-difference scheme known as the Keller-box method. The effects of Falkner-Skan power-law parameter (m), suction/injection parameter (f0) and the ratio of free stream velocity to boundary velocity parameter (λ) are discussed in detail. The numerical results for velocity distribution and skin friction coefficient are given for several values of these parameters. Comparisons with the existing results obtained by other researchers under certain conditions are made. The critical values off ₀,m and λ are obtained numerically and their significance on the skin friction and velocity profiles is discussed. The numerical evidence would seem to indicate the onset of reverse flow as it has been found by Riley and Weidman in 1989 for the Falkner-Skan equation for flow past an impermeable stretching boundary.     

An Asymptotic-Numerical Method for a Class of Weakly Coupled System of Singularly Perturbed Convection-Diffusion Equations

A weakly coupled convection dominated system of m-equations is analyzed. A higher order accurate asymptotic-numerical method is presented. The solutions of convection dominated problem are known to exhibit multi-scale character. There exist narrow region across the boundary of the domain where the solution exhibit steep gradient. This region is termed as boundary layer region and the solution of problem is said to have a boundary layer. Outside of this region, the solution of system behaves smoothly. To capture this multi-scale nature given system is factorized into two explicit systems. The degenerate system of initial value problems (IVPs), obtained by setting ??=?0, corresponds to the smooth solution, which lies outside of boundary layers. For solution inside boundary layers, a system of boundary value problems (BVPs) is obtained using stretching transformation. Regardless of this simple factorization, solutions of these systems preserve the key features of the given coupled system. Runge–Kutta method is used to solve the degenerate system of IVPs, whereas the system of BVPs is solved analytically. Stability and consistency of the proposed method is established. A uniform convergence of higher order is obtained. Possible extension to differential difference equations are also brought to attention. A comparative study of the present method with some state of art existing numerical schemes is carried out by means of several test problems. The results so obtained demonstrate the effectiveness and potential of present approach.     

The recovery of even polynomial potentials

We consider the problem of reconstructing an even polynomial potential from one set of spectral data of a Sturm-Liouville problem. We show that we can recover an even polynomial of degree 2m from m+1 given Taylor coefficients of the characteristic function whose zeros are the eigenvalues of one spectrum. The idea here is to represent the solution as a power series and identify the unknown coefficients from the characteristic function. We then compute these coefficients by solving a nonlinear algebraic system, and provide numerical examples at the end. Because of its algebraic nature, the method applies also to non self-adjoint problems.     

Error analysis of the DtN-FEM for the scattering problem in acoustics via Fourier analysis

In this paper, we are concerned with the error analysis for the finite element solution of the two-dimensional exterior Neumann boundary value problem in acoustics. In particular, we establish explicit priori error estimates in H¹ and L²- norms including both the effect of the truncation of the DtN mapping and that of the numerical discretization. To apply the finite element method (FEM) to the exterior problem, the original boundary value problem is reduced to an equivalent nonlocal boundary value problem via a Dirichlet-to-Neumann (DtN) mapping represented in terms of the Fourier expansion series. We discuss essential features of the corresponding variational equation and its modification due to the truncation of the DtN mapping in appropriate function spaces. Numerical tests are presented to validate our theoretical results.     

Numerical Solution of the Incompressible Navier-Stokes Equations with Singular Perturbation Considerations

A numerical method for coarse grids is proposed for the numerical solution of the incompressible Navier-Stokes equations. From singular perturbation considerations, we obtain partial differential equations and boundary conditions for the outer solution and the boundary layer correction. The former problem is solved with the finite difference method and the latter with the approximate method. Numerical experiments show that accurate outer flow and boundary flux result with little computational effort.     

Stability of boundary layer and rarefaction wave to an outflow problem for compressible Navier-Stokes equations under large perturbation

In this paper, we investigate the large-time behavior of solutions to an outflow problem for compressible Navier-Stokes equations. In 2003, Kawashima, Nishibata and Zhu [S. Kawashima, S. Nishibata, P. Zhu, Asymptotic stability of the stationary solution to the compressible Navier-Stokes equations in the half space, Comm. Math. Phys. 240 (2003) 483-500] showed there exists a boundary layer (i.e., stationary solution) to the outflow problem and the boundary layer is nonlinearly stable under small initial perturbation. In the present paper, we show that not only the boundary layer above but also the superposition of a boundary layer and a rarefaction wave are stable under large initial perturbation. The proofs are given by an elementary energy method.     

Struktur der Grenzschicht zwischen einem feldfreien Plasma und einem Magnetfeld im Vakuum

Self-consistent solutions of Vlasov's equations describing a plane boundary layer without shear are given. The currents within the layer are due to free and trapped electrons. The ion current is neglected sincem _e ?m _i . The free electrons and the ions are assumed Maxwellian at infinity. For the trapped electrons a one parameter family of distribution functions is assumed. Profiles are computed for potentials, fields, densities and currents for different numbers of trapped electrons. Numerical solutions are obtained by applying the method of finite elements to a variational formulation of the problem.     

On the Parameter-Uniform Convergence of Exponential Spline Interpolation in the Presence of a Boundary Layer

The paper is concerned with the problem of generalized spline interpolation of functions having large-gradient regions. Splines of the class C², represented on each interval of the grid by the sum of a second-degree polynomial and a boundary layer function, are considered. The existence and uniqueness of the interpolation L-spline are proven, and asymptotically exact two-sided error estimates for the class of functions with an exponential boundary layer are obtained. It is established that the cubic and parabolic interpolation splines are limiting for the solution of the given problem. The results of numerical experiments are presented.     

Numerical study of forced and free convective boundary layer flow of a magnetic fluid over a flat plate under the action of a localized magnetic field

The two-dimensional, steady, laminar, forced and free convective boundary layer flow of a magnetic fluid over a semi-infinite vertical plate, under the action of a localized magnetic field, is numerically studied. The magnetic fluid is considered to be water-based with temperature dependent viscosity and thermal conductivity. The study of the boundary layer is separated into two cases. In case I the boundary layer is studied near the leading edge, where it is dominated by the large viscous forces, whereas in case II the boundary layer is studied far from the leading edge of the plate where the effects of buoyancy forces increase. The numerical solution, for these two different cases, is obtained by an efficient numerical technique based on the common finite difference method. Numerical calculations are carried out for the value of Prandl number Pr =  49.832 (water-based magnetic fluid) and for different values of the dimensionless parameters entering into the problem and especially for the magnetic parameter Mn, the viscosity/temperature parameter Θ _r and the thermal/conductivity parameter S*. The analysis of the obtained results show that the flow field is influenced by the application of the magnetic field as well as by the variation of the viscosity and the thermal conductivity of the fluid with temperature. It is hoped that they could be interesting for engineering applications.     

The iterative solutions of 2nth-order nonlinear multi-point boundary value problems

The aim of this paper is to investigate the existence of iterative solutions for a class of 2nth-order nonlinear multi-point boundary value problems. The multi-point boundary condition under consideration includes various commonly discussed boundary conditions, such as three- or four-point boundary condition, (n + 2)-point boundary condition and 2(n − m)-point boundary condition. The existence problem is based on the method of upper and lower solutions and its associated monotone iterative technique. A monotone iteration is developed so that the iterative sequence converges monotonically to a maximal solution or a minimal solution, depending on whether the initial iteration is an upper solution or a lower solution. Two examples are presented to illustrate the results.     

The harmonic dirichlet problem for a cracked domain with jump conditions on cracks

The harmonic problem in a cracked domain is studied in R ^m , m?>?2. The boundary of the domain is assumed to be nonsmooth, while cracks are smooth. The Dirichlet condition is specified on the boundary of the domain. Jumps of the unknown function and its normal derivative are specified on the cracks. Uniqueness and solvability results are obtained. The problem is reduced to the uniquely solvable integral equation, its solution is given explicitely in the form of a series. The estimates of the solution of the problem depending on the boundary data are obtained.     

A new method for exact integration of some perturbed stiff linear systems of oscillatory type

This article presents a new family of real functions with values within the ring of M(m,R) matrices, Φ-functions for perturbed linear systems and a numerical method adapted for integration of this type of problem. This method permits the system solution to be expressed as a series of Φ-functions. The coefficients of this series are obtained through recurrences in which the perturbation intervenes.The Φ-functions series method has the advantage of being exactly integrated in the perturbed problem. For this purpose an appropriate B matrix is selected and used to construct the operator described in this article, thus annihilating the disturbance terms, transforming the system into a homogenous second-order system, which is exactly integrated with the two first Φ-functions.The article ends with a detailed study of four perturbed systems which illustrate how the method is used in stiff problems or in highly oscillatory problems, contrasting its behaviour by studying its accuracy in comparison with other well-known codes.     

B-spline collocation for solution of two-point boundary value problems

A numerical method based on B-spline is developed to solve the general nonlinear two-point boundary value problems up to order 6. The standard formulation of sextic spline for the solution of boundary value problems leads to non-optimal approximations. In order to derive higher orders of accuracy, high order perturbations of the problem are generated and applied to construct the numerical algorithm. The error analysis and convergence properties of the method are studied via Green’s function approach. O(h⁶) global error estimates are obtained for numerical solution of these classes of problems. Numerical results are given to illustrate the efficiency of the proposed method. Results of numerical experiments verify the theoretical behavior of the orders of convergence.     

Laguerre polynomial approach for solving linear delay difference equations

This paper presents a numerical method for the approximate solution of mth-order linear delay difference equations with variable coefficients under the mixed conditions in terms of Laguerre polynomials. The aim of this article is to present an efficient numerical procedure for solving mth-order linear delay difference equations with variable coefficients. Our method depends mainly on a Laguerre series expansion approach. This method transforms linear delay difference equations and the given conditions into matrix equation which corresponds to a system of linear algebraic equation. The reliability and efficiency of the proposed scheme are demonstrated by some numerical experiments and performed on the computer algebraic system Maple.     

Multipoint approach to the two-point boundary value problem

This paper treats nonlinear, two-point boundary value problems of the form $\text{x}\text{?}\text{? ?(x, t) = 0}$, in which the Jacobian matrix ?_x(x, t) is characterized by large positive eigenvalues. The resulting numerical difficulties are reduced by treating the two-point boundary value problem as a multipoint boundary value problem. A totally finite-difference approach is employed, thus bypassing the integration of the nonlinear equations, which characterizes shooting methods.The approach employed consists of extending to multipoint boundary value problems the modified-quasilinearization method developed by Miele and lyer for two-point boundary value problems. Basic to the method is the consideration of the performance index P, which measures the cumulative error in the differential equations, the boundary conditions, and the interface conditions.A modified-quasilinearization algorithm is generated by requiring the first variation of the performance index δP to be negative. This algorithm differs from the ordinary-quasilinearization algorithm because of the inclusion of the scaling factor or stepsize α in the system of variations. The main property of the modified-quasilinearization algorithm is the descent property: if the stepsize α is sufficiently small, the reduction in P is guaranteed. Convergence to the desired solution is achieved when the inequality P ? ? is met, where ? is a small, preselected number.The variations per unit stepsize $\text{Δx(t)}\text{α}\text{= A(t)}$ are governed by a system of mn nonhomogeneous, linear differential equations subjected to p initial conditions, q final conditions, and (m ? 1)n interface conditions, with p + q = n, where n is the dimension of the vector x and m is the number of subintervals. Therefore, the total number of boundary conditions and interface conditions is mn. The above system is solved employing the method of particular solutions: m(n + 1) particular solutions are combined linearly, and the coefficients of the combination are determined so that the linear system is satisfied.Two numerical examples are presented, one dealing with a linear system and one dealing with a nonlinear system. The examples illustrate the effectiveness as well as the rapidity of convergence of the present method.