Numerical investigation of transient heat transfer to hydromagnetic channel flow with radiative heat and convective cooling

This present study consists of a numerical investigation of transient heat transfer in channel flow of an electrically conducting variable viscosity Boussinesq fluid in the presence of a magnetic field and thermal radiation. The temperature dependent nature of viscosity is assumed to follow an exponentially model and the system exchanges heat with the ambient following Newton’s law of cooling. The governing nonlinear equations of momentum and energy transport are solved numerically using a semi-implicit finite difference method. Solutions are presented in graphical form and given in terms of fluid velocity, fluid temperature, skin friction and heat transfer rate for various parametric values. Our results reveal that combined effect of thermal radiation, magnetic field, viscosity variation and convective cooling have significant impact in controlling the rate of heat transfer in the boundary layer region.     

Ⅰ型拉伸试样金属注射成形多相流动的三维数值模拟

将粉末注射成形充模过程视为粘结剂、粉末和空气的三相流动过程,基于多相流理论给出金属注射成形的多相流动控制方程。根据金属注射成形的工艺特点确定多相流动控制方程的初边界条件,用多相流数值分析软件CFX对Ⅰ型拉伸试样的铁粉注射成形充模流动过程实现了三维数值模拟,分析了模腔中不同位置处粘结剂一粉末流动的速度曲线。数值模拟结果表明在Ⅰ型拉伸试样金属注射成形充模流动过程中发生过两相分离的现象,数值模拟的瞬态信息可用于金属注射成形产品缺陷的分析与控制。     

Finite element modelling of surface waves in a channel

A variational formulation of the vertically-integrated differential equations for free surface wave motion is presented. A finite element model is derived for solving this nonlinear system of hydrodynamic equations. The time integration scheme employed is discussed and the results obtained demonstrate its good stability and accuracy.Several applications of the model are considered: the first problem is an open channel of uniform depth and the second an open channel of linearly varying depth. The ‘inflow’ boundary condition is prescribed in terms of the velocity which represents a wavemaker and/or a flow source, while the ‘outflow’ boundary condition is specified in terms of the water elevation. The outflow condition is adjusted for two cases, a reflecting boundary (finite channel) and a non-reflecting boundary (open-ended channel). The latter boundary condition is examined in some detail and the results obtained show that the numerical model can produce the non-reflecting boundary that is similar to the analytical radiation condition for waves. Computational results for a third problem, involving wave reflection from a submerged cylinder, are also presented and compared with both experimental data and analytical predictions.The simplicity and the performance of the computational model suggest that free surface waves can be simulated without excessively complicated numerical schemes. The ability of the model to simulate outflow boundary conditions properly is of special importance since these conditions present serious problems for many numerical algorithms.     

Creeping flow analysis of slightly non-Newtonian fluid in a uniformly porous slit

This paper provides the analysis of the steady, creeping flow of a special class of slightly viscoelastic, incompressible fluid through a slit having porous walls with uniform porosity. The governing two dimensional flow equations along with non-homogeneous boundary conditions are non-dimensionalized. Recursive approach is used to solve the resulting equations. Expressions for stream function, velocity components, volumetric flow rate, pressure distribution, shear and normal stresses in general and on the walls of the slit, fractional absorption and leakage flux are derived. Points of maximum velocity components are also identified. A graphical study is carried out to show the effect of porosity and non-Newtonian parameter on above mentioned resulting expressions. It is observed that axial velocity of the fluid decreases with the increase in porosity and non-Newtonian parameter. The outcome of this theoretical study has significant importance both in industry and biosciences.     

Lattice Boltzmann equation model for multi-component multi-phase flow with high density ratios

Multi-component multi-phase (MCMP) flows are very common in engineering or industrial problems, as well as in nature. Because the lattice Boltzmann equation (LBE) model is based on microscopic models and mesoscopic kinetic equations, it offers many advantages for the study of multi-component or multi-phase flow problems. While the original formulation of Shan and Chen’s (SC) model can incorporate some MCMP flow scenarios, the density ratio of the different components is greatly restricted to less than approximately 2.0. This obviously limits the applications of this MCMP LBE model. Hence, based on the original SC MCMP model and the improvements in the single-component multi-phase (SCMP) flow model reported by Yuan and Schaefer, we have developed a new model that can simulate a MCMP system with a high density ratio.     

Group theoretic method for unsteady free convection flow of a micropolar fluid along a vertical plate in a thermally stratified medium

The group theoretic approach is applied for solving the problem of unsteady natural convection flow of micropolar fluid along a vertical flat plate in a thermally stratified medium. The application of two-parameter transformation group reduces the number of independent variables in the governing system consisting of partial differential equations and a set of auxiliary conditions from three to only one independent variable, and consequently the system of governing partial differential equations with boundary conditions reduces to a system of ordinary differential equations with appropriate boundary conditions. Numerical solution of the velocity, microrotation and heat transfer have been obtained. The possible forms of the ambient temperature variation with position and time are derived.     

Group solution for unsteady free-convection flow from a vertical moving plate subjected to constant heat flux

The problem of heat and mass transfer in an unsteady free-convection flow over a continuous moving vertical sheet in an ambient fluid is investigated for constant heat flux using the group theoretical method. The nonlinear coupled partial differential equation governing the flow and the boundary conditions are transformed to a system of ordinary differential equations with appropriate boundary conditions. The obtained ordinary differential equations are solved numerically using the shooting method. The effect of Prandlt number on the velocity and temperature of the boundary-layer is plotted in curves. A comparison with previous work is presented.     

Magnetohydrodynamic non-Darcy mixed convection heat transfer from a vertical heated plate embedded in a porous medium with variable porosity

A numerical model is developed to study magnetohydrodynamics (MHD) mixed convection from a heated vertical plate embedded in a Newtonian fluid saturated sparsely packed porous medium by considering the variation of permeability, porosity and thermal conductivity. The boundary layer flow in the porous medium is governed by Forchheimer–Brinkman extended Darcy model. The conservation equations that govern the problem are reduced to a system of non-linear ordinary differential equations by using similarity transformations. Because of non-linearity, the governing equations are solved numerically. The effects of magnetic field on velocity and temperature distributions are studied in detail by considering uniform permeability (UP) and variable permeability (VP) of the porous medium and the results are discussed graphically. Besides, skin friction and Nusselt number are also computed for various physical parameters governing the problem under consideration. It is found that the inertial parameter has a significant influence in increasing the flow field and the rate of heat transfer for variable permeability case. The important finding of the present work is that the magnetic field has considerable effects on the boundary layer velocity and on the rate of heat transfer for variable permeability of the porous medium. Further, the results obtained under the limiting conditions were found to be in good agreement with the existing ones.     

Thermal radiation effects on magnetohydrodynamic free convection heat and mass transfer from a sphere in a variable porosity regime

A mathematical model is presented for multiphysical transport of an optically-dense, electrically-conducting fluid along a permeable isothermal sphere embedded in a variable-porosity medium. A constant, static, magnetic field is applied transverse to the cylinder surface. The non-Darcy effects are simulated via second order Forchheimer drag force term in the momentum boundary layer equation. The surface of the sphere is maintained at a constant temperature and concentration and is permeable, i.e. transpiration into and from the boundary layer regime is possible. The boundary layer conservation equations, which are parabolic in nature, are normalized into non-similar form and then solved numerically with the well-tested, efficient, implicit, stable Keller-box finite difference scheme. Increasing porosity (ε) is found to elevate velocities, i.e. accelerate the flow but decrease temperatures, i.e. cool the boundary layer regime. Increasing Forchheimer inertial drag parameter (Λ) retards the flow considerably but enhances temperatures. Increasing Darcy number accelerates the flow due to a corresponding rise in permeability of the regime and concomitant decrease in Darcian impedance. Thermal radiation is seen to reduce both velocity and temperature in the boundary layer. Local Nusselt number is also found to be enhanced with increasing both porosity and radiation parameters.     

Analytical approximate solution of the cooling problem by Adomian decomposition method

The Adomian decomposition method (ADM) can provide analytical approximation or approximated solution to a rather wide class of nonlinear (and stochastic) equations without linearization, perturbation, closure approximation, or discretization methods. In the present work, ADM is employed to solve the momentum and energy equations for laminar boundary layer flow over flat plate at zero incidences with neglecting the frictional heating. A trial and error strategy has been used to obtain the constant coefficient in the approximated solution. ADM provides an analytical solution in the form of an infinite power series. The effect of Adomian polynomial terms is considered and shows that the accuracy of results is increased with the increasing of Adomian polynomial terms. The velocity and thermal profiles on the boundary layer are calculated. Also the effect of the Prandtl number on the thermal boundary layer is obtained. Results show ADM can solve the nonlinear differential equations with negligible error compared to the exact solution.     

Stochastic immersed boundary method incorporating thermal fluctuations

We give a brief introduction to the stochastic immersed boundary method which allows for simulation of small length-scale physical systems in which elastic structures interact with a fluid flow in the presence of thermal fluctuations. The conventional immersed boundary method is extended to account for thermal fluctuations by introducing stochastic forcing terms in the fluid equations. This gives a system of stiff SPDE's for which standard numerical approaches perform poorly. We discuss a numerical method derived using stochastic calculus to overcome the stiff features of the equations. We then discuss results which indicate that the method captures physical features predicted by statistical mechanics for small length-scale systems. The stochastic immersed boundary method holds promise as a numerical approach in simulating microscopic mechanical systems in which thermal fluctuations play a fundamental role. A more detailed discussion of this work is given in [1, 2, 3]. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Heat transfer in granular materials: effects of nonlinear heat conduction and viscous dissipation

In this paper, we study the heat transfer in a one‐dimensional fully developed flow of granular materials down a heated inclined plane. For the heat flux vector, we use a recently derived constitutive equation that reflects the dependence of the heat flux vector on the temperature gradient, the density gradient, and the velocity gradient in an appropriate frame invariant formulation. We use two different boundary conditions at the inclined surface: a constant temperature boundary condition and an adiabatic condition. A parametric study is performed to examine the effects of the material dimensionless parameters. The derived governing equations are coupled nonlinear second‐order ordinary differential equations, which are solved numerically, and the results are shown for the temperature, volume fraction, and velocity profiles. Copyright © 2013 John Wiley & Sons, Ltd.     

An application of a boundary element method to natural convection

The direct boundary element method is applied to the numerical modelling of thermal fluid flow in a transient state. The Navier-Stokes equations are considered under the Boussinesq approximation and the viscous thermal flow equations are expressed in terms of stream function, vorticity, and temperature in two dimensions. Boundary integral equations are derived using logarithmic potential and time-dependent heat potential as fundamental solutions. Boundary unknowns are discretized by linear boundary elements and flow domains are divided into a series of triangular cells. Charged points are translated upstream in the numerical evaluation of convective terms. Unknown stream function, vorticity, and temperature are staggered in the computational scheme.

Simple iteration is found to converge to the quasi steady-state flow. Boundary solutions for two-dimensional examples at a Reynolds number 100 and Grashoff number 10⁷ are obtained.     

Two-parameter extremum problems of boundary control for stationary thermal convection equations

Two-parameter extremum problems of boundary control are formulated for the stationary thermal convection equations with Dirichlet boundary conditions for velocity and with mixed boundary conditions for temperature. The cost functional is defined as the root mean square integral deviation of the desired velocity (vorticity, or pressure) field from one given in some part of the flow region. Controls are the boundary functions involved in the Dirichlet condition for velocity on the boundary of the flow region and in the Neumann condition for temperature on part of the boundary. The uniqueness of the extremum problems is analyzed, and the stability of solutions with respect to certain perturbations in the cost functional and one of the functional parameters of the original model is estimated. Numerical results for a control problem associated with the minimization of the vorticity norm aimed at drag reduction are discussed.     

Systematic derivation of an asymptotic model for the dynamics of curved viscous fibers

This paper presents a slender body theory for the dynamics of a curved inertial viscous Newtonian fiber. Neglecting surface tension and temperature dependence, the fiber flow is modeled as a three‐dimensional free boundary value problem in terms of instationary incompressible Navier–Stokes equations. From regular asymptotic expansions in powers of the slenderness parameter, leading‐order balance laws for mass (cross‐section) and momentum are derived that combine the unrestricted motion of the fiber centerline with the inner viscous transport. The physically reasonable form of the one‐dimensional fiber model results thereby from the introduction of the intrinsic velocity that characterizes the convective terms. For the numerical investigation of the viscous, gravitational and rotational effects on the fiber dynamics, a finite volume approach on a staggered grid with implicit upwind flux discretization is applied. Copyright © 2007 John Wiley & Sons, Ltd.     

Natural convection flow from a horizontal circular cylinder with uniform heat flux in presence of heat generation

Natural convection boundary layer laminar flow from a horizontal circular cylinder with uniform heat flux in presence of heat generation has been investigated. The governing boundary layer equations are transformed into a non-dimensional form and the resulting non-linear systems of partial differential equations, which are solved numerically by two distinct methods namely: (i) implicit finite difference method together with the Keller-box scheme and (ii) perturbation solution technique. The results of the surface shear stress in terms of local skin-friction and the rate of heat transfer in terms of local Nusselt number, velocity distribution, velocity vectors, temperature distribution as well as streamlines, isotherms and isolines of pressure are shown by graphically for a selection of parameter set consisting of heat generation parameter.     

A Mixed Analytical/Numerical Method for Velocity and Heat Transfer of Laminar Power-Law Fluids

This paper presents a relatively simple numerical method to investigate the flow and heat transfer of laminar power-law fluids over a semi-infinite plate in the presence of viscous dissipation and anisotropy radiation. On one hand, unlike most classical works, the effects of power-law viscosity on velocity and temperature fields are taken into account when both the dynamic viscosity and the thermal diffusivity vary as a power-law function. On the other hand, boundary layer equations are derived by Taylor expansion, and a mixed analytical/numerical method (a pseudo-similarity method) is proposed to effectively solve the boundary layer equations. This method has been justified by comparing its results with those of the original governing equations obtained by a finite element method. These results agree very well especially when the Reynolds number is large. We also observe that the robustness and accuracy of the algorithm are better when thermal boundary layer is thinner than velocity boundary layer.     

The effects of MHD and temperature dependent viscosity on the flow of non-Newtonian nanofluid in a pipe: Analytical solutions

This article examines the magnetohydrodynamic (MHD) flow of non-Newtonian nanofluid in a pipe. The temperature of the pipe is assumed to be higher than the temperature of the fluid. In particular two temperature dependent viscosity models, have been considered. The nonlinear partial differential equations along with the boundary conditions are first cast into a dimensionless form and then the equations are solved by homotopy analysis method (HAM). Explicit analytical expressions for the velocity field, the temperature distribution and nano concentration have been derived analytically. The effects of various physical parameters on velocity, temperature and nano concentration are discussed by using graphical approach.     

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.     

Stability of solutions to extremal problems of boundary control for stationary heat convection equations

We study extremal problems of boundary control for stationary heat convection equations with Dirichlet boundary conditions on velocity and temperature. As the cost functional we choose the mean square integral deviation of the required temperature field from a given temperature field measured in some part of the flow region. The controls are functions appearing in the Dirichlet conditions on velocity and temperature. We prove the stability of solutions to these problems with respect to certain perturbations of both the quality functional and one of the known functions appearing in the original equations of the model.