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 613 records · Page 34

Simulation studies for surfaces and materials strength

Computer simulation studies were conducted to investigate materials properties. All the calculations were carried out using atomic level simulation techniques which are based on semiempirical or model functions. During this project, functions with varying degrees of complexity were derived and employed in simulations. Investigations include covalently bonded materials as well as systems involving metal atoms. For small clusters calculated results provided information on various energy- and structure-related properties along with vibrational characteristics. Also, energy barriers for configurational transitions were calculated for selected cases. Simulation calculations for surfaces produced new results in areas related to surface energies, reconstructions and relaxations, surface defects, surface stresses as well as adsorption and nucleation processes. Simulations were also performed on bulk materials. Calculated results provided an atomic level understanding on energetics and structures of point defects, crystal stability, elastic properties, and materials strength for various systems. Calculations involving polymeric materials include studies of polar polymer melts and polymer/solid interfaces. Simulations employing specially developed codes provided significant information about energetics and conformational characteristics of different polymeric chains.

Halicioglu, Timur↗

Spatial distribution and yield of DNA double-strand breaks induced by 3-7 MeV helium ions in human fibroblasts

Accelerated helium ions with mean energies at the target location of 3-7 MeV were used to simulate alpha-particle radiation from radon daughters. The experimental setup and calibration procedure allowed determination of the helium-ion energy distribution and dose in the nuclei of irradiated cells. Using this system, the induction of DNA double-strand breaks and their spatial distributions along DNA were studied in irradiated human fibroblasts. It was found that the apparent number of double-strand breaks as measured by a standard pulsed-field gel assay (FAR assay) decreased with increasing LET in the range 67-120 keV/microm (corresponding to the energy of 7-3 MeV). On the other hand, the generation of small and intermediate-size DNA fragments (0.1-100 kbp) increased with LET, indicating an increased intratrack long-range clustering of breaks. The fragment size distribution was measured in several size classes down to the smallest class of 0.1-2 kbp. When the clustering was taken into account, the actual number of DNA double-strand breaks (separated by at least 0.1 kbp) could be calculated and was found to be in the range 0.010-0.012 breaks/Mbp Gy(-1). This is two- to threefold higher than the apparent yield obtained by the FAR assay. The measured yield of double-strand breaks as a function of LET is compared with theoretical Monte Carlo calculations that simulate the track structure of energy depositions from helium ions as they interact with the 30-nm chromatin fiber. When the calculation is performed to include fragments larger than 0.1 kbp (to correspond to the experimental measurements), there is good agreement between experiment and theory.

NASA Program Biomedical Research and Countermeasur↗

Axial focusing of impact energy in the Earth's interior: Proof-of-principle tests of a new hypothesis

A causal link between major impact events and global processes would probably require a significant change in the thermal state of the Earth's interior, presumably brought about by coupling of impact energy. One possible mechanism for such energy coupling from the surface to the deep interior would be through focusing due to axial symmetry. Antipodal focusing of surface and body waves from earthquakes is a well-known phenomenon which has previously been exploited by seismologists in studies of the Earth's deep interior. Antipodal focusing from impacts on the Moon, Mercury, and icy satellites has also been invoked by planetary scientists to explain unusual surface features opposite some of the large impact structures on these bodies. For example, 'disrupted' terrains have been observed antipodal to the Caloris impact basis on Mercury and Imbrium Basin on the Moon. Very recently there have been speculations that antipodal focusing of impact energy within the mantle may lead to flood basalt and hotspot activity, but there has not yet been an attempt at a rigorous model. A new hypothesis was proposed and preliminary proof-of-principle tests for the coupling of energy from major impacts to the mantle by axial focusing of seismic waves was performed. Because of the axial symmetry of the explosive source, the phases and amplitudes are dependent only on ray parameter (or takeoff angle) and are independent of azimuthal angle. For a symmetric and homogeneous Earth, all the seismic energy radiated by the impact at a given takeoff angle will be refocused (minus attenuation) on the axis of symmetry, regardless of the number of reflections and refractions it has experienced. Mantle material near the axis of symmetry will experience more strain cycles with much greater amplitude than elsewhere and will therefore experience more irreversible heating. The situation is very different than for a giant earthquake, which in addition to having less energy, has an asymmetric focal mechanism and a larger area. Two independent proof-of-principle approaches were used. The first makes use of seismic simulations, which are being performed with a realistic Earth model to determine the degree of focusing along the axis and to estimate the volume of material, if any, that experiences significant irreversible heating. The second involves two-dimensional hydrodynamic code simulations to determine the stress history, internal energy, and temperature rise as a function of radius along the axis.

Boslough, M. B.↗

New approach in direct-simulation of gas mixtures

Results are reported for an investigation of a new direct-simulation Monte Carlo method by which energy transfer and chemical reactions are calculated. The new method, which reduces to the variable cross-section hard sphere model as a special case, allows different viscosity-temperature exponents for each species in a gas mixture when combined with a modified Larsen-Borgnakke phenomenological model. This removes the most serious limitation of the usefulness of the model for engineering simulations. The necessary kinetic theory for the application of the new method to mixtures of monatomic or polyatomic gases is presented, including gas mixtures involving chemical reactions. Calculations are made for the relaxation of a diatomic gas mixture, a plane shock wave in a gas mixture, and a chemically reacting gas flow along the stagnation streamline in front of a hypersonic vehicle. Calculated results show that the introduction of different molecular interactions for each species in a gas mixture produces significant differences in comparison with a common molecular interaction for all species in the mixture. This effect should not be neglected for accurate DSMC simulations in an engineering context.

Chung, Chan-Hong↗

Numerical investigation of the control of the secondary instability process in boundary layers

A numerical model has been developed for investigating boundary layer transition control for a flat plate boundary layer. Active control of a periodically forced boundary layer in an incompressible fluid is studied using surface heating techniques. The spatially evolving boundary layer is simulated. The Navier-Stokes and energy equations are integrated using a fully implicit finite difference/spectral method. Temperature perturbations are introduced locally along finite heater strips to directly attenuate the instability waves in the flow. A feedback control loop is employed in which a downstream sensor is used to monitor wall shear stress fluctuations. Active control of small amplitude two-dimensional and three-dimensional disturbances is numerically simulated. With proper phase control, in-phase reinforcement and out-of-phase attenuation are demonstrated. A receptivity study of the localized temperature perturbations is made. It is shown that narrow heater strips are more receptive in that they maximize the amplitude level of the disturbances in the flow. Active control of the early stages of the fundamental breakdown process is also numerically simulated. Control is achieved with either two-dimensional or three-dimensional control inputs.

Kral, L. D.↗

Near-wall k-epsilon turbulence modeling

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

Mansour, N. N.↗

Free surfaces and multilayer interfaces in the GaAs/AlAs system

Semiempirical potential energy functions have been utilized for a variety of calculations in the Ga-Al-As system. Surface energies have been calculated for several orientations of GaAs. Ledge energies for the GaAs (0 0 -1) (As-terminated) surface show long-range interaction effects with the ledge energy increasing with spacing. GaAs (001)/AlAs (001) superlattices have been simulated for a range of interlayer spacings with the excess interfacial energy per interlayer increasing from 5 erg/sq cm at an interlayer spacing of 1 molecular layer (5.8 A) to 50 erg/sq cm at an interlayer spacing of 18 molecular layers (103.8 A).

Choi, D. K.↗

Near-wall k-epsilon turbulence modeling

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

Mansour, N. N.↗

Intensity variations in the interplanetary magnetic field measured by Voyager 2 and the 11-year solar cycle modulation of galactic cosmic rays

New evidence is presented to support the hypothesis that the 11-yr solar cycle modulation of galactic cosmic rays is caused by strong diffusion inside long-lived merged interaction regions. To test this hypothesis, the 1D force-field approximation of the cosmic ray modulation equation is solved. It is assumed that a constant solar wind speed convects magnetic field compressions and rarefactions unchanged through a model heliosphere. The result is a reasonable simulation of the integrated high-energy cosmic ray intensity profile from about 1982 to mid-1989. This period encompasses both the full recovery portion of the last profile from about 1982 to mid-1989. This model responds to the Voyager 2 magnetic field data by correctly timing the beginning of the new modulation cycle in late 1987. It is concluded that the present hypothesis is consistent with the results of this simulation.

Perko, J. S.↗

Elastic Model Transitions: A Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes a hybrid approach for determining physical displacements by augmenting the original quadratically constrained least squares (LSQI) algorithm with Direct Shape Mapping (DSM) and modifying the energy constraints. The approach presented is applicable to simulation of the elastic behavior of launch vehicles and other structures that utilize discrete LTI finite element model (FEM) derived mode sets (eigenvalues and eigenvectors) that are propagated throughout time. The time invariant nature of the elastic data presents a problem of how to properly transition elastic states from the prior to the new model while preserving motion across the transition and ensuring there is no truncation or excitation of the system. A previous approach utilizes a LSQI algorithm with an energy constraint to effect smooth transitions between eigenvector sets with no requirement that the models be of similar dimension or have any correlation. This approach assumes energy is conserved across the transition, which results in significant non-physical transients due to changing quasi-steady state energy between mode sets, a phenomenon seen when utilizing a truncated mode set. The computational burden of simulating a full mode set is significant so a subset of modes is often selected to reduce run time. As a result of this truncation, energy between mode sets may not be constant and solutions across transitions could produce non-physical transients. In an effort to abate these transients an improved methodology was developed based on the aforementioned approach, but this new approach can handle significant changes in energy across mode set transitions. It is proposed that physical velocities due to elastic behavior be solved for using the LSQI algorithm, but solve for displacements using a two-step process that independently addresses the quasi-steady-state and non-steady-state contributions to the elastic displacement. For structures subject to large external forces, such as thrust or atmospheric drag, it is imperative to capture these forces when solving for elastic displacement. To simplify the mathematical formulation, assumptions are made regarding mass matrix normalization, constant external forcing, and constant viscous damping. These simplifications allow for direct solutions to the quasi-steady-state displacements through a process titled Direct Shape Mapping. DSM solves for the displacements using the eigenvalues of the elastic modes and the external forcing and returns a set of elastic displacements dictated by the eigenvectors of the post-transition mode set. For the non-steady-state contributions to displacement we formulate a LSQI problem that is constrained by energy of the non-steady state terms. The contributions from the quasi-steady-state and non-steady state solutions are then combined to obtain the physical displacements associated with the new set of eigenvectors. Results for the LSQI-DSM approach show significant reduction/complete removal of transients across mode set transitions while maintaining elastic motion from the prior state. For time propagation applications employing discrete elastic models that need to be transitioned in time and where running with full a full mode set is not feasible, the method developed offers a practical solution to simulating vehicle elasticity.

Hannan, Mike R.↗

Combustion efficiency and altitude operational limits of three liquid hydrocarbon fuels having high volumetric energy content in a J33 single combustor

Combustion efficiency and altitude operational limits were determined in a J33 single combustor for AN-F-58 fuel and three liquid hydrocarbon fuels having high volumetric energy content (decalin, tetralin, and monomethylnaphthalene) at simulated altitude and combustor inlet-air conditions. At the conditions investigated, the combustion efficiency for the four fuels generally decreased with an increase in volumetric energy content. The altitude operational limits for decalin and tetralin fuels were higher than for AN-F-58 fuel; monomethylnaphthalene fuel gave the lowest altitude operational limit.

Stricker, Edward G↗

Ballistic Limit Equation for Single Wall Titanium

Hypervelocity impact tests and hydrocode simulations were used to determine the ballistic limit equation (BLE) for perforation of a titanium wall, as a function of wall thickness. Two titanium alloys were considered, and separate BLEs were derived for each. Tested wall thicknesses ranged from 0.5mm to 2.0mm. The single-wall damage equation of Cour-Palais [ref. 1] was used to analyze the Ti wall's shielding effectiveness. It was concluded that the Cour-Palais single-wall equation produced a non-conservative prediction of the ballistic limit for the Ti shield. The inaccurate prediction was not a particularly surprising result; the Cour-Palais single-wall BLE contains shield material properties as parameters, but it was formulated only from tests of different aluminum alloys. Single-wall Ti shield tests were run (thicknesses of 2.0 mm, 1.5 mm, 1.0 mm, and 0.5 mm) on Ti 15-3-3-3 material custom cut from rod stock. Hypervelocity impact (HVI) tests were used to establish the failure threshold empirically, using the additional constraint that the damage scales with impact energy, as was indicated by hydrocode simulations. The criterion for shield failure was defined as no detached spall from the shield back surface during HVI. Based on the test results, which confirmed an approximately energy-dependent shield effectiveness, the Cour-Palais equation was modified.

Ratliff, J. M.↗

Effects of a giant impact on Uranus

The effects of a giant impact on Uranus with respect to the axis tilt of Uranus and its satellites are discussed. The simulations of possible giant impacts were carried out using Cray supercomputers. The technique used is called smooth particle hydrodynamics (SPH). In this technique, the material in the proto-Uranus planet and in the impactor is divided into a large number of particles which can overlap one another so that local averages over these particles determine density and pressure in the problem, and the particles themselves have their own temperatures and internal energies. During the course of the simulation, these particles move around under the influence of the forces acting on them: gravity and pressure gradients. The results of model simulations are presented.

Slattery, W. L.↗

SIBatt-3D: In-Space/On-Surface 3D Printing of Sodium Ion Batteries from ISRU Materials

Constructed more than 20 years ago, the International Space Station’s primary power system originally used nickel-hydrogen batteries with a lifetime of 6.5 years, until NASA began the process of replacing them in 2016 with lithium-ion batteries with a lifetime of 10 years. The demanding and costly process was accomplished after four flights of the Japanese H-II Transfer Vehicle cargo spacecraft (with a cost of about $10,000 per pound of payload), and 13 different astronauts conducting 14 spacewalks. Besides utilization in the ISS, rechargeable batteries are present in many space applications: they are installed in exploration robots, life support systems and in portable communication devices, to mention some. In this context, this project is focused on the in-space manufacturing of shape-conformable batteries using in-situ resources, and aims to address the NASA’s gaps related to the development of next generation of energy storage devices (TX03), as well as in-space manufacturing and in-situ resource utilization (TX07). The proposed work also tackles the HEOMD’s objectives targeting the in-space additive manufacturing (AM) from Lunar/Martian materials (regolith as AM feedstock) to reinvigorate America’s Human Space Exploration Program (SPD-1). This project is in direct alignment with the STMD’s objectives to demonstrate in-space autonomous manufacturing and assembly of complete systems by 2030, and to enable humans to live and explore in space and on planetary surfaces by 2040 thanks to in-space habitation, infrastructure development and in-situ resource utilization (ST1 and ST5). Manufacturing of shape conformable batteries directly in-space and using in-situ resources would also contribute to reducing the payload weight and volume (TX12) for future missions, thus reducing risk for long term Mars missions where rapid resupply is logistically infeasible. Nowadays, commercial batteries consist of stacked two-dimensional (2D) sheets, which are only manufactured in restricted geometries (cylindrical and coin cell). Evolving from conventional 2D, complex 3D battery architectures have been proven to increase the electrochemical active surface area and ion diffusion path, leading to improved areal energy density and power performance. This tendency was illustrated in our recent in-depth modeling studies by simulating a classical Ragone plot exhibiting the energy-power relationship. Our team demonstrated through modeling that a complex gyroidal 3D printed battery architecture exhibits significantly improved power performances (>150% at the current density of 6C; full discharge in 10 minutes) in comparison to a traditional 3D printed planar geometry. Motivated by these results and as the fabrication of intricate 3D battery design is only possible experimentally thanks to the geometric freedom offered by additive manufacturing (AM), our team has already initiated leveraging thermoplastic material extrusion at the laboratory scale. While 3D printing of batteries is relatively recent (2013), it has witnessed a growing interest during the last recent years, as next-generation shape-conformable 3D batteries can be co-designed with the system. Consequently, dead-volume and mass brought from Earth are minimized, in addition to improved battery performance, in alignment with the aforementioned NASA’s objectives. Further, while this project is specifically dedicated to batteries, it lends itself towards the maturation of in-space manufacturing via 3D printing using in-situ resources, stated in HEOMD and STMD goals.

In-Space Manufacturing↗

Modeling and Simulation of a Parallel Hybrid-Electric Propulsion System - Electrified Powertrain Flight Demonstration (EPFD) Program

Electrified aircraft propulsion concepts have been proposed to meet aggressive future performance and environmental goals for the next generation of aircraft. However, electrified aircraft present a unique modeling and simulation challenge as they introduce multiple energy sources to the propulsion system, providing various means to meet thrust requirements, compared to conventional gas turbine propulsion architectures where only fuel is available. Additionally, the introduction of an electric powertrain to the existing system enables multiple electrified flight modes to exist (i.e. eTaxi, climb boost, takeoff boost, etc.), further increasing the complexity of the modeling environment. As part of the Electrified Powertrain Flight Demonstration program, this paper presents a modeling and simulation framework for a parallel hybrid-electric propulsion concept using the Environmental Design Space simulation tool. Electrical components are modeled in NPSS, and an overall sizing methodology is introduced. Finally, various operational modes of the electric powertrain are modeled and tested and their impact on key performance parameters is evaluated.

hybrid electric propulsion↗

Spacecraft Solar Particle Event (SPE) Shielding: Shielding Effectiveness as a Function of SPE model as Determined with the FLUKA Radiation Transport Code

Analysis of both satellite and surface neutron monitor data demonstrate that the widely utilized Exponential model of solar particle event (SPE) proton kinetic energy spectra can seriously underestimate SPE proton flux, especially at the highest kinetic energies. The more recently developed Band model produces better agreement with neutron monitor data ground level events (GLEs) and is believed to be considerably more accurate at high kinetic energies. Here, we report the results of modeling and simulation studies in which the radiation transport code FLUKA (FLUktuierende KAskade) is used to determine the changes in total ionizing dose (TID) and single-event environments (SEE) behind aluminum, polyethylene, carbon, and titanium shielding masses when the assumed form (i. e., Band or Exponential) of the solar particle event (SPE) kinetic energy spectra is changed. FLUKA simulations have fully three dimensions with an isotropic particle flux incident on a concentric spherical shell shielding mass and detector structure. The effects are reported for both energetic primary protons penetrating the shield mass and secondary particle showers caused by energetic primary protons colliding with shielding mass nuclei. Our results, in agreement with previous studies, show that use of the Exponential form of the event

Koontz, Steve↗

EAS accompanied by gamma-families at Mt. Norikura and comparison with Monte Carlo simulation

The experimental data of extensive air showers (EAS) accompanied by gamma-families, with total energy greater than 10 TeV, were compared with a Monte Carlo simulation with a rising cross section proportional to E sup 0.04 for the p-air inelastic cross section. It is found that the absolute intensity of size spectrum of such EAS is strongly affected by the primary protons intensity at 10 to the 15th power approx. 10 to the 16th power eV region and the experimental size spectrum agrees with the simulated spectra for the p-poor primary composition better than the p-rich one.

Shima, M.↗

Simulation and Measurement of Absorbed Dose from 137 Cs Gammas Using a Si Timepix Detector

The TimePix readout chip is a hybrid pixel detector with over 65k independent pixel elements. Each pixel contains its own circuitry for charge collection, counting logic, and readout. When coupled with a Silicon detector layer, the Timepix chip is capable of measuring the charge, and thus energy, deposited in the Silicon. Measurements using a NIST traceable 137Cs gamma source have been made at Johnson Space Center using such a Si Timepix detector, and this data is compared to simulations of energy deposition in the Si layer carried out using FLUKA.

Stoffle, Nicholas↗