Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixed convection”

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.

At least 19 records

A numerical study of Li-SF6 wick combustion - Forced and mixed convective burning

A numerical study is conducted to study Li-SF6 wick diffusion flames under mixed convective burning conditions at a pressure of 0.01 MPa. Both planar and cylindrical wicks are considered. The model is based on a conserved scalar approach. The objective of this study is to assess the effects of particular parameters on the burning rate and heat transfer. The flat-plate solution yields a fuel mass burning rate per unit surface area following the x exp -1/2 dependence of the classical similarity solution, where x is the streamwise distance. Cylindrical wick geometries yield enhanced burning rates over planar wicks. For the case of mixed convective burning, the burning rate results approach either the forced or natural convective burning limits as ambient streamwise velocity is changed. Critical Richardson numbers specifying these burning limits are determined for a given condition. Reducing gravity results in a lower burning rate because the influence of natural convection is diminished. Under reduced gravity of 1/1000 of the sea-level value, mixed convective burning nearly resembles forced convection.

Damaso, R. C.↗

Radiative heat transfer in FLiBe molten salt participating medium in a vertical heated tube under forced and mixed convection laminar flows

The contribution of radiative heat transfer (RHT) to convective heat transfer in a heated vertical pipe with laminar flow and constant wall temperature, for FLiBe molten salt, is investigated computationally. The combined effects of conduction in the fluid, forced convection, buoyancy, temperature-dependent physical properties, and thermal radiation are investigated. The P1 approximation is employed in the discretization of the Radiative Transfer Equation (RTE). The COMSOL Multiphysics software is used to generate the numerical solutions. The theoretical analysis developed here demonstrates that only for intermediate values of the optical thickness, i.e. τ D ~ O 1 , the participating media effects are expected to be important. This analysis is confirmed computationally and a value of τ D ≈ 4 is shown to lead to the highest increase in overall heat transfer behavior. Under forced convection, the Nusselt ratio N u total / N u no - rad reaches a peak value of 1.76 at τ D = 4 and z/D = 200 , and is less than 1.1 for τ D < 0 . 1 and τ D > 60 . Under aiding-flow mixed convection, N u total / N u no - rad reaches a peak of 1.34 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 4 and τ D > 36 . Under opposing-flow mixed convection, Nu total / Nu no - rad reaches a peak of 1.81 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 2 and τ D > 50 . RHT effects on the overall heat transfer are most pronounced in opposing mixed convection, where N u total / N u no - rad is 2.26 observed at z/D = 50 . The sensitivity to wall emissivity is evaluated; in forced convection, the peak N u total / N u no - rad is 1.66 at ε = 0 (reflective wall) and 2.04 at ε = 1 (absorptive wall). Finally, the sensitivity to pipe diameter is also discussed. Overall, for a vertical heated tube, radiative heat transfer effects lead to enhancement of heat transfer by as high as a factor of two, and they depend on the optical thickness of the flow, mixed convection environment ( Gr / R e 2 and direction of flow relative to the gravitational force), surface emissivity, and entrance effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantitative assessment of eddy viscosity rans models for turbulent mixed convection in a differentially heated plane channel

Turbulent mixed convection between two vertical, infinite parallel plates at different temperatures is studied using various two-equation turbulence models. The numerical simulations are performed at a turbulent Reynolds number of Re τ = 150 and a Grashof number of Gr = 9.6 × 10 5 . Comparisons are made against the highly trusted DNS results. Consistent with the DNS approach, the current simulations are performed using constant properties and the Boussinesq approximation to predict the influence of buoyancy. Previous studies have provided assessments of two-equation turbulence models for various scenarios, but often rely on a qualitative “eye” test in order to determine the most appropriate model to predict a given flow. This study aims to provide a new form of quantitative assessment that accounts for both the physics captured by the turbulence model as well as the magnitude of the system response quantities (SRQ) using a modified symmetric mean absolute percent error (SMAPE) method. This method is designed to be approachable to researchers at any level and can be applied to system response quantities from multiple research fields. Uncertainty quantification is also performed to determine the discretization error for each turbulence model. Recommendations are made as to which turbulence models best capture the physics – hydrodynamically and thermally – using both local and global validation metrics. Lastly, a sensitivity analysis is performed on the damping functions used in the most accurate models. This underpins the potential of model developments and adjustments most worth pursuing for buoyant flows. Finally, this framework provides a more physics-based comparative analysis of the selected turbulence models.

42 ENGINEERING↗

Opposed-flow flame spread and extinction in mixed-convection boundary layers

Experimental data for flame spread down thin fuel samples in an opposing, mixed-convection, boundary-layer flow are analyzed to determine the gas-phase velocity that characterizes how the flame reacts as it spreads toward the leading edge of the fuel sample into a thinning boundary layer. In the forced-flow limit where the cube of the Reynolds number divided by the Grashof number, Re exp 3/Gr, is large, L(q)/L(e), where L(q) is a theoretical flame standoff distance at extinction and L(e) is the measured distance from the leading edge of the sample where extinction occurs, is found to be proportional to Re exp n with n = -0.874 and Re based on L(e). The value of n is established by the character of the flow field near the leading edge of the flame. The Re dependence is used, along with a correction for the mixed-convection situation where Re exp 3/Gr is not large, to construct a Damkohler number with which the measured spread rates correlate for all values of Re exp 3/Gr.

Altenkirch, R. A.↗

Aerosols in the study of convective acinar mixing

Convective mixing (CM) refers to the different transport mechanisms except Brownian diffusion that irreversibly transfer inspired air into resident air and can be studied using aerosol bolus inhalations. This paper provides a review of the present understanding of how each of these mechanisms contributes to CM. Original data of the combined effect of stretch and fold and gravitational sedimentation on CM are also presented. Boli of 0.5 microm-diameter particles were inhaled at penetration volumes (V(p)) of 300 and 1200 ml in eight subjects. Inspiration was followed by a 10-s breath hold, during which small flow reversals (FR) were imposed, and expiration. There was no physiologically significant dependence in dispersion and deposition with increasing FR. The results were qualitatively similar to those obtained in a previous study in microgravity in which it was speculated that the phenomenon of stretch and fold occurred during the first breathing cycle without the need of any subsequent FR.

Convection↗

Mixed convection in a horizontal rectangular channel - Experimental and numerical velocity distributions

Mixed convection in a rectangular channel (width/height = 2) with bottom-heated and top-cooled sections is studied by laser Doppler anemometry in nitrogen at Ra = 22,200 and Re = 18.75, 36, and 54. At the lower Re values, symmetry breaking is observed in steady but spatially oscillating flows that prevail over a certain distance from the leading edge of the differentially heated section. Further downstream, unsteady flows are found even for Re = 18.75. Numerical models are used to investigate the effects of adiabatic, conducting (with a conductive-convective heat transfer coefficient), and perfectly conducting side walls; channel tilts and Prandtl number dependence. Good agreement between calculations and experiment is obtained for longitudinal convective roll velocities. The transverse velocities are found to be independent of Re.

Nyce, Thomas A.↗

Computation of three-dimensional mixed convective boundary layer flow

The paper presents the numerical solution of heat and mass transfer during cross-flow (orthogonal) mixed convection. In this class of flow, a buoyancy-driven transport in the vertical direction and a forced convective flow in the horizontal direction results in a three-dimensional boundary layer structure adjacent to the plate. The rates of heat and mass transfer are determined by a combined influence of the two transport processes. The equations for the conservation of mass, momentum, energy, and species concentration were solved along with appropriate boundary conditions to determine the distributions of velocity components, temperature, and concentration across the thickness of the boundary layer at different locations on the plate. Results were expressed in dimensionless form using Reynolds number, Richardson number for heat transfer, Richardson number for mass transfer, Prandtl number, and Schmidt number as parameters. It was found that the transport is dominated by buoyancy at smaller vertical locations and at larger distances away from the forced convection leading edge. Effects of forced convection appeared to be very strong at smaller horizontal distances from the leading edge. The cross stream forced convection enhanced the rate of heat and mass transfer by a very significant amount.

Gadepalli, Prashandt↗

3D modelling of MHD mixed convection flow in a vertical duct with transverse magnetic field and volumetric or surface heating

Understanding phenomena associated with the multiple effects/interactions of the fusion nuclear environment on liquid metal flow is required to correctly design liquid metal (LM) blankets for fusion facilities. These effects are investigated in the present work by numerically simulating 3D LM MHD flow. The simulated geometry consists of a straight, vertical duct which runs perpendicular to a strong, fringing applied magnetic field. There is also a region of applied heating as the primary goal is to explore buoyancy effects in MHD duct flows. Results are presented for both buoyancy assisted (upwards) and buoyancy opposed (downwards) flows in conducting and insulating ducts for a range of Hartmann numbers (Ha) up to 100, Reynolds numbers (Re) from 10 3 to 10 4 and Grashof (Gr) numbers from 10 7 to 10 8 . While increasing Gr or decreasing Re increases buoyancy effects, increasing Ha was shown to increase maximum temperature through turbulence reduction. The extent to which the MHD mixed convection flows are quasi-2D is analyzed and buoyant effects, in competition with electromagnetic forces, are shown to bring about 3D flow features not seen in purely MHD flows. Volumetric nuclear heating with steep gradients is applied to the vertical MHD flows for comparison to flows with surface heating only. Finally, surface heating generates stronger buoyancy effects than volumetric heating of the same total power; however, many of the same phenomena occur. Therefore, surface heating, the only option for lab experiments, can provide indication of the effects of volumetric heating in MHD flows.

42 ENGINEERING↗

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗

Heat transfer and pseudo phase transition for low-Reynolds, mixed-convection channel flow in the supercritical thermodynamic regime

Fluids at supercritical thermodynamic conditions are inherently complex due to the large variations in thermodynamic and transport properties. Multiple recent numerical and experimental investigations indicate ongoing interest for these fluids, especially supercritical CO and supercritical water, for a variety of applications. The potential heat transfer benefits in this regime are directly influenced by the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. Limited studies have been conducted for low-Reynolds/intermediate-Rayleigh numbers mixed-convection channel flow of supercritical water, which is the focus of this study. To investigate the thermally driven hydrodynamic instabilities in this regime, we use a high-order fully-implicit numerical method. As expected, in this channel configuration and flow regime, a competition arises between the forced and the natural convection. Buoyancy, driven by the thermal gradient at the wall, pushes the fluid upwards in the channel, while inertial forces drive the fluid from inlet to outlet. We observe significant differences in the heat transfer coefficient between the top and bottom walls of the channel flow, due to natural convection, and unstable plumes. We further quantify the buoyancy forces along the channel flow through evaluation of the Richardson number. Simulations show that with a higher negative Richardson number, there are more unstable structures along the channel, resulting in enhanced mixing. Furthermore, we demonstrate the existence of gas-like, liquid-like, and the possibility of pseudo-two-phase regimes in the supercritical fluid. The pseudo-two-phase regime occurs when the range of fluid temperatures is crossing the pseudo-critical line, also known as the Widom line. The density of the fluid varies significantly near the Widom line, and further influences the natural convection and heat transfer in the flow.

42 ENGINEERING↗

Time-dependent mixed convection heat transfer from a sphere in a micro-gravity environment

A fundamental problem of interest for crystal growth in micro-gravity applications involves the mixed convection heat transfer from a sphere in a uniform flow of fluid at a differing temperature. Under the combined influence of the imposed free stream as well as an induced buoyancy force due to thermal expansion of the fluid, the heat transfer from the sphere will be different from that of either the pure forced convection flow or the pure free convection flow. For the present study, the method of matched asymptotic expansions is applied to the laminar flow problem of an impulsively heated, impulsively started sphere in an originally quiescent fluid. Time series expansions are developed for the dependent variables by acknowledging the existence of two distinct regions: one, an inner region, near the sphere, in which viscous effects are significant; and two, an outer region in which the fluid may be treated as inviscid. The time series expansions are developed in terms of the Reynolds number and Richardson number (Buoyancy Parameter), and the relevant heat transfer and drag coefficients are calculated and plotted.

Hommel, Mark J.↗

From hydrogen to helium: the spectral evolution of white dwarfs as evidence for convective mixing

ABSTRACT We present a study of the hypothesis that white dwarfs undergo a spectral change from hydrogen- to helium-dominated atmospheres using a volume-limited photometric sample drawn from the Gaia-DR2 catalogue, the Sloan Digital Sky Survey (SDSS), and the Galaxy Evolution Explorer (GALEX). We exploit the strength of the Balmer jump in hydrogen-atmosphere DA white dwarfs to separate them from helium-dominated objects in SDSS colour space. Across the effective temperature range from 20 000 to 9000 K, we find that 22 per cent of white dwarfs will undergo a spectral change, with no spectral evolution being ruled out at 5σ. The most likely explanation is that the increase in He-rich objects is caused by the convective mixing of DA stars with thin hydrogen layers, in which helium is dredged up from deeper layers by a surface hydrogen convection zone. The rate of change in the fraction of He-rich objects as a function of temperature, coupled with a recent grid of 3D radiation-hydrodynamic simulations of convective DA white dwarfs – which include the full overshoot region – lead to a discussion on the distribution of total hydrogen mass in white dwarfs. We find that 60 per cent of white dwarfs must have a hydrogen mass larger than MH/MWD = 10−10, another 25 per cent have masses in the range MH/MWD = 10−14–10−10, and 15 per cent have less hydrogen than MH/MWD = 10−14. These results have implications for white dwarf asteroseismology, stellar evolution through the asymptotic giant branch and accretion of planetesimals on to white dwarfs.

Cunningham, Tim↗

Solution of mixed convection heat transfer from isothermal in-line fins

Transient and steady state combined natural and forced convective flows over two in-line finite thickness fins (louvers) in a vertical channel are numerically solved using two methods. The first method of solution is based on the 'Simple Arbitrary Lagrangian Eulerian' (SALE) technique which incorporates mainly two computational phases: (1) a Lagrangian phase in which the velocity field is updated by the effects of all forces, and (2) an Eulerian phase that executes all advective fluxes of mass, momentum and energy. The second method of solution uses the finite element code entitled FIDAP. In the first part, comparison of the results by FIDAP, SALE, and available experimental work were done and discussed for steady state forced convection over louvered fins. Good agreements were deduced between the three sets of results especially for the flow over a single fin. In the second part and in the absence of experimental literature, the numerical predictions were extended to the transient transports and to the opposing flow where pressure drop is reversed. Results are presented and discussed for heat transfer and pressure drop in assisting and opposing mixed convection flows.

Khalilollahi, Amir↗

On the Departure from Monin–Obukhov Surface Similarity and Transition to the Convective Mixed Layer

Large-eddy simulations are used to evaluate mean profile similarity in the convective boundary layer (CBL). Particular care is taken regarding the grid sensitivity of the profiles and the mitigation of inertial oscillations in the simulation spin-up. The nondimensional gradients Φ for wind speed and air temperature generally align with Monin–Obukhov similarity across cases but have a steeper slope than predicted within each profile. The same trend has been noted in several other recent studies. The Businger-Dyer relations are modified here with an exponential cutoff term to account for the decay in Φ to first-order approximation, yielding improved similarity from approximately 0.05z i to above 0.3z i , where z i is the CBL depth. The necessity for the exponential correction is attributed to an extended transition from surface scaling to zero gradient in the mixed layer, where the departure from Monin–Obukhov similarity may be negligible at the surface but becomes substantial well below the conventional surface layer height of 0.1 z i .

54 ENVIRONMENTAL SCIENCES↗

A Clue to the Extent of Convective Mixing Inside Massive Stars: The Surface Hydrogen Abundances of Luminous Blue Variables and Hydrogen-Poor Wolf-Rayet Stars

Interior layers of stars that have been exposed by surface mass loss reveal aspects of their chemical and convective histories that are otherwise inaccessible to observation. It must be significant that the surface hydrogen abundances of luminous blue variables (LBVs) show a remarkable uniformity, specifically X(sub surf) = 0.3 - 0.4, while those of hydrogen-poor Wolf-Rayet (WN) stars fall, almost without exception, below these values, ranging down to X(sub surf) = 0. According to our stellar model calculations, most LBVs are post-red-supergiant objects in a late blue phase of dynamical instability, and most hydrogen-poor WN stars are their immediate descendants. If this is so, stellar models constructed with the Schwarzschild (temperature-gradient) criterion for convection account well for the observed hydrogen abundances, whereas models built with the Ledoux (density-gradient) criterion fail. At the brightest luminosities, the observed hydrogen abundances of LBVs are too large to be explained by any of our highly evolved stellar models, but these LBVs may occupy transient blue loops that exist during an earlier phase of dynamical instability when the star first becomes a yellow supergiant. Independent evidence concerning the criterion for convection, which is based mostly on traditional color distributions of less massive supergiants on the Hertzsprung-Russell diagram, tends to favor the Ledoux criterion. It is quite possible that the true criterion for convection changes over from something like the Ledoux criterion to something like the Schwarzschild criterion as the stellar mass increases.

Stothers, Richard B.↗

Long‐Term Observations of Turbulence Vertical Velocity Spectra in a Convective Mixed Layer: Dependence on Land‐Surface Forcing in the U.S. Southern Great Plains

Doppler lidar vertical velocity retrievals were analyzed for the scale and structure of mixed-layer turbulence over a 7-year period, on fair-weather warm season days (442 cases) in the U.S. Southern Great Plains. Based on the hypothesized influence of land-surface forcing and updraft size on convective boundary layer clouds, a spectral analysis was performed to quantify effects of surface forcing (surface buoyancy flux, friction velocity, and evaporative fraction) on turbulence scale. Significant (order-of-magnitude) variations in spectral density were found in the energy-production subrange and mesoscale regimes. Integral scale decreased with increasing buoyancy flux (and Monin-Obukhov stability parameter), while spectral density in the energy-production subrange increased, implying a transition to buoyancy-driven cellular structures with narrower updrafts. The influence of stability parameter was limited to the neutral to convective transition, and could not explain the wide variation in spectral density in the mesoscale regime. However, high friction velocity (u* > 0.5 m s –1 ) was associated with larger integral scale (updraft width), and a greater portion of spectral density in the mesoscale regime. Lidar profiles and radar reflectivity for individual cases revealed evidence of roll and wave-like structures contributing to mesoscale variability on high friction velocity days, suggesting shear instabilities on days with wetter land surface conditions and smaller buoyancy flux. The role of friction velocity emphasizes the multivariate nature of surface influences on turbulence, and implies that variance-based turbulence closures may not be adequate for capturing effects of surface forcing on updraft size.

54 ENVIRONMENTAL SCIENCES↗