Stability of a nonuniformly heated fluid in a porous horizontal layer

We investigate the stability of a nonuniformly heated fluid in the gravitational field in a plane horizontal porous layer through which vertical forced motion is effected. A similar system was studied in [1, 2]. In the present paper, the nonuniformity of the permeability of the porous layer with respect to the depth and the dependence of the viscosity of the saturating fluid on the temperature are taken into account in addition. As a result of the application of the linear stability theory, an eigenvalue problem arises, which is solved numerically. A family of curves representing the dependence of the critical modified Rayleigh number Ra _k ^⋆ on the injection parameter (the Péclet number Pe) for different degrees of inhomogeneity of the permeability and the viscosity is obtained. It is found that although Pe=0 corresponds to Ra _k ^⋆ for uniform permeability and viscosity and the stability increases monotonically as Pe increases, the presence of nonuniformity of the permeability or the viscosity leads to the appearance of a stability minimum in the region Pe≈1, while under the simultaneous influence of these two factors, the minimum is shifted into the region Pe≈2. The results of the paper can be used, for example, in the investigation of heat transfer in the case of forced fluid motion in the fissures of a permeable rock mass, when, in the case of pumping through a horizontal fissure, the fluid penetrates vertically across its permeable walls into the stratum. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 3–7, November–December, 1986.     

Convective wall plume in power-law fluid: Second-order correction for the adiabatic wall

This paper considers the steady-state free convection flow arising from an infinitely long horizontal line source of heat embedded in the base of a vertical adiabatic surface when the ambient fluid is a non-Newtonian fluid for moderately large values of the generalized Grashof numbers by the method of matched asymptotic expansions. In particular, the second-order corrections to account for the non-boundary layer effects have been predicted. A family of numerical solutions for the power-law fluid behavior indexn ranging from 0.4 to 2.0 and for the Prandtl numberPr=10 and 100 are reported.     

Mixed convection from a cylinder oscillating vertically in a quiescent fluid

 The problem of heat convection from a vertically oscillating cylinder in a quiescent fluid is investigated. The governing equations of motion and energy are solved numerically in a non-inertial frame of references to determine the flow field and heat transfer characteristics under different conditions. The main dominating parameters are Keulegan–Carpenter number, KC, frequency parameter, β, Grashof number, Gr and Prandtl number, Pr. The ranges considered for these parameters are KC ≤ 10, β≤40 and Gr ≤ 10⁵ while Prandtl number is kept constant. The study revealed that the effect of amplitude and frequency of oscillation on heat transfer is strongly influenced by the Grashof number range. In the forced convection regime (Gr = 0), the increase of KC creates extensive vortex motion at all cylinder positions that leads to a significant increase in heat transfer. A similar trend, but with a lesser extent, is also observed for the increase of β. However, at high Grashof numbers, the effect of oscillation on heat convection is only significant at large values of KC. Received on 5 June 2000 / Published online: 29 November 2001     

Unsteady natural convection in an enclosure with vertical wavy walls

In this paper, unsteady heat transfer and fluid flow characteristics in an enclosure are investigated. The enclosure consists of two vertical wavy and two horizontal straight walls. The top and the bottom walls are considered adiabatic. Two wavy walls are kept isothermal and their boundaries are approximated by a cosine function. Governing equations including continuity, momentum and energy were discretized using the finite-volume method and solved by SIMPLE method in curvilinear coordinate. Simulation was carried out for a range of Grashof number Gr = 10³–10⁶, Prandtl number Pr = 0.5–4.0, wave ratio A (defined by amplitude/wavelength) 0.0–0.35 and aspect ratio W (defined by average width/wavelength) 0.5–1.0. Streamlines and isothermal lines are presented to corresponding flow and thermal fields. Local and average Nusselt number distributions are presented. The obtained results are in good agreement with available numerical and experimental data.     

Natural convection in an infinite square under low-gravity conditions

In order to study natural convection effects on fluid flows under low-gravity in space, we have expanded variables into a power series of Grashof number by using perturbation theory to reduce the Navier-Stokes equations to the Poisson equation for temperature T and biharmonic equation for stream function φ. Suppose that a square infinite closed cylinder horizontally imposes a specified temperature of linear distribution on the boundaries, we investigate the two dimensional steady flows in detail. The results for stream function φ, velocity u and temperature T are gained. The analysis of the influences of some parameters such as Grashof number G_r and Prandtl number P_r on the fluid motion lead to several interesting conclusions. At last, we make a comparison between two results, one from approximate equations, the other from the original version. It shows that the approximate theory correctly simplifies the physical problem, so that we can expect the theory will be applied to unsteady or three-dimensi     

Convective magnetohydrodynamic two fluid flow and heat transfer in an inclined channel

 The problem of fully developed free convection two fluid magnetohydrodynamic flow in an inclined channel is investigated. The governing momentum and energy equations are coupled and highly nonlinear due to dissipation terms, solutions are found employing perturbation technique for small values of Pr · Ec (=ɛ) the product of Prandtl number and Eckert number. Effects of Grashof number, Hartmann number, inclination angle, the ratios of electrical conductivities, viscosities and heights of two fluids on the flow are explored. It is observed that the flow can be controlled effectively by suitable adjustment of the values for the ratios of heights, electrical conductivities and viscosities of the two fluids. Received on 10 December 1999     

Effect of rotation on the formation of longitudinal vortices in mixed convection flow over a flat plate

This paper presents a numerical study of the effect of rotation on the formation of longitudinal vortices in mixed convection flow over a flat plate. The criterion on the position of marking the onset of longitudinal vortices is defined in this paper. The onset position characterized by the Goertler number G _δ depends on the Grashof number, the rotation number Ro, the Prandtl number Pr and the wave number. The results show that negative rotation stabilizes the boundary layer flow on the surface. On the contrary, positive rotation destabilizes the flow. The numerical data are compared with the experimental results.     

Effect of water density inversion on the plane-parallel flow and heat transfer in a channel of constant width

The results of a numerical study of the effect of cold-water density inversion (Prandtl number Pr=11.59) on the flow and heat transfer in a horizontal plane-parallel channel with isothermal top and bottom walls are presented. The calculations were performed for the Grashof number Gr=3·10⁴, the Reynolds number Re=10, and the channel segment length-to-height ratiol/d=40. The wall temperature was so varied that the temperature difference between the top and bottom walls remained constant. Surgut. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 72–78, January–February, 2000.     

Heat transfer in a hydromagnetic flow over a stretching sheet

Exact solutions are obtained for the heat transfer in an electrically conducting fluid past a stretching sheet subjected to the thermal boundary with either a prescribed temperature or a prescribed heat flux in the presence of a transverse magnetic field. The solutions for the heat transfer characteristics are evaluated numerically for different parameters, such as the magnetic parameterN, the Prandtl numberPr, the surface temperature indexs, and the surface heat flux indexd. It is observed that for the prescribed surface temperature case the fluid temperature increases due to the existance of the magnetic field, and decreases as the Prandtl number or the surface temperature index increases; for the prescribed surface heat flux case, the surface temperature decreases as the Prandtl number of the surface heat flux index increases, and the magnetic parameter decreases. In addition, varying the prescribed surface temperature indexs affects the mechanism of heat transfer.     

Determining modes and Grashof number in 2D turbulence: a numerical case study

We study how the number of numerically determining modes in the Navier–Stokes equations depends on the Grashof number. Consider the two-dimensional incompressible Navier–Stokes equations in a periodic domain with a fixed time-independent forcing function. We increase the Grashof number by rescaling the forcing and observe through numerical computation that the number of numerically determining modes stabilizes at some finite value as the Grashof number increases. This unexpected result implies that our theoretical understanding of continuous data assimilation is incomplete until an analytic proof which makes use of the non-linear term in the Navier–Stokes equations is found.      

Charging and coagulation of water aerosols with a small admixture of highly radioactive droplets

The mechanism of electrocoagulation of water aerosols with a small admixture of highly radioactive droplets is examined. A corresponding mathematical model describing the processes of ionization, electrization and coagulation of radioactive water aerosols is developed. The time dependence of the ion concentration and the charge and concentration of the nonradioactive droplets and of the charge and radius of the radioactive droplets is numerically investigated for a number of typical aerosols. It is shown that the electrocoagulation process may lead to an increase in the radius of the droplets from 5–10 to 30–40μm in ≃10⁴ sec ≃3 h and, consequently, may play a significant part in the development of aerosols with a droplet radius of up to 20μm, when gravitational coagulation is unimportant. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.2, pp. 90–96, March–April, 1994.     

The Linearized Crocco Equation

In this paper, we study the existence and uniqueness of a degenerate parabolic equation, with nonhomogeneous boundary conditions, coming from the linearization of the Crocco equation [12]. The Crocco equation is a nonlinear degenerate parabolic equation obtained from the Prandtl equations with the so-called Crocco transformation. The linearized Crocco equation plays a major role in stabilization problems of fluid flows described by the Prandtl equations [5]. To study the infinitesimal generator associated with the adjoint linearized Crocco equation – with homogeneous boundary conditions – we first study degenerate parabolic equations in which the x-variable plays the role of a time variable. This equation is doubly degenerate: the coefficient in front of ∂_x vanishes on a part of the boundary, and the coefficient of the elliptic operator vanishes in another part of the boundary. This makes very delicate the proof of uniqueness of solution. To overcome this difficulty, a uniqueness result is first obtained for an equation in which the elliptic operator is symmetric, and it is next extended to the original equation by combining an iterative process and a fixed point argument (see Th. 4.9). This kind of argument is also used to prove estimates, which cannot be obtained in a classical way.     

Boundary effects on the Bénard-Marangoni instability under an electric field

This article theoretically studies the Bénard-Marangoni instability problem for a liquid layer with a free upper surface, which is heated from below by a heating coil through a solid plate in ana.c. electric field. The boundary effects of the solid plate, which include its thermal conductivity, electric conductivity and thickness, have great influences on the onset of convective instability in the liquid layer. The stability analysis in this study is based on the linear stability theory. The eigenvalue equations obtained from the analysis are solved by using the fourth order Runge-Kutta-Gill's method with the shooting technique. The results indicate that the solid plate with a higher thermal or electric conductivity and a bigger thickness tends to stabilize the system. It is also found that the critical Rayleigh numberR _c, the critical Marangoni numberM _c, and the criticala.c. Rayleigh numberE _ac become smaller as the intensity of thea.c. electric field increases.     

The influence of the continuous phase Pe numbers on thermal wake phenomenon

The reverse of the transfer direction in the unsteady conjugate heat transfer between a spherical particle and a surrounding fluid flow has been analysed. The aspect this work is focused on is the influence of the continuous phase convection on the occurrence and development of this phenomenon. The energy equations are solved by the ADI finite difference method. The range of the Pe numbers investigated is between 0 and 10. The ratios of the thermal conductivity and volume heat capacity between the particle and its ambient flow belong to the interval 0.01–100. It was found that, in creeping flow, the thermal wake occurs at Pe=0.690·10⁻³. Increasing the Pe number up to 1 the dimension of thermal wake increases. For Pe>1, the increase in Pe decreases thermal wake. Received on 13 January 1998     

Mixed convection in a vertical channel with discrete heat sources at the wall

The mixed convection in a vertical plane-parallel channel with two heat sources of finite dimensions located at the wall is analyzed on the basis of a two-dimensional numerical simulation. The effect of the distance between the heat sources on the flow pattern and the temperature field is studied. Calculations are performed on the Grashof and Reynolds number ranges from 0–10⁵ and 0–10, respectively, at a Prandtl number of 0.7. The mathematical model is based on the time-dependent Navier-Stokes equations in the Boussinesq approximation. The problem is solved by the finite element method.     

Thermal instability in natural convection flow over a boundary layer subject to external electrical and magnetic fields

This paper presents a numerical study of external electrical and magnetic effects on the formation of longitudinal vortices in natural convection flow over a heated horizontal plate. The criterion on the position marking on the onset of longitudinal vortices is defined in the present paper. The onset position characterized by the Grashof number depends on the Prandtl number, the wave number, the electric field parameter, and the magnetic field parameter. The flow is found more stable as the value of the magnetic field parameter increases. Moreover, the stabilizing effect is also found on the flow when the positive electric field parameter Ec is applied. The results of the present numerical prediction show reasonable agreement with the experimental data with zero magnetic field and electric field parameters in literature.     

Forced convection of gaseous slip-flow in porous micro-channels under Local Thermal Non-Equilibrium conditions

Steady laminar forced convection gaseous slip-flow through parallel-plates micro-channel filled with porous medium under Local Thermal Non-Equilibrium (LTNE) condition is studied numerically. We consider incompressible Newtonian gas flow, which is hydrodynamically fully developed while thermally is developing. The Darcy–Brinkman–Forchheimer model embedded in the Navier–Stokes equations is used to model the flow within the porous domain. The present study reports the effect of several operating parameters on velocity slip and temperature jump at the wall. Mainly, the current study demonstrates the effects of: Knudsen number (Kn), Darcy number (Da), Forchheimer number (Γ), Peclet number (Pe), Biot number (Bi), and effective thermal conductivity ratio (K _R) on velocity slip and temperature jump at the wall. Results are given in terms of skin friction (C _f Re ^*) and Nusselt number (Nu). It is found that the skin friction: (1) increases as Darcy number increases; (2) decreases as Forchheimer number or Knudsen number increases. Heat transfer is found to (1) decreases as the Knudsen number, Forchheimer number, or K _R increases; (2) increases as the Peclet number, Darcy number, or Biot number increases.     

Transient natural convection of low-Prandtl-number fluids in a differentially heated cavity

The transient convective motion in a two-dimensional square cavity driven by a temperature gradient is analysed. The cavity is filled with a low-Prandtl-number fluid and the vertical walls are maintained at constant but different temperatures, while the horizontal boundaries are adiabatic. A control volume approach with a staggered grid is employed to formulate the finite difference equations. Numerically accurate solutions are obtained for Prandtl numbers of 0·001, 0·005 and 0·01 and for Grashof numbers up to 1 × 10⁷. It was found that the flow field exhibits periodic oscillation at the critical Grashof numbers, which are dependent on the Prandtl number. As the Prandtl number is decreased, the critical Grashof number and the frequency of oscillation decrease. Prior to the oscillatory flow, steady state solutions with an oscillatory transient period were predicted. In addition to the main circulation, four weak circulations were predicted at the corners of the cavity.     

Numerical study of transient laminar natural convection heat transfer over a sphere subjected to a constant heat flux

Transient laminar natural convection over a sphere which is subjected to a constant heat flux has been studied numerically for high Grashof numbers (10⁵ ≤ Gr ≤ 10⁹) and a wide range of Prandtl numbers (Pr = 0.02, 0.7, 7, and 100). A plume with a mushroom-shaped cap forms above the sphere and drifts upward continuously with time. The size and the level of temperature of the transient cap and plume stem decrease with increasing Gr and Pr. Flow separation and an associated vortex may appear in the wake of the sphere depending on the magnitude of Gr and Pr. A recirculation vortex which appears and grows until “steady state” is attained was found only for the very high Grashof numbers (10⁵ ≤ Gr ≤ 10⁹) and the lowest Prandtl number considered (Pr = 0.02). The appearance and subsequent disappearance of a vortex was observed for Gr = 10⁹ and Pr = 0.7. Over the lower hemisphere, the thickness of both the hydrodynamic (δ_H) and the thermal (δ_T) boundary layers remain nearly constant and the sphere surface is nearly isothermal. The surface temperature presents a local maximum in the wake of the sphere whenever a vortex is established in the wake of the sphere. The surface pressure recovery in the wake of the sphere increases with decreasing Pr and with increasing Gr. For very small Pr, unlike forced convection, the ratio δ_T/δ_H remains close to unity. The results are in good agreement with experimental data and in excellent agreement with numerical results available in the literature. A correlation has also been presented for the overall Nusselt number as a function of Gr and Pr.     

Mass transfer effects on flow past an impulsively started infinite vertical plate with constant mass flux — an exact solution

An exact solution to the problem of flow past an impulsively started infinite vertical plate in the presence of a foreign mass and constant mass flux at the plate is presented by the Laplace-transform technique. The velocity, the temperature and the concentration profiles are shown on graphs. The skin-friction and the Sherwood number are also shown on graphs. The effects of different parameters likeG (the Grashof number),Gc (the modified Grashof number),Pr (the Prandtl number) andSc (the Schmidt number) are discussed.