Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NATURAL 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 145 records · Page 8

Simulation of fluid flows during growth of organic crystals in microgravity

Several counter diffusion type crystal growth experiments were conducted in space. Improvements in crystal size and quality are attributed to reduced natural convection in the microgravity environment. One series of experiments called DMOS (Diffusive Mixing of Organic Solutions) was designed and conducted by researchers at the 3M Corporation and flown by NASA on the space shuttle. Since only limited information about the mixing process is available from the space experiments, a series of ground based experiments was conducted to further investigate the fluid dynamics within the DMOS crystal growth cell. Solutions with density differences in the range of 10 to the -7 to 10 to the -4 power g/cc were used to simulate microgravity conditions. The small density differences were obtained by mixing D2O and H2O. Methylene blue dye was used to enhance flow visualization. The extent of mixing was measured photometrically using the 662 nm absorbance peak of the dye. Results indicate that extensive mixing by natural convection can occur even under microgravity conditions. This is qualitatively consistent with results of a simple scaling analysis. Quantitave results are in close agreement with ongoing computational modeling analysis.

Roberts, Gary D.↗

Flow Effects on the Flammability Diagrams of Solid Fuels: Microgravity Influence on Ignition Delay

The possibility of an accidental fire in space-based facilities is a primary concern of space exploration programs. Spacecraft environments generally present low velocity air currents produced by ventilation and heating systems (of the order of 0.1 m/s), and fluctuating oxygen concentrations around that of air due to CO2 removal systems. Recent experiments of flame spread in microgravity show the spread rate to be faster and the limiting oxygen concentration lower than in normal-gravity. To date, there is not a material flammability-testing protocol that specifically addresses issues related to microgravity conditions. The present project (FIST) aims to establish a testing methodology that is suitable for the specific conditions of reduced gravity. The concepts underlying the operation of the LIFT apparatus, ASTM-E 1321-93, have been used to develop the Forced-flow Ignition and flame-Spread Test (FIST). As in the LIFT, the FIST is used to obtain the flammability diagrams of the material, i.e., graphs of ignition delay time and flame spread rate as a function of the externally applied radiant flux, but under forced flow rather than natural convection conditions, and for different oxygen concentrations. Although the flammability diagrams are similar, the flammability properties obtained with the FIST are found to depend on the flow characteristics. A research program is currently underway with the purpose of implementing the FIST as a protocol to characterize the flammability performance of solid materials to be used in microgravity facilities. To this point, tests have been performed with the FIST apparatus in both normal-gravity and microgravity conditions to determine the effects of oxidizer flow characteristics on the flammability diagrams of polymethylmethacrylate (PMMA) fuel samples. The experiments are conducted at reduced gravity in a KC- 135 aircraft following a parabolic flight trajectory that provides up to 25 seconds of low gravity. The objective of the experiments is to obtain data of ignition delay and flame spread rate at low flow velocities (0.1 to 0.2 m/s), which cannot be obtained under normal gravity because of the natural convection induced flows (approx. 0.5 m/s). Due to the limited reduced gravity time, the data can only be obtained for high radiant fluxes, and are consequently limited in scope. These tests do, however, provide insight into the flammability diagram characteristics at low velocity and reduced gravity, and also into the implications of the flow-dependence of the flammability properties under environments similar to those encountered in space facilities.

Cordova, J. L.↗

Thermovelocimetric characterization of liquid metal convection in a rotating slender cylinder

Rotating turbulent convection occurs ubiquitously in natural convective systems encompassing planetary cores, oceans and atmospheres, as well as in many industrial applications. While the global heat and mass transfer of water-like rotating Rayleigh–Bénard convection is well-documented, the dynamics in low-Prandtl-number liquid metals remain less understood. In this study, we experimentally investigate rotating Rayleigh-Bénard convection in liquid gallium (Prandtl number Pr ≈ 0.027) within a slender cylinder (diameter-to-height aspect ratio Γ = D/H =1/2) using novel thermovelocimetric diagnostic techniques that integrate simultaneous multi-point thermometry and ultrasonic Doppler velocity measurements. Our results reveal the formation of a stable, global-scale azimuthal wavenumber m = 2 quadrupolar vortex at low supercriticality. We propose that enhanced wall modes facilitated by the slender cylinder geometry interact with the bulk flow to create these large-scale axialized vortices. Furthermore, our findings imply a distinct scaling behavior for the wall-mode precession frequency in liquid metals, extending previous results obtained for moderate-Pr fluids. In conclusion, this provides new insights into wall-bulk coupling mechanisms of low-Pr rotating convective turbulence.

Magnetohydrodynamics↗

Heat Entrapment Effects Within Liquid Acquisition Devices

We introduce a model problem to address heat entrapment effects or the local accumulation of thermal energy within liquid acquisition devices. We show that the parametric space consists of six parameters, namely the Rayleigh and Prandtl numbers, the aspect ratio, and heat flux ratios for the bottom, side, and top boundaries of the enclosure. For the range of Ra considered 1 to 10(sup 9), beyond Ra on the order of 10(sup 5), convective instability is the dominant mode of convection in comparison to natural convection. The flow field transitions to asymmetric modes at Ra on the order of 10(sup 7). Direct numerical simulation of a large geometric length scale prototype for Ra on the order of 10(sup 9) shows that the flow field evolves from small wavelength instability which gives rise to nonlinear growth of thermals, propagation of the instability occurs via growth of secondary and tertiary modes, and a travelling wave mode occurs prior to asymmetry. The effect of a large aspect ratio is to increase the number of modes in the vertical direction. Due to the slow diffusion of heat in the prototype, asymptotic states are not readily attained, we show that dynamical similarity can be used for a model which allows the attainment of asymptotic states and that transition to a chaotic state occurs for Ra on the order of 10(sup 9) via a broadband power spectrum. These dynamical events show that for the baseline condition in which heat is absorbed from background laboratory environment, higher heat flux is absorbed at the top and bottom boundaries of the enclosure than a nominal value of 34.9 ergs per square centimeter -second.

Duval, W. M. B.↗

Combustion of Interacting Droplet Arrays in a Microgravity Environment

This research program involves the study of single droplets and linear arrays of droplets in weakly-buoyant and non-bouyant environments. The primary purpose of the single droplet work was to (1) provide a data base from which to compare droplet array results and (2) to correlate the effects of buoyancy on flame shape. Traditionally convective effects in droplet combustion are represented in terms of the Reynolds number, Re, for forced convection and the Grashof number, Gr, for natural convection. Typically, corrections to the burning rate constant for convective effects are written in terms of Re or Gr(exp 1). The Stefan velocity is not included in these correlations, even though from purely physical reasons, one would expect it to be important, especially at higher burning rates. The flame distortion due to convective effects is less documented quantitatively. Kumagai and Isoda do predict flame shape in natural and forced convective flow fields. Their focus, however, was to predict the actual flame dimensions. Law and co-workers used reduced pressure, high oxidizer ambients to obtain spherical flames. This implies that buoyant flows were reduced at the low pressures, as indicated by a very small Grashof number. Ross et al, however, using scaling arguments showed that reducing the pressure does not have a large effect on the magnitude of the buoyant velocity. Struk et al showed elongated flame shapes during simulated (porous sphere) droplet combustion. The elongation of the flames was due to residual gravity levels aboard the reduced gravity aircraft on which the experiments were conducted. These flame shapes, as well as some data from the literature were interpreted based on a dimensionless grouping called the sphericity parameter, Sp. Sp is the ratio of a characteristic computed buoyant velocity to the Stefan velocity at the flame front. One purpose of the droplet arrays work is to extend the database and theories that exist for single droplets into the regime where droplet interactions are important. The eventual goal being to use the results of this work as inputs to models on spray combustion where droplets seldom burn individually; instead the combustion history of a droplet is strongly influenced by the presence of the neighboring droplets. Recently, Annamali and Ryan have summarized he current status of droplet array, cloud and spray combustion. A number of simplified theories led numerical studies of droplet vaporization/combustion where multiple droplet effects are present are now available. These theories all neglect the effect of buoyancy. Experimentally, most studies to date suffer the effects of buoyancy. It is the dominant transport mechanism in the problem. Only the works of Law and co-worker and more recently by Mikami et al were performed in an environment where buoyancy effects were small. Law and co-workers were limited to high oxygen index, low pressure ambient environments since there studies were conducted in normal gravity.

Dietrich, Daniel L.↗

Influence of Reservoir Convection on Heat Extraction with Closed-Loop Geothermal Systems

Over the past few years, advancements in closed-loop geothermal systems (CLGS), also called advanced geothermal systems (AGS), have sparked a renewed interest in these types of designs. CLGS have certain advantages over traditional and enhanced geothermal systems (EGS), including not requiring in-situ reservoir permeability, conservation of the circulating fluid, and allowing for different fluids, including working fluids directly driving a turbine at the surface. CLGS may be attractive in environments where water resources are limited, rock contaminants must be avoided, and stimulation treatments are not available (e.g., due to regulatory or technical reasons). Despite these advantages, CLGS have some challenges, including limited surface area for heat transfer and requiring long wellbores and laterals to obtain multi-MW output in conduction-only reservoirs. CLGS have been investigated in conduction-only systems. In this paper, we explore the impact of both forced and natural convection on the levels of heat extraction with a CLGS deployed in a hot wet rock reservoir. We bound potential benefits of convection by investigating liquid reservoirs over a range of natural and forced convective coefficients. Additionally, we investigate the effects of permeability, porosity, and geothermal temperature gradient in the reservoir on CLGS outputs. Reservoir simulations indicate that reservoir permeabilities of at least ~100 mD are required for natural convection to increase the heat output with respect to a conduction-only scenario. The impact increases with increasing reservoir temperature. When subject to a forced convection flow field, Darcy velocities of at least 10-7 m/s are required to obtain an increase in heat output.

advanced geothermal systems↗

Experimental and Numerical Investigation of Buoyancy Driven Convection During PDAMNA Thin Film Growth

This paper presents results from numerical simulations as well as laboratory experiments of buoyancy driven convection in an ampoule under varying heating and gravitational acceleration loadings. The modeling effort in this work resolves the large scale natural convective motion that occurs in the fluid during photodeposition of polydiacetelene films which is due to energy absorbed by the growth solution from a UV source. Consequently, the growth kinetics of the film are ignored in the model discussed here, and also a much simplified ampoule geometry is considered. The objective of this work is to validate the numerical prediction on the strength and structure of buoyancy driven convection that could occur under terrestrial conditions during nonlinear optical film growth. The validation is used to enable a reliable predictive capability on the nature and strength of the convective motion under low gravity conditions. The ampoule geometry is in the form of a parallelepiped with rectangular faces. The numerical results obtained from the solution to the Boussinesq equations show that natural convection will occur regardless of the orientation of the UV source with respect to the gravity vector. The least strong convective motion occurred with the UV beam directed at the top face of the parallelepiped. The strength of the convective motion was found to be almost linearly proportional to the total power of the UV source. Also, it was found that the strength of the convective motion decreased linearly with the gravity due to acceleration. The pattern of the convective flow on the other hand, depended on the source location.

Antar, Basil N.↗

Convective heat transfer to low-temperature fluids

Research into forced and natural convection processes in low-temperature (cryogenic) fluids is reviewed with primary emphasis on forced convection. Boundaries of the near-critical region are defined, fluid properties near the critical state are discussed, and heat-transfer processes around the critical point are described. The thermodynamics of the critical point is analyzed together with transport properties of a near-critical fluid, and the quantum states of low-temperature molecular hydrogen (para and ortho) are discussed. Experimental work on heat transfer in free, natural, and forced convection systems is briefly summarized. Graham's (1969) penetration model for near-critical fluids is outlined, near-critical heat transfer is discussed in relation to conventional geometric effects, and the effects of curvature on the properties of near-critical hydrogen are noted. Theoretical considerations in free and forced convection are examined.

Graham, R. W.↗

Numerical studies of convective heat transfer in an inclined semiannular enclosure

Natural convection heat transfer in a two-dimensional differentially heated semiannular enclosure is studied. The enclosure is isothermally heated and cooled at the inner and outer walls, respectively. A commercial software based on the SIMPLER algorithm was used to simulate the velocity and temperature profiles. Various parameters that affect the momentum and heat transfer processes were examined. These parameters include the Rayleigh number, Prandtl number, radius ratio, and the angle of inclination. A flow regime extending from conduction-dominated to convection-dominated flow was examined. The computed results of heat transfer are presented as a function of flow parameter and geometric factors. It is found that the heat transfer rate attains a minimum when the enclosure is tilted about +50 deg with respect to the gravitational direction.

Wang, Lin-Wen↗

Single-phase, natural circulation annular flow measurements for cartridge loop irradiation experiments

The nuclear industry is increasingly considering cartridge-style experiments for irradiation testing of advanced reactor fuels and materials under flowing conditions. Cartridge loops do not require the extensive support infrastructure that are necessary for external flow loops and minimize the possibility of coolant solidification over the long distance from the reactor to the external facilities. However, there is a general lack of quality flow data for internally heated fluids in an annular configuration representative of a cartridge-type irradiation experiment, particularly one with natural circulation. To address this data need, a series of experiments was conducted to measure the natural circulation flow rates of pressurized water in a sealed, internally heated vessel with annular flow conditions that represent a molten salt or sodium cartridge loop. Temperatures and flow rates were measured under steady-state and transient conditions. Here, this paper describes the facility, methods, and results of the experiments, including the determination of nondimensional parameters. A simple 1D model of the natural convection flow rates agrees well with the experimental results. Applying this model to simulate a liquid salt cartridge experiment predicts that natural circulation flow might be able to provide liquid salt Reynolds numbers similar to those of some molten salt reactor concepts at relevant power densities.

42 ENGINEERING↗

Solidification Using a Baffle in Sealed Ampoules (SUBSA)

Solidification Using a Baffle in Sealed Ampoules (SUBSA) will be the first materials science experiment conducted in the Microgravity Science Glovebox (MSG) Facility at the International Space Station (ISS) Alpha. The launch is schedule for May 31, 2002. Using the specially developed furnace, 10 Te and Zn-doped single crystals of InSb will be directionally solidified in microgravity. A key goal of the SUBSA investigation is to (i) clarify the origin of the melt motion in space laboratories and (ii) to reduce the magnitude of the melt motion to the point that it does not interfere with the transport phenomena. These goals will be accomplished through a special ampoule and furnace design. A disk-shaped baffle, positioned close to the freezing front, is used to reduce melt motion. Furthermore, the solidification will be visualized by using a transparent furnace, with a video camera, continuously sending images to the earth. This allows detection of bubbles and melt de-wetting that could cause surface tension driven convection. In preparation for the space experiments, 30 ground-based experiments were conducted. The results of ground based tests and numerical modeling will be presented. Based on numerical modeling, 12 mm 1D silica ampoules were selected. The small diameter ampoule favors closer placement of the baffle to the interface, without excessive radial segregation caused by forced convection while providing more damping of natural convection. The parts in the silica ampoule include 2 carbon springs made by Energy Science Laboratories, Inc., a pyrocarbon-coated graphite cylinder, pyrocarbon-coated graphite a baffle with the shaft and the InSb charge with the seed crystal grown by W.A. Bonner of Crystallod Inc.

Marin, C.↗

Numerical modeling of enclosure convection

A numerical study on the steady and unsteady natural convection in two-dimensional rectangular enclosures has been performed by a time-accurate ADI finite difference scheme. The study covered a range of Rayleigh numbers between 1000 and 10 to the 7th, aspect ratios between 0.2 and 10.0, and tilt angles between -90 (heating from bottom) and +90 deg (heating from top). Various Prandtl numbers have been studied, but only the results of water (Pr = 7.0) are reported here due to space limitations. The physics revealed, however, includes the convection phenomena and the Rayleigh-Benard stability, as well as the combined mechanism of these two. The onset of secondary cells is determined by using a velocity map, which is simpler and cleaner, instead of a streamline plot. The critical Ra number for the occurrence of these secondary cells is shown to be lower than can be detected by experimental studies. On the Rayleigh-Benard stability part, a second transition from stable single-cell convection to periodic multicellular convection is disclosed.

Duh, J. C.↗

Flow dynamics and heat transfer in simplified battery energy storage systems with heated battery modules

Large-scale energy storage systems (ESSs) composed of batteries show promise in addressing current energy challenges, but dissipation of generated heat is important. Here, this paper focuses on buoyant convective flows in simplified ESS battery racks. Natural convection is not generally the primary cooling strategy but can be important in abnormal scenarios where there is module overheat or potentially thermal runaway. We use computational fluid dynamics to investigate the flow dynamics and heat transfer mechanisms in a simplified parameterized rack design. Despite its simplicity, this configuration produces many of the relevant features expected in real ESSs without details of module geometry or hardware, allowing broad conclusions independent of manufacture-specific designs. We start by providing visualizations of the flowfield and measurements of entrainment, heat flux, and pressure. To characterize the dependence on the system parameters, we develop an integral-scale analysis of the average temperature equation to highlight the dominant source terms. We use results from this analysis to derive a steady network model composed of simple algebraic expressions to provide first-order predictions of entrainment through the rack. The network model leads to a linear scaling of the Reynolds number based on convective mass flux with respect to the Grashof number based on the heat source. We deduce empirical relationships that relate the heat exchanged between modules using a surface-averaged Nusselt number as a function of the local Reynolds and Rayleigh numbers. Lastly, we investigate how space between the modules and rack in the spanwise direction creates flow bypass, resulting in different flow pathways.

Battery thermal management↗

Numerical studies of convective heat transfer in an inclined semi-annular enclosure

Natural convection heat transfer in a two-dimensional differentially heated semiannular enclosure is studied. The enclosure is isothermally heated and cooled at the inner and outer walls, respectively. A commercial software based on the SIMPLER algorithm was used to simulate the velocity and temperature profiles. Various parameters that affect the momentum and heat transfer processes were examined. These parameters include the Rayleigh number, Prandtl number, radius ratio, and the angle of inclination. A flow regime extending from conduction-dominated to convection-dominated flow was examined. The computed results of heat transfer are presented as a function of flow parameter and geometric factors. It is found that the heat transfer rate attains a minimum when the enclosure is tilted about +50 deg with respect to the gravitational direction.

Wang, L.-W.↗

Determination of blade-to-coolant heat-transfer coefficients on a forced-convection, water-cooled, single-stage turbine

Blade-to-coolant convective heat-transfer coefficients were obtained on a forced-convection water-cooled single-stage turbine over a large laminar flow range and over a portion of the transition range between laminar and turbulent flow. The convective coefficients were correlated by the general relation for forced-convection heat transfer with laminar flow. Natural-convection heat transfer was negligible for this turbine over the Grashof number range investigated. Comparison of turbine data with stationary tube data for the laminar flow of heated liquids showed good agreement. Calculated average midspan blade temperatures using theoretical gas-to-blade coefficients and blade-to-coolant coefficients from stationary-tube data resulted in close agreement with experimental data.

COOLING - GAS TURBINE SYSTEMS↗

Progress in incompressible Navier-Stokes computations for propulsion flows and its dual-use applications

Development of an incompressible Navier-Stokes solution procedure was performed for the analysis of a liquid rocket engine pump components and for the mechanical heart assist devices. The solution procedure for the propulsion systems is applicable to incompressible Navier-Stokes flows in a steadily rotating frame of reference for any general complex configurations. The computer codes were tested on different complex configurations such as liquid rocket engine inducer and impellers. As a spin-off technology from the turbopump component simulations, the flow analysis for an axial heart pump was conducted. The baseline Left Ventricular Assist Device (LVAD) design was improved by adding an inducer geometry by adapting from the liquid rocket engine pump. The time-accurate mode of the incompressible Navier-Stokes code was validated with flapping foil experiment by using different domain decomposition methods. In the flapping foil experiment, two upstream NACA 0025 foils perform high-frequency synchronized motion and generate unsteady flow conditions for a downstream larger stationary foil. Fairly good agreement was obtained between unsteady experimental data and numerical results from two different moving boundary procedures. Incompressible Navier-Stokes code (INS3D) has been extended for heat transfer applications. The temperature equation was written for both forced and natural convection phenomena. Flow in a square duct case was used for the validation of the code in both natural and forced convection.

Kiris, Cetin↗

Uncertainty estimation of bifurcated solutions in the Rayleigh–Bénard problem for advanced nuclear reactors applications

Multiphysics models of nuclear reactors frequently comprise nonlinear systems of equations. The nonlinear nature of these models could lead to solution bifurcations, where a small change in a certain parameter, e.g., the thermophysical properties of the coolant, can lead to a sudden change in the system’s behavior. At the point in parameter space where this happens, called a critical point, the Jacobian matrix of the model’s nonlinear operator becomes singular potentially permitting multiple solutions to coexist. In this paper, we perform uncertainty estimation (UE) in a parameter range that includes bifurcated solutions within the context of Rayleigh–Bénard problem. We perform this analysis assuming uncertain temperature difference, and tilt angle for the iterative solution algorithm with a unit Prandtl number (Pr = 1). Also, we perform this analysis under uncertain thermophysical properties for both FLiBe molten salt and liquid sodium as working fluid. We deploy two approaches to compute statistical moments for the resulting distributions of selected flow-field variables. The first approach is the blind computation of the mean and the standard deviation without any consideration of solution bifurcation, while the second approach utilizes k-means clustering to cluster each branch’s solutions together and compute separate statistical moments for each branch. The statistical distributions are obtained by perturbing the selected parameters about nominal values that correspond to a solution on one of the valid branches, and that solution is used as initial guess for the iterative solution algorithm. We found that perturbation of any parameter when its nominal value is close to its critical point always leads to branch jumping, i.e., the iterations converge to a solution on a branch different from the branch of the initial guess. This produces a statistical ensemble comprised of fundamentally different solutions leading to wrong mean values and uncertainty estimates, whereas clustering provides an efficient way to deal with this type of computation. This work is important for developing Gen IV nuclear systems because many of these systems rely on natural convection for cooling especially in accident conditions.

97 - MATHEMATICS AND COMPUTING↗