High Rayleigh number variational multiscale large eddy simulations of Rayleigh-Bénard convection
Not Available
Sondak, David↗
SEARCH · Search NASA
Results for “RAYLEIGH NUMBER”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Reduced modeling of porous media convection in a minimal flow unit at large Rayleigh number
Not Available
Wen, Baole (ORCID:0000000320731508)↗
Inclined porous medium convection at large Rayleigh number
High-Rayleigh-number ( $Ra$ ) convection in an inclined two-dimensional porous layer is investigated using direct numerical simulations (DNS) and stability and variational upper-bound analyses. When the inclination angle $\unicode[STIX]{x1D719}$ of the layer satisfies $0^{\circ }<\unicode[STIX]{x1D719}\lesssim 25^{\circ }$ , DNS confirm that the flow exhibits a three-region wall-normal asymptotic structure in accord with the strictly horizontal ( $\unicode[STIX]{x1D719}=0^{\circ }$ ) case, except that as $\unicode[STIX]{x1D719}$ is increased the time-mean spacing between neighbouring interior plumes also increases substantially. Both DNS and upper-bound analysis indicate that the heat transport enhancement factor (i.e. the Nusselt number) $Nu\sim CRa$ with a $\unicode[STIX]{x1D719}$ -dependent prefactor $C$ . When $\unicode[STIX]{x1D719}>\unicode[STIX]{x1D719}_{t}$ , however, where $30^{\circ }<\unicode[STIX]{x1D719}_{t}<32^{\circ }$ independently of $Ra$ , the columnar flow structure is completely broken down: the flow transitions to a large-scale travelling-wave convective roll state, and the heat transport is significantly reduced. To better understand the physics of inclined porous medium convection at large $Ra$ and modest inclination angles, a spatial Floquet analysis is performed, yielding predictions of the linear stability of numerically computed, fully nonlinear steady convective states. The results show that there exist two types of instability when $\unicode[STIX]{x1D719}\neq 0^{\circ }$ : a bulk-mode instability and a wall-mode instability, consistent with previous findings for $\unicode[STIX]{x1D719}=0^{\circ }$ (Wen et al. , J. Fluid Mech. , vol. 772, 2015, pp. 197–224). The background flow induced by the inclination of the layer intensifies the bulk-mode instability during its subsequent nonlinear evolution, thereby favouring increased spacing between the interior plumes relative to that observed in convection in a horizontal porous layer.
Mechanics↗
THE SPECTRAL AMPLITUDE OF STELLAR CONVECTION AND ITS SCALING IN THE HIGH-RAYLEIGH-NUMBER REGIME
Convection plays a central role in the dynamics of any stellar interior, and yet its operation remains largely hidden from direct observation. As a result, much of our understanding concerning stellar convection necessarily derives from theoretical and computational models. The Sun is, however, exceptional in that regard. The wealth of observational data afforded by its proximity provides a unique test bed for comparing convection models against observations. When such comparisons are carried out, surprising inconsistencies between those models and observations become apparent. Both photospheric and helioseismic measurements suggest that convection simulations may overestimate convective flow speeds on large spatial scales. Moreover, many solar convection simulations have difficulty reproducing the observed solar differential rotation owing to this apparent overestimation. We present a series of three-dimensional stellar convection simulations designed to examine how the amplitude and spectral distribution of convective flows are established within a star’s interior. While these simulations are nonmagnetic and nonrotating in nature, they demonstrate two robust phenomena. When run with sufficiently high Rayleigh number, the integrated kinetic energy of the convection becomes effectively independent of thermal diffusion, but the spectral distribution of that kinetic energy remains sensitive to both of these quantities. A simulation that has converged to a diffusion-independent value of kinetic energy will divide that energy between spatial scales such that low-wavenumber power is overestimated and high-wavenumber power is underestimated relative to a comparable system possessing higher Rayleigh number. We discuss the implications of these results in light of the current inconsistencies between models and observations.
79 ASTRONOMY AND ASTROPHYSICS↗
Parametric model-order-reduction development for unsteady convection
A time-averaged error indicator with POD- h Greedy is developed to drive parametric model order reduction (pMOR) for 2D unsteady natural convection in a high-aspect ratio slot parameterized with the Prandtl number, Rayleigh number, and slot angle with respect to the gravity. The error indicator is extended to accommodate the energy equation and Leray regularization. Despite being two-dimensional and laminar, the target flow regime presents several challenges: 1) there is a bifurcation in the angle parameter space; 2) the solution can be multivalued, even at steady state; and 3) the solution exhibits spatio-temporal chaos at several points in the parameter space. The authors explore several reduced-order models (ROMs) and demonstrate that Leray-regularized Galerkin ROMs provide a robust solution approach for this class of flows. They further demonstrate that error-indicated pMOR can efficiently predict several QOIs, such as mean flow, mean Nusselt number and mean turbulent kinetic energy, even in the presence of a bifurcation. Finally, they show that spatio-temporal chaos can lead to lack of reproducibility in both the full-order model and the reduced-order model and that the variance in the full-order model provides a lower bound on the pMOR error in these cases.
leray regularization↗
Experimental Investigation of Natural Convection Heat Transfer in Heated Vertical Tubes with Discrete Rings
Natural convection heat transfer from heated vertical pipes has been investigated experimentally. The pipes are open-ended and circular in cross section. The test sections are electrically heated imposing the circumferentially and axially constant wall heat flux. The outer surface of the test section is insulated and heat is allowed to dissipate from the internal surface. Experiments are carried out for four different channels of 45 mm internal diameter and 3.8 mm thickness with length varying from 450 mm to 850 mm. Ratios of length to diameter of the channel are taken as L/D = 10, 12.22, 15.56, and 18.89. Wall heat fluxes are maintained at q{sup //} = 250 to 3341 W/m{sup 2}. Experiments are also carried out on intensified tubes of the same geometrical sizes with the discrete rings of rectangular cross section provided on the internal surface. Thickness of the rings are taken as t = 4, 6 and 8 mm; step size (s) of the rings are varying from 75 mm to 283.3 mm and the other ratios are taken as s/D = 1.67 to 6.29, t/D = 0.088 to 0.178 and s / t = 9.38 to 70.82. The effects of L/D ratio, ring thickness, ring spacing and wall heating conditions on local steady-state heat transfer phenomena are studied. A correlation is developed for average Nusselt number and modified Rayleigh number and another correlation is also developed for modified Rayleigh number and modified Reynolds number. These correlations predict the data accurately within ±10% error.
42 ENGINEERING↗
Improved natural convection heat transfer correlations for reactor cavity cooling systems of high-temperature gas-cooled reactors: From computational fluid dynamics to Pronghorn
The Reactor Cavity Cooling System (RCCS) is a common reactor safety system in High Temperature Gas Cooled Reactors (HTGR) that removes heat from the Reactor Pressure Vessel (RPV) by radiation ($\sim 80\%$) and natural convection ($\sim 20\%$). For simulation of accident scenarios of HTGRs, intermediate fidelity and system codes models must be employed for limiting the models' execution time. While accurate quantification of the radiative heat transfer is available in these models, quantification of natural convection must rely on correlations of questionable accuracy for the Nusselt number. Commonly used correlations are based in experiments performed at low Rayleigh numbers and/or using isothermal walls in simplified geometries. Here, this work improves on the accuracy of natural convection heat transfer correlations in support for HTGR designs. These correlations include both local and average Nusselt numbers as a function of the global Rayleigh number, the local Rayleigh number, and the temperature profile at the hot wall of the RCCS. In the absence of dedicated experiments and the difficulty of performing high-fidelity simulations at realistic Rayleigh numbers, the data to fit the correlations are generated with Computational Fluid Dynamics (CFD) using Reynolds Averaged Navier-Stokes (RANS) models. First, a careful selection of the RANS turbulence model is performed by comparing the results obtained with different RANS turbulence models against high fidelity simulations of natural convection at $Ra \ 1 \times 10^{11}$ in a rectangular cavity. Next, the selected model is used to perform simulations of an HTGR cavity at different high Rayleigh numbers $\in [6.1 \times 10^{11},2.9 \times 10^{13}]$ to encompass several HTGR designs, assuming an isothermal RPV wall. The results obtained are used to fit a correlation for the average and space-varying Nusselt number as a function of the global and local Rayleigh numbers via a sparsity-promoting least-squares method. The selected RANS model is then used to perform simulations of a PBMR 400 HTGR cavity with the temperature profiles at the RPV wall obtained during a PLOFC transient. We use the results obtained to fit a temperature-dependent correction to the space-varying Nusselt number with the sparsity-promoting least-squares method. The results obtained in this work, enable system-level codes, such as Pronghorn, to perform higher-fidelity simulations of the heat exchange process in the RCCS while still maintaining a low computational cost.
42 ENGINEERING↗
Fluxes across double-diffusive interfaces : a one-dimensional turbulence study.
The parameterization of the fluxes of heat and salt across double-diffusive interfaces is of interest in geophysics, astrophysics, and engineering. The present work is a parametric study of these fluxes using one-dimensional-turbulence (ODT) simulations. Its main distinction is that it considers a parameter space larger than previous studies. Specifically, this work considers the effect on the fluxes of the stability parameter R{sub {rho}}, Rayleigh number Ra, Prandtl number, Lewis number, and Richardson number. The ratio Ra/R{sub {rho}} is found to be a dominant parameter. Here Ra/R{sub {rho}} can be seen as a ratio of destabilizing and stabilizing effects. Trends predicted by the simulations are in good agreement with previous models and available measurements.
58 GEOSCIENCES↗
Effect of adiabatic wall on the natural convection heat transfer from a wavy surface created by attached horizontal cylinders
The effect of a vertical adiabatic wall on the natural convection heat transfer from vertical array of attached cylinders, which can be considered as wavy surface, was investigated experimentally and numerically. The experiments were carried out using Mach-Zehnder interferometer and the commercial FLUENT code was used for numerical study. This paper focuses on the effect of wall-wavy surface spacing and Rayleigh number variation on the local and average free convection heat transfer coefficients from the each cylinder and the wavy surface. Rayleigh number ranges from 2400 to 10,000 and from 300,000 to 1,250,000 based on cylinder diameter and wavy surface height respectively. The local and average Nusselt numbers were determined for the different Rayleigh numbers, and the ratio of wall- wavy surface spacing to cylinder diameter 0.75, 1, 1.5, 2, 3, 4, 5, and {infinity}. Results are indicated with a single correlation which gives the average Nusselt number as a function of the ratio of the wall-wavy surface spacing to cylinder diameter and the Rayleigh numbers. There is an optimum distance between the wall and wavy surface in which the Nusselt number attain its maximum value. This optimum distance depends on the Rayleigh number. (author)
14 SOLAR ENERGY↗
Turbulent Natural Convection in a Square Cavity with a Circular Cylinder
In this paper, numerical simulations of high Rayleigh number flows (10 8 -10 10 ) were conducted to investigate the turbulent fluid flow and thermal characteristics of natural convection induced by a centrally placed hot cylinder in a cold square enclosure. The effect of the aspect ratio (radius of the cylinder to the side of the cavity) was investigated for three values (0.1, 0.2, and 0.3) for each Rayleigh number. Effects of turbulence induced by the high Rayleigh number (>10 7 ) were computed by using the unsteady k-ω model. A spectral-element method with high polynomial order (high resolution) was used to solve the system of unsteady time-averaged equations of continuity, momentum, and energy, along with the turbulence equations. Detailed comparison with other numerical work is presented. Contours of velocity, temperature, and turbulence quantities are presented for various high Rayleigh numbers. Also presented is the influence of the Rayleigh number on the local Nusselt number on the centrally placed hot cylinder and the cold enclosure walls. Time-marching results show that the steady-state solutions can be obtained even for high Rayleigh numbers considered in this study. The results also show that the average and peak Nusselt numbers roughly double for each order of magnitude increase of the Rayleigh number for all radii considered. Finally, a correlation for the average Nusselt number as a function of Rayleigh number and aspect ratio is also presented.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Scaling and excitation of combined convection in a rapidly rotating plane layer
The optimum (to my mind) scaling of the combined thermal and compositional convection in a rapidly rotating plane layer is proposed.This scaling follows from self-consistent estimates of typical physical quantities. Similarity coefficients are introduced for the ratio convection dissipation/convection generation (s) and the ratio thermal convection/compositional convection (r). The third new and most important coefficient δ is the ratio of the characteristic size normal to the axis of rotation to the layer thickness. The faster the rotation, the lower δ. In the case of the liquid Earth core, δ ~ 10{sup –3} substitutes for the generally accepted Ekman number (E ~ 10{sup –15}) and s ~ 10{sup –6} substitutes for the inverse Rayleigh number 1/Ra ~ 10{sup –30}. It is found that, at turbulent transport coefficients, number s and the Prandtl number are on the order of unity for any objects and δ is independent of transport coefficients. As a result of expansion in powers of δ, an initially 3D system of six variables is simplified to an almost 2D system of four variables without δ. The problem of convection excitation in the main volume is algebraically solved and this problem for critical values is analytically solved. Dispersion relations and general expressions for critical wavenumbers, numbers s (which determine Rayleigh numbers), other critical parameters, and asymptotic solutions are derived. Numerical estimates are made for the liquid cores in the planets that resemble the Earth. Further possible applications of the results obtained are proposed for the interior of planets, moons, their oceans, stars, and experimental objects.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
RANS simulation of natural convection in an enclosure with curved surface heating from below
Prediction accuracy of natural convection is strongly dependent on the choice of turbulence model. In this paper, different turbulence models as well as turbulent heat flux models are employed to predict the natural convection in an enclosure with curved surface heating from below and investigate their performance. Based on the simulation results, it is found that k-ε model which is well known to overestimate the production of kinetic energy provides a higher prediction of Nusselt number while k-ε-ζ-f model provides a relatively smaller estimation. A systematic CFD simulation is carried out for Rayleigh numbers in the range 10{sup 9} to 10{sup 11}, with which the influence from aspect ratio is also investigated. The effect of Rayleigh number on heat transfer is found to be significant and the dependency of Nusselt number (Nu) is obtained as a function of Rayleigh number (Ra). The influence from the aspect ratio of the simulation domain is found to be dependent on the choice of turbulence model. For k-ε model, Nu increases as aspect ratio increases, while an opposite dependency is observed when k-ε-ζ-f model is used. For this discrepancy, arguments are provided based on the streamlines and the profiles of Nu over the heating surface. (authors)
42 ENGINEERING↗
Equatorial symmetry of Boussinesq convective solutions in a rotating spherical shell allowing rotation of the inner and outer spheres
We investigate properties of convective solutions of the Boussinesq thermal convection in a moderately rotating spherical shell allowing the respective rotation of the inner and outer spheres due to the viscous torque of the fluid. The ratio of the inner and outer radii of the spheres, the Prandtl number, and the Taylor number are fixed to 0.4, 1, and 500{sup 2}, respectively. The Rayleigh number is varied from 2.6 × 10{sup 4} to 3.4 × 10{sup 4}. In this parameter range, the behaviours of obtained asymptotic convective solutions are almost similar to those in the system whose inner and outer spheres are restricted to rotate with the same constant angular velocity, although the difference is found in the transition process to chaotic solutions. The convective solution changes from an equatorially symmetric quasi-periodic one to an equatorially symmetric chaotic one, and further to an equatorially asymmetric chaotic one, as the Rayleigh number is increased. This is in contrast to the transition in the system whose inner and outer spheres are assumed to rotate with the same constant angular velocity, where the convective solution changes from an equatorially symmetric quasi-periodic one, to an equatorially asymmetric quasi-periodic one, and to equatorially asymmetric chaotic one. The inner sphere rotates in the retrograde direction on average in the parameter range; however, it sometimes undergoes the prograde rotation when the convective solution becomes chaotic.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Natural convection in a cubical cavity with a coaxial heated cylinder
High-resolution three-dimensional simulations were conducted to investigate the velocity and temperature fields in a cold cubical cavity due to natural convection induced by a centrally placed hot cylinder. Unsteady, incompressible Navier-Stokes equations were solved by using a spectral- element method for Rayleigh numbers ranging from 10 3 to 10 9 . The effect of spanwise thermal boundary conditions, aspect ratio (radius of the cylinder to the side of the cavity), and spanwise temperature distribution of the inner cylinder on the velocity and thermal fields were investigated for each Rayleigh number. Results from two-dimensional calculations were compared with three-dimensional simulations. The 3D results indicate a complex flow structure in the vicinity of the spanwise walls. The results also show that the imposed thermal wall boundary condition impacts the flow and temperature fields strongly near the spanwise walls. The variation of the local Nusselt number on the cylinder surface and enclosure walls at various spanwise locations was also investigated. The local Nusselt number on the cylinder surface and enclosure walls at the cavity mid-plane (Z = 0) is close to 2D simulations for 10 3 ≤ Ra ≤ 10 8 . Simulations also show a variation in the local Nusselt number, on both the cylinder surface and the enclosure walls, in the spanwise direction, for all Rayleigh numbers studied in this work. The results also indicate that if the enclosure walls are insulated in the spanwise direction (as opposed to a constant temperature), the peak Nusselt number on the enclosure surface occurs near the spanwise walls and is about 20% higher than the peak Nusselt number at the cavity mid-plane. The temporal characteristics of 3D flows are also different from 2D results for Ra > 10 8 . Furthermore, these results suggest that 3D simulations would be more appropriate for flows with Ra > 10 8 .
42 ENGINEERING↗
EFFECTS OF STELLAR FLUX ON TIDALLY LOCKED TERRESTRIAL PLANETS: DEGREE-1 MANTLE CONVECTION AND LOCAL MAGMA PONDS
We model the geodynamical evolution of super-Earth exoplanets in synchronous rotation about their star. While neglecting the effects of a potential atmosphere, we explore the parameter spaces of both the Rayleigh number and intensity of incoming stellar flux, and identify two main stages of mantle convection evolution. The first is a transient stage in which a lithospheric temperature and thickness dichotomy emerges between the substellar and the antistellar hemispheres, while the style of mantle convection is dictated by the Rayleigh number. The second stage is the development of degree-1 mantle convection. Depending on mantle properties, the timescale of onset of this second stage of mantle evolution varies from order 1 to 100 billion years of simulated planetary evolution. Planets with higher Rayleigh numbers (due to, for instance, larger planetary radii than the Earth) and planets whose incoming stellar flux is high (likely for most detectable exoplanets) will develop degree-1 mantle convection most quickly, on the order of 1 billion years, which is within the age of many planetary systems. Surface temperatures range from 220 K to 830 K, implying the possibility of liquid water in some regions near the surface. These results are discussed in the context of stable molten magma ponds on hotter planets, and the habitability of super-Earths which may lie outside the Habitable Zone.
79 ASTRONOMY AND ASTROPHYSICS↗
Convective carbon dioxide dissolution in a closed porous medium at low pressure
Motivated by the persistence of natural carbon dioxide ( $\text{CO}_{2}$ ) fields, we investigate the convective dissolution of $\text{CO}_{2}$ at low pressure (below 1 MPa) in a closed system, where the pressure in the gas declines as convection proceeds. This introduces a negative feedback that reduces the convective dissolution rate even before the brine becomes saturated. We analyse the case of an ideal gas with a solubility given by Henry’s law, in the limits of very low and very high Rayleigh numbers. The equilibrium state in this system is determined by the dimensionless dissolution capacity, $\unicode[STIX]{x1D6F1}$ , which gives the fraction of the gas that can be dissolved into the underlying brine. Analytic approximations of the pure diffusion problem with $\unicode[STIX]{x1D6F1}>0$ show that the diffusive base state is no longer self-similar and that diffusive mass transfer declines rapidly with time. Direct numerical simulations at high Rayleigh numbers show that no constant flux regime exists for $\unicode[STIX]{x1D6F1}>0$ ; nevertheless, the quantity $F/C_{s}^{2}$ remains constant, where $F$ is the dissolution flux and $C_{s}$ is the dissolved concentration at the top of the domain. Simple mathematical models are developed to predict the evolution of $C_{s}$ and $F$ for high-Rayleigh-number convection in a closed system. The negative feedback that limits convection in closed systems may explain the persistence of natural $\text{CO}_{2}$ accumulations over millennial time scales.
Mechanics↗
Opposing flow in square porous annulus: Influence of Dufour effect
Heat and mass transfer in porous medium is very important area of research which is also termed as double diffusive convection or thermo-solutal convection. The buoyancy ratio which is the ratio of thermal to concentration buoyancy can have negative values thus leading to opposing flow. This article is aimed to study the influence of Dufour effect on the opposing flow in a square porous annulus. The partial differential equations that govern the heat and mass transfer behavior inside porous medium are solved using finite element method. A three node triangular element is used to divide the porous domain into smaller elements. Results are presented with respect to geometric and physical parameters such as duct diameter ratio, Rayleigh number, radiation parameter etc. It is found that the heat transfer increase with increase in Rayleigh number and radiation parameter. It is observed that Dufour coefficient has more influence on velocity profile.
42 ENGINEERING↗