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.
The forced-convection heat-transfer characteristics of sodium hydroxide were experimentally investigated. The heat-transfer data for heating fall slightly above the McAdams correlation line, and the heat-transfer data for cooling are fairly well represented by the McAdams correlation line.
A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.
GFSSP has been extended to model conjugate heat transfer Fluid Solid Network Elements include: a) Fluid nodes and Flow Branches; b) Solid Nodes and Ambient Nodes; c) Conductors connecting Fluid-Solid, Solid-Solid and Solid-Ambient Nodes. Heat Conduction Equations are solved simultaneously with Fluid Conservation Equations for Mass, Momentum, Energy and Equation of State. The extended code was verified by comparing with analytical solution for simple conduction-convection problem The code was applied to model: a) Pressurization of Cryogenic Tank; b) Freezing and Thawing of Metal; c) Chilldown of Cryogenic Transfer Line; d) Boil-off from Cryogenic Tank.
A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.
Current liquid oxygen feed systems waste propellant and use hardware, unnecessary during flight, to condition the propellant at the engine turbopumps prior to launch. Simplified liquid oxygen propellant conditioning concepts are being sought for future launch vehicles. During a joint program, four alternative propellant conditioning options were studied: (1) passive recirculation; (2) low bleed through the engine; (3) recirculation lines; and (4) helium bubbling. The test configuration for this program was based on a vehicle design which used a main recirculation loop that was insulated on the downcomer and uninsulated on the upcomer. This produces a natural convection recirculation flow. The test article for this program simulated a feedline which ran from the main recirculation loop to the turbopump. The objective was to measure the temperature profile of this test article. Several parameters were varied from the baseline case to determine their effects on the temperature profile. These parameters included: flow configuration, feedline slope, heat flux, main recirculation loop velocity, pressure, bleed rate, helium bubbling, and recirculation lines. The heat flux, bleed rate, and recirculation configurations produced the greatest changes from the baseline temperature profile. However, the temperatures in the feedline remained subcooled. Any of the options studied could be used in future vehicles.
A determination of nominal flow phenomena in liquid metal fast reactor (LMFR) fuel assemblies is critical toward generation-IV reactor development. Axially positioned spacer grids are used to maintain the geometry of hexagonal rod bundles and simultaneously introduce perturbations in the flow. Three-dimensional (3D) printed asymmetric honeycomb spacer grids were installed in a prototypical 127-pin LMFR fuel assembly model to study complex fluid dynamics interactions induced by the spacer grid and rods. To characterize flow dynamics in this intricate geometry, time-resolved particle image velocimetry (TR-PIV) using the matched-index-of-refraction method was employed to obtain non-intrusive velocity measurements for three axial planes (one near-wall and two interior planes) at a Reynolds number of 6000. The statistical TR-PIV results compared sub-channel-dependent normalized time-averaged velocity, velocity fluctuations, Reynolds stress, vorticity, and turbulence kinetic energy distributions. TR-PIV line profiles characterized downstream spacer grid flow dynamics. Two-point spatial and spatial–temporal cross-correlation fields revealed local coherent structures and quantified convection velocities of traveling vortices. Spatial–temporal decomposition using dynamic mode decomposition (DMD) applied to the near-wall vorticity fields extracted turbulent structures and flow instabilities in the wake region of the spacer grid, along with their decay and frequency rates. Reduced-order velocity fields from DMD reconstructions identified the most energy-containing coherent structures persistent in the near-wall region. This research provides experimental data sets and analyses of flow behavior in rod bundles with hexagonal spacer grids. In conclusion, the results are critical toward LMFR design and geometry optimization, crucial for the validation of computational fluid dynamics and reduced-order flow models.
Frequency measurements for the Delta V = 2 transitions of CO in the integrated light spectrum of the sun are presented. The nature and magnitude of systematic errors which typically arise in absolute velocity measurements of integrated sunlight are explored in some detail, and measurements believed accurate at the level of about 5 m/s or less are presented. It is found that the integrated light velocity varies by about 3 m/s or less over a one-day period. Over the long term, the data indicate an increasing blue-shift in these weak infrared lines amounting to 30 m/s from 1983 to 1985. The sense of the drift is consistent with a lessening in the magnetic inhibition of granular convection at solar minimum. Such an effect has implications for the spectroscopic detectability of planetary-mass companions to solar-type stars.
This study extends the reduced Navier-Stokes (RNS) global pressure relaxation procedure developed by Rubin and co-workers for external flow to internal flow applications. The streamwise pressure gradient is split into a backward-differenced or initial value component, as in boundary layer marching, and a forward-differenced or boundary value component that represents the elliptic downstream effects. The streamwise convection terms are upwind-differenced and all other streamwise derivatives are backward-differenced. A standard boundary layer marching technique imbedded in a conventional line relaxation technique is obtained. For compressible flow the pressure iteration determines the interior flow interaction as well as the inlet mass flux that is consistent with the outflow pressure boundary condition. Results have been computed for incompressible flow in both rectangular and curved channels, and for subsonic compressible flow in the simulation of an aerofoil in a wind tunnel. Converged solutions were obtained over a range of Reynolds numbers generating small to moderately large separation bubbles.
A full 3D Navier-Stokes numerical investigation has been conducted of the shock-wave/boundary-layer flow interactions caused by four diamond-shaped struts, of varying thickness, in an annular duct with Mach 3 core flow and turbulent boundary-layers on both walls. Secondary flows caused by weak-to-strong interactions are examined in the vicinity of a strut which is bounded by curved endwalls. The duct endwall boundary-layer separated for the strongest interaction. The struts studied had maximum thickness-to-chord ratios of 0.125, 0.188, 0.250, and 0.500. The duct gap height is 0.7 strut chords, the duct inner-to-outer wall radius ratio is 0.7, and the Reynolds number is 3 x 10 exp 5 based on the strut chord length which was held constant for all interactions considered. The effects of strut thickness on the secondary flows are discussed, including: trajectories for the leading and trailing edge horseshoe vortices, strut/endwall corner vortices, and boundary-layer separation. The line of coalescence discussed in the literature, previously ascribed to boundary-layer separation, is shown to be caused by the leading edge horseshoe vortex convecting along the shock front.
A high resolution global magnetohydrodynamic simulation model is used to investigate magnetospheric dynamics during intervals with southward interplanetary magnetic field (IMF). When the southward IMF reaches the dayside magnetopause reconnection begins and magnetic flux is convected into the tail lobes. After about 35 m, reconnection begins within the plasma sheet near midnight at x = -14 RE. Later the x-line moves towards the magnetopause. The reconnection occurs just tailward of the region where the tail attaches onto the dipole dominated inner magnetosphere. Later when all the plasma sheet field lines have reconnected a plasmoid moves down the tail. The region of the ionosphere where the energy flux from the magnetosphere is greatest is calculated. The energy flux is confined to a region which approximates the auroral oval.
Concerning our knowledge of convective clouds and their important processes, GATE marks a watershed between "seeing through a glass darkly" and "then face to face" in some of the most vital aspects of convective systems. Two main areas will be discussed here, namely the interactions of observations and modeling in documentation of slow and fast-moving squall lines, for which the publications began in the 70's, and the use of the GATE ships' radar data set as the basis for demonstrating that tropical rainfall could be measured from space, resulting in the Tropical Rainfall Measuring Mission (TRMM). The TRMM satellite was launched in 1997, with the enabling analyses reaching publication in the 1980's.
Cold pools are increasingly being recognized as important players in the evolution of both shallow and deep convection; hence, the incorporation of cold pool processes into a number of recently developed convective parameterizations. Unfortunately, observations serving to inform cold pool parameterization development are limited to select field programs and limited radar domains. However, a number of recent studies have noted that cold pools are often associated with arcs-lines of shallow clouds traversing 10 100 km in visible satellite imagery. Boundary layer thermodynamic perturbations are plausible at such scales, coincident with such mesoscale features. Atmospheric signatures of features at these spatial scales are potentially observable from satellites. In this presentation, we discuss recent work that uses multi-sensor, high-resolution satellite products for observing mesoscale wind vector fluctuations and boundary layer temperature depressions attributed to cold pools produced by antecedent convection. The relationship to subsequent convection as well as convective system longevity is discussed. As improvements in satellite technology occur and efforts to reduce noise in high-resolution orbital products progress, satellite pixel level (10 km) thermodynamic and dynamic (e.g. mesoscale convergence) parameters can increasingly serve as useful benchmarks for constraining convective parameterization development, including for regimes where organized convection contributes substantially to the cloud and rainfall climatology.
We have calculated a grid of more than 700 model atmospheres valid for a wide range of parameters encompassing the coolest known M dwarfs, M subdwarfs, and brown dwarf candidates: 1500 less than or equal to T(sub eff) less than or equal to 4000 K, 3.5 less than or equal to log g less than or equal to 5.5, and -4.0 less than or equal to (M/H) less than or equal to +0.5. Our equation of state includes 105 molecules and up to 27 ionization stages of 39 elements. In the calculations of the base grid of model atmospheres presented here, we include over 300 molecular bands of four molecules (TiO, VO, CaH, FeH) in the JOLA approximation, the water opacity of Ludwig (1971), collision-induced opacities, b-f and f-f atomic processes, as well as about 2 million spectral lines selected from a list with more than 42 million atomic and 24 million molecular (H2, CH, NH, OH, MgH, SiH, C2, CN, CO, SiO) lines. High-resolution synthetic spectra are obtained using an opacity sampling method. The model atmospheres and spectra are calculated with the generalized stellar atmosphere code PHOENIX, assuming LTE, plane-parallel geometry, energy (radiative plus convective) conservation, and hydrostatic equilibrium. The model spectra give close agreement with observations of M dwarfs across a wide spectral range from the blue to the near-IR, with one notable exception: the fit to the water bands. We discuss several practical applications of our model grid, e.g., broadband colors derived from the synthetic spectra. In light of current efforts to identify genuine brown dwarfs, we also show how low-resolution spectra of cool dwarfs vary with surface gravity, and how the high-regulation line profile of the Li I resonance doublet depends on the Li abundance.
This study investigates the evolution of the magnetotail's magnetic field with the aid of a self-consistent two-dimensional model. In this model the plasma mantle continuously supplies particles to the magnetotail, the ion current periodically updates the magnetic field using the Biot-Savart law. The simulated magnetotail evolves into a quasi-steady state, characterized by the periodic motion of the model's near-Earth X-line. This variability results from the nonadiabatic acceleration of ions in the current sheet and their rapid loss from the tail. The characteristic time scale of variability in the magnetotail is on the order of 4 - 5 minutes. We also investigate how the magnetotail's topology responds to increased convection electric fields, and show examples of observations of variability in the magnetotail.
Spacecraft measurements of the plasma populations and magnetic fields near Jupiter and Saturn have revealed that large magnetospheres surround both planets. Magnetic field measurements have indicated closed field line topologies in the dayside magnetospheres of both planets while plasma instruments have shown these regions to be populated by both hot and cold plasma components convected azimuthally in the sense of planetary rotation. By using published data from the Voyager Plasma Science (PLS), Low Energy Charged Particle (LECP), and Magnetometer (MAG) instruments, it is possible to investigate the validity of the time stationary MHD momentum equation in the middle magnetospheres of Jupiter and Saturn. At Saturn, the hot plasma population is negligible in the dynamic sense and the centrifugal force of the cold rotating plasma appears to balance the Lorentz force. At Jupiter, the centrifugal force balances about 25 percent of the Lorentz force. The remaining inward Lorentz force is balanced by pessure gradients in the hot, high-beta plasma of the Jovian magnetodisk.
The first simultaneous study of dc and ac electric and magnetic fields, E x B velocity, plasma flows, ratio of plasma to magnetic field pressure, total energy density, energetic particles, and ion composition from the ISEE satellites and ground and interplanetary magnetic fields has been made to determine (1) the relationship of the previously observed electric fields at the plasma sheet boundary and at the neutral sheet to plasma parameters, and (2) whether the phenomena occurring during quiet and active times were consistent with the formation of a near-earth neutral line during substorms or with the boundary layer model. Five observations made during the study of two substorms were seen to be in agreement with the neutral-line model. The observations are consistent with the satellite being located at varying distances from the neutral line and diffusion region where reconnection and plasma acceleration were occurring. Although the z component (into or out of the ecliptic plane) of E x B convection was generally toward the neutral sheet, there were examples when it was consistent with the inferred motion of the plasma sheet past the satellite. A synthesis of previous reports on large electric fields at the plasma sheet boundary and variable fields at the neutral sheet including the associated plasma flows is also described.
A numerical model was developed to examine the flow behavior of He II under forced flow conditions. The model basically solves the governing equations in the two-fluid model of He II, and can be used to assess the pump requirement for transferring He II in spaceborne systems. Predictions of liquid temperature and pressure drop agree well with empirical data obtained from small scale tests. At high transfer rates, the model predicts that the forced convection dominates the Gorter-Mellink conduction in heat transfer, and the classical turbulent flow correlation is sufficient for calculating the pressure drop across the transfer line.
Magnetic field reconnection is a fundamental process that occurs in the magnetotail during geomagnetic substorms. Some 2D reconnection models predict the formation of a plasmoid, or closed loop of magnetic field lines, in the noon-midnight meridional plane at those times. When the 3D magnetotail magnetic field is considered, it becomes clear that reconnection produces a flux rope with an axis transverse to the earth-sun line. Three signatures mark both 2D plasmoids and 3D flux ropes: (1) a bipolar magnetic field signature, (2) tailward flow of a hot plasma, and (3) convecting isotropic energetic particle distributions. Plasmoids and flux ropes may be distinguished by (4) the axial magnetic field that only flux ropes possess. All four signatures have been identified in near-earth, middle, and distant magnetotail observations, but their interpretation is disputed. Thus, the existence of magnetotail flux ropes remains a controversial subject.