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.
Contains data from the model validation in the 1D Heat Loss Models to Predict the Aquifer Temperature Profile during Hot/Cold Water Injection Project. The data include two COMSOL models (2D axisymmetric benchmark model and 2D Vinsome model), one python code (1D Vinsome based FEM numerical simulation), one matlab main code (1D Newton analytical solution and all results comparison visualization), and output files generated from the above models.
A simple eddy-viscosity model is shown to make it possible to calculate numerically the mean properties of a turbulent wake. Although the structure of the Reynolds stress terms is not resolved, the results obtained are adequate for predicting velocity profiles and displacement thicknesses.
Development of a simple model of the accretion process of the moon to explain its initial temperature distribution. The model assumes that the moon was formed from the accumulation of the solid particles or gases in the isolated, closed cloud. Two equations are derived to calculate the accretion rate and surface temperature of the accreting moon. Numerical calculations are made for a wide range of the particle concentrations and particle velocities in the cloud. A limited set of parameters gives the initial temperature profiles as required by geochemical and geophysical data. These models of the protomoon cloud indicate that the lunar outershell, about 400 km thick, was partially or completely molten just after the accretion of the moon and that the moon should have been formed in a period shorter than 1000 yr. If the moon formed at a position nearer to the earth than its present one, the moon might have been formed in a period of less than one year.
The paper discusses the results of an experimental and numerical investigation of tangential slot injection film cooling with zero pressure gradients in subsonic boundary layers at freestream Mach numbers of 0.4, 0.6, and 0.8. The results are compared with the predictions obtained from a finite-difference boundary-layer program developed by NASA for slot injection into turbulent supersonic boundary layers. The two sets of results are found to compare favorably. The numerical results point to the existence of a unique relation between isothermal effectiveness and adiabatic effectiveness, thereby confirming the existence of similarity conditions for temperature and velocity profiles.
A discrete time model for simulation of the dynamics of samarium cobalt-type permanent magnet brushless dc machines is presented. The simulation model includes modeling of the interaction between these machines and their attached power conditioners. These are transistorized conditioner units. This model is part of an overall discrete-time analysis of the dynamic performance of electromechanical actuators, which was conducted as part of prototype development of such actuators studied and built for NASA-Johnson Space Center as a prospective alternative to hydraulic actuators presently used in shuttle orbiter applications. The resulting numerical simulations of the various machine and power conditioner current and voltage waveforms gave excellent correlation to the actual waveforms collected from actual hardware experimental testing. These results, numerical and experimental, are presented here for machine motoring, regeneration and dynamic braking modes. Application of the resulting model to the determination of machine current and torque profiles during closed-loop actuator operation were also analyzed and the results are given here. These results are given in light of an overall view of the actuator system components. The applicability of this method of analysis to design optimization and trouble-shooting in such prototype development is also discussed in light of the results at hand.
The program for Regional Observing and Forecasting Services (PROFS) has been using wind profile data in experimental forecast applications for over two years, mostly in the form of real-time color displays on the PROFS forecast workstation. The most ambitious test of the workstation to date, the 1985 PROFS Real-Time Experiment (RT-85), ran from 15 May to 23 August, 1985. The use of wind profiler products during this and previous experiments is described. Data from the experimental profiler network in Colorado and from the profiler in Oklahoma are in the form of hourly averages. Arriving data frequently contain errors whose origins range from interference by aircraft in the beams to highway truck traffic. Most of the irregularities are apparent through visual inspection of profiler wind observations plotted on a time-height cross section, but this method of quality control is inadequate if the intended uses of the data involve numerical calculations.
The charge of cosmic grains could play an important role in many astrophysical phenomena. It probably has an influence on the coagulation of grains and more generally on grain-grain collisions, and on interaction between charged particles and grains which could lead to the formation of large grains or large molecules. The electrostatic charge of grains depends mainly on the nature of constitutive material of the grain and on the physical properties of its environment: it results from a delicate balance between the plasma particle collection and the photoelectron emission, both of them depending on each other. The charge of the grain is obtained in two steps: (1) using the numerical model the characteristics of the environment of the grain are computed; (2) the charge of a grain which is embedded in this environment is determined. The profile of the equilibrium charge of some typical grains through different types of interstellar clouds is obtained as a function of the depth of the cloud. It is shown that the grain charge can reach high values not only in hot diffuse clouds, but also in clouds with higher densities. The results are very sensitive to the mean UV interstellar radiation field. Three parameters appear to be essential but with different levels of sensitivity of the charge: the gas density, the temperature, and the total thickness of the cloud.
Potential uses of specific extremely High Frequency (EHF) and Sub-Millimeter-Wave (SMMW) channels at 90, 166, 183, 220, 325, 340, and 410 GHz for passive spaceborne remote sensing of the troposphere and lower stratosphere are investigated using an iterative numerical radiative transfer model. Collectively, these channels offer potential for high spatial resolution imaging using diffraction-limited apertures of practical size, along with the ability to profile water vapor, map precipitation beneath optically opaque cloud cover, and to measure nonprecipitating cloud (e.g., cirrus) parameters. A widely-spaced set of EHF and SMMW channels can yield observable degrees of freedom related to clouds and precipitation not available by exclusively using the more thoroughly studied microwave channels below 183 GHz. A new passive airborne imaging instrument for tropospheric meteorological sensing is described.
The pressure-driven flow of a rarefied monatomic gas through a two-dimensional slit is simulated using the direct simulation Monte Carlo technique. Of particular interest is the change in flow field structure as pressure ratio and Knudsen number are varied. Comparisons are made to quantify the limits of validity of free-molecular theory and approximate, nearly free-molecular iterative methods. Also addressed is the sensitivity of the numerical solutions to grid structure and boundary conditions. The free-molecular theory is found to predict quantitative flow field properties (e.g., centerline velocities or downstream flux profiles) reasonably well for large finite Knudsen number with the error dependent on the pressure ratio. The nearly free-molecular corrections are shown to have limited range of applicability. A previously derived parameter is found to correlate total mass flux well as a function of pressure ratio and Knudsen number over a large portion of the transitional regime.
A model of the interaction of the noise from a spreading subsonic jet with a 4 panel assembly is studied numerically in two dimensions. The effect of forward motion of the jet is accounted for by considering a uniform flow field superimposed on a mean jet exit profile. The jet is initially excited by a pulse-like source inserted into the flow field. The pulse triggers instabilities associated with the inviscid instability of the jet shear layer. These instabilities generate sound which in turn serves to excite the panels. We compare the sound from the jet, the responses of the panels and the resulting acoustic radiation for the static jet and the jet in forward motion. The far field acoustic radiation, the panel response and sound radiated from the panels are all computed and compared to computations of a static jet. The results demonstrate that for a jet in forward motion there is a reduction in sound in downstream directions and an increase in sound in upstream directions in agreement with experiments. Furthermore, the panel response and radiation for a jet in forward motion exhibits a downstream attenuation as compared with the static case.
The goal of a coordinated effort between groups at GWU and NASA GSFC is the development of a low-cost, global, surface instrument network that continuously monitors three key carbon cycle gases in the atmospheric column: carbon dioxide (CO2), methane (CH4), carbon monoxide (CO), as well as oxygen (O2) for atmospheric pressure profiles. The network will implement a low-cost, miniaturized, laser heterodyne radiometer (mini-LHR) that has recently been developed at NASA Goddard Space Flight Center. This mini-LHR is designed to operate in tandem with the passive aerosol sensor currently used in AERONET (a well established network of more than 450 ground aerosol monitoring instruments worldwide), and could be rapidly deployed into this established global network. Laser heterodyne radiometry is a well-established technique for detecting weak signals that was adapted from radio receiver technology. Here, a weak light signal, that has undergone absorption by atmospheric components, is mixed with light from a distributed feedback (DFB) telecommunications laser on a single-mode optical fiber. The RF component of the signal is detected on a fast photoreceiver. Scanning the laser through an absorption feature in the infrared, results in a scanned heterodyne signal io the RF. Deconvolution of this signal through the retrieval algorithm allows for the extraction of altitude contributions to the column signal. The retrieval algorithm is based on a spectral simulation program, SpecSyn, developed at GWU for high-resolution infrared spectroscopies. Variations io pressure, temperature, composition, and refractive index through the atmosphere; that are all functions of latitude, longitude, time of day, altitude, etc.; are modeled using algorithms developed in the MODTRAN program developed in part by the US Air Force Research Laboratory. In these calculations the atmosphere is modeled as a series of spherically symmetric shells with boundaries specified at defined altitudes. Temperature, pressure, and species mixing ratios are defined at these boundaries. Between the boundaries, temperature is assumed to vary linearly with altitude while pressure (and thus gas density) vary exponentially. The observed spectrum at the LHR instrument will be the integration of the contributions along this light path. For any absorption measurement the signal at a particular spectral frequency is a linear combination of spectral line contributions from several species. For each species that might absorb in a spectral region, we have pre-calculated its contribution as a function of temperature and pressure. The integrated path absorption spectrum can then by calculated using the initial sun angle (from location, date, and time) and assumptions about pressure and temperature profiles from an atmospheric model. The modeled spectrum is iterated to match the experimental observation using standard multilinear regression techniques. In addition to the layer concentrations, the numerical technique also provides uncertainty estimates for these quantities as well as dependencies on assumptions inherent in the atmospheric models.
Recap: NASA-Specific Objectives for ICE-POP: Provide real-time observational and NWP data in support of ICE-POP, participate in significant international science effort; GPM (Global Precipitation Measurement) Ground Validation and NASA Weather Program -Direct/physical validation of active/passive satellite-based snowfall retrieval algorithms over coastline and mountains; melting layer interaction with terrain -Physics of snow, coupling to snow water equivalent rate and satellite remote sensor retrieval algorithm assumptions - -Size distributions, types/habit, water equivalent, profiles -NU-WRF (NASA-Unified Weather Research and Forecasting) Model plus Observational analyses: Movement toward “level IV products” leverage intensive and multi-faceted NWP (Numerical Weather Prediction) component -Model precipitation processes (liquid, mixed phase and frozen); Build model testing database for further active/passive remote sensing algorithm development (e.g., satellite data simulators) -"Integrated" validation of products in operational context.
This paper describes the one-node P{sub N}-Source Expansion Nodal Method (SENM) axial solver for the 2D/1D method recently implemented in MPACT to support Boiling Water Reactors (BWR) analysis. Since the BWR has a more complicated design and strong burnable absorber, the existing PN-Nodal Expansion Method (NEM) axial solver in MPACT may not be sufficient to accurately represent the intranodal flux and source profiles for BWRs. The one-node P{sub N}-SENM has been implemented in this work to reduce the axial spatial discretization error for BWRs. From numerical results, we verify that the P{sub N}-SENM can improve the accuracy of the pin power prediction, and confirm that P{sub N}-SENM can use more than a 1.5 times larger axial mesh size than P{sub N}-NEM to have similar accuracy for the BWR GE14 3D assembly problem. (authors)
Solid-phase diffusion in active materials of lithium-ion batteries significantly affects charging and safety-related behavior of lithium-ion batteries. Therefore, it is essential to develop an efficient and robust numerical algorithm for solving solid-phase diffusion equations in physics-based battery models. In this work, we discuss the origins of numerical instabilities that can occur when solving the solid-phase diffusion equations using iterative methods. Then, in order to resolve such issues, we propose a simple numerical treatment to the surface flux term of discretized solid-phase diffusion equations. To demonstrate its numerical robustness, the proposed method is implemented into a pseudo two-dimensional (P2D) physics-based battery model and simulations are conducted at wide ranges of operating conditions. Even with extremely poor initial guesses for the Li+ concentrations of the active materials, computations using the proposed method do not diverge and the their computational speeds are comparable to those with conventional initial guesses. Comprehensive tests of the proposed method are also performed with a dynamic current profile based on US06 driving profile and a multi-stage charging profile with very high initial C-rate (12C).
Comparisons of radial dose calculation are performed, as well as the introduction of important physics to improve the calculation techniques. Also, the consequences to device performance are explored via numerical simulations.
Microorganisms efficiently coordinate phenotype expressions through a decision-making process known as quorum sensing (QS). We investigated QS amongst distinct, spatially distributed microbial aggregates under various flow conditions using a process-driven numerical model. Model simulations assess the conditions suitable for QS induction and quantify the importance of advective transport of signaling molecules. In addition, advection dilutes signaling molecules so that faster flow conditions require higher microbial densities, faster signal production rates, or higher sensitivities to signaling molecules to induce QS. However, autoinduction of signal production can substantially increase the transport distance of signaling molecules in both upstream and downstream directions. We present empirical approximations to the solutions of the advection–diffusion–reaction equation that describe the concentration profiles of signaling molecules for a wide range of flow and reaction rates. These empirical relationships, which predict the distribution of dissolved solutes along pore channels, allow to quantitatively estimate the effective communication distances amongst multiple microbial aggregates without further numerical simulations.
We have imaged the planetary nebula (PN) NGC 246 in the near-ultraviolet wavelengths [Ne v] 342.6 nm, the Bowen fluorescence line of 0 111 at 344.4 nm, and a nearby line-free region centered on 338.6 nm, as well as H(alpha), [O III] 500.7 nm, and [S II] 673.0 and 671.5 nm. Imaging in the 344.4 nm line is necessary to deconvolve contamination of the [Ne v] images by O III 342.9 nm. The emission from the shell and inner parts of the nebula is detected in [Ne v]. The radial profiles of the [Ne v] brightness decrease with radius from the exciting star, indicating that the bulk of the emission from this ion is due to the hard UV stellar radiation field, with a (probably) marginal contribution from collisional ionization in a shock between the PN shell and the interstellar medium (ISM). In contrast, the radial profiles of the emission in H(alpha), [0 III] 500.7 nm, and [S II] are flatter and peak at the location of the shell. The emission of [S II] probably traces the interaction of the PN with the ambient ISM. We also present two-dimensional numerical simulations for this PN-ISM interaction. The simulations consider the stellar motion with respect to the ambient ISM, with a velocity of 85 km/s , and include the time evolution of the wind parameters and UV radiation field from the progenitor star.
Multiscale experiments in heterogeneous materials and the knowledge of their physics under shock compression are limited. Here, this study examines the multiscale shock response of particulate composites comprised of soda-lime glass particles in a PMMA matrix using full-field high speed digital image correlation (DIC) for the first time. Normal plate impact experiments, and complementary numerical simulations, are conducted at stresses ranging from 1.1 - 3.1 GPa to elucidate the mesoscale mechanisms responsible for the distinct shock structure observed in particulate composites. The particle velocity from the macroscopic measurement at continuum scale shows a relatively smooth velocity profile, with shock thickness decreasing with an increase in shock stress, and the composite exhibits strain rate scaling as the second power of the shock stress. In contrast, the mesoscopic response was highly heterogeneous, which led to a rough shock front and the formation of a train of weak shocks traveling at different velocities. Additionally, the normal shock was seen to diffuse the momentum in the transverse direction, affecting the shock rise and the rounding-off observed at the continuum scale measurements. The numerical simulations indicate that the reflections at the interfaces, wave scattering, and interference of these reflected waves are the primary mechanisms for the observed rough shock fronts.