SEARCH · Engineering Papers
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.
Core hydrogen-burning models are discussed together with the associated physical instabilities. A mathematical relation for dynamical instability is presented. Conditions leading to convective instability are analyzed, giving attention to the nature of the convectively unstable, chemically inhomogeneous, zone that is left behind by the retreating convective core. Questions of pulsational instability are also investigated.
Advances in computing capabilities and scientific development have come together to evolve general circulation models into multi-scale Earth system modeling tools. The Goddard Earth Observing System (GEOS) model is one such example of this evolution. The GEOS model is driven by the finite-volume cubed-sphere (FV3) non-hydrostatic dynamical core. Surrounding FV3 is a scale-aware physics package and data assimilation capability permitting multi-scale application of GEOS for sub-seasonal to seasonal climate prediction, medium range weather prediction, and global mesoscale modeling at convection allowing resolutions. GEOS also includes a comprehensive chemistry package representing a range of capabilities from basic chemistry and interactve aerosols and gaseous species, to carbon emissions and uptake, and complex ozone photochemistry. In this study, we apply the GEOS model to study the fidelity of these processes with increasing horizontal resolution and scale-aware processes in GEOS on extreme events. The GEOS model is run for 40-days beginning in August 2016 at three uniform global resolutions of 13-km (c768), 6-km (c1536) and 3-km (c3072) with 72 vertical levels up to 0.01mb. The model physics use the Grell-Freitas scale-aware convection scheme to dynamically reduce the role of parameterized deep convection as resolved scale processes in the model take over at higher resolutions. We will include high-resolution global emissions and fluxes of aerosols and carbon downscaled from recent satellite observations. We will compare these simulations with reanalyses and observations, focusing on rainfall, clouds and radiative forcing at hourly to monthly timescales. We will closely examine the probability distribution of precipitation intensities and radiative properties of clouds on a daily time scale. In addition, we will focus on extreme events, in particular the diurnal cycle of convection over the US and the frequency and physical nature of organized convection and heavy rain events across the globe.
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.
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.
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.
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.
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.
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.
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.
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.
Tiros i photography of cellular cloud pattern related to conventional meteorological analyses, physical nature of cellular convection
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.
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.
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.
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.
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.
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.
The burning of liquid hydrocarbon fuels accounts for a significant portion of global energy production. With predicted future increases in demand and limited reserves of hydrocarbon fuel, it is important to maximize the efficiency of all processes that involve conversion of fuel. With the exception of unwanted fires, most applications involve introduction of liquid fuels into an oxidizing environment in the form of sprays which are comprised of groups of individual droplets. Therefore, tremendous benefits can result from a better understanding of spray combustion processes. Yet, theoretical developments and experimental measurements of spray combustion remains a daunting task due to the complex coupling of a turbulent, two-phase flow with phase change and chemical reactions. However, it is recognized that individual droplet behavior (including ignition, evaporation and combustion) is a necessary component for laying the foundation for a better understanding of spray processes. Droplet combustion is also an ideal problem for gaining a better understanding of non-premixed flames. Under the idealized situation producing spherically-symmetric flames (produced under conditions of reduced natural and forced convection), it represents the simplest geometry in which to formulate and solve the governing equations of mass, species and heat transfer for a chemically reacting two phase flow with phase change. The importance of this topic has promoted extensive theoretical investigations for more than 40 years.