Computing Viscous Flow in an Elastic Tube

We have developed a numerical method for simulating viscous flow through a compliant closed tube, driven by a pair of fluid source and sink. As is natural for tubular flow simulations, the problem is formulated in axisymmetric cylindrical coordinates, with fluid flow described by the Navier-Stokes equations. Because the tubular walls are assumed to be elastic, when stretched or compressed they exert forces on the fluid. Since these forces are singularly supported along the boundaries, the fluid velocity and pressure fields become unsmooth. To accurately compute the solution, we use the velocity decomposition approach, according to which pressure and velocity are decomposed into a singular part and a remainder part. The singular part satisfies the Stokes equations with singular boundary forces. Because the Stokes solution is unsmooth, it is computed to second-order accuracy using the immersed interface method, which incorporates known jump discontinuities in the solution and derivatives into the finite difference stencils. The remainder part, which satisfies the Navier-Stokes equations with a continuous body force, is regular. The equations describing the remainder part are discretized in time using the semi-Lagrangian approach, and then solved using a pressure-free projection method. Numerical results indicate that the computed overall solution is second-order accurate in space, and the velocity is second-order accurate in time.     

The study of Heat and Mass Transfer in a Visco elastic fluid due to a continuous stretching surface using Homotopy Analysis Method

In this paper, an approximate analytical solution is derived for the flow velocity and temperature due to the laminar, two-dimensional flow of non-Newtonian incompressible visco elastic fluid due to a continuous stretching surface. The surface is stretched with a velocity proportional to the distance $x$ along the surface. The surface is assumed to have either power-law heat flux or power-law temperature distribution. The presence of source/sink and the effect of uniform suction and injection on the flow are considered for analysis. An approximate analytical solution has been obtained using Homotopy Analysis Method(HAM) for various values of visco elastic parameter, suction and injection rates. Optimal values of the convergence control parameters are computed for the flow variables. It was found that the computational time required for averaged residual error calculation is very very small compared to the computation time of exact squared residual errors. The effect of mass transfer parameter, visco elastic parameter, source/sink parameter and the power law index on flow variables such as velocity, temperature profiles, shear stress, heat and mass transfer rates are discussed.     

Convective heat and mass transfer in three-dimensional mixed convection flow of viscoelastic fluid in presence of chemical reaction and heat source/sink

Heat and mass transfer effects in the three-dimensional mixed convection flow of a viscoelastic fluid with internal heat source/sink and chemical reaction have been investigated in the present work. The flow generation is because of an exponentially stretching surface. Magnetic field normal to the direction of flow is considered. Convective conditions at the surface are also encountered. Appropriate similarity transformations are utilized to reduce the boundary layer partial differential equations into the ordinary differential equations. The homotopy analysis method is used to develop the solution expressions. Impacts of different controlling parameters such as ratio parameter, Hartman number, internal heat source/sink, chemical reaction, mixed convection, concentration buoyancy parameter and Biot numbers on the velocity, temperature and concentration profiles are analyzed. The local Nusselt and Sherwood numbers are sketched and examined.     

Combined effects of non-uniform heat source/sink and thermal radiation on heat transfer over an unsteady stretching permeable surface

The present paper is concerned with the study of flow and heat transfer characteristics in the unsteady laminar boundary layer flow of an incompressible viscous fluid over continuously stretching permeable surface in the presence of a non-uniform heat source/sink and thermal radiation. The unsteadiness in the flow and temperature fields is because of the time-dependent stretching velocity and surface temperature. Similarity transformations are used to convert the governing time-dependent nonlinear boundary layer equations for momentum and thermal energy are reduced to a system of nonlinear ordinary differential equations containing Prandtl number, non-uniform heat source/sink parameter, thermal radiation and unsteadiness parameter with appropriate boundary conditions. These equations are solved numerically by applying shooting method using Runge–Kutta–Fehlberg method. Comparison of numerical results is made with the earlier published results under limiting cases. The effects of the unsteadiness parameter, thermal radiation, suction/injection parameter, non-uniform heat source/sink parameter on flow and heat transfer characteristics as well as on the local Nusselt number are shown graphically.     

Analytical solutions for a time-fractional axisymmetric diffusion–wave equation with a source term

In this paper, a time-fractional axisymmetric diffusion–wave equation with a source term is considered in cylindrical coordinates. The analytical solution is obtained with the help of an integral transform method and some properties of special functions. In addition, we discuss two kinds of different boundary conditions and different forms of the source term. Finally, we analyze the effects of the fractional derivative on the solutions by using numerical results and find that sub-diffusion phenomena and oscillations exist.     

Analytical vortex and source surface distributions for flow about elliptical cylinders

Summary A cylindrical body is introduced into a potential flow and the resulting disturbance is treated using only an analytical singularity distribution on the surface of the body. Conformai transformation methods, using previously published results for singularity distributions on a circular cylinder, enable the derivation of the required form of vorticity distribution on the transformed body. This approach, however, does not facilitate the derivation of the analogous source/sink distribution. For an elliptic cylinder, examples of vorticity distributions are presented and a new solution is obtained giving the source/sink distribution on the surface of the cylinder in a uniform stream at incidence.     

Propagation of harmonic waves in prestressed fiber viscoelastic thick tubes filled with a viscous dusty fluid

In this work, propagation of harmonic waves in initially stressed cylindrical viscoelastic thick tubes filled with a Newtonian fluid is studied. The tube, subjected to a static inner pressure P_i and a positive axial stretch λ, will be considered as an incompressible viscoelastic and fibrous material. The fluid is assumed as an incompressible, viscous and dusty fluid. The field equations for the fluid are obtained in the cylindrical coordinates. The governing differential equations of the tube’s viscoelastic material are obtained also in the cylindrical coordinates utilizing the theory of small deformations superimposed on large initial static deformations. For the axially symmetric motion the field equations are solved by assuming harmonic wave solutions. A closed form solution can be obtained for equations governing the fluid body, but due to the variability of the coefficients of resulting differential equations of the solid body, such a closed form solution is not possible to obtain. For that reason, equations for the solid body and the boundary conditions are treated numerically by the finite-difference method to obtain the effects of the thickness of the tube on the wave characteristics. Dispersion relation is obtained using the long wave approximation and, the wave velocities and the transmission coefficients are computed.     

Geometry of singularities for the steady Boussinesq equations

Analysis and computations are presented for singularities in the solution of the steady Boussinesq equations for two-dimensional, stratified flow. The results show that for codimension 1 singularities, there are two generic singularity types for general solutions, and only one generic singularity type if there is a certain symmetry present. The analysis depends on a special choice of coordinates, which greatly simplifies the equations, showing that the type is exactly that of one dimensional Legendrian singularities, generalized so that the velocity can be infinite at the singularity. The solution is viewed as a surface in an appropriate compactified jet space. Smoothness of the solution surface is proved using the Cauchy-Kowalewski Theorem, which also shows that these singularity types are realizable. Numerical results from a special, highly accurate numerical method demonstrate the validity of this geometric analysis. A new analysis of general Legendrian singularities with blowup, i.e., at which the derivative may be infinite, is also presented, using projective coordinates.Research supported in part by the ARPA under URI grant number #N00014092-J-1890.Research supported in part by the NSF under grant number #DMS93-02013.Research supported in part by the NSF under grant #DMS-9306488.     

A numerical study of flows driven by a rotating magnetic field in a square container

The influence of the geometry on the magnetically driven flow is studied by means of numerical simulations. Low–frequency, low–induction and low–interaction conditions are assumed. The rotating magnetic field (RMF) gives rise to a time–independent magnetic body force, computed via the electrical potential equation and Ohm's law and a time–dependent part that is neglected due to the low interaction parameter. Flow results of the cylindrical and square container are compared with respect to the magnetic body force, time–averaged velocity fields, first flow instabilities and Reynolds stress tensors. The dependency of the maximal velocity magnitude and the intensity of the magnetic induction is identical in axisymmetric and non–axisymmetric containers and in good agreement with Davidson's theory. However, significant differences are recognized, for instance, in the distribution of the Reynolds stress tensors. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Lie group analysis for the effect of temperature-dependent fluid viscosity with thermophoresis and chemical reaction on MHD free convective heat and mass transfer over a porous stretching surface in the presence of heat source/sink

This paper concerns with a steady two-dimensional flow of an electrically conducting incompressible fluid over a vertical stretching sheet. The flow is permeated by a uniform transverse magnetic field. The fluid viscosity is assumed to vary as a linear function of temperature. A scaling group of transformations is applied to the governing equations. The system remains invariant due to some relations among the parameters of the transformations. After finding three absolute invariants a third-order ordinary differential equation corresponding to the momentum equation and two second-order ordinary differential equation corresponding to energy and diffusion equations are derived. The equations along with the boundary conditions are solved numerically. It is found that the decrease in the temperature-dependent fluid viscosity makes the velocity to decrease with the increasing distance of the stretching sheet. At a particular point of the sheet the fluid velocity decreases with the decreasing viscosity but the temperature increases in this case. It is found that with the increase of magnetic field intensity the fluid velocity decreases but the temperature increases at a particular point of the heated stretching surface. Impact of thermophoresis particle deposition with chemical reaction in the presence of heat source/sink plays an important role on the concentration boundary layer. The results thus obtained are presented graphically and discussed.     

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.     

Selfsimilar solutions of the problem of heat and mass transfer in a saturated porous medium with a volume heat source

Selfsimilar solutions of a system of stationary equations of heat condunction and filtration of molten material in the presence of a volume heat source generated by absorption of the energy of electromagnetic radiation, are considered. The possibility of the existence of a self-similar solution in the case of various (plane, cylindrical and spherical) spatial symmetries is studied. The existence of a selfsimilar solution is shown for the axisymmetric case when the radiation obeys a prescribed law. The influence of the surface volume heating and convective heat transfer due to filtration is studied. A solution for the case when the filtration of the molten phase is quasistationary is also investigated.     

A general solution of equations of equilibrium in linear elasticity

A general solution of equations of equilibrium in linear elasticity is presented in cylindrical coordinates in terms of three harmonic functions describing an arbitrary displacement field. The structure of this solution is similar to the general solution given by Love (Kelvin’s solution) in spherical coordinates. Galerkin vector representation of our solution leads to an integral connecting the harmonic functions. The connections to Papkovich–Neuber and Muki’s general representations are also provided. Suitable choices of the harmonic functions in our new representation yield general solutions for axisymmetric deformations due to Love, Boussinesq and Michell. Some unbounded deformations induced by singular forces are tabulated in terms of the scalar harmonic functions to justify the simple nature of our representation. Exact solution of the half-space boundary value problem is also provided to demonstrate the power of our approach. The stress components computed via our solution are also listed (see the Appendix).     

光滑区域上二维无黏性无热传导Boussinesq方程组与三维轴对称不可压Euler方程组的指数增长全局光滑解

研究二维无黏性无热传导Boussinesq方程组和三维轴对称不可压Euler方程组光滑解的增长情况,找各种区域使其上的方程组有快增长的解。对Boussinesq方程组,通过选取初始温度和速度的一个分量,可以把方程去耦为两部分。从关于涡量的部分求出涡量、速度场和使结论成立的区域,从关于温度的部分,可见温度的高阶导的增长仅依赖于速度场的一个分量。通过适当选取该分量,得到温度高阶导有指数增长的全局光滑解。对轴对称Euler方程组做类似的处理,适当选取速度场的径向分量,可把方程组去耦,最终得到一类光滑区域,在其上方程组有指数增长全局光滑解。该研究把Chae、Constantin、Wu对一个二维锥形区域上无黏性无热传导Boussinesq方程的结果,推广到一类光滑区域上, 并把他们的方法应用到三维轴对称不可压Euler方程组, 得到了类似的结果。     

The effect of variable viscosity on MHD viscoelastic fluid flow and heat transfer over a stretching sheet

An analysis has been carried out to study the momentum and heat transfer characteristics in an incompressible electrically conducting non-Newtonian boundary layer flow of a viscoelastic fluid over a stretching sheet. The partial differential equations governing the flow and heat transfer characteristics are converted into highly non-linear coupled ordinary differential equations by similarity transformations. The effect of variable fluid viscosity, Magnetic parameter, Prandtl number, variable thermal conductivity, heat source/sink parameter and thermal radiation parameter are analyzed for velocity, temperature fields, and wall temperature gradient. The resultant coupled highly non-linear ordinary differential equations are solved numerically by employing a shooting technique with fourth order Runge–Kutta integration scheme. The fluid viscosity and thermal conductivity, respectively, assumed to vary as an inverse and linear function of temperature. The analysis reveals that the wall temperature profile decreases significantly due to increase in magnetic field parameter. Further, it is noticed that the skin friction of the sheet decreases due to increase in the Magnetic parameter of the flow characteristics.     

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.     

Natural convection on a vertical stretching cylinder

The flow and natural (or mixed) convection due to a vertical stretching cylinder is studied. Using similarity transforms, the Navier-Stokes and energy equations reduce to a set of nonlinear ordinary differential equations. Asymptotic analysis for large Reynolds numbers shows the relation between axisymmetric flow and two-dimensional flow. Due to the algebraic decay of the similarity functions, numerical integration is performed using a compressed coordinate. The axial velocity is composed of forced convection due to stretching and natural convection from the heated cylinder. The heat transfer increases with both the Reynolds number and the Prandtl number. The result is also a rare similarity solution of the free convection and Navier-Stokes equations.     

Effect of variable viscosity on MHD non-Darcy mixed convective heat transfer over a stretching sheet embedded in a porous medium with non-uniform heat source/sink

An analysis has been presented to investigate the effect of temperature-dependent viscosity on non-Darcy MHD mixed convective heat transfer past a porous medium by taking into account of Ohmic dissipation and non-uniform heat source/sink. Thermal boundary layer equation takes into account of viscous dissipation and Ohmic dissipation due to transverse magnetic field and electric field. The governing fundamental equations are first transformed into system of ordinary differential equations using self-similarity transformation and are solved numerically by using the fifth-order Runge–Kutta–Fehlberg method with shooting technique for various values of the physical parameters. The effects of variable viscosity, porosity, Eckert number, Prandtl number, magnetic field, electric field and non-uniform heat source/sink parameters on velocity and temperature profiles are analyzed and discussed. Favorable comparisons with previously published work on various special cases of the problem are obtained. Numerical results on the development of the local skin-friction co-efficient and local Nusselt number with non-uniform heat source/sink are tabulated for various physical parameters to show the interesting aspects of the solution.     

Genesis of a swirling source/sink line into an updraft

Several mathematical models applied to tornadoes consist of exact and axisymmetric solutions of the steady and incompressible Navier—Stokes equations. These models studied by J. Serrin [5] and M.A. Gol'dshtik and V.N. Sthern (1995) describe families of fluid motions vanishing at the ground and are restricted for not developing neither a source nor a sink near the vortex line. Therefore, J. Serrin showed that the flow patterns of the resulting velocity field may have some realistic characteristics to model the mature phase of the lifetime of a tornado, in comparison with atmospheric observations. On another hand, no reason has been explicated to motivate the absence of a source/sink vortex line. We present here a result of local existence and uniqueness of a family of fluid motions, leading to the genesis of such lines inside a non-rotating updraft which do not develops neither a source nor a sink.     

Regularity of the 3D Navier–Stokes equations with viewpoint of 2D flow

The regularity of 2D Navier–Stokes flow is well known. In this article we study the relationship of 3D and 2D flow, and the regularity of the 3D Naiver–Stokes equations with viewpoint of 2D equations. We consider the problem in the Cartesian and in the cylindrical coordinates.