Effect of the adsorption-desorption process intensity on solutal convection near a drop in a horizontal channel

The interaction between gravity convection and Marangoni convection in a horizontal rectangular channel filled with a liquid containing a surfactant and a drop of another liquid is numerically investigated. For large Schmidt numbers the occurring oscillatory regime of solutal convection is analyzed. In the model with a surface phase the effect of the adsorption and desorption processes on the convective flow structure is determined. The corresponding initial and boundary value problem is solved using a difference method.     

Mixed convection boundary layer flow of a viscoelastic fluid over a horizontal circular cylinder

The steady mixed convection boundary layer flow of a viscoelastic fluid over a horizontal circular cylinder in a stream flowing vertically upwards is numerically studied for both cases of heated and cooled cylinders. The governing partial differential equations are transformed into dimensionless forms using an appropriate transformation and then solved numerically using the Keller-box method. The comparison between the solutions obtained and those for a Newtonian fluid is found to be very good. Effects of the mixed convection and elasticity parameters on the skin friction and heat transfer coefficients for a fluid having the Prandtl number equal to one are also discussed. It is found that for some values of the viscoelastic parameter and some negative values of the mixed convection parameter (opposing flow) the boundary layer separates from the cylinder. Heating the cylinder delays separation and can, if the cylinder is warm enough, suppress the separation completely. Similar to the case of a Newtonian fluid, cooling the cylinder brings the separation point nearer to the lower stagnation point. However, for a sufficiently cold cylinder there will not be a boundary layer.     

Low-Prandtl-number Bénard–Marangoni convection in a vertical magnetic field

The effect of a homogeneous magnetic field on surface-tension-driven Bénard convection is studied by means of direct numerical simulations. The flow is computed in a rectangular domain with periodic horizontal boundary conditions and the free-slip condition on the bottom wall using a pseudospectral Fourier–Chebyshev discretization. Deformations of the free surface are neglected. Two- and three-dimensional flows are computed for either vanishing or small Prandtl number, which are typical of liquid metals. The main focus of the paper is on a qualitative comparison of the flow states with the non-magnetic case, and on the effects associated with the possible near-cancellation of the nonlinear and pressure terms in the momentum equations for two-dimensional rolls. In the three-dimensional case, the transition from a stationary hexagonal pattern at the onset of convection to three-dimensional time-dependent convection is explored by a series of simulations at zero Prandtl number.     

Resonance Driven by Horizontal Oscillations of a Vessel Occupied by a Fluid Heated from Above

Unsteady convection driven by periodic horizontal displacements of a rectangular cavity heated from below is investigated numerically. The resonance effects clearly manifest themselves in the frequency dependence of the induced convection flow. These effects appear at a sufficient temperature difference (Grashof number G>10⁴) and a considerable amplitude of the horizontal displacements of the cavity. The dependence of the parameters (resonance frequency and amplitude under resonance conditions) on the problem parameters (Grashof and Prandtl numbers and the cavity length/height ratio) and the method are found. Estimates show that the effects in question can be detected under laboratory conditions.     

Natural convections in conjugated single and double enclosures

The natural convection in single and double conjugated enclosures are numerically investigated. The single and double enclosures are formed by low conductance walls with finite thickness. The outside vertical surfaces of the conducting walls are of the third kind of boundary condition while the top and bottom outside surfaces are adiabatic. The problem studied is characterized by a dominant horizontal temperature gradient and the thermal boundary conditions at the cavity surfaces can not be specified in priori. Numerical results reveal the characteristics in such kind of enclosures and show the importance of the thermal boundary conditions on the natural convection in enclosures. It is also found that the natural convections in the conjugated double enclosures are basically the same, with a major difference in their fluid temperature levels.     

Transient mixed convection flow arising due to thermal and mass diffusion over porous sensor surface inside squeezing horizontal channel

The double diffusion effect on the mixed convection flow over a horizontal porous sensor surface placed inside a horizontal channel is analyzed.With the appropriate transformations,the unsteady equations governing the flow are reduced to non-similar boundary layer equations which are solved numerically for the time-dependent mixed convection parameter.The asymptotic solutions are obtained for small and large values of the time-dependent mixed convection parameter.The results are discussed in terms of the skin friction,the heat transfer coefficient,the mass transfer coefficient,and the velocity,temperature,and concentration profiles for different values of the Prandtl number,the Schmidt number,the squeezing index,and the mixed convection parameter.     

A numerical study of 3D natural convection in a cube: effects of the horizontal thermal boundary conditions

A high-resolution, three-dimensional finite-difference numerical study of natural convection flows of a viscous fluid in a differentially heated cubical box is reported. The vertical sidewalls of the enclosure are maintained at constant temperatures of different values. The other vertical walls (the end walls) are thermally insulated. For the horizontal walls, two kinds of thermal boundary conditions are specified: adiabatic and perfectly conducting. Computations have been performed for an air-filled cavity for Rayleigh numbers of 10⁵ and 10⁶. The specific effects of the horizontal thermal boundary conditions on the flow structure are examined in detail. In the case of conducting walls, heat transfer through the horizontal walls enhances the convective flow activities. The numerically predicted velocity and temperature profiles in the symmetry planes are consistent with previous experimental measurements and computations.     

NUMERICAL INVESTIGATION ON FLOW STABILITY OF RAYLEIGH-B?NARD CONVECTION OF COLD WATER IN A RECTANGULAR CAVITY COOLED AND HEATED SYMMETRICALLY RELATIVE TO THE TEMPERATURE OF DENSITY MAXIMUM

为了解具有密度极值流体瑞利-贝纳德对流特有现象和规律,利用有限容积法对长方体腔内关于密度极值温度对称加热-冷却时冷水瑞利-贝纳德对流的分岔特性进行了三维数值模拟,得到了不同条件下的对流结构型态及其分岔序列,分析了密度极值特性、瑞利数、热边界条件以及宽深比对瑞利-贝纳德对流的影响. 结果表明:具有密度极值冷水瑞利-贝纳德对流系统较常规流体更加稳定,且流动型态及其分岔序列更加复杂;相同瑞利数下多种流型可以稳定共存,各流型在相互转变中存在滞后现象;随着宽深比的增加,流动更易失稳,对流传热能力增强;系统在导热侧壁时比绝热侧壁更加稳定,对流传热能力有所减弱;基于计算结果,采用线性回归方法,得到了热壁传热关联式.     

Modelling far-field entrainment in compressible flows

A procedure that generates steady-states with accurate far-field entrainment is described here. It can be applied to any existing code originally designed for unsteady simulations. Such steady-states are obtained with physical-time damping, introduced through a dual time stepping methodology to effectively eliminate all transient fluctuations without affecting spatial resolution. Far-field entrainment is guaranteed by using properly selected local boundary conditions. Hence, reference profiles with accurate far-field entrainment are generated from the very same code that will eventually use them in unsteady simulations. This new procedure is tested on an unsteady code designed for compressible planar mixing-layer simulations at arbitrary Mach numbers. It is numerically stable for a wide range of Mach numbers and velocity ratios.     

Natural Convection Induced by a Centrifugal Force Field in a Horizontal Annular Porous Layer Saturated with a Binary Fluid

The Darcy Model with the Boussinesq approximation is used to study natural convection in a horizontal annular porous layer filled with a binary fluid, under the influence of a centrifugal force field. Neumann boundary conditions for temperature and concentration are applied on the inner and outer boundary of the enclosure. The governing parameters for the problem are the Rayleigh number, Ra, the Lewis number, Le, the buoyancy ratio, j{\varphi } , the radius ratio of the cavity, R, the normalized porosity, e{\varepsilon } , and parameter a defining double-diffusive convection (a = 0) or Soret induced convection (a = 1). For convection in a thin annular layer (R → 1), analytical solutions for the stream function, temperature and concentration fields are obtained using a concentric flow approximation and an integral form of the energy equation. The critical Rayleigh number for the onset of supercritical convection is predicted explicitly by the present model. Also, results are obtained from the analytical model for finite amplitude convection for which the flow and heat and mass transfer are presented in terms of the governing parameters of the problem. Numerical solutions of the full governing equations are obtained for a wide range of the governing parameters. A good agreement is observed between the analytical model and the numerical simulations.     

Analytical and numerical studies of the boundary slip in the immersed boundary‐thermal lattice Boltzmann method

We analytically and numerically investigate the boundary slip, including the velocity slip and the temperature jump, in immersed boundary‐thermal lattice Boltzmann methods (IB‐TLBMs) with the two‐relaxation‐time collision operator. We derive the theoretical equation for the relaxation parameters considering the effect of the advection velocity on the temperature jump of the IB‐TLBMs. The analytical and numerical solutions demonstrate that the proposed iterative correction methods without the computational cost of the sparse matrix solver reduce the boundary slip and boundary‐value deviation as effectively as the implicit correction method for any relaxation time. Because the commonly used multi‐direct forcing method does not consider the contributions of the body force to the momentum flux, it cannot completely eliminate the boundary slip because of the numerical instability for a long relaxation time. Both types of proposed iterative correction methods are more numerically stable than the implicit correction method. In simulations of flow past a circular cylinder and of natural convection, the present iterative correction methods yield adequate results without the errors of the velocity slip, the temperature jump, and the boundary‐value deviation for any relaxation time parameters and for any number of Lagrangian points per length. The combination of the present methods and the two‐relaxation‐time collision operator is suitable for simulating fluid flow with thermal convection in the multiblock method in which the relaxation time increases in inverse proportion to the grid size.     

Large eddy simulation of a confined square cavity with natural convection based on compressible flow equations

Large eddy simulation of natural convection in a confined square cavity is described. The use of a complex compressible code with an artificial acoustic stiffness correction method, allows the use of higher time steps for a faster time and statistical convergence. We consider a broadly studied experimental case, consisting of a natural convective flow in a confined square cavity, with vertical walls heated at different rates (active walls), set at Ra = 1.58 × 10⁹. Turbulent boundary layers developing on the active walls and a vertical stable stratification characterize the mean flow. It is shown here that the results of this study match the experimental results reported in literature; for instance, mean velocity results. Although results for rms velocity fluctuations are barely over-predicted, the peak region is properly represented, while the greatest disagreements are found in the turbulent heat flow rate (velocity–temperature correlations). Turbulent structures were identified using different visualization methods and statistical studies. The authors found that the boundary layers on the active walls almost reach the fully turbulent regime, tending toward the laminar regime along the horizontal walls.     

Non-darcy double-diffusive mixed convection from heated vertical and horizontal plates in saturated porous media

The non-darcy mixed convection flows from heated vertical and horizontal plates in saturated porous media have been considered using boundary layer approximations. The flows are considered to be driven by multiple buoyancy forces. The similarity solutions for both vertical and horizontal plates have been obtained. The governing equations have been solved numerically using a shooting method. The heat transfer, mass transfer and skin friction are reduced due to inertial forces. Also, they increase with the buoyancy parameter for aiding flow and decrease for the opposing flow. For aiding flow, the heat and mass transfer coefficients are found to approach asymptotically the forced or free convection values as the buoyancy parameter approaches zero or infinity.     

Time‐domain BEM solution of convection–diffusion‐type MHD equations

The two‐dimensional convection–diffusion‐type equations are solved by using the boundary element method (BEM) based on the time‐dependent fundamental solution. The emphasis is given on the solution of magnetohydrodynamic (MHD) duct flow problems with arbitrary wall conductivity. The boundary and time integrals in the BEM formulation are computed numerically assuming constant variations of the unknowns on both the boundary elements and the time intervals. Then, the solution is advanced to the steady‐state iteratively. Thus, it is possible to use quite large time increments and stability problems are not encountered. The time‐domain BEM solution procedure is tested on some convection–diffusion problems and the MHD duct flow problem with insulated walls to establish the validity of the approach. The numerical results for these sample problems compare very well to analytical results. Then, the BEM formulation of the MHD duct flow problem with arbitrary wall conductivity is obtained for the first time in such a way that the equations are solved together with the coupled boundary conditions. The use of time‐dependent fundamental solution enables us to obtain numerical solutions for this problem for the Hartmann number values up to 300 and for several values of conductivity parameter. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical simulation of two-dimensional convection: Role of boundary conditions

The problem of two-dimensional, periodic in the horizontal coordinate, convection of an incompressible fluid heated from below between two horizontal planes is considered. The problem is solved in two formulations: with (stress-)free and hard (no-slip) boundary conditions on the horizontal planes. It is shown that at small supercriticalities the two-dimensional convection calculation leads to more correct results with hard than with free boundary conditions. It is established that the difference between the free and hard conditions is most strongly manifested in the pulsations of the vertical velocity component, whereas the dependence of the Nusselt number and the pulsations of the horizontal velocity component on the boundary conditions is more weakly expressed.     

Quasiperiodic structures in the problem of fluid convection between horizontal planes

A steady flow induced by quasiperiodic temperature distributions at the boundaries is considered with reference to the example of fluid convection in a plane horizontal layer heated from below. A spectral-finite-difference method is employed to obtain convective vortex structures varying quasiperiodically along the horizontal coordinate. The properties of the quasiperiodic solutions, together with their spectra and integral characteristics (Nusselt number) are studied as functions of the Grashof number and the boundary conditions of the problem.     

Mixed convection boundary layer flow over a horizontal circular cylinder with Newtonian heating

The steady mixed convection boundary layer flow over a horizontal circular cylinder, generated by Newtonian heating in which the heat transfer from the surface is proportional to the local surface temperature, is considered in this study. The governing boundary layer equations are first transformed into a system of non-dimensional equations via the non-dimensional variables, and then into non-similar equations before they are solved numerically using a numerical scheme known as the Keller-box method. Numerical solutions are obtained for the skin friction coefficient Re ^1/2 C _f and the local wall temperature θ _w (x) as well as the velocity and temperature profiles with two parameters, namely the mixed convection parameter λ and the Prandtl number Pr.     

Effect of Surface Corrugation on Convection in a Three-Dimensional Finite Box of Fluid Saturated Porous Material

The problem of finite amplitude thermal convection in a three-dimensional finite box of fluid saturated porous material is investigated, when the lower boundary of the fluid is corrugated. The nonlinear problem of three-dimensional convection in the box for the values of the Rayleigh number close to the classical critical value and for small values of the amplitude of the corrugations is solved by a perturbation technique. The preferred mode of convection is determined by stability analysis. In the absence of corrugation three-dimensional modes of convection can be either stable or unstable depending on the values of the aspect ratios of the box, while two-dimensional rolls are always stable, provided that the box aspect ratios allow the existence of such modes of convection. In the presence of boundary corrugation with the appropriate form, different three-dimensional or two-dimensional modes of corrugation can be stable or unstable. For a rough boundary with local roughness sites, the location, size, and number of the roughness elements plus the wave numbers of the convection modes and the box aspect ratios can all play a role leading to either stable or unstable particular three- or two-dimensional flow patterns. For a wavy boundary, resonant wave-vector excitation can lead to the preference of stable two- or three-dimensional flow patterns whose wave vectors are in a subset of those due to the wavy boundary, while nonresonant wave-vector excitation can lead to the preference of stable flow patterns whose wave vectors are not generally in a subset of those due to the wavy boundary. Heat transported by convection can either be enhanced or be reduced by certain proper forms of the corrugations and by appropriate values of the box aspect ratios. Due to the surface corrugation highly subcritical modes of convection are stable, while highly supercritical modes of convection are unstable. Received 24 July 1998 and accepted 11 April 1999     

Laminar natural convection of Cu-water nanofluid in concentric annuli with radial fins attached to the inner cylinder

Laminar natural convection of Cu-water nano-fluid between two horizontal concentric cylinders with radial fins attached to the inner cylinder is studied numerically. The inner and outer cylinders are maintained at constant temperature. The governing equations in the polar two-dimensional space with the respective boundary conditions are solved using the finite volume method. The hybrid-scheme is used to discretize the convection terms. In order to couple the velocity field and the pressure in the momentum equations, the well known semi-implicit method for pressure linked equation reformed algorithm is adopted. Using the developed code, a parametric study is undertaken, and the effects of the Rayleigh number, Number of fins, length of the fins and the volume fraction of nano-particles on the fluid flow and heat transfer inside the annuli are investigated. In this study, two cases with different number of fins are considered. It is observed from the results that the average Nusselt number increases with increasing both the Rayleigh number and the volume fraction of the nano-particles. Moreover, the average Nusselt number decreases by increasing the fins’ length and the number of fins. Heat transfer rate increases by increasing the fins’ length at all Rayleigh numbers, but it increases by increasing the number of fins at high Rayleigh numbers.     

Mixed convection boundary layer flow past a wedge with permeable walls

The effects of suction/injection on steady laminar mixed convection boundary layer flow over a permeable horizontal surface of a wedge in a viscous and incompressible fluid is considered in this paper. The similarity solutions of the governing boundary layer equations are obtained for some values of the suction/injection parameter f ₀, the constant exponent m of the wall temperature as well as the mixed convection parameter λ. The resulting system of nonlinear ordinary differential equations is solved numerically for both assisting and opposing flow regimes using an implicit finite-difference scheme known as the Keller-box method. Numerical results for the reduced skin friction coefficient, the local Nusselt number, and the velocity and temperature profiles are obtained for various values of parameters considered. Dual solutions are found to exist for the case of opposing flow.