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 145 records · Page 8

Structural dynamic model obtained from flight use with piloted simulation and handling qualities analysis

The ability to use flight data to determine an aircraft model with structural dynamic effects suitable for piloted simulation. and handling qualities analysis has been developed. This technique was demonstrated using SR-71 flight test data. For the SR-71 aircraft, the most significant structural response is the longitudinal first-bending mode. This mode was modeled as a second-order system, and the other higher order modes were modeled as a time delay. The distribution of the modal response at various fuselage locations was developed using a uniform beam solution, which can be calibrated using flight data. This approach was compared to the mode shape obtained from the ground vibration test, and the general form of the uniform beam solution was found to be a good representation of the mode shape in the areas of interest. To calibrate the solution, pitch-rate and normal-acceleration instrumentation is required for at least two locations. With the resulting structural model incorporated into the simulation, a good representation of the flight characteristics was provided for handling qualities analysis and piloted simulation.

Powers, Bruce G.↗

Modeling human vestibular responses during eccentric rotation and off vertical axis rotation

A mathematical model has been developed to help explain human multi-sensory interactions. The most important constituent of the model is the hypothesis that the nervous system incorporates knowledge of sensory dynamics into an "internal model" of these dynamics. This internal model allows the nervous system to integrate the sensory information from many different sensors into a coherent estimate of self-motion. The essence of the model is unchanged from a previously published model of monkey eye movement responses; only a few variables have been adjusted to yield the prediction of human responses. During eccentric rotation, the model predicts that the axis of eye rotation shifts slightly toward alignment with gravito-inertial force. The model also predicts that the time course of the perception of tilt following the acceleration phase of eccentric rotation is much slower than that during deceleration. During off vertical axis rotation (OVAR) the model predicts a small horizontal bias along with small horizontal, vertical, and torsional oscillations. Following OVAR stimulation, when stopped right- or left-side down, a small vertical component is predicted that decays with the horizontal post-rotatory response. All of the predictions are consistent with measurements of human responses.

NASA Discipline Neuroscience↗

On the Onset Delays of Solar Energetic Electrons and Protons: Evidence for a Common Accelerator

The processes responsible for the acceleration of solar energetic particles (SEPs) are still not well understood, including whether SEP electrons and protons are accelerated by common or separate processes. Using a numerical particle transport model that includes both pitch-angle and perpendicular spatial diffusion, we simulate, among other quantities, the onset delay for MeV electrons and protons and compare the results to observations of SEPs from widely separated spacecraft. Such observations have previously been interpreted, in a simple scenario assuming no perpendicular diffusion, as evidence for different electron and proton sources. We show that, by assuming a common particle source together with perpendicular diffusion, we are able to simultaneously reproduce the onset delays for both electrons and protons. We argue that this points toward a common accelerator for these particles. Moreover, a relatively broad particle source is required in the model to correctly describe the observations. This is suggestive of diffusive shock acceleration occurring at large shock structures playing a significant role in the acceleration of these SEPs.

R D Strauss↗

Assessment of Artemis-1 Pogo Flight Instrumentation System

As the space industry continues to strive for more efficient launch vehicles, it must relyon increasingly accurate predictive models. Verification of models typically requires physical testing. Flight data measurements offer the most real and therefore the most accurate data for model correlation. As NASA prepares for the inaugural launch of Space Launch System(SLS), Artemis-1, they must rely heavily on predictive system models to ensure flight safety.NASA has implemented a Development Flight Instrumentation (DFI) system in hopes of recovering useful flight data to aid in model correlation. Historically, some of the most important flight measurements are those that monitor the potentially destructive dynamic interaction of the structural and propellant modes – a phenomenon known as pogo. Pogo is a dynamic instability that can occur on a launch vehicle during any phase of ascent. During this investigation an end-to-end assessment of the Artemis-I Pogo-related DFI was performed to identify any obstacles inherent in the current instrumentation system which may prevent successful measurement of the data necessary to validate the current predictive fluid-structural models. Input drive signals were approximated and applied to a state-space system model to derive predicted pressure and acceleration responses. These predicted responses were then fed through a simulation of the data acquisition process in order to recover predicted measurements. Finally, a mock Flight Data Analysis (FDA) was performed to assess the ability of these measurements to meet the Flight Test Objectives(FTO).

Adam Johnson↗

Assessment of Artemis-1 Pogo Flight Instrumentation System

As the space industry continues to strive for more efficient launch vehicles, it must rely on increasingly accurate predictive models. Verification of models typically requires physical testing. Flight data measurements offer the most real and therefore the most accurate data for model correlation. As NASA prepares for the inaugural launch of Space Launch System (SLS), Artemis-1, they must rely heavily on predictive system models to ensure flight safety. NASA has implemented a Development Flight Instrumentation (DFI) system in hopes of recovering useful flight data to aid in model correlation. Historically, some of the most important flight measurements are those that monitor the potentially destructive dynamic interaction of the structural and propellant modes – a phenomenon known as pogo. Pogo is a dynamic instability that can occur on a launch vehicle during any phase of ascent. During this investigation an end-to-end assessment of the Artemis-I Pogo-related DFI was performed to identify any obstacles inherent in the current instrumentation system which may prevent successful measurement of the data necessary to validate the current predictive models of the Main Propulsion System (MPS). Input drive signals were approximated and applied to a system-level SLS state-space model to derive predicted pressure and acceleration responses. These predicted responses were then fed through a simulation of the data acquisition process in order to recover predicted measurements. Finally, a mock Flight Data Analysis (FDA) was performed to assess the ability of these measurements to meet the Flight Test Objectives (FTO).

Adam Johnson↗

Preferential Acceleration of Heavy Ions in a Spontaneously Fragmenting Flare Current Sheet

We study the ion acceleration in a mesoscale, spontaneously fragmenting flare current sheet (SFCS) characterized by the presence of a plasmoid cascade. The main subject of our investigation is to determine whether and how plasmoid cascades at intermediate scales in a fragmented current sheet of a solar flare can impact the (preferential) acceleration of specific ions. The time evolution of the SFCS is obtained from high-resolution 2.5D MHD simulations. The ion trajectories (in the background fields resulting from the MHD model), energies, and pitch angles are calculated using a relativistic test-particle code based on the half-acceleration–rotation–half-acceleration method. For light ions, the main acceleration effects of electromagnetic fields within the SFCS are analyzed using the guiding center approximation. We identify regions with the most-efficient ion acceleration within the SFCS, the accelerator efficiency, and spectra of the accelerated ions. The influence of the charge-to-mass ratio on ion behavior is also studied and resulting ion abundances are compared with observational data. The main ion acceleration takes place in the regions with a strong polarization term, which is part of the first-order Fermi acceleration. Because the term is mass dependent, heavier ions undergo preferential acceleration. The ion energy spectra, abundance-enhancement factors, and differential fluxes, obtained from the model, exhibit power-law profiles, in agreement with observed solar energetic particle events. Nonetheless, the obtained slopes for the abundance-enhancement factor do not exactly match the observed data. The computed slopes and profiles are not sensitive to changes in the initial plasma temperature.

David Kramoliš↗

Frostwing Co-Operation in Aircraft Icing Research

e aerodynamic effects of Cold Soaked Fuel Frost have become increasingly significant as airworthiness authorities have been asked to allow it during aircraft take-off. The Federal Aviation Administration and the Finnish Transport Safety Agency signed a Research Agreement in aircraft icing research in 2015 and started a research co-operation in frost formation studies, computational fluid dynamics for ground de/anti-icing fluids, and de/anti-icing fluids aerodynamic characteristics. The main effort has been so far on the formation and aerodynamic effects of CSFF. To investigate the effects, a generic high-lift common research wind tunnel model and DLR-F15 airfoil, representing the wing of a modern jet aircraft, was built including a wing tank cooling system. Real frost was generated on the wing in a wind tunnel test section and the frost thickness was measured with an Elcometer gauge. Frost surface geometry was measured with laser scanning and photogrammetry. The aerodynamic effect of the frost was studied in a simulated aircraft take-off sequence, in which the speed was accelerated to a typical rotation speed and the wing model was then rotated to an angle of attack used at initial climb. Time histories of the lift coefficient were measured with a force balance. The experiments showed that depending on the ambient temperature the frost may evaporate/melt during the take-off sequence. Lift losses after rotation with CSFF contamination at ambient temperatures of 4˚ to 7˚C above freezing point were measured to be 4 to 5 % for roughness values, k/c, below 10(exp -3). For comparison, lift loss tests with typical anti-icing fluids were to roughly equal lift losses. This paper gives an overview of the performed activities.

Koivisto, Pekka↗

Evaluation of Mid-Size Male Hybrid III Models for use in Spaceflight Occupant Protection Analysis

In an effort to improve occupant safety during dynamic phases of spaceflight, the National Aeronautics and Space Administration (NASA) has worked to develop occupant protection standards for future crewed spacecraft. One key aspect of these standards is the identification of injury mechanisms through anthropometric test devices (ATDs). Within this analysis, both physical and computational ATD evaluations are required to reasonably encompass the vast range of loading conditions any spaceflight crew may encounter. In this study the accuracy of publically available mid-size male HIII ATD finite element (FE) models are evaluated within applicable loading conditions against extensive sled testing performed on their physical counterparts. Methods: A series of sled tests were performed at the Wright Patterson Air force Base (WPAFB) employing variations of magnitude, duration, and impact direction to encompass the dynamic loading range for expected spaceflight. FE simulations were developed to the specifications of the test setup and driven using measured acceleration profiles. Both fast and detailed FE models of the mid-size male HIII were ran to quantify differences in their accuracy and thus assess the applicability of each within this field. Results: Preliminary results identify the dependence of model accuracy on loading direction, magnitude, and rate. Additionally the accuracy of individual response metrics are shown to vary across each model within evaluated test conditions. Causes for model inaccuracy are identified based on the observed relationships. Discussion: Computational modeling provides an essential component to ATD injury metric evaluation used to ensure the safety of future spaceflight occupants. The assessment of current ATD models lays the groundwork for how these models can be used appropriately in the future. Identification of limitations and possible paths for improvement aid in the development of these effective analysis tools.

Putnam, Jacob B.↗

Evaluation of Mid-Size Male Hybrid III Models for use in Spaceflight Occupant Protection Analysis

Introduction: In an effort to improve occupant safety during dynamic phases of spaceflight, the National Aeronautics and Space Administration (NASA) has worked to develop occupant protection standards for future crewed spacecraft. One key aspect of these standards is the identification of injury mechanisms through anthropometric test devices (ATDs). Within this analysis, both physical and computational ATD evaluations are required to reasonably encompass the vast range of loading conditions any spaceflight crew may encounter. In this study the accuracy of publically available mid-size male HIII ATD finite element (FE) models are evaluated within applicable loading conditions against extensive sled testing performed on their physical counterparts. Methods: A series of sled tests were performed at the Wright Patterson Air force Base (WPAFB) employing variations of magnitude, duration, and impact direction to encompass the dynamic loading range for expected spaceflight. FE simulations were developed to the specifications of the test setup and driven using measured acceleration profiles. Both fast and detailed FE models of the mid-size male HIII were ran to quantify differences in their accuracy and thus assess the applicability of each within this field. Results: Preliminary results identify the dependence of model accuracy on loading direction, magnitude, and rate. Additionally the accuracy of individual response metrics are shown to vary across each model within evaluated test conditions. Causes for model inaccuracy are identified based on the observed relationships. Discussion: Computational modeling provides an essential component to ATD injury metric evaluation used to ensure the safety of future spaceflight occupants. The assessment of current ATD models lays the groundwork for how these models can be used appropriately in the future. Identification of limitations and possible paths for improvement aid in the development of these effective analysis tools.

Putnam, J.↗

Generation and propagation of an electromagnetic pulse in the Trigger experiment and its possible role in electron acceleration

Instruments onboard the Trigger payload detected a large-amplitude, low-frequency, electric field pulse which was observed with a time delay consistent only with an electromagnetic wave. A model for this perturbation is constructed, and the associated field-aligned current is calculated as a function of altitude. This experiment may simulate the acceleration mechanism which results in the formation of auroral arcs, and possibly even other events in cosmic plasmas.

Kelley, M. C.↗

The Consequences of Alfven Waves and Parallel Potential Drops in the Auroral Zone

The goal of this research is to examine the causes of field-aligned plasma acceleration in the auroral zone using satellite data and numerical simulations. A primary question to be addressed is what causes the field-aligned acceleration of electrons (leading to precipitation) and ions (leading to upwelling ions) in the auroral zone. Data from the Fast Auroral SnapshoT (FAST) and Polar satellites is used when the two satellites are in approximate magnetic conjunction and are in the auroral region. FAST is at relatively low altitudes and samples plasma in the midst of the auroral acceleration region while Polar is at much higher altitudes and can measure plasmas and waves propagating towards the Earth. Polar can determine the sources of energy streaming earthward from the magnetotail, either in the form of field-aligned currents, electromagnetic waves or kinetic particle energy, that ultimately leads to the acceleration of plasma in the auroral zone. After identifying and examining several events, numerical simulations are run that bridges the spatial region between the two satellites. The code is a one-dimensional, long system length particle in cell simulation that has been developed to model the auroral region. A main goal of this research project is to include Alfven waves in the simulation to examine how these waves can accelerate plasma in the auroral zone.

Schriver, David↗

Adaptive sequential estimation with applications

A sequential estimator is presented and demonstrated which successfully tracks the system state and model errors in the presence of significant and unpredictable system or environmental variations. This adaptive estimation concept is shown to lead to a new and attractive approach to parameter identification problems. Simulation results are presented in an orbit determination problem, where the estimator tracks the orbit and unmodeled accelerations due to errors in the geopotential model. Some results in a re-entry trajectory estimation problem are also summarized.

Jazwinski, A. H.↗

What We Don't Understand About Ion Acceleration Flares

There are now strong associations between the (3)He-rich, Fe-rich ions in "impulsive" solar energetic particle (SEP) events and the similar abundances derived from gamma-ray lines from flares. Compact flares, where wave energy can predominate, are ideal sites for the study of wave-particle physics. Yet there are nagging questions about the magnetic geometry, the relation between ions that escape and those that interact, and the relative roles of cascading Alfven waves and the EMIC waves required to enhance He-3. There are also questions about the relative timing of ion and electron acceleration and of heating; these relate to the variation of ionization states before and during acceleration and during transport out of the corona. We can construct a model that addresses many of these issues, but problems do remain. Our greatest lack is realistic theoretical simulations of element abundances, spectra, and their variations. By contrast, we now have a much better idea of the acceleration at CME-driven shock waves in the rare but large "gradual" SEP events, largely because of their slow temporal evolution and great spatial extent.

Reames, Donald V.↗

Time-dependent simulation of oblique MHD cosmic-ray shocks using the two-fluid model

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic-ray (CR)-modified plane shock evolution. Most of the calculations are done with a two-fluid model for diffusive shock acceleration, but we provide also comparisons between a typical shock computed that way against calculations carried out using the more complete, momentum-dependent, diffusion-advection equation. We also illustrate a test showing that these simulations evolve to dynamical equilibria consistent with previously published steady state analytic calculations for such shocks. In order to improve understanding of the dynamical role of magnetic fields in shocks modified by CR pressure we have explored for time asymptotic states the parameter space of upstream fast mode Mach number, M(sub f), and plasma beta. We compile the results into maps of dynamical steady state CR acceleration efficiency, epsilon(sub c). We have run simulations using constant, and nonisotropic, obliquity (and hence spatially) dependent forms of the diffusion coefficient kappa. Comparison of the results shows that while the final steady states achieved are the same in each case, the history of CR-MHD shocks can be strongly modified by variations in kappa and, therefore, in the acceleration timescale. Also, the coupling of CR and MHD in low beta, oblique shocks substantially influences the transient density spike that forms in strongly CR-modified shocks. We find that inside the density spike a MHD slow mode wave can be generated that eventually steepens into a shock. A strong layer develops within the density spike, driven by MHD stresses. We conjecture that currents in the shear layer could, in nonplanar flows, results in enhanced particle accretion through drift acceleration.

Frank, Adam↗

Direct Numerical Simulations of Turbulent Convection: Variable Gravity and Uniform Rotation - Part 1

Turbulent convection may have played a major role in determining the structure and evolution of the primordial solar nebula, but current, incomplete models of convection and turbulence give very different results and remain largely untested in the absence of detailed astronomical observations. Numerical simulations provide an 'experimental' database for comparison with these models, and, to this end, direct numerical simulations of turbulent convection were performed with modifications intended to mimic some of the unique physical features of thin accretion disks, such as the primordial solar nebula: internal heating, a gravitational acceleration that is linearly proportional to the distance from midplane of the nebula, and rapid rotation. Peclet numbers in the simulations are comparable to those in solar nebula models; Rossby numbers in the simulations are an order of magnitude larger than those in solar nebula models because of the unrealistically high Prandtl and low, Reynolds numbers required to resolve all scales of the convective flow. We find that, despite the loss of buoyancy at midplane, turbulent motions easily penetrate the midplane region with little loss of intensity, providing efficient transport of heat and turbulent kinetic energy throughout the interior. A simple mixing length model modified to include rotation is found to give convective heat fluxes for the interior flow in rough agreement with the numerical simulations. We discuss these preliminary results with regard to assumptions about heating distributions and convective heat fluxes made in standard solar nebula models. More definitive comparisons with solar nebula modelling will become possible when more realistic effects of shear, density, stratification, and compression are included.

Cabot, W.↗

Multibody Based Digital Astronaut Dynamics Simulation

BACKGROUND: This study provides the Software, Robotics, & Simulation Division at the NASA Johnson Space Center with a verification tool for multibody dynamics simulation requiring human motion. The motivation stems from current studies of several Vibration Isolation & Stabilization(VIS)system designs that attenuate the moments and forces which would be transmitted to a spacecraft during an exercise. A multibody dynamics model for a proposed VIS was available previously[1], therefore modeling of the VIS was not needed for this work. The interest here is in creating the multibody dynamics model of an astronaut in motion which may be utilized independently or while attached to a mechanism. An existing simulation [2] that utilizes OpenSim [3,4] and an in-house multibody dynamics package (MBDyn) [5] is used in order to verify the astronaut model. The main advantage this model will have over the existing simulation is that everything will be processed in one tool. METHODS AND RESULTS: Creating the simulation required; estimation of Body Segment Inertial Parameters (BSIP),a multibody model of the human-VIS system, joint acceleration profiles, and input files for MBDyn, which is used for this analysis. The scaling factors provided by Dumas et al. [6] are utilized in estimating the BSIP. Anthropometric data are used for estimating these parameters, the Anthropometric Survey of US Army Personnel (ANSUR II) [7] was the source. The astronaut model consists of 15 bodies, 14 joints and 32 degrees of freedom, with the dynamics topology generated using the center of mass locations and anthropometric data. MBDyn has an option for prescribed joint motion (PJM), which requires joint acceleration data as input. The joint angle data is first obtained from a motion capture system and then processed through code that has been created to generate approximate joint acceleration profiles. The topology tree for the astronaut model begins at the right foot up to the pelvis where there is one branch for going down the left leg and another for the torso. The torso branch leads to branches for the arms and a leaf body for the head/neck segment. For attachment to the VIS, the heel of the right foot is connected to the VIS platform through a fixed joint, resembling a foot restraint. The left foot does not attach to the platform in order to prevent a system with a closed loop. Topology and symmetry of the astronaut model were verified through kinematic analysis. Further verification of the forces and moments transmitted to the VIS were verified against the existing simulation. There was a satisfactory level of agreement when testing a simple motion, for example, rocking back and forth. Full exercise motions are to be tested soon. The main outcome has been a novel application of MBDyn for biomechanics modeling that is now available for dynamic simulations involving human motion. The estimation of BSIP was another useful result of this study, requiring only 15 inputs for generating mass properties of a theoretical astronaut model. Expansion on this work is possible by going through an alternative route in obtaining the joint motion data. Instead of high-tech and often expensive motion capture systems, an individual may watch videos with high focus and at a slow motion for each individual segment in order to determine the initial and final time and angle for that specific degree of freedom. Synthetic trajectories may also be created if there is no video reference available.

F N Matari↗

Magnetosphere-ionosphere coupling during plasmoid evolution - First results

The influence of magnetosphere-ionosphere coupling on the dynamic evolution of the magnetotail is investigated by a three-dimensional resistive MHD code that includes the effects of the closure of field-aligned currents in a simple resistive model ionosphere. Particular emphasis is on the role of this coupling during substorm evolution and the modification of the latter by the convection driven by the ionospheric electric fields. For comparison, results are presented from a simulation which uses an infinitely conducting ionosphere but is otherwise identical. Comparison of the two simulations shows that the major impact of magnetosphere-ionospheric communication is an acceleration of magnetotail evolution. Otherwise, phenomena in the two models are qualitatively similar. It is concluded that ionospheric effects do not significantly affect substorm associated magnetotail dynamics.

Hesse, Michael↗

Graphics Processing Unit (GPU) Acceleration of the Goddard Earth Observing System Atmospheric Model

The Goddard Earth Observing System 5 (GEOS-5) is the atmospheric model used by the Global Modeling and Assimilation Office (GMAO) for a variety of applications, from long-term climate prediction at relatively coarse resolution, to data assimilation and numerical weather prediction, to very high-resolution cloud-resolving simulations. GEOS-5 is being ported to a graphics processing unit (GPU) cluster at the NASA Center for Climate Simulation (NCCS). By utilizing GPU co-processor technology, we expect to increase the throughput of GEOS-5 by at least an order of magnitude, and accelerate the process of scientific exploration across all scales of global modeling, including: The large-scale, high-end application of non-hydrostatic, global, cloud-resolving modeling at 10- to I-kilometer (km) global resolutions Intermediate-resolution seasonal climate and weather prediction at 50- to 25-km on small clusters of GPUs Long-range, coarse-resolution climate modeling, enabled on a small box of GPUs for the individual researcher After being ported to the GPU cluster, the primary physics components and the dynamical core of GEOS-5 have demonstrated a potential speedup of 15-40 times over conventional processor cores. Performance improvements of this magnitude reduce the required scalability of 1-km, global, cloud-resolving models from an unfathomable 6 million cores to an attainable 200,000 GPU-enabled cores.

Putnam, Williama↗