Search results

We introduce a second‐order accurate time‐marching pressure‐correction algorithm to accommodate weakly‐compressible highly‐viscous liquid flows at low Mach number. As the incompressible limit is approached (Ma ≈ 0), the consistency of the compressible scheme is highlighted in recovering equivalent incompressible solutions. In the viscous‐dominated regime of low Reynolds number (zone of interest), the algorithm treats the viscous part of the equations in a semi‐implicit form. Two discrete representations are proposed to interpolate density: a piecewise‐constant form with gradient recovery and a linear interpolation form, akin to that on pressure. Numerical performance is considered on a number of classical benchmark problems for highly viscous liquid flows to highlight consistency, accuracy and stability properties. Validation bears out the high quality of performance of both compressible flow implementations, at low to vanishing Mach number. Neither linear nor constant density interpolations schemes degrade the second‐order accuracy of the original incompressible fractional‐staged pressure‐correction scheme. The piecewise‐constant interpolation scheme is advocated as a viable method of choice, with its advantages of order retention, yet efficiency in implementation.

The purpose of this paper is to consider the unsteady natural convection flow of a viscous incompressible fluid, which is induced by differential heating on the solid vertical boundary of an open‐ended rectangular cavity with the two horizontal surfaces which are permeable and maintained at the temperature of ambient fluid. Attention is focused on how the flow and heat transfer is affected by variations of the buoyancy force, as well as by the permeability of the surfaces.

An upwind finite‐difference method in conjunction with a successive over‐relaxation iteration technique is used to solve the governing boundary layer equations. To do this, the first and second derivatives were approximated by central differences and were used in the vorticity, energy and Poisson equations. To preserve the conservative property, the finite‐difference forms of the vorticity and energy equations were written in conservative form for the convective terms.

Local rate of heat transfer from the heated surface increases owing to an increase in the value of Ra. In the region near the bottom surface, the heat transfer from the left vertical surface decreases, but that increases in the region near the upper surface. Due to blowing of fluid through the permeable surfaces, the rate of heat transfer is higher than the situation where fluid is being withdrawn. This difference was found to be higher in the case of larger value of Ra.

The analysis is valid for unsteady, two‐dimensional natural convection flow of a viscous fluid filled in an open‐ended rectangular enclosure. An extension to three‐dimensional flow case is left for future work.

The method is very useful to analyze solar receiver systems, fire research, electronic cooling, brake housing of an aircraft and many environmental geothermal processes.

The results of this study may be of some interest to engineers interested in heat transfer in ventilated rooms or enclosures.

The paper's aim is to investigate the mixed convection flow of an electrically conducting and viscous incompressible fluid past an isothermal vertical surface with Joule heating in the presence of a uniform transverse magnetic field fixed relative to the surface. It was assumed that the electrical conductivity of the fluid varies linearly with the transverse velocity component.

The governing boundary layer equations were solved numerically. The boundary layer equations were first reduced to a convenient form by using two different formulations, namely, (i) the stream function formulation (SFF) and (ii) primitive variable formulation (PVF).

It was observed that both the local shear‐stress and Nusselt number increase with increasing value of local magnetic parameter, ξ.

In the present investigation, we investigated the effects of Joule heating on MHD mixed convection boundary layer flow of an electrically conducting viscous incompressible fluid past an isothermal vertical flat plate in the presence of a transverse magnetic field fixed relative to the surface of the plate. The analysis was valid for a steady, two dimensional laminar flow. An extension to three dimensional flow case is left for future work.

Here we have analyzed the problem of mixed convection flow of electrically conducting and viscous incompressible fluid past an isothermal vertical surface with viscous and Joule heating in presence of a uniform transverse magnetic field fixed relative to the surface. The work would be useful in the thermal management of heat transfer devices.

The results of this study may be of interest to engineers interested in heat exchanger design.

The purpose of this paper is to investigate the effects of thermal radiation and viscous dissipation on steady natural convection flow of a viscous incompressible fluid along a uniformly moving infinite vertical porous plate with Newtonian heating in the presence of transverse magnetic field. The governing non-linear boundary layer equations are solved by using homotopy analysis method (HAM). The effects of various system parameters on velocity and temperature fields are discussed graphically, and the numerical values for skin friction and Nusselt number are presented in tabular form.

The problem is formulated using the Boussinesq approximation under the effects of thermal radiation and transverse magnetic field. The resulting coupled system of non-linear differential equations is solved using HAM with appropriate boundary conditions for Newtonian heating of the plate. HAM is a powerful method which provides rapidly converging series solution for the velocity and temperature fields. The effects of Prandtl number, Grashof number, suction parameter, magnetic field parameter, radiation parameter and Eckert number on the fluid velocity, temperature, skin friction and Nusselt number have been investigated.

The HAM solution has been successfully applied to find the converging series solution for velocity and temperature fields in terms of pertinent system parameters. Comparison of the exact solution results agree well with the HAM solution results in the absence of Eckert number and this indicates that the HAM solutions are accurate. It is found that the velocity and temperature profiles decreases with the increase of thermal radiation and suction parameters. An increase in the magnetic field parameter leads to a rise in the fluid temperature and fall in the fluid velocity.

The present analysis is limited to steady state laminar natural convection flow only. Unsteady natural- /mixed-convection laminar flow in the presence of thermal radiation, chemical reaction and transverse magnetic field will be investigated in a future work.

The study provides very useful information for heat transfer engineers to understand the heat transfer rate when the moving vertical porous surface temperature is not known a prior. The present results have immediate relevance in the design of nuclear reactors where vertical moving porous plates are using as control rods.

The present research work is relatively original and illustrates the effects of thermal radiation, viscous dissipation and transverse magnetic field on natural convection flow past a uniformly moving infinite vertical porous plate with Newtonian heating.

The research is directed at development of an efficient and accurate technique for modelling incompressible free surface flows in which viscous effects may not be neglected. The paper describes the methodology, and gives illustrative results for simple geometries.

The pseudospectral matrix element method of discretisation is selected as the basis for the CFD technique adopted, because of its high spectral accuracy. It is implemented as a means of solving the Navier‐Stokes equations coupled with the modified compressibility method.

The viscous solver has been validated for the benchmark cases of uniform flow past a cylinder at a Reynolds number of 40, and 2D cavity flows. Results for sloshing of a viscous fluid in a tank have been successfully compared with those from a linearised analytical solution. Application of the method is illustrated by the results for the interaction of an impulsive wave with a surface piercing circular cylinder in a cylindrical tank.

The paper demonstrates the viability of the approach adopted. The limitation of small amplitude waves should be tackled in future work.

The results will have particular significance in the context of validating computations from more complex schemes applicable to arbitrary geometries.

The new methodology and results are of interest to the community of those developing numerical models of flow past marine structures.

The theoretical derivation of the start‐up laminar flow of incompressible viscous fluid in a long pipe as suggested by Szymanski, could not be verified experimentally. This leads to the checking of assumption of constant pressure gradient across the ends of the pipe, on the basis of which the theoretical development was made. Recently, the problem was again investigated for viscous fluid by Otis. In the present paper, the laminar start‐up flow of elastico‐viscous fluid in a pipe, without assuming constant pressure gradient across its ends, has been investigated. The non‐linear governing equations are solved numerically and the effects of start‐up flow parameters and elastico‐viscous parameter on the velocity distribution have been studied.

A software package is developed for the modelling and animation of viscous incompressible fluids. The full time‐dependent Navier‐Stokes equations are used to simulate 2D and 3D incompressible fluid phenomena which include shallow and deep fluid flow, transient dynamic flow, vorticity and splashing in simulated physical environments. The package also allows the inclusion of variously shaped and spaced static or moving obstacles that are fully submerged or penetrate the fluid surface. Stable numerical analysis techniques based on finite‐differences are used for the solution of the Navier‐Stokes equations. To model free‐surface fluids, a technique based on the Marker‐and‐Cell method is presented. Based on the fluid’s pressure and velocities obtained from the solution of the Navier‐Stokes equations this technique allows modelling of the fluid’s free surface either by solving a surface equation of by tracking the motion of marker particles. The latter technique is suitable for visualization of splashing and vorticity. Furthermore, an editing tool is developed for easy definition of a physical‐world which includes obstacles, boundaries and fluid properties such as viscosity, initial velocity and pressure. Using the editor, complex fluid simulations can be performed without prior knowledge of the underlying fluid dynamics equations. Finally, depending on the application fluid rendering techniques are developed using standard Silicon Graphics workstation hardware routines.

The purpose of this paper is to analyze the mangnetohydrodynamic boundary layer flow of a viscous, incompressible and electrically conducting non-Newtonian nanofluid obeying power-law model over a non-linear stretching sheet under the influence of thermal radiation with heat source/sink.

The transverse magnetic field is applied normal to the sheet. The model used for the nanofluid incorporates the effects of Brownian motion with thermophoresis in the presence of thermal radiation. On this regard, thermophoresis effect on convective heat transfer on nanofluids are investigated simultaneously. The governing partial differential equations are reduced to ordinary differential equations by suitable similarity transformations which are solved numerically by variational finite element method.

The computations carried out for some values of the power-law index, magnetic parameter, radiation parameter, Brownian motion and thermophoresis. The effect of these parameters on the velocity, temperature and nanoparticle volume fraction distribution are presented graphically. The skin friction coefficient, Nusselt number and Sherwood number for various values of the flow parameters of the problem are also presented.

To the best of the authors’ knowledge, no investigations has been reported regarding the study of non-Newtonian nanofluids which obeying power-law model over a nonlinear stretching sheet. The principal aim of this paper is to study the boundary layer MHD flow of a non-Newtonian power-law model over a non-linear stretching sheet on a quotient viscous incompressible electrically conducting with a nanofluid.

This paper aims to present a high-order algorithm implemented with the modal spectral element method and simulations of three-dimensional thermal convective flows by using the full viscous dissipation function in the energy equation. Three benchmark problems were solved to validate the algorithm with exact or theoretical solutions. The heated rotating sphere at different temperatures inside a cold planar Poiseuille flow was simulated parametrically at varied angular velocities with positive and negative rotations.

The fourth-order stiffly stable schemes were implemented and tested for time integration. To provide the hp-refinement and spatial resolution enhancement, a modal spectral element method using hierarchical basis functions was used to solve governing equations in a three-dimensional space.

It was found that the direction of rotation of the heated sphere has totally different effects on drag, lateral force and torque evaluated on surfaces of the sphere and walls. It was further concluded that the angular velocity of the heated sphere has more influence on the wall normal velocity gradient than on the wall normal temperature gradients and therefore, more influence on the viscous dissipation than on the thermal dissipation.

This paper concerns incompressible fluid flow at constant properties with up to medium temperature variations in the absence of thermal radiation and ignoring the pressure work.

This paper contributes a viable high-order algorithm in time and space for modeling convective heat transfer involving an internal heated rotating sphere with the effect of viscous heating.

Results of this paper could provide reference for related topics such as enhanced heat transfer forced convection involving rotating spheres and viscous thermal effect.

The merits include resolving viscous dissipation and thermal diffusion in stationary and rotating boundary layers with both h- and p-type refinements, visualizing the viscous heating effect with the full viscous dissipation function in the energy equation and modeling the forced advection around a rotating sphere with varied positive and negative angular velocities subject to a shear flow.

The purpose of this paper is to focus on the numerical modelling of transient natural convection flow of an incompressible viscous nanofluid past an impulsively started semi-infinite vertical plate with variable surface temperature.

The problem is governed by the coupled non-linear partial differential equations with appropriate boundary conditions. A robust, well-tested, Crank-Nicolson type of implicit finite-difference method, which is unconditionally stable and convergent, is used to solve the governing non-linear set of partial differential equations.

The local and average values of the skin-friction coefficient (viscous drag) and the average Nusselt number (the rate of heat transfer) decreased, while the local Nusselt number increased for all nanofluids, namely, aluminium oxide-water, copper-water, titanium oxide-water and silver-water with an increase in the temperature exponent m. Selecting aluminium oxide as the dispersing nanoparticles leads to the maximum average Nusselt number (the rate of heat transfer), while choosing silver as the dispersing nanoparticles leads to the minimum local Nusselt number compared to the other nanofluids for all values of the temperature exponent m. Also, choosing silver as the dispersing nanoparticles leads to the minimum skin-friction coefficient (viscous drag), while selecting aluminium oxide as the dispersing nanoparticles leads to the maximum skin-friction coefficient (viscous drag) for all values of the temperature exponent m.

The Brinkman model for dynamic viscosity and Maxwell-Garnett model for thermal conductivity are employed. The governing boundary layer equations are written according to The Tiwari-Das nanofluid model. A range of nanofluids containing nanoparticles of aluminium oxide, copper, titanium oxide and silver with nanoparticle volume fraction range less than or equal to 0.04 are considered.

The present simulations are relevant to nanomaterials thermal flow processing in the chemical engineering and metallurgy industries. This study also provides an important benchmark for further simulations of nanofluid dynamic transport phenomena of relevance to materials processing, with alternative computational algorithms (e.g. finite element methods).

This paper is relatively original and illustrates the influence of variable surface temperature on transient natural convection flow of a viscous incompressible nanofluid and heat transfer from an impulsively started semi-infinite vertical plate.