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.     

Bernstein Ritz‐Galerkin method for solving an initial‐boundary value problem that combines Neumann and integral condition for the wave equation

In this article, the Ritz‐Galerkin method in Bernstein polynomial basis is implemented to give an approximate solution of a hyperbolic partial differential equation with an integral condition. We will deal here with a type of nonlocal boundary value problem, that is, the solution of a hyperbolic partial differential equation with a nonlocal boundary specification. The nonlocal conditions arise mainly when the data on the boundary cannot be measured directly. The properties of Bernstein polynomial and Ritz‐Galerkin method are first presented, then Ritz‐Galerkin method is used to reduce the given hyperbolic partial differential equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique presented in this article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

流过半无限竖直可渗透壁面时的磁流体动力学自由对流

研究不可压缩粘性导电流体,流过半无限竖直可渗透平板时,将其偏微分形式的流动和传热的基本控制方程,应用适当的相似变换,简化为非线性的常微分方程组．对两种抽吸参数:大的和小的抽吸参数,采用摄动法得到变换后方程的近似解．数值结果表明,随着磁场参数和抽吸参数的增大,任意点的速度场在减小;磁场参数的影响,引起热边界层厚度的增大;速度和温度场随着热汇参数的增大而减小．     

Hopf-Cole transformation to some systems of partial differential equations

In this paper we study a system of nonlinear partial differential equations which we write as a Burgers equation for matrix and use the Hopf-Cole transformation to linearize it. Using this method we solve initial value problem and initial boundary value problems for some systems of parabolic partial differential equations. Also we study an initial value problem for a system of nonlinear partial differential equations of first order which does not have solution in the standard distribution sense and construct an explicit solution in the algebra of generalized functions of Colombeau. Received November 1999     

Group method analysis of melting effect on MHD mixed convection flow from radiate vertical plate embedded in a saturated porous media

The linear transformation group approach is developed to simulate problem of hydromagnetic heat transfer by mixed convection along vertical plate in a liquid saturated porous medium in the presence of melting and thermal radiation effects for opposing external flow. The application of a one-parameter transformation group reduces the number of independent variables by one so that the governing partial differential equations with the boundary conditions reduce to an ordinary differential equations with appropriate corresponding conditions. The Runge-Kutta shooting method is used to solve the determining equations of the set of nonlinear ordinary differential equations. are presented in the form of the temperature and flow fields in the melting region within the boundary layer for different parameters entering into the analysis. Also the effects of the pertinent parameters on the rate of the heat transfer in terms of the local Nusselt number at the solid–liquid interface are also discussed.     

Boundary element solution of the two-dimensional groundwater flow equation with stochastic free surface boundary condition

An innovative approach to the approximate solution of stochastic partial differential equations in groundwater flow is presented. The method uses a formulation of the Ito's lemma in Hilbert spaces to derive partial differential equations satisfying the moments of the solution process. Since the moments equations are deterministic, they could be solved by any analytical or numerical method existing in the literature. This permits the analysis and solution of stochastic partial differential equations occurring in two-dimensional or three-dimensional domains of any geometrical shape. The method is tested for the first time in the present paper through a practical application in a sandy phreatic aquifer at the Chalk River Nuclear Laboratories, Ontario, Canada. The equation solved is the two-dimensional LaPlace equation with a dynamic, randomly perturbed, free surface boundary condition. The moments equations are derived and solved by using the boundary integral equation method. A comparison is made with a previous analytical solution obtained by applying the randomly forced one-dimensional Boussinesq equation, and some observations on modeling procedures are given.     

MHD power-law fluid flow and heat transfer over a non-isothermal stretching sheet

This article presents a numerical solution for the magnetohydrodynamic (MHD) non-Newtonian power-law fluid flow over a semi-infinite non-isothermal stretching sheet with internal heat generation/absorption. The flow is caused by linear stretching of a sheet from an impermeable wall. Thermal conductivity is assumed to vary linearly with temperature. The governing partial differential equations of momentum and energy are converted into ordinary differential equations by using a classical similarity transformation along with appropriate boundary conditions. The intricate coupled non-linear boundary value problem has been solved by Keller box method. It is important to note that the momentum and thermal boundary layer thickness decrease with increase in the power-law index in presence/absence of variable thermal conductivity.     

Unsteady convective boundary layer flow of a viscous fluid at a vertical surface with variable fluid properties

In this paper we present numerical solutions to the unsteady convective boundary layer flow of a viscous fluid at a vertical stretching surface with variable transport properties and thermal radiation. Both assisting and opposing buoyant flow situations are considered. Using a similarity transformation, the governing time-dependent partial differential equations are first transformed into coupled, non-linear ordinary differential equations with variable coefficients. Numerical solutions to these equations subject to appropriate boundary conditions are obtained by a second order finite difference scheme known as the Keller-Box method. The numerical results thus obtained are analyzed for the effects of the pertinent parameters namely, the unsteady parameter, the free convection parameter, the suction/injection parameter, the Prandtl number, the thermal conductivity parameter and the thermal radiation parameter on the flow and heat transfer characteristics. It is worth mentioning that the momentum and thermal boundary layer thicknesses decrease with an increase in the unsteady parameter.     

利用Prüfer变换求解微分方程

针对一类二阶非线性常微分方程,利用Prüfer变换将其约化为特殊的一阶常微分方程组,从而使其求解过程得以简化.实例说明应用Prüfer变换求解一类偏微分方程边值问题的技巧.     

偶数阶非线性中立型偏微分方程系统的振动性

研究一类偶数阶非线性中立型偏泛函微分方程系统解的振动性,利用微分不等式方法和黎卡提变换,获得了该类系统在两类子同边值条件下振动的若干充分条件.     

Goursat problem for two-dimensional second-order hyperbolic operator-differential equations with variable domains

We develop a modification of the energy inequality method and use it to prove the well-posedness of the Goursat problem for linear second-order hyperbolic differential equations with operator coefficients whose domains depend on the two-dimensional time. An energy inequality for strong solutions of this Goursat problem is derived with the help of abstract smoothing operators, and we prove that the range of the problem is dense by using properties of a regularizing Cauchy problem whose inverse operator is a family of smoothing operators of a new type. We give an example of a well-posed boundary value problem for a two-dimensional complete second-order hyperbolic partial differential equation with Goursat conditions and with a boundary condition depending on the two-dimensional time.     

Numerical solution of two-dimensional reaction-diffusion Brusselator system

In this paper, polynomial based differential quadrature method (DQM) is applied for the numerical solution of a class of two-dimensional initial-boundary value problems governed by a non-linear system of partial differential equations. The system is known as the reaction-diffusion Brusselator system. The system arises in the modeling of certain chemical reaction-diffusion processes. In Brusselator system the reaction terms arise from the mathematical modeling of chemical systems such as in enzymatic reactions, and in plasma and laser physics in multiple coupling between modes. The numerical results reported for three specific problems. Convergence and stability of the method is also examined numerically.     

Three-layer finite difference method for solving linear differential algebraic systems of partial differential equations

A boundary value problem is examined for a linear differential algebraic system of partial differential equations with a special structure of the associate matrix pencil. The use of an appropriate transformation makes it possible to split such a system into a system of ordinary differential equations, a hyperbolic system, and a linear algebraic system. A three-layer finite difference method is applied to solve the resulting problem numerically. A theorem on the stability and the convergence of this method is proved, and some numerical results are presented.     

New computational method for solving some 2-dimensional nonlinear Volterra integro-differential equations

The aim of this paper is to present an efficient numerical procedure for solving the two-dimensional nonlinear Volterra integro-differential equations (2-DNVIDE) by two-dimensional differential transform method (2-DDTM). The technique that we used is the differential transform method, which is based on Taylor series expansion. Using the differential transform, 2-DNVIDE can be transformed to algebraic equations, and the resulting algebraic equations are called iterative equations. New theorems for the transformation of integrals and partial differential equations are introduced and proved. The reliability and efficiency of the proposed scheme are demonstrated by some numerical experiments.     

A Lyapunov function for a two-chemical species version of the chemotaxis model

Since the accuracy of finite element solutions of partial differential equations is generally mesh dependent, especially when solutions have singularities and discontinuities, a proper mesh generation is often important and sometimes crucial for an accurate numerical approximation of such problems. In this paper, the mesh transformation method is applied to the boundary value problems of elliptic partial differential equations, and it is proved that the method leads to the optimal finite element solutions. AMS subject classification (2000) 73C50, 65K10, 65N12, 65N30     

Heat transfer in MHD viscoelastic fluid flow over a stretching sheet with variable thermal conductivity,non-uniform heat source and radiation

An analysis has been carried out to study the magnetohydrodynamic boundary layer flow and heat transfer characteristics of a non-Newtonian viscoelastic fluid over a flat sheet with a linear velocity in the presence of thermal radiation and non-uniform heat source. The thermal conductivity is assumed to vary as a linear function of temperature. The basic equations governing the flow and heat transfer are in the form of partial differential equations, the same have been reduced to a set of non-linear ordinary differential equations by applying suitable similarity transformation. The transformed equations are solved analytically by regular perturbation method. Numerical solution of the problem is also obtained by the efficient shooting method, which agrees well with the analytical solution. The effects of various physical parameters such as viscoelastic parameter, Chandrasekhar number, Prandtl number, variable thermal conductivity parameter, Eckert number, thermal radiation parameter and non-uniform heat source/sink parameters which determine the temperature profiles are shown in several plots and the heat transfer coefficient is tabulated for a range of values of said parameters. Some important findings reported in this work reveals that combined effect of variable thermal conductivity, radiation and non-uniform heat source have significant impact in controlling the rate of heat transfer in the boundary layer region.     

Analytical eigenvalue bounds for matrices arising when discretizing the self-adjoint second order elliptic partial differential equations

This paper shows how to derive analytical expressions for the eigenvalue bounds of matrices arising when using a fast method for separable finite difference equations for the numerical solution of the first three boundary value problems for the two-dimensional self-adjoint second order elliptic partial differential equation in a rectangle.     

Numerical solution of the momentum and heat transfer equations for a hydromagnetic flow due to a stretching sheet of a non-uniform property micropolar liquid

A study of the hydromagnetic flow due to a stretching sheet and heat transfer in an incompressible micropolar liquid is made. Temperature-dependent thermal conductivity and a non-uniform heat source/sink render the problem analytically intractable and hence a numerical study is made using the shooting method based on Runge-Kutta and Newton-Raphson methods. The two problems of horizontal and vertical stretching are considered to implement the numerical method. The former problem involves one-way coupling between linear momentum and heat transport equations and the latter involves two-way coupling. Further, both the problems involve two-way coupling between the non-linear equations of conservation of linear and angular momentums. A similarity transformation arrived at for the problem using the Lie group method facilitates the reduction of coupled, non-linear partial differential equations into coupled, non-linear ordinary differential equations. The algorithm for solving the resulting coupled, two-point, non-linear boundary value problem is presented in great detail in the paper. Extensive computation on velocity and temperature profiles is presented for a wide range of values of the parameters, for prescribed surface temperature (PST) and prescribed heat flux (PHF) boundary conditions.     

Qualitative analysis of the invasion free boundary problem in biofilms

The work presents the qualitative analysis of the free boundary value problem related to the invasion model for multispecies biofilms. This model is based on the continuum approach for biofilm modeling and consists of a system of nonlinear hyperbolic partial differential equations for microbial species growth and spreading, a system of semilinear elliptic partial differential equations describing the substrate trends and a system of semilinear elliptic partial differential equations accounting for the diffusion and reaction of motile species within the biofilm. The free boundary evolution is regulated by a nonlinear ordinary differential equation. Overall, this leads to a free boundary value problem essentially hyperbolic. By using the method of characteristics, the partial differential equations constituting the invasion model are converted to Volterra integral equations. Then, the fixed point theorem is used for the uniqueness and existence result. The work is completed with numerical simulations describing the invasion of nitrite oxidizing bacteria in a biofilm initially constituted by ammonium oxidizing bacteria.     

A bidimensional phase-field model with convection for change phase of an alloy

The article analyzes a two-dimensional phase-field model for a non-stationary process of solidification of a binary alloy with thermal properties. The model allows the occurrence of fluid flow in non-solid regions, which are a priori unknown, and is thus associated to a free boundary value problem for a highly non-linear system of partial differential equations. These equations are the phase-field equation, the heat equation, the concentration equation and a modified Navier-Stokes equations obtained by the addition of a penalization term of Carman-Kozeny type which accounts for the mushy effects. A proof of existence of weak solutions for such system is given. The problem is firstly approximated and a sequence of approximate solutions is obtained by Leray-Schauder fixed point theorem. A solution is then found by using compactness argument.