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 649 records · Page 36

Analysis of SOFCAL calibration data

Ionization calorimeters determine the total energy of a particle by absorbing, within the mass of the calorimeter, the entire energy of the particle or a significant (and determinable) fraction of it. The JACEE group has used passive ionization calorimetery, which employs photographic emulsions as the detector medium, to measure the charge composition and energy spectrum of cosmic rays up to, and exceeding, 10(exp 14)eV. The SOFIC approach depends similarly on the use of three-dimensional nuclear-electromagnetic shower theory to relate the ionization deposits obtained during a shower to the energy of the particle causing the shower. Bundles of thousands of scintillating optical fibers are read out using image-intensified CCD's. Such an event may be triggered by a fast shower detector placed under the instrument or if only heavy particles are of interest, from a fast primary Cerention detector placed above the calorimeter. In the first study, fast Hammatsu photomultiplier tubes were purchased and tested for possible application for a triggering purpose. In the second study, some refinements have been made to the theoretical treatment of hadronic interactions in the central collision region. These will be helpful in improving the simulations necessary for observations of high energy cosmic ray nuclei with a SOFIC.

Gregory, John C.↗

Orbital simulations of satellite escape/capture and the origin of satellites such as Triton

We investigate satellite escape/capture in the context of the restricted, circular three body problem as applied to the Sun, Neptune, and Triton. We have computed a large number of coplanar prograde and retrograde orbital simulations over a range of initial distances and velocities. The satellite starts at superior conjunction within approximately 2 Hill radii of Neptune and has a velocity orthogonal to the Sun-planet line. Orbits with these initial conditions can be reflected with respect to time, so an escape is simply the reverse of a capture. We numerically integrate the equations of motion to compute the satellite's position until it escapes, collides with Neptune, or after 100 planetary years fails to escape, when computations cease. The initial distance x and velocity v in the restricted problem uniquely define the Jacobi constant C, a conserved energy-like quantity. Plots of the simulation outcomes in the prograde and retrograde C, x phase spaces reveal distinct zones in which temporary satellites approach the planet closely enough that permanent capture can be effected by gas drag with a protoplanetary nebula or by collision with a pre-existing satellite. Single and double close-flybys constitute the most common possible capture orbits. Long term multiple flyby orbits occur near the stability limits between bound and unbound orbits, and are more common among retrograde captures.

Benner, Lance A. M.↗

Simulation of the turbulent Rayleigh-Benard problem using a spectral/finite difference technique

The three-dimensional, incompressible Navier-Stokes and energy equations with the Bousinesq assumption have been directly simulated at a Rayleigh number of 3.8 x 10 to the 5th power and a Prandtl number of 0.76. In the vertical direction, wall boundaries were used and in the horizontal, periodic boundary conditions were used. A spectral/finite difference numerical method was used to simulate the flow. The flow at these conditions is turbulent and a sufficiently fine mesh was used to capture all relevant flow scales. The results of the simulation are compared to experimental data to justify the conclusion that the small scale motion is adequately resolved.

Eidson, T. M.↗

Momentum transport near a magnetic X line in collisionless reconnection

Plasma dynamics and momentum transport near an X line during time-dependent magnetic reconnection in a collisionless plasma are investigated based on two-dimensional particle simulations. We find that a weakly skewed velocity distribution is formed near the magnetic X line, leading to the presence of off-diagonal elements of the plasma pressure tensor. Let the reconnection electric field be in the y direction. The gradients of the off-diagonal elements of the pressure tensor can provide a transport of the y momentum. During the normal magnetic reconnection, the momentum transport associated with the off-diagonal terms of the pressure tensor mediates a transfer of the y momentum from the region near the X line to regions outside the X line. A period of 'reverse magnetic reconnection,' during which the plasma kinetic energy is converted into magnetic energy, is also observed in the simulation. When reverse reconnection occurs, the gradients of the off-diagonal pressure tensor elements can mediate a transfer of y momentum into the X line. It is found that the inertial term also plays a significant role in the force balance near the magnetic X line. An explanation for the origin of the off-diagonal pressure terms is also given in this paper.

Cai, H. J.↗

Space-based Ornithology-Studying Bird Migration and Environmental Change in North America

Natural fluctuations in the availability of critical stopover sites coupled with anthropogenic destruction of wetlands, land-use change, and anticipated losses due to climate change present migratory birds with a formidable challenge. We have developed an individual-based, spatially explicit bird migration model that simulates the migration routes, timing and energy budgets of individual birds under dynamic weather and land surface conditions. Our model incorporates biophysical constraints, individual bird energy status, bird behavior, and flight aerodynamics. We model the speed, direction, and timing of individual birds moving through a user specified Lagrangian grid. The model incorporates environmental properties including wind speed and direction, topography, dynamic hydrologic properties of the landscape, and environmental suitability. The model is driven by important variables estimated from satellite observations of the land surface, by data assimilation products from weather and climate models, and biological field data. We illustrate the use of the model to study the impact of both short- and long-term environmental variatios, e.g. climate, drought, anthropogenic, on migration timing (phenology), spatial pattern, and fitness (survival and reproductive success). We present several theoretical simulations of the spring migration of Pectoral Sandpiper (Calidris melanotos) in North America with emphasis on the Central flyway from the Gulf of Mexico to Alaska.

Smith, James↗

Robotic Specialization in Autonomous Robotic Structural Assembly

Robotic in-space assembly of large space structures is a long-term NASA goal to reduce launch costs and enable larger scale missions. Recently, researchers have proposed using discrete lattice building blocks and co-designed robots to build high-performance, scalable primary structure for various on-orbit and surface applications. These robots would locomote on the lattice and work in teams to build and reconfigure building-blocks into functional structure. However, the most reliable and efficient robotic system architecture, characterized by the number of different robotic 'species' and the allocation of functionality between species, is an open question. To address this problem, we decompose the robotic building-block assembly task into functional primitives and, in simulation, study the performance of the the variety of possible resulting architectures. For a set consisting of five process types (move self, move block, move friend, align bock, fasten block), we describe a method of feature space exploration and ranking based on energy and reliability cost functions. The solution space is enumerated, filtered for unique solutions, and evaluated against energy and reliability cost functions for various simulated build sizes. We find that a 2 species system, dividing the five mentioned process types between one unit cell transport robot and one fastening robot, results in the lowest energy cost system, at some cost to reliability. This system enables fastening functionality to occupy the build front while reducing the need for that functional mass to travel back and forth from a feed station. Because the details of a robot design affect the weighting and final allocation of functionality, a sensitivity analysis was conducted to evaluate the effect of changing mass allocations on architecture performance. Future systems with additional functionalities such as repair, inspection, and others may use this process to analyze and determine alternative robot architectures.

Bernus, Borbala↗

Evaluation of Material Models within LS-DYNA(Registered TradeMark) for a Kevlar/Epoxy Composite Honeycomb

The performance of material models to simulate a novel composite honeycomb Deployable Energy Absorber (DEA) was evaluated using the nonlinear explicit dynamic finite element code LS-DYNA(Registered TradeMark). Prototypes of the DEA concept were manufactured using a Kevlar/Epoxy composite material in which the fibers are oriented at +/-45 degrees with respect to the loading axis. The development of the DEA has included laboratory tests at subcomponent and component levels such as three-point bend testing of single hexagonal cells, dynamic crush testing of single multi-cell components, and impact testing of a full-scale fuselage section fitted with a system of DEA components onto multi-terrain environments. Due to the thin nature of the cell walls, the DEA was modeled using shell elements. In an attempt to simulate the dynamic response of the DEA, it was first represented using *MAT_LAMINATED_COMPOSITE_FABRIC, or *MAT_58, in LS-DYNA. Values for each parameter within the material model were generated such that an in-plane isotropic configuration for the DEA material was assumed. Analytical predictions showed that the load-deflection behavior of a single-cell during three-point bending was within the range of test data, but predicted the DEA crush response to be very stiff. In addition, a *MAT_PIECEWISE_LINEAR_PLASTICITY, or *MAT_24, material model in LS-DYNA was developed, which represented the Kevlar/Epoxy composite as an isotropic elastic-plastic material with input from +/-45 degrees tensile coupon data. The predicted crush response matched that of the test and localized folding patterns of the DEA were captured under compression, but the model failed to predict the single-cell three-point bending response.

Polanco, Michael A.↗

Aerosol-Ice Formation Closure Pilot Study (AEROICSTUDY) Field Campaign Report

The aerosol-ice formation closure pilot study (AEROICESTUDY) was a field campaign conducted at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility at the Southern Great Plains (SGP) observatory from October 7 to 28, 2019. The purpose of this campaign was to test a field observational approach for conducting an aerosol-ice formation closure study. In other words, the goal was to predict the number of ice nucleating particles (INPs) from the measurements of physical and chemical aerosol properties and compare those to measured INP number concentrations. This project was motivated by the fact that the aerosol community has widely conducted aerosol radiative closure, aerosol-cloud condensation nuclei (CCN), and CCN-droplet closure studies to test the physical models and parameterizations that cloud-resolving and climate models rely on to perform reliable simulations of the Earth system and energy budget. However, very few closure studies related to INPs have been conducted, and to our knowledge, none using robust, size-resolved, ambient aerosol composition measurements as model inputs.

aerosol-ice formation↗

Aerosol-Ice Formation Closure Pilot Study (AEROICESTUDY) Field Campaign Report

The aerosol-ice formation closure pilot study (AEROICESTUDY) was a field campaign conducted at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility at the Southern Great Plains (SGP) observatory from October 7 to 28, 2019. The purpose of this campaign was to test a field observational approach for conducting an aerosol-ice formation closure study. In other words, the goal was to predict the number of ice nucleating particles (INPs) from the measurements of physical and chemical aerosol properties and compare those to measured INP number concentrations. This project was motivated by the fact that the aerosol community has widely conducted aerosol radiative closure, aerosol-cloud condensation nuclei (CCN), and CCN-droplet closure studies to test the physical models and parameterizations that cloud-resolving and climate models rely on to perform reliable simulations of the Earth system and energy budget. However, very few closure studies related to INPs have been conducted, and to our knowledge, none using robust, size-resolved, ambient aerosol composition measurements as model inputs

AEROICESTUDY↗

Quantification of Energy Release in Composite Structures

Energy release rate is usually suggested as a quantifier for assessing structural damage tolerance. Computational prediction of energy release rate is based on composite mechanics with micro-stress level damage assessment, finite element structural analysis and damage progression tracking modules. This report examines several issues associated with energy release rates in composite structures as follows: Chapter I demonstrates computational simulation of an adhesively bonded composite joint and validates the computed energy release rates by comparison with acoustic emission signals in the overall sense. Chapter II investigates the effect of crack plane orientation with respect to fiber direction on the energy release rates. Chapter III quantifies the effects of contiguous constraint plies on the residual stiffness of a 90 deg ply subjected to transverse tensile fractures. Chapter IV compares ICAN and ICAN/JAVA solutions of composites. Chapter V examines the effects of composite structural geometry and boundary conditions on damage progression characteristics.

Minnetyan, Levon↗

Quantification of Energy Release in Composite Structures

Energy release rate is usually suggested as a quantifier for assessing structural damage tolerance. Computational prediction of energy release rate is based on composite mechanics with micro-stress level damage assessment, finite element structural analysis and damage progression tracking modules. This report examines several issues associated with energy release rates in composite structures as follows: Chapter I demonstrates computational simulation of an adhesively bonded composite joint and validates the computed energy release rates by comparison with acoustic emission signals in the overall sense. Chapter II investigates the effect of crack plane orientation with respect to fiber direction on the energy release rates. Chapter III quantifies the effects of contiguous constraint plies on the residual stiffness of a 90 ply subjected to transverse tensile fractures. Chapter IV compares ICAN and ICAN/JAVA solutions of composites. Chapter V examines the effects of composite structural geometry and boundary conditions on damage progression characteristics.

Minnetyan, Levon↗

Dynamic Simulation of Human Thermoregulation and Heat Transfer for Spaceflight Applications

Models of human thermoregulation and heat transfer date from the early 1970s and have been developed for applications ranging from evaluating thermal comfort in spacecraft and aircraft cabin environments to predicting heat stress during EVAs. Most lumped or compartment models represent the body as an assemblage cylindrical and spherical elements which may be subdivided into layers to describe tissue heterogeneity. Many existing models are of limited usefulness in asymmetric thermal environments, such as may be encountered during an EVA. Conventional whole-body clothing models also limit the ability to describe local surface thermal and evaporation effects in sufficient detail. A further limitation is that models based on a standard man model are not readily scalable to represent large or small subjects. This work describes development of a new human thermal model derived from the 41-node man model. Each segment is divided into four concentric, constant thickness cylinders made up of a central core surrounded by muscle, fat, and skin, respectively. These cylinders are connected by the flow of blood from a central blood pool to each part. The central blood pool is updated at each time step, based on a whole-body energy balance. Results show the model simulates core and surface temperature histories, sweat evaporation and metabolic rates which generally are consistent with controlled exposures of human subjects. Scaling rules are developed to enable simulation of small and large subjects (5th percentile and 95th percentile). Future refinements will include a clothing model that addresses local surface insulation and permeation effects and developing control equations to describe thermoregulatory effects such as may occur with prolonged weightlessness or with aging.

Miller, Thomas R.↗

Principles of magnetohydrodynamic simulation in space plasmas

Attention is given to the philosophical as well as physical principles that are essential to the establishment of MHD simulation studies for solar plasma research, assuming the capabilities of state-of-the-art computers and emphasizing the importance of 'local' MHD simulation. Solar-terrestrial plasma space is divided into several elementary regions where a macroscopic elementary energy conversion process could conceivably occur; the local MHD simulation is defined as self-contained in each of the regions. The importance of, and the difficulties associated with, the boundary condition are discussed in detail. The roles of diagnostics and of the finite difference method are noted.

Sato, T.↗

Dissociation and Internal Excitation of Molecular Nitrogen Due to N + N2 Collisions Using Direct Molecular Simulation

In this work we present a molecular level study of N2+N collisions, focusing on excitation of internal energy modes and non-equilibrium dissociation. The computation technique used here is the direct molecular simulation (DMS) method and the molecular interactions have been modeled using an ab−initio potential energy surface (PES) developed at NASA's Ames Research Center. We carried out vibrational excitation calculations between 5000K and 30000K and found that the characteristic vibrational excitation time for the N + N2 process was an order of magnitude lower than that predicted by the Millikan and White correlation. It is observed that during vibrational excitation the high energy tail of the vibrational energy distribution gets over populated first and the lower energy levels get populated as the system evolves. It is found that the non-equilibrium dissociation rate coefficients for the N + N2 process are larger than those for the N2 + N2 process. This is attributed to the non-equilibrium vibrational energy distributions for the N + N2 process being less depleted than that for the N2 +N2 process. For an isothermal simulation we find that the probability of dissociation goes as 1/T(sub tr) for molecules with internal energy (epsilon(sub int)) less than approximately 9.9eV, while for molecules with epsilon (sub int) greater than 9.9eV the dissociation probability was weakly dependent on translational temperature of the system. We compared non-equilibrium dissociation rate coefficients and characteristic vibrational excitation times obtained by using the ab-initio PES developed at NASA's Ames Research Center to those obtained by using an ab-initio PES developed at the University of Minnesota. Good agreement was found between the macroscopic properties and molecular level description of the system obtained by using the two PESs.

Grover, Maninder S.↗

The Influence of surface binding energy on sputtering in models of the sodium exosphere of Mercury

We have simulated the sodium (Na) exosphere of Mercury to show how the exosphere is affected by the assumed surface binding energy (SBE) of Na in the sputtered component. We constrained ion precipitation onto the surface using distributions for the cusp regions that are consistent with measurements by the MErcury Surface, Space ENvironment, GEochemistry, and Ranging Fast Imaging Plasma Spectrometer instrument. We have simulated sputtering with SBEs of 0.27, 2.6, 4.4, and 7.9 eV, with the lowest value commonly used in exosphere models and the highest from recent molecular dynamics calculations for the Na-bearing feldspar end-member, albite. A gradual change in the exosphere is seen as the yield decreases and the ejecta energy increases with increasing SBE. We describe the corresponding exosphere source functions for ion sputtering (IS), as well as for the previously studied processes of micrometeoroid impact vaporization and photon-stimulated desorption (PSD), along with their release energy distributions and spatial distributions. We have summed the contributions of the various source processes to explain how and when the different sources can be distinguished by observations. The modeled exosphere scale heights range from 72 km for PSD to over 1000 km for IS using a SBE of 7.9 eV. We find that the processes responsible for generating Mercury's Na exosphere are separable by measuring line-of-sight column densities tangent to the planet at various altitudes and positions around the planet. Our initial results are consistent with the Na being sputtered from a high-SBE material such as feldspar, which has been predicted to be abundant on the Mercury's surface.

Rosemary M. Killen↗

Phase diagrams and microscopic structures of (Hg,Cd)Te, (Hg,Zn)Te, and (Cd,Zn)Te alloys

A cluster theory based on the quasi-chemical approximation has been applied to study the local correlation bond-length distribution, and phase diagrams of the II-VI pseudobinary alloys Hg(1 - x)Cd(x)Te, Hg(1 - x)Zn(x)Te, and Cd(1 - x)Zn(x)Te. The cluster energy is calculated by letting it relax in some effective alloy medium and then considering the contributions from the strain and chemical energies. Two different models are presented to simulate the alloy medium. While both models show that all three alloys have nearly random distributions, the signs of the local correlation prove to be sensitive to the alloy medium chosen for the energy calculation. Good agreement is found between experiment and the bond lengths and phase diagrams in both models.

Patrick, R. S.↗

High Energy Boundary Conditions for a Cartesian Mesh Euler Solver

Inlets and exhaust nozzles are common place in the world of flight. Yet, many aerodynamic simulation packages do not provide a method of modelling such high energy boundaries in the flow field. For the purposes of aerodynamic simulation, inlets and exhausts are often fared over and it is assumed that the flow differences resulting from this assumption are minimal. While this is an adequate assumption for the prediction of lift, the lack of a plume behind the aircraft creates an evacuated base region thus effecting both drag and pitching moment values. In addition, the flow in the base region is often mis-predicted resulting in incorrect base drag. In order to accurately predict these quantities, a method for specifying inlet and exhaust conditions needs to be available in aerodynamic simulation packages. A method for a first approximation of a plume without accounting for chemical reactions is added to the Cartesian mesh based aerodynamic simulation package CART3D. The method consists of 3 steps. In the first step, a components approach where each triangle is assigned a component number is used. Here, a method for marking the inlet or exhaust plane triangles as separate components is discussed. In step two, the flow solver is modified to accept a reference state for the components marked inlet or exhaust. In the third step, the flow solver uses these separated components and the reference state to compute the correct flow condition at that triangle. The present method is implemented in the CART3D package which consists of a set of tools for generating a Cartesian volume mesh from a set of component triangulations. The Euler equations are solved on the resulting unstructured Cartesian mesh. The present methods is implemented in this package and its usefulness is demonstrated with two validation cases. A generic missile body is also presented to show the usefulness of the method on a real world geometry.

Pandya, Shishir↗

Hybrid Particle-Element Simulation of Impact on Composite Orbital Debris Shields

This report describes the development of new numerical methods and new constitutive models for the simulation of hypervelocity impact effects on spacecraft. The research has included parallel implementation of the numerical methods and material models developed under the project. Validation work has included both one dimensional simulations, for comparison with exact solutions, and three dimensional simulations of published hypervelocity impact experiments. The validated formulations have been applied to simulate impact effects in a velocity and kinetic energy regime outside the capabilities of current experimental methods. The research results presented here allow for the expanded use of numerical simulation, as a complement to experimental work, in future design of spacecraft for hypervelocity impact effects.

Fahrenthold, Eric P.↗