Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy simulation”

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 235 records · Page 13

Particle acceleration by propagating interplanetary shocks

The process of charged particle acceleration in interplanetary shocks has been simulated for typical parameters. Since space probe observations of charged particle fluxes are the principle means of inferring the action of an acceleration mechanism, the simulation was designed to follow particles backwards in time from a given observing point, usually 1 AU. This allowed a full diagnosis of which observed particles had interacted with an oncoming shock, where the interacting occurred, and by how much the energy had been changed by the shock interaction. Simple assumptions about the preshock energy spectrum allow the construction of a full temporal profile of the expected intensity, anisotropy, and energy spectrum. The simulations are apparently capable of reproducing the main features of the less than 10 MeV/nuc ion flux enhancements observed at interplanetary shocks using only a laminar interplanetary magnetic field, Archimidean spiral, and a spherical oblique shock with typical speed, strength, and shock normal-magnetic field angle.

Armstrong, T. P.↗

Protecting Venus from Asteroids, Comets, and Meteors

It is well accepted that the dense, thick atmosphere of Venus prevents most small cosmic bodies from reaching the surface and forming craters. We have examined this atmospheric intervention in detail, incorporating the lessons learned from the extensive modeling of impactor deceleration and flattening motivated by the SL-9 impacts with Jupiter. We employ a "pancake" model, which best matches detailed code simulations of atmospheric energy deposition, and Schmidt-Holsapple crater scaling modified for complex (flattened) craters. We adopt the distributions of Venus-crossing asteroids and comets determined by E.M. Shoemaker and co-workers, as well as generalizations of these distributions. Our nominal simulation of the venusian crater record is shown below, calibrated to the total number of venusian craters (940). As nearly all craters on Venus are well-preserved and relatively uniformly distributed, such simulations constrain the age of the surface. The fit is reasonable, with a nominal crater retention age of approx. 700 Ma. The fit at the large-crater end is improved if the number of large asteroids is increased, which Shoemaker argues is in fact more representative of the long-term (over several 100 Ma) average, and if Halley-family comets are included. The ages we obtain under a variety of modeling choices that produce good fits (including using Shoemaker's preferred crater scaling) are approx. 700-900 Ma, substantially greater than the most widely cited age estimate in the literature (-300 Ma). The key difference is that we find very large depletions in the production of 20-30-km craters (see figure) compared with previous estimates, the size range at which atmospheric effects are often calibrated or assumed nearly negligible. As venusian global resurfacing recedes deeper into history, the likelihood that Venus is resting between bouts of activity diminishes. Venus, like Mars, may instead be dying or dead.

McKinnon, William B.↗

Numerical simulation of large-scale ocean-atmosphere coupling and the ocean's role in climate

The problem of reducing model generated sigma coordinate data to pressure levels is considered. A mass consistent scheme for performing budget analyses is proposed, wherein variables interpolated to a given pressure level are weighted according to the mass between a nominal pressure level above and either a nominal pressure level below or the Earth's surface, whichever is closer. The method is applied to the atmospheric energy cycle as simulated by the OSU two level atmospheric general circulation model. The results are more realistic than sigma coordinate analyses with respect to eddy decomposition, and are in agreement with the sigma coordinate evaluation of the numerical energy sink. Comparison with less sophisticated budget schemes indicates superiority locally, but not globally.

Gates, W. L.↗

Battery model for electrical power system energy balance

A model to simulate nickel-cadmium battery performance and response in a spacecraft electrical power system energy balance calculation was developed. The voltage of the battery is given as a function of temperature, operating depth-of-charge (DOD), and battery state-of-charge. Also accounted for is charge inefficiency. A battery is modeled by analysis of the results of a multiparameter battery cycling test at various temperatures and DOD's.

Hafen, D. P.↗

High Energy Electron Detection with ATIC

The ATIC (Advanced Thin Ionization Calorimeter) balloon-borne ionization calorimeter is well suited to record and identify high energy cosmic ray electrons. The instrument was exposed to high-energy beams at CERN H2 bean-dine in September of 1999. We have simulated the performance of the instrument, and compare the simulations with actual high energy electron exposures at the CERN accelerator. Simulations and measurements do not compare exactly, in detail, but overall the simulations have predicted actual measured behavior quite well.

Chang, J.↗

Monte Carlo simulation of the nuclear-electromagnetic cascade development and the energy response of ionization spectrometers

Modifications to the basic computer program for performing the simulations are reported. The major changes include: (1) extension of the calculations to include the development of cascades initiated by heavy nuclei, (2) improved treatment of the nuclear disintegrations which occur during the interactions of hadrons in heavy absorbers, (3) incorporation of accurate multi-pion final-state cross sections for various interactions at accelerator energies, (4) restructuring of the program logic so that calculations can be made for sandwich-type detectors, and (5) logic modifications related to execution of the program.

Jones, W. V.↗

Heating-related flows in cool solar loops

The effects of spatial and temporal variations in the heating of cool loop models are investigated in an attempt to explain the net redshifts that are observed on the sun. The response of initially static cool loops to changes in the energy input is simulated. For hot loops, it is found that spatially asymmetric changes produce a final steady state that is dynamic, and that spatially symmetric changes produce a final state that is static. Some general properties of cool loop equilibria are discussed, emphasizing the relationship between structure and energy input. The results are unable to explain the net redshifts observed in emission lines formed near 100,000 K on the sun.

Klimchuk, J. A.↗

Statistical Properties of Maximum Likelihood Estimators of Power Law Spectra Information

A simple power law model consisting of a single spectral index, a is believed to be an adequate description of the galactic cosmic-ray (GCR) proton flux at energies below 10(exp 13) eV, with a transition at the knee energy, E(sub k), to a steeper spectral index alpha(sub 2) greater than alpha(sub 1) above E(sub k). The Maximum likelihood (ML) procedure was developed for estimating the single parameter alpha(sub 1) of a simple power law energy spectrum and generalized to estimate the three spectral parameters of the broken power law energy spectrum from simulated detector responses and real cosmic-ray data. The statistical properties of the ML estimator were investigated and shown to have the three desirable properties: (P1) consistency (asymptotically unbiased). (P2) efficiency asymptotically attains the Cramer-Rao minimum variance bound), and (P3) asymptotically normally distributed, under a wide range of potential detector response functions. Attainment of these properties necessarily implies that the ML estimation procedure provides the best unbiased estimator possible. While simulation studies can easily determine if a given estimation procedure provides an unbiased estimate of the spectra information, and whether or not the estimator is approximately normally distributed, attainment of the Cramer-Rao bound (CRB) can only he ascertained by calculating the CRB for an assumed energy spectrum-detector response function combination, which can be quite formidable in practice. However. the effort in calculating the CRB is very worthwhile because it provides the necessary means to compare the efficiency of competing estimation techniques and, furthermore, provides a stopping rule in the search for the best unbiased estimator. Consequently, the CRB for both the simple and broken power law energy spectra are derived herein and the conditions under which they are attained in practice are investigated. The ML technique is then extended to estimate spectra information from an arbitrary number of astrophysics data sets produced by vastly different science instruments. This theory and its successful implementation will facilitate the interpretation of spectral information from multiple astrophysics missions and thereby permit the derivation of superior spectral parameter estimates based on the combination of data sets.

Howell, L. W.↗

Statistical Properties of Maximum Likelihood Estimators of Power Law Spectra Information

A simple power law model consisting of a single spectral index, sigma(sub 2), is believed to be an adequate description of the galactic cosmic-ray (GCR) proton flux at energies below 10(exp 13) eV, with a transition at the knee energy, E(sub k), to a steeper spectral index sigma(sub 2) greater than sigma(sub 1) above E(sub k). The maximum likelihood (ML) procedure was developed for estimating the single parameter sigma(sub 1) of a simple power law energy spectrum and generalized to estimate the three spectral parameters of the broken power law energy spectrum from simulated detector responses and real cosmic-ray data. The statistical properties of the ML estimator were investigated and shown to have the three desirable properties: (Pl) consistency (asymptotically unbiased), (P2) efficiency (asymptotically attains the Cramer-Rao minimum variance bound), and (P3) asymptotically normally distributed, under a wide range of potential detector response functions. Attainment of these properties necessarily implies that the ML estimation procedure provides the best unbiased estimator possible. While simulation studies can easily determine if a given estimation procedure provides an unbiased estimate of the spectra information, and whether or not the estimator is approximately normally distributed, attainment of the Cramer-Rao bound (CRB) can only be ascertained by calculating the CRB for an assumed energy spectrum- detector response function combination, which can be quite formidable in practice. However, the effort in calculating the CRB is very worthwhile because it provides the necessary means to compare the efficiency of competing estimation techniques and, furthermore, provides a stopping rule in the search for the best unbiased estimator. Consequently, the CRB for both the simple and broken power law energy spectra are derived herein and the conditions under which they are stained in practice are investigated.

Howell, L. W., Jr.↗

Evaluation of Stirling engine appendix gap losses

The efficiency of a Stirling engine can be strongly influenced by the heat transfer losses associated with its appendix gap region. The cyclic energy flows in this region are related to the temperature gradient along the piston and cylinder partition walls, the reciprocating motion of the piston, the pressure variation in the Stirling cycle and the leakage flow across the cold-end seal. This paper outlines a numerical model to simulate these cyclic energy flows and the subsequent effort to correlate it with engine test data. The sensitivity of the appendix gap loss to selected parameters and comparisons with test results are presented.

Huang, S. C.↗

Non-Equilibrium Properties from Equilibrium Free Energy Calculations

Calculating free energy in computer simulations is of central importance in statistical mechanics of condensed media and its applications to chemistry and biology not only because it is the most comprehensive and informative quantity that characterizes the eqUilibrium state, but also because it often provides an efficient route to access dynamic and kinetic properties of a system. Most of applications of equilibrium free energy calculations to non-equilibrium processes rely on a description in which a molecule or an ion diffuses in the potential of mean force. In general case this description is a simplification, but it might be satisfactorily accurate in many instances of practical interest. This hypothesis has been tested in the example of the electrodiffusion equation . Conductance of model ion channels has been calculated directly through counting the number of ion crossing events observed during long molecular dynamics simulations and has been compared with the conductance obtained from solving the generalized Nernst-Plank equation. It has been shown that under relatively modest conditions the agreement between these two approaches is excellent, thus demonstrating the assumptions underlying the diffusion equation are fulfilled. Under these conditions the electrodiffusion equation provides an efficient approach to calculating the full voltage-current dependence routinely measured in electrophysiological experiments.

Pohorille, Andrew↗

Parametric Study of Decay of Homogeneous Isotropic Turbulence Using Large Eddy Simulation

Numerical simulations of decaying homogeneous isotropic turbulence are performed with both low-order and high-order spatial discretization schemes. The turbulent Mach and Reynolds numbers for the simulations are 0.2 and 250, respectively. For the low-order schemes we use either second-order central or third-order upwind biased differencing. For higher order approximations we apply weighted essentially non-oscillatory (WENO) schemes, both with linear and nonlinear weights. There are two objectives in this preliminary effort to investigate possible schemes for large eddy simulation (LES). One is to explore the capability of a widely used low-order computational fluid dynamics (CFD) code to perform LES computations. The other is to determine the effect of higher order accuracy (fifth, seventh, and ninth order) achieved with high-order upwind biased WENO-based schemes. Turbulence statistics, such as kinetic energy, dissipation, and skewness, along with the energy spectra from simulations of the decaying turbulence problem are used to assess and compare the various numerical schemes. In addition, results from the best performing schemes are compared with those from a spectral scheme. The effects of grid density, ranging from 32 cubed to 192 cubed, on the computations are also examined. The fifth-order WENO-based scheme is found to be too dissipative, especially on the coarser grids. However, with the seventh-order and ninth-order WENO-based schemes we observe a significant improvement in accuracy relative to the lower order LES schemes, as revealed by the computed peak in the energy dissipation and by the energy spectrum.

Swanson, R. C.↗

Experimental solar heating-cooling system model tests of a full-scale building system

An experimental solar heating and cooling system model has been built and operated, combining elements that are programmable - heating and/or cooling load of a building, collected solar energy - with experimental equipment. The system model was based on the loads and components used in the Solar Building Test Facility (SBTF) which includes a 1394 square meter solar collector field, at NASA-Langley. Operations covered 5 continuous days under summer conditions. For the system model, up to 55% of the simulated collected solar energy was used for the building load. This amount constituted 35% of the building cooling load. Heat loss was significant. If tank heat loss were eliminated, 75% of the collected solar energy would be used. This amount would supply approximately 50% of the building cooling load.

Namkoong, D.↗

Tests of a reduced-scale experimental model of a building solar heating-cooling system

An experimental solar heating and cooling system model has been built and operated, combining elements that are programmable (e.g., heating and cooling load of a building and collected solar energy) with experimental equipment. The experimental system model was based on the loads and components used in the Solar Building Test Facility (SBTF), which includes a 1394 sq m solar collector field at NASA Langley. These tests covered 5 continuous days under summer conditions. For the system model up to 55 percent of the simulated collected solar energy was used for the building load. This amount of solar energy supplied 35 percent of the building cooling load. Heat loss was significant. If tank heat loss were eliminated, which would make it similar to the actual SBTF, 75 percent of the collected solar energy would be used. This amount would supply approximately 50 percent of the building cooling load. A higher fraction of solar energy is possible with a more performance-optimized system.

Namkoong, D.↗

Simulation studies for surfaces and materials strength

Three different investigations based on computer simulations were carried out. Simulation calculations were performed using model potentials with two- and three-body interactions which were represented by the Mie and the Axilrod-Teller potentials, respectively. In the first part, to analyze stability criteria for simple crystals, a parametrical investigation was performed and stability regions for five different crystalline structures were determined with respect to potential energy parameters. To simulate some of the basic surface processes such as diffusion and nucleation, in the second and third parts, silicon surfaces were considered as a specific example. In these studies diffusion of adatoms on stepped surfaces were simulated and energetics associated with the formation of kink sites were calculated.

Halicioglu, Timur↗