Computation of heat transfer enhancement in a plate‐fin heat exchanger with triangular inserts and delta wing vortex generator

Longitudinal vortices disrupt the growth of the thermal boundary layer, thereby the vortex generators producing the longitudinal vortices are well known for the enhancement of heat transfer in compact heat exchangers. The present investigation determines the heat transfer characteristics with secondary flow analysis in plate fin triangular ducts with delta wing vortex generators. This geometrical configuration is investigated for various angles of attack of the wing i.e. 15°, 20°, 26° and 37° and Reynolds numbers 100 and 200. The constant wall temperature boundary condition is used. The solution of the complete Navier Stokes equation and the energy equation is carried out using the staggered grid arrangement. The performance of the combination of triangular secondary fins and delta wing with stamping on slant surfaces has also been studied. Copyright © 2009 John Wiley & Sons, Ltd.     

The laminar free‐convection boundary‐layer flow about a heated and rotating down‐pointing vertical cone in the presence of a transverse magnetic field

This paper deals with the study of the laminar free‐convection boundary‐layer flow about a heated and rotating down‐pointing vertical cone in the presence of a transverse magnetic field. Two cases of heat transfer analysis are discussed. These are: (i) the rotating cone with prescribed surface temperature and (ii) the rotating cone with prescribed surface heat flux. By means of similarity transformation, the governing partial differential equations are reduced into highly non‐linear ordinary differential equations. The resulting non‐linear system has been solved analytically using a very efficient technique, namely homotopy analysis method. Expressions for velocity and temperature fields are developed in a series form. The influence of various pertinent parameters is also seen on the velocity and temperature fields. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical analysis on the propulsive performance and vortex shedding of fish‐like travelling wavy plate

Numerical analysis is carried out to investigate viscous flow over a travelling wavy plate undergoing lateral motion in the form of a streamwise travelling wave, which is similar to the backbone undulation of swimming fish. The two‐dimensional incompressible Navier–Stokes equations are solved using the finite element technique with the deforming‐spatial‐domain/stabilized space–time formulation. The objective of this study is to elucidate hydrodynamic features of flow structure and vortex shedding near the travelling wavy plate and to get into physical insights to the understanding of fish‐like swimming mechanisms in terms of drag reduction and optimal propulsive performance. The effects of some typical parameters, including the phase speed, amplitude, and relative wavelength of travelling wavy plate, on the flow structures, the forces, and the power consumption required for the propulsive motion of the plate are analysed. These results predicted by the present numerical analysis are well consistent with the available data obtained for the wave‐like swimming motion of live fish in nature. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical analysis of mixed convection in three‐dimensional rectangular channel with flush‐mounted heat sources based on field synergy principle

In this paper the fluid flow and heat transfer characteristics of mixed convection in three‐dimensional rectangular channel with four heat sources are investigated numerically. The SIMPLEC algorithm is applied to deal with the coupling between pressure and velocity, and a new high‐order stability‐guaranteed second‐order difference (SGSD) scheme is adopted to discretize the convection term. The influence of four parameters is studied: Richardson number, heat source distribution, channel height and inclination angle. The numerical results are analysed from the viewpoint of the field synergy principle, which says that the enhanced convective heat transfer is related not only to the velocity field and temperature field, but also to the synergy between them. It is found that the effects of the four parameters on the thermal performance can all be explained with the field synergy principle. To obtain better electronic cooling, the synergy between the velocity and temperature gradient should be increased when other conditions are unchanged. Copyright © 2006 John Wiley & Sons, Ltd.     

Effects of inclination angle on natural convection in enclosures filled with Cu–water nanofluid

Effects of inclination angle on natural convection heat transfer and fluid flow in a two-dimensional enclosure filled with Cu-nanofluid has been analyzed numerically. The performance of nanofluids is tested inside an enclosure by taking into account the solid particle dispersion. The angle of inclination is used as a control parameter for flow and heat transfer. It was varied from  = 0° to  = 120°. The governing equations are solved with finite-volume technique for the range of Rayleigh numbers as 10³  Ra  10⁵. It is found that the effect of nanoparticles concentration on Nusselt number is more pronounced at low volume fraction than at high volume fraction. Inclination angle can be a control parameter for nanofluid filled enclosure. Percentage of heat transfer enhancement using nanoparticles decreases for higher Rayleigh numbers.     

Internal flow numerical simulation of a cross‐flow fan in refrigerant operating condition using user‐defined functions

In the paper, a cross‐flow fan in refrigerant operating condition is systematically simulated using user‐defined functions. Three‐dimensional simulations are acquired with Navier–Stokes equations coupled with k–ε turbulence model, and internal flow characteristics of an indoor split‐type air conditioner are obtained, which is mainly composed of cross‐flow fan and heat exchanger. It has systematically been simulated in the isothermal flow condition that the performance of cross‐flow fan may be reduced easily with dry or humid air, and in the refrigerant operating condition in which user‐defined functions are applied to the humid air, considered as a mixture of dry air and vapor. A density‐modulated function is adopted to deal with the condensation of the vapor at the heat‐transfer region approximately. The results show flow mechanism of the two gas‐phase flow, including phase‐vary process. The distribution of the parameters is not uniform at the inlet of the machine, the intensity and position of pressure and velocity vary along the axial direction of the fan, the distribution of vapor volume fraction and turbulent intensity in heat‐transfer region is obtained, and the external characteristic data of the indoor machine are obtained and analyzed. Compared with the experimental data, the calculated characteristic curves and designed parameters are on target. © British Crown Copyright 2010/MOD. Reproduced with permission. Published by John Wiley & Sons, Ltd.     

On the performance of sudden‐expansion pipe without/with cross‐flow injection: experimental and numerical studies

The paper explores the possibilities that different turbulence closures offer, for in‐depth analysis of a complex flow. The case under investigation is steady, turbulent flow in a pipe with sudden expansion without/with normal‐to‐wall injection through jets. This is a typical geometry where generation of turbulence energy takes place, due to sudden change in boundary conditions. This study is aimed at investigating the capability of a developed computational program by the present authors with three different turbulence models to calculate the mean flow variables. Three two‐equation models are implemented, namely the standard linear k ? ε model, the low Reynolds number k ? ε model and the cubic nonlinear eddy viscosity (NLEV) k ? ε model. The performance of the chosen turbulence models is investigated with regard to the available data in the literature including velocity profiles, turbulent kinetic energy and reattachment position in a pipe expansion. In order to further assess the reliability of the turbulence models, an experimental program was conducted by the present authors also at the fluid mechanics laboratory of Menoufiya University. Preliminary measurements, including the surface pressure along the two walls of the expansion pipe and the pressure drop without and with the presence of different arrangements of wall jets produced by symmetrical or asymmetrical fluid cross‐flow injection, are introduced. The results of the present studies demonstrate the superiority of the cubic NLEV k ? ε model in predicting the flow characteristics over the entire domain. The simple low Reynolds number k ? ε model also gives good prediction, especially when the reattachment point is concerned. The evaluation of the reattachment point and the pressure‐loss coefficient is numerically addressed in the paper using the cubic NLEV k ? ε model. The results show that the injection location can control the performance of the pipe‐expansion system. It is concluded that the introduction of flow injection can increase the energy loss in the pipe expansion. The near‐field turbulence structure is also considered in the present study and it is noticed that the turbulence level is strongly affected by the cross‐flow injection and the jet location. Copyright © 2011 John Wiley & Sons, Ltd.     

Turbulence in a skewed three‐dimensional wall‐bounded shear flow: effect of mean vorticity on structure modification

Mean‐flow three‐dimensionalities affect both the turbulence level and the coherent flow structures in wall‐bounded shear flows. A tailor‐made flow configuration was designed to enable a thorough investigation of moderately and severely skewed channel flows. A unidirectional shear‐driven plane Couette flow was skewed by means of an imposed spanwise pressure gradient. Three different cases with 8°, 34°and 52°skewing were simulated numerically and the results compared with data from a purely two‐dimensional plane Couette flow. The resulting three‐dimensional flow field became statistically stationary and homogeneous in the streamwise and spanwise directions while the mean velocity vector V and the mean vorticity vector Ω remained parallel with the walls. Mean flow profiles were presented together with all components of the Reynolds stress tensor. The mean shear rate in the core region gradually increased with increasing skewing whereas the velocity fluctuations were enhanced in the spanwise direction and reduced in the streamwise direction. The Reynolds shear stress is known to be closely related to the coherent flow structures in the near‐wall region. The instantaneous and ensemble‐averaged flow structures were turned by the skewed mean flow. We demonstrated for the medium‐skewed case that the coherent structures should be examined in a coordinate system aligned with V to enable a sound interpretation of 3D effects. The conventional symmetry between Case 1 and Case 2 vortices was broken and Case 1 vortices turned out to be stronger than Case 2. This observation is in conflict with the common understanding on the basis of the spanwise (secondary) mean shear rate. A refined model was proposed to interpret the structure modifications in three‐dimensional wall‐flows. What matters is the orientation of the mean vorticity vector Ω relative to the vortex vorticity vector ω _v, that is, the sign of Ω · ω _v. In the present situation, Ω · ω _v > 0 for the Case 1 vortices causing a strengthening relative to the Case 2 vortices. Copyright © 2011 John Wiley & Sons, Ltd.     

On the relative impact of subgrid‐scale modelling and conjugate heat transfer in LES of hot jets in cross‐flow over cold plates

This work describes numerical simulations of a hot jet in cross‐flow with applications to anti‐ice systems of aircraft engine nacelles. Numerical results are compared with experimental measurements obtained at ONERA to evaluate the performances of LES in this industrial context. The combination of complex geometries requiring unstructured meshes and high Reynolds number does not allow the resolution of boundary layers so that wall models must be employed. In this framework, the relative influence of subgrid‐scale modelling and conjugate heat transfer in LESs of aerothermal flows is evaluated. After a general overview of the transverse jet simulation results, a LES coupled with a heat transfer solver in the walls is used to show that thermal boundary conditions at the wall have more influence on the results than subgrid scale models. Coupling fluid flow and heat transfer in solids simulations is the only method to specify their respective thermal boundary conditions. Copyright © 2011 John Wiley & Sons, Ltd.     

Effects of Reynolds number on physiological‐type pulsatile flows in a pipe with ring‐type constrictions

The effects of Reynolds number on the physiological‐type of laminar pulsatile flow fields within the vicinity of mechanical ring‐type constriction in small pipes were studied numerically. The parameters considered are: the Reynolds number (Re) in the range of 50–1500; Strouhal number (St) in the range of 0.00156–3.98; Womersley number (Nw) from 0.0 to 50.0. The pulsatile flows considered were physiological‐type of simulated flows. Within a pulsating cycle, detailed flow characteristics were studied through the pulsating contours of streamline (ψ), vorticity (Ω), shear stress (τ) and isobar. The relations between the instantaneous flow rate (Q) and instantaneous pressure gradients (dp/dz) are observed to be elliptic. The relations between the instantaneous flow rate (Q) and pressure loss (P_loss) are quadratic. Linear relations were observed between the instantaneous flow rate (Q) and the maximum velocity, maximum vorticity and maximum shear stress. The Reynolds number of the flow in a pulsating cycle was found to have significant effects on the recirculation length and the pressure gradient within the pulsatile flow regime. Copyright © 1999 John Wiley & Sons, Ltd.     

Effect of duct velocity profile and buoyancy‐induced flow on efficiency of transient hydrodynamic removal of a contaminant from a cavity

This paper describes a preliminary numerical analysis of the effect of duct velocity profile and buoyancy‐induced flow generated by the heat source on hydrodynamic removal of contaminants contained in cavities. The process of fluid renewal in a cavity is modelled via a numerical solution of the Navier–Stokes equations coupled with the energy equation for transient flows. The foulant has the same density as the fluid in the duct and the duct velocity profile is considered to be Poiseuille flow and Couette flow, respectively. The results show that the change in Grashof number and duct flow velocity profile causes a dramatic difference in the observed flow patterns and cleaning efficiency. From a cleaning perspective, the results suggest that Couette flow at higher value of Grashof number becomes more effective in further purging of contaminated fluid from a cavity. Copyright © 2004 John Wiley & Sons, Ltd.     

Accelerated convergence of the numerical simulation of incompressible flow in general curvilinear co‐ordinates by discretizations on the double‐staggered grids

The convergence rate of a methodology for solving incompressible flow in general curvilinear co‐ordinates is analyzed. Double‐staggered grids (DSGs), each defined by the same boundaries as the physical domain, are used for discretization. Both grids are MAC quadrilateral meshes with scalar variables (pressure, temperature, etc.) arranged at the center and the Cartesian velocity components at the middle of the sides of the mesh cells. The problem was checked against benchmark solutions of natural convection in a squeezed cavity, heat transfer in concentric horizontal cylindrical annuli, and a hot cylinder in a duct. Poisson's pressure‐correction equations that arise from the SIMPLE‐like procedure are solved by several methods: successive overrelaxation, symmetric overrelaxation, modified incomplete factorization preconditioner, conjugate gradient (CG), and CG with preconditioner. A genetic algorithm was developed to solve problems of numerical optimization of SIMPLE‐like calculation time in a space of iteration numbers and relaxation parameters. The application provides a means of making an unbiased comparison between the DSGs method and the widely used interpolation method. Furthermore, the convergence rate was demonstrated by application to the calculation of natural convection heat transfer in concentric horizontal cylindrical annuli. Calculation times when DSGs were used were 2–10 times shorter than those achieved by interpolation. With the DSGs method, calculation time increases slightly with increasing non‐orthogonality of the grids, whereas an interpolation method calls for very small iteration parameters that lead to unacceptable calculation times. Copyright © 2007 John Wiley & Sons, Ltd.     

Simulation of transient flow in pipelines for computer‐based operations monitoring

Computer models can provide the basis for real‐time monitoring and control of fluid flow in pipelines. Problems of fluid flow in pipelines are mathematically represented by a non‐linear system of coupled partial differential equations. In this paper, several numerical techniques are evaluated with respect to their suitability for the purpose of real‐time monitoring of fluid flow in pipelines. The proposed techniques are evaluated in terms of the L₁, the L₂, and the L_∞ error norms. Moreover, the developed simulators will be compared in terms of their speed of response and settling time which are essential factors for an effective real‐time monitoring scheme. Finally, the selected simulation scheme is further tested under assumed pipeline leak conditions. Copyright © 2004 John Wiley & Sons, Ltd.     

Evaluation of linear and nonlinear convection schemes on multidimensional non‐orthogonal grids with applications to KVLCC2 tanker

Almost all evaluations of convection schemes reported in the literature are conducted using simple problems on uniform orthogonal grids; thus, having limited contribution when solving industrial computational fluid dynamics (CFD), where the grids are usually non‐orthogonal with distortions. Herein, several convection schemes are assessed in uniform and distorted non‐orthogonal grids with emphasis on industrial applications. Linear and nonlinear (TVD) convection schemes are assessed on analytical benchmarks in both uniform and distorted grids. To evaluate the performance of the schemes, four error metrics are used: dissipation, phase and L₁ errors, and the schemes' effective order of accuracy. Qualitative and quantitative deterioration of these error metrics as a function of the grid distortion metrics are investigated, and rigorous verifications are performed. Recommendations for effective use of the convection schemes based on the range of grid aspect ratio (AR), expansion ratio (ER) and skewness (Q) are included. A ship hydrodynamics case is studied, involving a Reynolds averaged Navier–Stokes simulation of a bare‐hull KVLCC2 tanker using linear and nonlinear convection schemes coupled with isotropic and anisotropic Reynolds‐stress (ARS) turbulence models using CFDShip‐Iowa v4. Predictions of local velocities and turbulent quantities from the midships to the nominal wake plane are compared with experimental fluid dynamics (EFD), and rigorous verification and validation analyses for integral forces and moments are performed for 0° and 12° drift angles. Best predictions are observed when coupling a second‐order TVD scheme with the anisotropic turbulence model. Further improvements are observed in terms of prediction of the vortical structures for 30° drift when using TVD2S‐ARS coupled with DES. Copyright © 2009 John Wiley & Sons, Ltd.     

Control of hydrothermal waves in a thermocapillary flow using a gradient‐based control strategy

The problem of controlling the hydrothermal waves in a thermocapillary flow is addressed using a gradient‐based control strategy. The state equations are the two‐dimensional unsteady incompressible Navier–Stokes and energy equations under the Boussinesq approximation. The modeled problem is the ‘open boat’ process of crystal growth, the flow which is driven by Marangoni and buoyancy effects. The control is a spatially and temporally varying heat flux boundary condition at the free surface. The control that minimizes the hydrothermal waves is found using a conjugate gradient method, where the gradient of the objective function with respect to the control variables is obtained from solving a set of adjoint equations. The effectiveness of choices of the parameters governing the control algorithm is examined. Almost complete suppression of the hydrothermal waves is obtained for certain choices of the parameters governing the control algorithm. The numerical issues involved with finding the control using the optimizer are discussed, and the features of the resulting control are analyzed with the goal of understanding how it affects the flow.Copyright © 2012 John Wiley & Sons, Ltd.     

Preserving bounded and conservative solutions of transport in one‐dimensional shallow‐water flow with upwind numerical schemes: Application to fertigation and solute transport in rivers

This work intends to show that conservative upwind schemes based on a separate discretization of the scalar solute transport from the shallow‐water equations are unable to preserve uniform solute profiles in situations of one‐dimensional unsteady subcritical flow. However, the coupled discretization of the system is proved to lead to the correct solution in first‐order approximations. This work is also devoted to show that, when using a coupled discretization, a careful definition of the flux limiter function in second‐order TVD schemes is required in order to preserve uniform solute profiles. The work shows that, in cases of subcritical irregular flow, the coupled discretization is necessary but nevertheless not sufficient to ensure concentration distributions free from oscillations and a method to avoid these oscillations is proposed. Examples of steady and unsteady flows in test cases, river and irrigation are presented. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical analysis of turbulent flow separation in a rectangular duct with a sharp 180‐degree turn by algebraic Reynolds stress model

Turbulent flow in a rectangular duct with a sharp 180‐degree turn is difficult to predict numerically because the flow behavior is influenced by several types of forces, including centrifugal force, pressure‐driven force, and shear stress generated by anisotropic turbulence. In particular, this type of flow is characterized by a large‐scale separated flow, and it is difficult to predict the reattachment point of a separated flow. Numerical analysis has been performed for a turbulent flow in a rectangular duct with a sharp 180‐degree turn using the algebraic Reynolds stress model. A boundary‐fitted coordinate system is introduced as a method for coordinate transformation to set the boundary conditions next to complicated shapes. The calculated results are compared with the experimental data, as measured by a laser‐Doppler anemometer, in order to examine the validity of the proposed numerical method and turbulent model. In addition, the possibility of improving the wall function method in the separated flow region is examined by replacing the log‐law velocity profile for a smooth wall with that for a rough wall. The analysis results indicated that the proposed algebraic Reynolds stress model can be used to reasonably predict the turbulent flow in a rectangular duct with a sharp 180‐degree turn. In particular, the calculated reattachment point of a separated flow, which is difficult to predict in a turbulent flow, agrees well with the experimental results. In addition, the calculation results suggest that the wall function method using the log‐law velocity profile for a rough wall over a separated flow region has some potential for improving the prediction accuracy. Copyright © 2007 John Wiley & Sons, Ltd.     

Exact and numerical simulation of peristaltic flow of a non‐Newtonian fluid with inclined magnetic field in an endoscope

In the present article, we have studied the effects of inclined magnetic field on the peristaltic flow of Jeffrey fluid through the gap between two coaxial inclined tubes. The inner tube is rigid, whereas the outer tube has sinusoidal wave traveling down its wall. The governing equations are simplified using long wave length and low Reynolds number approximations. Exact and numerical solutions have been derived for velocity profile. The expressions for pressure rise and friction force are calculated using numerical integration. Graphical results and trapping phenomenon is presented at the end of the article to see the physical behavior of different parameters. Copyright © 2010 John Wiley & Sons, Ltd.     

Effect of heat transfer on the flow of a second‐grade fluid in divergent/convergent channel

This paper studies the effects of a second‐grade fluid on the flow and heat transfer characteristics in a divergent/convergent channel. The momentum and energy equations are first given in a nondimensional form and then solved analytically using the method of homotopy analysis method. Convergence of derived series solutions is shown. Graphical results for the velocity and the temperatures are presented and discussed for various emerging parameters. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical simulation of 2D lid‐driven cavity flow with CLEARER algorithm on extremely highly skewed grids at high Reynolds numbers

It is crucial to deal with the grid non‐orthogonality effectively in solving the flow in complex geometries, especially at high Reynolds numbers. In this study, the newly proposed Coupled and Linked Equations Algorithm Revised‐ER (CLEARER) algorithm is adopted to solve this problem successfully. In CLEARER algorithm the second relaxation factor is introduced in constructing the contravariant interface velocities, by setting it to a low value. CLEARER algorithm can overcome the severe grid non‐orthogonality and non‐linearity of equations effectively. After the numerical results with CLEARER are validated with the benchmark solutions, this algorithm is used to solve the lid‐driven flow in inclined cavity with inclination angles varying from 10 to 170^°, and Reynolds numbers varying from 5000 to 15 000. The streamlines and the centerline velocity distributions are provided in detail for all cases, which may offer some guidance for the study in this area. Copyright © 2010 John Wiley & Sons, Ltd.