Magnetoconvection of two-immiscible fluids in vertical enclosure

Two-fluid magnetoconvection in a vertical enclosure consisting of two regions, one electrically conducting and another electrically non-conducting is investigated analytically. The transport properties of both fluids are assumed constant. The resulting coupled non-linear equations are solved using regular perturbation method. It is found that the effect of negative electric load parameter E is to aid the flow while the effect of positive E is to oppose the flow as compared to the case when E=0. The effects of various physical parameters such as Grashof number, viscosity ratio and width ratio are presented for heat absorption and heat generation case.     

Transient three-dimensional flow in an enclosure with a hot spot on a vertical wall

A numerical method for the solution of the vector potential/vorticity vector formulation of the transient, fully three-dimensional Navier-Stokes energy and continuity equations has been applied to simulate the development of natural convective flow within a ‘box’ after a sudden temperature change on a vertical portion of the wall. Only one cavity size has been considered, this having a vertical height of three times its width and a horizontal length of six times its width. A single heated rectangular hot spot or ‘element’ on an otherwise adiabatic wall is centred between the vertical end walls. The opposite vertical wall is held at the intial fluid temperature, and all other walls are assumed to be adiabatic. Fluid properties have been assumed constant except for the density change with temperature that gives rise to the buoyancy force. The numerical method is an underrelaxation Gauss-Seidel method using finite differencing at each time step. Solutions have been obtained for a Prandtl number of 0.71, for Rayleigh numbers, based on the width, of between 0 and 100000 and for a number of heated element locations and sizes.     

Numerical analysis of natural convection around a vertical channel in a rectangular enclosure

Numerical computations were performed for the heat transfer and fluid flow characteristics of an internal vertical channel composed by a pair of parallel plates situated in a rectangular enclosure, with the inner plates and the bounding wall of the enclosure maintained at uniform but different temperatures. Natural convection occurred in the air which occupied the enclosure space. The plates were symmetrically arranged. The dimensionless channel widthS was varied parametrically. The Rayleigh numbers ranged from 10² to 10⁷. Static bifurcation was found in this configuration. The bifurcation is related to the flow pattern transition from single-vortex structure to double-vortex structure or vice versa. Comparison with the empirical correlations obtained for a vertical plate and a channel in an infinite space showed that the heat transfer process of the plates and the channel was deteriorated by the existence of the enclosure.     

Analytical and numerical study of double diffusive convection in a vertical enclosure

This paper presents an analytical and numerical study of natural convection of a double-diffusive fluid contained in a rectangular slot subject to uniform heat and mass fluxes along the vertical sides. Governing parameters of the problem under study are the thermal Rayleigh number, Ra _T ; buoyancy ratio, N; Lewis number, Le; Prandtl number, Pr and aspect ratio of the cavity, A. In the first part of the analytical study a scale analysis is applied to the two extreme cases of heat-transfer and mass-transfer-driven flows. In the second part, an analytical solution, based on the parallel flow approximation, is reported for tall enclosures (A?1). Solutions for the flow fields, temperature and concentration distributions and Nusselt and Sherwood numbers are obtained in terms of the governing parameters of the problem. In the limits of heat-driven and solute-driven flows a good agreement is obtained between the prediction of the scale analysis and those of the analytical solution. The numerical solutions are based on the complete governing equations for two-dimensional flows, and cover the range 1≤Ra _T ≤10⁷, 0≤N≤10⁵, 10^-3≤Le≤10³, 1≤A≤20 and Pr=7. A good agreement is found between the analytical predictions and the numerical simulation.     

Natural convection in a square enclosure with an internal isolated vertical plate

Numerical computations were performed for the average Nusselt number at an internal vertical plate situated in a square enclosure, with the inner plate and the bounding wall of the enclosure maintained at uniform but different temperatures. Natural convection occurred in the air which occupied the enclosure space. The position of the inner vertical plate within the enclosure was varied parametrically. The plate height-cavity height ratio was 0.513. For narrow distance between the inner plate and the bounding wall the inner plate Nusselt number was enhanced. Aside from this, the plate average Nusselt number was remarkably insensitive to the plate position. The effect of the Rayleigh number on the velocity and temperature fields and local Nusselt numbers are also discussed. The agreement between the predicted flow pattern forRa=1.1×10⁶ and the flow visualization result was reasonably good.

---

Natürliche Konvektion in einem quadratischen Horizontalschacht, der eine freistehende, senkrechte Platte enthält

Zusammenfassung Eine numerische Untersuchung liefert mittlere Nußelt-Zahlen an einer, in einem quadratischen Horizontalschacht freistehenden, senkrechten Platte, wobei deren Temperatur und die der umgebenden Wände jeweils konstant gehalten werden. Im Luftraum dazwischen stellte sich freie Konvektion ein. Die Position der Platte war veränderlich, ihre Höhe blieb mit 51.3% der Schachthöhe konstant. Rückte die Platte nahe an eine Schachtwand, so erhöhte sich die Nußelt-Zahl auf der dieser zugewandten Seite, während die Gesamt-Nußelt-Zahl bezüglich der Platte fast konstant bleibt. Es wird auch der Einfluß der Rayleigh-Zahl auf das Geschwindigkeitsund Temperaturfeld diskutiert. BeiRa=1.1·10⁶ stimmten die Ergebnisse aus der Berechnung gut mit den experimentellen Befunden einer Strömungsvisualisation überein.

---

Nomenclature a distance between vertical plate and side-wall of enclosure thermal diffusivity (in definition ofu _r) - b distance between vertical plate and bottom of enclosure - g gravitational acceleration - G characteristic flow rate - H height of vertical plate - k thermal conductivity - k _f fluid thermal conductivity - K relative thermal conductivity,k/k _f - L width of square enclosure - M _res mass residual - Nu local Nusselt number - Nu _m average Nusselt number - Nu _L local Nusselt number of left side of vertical plate - Nu _R local Nusselt number of right side of vertical plate - Nu _B local Nusselt number of bottom side of vertical plate - Nu _T local Nusselt number of top side of vertical plate - p effective pressure - P dimensionless pressure,P=p/[(Ra Pr)(a/H)²] - Pr Prandtl number - Ra Rayleigh number,Ra=gTH ³ Pr/ ² - T temperature - T _i temperature of internal plate - T _o temperature of enclosure surface - u, v velocity components inx-, y-direction - U, V dimensionless velocities,U=u/u _r, V=v/u_r - u _r reference velocity,u _r=(Ra Pr)^1/2/(a/H) - X, {iyY} dimensionless coordinates,X=x/H, Y=y/H Greek symbols heat transfer coefficient - volume expansion coefficient - thickness of plate - kinematic viscosity - density - dimensionless temperature, (T _i–T)/(T _i–T_o)     

Experimental and numerical analysis of thermo-magnetic convection in a vertical annular enclosure

In the present paper an experimental and numerical analysis of a thermo-magnetic convective flow of paramagnetic fluid in an annular enclosure with a round rod core and a cylindrical outer wall is presented. It is complemented by an experimental analysis of natural convection depending on the inclination angle to show the stability of the present configuration. Convection in an annulus between two vertical co-axial cylinders resulting from gravitational and magnetic environments has been investigated. A strong magnetic field can be an alternative to heat transfer enhancement. The effect of the magnetic field on the convection of the paramagnetic fluid in the annular vessel in various positions was compared. The numerical analysis was done based on the continuity, momentum and energy equations. A term related to the magnetic buoyancy force was added to the momentum equation. The distributions of Nusselt number present minima in two positions of the enclosure, which depends on the reciprocal relationship between the gravitational and magnetising forces.     

Entropy generation in a square enclosure with partial heating from a vertical lateral wall

The entropy generation during the transient laminar natural convection in a square enclosure that is partially heated from a vertical lateral wall is numerically investigated. The active sites referring to the main irreversibility locations are determined. The Boussinesq approximation is used in the natural convection modelling. The effects of Prandtl (Pr) and Rayleigh (Ra) number combinations on the entropy generation are investigated. The study is restricted to the fluids of Prandtl number from 0.01 to 1.0, and Rayleigh numbers in the range of 10²–10⁸. It is found that the upper corner of the heated part of the side wall is the active site where the entropy generation initiates due to irreversibilities representing the energy loss.     

Solution of three-dimensional natural convection about a vertical square rod in a cylindrical enclosure

The problem of three-dimensional laminar natural convection in a vertical enclosure with an inner square rod is treated by a numerical method in boundary-fitted co-ordinates. The inner and outer cylinders are heated and cooled, respectively, to maintain different constant surface temperatures. The horizontal enclosure surfaces are maintained at adiabatic conditions. The Prandtl number is that of air, 0·703, and the Rayleigh numbers span the conduction, transition and boundary layer regimes of flow. The radius ratio is 1, and the aspect ratio (cylinder length divided by maximum annular gap) is 1. The results of the study provide data useful in the design and performance assessment of nuclear reactor spent fuel shipping casks.     

Natural convection in a vertical rectangular enclosure with localized heating and cooling zones

Experimental and numerical studies of natural convection in a single phase, closed thermosyphon were carried out using a vertical, rectangular enclosure model. Only one vertical plate plays the role of heat transfer surface having 100 mm height and 100 mm width, and others act as the adiabatic wall made of transparent plexi-glass. The heat transfer surface is separated into three horizontal zones with an equal height; top 1/3 and bottom 1/3 of the surface are cooling and heating zones, respectively and intermediate section is an adiabatic zone. Water is used as the working fluid. Variable parameters are distance D between the heat transfer surface and an adiabatic plate opposite to the heat transfer plate, and temperature difference ΔT between heating and cooling zones. By changing both D and ΔT, three regimes of the natural convection flow; quasi-two-dimensional steady, three-dimensional steady and unsteady flows are observed by means of thermo-sensitive liquid crystal powder and numerically simulated very well by solving a set of governing equations. Received on 17 January 2000     

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.     

Analytical and numerical study of buoyancy-driven convection in a vertical enclosure filled with nanofluids

This paper presents an analytical and numerical study of natural convection of nanofluids contained in a rectangular enclosure subject to uniform heat flux along the vertical sides. Governing parameters of the problem under study are the thermal Rayleigh number Ra, the Prandtl number Pr, the aspect ratio of the cavity A and the solid volume fraction of nanoparticles, Φ. Three types of nanoparticles are taken into consideration: Cu, Al₂O₃ and TiO₂. Various models are used for calculating the effective viscosity and thermal conductivity of nanofluids. In the first part of the analytical study, a scale analysis is made for the boundary layer regime situation. In the second part, an analytical solution based on the parallel flow approximation is reported for tall enclosures (A ≫ 1). In the boundary layer regime a good agreement is obtained between the predictions of the scale analysis and those of the analytical solution. Solutions for the flow fields, temperature distributions and Nusselt numbers are obtained explicitly in terms of the governing parameters of the problem. A numerical study of the same phenomenon, obtained by solving the complete system of the governing equations, is also conducted. A good agreement is found between the analytical predictions and the numerical simulations.     

A numerical study of the melting of phase change material heated from a vertical wall of a rectangular enclosure

This article presents the first research effort of our group to formulate, implement and validate a numerical method in order to optimise the design of solar passive walls involving phase change materials (PCMs). The fusion of ice, gallium and the commercially available PCM 27 (hydrated salt), engineered by Cristopia and later embedded within an experiment unit, was studied. Comparisons against other prediction methods and experimental data for the fusion of gallium were carried out with good agreement of the solutions. The proposed enthalpy-based method is found to be excellent to predict the fusion, but still fails to reproduce adequately the exact solidification pattern measured for the PCM 27. Further research is going on to improve the model.     

A note on the free convection boundary-layer flow over a vertical surface with spatial internal heat generation

Meccanica - The model proposed by Postelnicu et al. [17, 18] for the natural convection boundary-layer flow on a vertical surface in a porous medium driven by spatially-dependent localised internal...     

Onset of unsteady axi-symmetric laminar natural convection in a vertical cylindrical enclosure heated at the wall

In the present study laminar transition to oscillatory convection of fluids having different Prandtl numbers in a laterally heated vertical cylindrical enclosure for different aspect ratios (melt height to crucible radius) of 2–4 is investigated numerically for 0.01 ≤ Pr ≤ 10. Numerical solution to two-dimensional axisymmetric transient Navier Stokes equations and energy equation were solved by finite volume method using SIMPLE algorithm. Numerical results illustrate that there exists a critical Rayleigh number for each Prandtl number beyond which sustained laminar oscillatory flow sets in. The oscillatory regime was characterised by the oscillation of the average kinetic energy and average thermal energy of the melt. For a given aspect ratio, critical Rayleigh number increases with Pr upto 1 and then flattens. It was observed that for low Prandtl number fluids, Pr < 1.0, critical Rayleigh number is found to increase with increase in aspect ratio while for high Prandtl number fluids, Pr ≥ 1.0, it is found to decrease with increase in aspect ratio. The influence of aspect ratio on the transient behaviour of the melt volume below and above the critical Rayleigh number was studied.     

A note on the diffraction of an obliquely incident surface wave by a partially immersed fixed vertical barrier

The problem of the diffraction of surface waves, obliquely incident on a partially immersed fixed vertical barrier in deep water, is solved approximately by reducing it to the solution of an integral equation, for small angle of incidence of the incident wave. The corrections to the reflection and transmission coefficients over their normal incidence values for small angle of incidence are obtained and presented graphically for some intermediate values of wave numbers.     

Numerical simulation of laminar natural convection in a laterally heated vertical cylindrical enclosure: application to crystal growth

Laminar natural convection has been studied in a laterally heated vertical cylindrical enclosure with a free insulated surface and a centrally located constant temperature wall at the top. These conditions are a simplification of the conditions existing in a Czochralski crystal pulling system. The laminar, axisymmetric flow of a Newtonian, constant physical properties fluid under Boussinesq’s approximation has been considered. Governing equations in primitive variable form are solved numerically by control volume method. SIMPLE algorithm due to Patankar has been used for the numerical simulation. The effects of the constant wall heat flux boundary condition at the side wall have been investigated whereas the bottom wall is considered to be insulated. Streamlines and isotherms are presented for various Rayleigh numbers and Prandtl numbers. Heat flux vectors through the melt are plotted for selected cases. The axial velocity and temperature variations at different horizontal sections of the crucible have been presented graphically to explain the transport processes inside the crucible. It has been observed that in case of low Pr and high Ra, flow separation occurs at the vertical wall of the crucible which leads to an oscillatory flow as Ra increases. The investigation has been extended to the oscillatory regime of flow in the zone of supercritical Rayleigh numbers and some unsteady results are also presented. Finally a heat transfer correlation has been developed for steady-state case.     

On cooling behavior of a vertical plate in a phase change material/water composite enclosure under pulsating heat load

This paper presents a numerical and experimental study concerning cooling characteristics of a pulsating heated vertical plate sandwiched between a substrate of phase change material (PCM) and an enclosure of water, forming a composite vertical rectangular enclosure. The vertical plate is assumed to have a uniform pulsating (on/off) volumetric heat source. The PCM considered in the present study is n-Octadecane. In the finite-difference simulation, the two-dimensional buoyancy-driven fluid flows developed in both the water-filled subenclosure and the molten PCM region of the PCM-filled subenclosure were modeled as laminar Newtonian fluid flow adhering to the Boussinesq approximation. Meanwhile, two-dimensional conduction is accounted for the plate heater as well as the solid PCM zone. Numerical results are presented to unveil the cooling behavior of the pulsating heat-generating plate through the PCM substrate and water-filled enclosure. Results of the parametric simulations reveal that the water layer has the better capability of heat dissipation than the PCM substrate. Heat dissipation from the plate through the PCM substrate is mainly via the latent heat absorption as associated with melting phenomenon. Moreover, numerical results obtained are compared with the corresponding experiments.     

A note on a new theory of turbulence

In this paper,the creation and annihilaton of turbulent eddies are described aselementary particles in the quantum field theory.An elementary particle may be consideredas a solid entity as it exists in quantum theory,but a turbulent eddy is often changed in sizeand shape with time due to its energy dissipation in a turbulent field.Therefore,in order toapply the method of the quantum field theory to the turbulent field by analogy,the entity ofthe same eddy should be defined firstly.According to the linearized theory,the turbulenteddies with the similarity character in lime duration may be considered as the entity of thesame eddy,and the creation and annihilation of turbulent eddies without the similarcharacters are related to the interaction term φ_2 in equation(2.6).Then,the creationoperator and annihilation operator similar to those in the quantum field theory are used todescribe the state of turbulent eddy field.Finally,a“Schrǒdinger”equation of turbulenteddies is formulated based upon the nonline