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 469 records · Page 26

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↗

The relationships between high latitude convection reversals and the energetic particle morphology observed by the Atmosphere Explorer

Simultaneous measurements of the auroral zone particle precipitation and the ion convection velocity by Atmosphere Explorer show a consistent difference between the location of the poleward boundary of the auroral particle precipitation and the ion convection reversal. The difference of about 1.5 degrees of invariant latitude is such that some part of the antisunward convection lies wholly within the auroral particle precipitation region. The nature of the convection reversals within the precipitation region suggests that in this region the convection electric field is generated on closed field lines that connect in the magnetosphere to the low latitude boundary layer.

Heelis, R. A.↗

Topographic forcing of supercritical convection in a porous medium such as the oceanic crust

Experiments carried out with a Hele-Shaw cell indicate that above the critical Rayleigh number, topography does not control the convection pattern, except when the topographic wavelength is comparable to the depth of water penetration. Scaled to the depth of the layer, the convective wavenumbers are restricted to values between 2.5 and 4.8 - a range that brackets pi, the natural wavenumber for convection in a porous slab with planar, isothermal, impermeable boundaries. Topographies within this range are found to control the circulation pattern perfectly, with downwelling under valleys and upwelling aligned with topographic height. Other topographies do not force the pattern, although it is noted that in some cases the convection wavenumber may be a harmonic of the topographic wavenumber. It is pointed out that unforced circulation cells wander and vary in size because they are not locked to the topography.

Hartline, B. K.↗

The convective liquidus in a solidifying magma chamber - A fluid dynamic investigation

An experimental study of the solidification of a layer of paraffin cooled from above is reported. When the fluid is superheated, convection sets in at the beginning of cooling but rapidly decreases in intensity. Once the layer has lost its superheat, convection ceases and further cooling is by conduction. The interior temperature then remains constant until encountered by the upper crust. The convective liquidus is defined as the temperature threshold below which convection is very weak or nonexistent. In natural basaltic systems the convective liquidus, although strictly unknown, may be very close to the true liquidus. Rapid convection may therefore not be a dominant process during the crystallization of many magmna chambers; instead, convection is part of an overall intimate balance between phase equilibria, crystal growth, and heat transfer.

Brandeis, Genevieve↗

Volatiles, rheology, and mantle convection: Comparing Earth, Venus, and Mars

Silicate rheology is controlled in part by volatile content. The variation of viscosity with position in the mantle will influence the nature of mantle convection; hence, modeling mantle convection and its effect on surface observables such as the geoid places constraints on the viscosity structure of a planet's mantle and may indirectly constrain the volatile distribution. Models of viscous mantle flow and the Earth's geoid indicate that there is roughly a two order of magnitude variation in viscosity between the upper and lower mantles, although there is some disagreement over the depth of the viscosity minimum in the upper mantle. Some studies of post-glacial rebound also support such a viscosity contrast between the upper and lower mantles. On Venus, several highland regions appear to be supported by mantle plumes. Modeling of the geoid and topography of these regions indicates that if these features are plume-related, then the mantle of Venus can not have an Earth-like low viscosity zone in its upper mantle. On Mars, the Tharsis volcanic province has alternatively been explained as supported either by mantle convection or by flexure of a thick lithosphere. If the convective model is correct, then the large geoid anomaly requires that Mars can not have a low viscosity zone in its upper mantle.

Kiefer, Walter S.↗

High Marangoni number convection in a square cavity

Steady thermocapillary flows in a square cavity are computed using a finite difference procedure. Accurate numerical solutions are obtained for as high a Reynolds number as possible, in order to characterize the nature of strongly convective flows of this general class. Only the case of capillary number going to zero is considered, so that the free surface is assumed flat at the leading order. Surface deflections are computed by domain perturbation. Boundary layer formation at different values of the Prandtl number is observed at large values of the Marangoni number and Re. The numerical results are used to infer the relevant scalings of a consistent boundary layer picture of the flow, valid asymptotically as Re goes to infinity.

Zebib, A.↗

Direct numerical simulations of turbulent convection with a variable gravity and Keplerian rotation

Thermal convection was proposed as a possible mechanism for generation and maintenance of turbulence in the inner accretion disk regime of the primordial solar nebula. It is of fundamental interest to design experiments with the basic physical features of the solar nebula conditions cannot be produced in the laboratory, numerical simulations of hydrodynamic flows, which have been very successful in describing aerodynamic flows, can be suitable modified to provide experimental data for solar nebula modelling. The goals are to modify an extant, proven hydrodynamics code with the most important features of the solar nebula and other thin accretion disks: bouyancy terms to generate convection, internal heating representing the release of gravitational potential energy, a variable gravity linearly proportional the the distance from the vertical midplane due to centrifugal balance, rapid rotation with axis aligned with gravity, and Keplerian rotational shear; to determine the effect that these features have on the turbulent convection by introducing them individually and to determine the cumulative nature of the turbulent convection for accretion disk conditions; and to model the convection and the turbulence. In this manner, prior solar nebula models can be tested and their deficiencies rectified.

Cabot, William H.↗

Exploiting the Temperature Dependence of Magnetic Susceptibility to Control Convection in Fundamental Studies of Solidification Phenomena

It is well known that convection is a dominant mass transport mechanism when materials are solidified on Earth's surface. This convection is caused by gradients in density (and therefore gravitational force) that are brought about by gradients in temperature, composition or both. Diffusion of solute is therefore dwarfed by convection and the study of fundamental parameters, such as dendrite tip shape and growth velocity in the absence of convection is nearly impossible. Significant experimental work has therefore been carried out in orbiting laboratories with the intent of minimizing convection by minimizing gravity. One of the best known experiments of this kind is the Isothermal Dendritic Growth Experiment (IDGE), supported by NASA. Naturally such experiments are costly and one objective of the present investigation is to develop an experimental method whereby convection can be- halted, in solidification and other experiments, on the surface. A second objective is to use the method to minimize convection resulting from the residual accelerations suffered by experiments in microgravity. The method to be used to minimize convection relies on the dependence of the magnetic susceptibility of a fluid on temperature or composition (whichever is driving convection). All materials experience a force when placed in a magnetic field gradient. The direction and magnitude of that force depend on the magnetic susceptibility of the material. Consequently the force will vary if the susceptibility varies with temperature or composition. With a magnetic field gradient in the right direction (typically upward) and of the right magnitude, this variation in the magnetic force can be made to exactly cancel the variation in the gravitational force. Expressed another way, normal buoyancy is exactly countered by a "magnetic buoyancy". To demonstrate the principle, a solution of MnC12 in water has been used. First the variation of the susceptibility of this paramagnetic solution with temperature and concentration was measured. Then a "cell", containing this solution and 50mm long by 15mm high by 155mm wide, was placed in a superconducting magnet at Marshall Space Flight Center. The magnetic field was measured at various positions within the bore of the magnet using a Hall effect probe. In this way, a position was found where the magnetic field gradient was predominantly upward; the magnitude of the gradient could then be adjusted by adjusting the current of the magnet. The ends of the cell consisted of machined copper blocks maintained at controlled temperatures by circulating water from constant temperature baths. The walls of the cell were of rectangular section glass tubing so that the cell contents could be seen. Velocities arising from thermal gradients within the cell were measured by particle image velocimetry (PIV). Particles used for this purpose were silver-coated hollow glass spheres of micrometers diameter and nearly the same density as the solution. A central vertical plane of the cell was illuminated by a laser beam passing through a cylindrical lens. Digital images of the particles were captured on a CCD camera and fed to a computer so that frame-to-frame movements of particles traveling with the fluid were captured. These images were employed to compute velocity maps using commercial PIV software. In a typical experiment the cold end of the cell was maintained at 10C and the warm end at 30 C. With no current in the magnet, i.e.- with natural -convection allowed to occur, the fluid was observed to circulate with an average speed of approximately 0.3 millimeters per second. It was visually apparent that this circulation was diminished as the current was increased. At currents of approximately 20A the flow was halted, to within the precision of the PIV measurements. At yet higher currents the convection was reversed with the hotter solution sinking and the cooler solution rising. At 40A this reversed convection had speeds averaging 0.43 millimeters per second. The measurements of susceptibility and density allow an estimate of the field gradient necessary to halt convection in the experiment. That estimate was 7.8T (squared) per meter and the convection was observed to halt in the magnet at a current giving 7.21T (squared) per meter from the magnetic field measurements. Calculations of the flow have been carried out using the computational fluid dynamics software FLUENT and show good agreement with the measurements.

Seybert, C. D.↗

Structural modeling and dynamic analysis of condensate storage tanks in nuclear power plants

Condensate storage tanks (CSTs) in nuclear power plants (NPPs) are classified as critical equipment capable of surviving strong shaking in a design basis earthquake to assure the ability to subsequently provide cooling water. Structural modeling and dynamic analysis of CSTs are complicated due to fluid-structure interaction (FSI) and coupling issues, which could have significant effects on seismic response. The capability of a number of approximate models to capture seismic response is investigated using dynamic analysis results from detailed finite element models. Modal analyses and time history analyses are carried out using both 2D and 3D models to predict the dynamic behavior of CSTs as a function of ground motion intensity. Although 2D simplified models can be used to quickly evaluate the dynamic response of CSTs when there is substantial margin to failure, a detailed 3D model is required when it is necessary to examine a limit state associated with the failure mode of the tank. The findings show good agreement between the natural frequencies of convective modes determined from the 2D and 3D CST models. The difference between the frequencies of impulsive modes were found to be equal to or less than 25%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

High-Fidelity Velocity and Concentration Measurements of Turbulent Buoyant Jets

Accurate models of turbulent buoyant flows are essential for the design of the cooling circuit of nuclear reactors and passive safety systems. However, available models fail to fully capture the physics of turbulent mixing when buoyancy becomes predominant with respect to momentum. Therefore, high-fidelity experiments of well-controlled fundamental flows are needed to develop and validate more accurate models. We analyze experiments of positive and negative turbulent buoyant jets, both in uniform and stratified environments, with the aim of understanding the thermal hydraulics of turbulent mixing with variable density and providing high-fidelity data for the development and validation of turbulence models. Non-intrusive, simultaneous particle image velocimetry and laser-induced fluorescence measurements were carried out to acquire instantaneous velocity and concentration fields on a vertical section parallel to the axis of a jet in the self-similar region. The refractive index matching method was applied to measure high-resolution buoyant jets with up to 8.6% density difference. These data are free of the typical errors that characterize optical measurements of buoyancy-driven flows (e.g. natural and mixed convection) where the refractive index of the fluid is inhomogeneous throughout the measurement domain. Turbulent statistics and entrainment of buoyant jets in uniform and stratified environments are presented. These data are compared with non-buoyant jets in a uniform environment, as a reference to investigate the effects of buoyancy and stratification on turbulent mixing. The results will be used for the assessment of current turbulence models and as a basis for the development of a new model that captures turbulent mixing.

laser induced fluorescence↗

FY22 Status Report on the ART-GCR CMVB and CNWG International Collaborations

This work presents the numerical model for the high temperature test reactor (HTTR) loss of forced cooling (LOFC) experiment with the INL codes Griffin, BISON, and RELAP-7, based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Promising results were obtained, with the overall behavior of the reactor successfully captured. Changes in the heat transfer coupling, as compared to the fiscal year (FY)-21 model, enabled drastic improvement of the steady-state solution, both in terms of computational time and global energy discrepancies. The former was reduced by a factor of roughly 60, whereas the latter decreased from 11% to less than 2%. Furthermore, the discrepancy in the steady-state multiplication factor was improved from +2,300 and +2,900 pcm (for the 30 and 9 MW cases) to -700 and +1,200 pcm, respectively, and now falls well within the large measurement uncertainties stemming from graphite impurities. Validation of the Monte Carlo model used to generate cross sections was also performed against available measurements. Though significant, the discrepancies remain acceptable overall in light of the large uncertainty. Numerous improvements are still needed to better compare with the experiment involving the 9 MW case and to instill greater confidence in the model’s ability to accurately predict the 30 MW behavior. Specifically, the power levels predicted by the 9 MW transient simulation following re-criticality remain low, pointing to an underestimation of the passive cooling of the core. A key aspect of future work will be to better understand the flow pattern during the LOFC event, particularly to determine if natural or forced convection is occurring inside the reactor pressure vessel (RPV). More generally, additional validation data would be immensely useful for further enhancing the numerical model and better matching the experiments. In addition, a more sophisticated thermal-hydraulics model that simulates all the channels as a single system model should be considered to take into account the rest of the primary loop. Finally, even if the results are in better agreement with the experiments, sensitivity analysis and uncertainty quantification will be necessary to evaluate the model uncertainty.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Survey of heat transfer to near critical fluids

Survey topics include - heat transfer boundaries of near critical region, free, natural, and forced convection experiments, oscillations, geometric effects, parameters which appear to be significant to heat transfer in critical region, and theories which have been proposed for region.

Hendricks, R. C.↗

Heat Transfer and Fluid Mechanics Institute, Meeting, 26th, Washington State University, Pullman, Wash., June 26-28, 1978, Proceedings

Papers are presented on the heat transfer and fluid mechanics in energy systems including diesel combustion analysis using rapid sampling and the numerical simulation of thermal response of architectural structures. Experiments associated with hypersonic boundary layers, the prediction of nozzle-wall pressure distribution, and the effects of nozzle geometry on an axisymmetric thrust augmentor are outlined. Predictions are made of film thickness on a roll-metering device and the thermal characteristics of multiple-surface submerged condenser modules are described. Numerical studies are made of the interaction between natural and forced convection in a complex geometry and of pressure drop and holdup in two-phase stratified flow. A numerical solution is given for viscous flows at high Reynolds numbers and for compressible boundary layer flows near separation or reattachment. The structure and transport properties in turbulent shear flows are reviewed along with the spectra of Reynolds stresses in a developing pipe flow.

Crowe, C. T.↗

Simulation of the temperature distribution in crystals grown by Czochralski method

Production of perfect crystals, free of residual strain and dislocations and with prescribed dopant concentration, by the Czochralski method is possible only if the complex, interacting phenomena that affect crystal growth in a Cz-puller are fully understood and quantified. Natural and forced convection in the melt, thermocapillary effect and heat transfer in and around the crystal affect its growth rate, the shape of the crystal-melt interface and the temperature gradients in the crystal. The heat transfer problem in the crystal and between the crystal and all other surfaces present in the crystal pulling apparatus are discussed at length. A simulation and computer algorithm are used, based on the following assumptions: (1) only conduction occurs in the crystal (experimentally determined conductivity as a function of temperature is used), (2) melt temperature and the melt-crystal heat transfer coefficient are available (either as constant values or functions of radial position), (3) pseudo-steady state is achieved with respect to temperature gradients, (4) crystal radius is fixed, and (5) both direct and reflected radiation exchange occurs among all surfaces at various temperatures in the crystal puller enclosure.

Dudokovic, M. P.↗

Theoretical modelling of the fine structures in sunspots

Until a decade ago most solar physicists thought of a sunspot as the upper end of a giant flux tube floating vertically. The existence of umbral dots and penumbral grains has been known for several decades. On the basis of available observations, they seem to be regions of photospheric intensity with upflowing gas motion and magnetic fields much weaker than in the surrounding sunspot surface. It has also been suggested that the differences in the appearances of umbral dots and granular cells are caused by the highly nonlinear nature of the convection problem in the presense of strong magnetic fields. The main ideas are presented here without any equations. It can be shown that a pocket of field free gas surrounded by a vertical magnetic field in the presence of gravity takes up the shape of a tapering column ending at a vertex at the top. Some convection is expected to take place in the trapped field free gas, whereas the magnetic field around it makes those regions stable against convection. Eventually the apex of the tapering column reaches the photospheric surface where the bulging of the magnetic field makes the field no longer able to close on the field free gas and trap it underneath.

Choudhuri, Arnab Rai↗

Survey Of Cooling Techniques

Methods of maintaining electronics at best operating temperatures reviewed. Paper surveys cooling techniques for electronic components and sensors that must be kept at temperatures ranging from thousandths of kelvin to room temperature. Discusses characteristics of various coolers in terms of ground and space applications, cooling capacity, reliability, and temperature range. Also discusses briefly, cooling of conventional electronic circuitry operating few degrees above environmental temperature by natural or forced convection of air or liquid. At greater length, discusses passive and active refrigeration of "cold" electronic circuitry operating below environmental temperature.

Chan, Chung K.↗

Flow Effects on the Flammability Diagrams of Solid Fuels

A research program is currently underway with the final objective of developing a fundamental understanding of the controlling mechanisms underlying the flammability diagrams of solid combustible materials and their derived fire properties. Given that there is a high possibility of an accidental fire occurring in a space-based facility, understanding the fire properties of materials that will be used in such facilities is of critical importance. With this purpose, the flammability diagrams of the materials, as those produced by the Lateral Ignition and Flame Spread Test (LIFT) apparatus and by a new forced flow device, the Forced Flow Ignition and Flame Spread Test (FIST) apparatus, will be obtained. The specific objective of the program is to apply the new flammability apparatus, which will more accurately reflect the potential ambient conditions of space-based environments, to the characterization of the materials for space applications. This paper presents a parametric study of oxidizer flow effects on the ignition curve of the flammability diagrams of PMMA. The dependence of the ignition delay time on the external radiant flux and either the sample width (LIFT) or the flow velocity (FIST) has been studied. Although preliminary, the results indicate that natural and forced convection flow changes, affect the characteristics of the ignition curves of the flammability diagrams. The major effect on the ignition time appears to be due to convective transfer variations at the fuel surface. At high radiant fluxes or high flow velocities, however, it appears that gas phase processes become increasingly important, affecting the overall ignition delay time. A numerical analysis of the solid fuel heating and pyrolysis has also been developed. The theoretical predictions approximate the experiments well for conditions in which the gas phase induction time is negligible.

Cordova, J. L.↗