Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerator modelling and simulations”

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 271 records · Page 15

Simulating Sources of Superstorm Plasmas

We evaluated the contributions to magnetospheric pressure (ring current) of the solar wind, polar wind, auroral wind, and plasmaspheric wind, with the surprising result that the main phase pressure is dominated by plasmaspheric protons. We used global simulation fields from the LFM single fluid ideal MHD model. We embedded the Comprehensive Ring Current Model within it, driven by the LFM transpolar potential, and supplied with plasmas at its boundary including solar wind protons, polar wind protons, auroral wind O+, and plasmaspheric protons. We included auroral outflows and acceleration driven by the LFM ionospheric boundary condition, including parallel ion acceleration driven by upward currents. Our plasmasphere model runs within the CRCM and is driven by it. Ionospheric sources were treated using our Global Ion Kinetics code based on full equations of motion. This treatment neglects inertial loading and pressure exerted by the ionospheric plasmas, and will be superceded by multifluid simulations that include those effects. However, these simulations provide new insights into the respective role of ionospheric sources in storm-time magnetospheric dynamics.

Fok, Mei-Ching↗

Simulation Accelerator

Under a NASA SBIR (Small Business Innovative Research) contract, (NAS5-30905), EAI Simulation Associates, Inc., developed a new digital simulation computer, Starlight(tm). With an architecture based on the analog model of computation, Starlight(tm) outperforms all other computers on a wide range of continuous system simulation. This system is used in a variety of applications, including aerospace, automotive, electric power and chemical reactors.

Source record↗

Impact Testing and Simulation of a Sinusoid Foam Sandwich Energy Absorber

A sinusoidal-shaped foam sandwich energy absorber was developed and evaluated at NASA Langley Research Center through multi-level testing and simulation performed under the Transport Rotorcraft Airframe Crash Testbed (TRACT) research project. The energy absorber, designated the "sinusoid," consisted of hybrid carbon- Kevlar® plain weave fabric face sheets, two layers for each face sheet oriented at +/-45deg with respect to the vertical or crush direction, and a closed-cell ELFOAM(TradeMark) P200 polyisocyanurate (2.0-lb/ft3) foam core. The design goal for the energy absorber was to achieve an average floor-level acceleration of between 25- and 40-g during the full-scale crash test of a retrofitted CH-46E helicopter airframe, designated TRACT 2. Variations in the design were assessed through quasi-static and dynamic crush testing of component specimens. Once the design was finalized, a 5-ft-long subfloor beam was fabricated and retrofitted into a barrel section of a CH-46E helicopter. A vertical drop test of the barrel section was conducted onto concrete to evaluate the performance of the energy absorber prior to retrofit into TRACT 2. Finite element models were developed of all test articles and simulations were performed using LSDYNA ®, a commercial nonlinear explicit transient dynamic finite element code. Test analysis results are presented for the sinusoid foam sandwich energy absorber as comparisons of load-displacement and acceleration-time-history responses, as well as predicted and experimental structural deformations and progressive damage for each evaluation level (component testing through barrel section drop testing).

Jackson, Karen E.↗

Three-dimensional calculation of shuttle charging in polar orbit

The charged particles environment in polar orbit can be of sufficient intensity to cause spacecraft charging. In order to gain a quantitative understanding of such effects, the Air Force is developing POLAR, a computer code which simulates in three dimensions the electrical interaction of large space vehicles with the polar ionospheric plasma. It models the physical processes of wake generation, ambient ion collection, precipitating auroral electron fluxes, and surface interactions, including secondary electron generation and backscattering, which lead to vehicle charging. These processes may be followed dynamically on a subsecond timescale so that the rapid passage through intense auroral arcs can be simulated. POLAR models the ambient plasma as isotropic Maxwellian electrons and ions (0+, H+), and allows for simultaneous precipitation of power-law, energetic Maxwellian, and accelerated Gaussian distributions of electrons. Magnetic field effects will be modeled in POLAR but are currently ignored.

Cooke, D. L.↗

Analytical determination of space station response to crew motion and design of suspension system for microgravity experiments

The objective of this investigation is to make analytical determination of the acceleration produced by crew motion in an orbiting space station and define design parameters for the suspension system of microgravity experiments. A simple structural model for simulation of the IOC space station is proposed. Mathematical formulation of this model provides the engineers a simple and direct tool for designing an effective suspension system.

Liu, F. C.↗

Analytical determination of space station response to crew motion and design of suspension system for microgravity experiments

The objective of this investigation is to make analytical determination of the acceleration produced by crew motion in an orbiting space station and define design parameters for the suspension system of microgravity experiments. A simple structural model for simulation of the IOC space station is proposed. Mathematical formulation of this model provides the engineers a simple and direct tool for designing an effective suspension system.

Liu, Frank C.↗

Nonlinear maneuver autopilot for the F-15 aircraft

A methodology is described for the development of flight test trajectory control laws based on singular perturbation methodology and nonlinear dynamic modeling. The control design methodology is applied to a detailed nonlinear six degree-of-freedom simulation of the F-15 and results for a level accelerations, pushover/pullup maneuver, zoom and pushover maneuver, excess thrust windup turn, constant thrust windup turn, and a constant dynamic pressure/constant load factor trajectory are presented.

Menon, P. K. A.↗

Nonaxisymmetric evolution in protostellar disks

We present a two-dimensional, multigridded hydrodynamical simulation of the collapse of an axisymmetric, rotating, 1 solar mass protostellar cloud, which forms a resolved, hydrotastic disk. The code includes the effects of physical viscosity, radiative transfer and radiative acceleration but not magnetic fields. We examine how the disk is affected by the inclusion of turbulent viscosity by comparing a viscous simulation with an inviscid model evolved from the same initial conditions, and we derive a disk evolutionary timescale on the order of 300,000 years if alpha = 0.01. Effects arising from non-axisymmetric gravitational instabilities in the protostellar disk are followed with a three-dimensional SPH code, starting from the two-dimensional structure. We find that the disk is prone to a series of spiral instabilities with primary azimulthal mode number m = 1 and m = 2. The torques induced by these nonaxisymmetric structures elicit material transport of angular momentum and mass through the disk, readjusting the surface density profile toward more stable configurations. We present a series of analyses which characterize both the development and the likely source of the instabilities. We speculate that an evolving disk which maintains a minimum Toomre Q-value approximately 1.4 will have a total evolutionary span of several times 10(exp 5) years, comparable to, but somewhat shorter than the evolutionary timescale resulting from viscous turbulence alone. We compare the evolution resulting from nonaxisymmetric instabilities with solutions of a one-dimensional viscous diffusion equation applied to the initial surface density and temperature profile. We find that an effective alpha-value of 0.03 is a good fit to the results of the simulation. However, the effective alpha will depend on the minimum Q in the disk at the time the instability is activated. We argue that the major fraction of the transport characterized by the value of alpha is due to the action of gravitational torques, and does not arise from inherent viscosity within the smoothed particle hydrodynamics method.

Laughlin, Gregory↗

Mechanisms for the Dissipation of Alfven Waves in Near-Earth Space Plasma

Alfven waves are a major mechanism for the transport of electromagnetic energy from the distant part of the magnetosphere to the near-Earth space. This is especially true for the auroral and polar regions of the Earth. However, the mechanisms for their dissipation have remained illusive. One of the mechanisms is the formation of double layers when the current associated with Alfven waves in the inertial regime interact with density cavities, which either are generated nonlinearly by the waves themselves or are a part of the ambient plasma turbulence. Depending on the strength of the cavities, weak and strong double layers could form. Such double layers are transient; their lifetimes depend on that of the cavities. Thus they impulsively accelerate ions and electrons. Another mechanism is the resonant absorption of broadband Alfven- wave noise by the ions at the ion cyclotron frequencies. But this resonant absorption may not be possible for the very low frequency waves, and it may be more suited for electromagnetic ion cyclotron waves. A third mechanism is the excitation of secondary waves by the drifts of electrons and ions in the Alfven wave fields. It is found that under suitable conditions, the relative drifts between different ion species and/or between electrons and ions are large enough to drive lower hybrid waves, which could cause transverse accelerations of ions and parallel accelerations of electrons. This mechanism is being further studied by means of kinetic simulations using 2.5- and 3-D particle-in-cell codes. The ongoing modeling efforts on space weather require quantitative estimates of energy inputs of various kinds, including the electromagnetic energy. Our studies described here contribute to the methods of determining the estimates of the input from ubiquitous Alfven waves.

Singh, Nagendra↗

Universal Responses of Cyclic-Oxidation Models Studied

Oxidation is an important degradation process for materials operating in the high-temperature air or oxygen environments typical of jet turbine or rocket engines. Reaction of the combustion gases with the component material forms surface layer scales during these oxidative exposures. Typically, the instantaneous rate of reaction is inversely proportional to the existing scale thickness, giving rise to parabolic kinetics. However, more realistic applications entail periodic startup and shutdown. Some scale spallation may occur upon cooling, resulting in loss of the protective diffusion barrier provided by a fully intact scale. Upon reheating, the component will experience accelerated oxidation due to this spallation. Cyclic-oxidation testing has, therefore, been a mainstay of characterization and performance ranking for high-temperature materials. Models simulate this process by calculating how a scale spalls upon cooling and regrows upon heating (refs. 1 to 3). Recently released NASA software (COSP for Windows) allows researchers to specify a uniform layer or discrete segments of spallation (ref. 4). Families of model curves exhibit consistent regularity and trends with input parameters, and characteristic features have been empirically described in terms of these parameters. Although much insight has been gained from experimental and model curves, no equation has been derived that can describe this behavior explicitly as functions of the key oxidation parameters.

Smialek, James L.↗

Dynamic Model Investigation of the Rough-Water Landing Characteristics of a Spacecraft

The investigation was made to study the rough-water landing characteristics of a Gemini type of spacecraft. The investigations were made with a 1/6-scale dynamic model in a simulated sea state 4 rough water. Parachute letdown landings were simulated with the model at various yaw angles and horizontal velocities. The vertical velocity and landing attitude remained constant. The range of maximum lateral and longitudinal acceleration was from about 3-1/2g to 16g while that for the maximum normal acceleration was from lg to 15g. The range of maximum angular acceleration was from about 0 to 190 radians per second(exp 2). The smoothest behavior and the lowest angular acceleration occurred at the 90 degree yaw angle. The normal acceleration was near minimum at this condition.

[Division] Imp. Str↗

Effects of Cavities and Protuberances on Transition over Hypersonic Vehicles

Surface protuberances and cavities on a hypersonic vehicle are known to cause several aerodynamic or aerothermodynamic issues. Most important of all, premature transition due to these surface irregularities can lead to a significant rise in surface heating. To help understand laminar-turbulent transition induced by protuberances or cavities on a Crew Exploration Vehicle (CEV) surface, high-fidelity numerical simulations are carried out for both types of trips on a CEV wind tunnel model. Due to the large bluntness, these surface irregularities reside in an accelerating subsonic boundary layer. For the Mach 6 wind tunnel conditions with a roughness Reynolds number Re(sub kk) of 800, it was found that a protuberance with a height to boundary layer thickness ratio of 0.73 leads to strong wake instability and spontaneous vortex shedding, while a cavity with identical geometry only causes a rather weak flow unsteadiness. The same cavity with a larger Reynolds number also leads to similar spontaneous vortex shedding and wake instability. The wake development and the formation of hairpin vortices for both protuberance and cavity were found to be qualitatively similar to that observed for an isolated hemisphere submerged in a subsonic, low speed flat-plate boundary layer. However, the shed vortices and their accompanying instability waves were found to be slightly stabilized downstream by the accelerating boundary layer along the CEV surface. Despite this stabilizing influence, it was found that the wake instability spreads substantially in both wall-normal and azimuthal directions as the flow is evolving towards a transitional state. Similarities and differences between the wake instability behind a protuberance and a cavity are investigated. Computations for the Mach 6 boundary layer over a slender cylindrical roughness element with a height to the boundary layer thickness of about 1.1 also shows spontaneous vortex shedding and strong wake instability. Comparisons of detailed flow structures associated with protuberances at subsonic and supersonic edge Mach numbers indicate distinctively different instability mechanisms.

Chang, Chau-Lyan↗

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, Joan↗

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. [1995] that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, M. Joan↗

Development of Multi-Physics Dynamics Models for High-Frequency Large-Amplitude Structural Response Simulation

An analytic approach is demonstrated to reveal potential pyroshock-driven dynamic effects causing power losses in the Thermo-Electric (TE) module bars of the Mars Science Laboratory (MSL) Multi-Mission Radioisotope Thermoelectric Generator (MMRTG). This study utilizes high-fidelity finite element analysis with SIERRA/PRESTO codes to estimate wave propagation effects due to large-amplitude suddenly-applied pyro shock loads in the MMRTG. A high fidelity model of the TE module bar was created with approximately 30 million degrees-of-freedom (DOF). First, a quasi-static preload was applied on top of the TE module bar, then transient tri-axial acceleration inputs were simultaneously applied on the preloaded module. The applied input acceleration signals were measured during MMRTG shock qualification tests performed at the Jet Propulsion Laboratory. An explicit finite element solver in the SIERRA/PRESTO computational environment, along with a 3000 processor parallel super -computing framework at NASA AMES, was used for the simulation. The simulation results were investigated both qualitatively and quantitatively. The predicted shock wave propagation results provide detailed structural responses throughout the TE module bar, and key insights into the dynamic response (i.e., loads, displacements, accelerations) of critical internal spring/piston compression systems, TE materials, and internal component interfaces in the MMRTG TE module bar. They also provide confidence on the viability of this high-fidelity modeling scheme to accurately predict shock wave propagation patterns within complex structures. This analytic approach is envisioned for modeling shock sensitive hardware susceptible to intense shock environments positioned near shock separation devices in modern space vehicles and systems.

Shock Waves↗

Iterative Nonlinear Tikhonov Algorithm with Constraints for Electromagnetic Tomography

Low frequency electromagnetic tomography such as the capacitance tomography (ECT) has been proposed for monitoring and mass-gauging of gas-liquid two-phase system under microgravity condition in NASA's future long-term space missions. Due to the ill-posed inverse problem of ECT, images reconstructed using conventional linear algorithms often suffer from limitations such as low resolution and blurred edges. Hence, new efficient high resolution nonlinear imaging algorithms are needed for accurate two-phase imaging. The proposed Iterative Nonlinear Tikhonov Regularized Algorithm with Constraints (INTAC) is based on an efficient finite element method (FEM) forward model of quasi-static electromagnetic problem. It iteratively minimizes the discrepancy between FEM simulated and actual measured capacitances by adjusting the reconstructed image using the Tikhonov regularized method. More importantly, it enforces the known permittivity of two phases to the unknown pixels which exceed the reasonable range of permittivity in each iteration. This strategy does not only stabilize the converging process, but also produces sharper images. Simulations show that resolution improvement of over 2 times can be achieved by INTAC with respect to conventional approaches. Strategies to further improve spatial imaging resolution are suggested, as well as techniques to accelerate nonlinear forward model and thus increase the temporal resolution.

Xu, Feng↗

A simple propulsion system model for the simulation of nonlinear dynamic thrust response

A simple mathematical jet engine description is presented where the measured transition functions of the engine thrust can be simulated in a quasi-stationary operation as well as in acceleration or deceleration schedules. Because of its simplicity and high fidelity, it is especially suited for representing jet engines in digital simulation programs.

Schaenzer, G.↗

Analysis of the Rim Inertial Measuring System (RIMS)

Nonlinear equations of motion are derived for the rim inertial measuring system, which is mounted in a strapped down configuration on a carrier vehicle. The measurement and equations derived from angular rates and linear accelerations were interfaced with the dynamic model of a user defined carried vehicle. Two methods are presented for rate measurement, and the results of the nonlinear computer simulation are presented.

Joshi, S. M.↗