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↗
High Rayleigh number heat transfer in a horizontal cylinder with adiabatic wall
The present investigation is concerned with an experimentally guided approach to the estimation of Nusselt numbers (Nu) at high Rayleigh numbers (Ra) for a cylinder with an adiabatic side wall. The Rayleigh number dependence of the Nusselt number for a horizontal cylinder with an adiabatic wall is presented in a graph. The obtained data are compared with results reported by Shih (1981). Shih has extended a three-term expansion for velocity and temperature distributions reported by Bejan and Tien (1978).
Schiroky, G. H.↗
Stellar convection. III - Convection at large Rayleigh numbers
The results of a numerical, three-dimensional study of convection in a self-gravitating sphere of Boussinesq fluid with a Rayleigh number of 10-billion and a Prandtl number of 1 are presented. The velocity and temperature are computed by using spectral methods in the horizontal and finite-differencing in the radial directions. An eddy viscosity and diffusivity are needed to model the subresolution flow. For Rayleigh numbers much less than 10-billion the flows do not have well-defined inertial ranges, and an eddy viscosity and diffusivity cannot be assigned in this self-consistent manner. By computing the energy spectra as well as the detailed energy budgets as a function of wavenumber, it is shown that for Rs = 10-billion there is an inertial range for the modes corresponding to spherical harmonics with l greater than 6.
Marcus, P. S.↗
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↗
Localization of toroidal motion and shear heating in 3-D high Rayleigh number convection with temperature-dependent viscosity
We have applied spectral-transform methods to study three-dimensional thermal convection with temperature-dependent viscosity. The viscosity varies exponentially with the form exp(-BT), where B controls the viscosity contrast and T is temperature. Solutions for high Rayleigh numbers, up to an effective Ra of 6.25 x 10(exp 6), have been obtained for an aspect-ratio of 5x5x1 and a viscosity contrast of 25. Solutions show the localization of toroidal velocity fields with increasing vigor of convection to a coherent network of shear-zones. Viscous dissipation increases with Rayleigh number and is particularly strong in regions of convergent flows and shear deformation. A time-varying depth-dependent mean-flow is generated because of the correlation between laterally varying viscosity and velocity gradients.
Balachandar, S.↗
Stellar convection 3: Convection at large Rayleigh numbers
A three dimensional study of convection in a self gravitating sphere of Boussinesq fluid with a Rayleigh number of 10 to the 10th power and a Prandtl of 1 is presented. The velocity and temperature of the fluid are computed at the largest wavelengths using spectral methods. A confirmation that the fluid is anisotropic and that the energy spectra are not smooth functions of wavelength but have a large amount of fine structure is discussed. The parameterization of the transport properties of the unresolvable inertial subrange with eddy viscosities and diffusivities is described. The time dependent fluctuations in the energy spectra and how they cascade from large to small wavelengths is examined.
Marcus, P. S.↗
High Rayleigh number convection in rectangular enclosures with differentially heated vertical walls and aspect ratios between zero and unity
High Rayleigh number convection in a rectangular cavity with insulated horizontal surfaces and differentially heated vertical walls was analyzed for an arbitrary aspect ratio smaller than or equal to unity. Unlike previous analytical studies, a systematic method of solution based on linearization technique and analytical iteration procedure was developed to obtain approximate closed-form solutions for a wide range of aspect ratios. The predicted velocity and temperature fields are shown to be in excellent agreement with available experimental and numerical data.
Kassemi, Siavash A.↗
Surface deformation, gravity and the geoid from a three-dimensional convection model at low Rayleigh numbers
The two-dimensional surface deformation, gravity field, and geoid are calculated from the temperature fields of a number of numerical models of constant-viscosity three-dimensional convective flows, heated from within and from below. The surface deformation and the geoid are insensitive to the short-wavelength features of the temperature variation; the gravity field is less smooth, although it still does not clearly indicate the narrowness of the upwelling and downwelling regions at large Rayleigh numbers. The results indicate that the convective upwelling beneath regions showing positive geoid and residual depth anomalies is more localized than the horizontal extent of these features would suggest. The results are relevant to attempts to understand the behavior of the earth mantle.
Craig, Claire Harvey↗
A spectral method determination of the first critical Rayleigh number for a low-Prandtl number crystal melt in a cylindrical container
The onset of laminar axisymmetric Rayleigh-Benard convection is investigated for a low-Prandtl number liquid metal in a cylindrical container. All surfaces are considered to be solid and no-slip. Two separate cases are examined for the thermal boundary conditions at the side wall, one with conducting and the other with insulated surface. The governing Boussinesq system is first perturbed and then simplified by introducing a Stokes stream function. Subsequently, a Chebyshev Galerkin spectral model is employed to reduce the simplified system to a system of first-order nonlinear ordinary differential equations. A local stability analysis determines the two values of the first critical Rayleigh number, Ra(sub cl), for the insulated and conducting side walls. As expected, the conducting Ra(sub cl) value of 2882.5 obtained from the present approach exceeded the corresponding insulated Ra(sub cl) value of 2331.6. For the insulated case, an earlier study using a different numerical approach suggests that Ra(sub cl) = 2261.9, while an experimental study measured Ra(sub cl) = 2700.
Dietz, C. M., Jr.↗
A Numerical Study of High Rayleigh Number Thermal Convection in
Two-dimensional convective flows in shallow cavities with adiabatic.
mathematics numerical study thermal convection Ray↗
Exchange of stability for a column of fluid with variable Rayleigh number under free boundary conditions
The exchange of stabilities is demonstrated for a system with harmonic boundary conditions. The motion of fluid in the presence of temperatures gradients is described. It is shown that this principle holds under free, but not rigid or semirigid, boundary conditions.
Korbly, L.↗
A direct interaction approximation treatment of high Rayleigh number convective turbulence and comparison with experiment
A method is proposed for treating steady-state, thermally driven convection using the full direct interaction approximation to treat the nonlinear energy transfer and a prescribed instability function to model the eenrgy input. The instability function used is the growth rate of the mechanism generating the turbulence. This methodology can be easily generalized and applied to other driving mechanisms. The 1/3 power law form of the N vs. R relation for water is duplicated here and the coefficient is computed, using a two-point closure, to be less than about 0.08.
Hartke, Gregory J.↗
Convection treatment for high Rayleigh number, laminar, natural convection calculation
The problem of natural convection in a square cavity of a Boussinesq incompressible fluid is computed using a consistent central-difference approximation, a consistent second-order upwind scheme, and the hybrid scheme for the energy equation. It was found that, for this problem, the second-order upwind scheme yields satisfactory solutions and that an apparently low-order combination of schemes, with the second-order upwind scheme for the momentum equations and the hybrid scheme for the energy equation, can produce satisfactory predictions for this problem.
Shyy, Wei↗
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↗
Chemical differentiation of a convecting planetary interior: Consequences for a one-plate planet such as Venus
Chemically depleted mantle forming a buoyant, refractory layer at the top of the mantle can have important implications for the evolution of the interior and surface. On Venus, the large apparent depths of compensation for surface topographic features might be explained if surface topography were supported by variations in the thickness of a 100-200 km thick chemically buoyant mantle layer or by partial melting in the mantle at the base of such a layer. Long volcanic flows seen on the surface may be explained by deep melting that generates low-viscosity MgO-rich magmas. The presence of a shallow refractory mantle layer may also explain the lack of volcanism associated with rifting. As the depleted layer thickens and cools, it becomes denser than the convecting interior and the portion of it that is hot enough to flow can mix with the convecting mantle. Time dependence of the thickness of a depleted layer may create episodic resurfacing events as needed to explain the observed distribution of impact craters on the venusian surface. We consider a planetary structure consisting of a crust, depleted mantle layer, and a thermally and chemically well-mixed convecting mantle. The thermal evolution of the convecting spherical planetary interior is calculated using energy conservation: the time rate of change of thermal energy in the interior is equated to the difference in the rate of radioactive heat production and the rate of heat transfer across the thermal boundary layer. Heat transfer across the thermal boundary layer is parameterized using a standard Nusselt number-Rayleigh number relationship. The radioactive heat production decreases with time corresponding to decay times for the U, Th, and K. The planetary interior cools by the advection of hot mantle at temperature T interior into the thermal boundary layer where it cools conductively. The crust and depleted mantle layers do not convect in our model so that a linear conductive equilibrium temperature distribution is assumed. The rate of melt production is calculated as the product of the volume flux of mantle into the thermal boundary layer and the degree of melting that this mantle undergoes. The volume flux of mantle into the thermal boundary layer is simply the heat flux divided by amount of heat lost in cooling mantle to the average temperature in the thermal boundary layer. The degree of melting is calculated as the temperature difference above the solidus, divided by the latent heat of melting. A maximum degree of melting is prescribed corresponding to the maximum amount of basaltic melt that the mantle can initially generate. As the crust thickens, the pressure at the base of the crust becomes high enough and the temperature remains low enough for basalt to transform to dense eclogite.
Parmentier, E. M.↗