scholarly journals From a steady plume to periodic puffs during confined carbon dioxide dissolution

2018 ◽  
Vol 855 ◽  
pp. 1-27
Author(s):  
P. Meunier ◽  
F. Nadal

In this paper, the stability of a laminar plume due to solutal convection is addressed from experimental, numerical and theoretical points of view. A topless vertical tube containing water is put in a pressure cell filled with carbon dioxide ( $\text{CO}_{2}$ ). The diffusion of $\text{CO}_{2}$ at the free surface creates a thin layer of heavy fluid underneath the surface. This unstable density gradient generates a steady laminar plume which goes downward through the entire tube. A quasi-steady flow settles in the tube, filling gradually the bottom of the tube with heavy fluid. During this laminar regime, the velocity of the plume slowly decreases due to the build-up of the background density gradient. Surprisingly, despite the decrease of the Reynolds number, the laminar plume suddenly destabilises via a varicose mode into periodic pulsed puffs after an onset time which depends on the height of the tube and on the solutal Rayleigh number $Ra$ . This periodic regime is followed by an aperiodic regime, which lasts until the complete saturation of the solution. The observed destabilisation is explained as a result of the interplay between the feedback of the global recirculating flow and the progressive density stratification of the background fluid. The wavelength, frequency, onset time and phase velocity of the instability are explored using particle image velocimetry (PIV) measurements over two decades of Rayleigh number. The characteristics of the instability appear to be almost independent of the Bond number but strongly dependent on the solutal Rayleigh number and the aspect ratio. The phase velocity is very close to the fluid velocity of the plume before the instability, which has been predicted in various works to scale as $Ra^{2/3}(\ln \,Ra)^{1/3}$ . The wavelength is close to 4.5 times the radius of the cylinder (independent of aspect ratio, Bond number and Rayleigh number) such that the frequency scales as the phase velocity. The onset time, which is proportional to the height of the cylinder, scales as $Ra^{-0.55}$ and depends on the Bond number. A simplified model inspired from Lorenz’ waterwheel is proposed to explain the destabilisation process after partial fill-up of the cylinder. Although very qualitative, the model captures the key features of the experimental observations.

Author(s):  
Salaika Parvin ◽  
Nepal Chandra Roy ◽  
Litan Kumar Saha ◽  
Sadia Siddiqa

A numerical study is performed to investigate nanofluids' flow field and heat transfer characteristics between the domain bounded by a square and a wavy cylinder. The left and right walls of the cavity are at constant low temperature while its other adjacent walls are insulated. The convective phenomena take place due to the higher temperature of the inner corrugated surface. Super elliptic functions are used to transform the governing equations of the classical rectangular enclosure into a system of equations valid for concentric cylinders. The resulting equations are solved iteratively with the implicit finite difference method. Parametric results are presented in terms of streamlines, isotherms, local and average Nusselt numbers for a wide range of scaled parameters such as nanoparticles concentration, Rayleigh number, and aspect ratio. Several correlations have been deduced at the inner and outer surface of the cylinders for the average Nusselt number, which gives a good agreement when compared against the numerical results. The strength of the streamlines increases significantly due to an increase in the aspect ratio of the inner cylinder and the Rayleigh number. As the concentration of nanoparticles increases, the average Nusselt number at the internal and external cylinders becomes stronger. In addition, the average Nusselt number for the entire Rayleigh number range gets enhanced when plotted against the volume fraction of the nanofluid.


2011 ◽  
Vol 71-78 ◽  
pp. 1187-1190
Author(s):  
Yan Lai Zhang ◽  
Zhong Hao Rao ◽  
Shuang Feng Wang ◽  
Hong Zhang ◽  
Li Jun Li ◽  
...  

This experiment is performed to investigate heat transfer characteristics with the PCM microcapsule slurry in a solid phase state at a horizontal rectangular enclosure heating from below and cooling from top. Some important parameters are taken into account such as the mass concentration of the PCM, the temperature difference between heating plate and cooling plate, Nusselt number Nu, Rayleigh number Ra and the aspect ratio (width/height) of the horizontal rectangular enclosure. Experiment is done under the thermal steady condition in the PCM microcapsule slurry. Heat transfer coefficient is measured under various temperature differences in PCM mass concentrations of 10% and 20%. And relationship with Nusselt number Nu and Rayleigh number Ra is summarized to various heights H or the aspect ratio (width/height) Ar of enclosure.


1978 ◽  
Vol 100 (2) ◽  
pp. 205-211 ◽  
Author(s):  
L. A. Clomburg

Laminar natural convection in a two-dimensional enclosure with both source (uniform heat flux density) and sink (temperature specified) located on the top horizontal boundary is investigated numerically. Temperature and velocity profiles are presented for a high Prandtl number fluid for length Rayleigh numbers in the range 107 to 109 for length to depth ratios of 1:1 to 4:1. To generalize the results, an order of magnitude analysis is used to determine the dependence of temperature, velocity, and boundary-layer thickness scales on aspect ratio and Rayleigh number. The numerical data are well correlated using these suggested scales. The analysis shows the Nusselt number and the maximum horizontal velocity to depend on the 1/6 and 1/3 powers of the Rayleigh number, independent of aspect ratio.


2019 ◽  
Vol 20 (1) ◽  
pp. 229-244
Author(s):  
Mehdi Ahmadi ◽  
Seyed Ali Agha Mirjalily ◽  
Seyed Amir Abbas Oloomi

ABSTRACT: This study is conducted to investigate turbulent natural convection flow in an enclosure with thermal sources using the low-Reynolds number (LRN) k-? model. This enclosure has a cold source with temperature Tc and a hot source with temperature Th as thermal sources, other walls of the enclosure are adiabatic. The aim of this study is to predict the effect of change in Rayleigh number, repositioning of cold and hot sources, and thermal sources aspect ratio on the flow field, temperature, and rate of heat transfer. To achieve this aim, the equations of continuity, momentum, energy, turbulent kinetic energy, and kinetic energy dissipation are employed in the case of 2D turbulence with constant thermo-physical properties except the density in the buoyancy term (Boussinesq approximation). To numerically solve these equations, the finite volume method and SIMPLE algorithm are used. According to the modeling results, the most optimal temperature distribution in the enclosure is seen when the hot source is below the cold source. With decreasing distance between hot and cold sources, heat transfer rate increases. The maximal heat transfer rate is derived via study of the heating sources aspect ratio. In constant positions of cold and hot sources on a wall, the heat transfer rate increases with increasing Rayleigh number (Ra=109-1011). ABSTAK: Kajian ini dijalankan bagi mengkaji perubahan semula jadi aliran perolakan dalam tempat tertutup dengan sumber haba menggunakan model k-? nombor Reynolds-rendah (LRN). Bekas tertutup ini mempunyai dua sumber haba iaitu sumber sejuk dengan suhu Tc dan sumber panas dengan suhu Th, manakala dinding lain bekas ini adalah adiabatik. Tujuan kajian ini adalah bagi mengesan perubahan nombor Rayleigh, mengubah sumber sejuk dan panas dan nisbah sumber haba kepada kawasan aliran, suhu dan halaju perubahan haba. Bagi mencapai tujuan tersebut, persamaan sambungan, momentum, tenaga, tenaga kinetik perolakan, dan pengurangan tenaga kinetik telah dilaksanakan dalam kes perolakan 2D dengan sifat fizikal-haba berterusan (malar) kecuali isipadu terma keapungan (anggaran Boussinesq). Bagi menyelesaikan persamaan ini secara berangka, kaedah isipadu terhad dan algorithma MUDAH telah digunakan. Berdasarkan keputusan model, suhu distribusi optimal dalam bekas tertutup dilihat apabila sumber panas adalah kurang daripada sumber sejuk. Dengan pengurangan jarak antara sumber panas dan sejuk, kadar pertukaran haba meningkat. Kadar pertukaran haba maksima telah diperoleh melalui kajian nisbah  aspek sumber pemanasan. Kadar pertukaran haba bertambah dengan bertambahnya nombor Rayleigh  (Ra=109-1011), pada posisi tetap sumber sejuk dan panas pada dinding bekas.


Author(s):  
C. Israel-Cookey ◽  
L. Ebiwareme ◽  
E. Amos

In this research article, the effect of Vadasz number on magnetoconvection in a Darcy Porous Layer with concentration based internal heating is studied. The flow is governed by the Oberbeck-Boussineq model for Newtonian fluid. The stability analysis method based on the perturbation of infinitesimal amplitude is carried out using the normal mode analysis. The onset criterion for both the stationary and oscillatory convection on the stability of system is obtained. The analysis examines the effects of pertinent parameters on the stability of the system: magnetic field parameter, solutal Rayleigh number, Lewis number and Vadasz number. The result show that, internal heat parameter,  and Lewis number, , hastens the onset of instability in the system, whereas magnetic field, , Vadasz number,  and solutal Rayleigh number,  delay the onset of instability.


Water ◽  
2020 ◽  
Vol 12 (3) ◽  
pp. 738
Author(s):  
Holger Class ◽  
Kilian Weishaupt ◽  
Oliver Trötschler

Carbon dioxide density-driven dissolution in a water-filled laboratory flume of the dimensions 60 cm length, 40 cm height, 1 cm thickness, was visualized using a pH-sensitive color indicator. We focus on atmospheric pressure conditions, like in caves where CO2 concentrations are typically higher. Varying concentrations of carbon dioxide were applied as boundary conditions at the top of the experimental setup, leading to the onset of convective fingering at differing times. The data were used to validate a numerical model implemented in the numerical simulator DuMux. The model solves the Navier–Stokes equations for density-induced water flow with concentration-dependent fluid density and a transport equation, including advective and diffusive processes for the carbon dioxide dissolved in water. The model was run in 2D, 3D, and pseudo-3D on two different grids. Without any calibration or fitting of parameters, the results of the comparison between experiment and simulation show satisfactory agreement with respect to the onset time of convective fingering, and the number and the dynamics of the fingers. Grid refinement matters, in particular, in the uppermost part where fingers develop. The 2D simulations consistently overestimated the fingering dynamics. This successful validation of the model is the prerequisite for employing it in situations with background flow and for a future study of karstification mechanisms related to CO2-induced fingering in caves.


2012 ◽  
Vol 713 ◽  
pp. 216-242 ◽  
Author(s):  
Jun Hu ◽  
Daniel Henry ◽  
Xie-Yuan Yin ◽  
Hamda BenHadid

AbstractThree-dimensional Rayleigh–Bénard instabilities in binary fluids with Soret effect are studied by linear biglobal stability analysis. The fluid is confined transversally in a duct and a longitudinal throughflow may exist or not. A negative separation factor $\psi = \ensuremath{-} 0. 01$, giving rise to oscillatory transitions, has been considered. The numerical dispersion relation associated with this stability problem is obtained with a two-dimensional Chebyshev collocation method. Symmetry considerations are used in the analysis of the results, which allow the classification of the perturbation modes as ${S}_{l} $ modes (those which keep the left–right symmetry) or ${R}_{x} $ modes (those which keep the symmetry of rotation of $\lrm{\pi} $ about the longitudinal mid-axis). Without throughflow, four dominant pairs of travelling transverse modes with finite wavenumbers $k$ have been found. Each pair corresponds to two symmetry degenerate left and right travelling modes which have the same critical Rayleigh number ${\mathit{Ra}}_{c} $. With the increase of the duct aspect ratio $A$, the critical Rayleigh numbers for these four pairs of modes decrease and closely approach the critical value ${\mathit{Ra}}_{c} = 1743. 894$ obtained in a two-dimensional situation, one of the mode (a ${S}_{l} $ mode called mode A) always remaining the dominant mode. Oscillatory longitudinal instabilities ($k\approx 0$) corresponding to either ${S}_{l} $ or ${R}_{x} $ modes have also been found. Their critical curves, globally decreasing, present oscillatory variations when the duct aspect ratio $A$ is increased, associated with an increasing number of longitudinal rolls. When a throughflow is applied, the symmetry degeneracy of the pairs of travelling transverse modes is broken, giving distinct upstream and downstream modes. For small and moderate aspect ratios $A$, the overall critical Rayleigh number in the small Reynolds number range studied is only determined by the upstream transverse mode A. In contrast, for larger aspect ratios as $A= 7$, different modes are successively dominant as the Reynolds number is increased, involving both upstream and downstream transverse modes A and even the longitudinal mode.


1991 ◽  
Vol 113 (1) ◽  
pp. 97-107 ◽  
Author(s):  
S. Ramanathan ◽  
R. Kumar

This paper presents the numerical results of natural convective flows between two vertical, parallel plates within a large enclosure. A parametric study has been conducted for various Prandtl numbers and channel aspect ratios. The results are in good agreement with the reported results in the literature for air for large aspect ratios. However, for small aspect ratios, the present numerical results do not agree with the correlations given in the literature. The discrepancy is due to the fact that the published results were obtained for channels where the diffusion of thermal energy in the vertical direction is negligible. The results obtained in this paper indicate that vertical conduction should be considered for channel aspect ratios less than 10 for Pr = 0.7. Correlations are presented to predict the maximum temperature and the average Nusselt number on the plate as explicit functions of the channel Rayleigh number and the channel aspect ratio for air. The plate temperature is a weak function of Prandtl number for Prandtl numbers greater than 0.7, if the channel Rayleigh number is chosen as the correlating parameter. For Prandtl numbers less than 0.1, the plate temperature is a function of the channel Rayleigh number and the Prandtl number. A correlation for maximum temperature on the plate is presented to include the Prandtl number effect for large aspect ratio channels.


2020 ◽  
Vol 143 (2) ◽  
Author(s):  
Basanta Kumar Rana ◽  
Bhajneet Singh ◽  
Jnana Ranjan Senapati

Abstract Numerical investigations are performed on natural and mixed convection around stationary and rotating vertical heated hollow cylinder with negligible wall thickness suspended in the air. The fluid flow and heat transfer characterization around the hollow cylinder are obtained by varying the following parameters, namely, Rayleigh number (Ra), Reynolds number (ReD), and cylindrical aspect ratio (L/D). The heat transfer quantities are estimated by varying the Rayleigh number (Ra) from 104 to 108 and aspect ratio (L/D) ranging from 1 to 20. Steady mixed convection with active rotation of hollow vertical cylinder is further studied by varying the Reynolds number (ReD) from 0 to 2100. The velocity vectors and temperature contours are shown in order to understand the fluid flow and heat transfer around the vertical hollow cylinder for both rotating and nonrotating cases. The surface average Nusselt number trends are presented for various instances of Ra, ReD, and L/D and found out that the higher rate of heat loss from the cylinder wall occurs at high Ra, low L/D (short cylinder) and high ReD.


Sign in / Sign up

Export Citation Format

Share Document