Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamics 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 1,333 records · Page 74

Numerical simulation of dynamics of brushless dc motors for aerospace and other applications. Volume 2: User's guide to computer EMA model

A description and user's guide of the computer program developed to simulate the dynamics of an electromechanical actuator for aerospace applications are presented. The effects of the stator phase currents on the permanent magnets of the rotor are examined. The voltage and current waveforms present in the power conditioner network during the motoring, regenerative braking, and plugging modes of operation are presented and discussed.

Demerdash, N. A. O.↗

Mars Tumbleweed Simulation Using Singular Perturbation Theory

The Mars Tumbleweed is a new surface rover concept that utilizes Martian winds as the primary source of mobility. Several designs have been proposed for the Mars Tumbleweed, all using aerodynamic drag to generate force for traveling about the surface. The Mars Tumbleweed, in its deployed configuration, must be large and lightweight to provide the ratio of drag force to rolling resistance necessary to initiate motion from the Martian surface. This paper discusses the dynamic simulation details of a candidate Tumbleweed design. The dynamic simulation model must properly evaluate and characterize the motion of the tumbleweed rover to support proper selection of system design parameters. Several factors, such as model flexibility, simulation run times, and model accuracy needed to be considered in modeling assumptions. The simulation was required to address the flexibility of the rover and its interaction with the ground, and properly evaluate its mobility. Proper assumptions needed to be made such that the simulated dynamic motion is accurate and realistic while not overly burdened by long simulation run times. This paper also shows results that provided reasonable correlation between the simulation and a drop/roll test of a tumbleweed prototype.

Raiszadeh, Behzad↗

Simulation of dynamic crystal plasticity with a Lagrangian discontinuous Galerkin hydrodynamic method

Here we present a new Lagrangian modal discontinuous Galerkin (DG) hydrodynamic method that supports a dynamic dislocation based crystal plasticity model for simulating the mechanical behavior of crystallographic materials, both single crystal and polycrystalline, under dynamic conditions. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated by Taylor series polynomials of varying degree. These polynomials describe macro-scale hydrodynamic behavior while their evolution is determined by evaluating the dynamic crystal plasticity model at material points within the element. The dynamic crystal plasticity model is sensitive to the time increment size, with too large of time increments leading to instability in the model. To mitigate this, the temporal evolution of the dynamic crystal plasticity model is achieved with the combination of a sub-incrementing scheme with Heun’s third-order time integration scheme, which is also used to temporally evolve the governing equations. The implementation of the dynamic crystal plasticity model within the DG framework is tested using a 2D approximation of the Taylor impact test with a single crystal material, using quadratic elements that have faces that can bend. In addition to the standard continuous material modeling, we propose a new simulation method that would represent the heterogeneous behavior of polycrystalline microstructures within an element by varying the position and material properties of the material points within the element. This method is demonstrated using random orientation distributions on materials points that are arranged in both structured and random configurations.

42 ENGINEERING↗

Ultrafast terahertz-field-driven ionic response in ferroelectric BaTiO 3

The dynamical processes associated with electric field manipulation of the polarization in a ferroelectric remain largely unknown but fundamentally determine the speed and functionality of ferroelectric materials and devices. Here in this paper we apply subpicosecond duration, single-cycle terahertz pulses as an ultrafast electric field bias to prototypical BaTiO 3 ferroelectric thin films with the atomic-scale response probed by femtosecond x-ray-scattering techniques. We show that electric fields applied perpendicular to the ferroelectric polarization drive large-amplitude displacements of the titanium atoms along the ferroelectric polarization axis, comparable to that of the built-in displacements associated with the intrinsic polarization and incoherent across unit cells. This effect is associated with a dynamic rotation of the ferroelectric polarization switching on and then off on picosecond time scales. These transient polarization modulations are followed by long-lived vibrational heating effects driven by resonant excitation of the ferroelectric soft mode, as reflected in changes in the c-axis tetragonality. The ultrafast structural characterization described here enables a direct comparison with first-principles-based molecular-dynamics simulations, with good agreement obtained.

36 MATERIALS SCIENCE↗

A data-driven multiscale model for reactive wetting simulations

Here, we describe a data-driven, multiscale technique to model reactive wetting of a silver–aluminum alloy on a Kovar™ (Fe-Ni-Co alloy) surface. We employ molecular dynamics simulations to elucidate the dependence of surface tension and wetting angle on the drop’s composition and temperature. A design of computational experiments is used to efficiently generate training data of surface tension and wetting angle from a limited number of molecular dynamics simulations. The simulation results are used to parameterize models of the material’s wetting properties and compute the uncertainty in the models due to limited data. The data-driven models are incorporated into an engineering-scale (continuum) model of a silver–aluminum sessile drop on a Kovar™ substrate. Model predictions of the wetting angle are compared with experiments of pure silver spreading on Kovar™ to quantify the model-form errors introduced by the limited training data versus the simplifications inherent in the molecular dynamics simulations. The paper presents innovations in the determination of “convergence” of noisy MD simulations before they are used to extract the wetting angle and surface tension, and the construction of their models which approximate physio-chemical processes that are left unresolved by the engineering-scale model. Together, these constitute a multiscale approach that integrates molecular-scale information into continuum scale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Validation of Fault Displacements from Dynamic Rupture Simulations against the Observations from the 1992 Landers Earthquake

ABSTRACT Coseismic fault displacements in large earthquakes have caused significant damage to structures and lifelines on and near fault lines. Coseismic displacements represent a real threat, especially to distributed infrastructure systems. For infrastructure systems that can not avoid active faults, engineering displacement demands are defined using probabilistic fault-displacement hazard analyses (PFDHA). However, PFDHA models are sparse and poorly constrained partly due to the scarcity of detailed fault-displacement observations. Advancements in dynamic rupture simulation methods make them an attractive approach to address this important issue. Because fault displacements can be simulated for various geologic conditions as constrained by current knowledge about earthquake processes, they can be used to supplement the observation datasets. In addition to providing on-fault displacements, when used with appropriate constitutive models for the bulk medium, they can capture off-fault distributed inelastic deformations as well. For viable extrapolation, simulations must first be validated against data. In this article, we summarize the calibration and validation of the dynamic rupture model against the observations of the well-documented 1992 Landers earthquake. We defined a preferred model that reproduces several first-order fault-displacement metrics such as the on-fault partition of the total displacement, the mean fault-zone width, and the location of the peak displacement. Simulated ground motions consistent with the observations ensure that all physics important to modeling have been properly parameterized. For the extrapolation, we generated a suite of dynamic rupture models to quantify expected fault-displacement metrics, their intercorrelations, and magnitude dependencies, which are in part supported by the Landers and other recent earthquakes. Our validation and extrapolation exercise paves the way for using dynamic rupture modeling to quantitatively address fault-displacement hazard on a broader scale. The results are promising and are expected to be useful to inform PFDHA model development.

Geochemistry & Geophysics↗

Numerical simulation and dynamical response of a moored hydrokinetic turbine operating in the wake of an upstream turbine for control design

Herein numerical simulation of a downstream hydrokinetic turbine operating in the wake of an upstream turbine for feedback control design is presented. Wake effects from an upstream turbine are quantified in terms of wake velocity and amplified turbulence levels. These effects are integrated in an in-stream hydrokinetic turbine numerical simulation that utilizes a Blade Element Momentum approach with a dynamic wake inflow model. Simulations are carried out on a fixed turbine model to simulate operation in river or tidal channels with conventional foundations, as well as on a compliantly moored turbine model such as those designed to operate in open ocean currents.

13 HYDRO ENERGY↗

PB-AM: An open-source, fully analytical linear poisson-boltzmann solver

Here, we present the open source distributed software package Poisson-Boltzmann Analytical Method (PB-AM), a fully analytical solution to the linearized PB equation, for molecules represented as non-overlapping spherical cavities. The PB-AM software package includes the generation of outputs files appropriate for visualization using visual molecular dynamics, a Brownian dynamics scheme that uses periodic boundary conditions to simulate dynamics, the ability to specify docking criteria, and offers two different kinetics schemes to evaluate biomolecular association rate constants. Given that PB-AM defines mutual polarization completely and accurately, it can be refactored as a many-body expansion to explore 2- and 3-body polarization. Additionally, the software has been integrated into the Adaptive Poisson-Boltzmann Solver (APBS) software package to make it more accessible to a larger group of scientists, educators, and students that are more familiar with the APBS framework.

97 MATHEMATICS AND COMPUTING↗

Estimates for the Aerodynamic Coefficients of Ringsail and Disk-Gap-Band Parachutes Operating on Mars

The creation of flight dynamics simulations is often part of the design, development, and flight operations of missions that include entry, descent, and landing (EDL) in their concept of operation. Such flight dynamics simulations include models for many of the EDL sub-systems. For EDL systems that include a parachute, models of the parachute’s aerodynamic characteristics must be created for inclusion in the flight dynamics simulations. Models are presented for the aerodynamic coefficients (i.e., C(D),C(Tot),C(T),C(N), and C(m,SLCP)) of Supersonic Ringsail (SSRS) and Disk-Gap-Band (DGB) parachutes as functions of total porosity, λ(T), Mach number, M, and total angle of attack, λ(T) (when necessary). The source aerodynamic coefficients data used for creating these models were obtained from a wind tunnel test of subscale parachutes. In this wind tunnel test, subscale parachutes of both parachute types were fabricated from two different fabrics with very different permeabilities. By varying the fabric permeability, while maintaining the parachute geometry constant, it was possible to vary λ(T). The fabric permeability test data necessary for the calculation of λ(T) were obtained from samples of the same fabrics used to fabricate the subscale parachutes. Although the models for the aerodynamic coefficients are simple polynomial functions of λ(T) and M, they are capable of producing good reproductions of the source data. The (λ(T), M) domains over which these models are applicable are clearly defined. The models and their domains are applicable to flight operations on Mars. The models presented in this report are not unique – other models could be created from the source data. To make the creation of alternative numerical models, the aerodynamic coefficients source data are included in this report.

Cruz, Juan R.↗

Developing a Reduced 240-Bus WECC Dynamic Model for Frequency Response Study of High Renewable Integration

The ongoing changes in the generation resource mix, driven by the rapid adoption of inverter-based resources (IBR) as well as the early retirement of synchronous generators, are bringing new challenges to the planning and operation of bulk electric power systems. Consequently, there is an increasing need to understand, design, and quantify the reliability service provision from IBRs by performing integrated scheduling and dynamic simulations. However, test systems that have consistent scheduling and dynamic models rarely exist largely because of the decoupled nature of the two simulations on a synchronous generator-dominated system. This paper develops the dynamic model of a reduced Western Electricity Coordinating Council (WECC) system. In conjunction with the existing scheduling model, it is suitable for integrating scheduling and dynamic simulations. The reduced 240-bus WECC model reflects the generation resource mix of the Western Interconnection as of 2018. Moreover, the developed dynamic model is validated against field frequency events measured by FNET/GridEye and preserves the dominant inter-area oscillation mode in WECC.

dynamic simulation↗

Validation of Deterministic Broadband Ground Motion and Variability from Dynamic Rupture Simulations of Buried Thrust Earthquakes

According to the article, we numerically model broadband ground motion (up to 5–7.5 Hz⁠) from blind‐thrust scenario earthquakes matching the fault geometry of the 1994 Mw 6.7 Northridge earthquake. Several realizations are modeled (by varying the hypocenter location in the dynamic rupture simulation) in a 1D‐layered velocity profile. In addition, we include Q(f )⁠, nonlinear effects from Drucker–Prager plasticity, and superimpose small‐scale medium complexity in both a 1D‐layered and 3D velocity model within the subsequent wave propagation. We investigate characteristics of the ground motion and its variability up to 50 km from the fault by comparing them with ground‐motion prediction equations (GMPEs), simple proxy metrics, as well as strong ground motion records from the Northridge event. We find that median ground motion closely follows the trend predicted by GMPEs and that the intraevent standard deviation, although varying with hypocenter location, lies near that of GMPE models. Plasticity affects ground‐motion amplitudes in regions near the source, reducing intraevent variability above ∼0.5 Hz⁠. Heterogeneity in the velocity structure on both the regional and small scales is needed for the simulated data to match two proxy metrics: the period‐to‐period correlation of spectral acceleration (SA) and the ratio of maximum‐to‐median SA. Although small‐scale heterogeneity has a negligible effect on median SA for this style of rupture, it serves to significantly increase the cumulative absolute velocity, better agreeing with observations. When compared with strong‐motion data, we find that long‐wavelength velocity structure within our deterministic simulations reduces bias at both short and long periods. Finally, synthetic ground motion at both footwall and hanging‐wall sites has no clear dependence on the distance to rupture (at both short and long periods); directivity is likely overpowering any hanging‐wall effect.

Geochemistry & Geophysics↗

Insight into the Li{sub 2}CO{sub 3}–K{sub 2}CO{sub 3} eutectic mixture from classical molecular dynamics: Thermodynamics, structure, and dynamics

We use molecular dynamics simulations to study the thermodynamics, structure, and dynamics of the Li{sub 2}CO{sub 3}–K{sub 2}CO{sub 3} (62:38 mol. %) eutectic mixture. We present a new classical non-polarizable force field for this molten salt mixture, optimized using experimental and first principles molecular dynamics simulations data as reference. This simple force field allows efficient molecular simulations of phenomena at long time scales. We use this optimized force field to describe the behavior of the eutectic mixture in the 900–1100 K temperature range, at pressures between 0 and 5 GPa. After studying the equation of state in these thermodynamic conditions, we present molecular insight into the structure and dynamics of the melt. In particular, we present an analysis of the temperature and pressure dependence of the eutectic mixture’s self-diffusion coefficients, viscosity, and ionic conductivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Imogolite in water: Simulating the effects of nanotube curvature on structure and dynamics

Imogolite is a fascinating inorganic nanotube that is found in nature or synthesized in a laboratory. The synthesis process is carried out in liquid media, and leads to the formation of almost monodisperse diameter nanotubes. Here we investigate, employing classical molecular dynamics simulations, the interaction of water and imogolite for nanotubes of several radii. We established that water penetrates the pores of N = 9 and larger nanotubes, and adopts a coaxial arrangement in it. Also, while water molecules can diffuse along the center of the nanotube, the molecules next to the inner imogolite walls have very low mobility. At the outer nanotube wall, an increase of water density is observed, this effect extends up to 1 nm, beyond which water properties are bulk-like. Both phenomena are affected by the imogolite curvature.

58 GEOSCIENCES↗

Man-in-the-control-loop simulation of manipulators

A method to achieve man-in-the-control-loop simulation is presented. Emerging real-time dynamics simulation suggests a potential for creating an interactive design workstation with a human operator in the control loop. The recursive formulation for multibody dynamics simulation is studied to determine requirements for man-in-the-control-loop simulation. High speed computer graphics techniques provides realistic visual cues for the simulator. Backhoe and robot arm simulations are implemented to demonstrate the capability of man-in-the-control-loop simulation.

Chang, J. L.↗

A new environment to simulate the dynamics in the close proximity of rubble-pile asteroids

This paper presents a new environment to simulate close-proximity dynamics around rubble-pile asteroids. The code provides methods for modeling the asteroid’s gravity field and surface through granular dynamics. It implements stateof-the-art techniques to model both gravity and contact interaction between particles: 1) mutual gravity as either direct N2 or Barnes-Hut GPU-parallel octree and 2) contact dynamics with a soft-body (force-based, smooth dynamics), hard-body (constraint-based, non-smooth dynamics), or hybrid (constraint-based with compliance and damping) approach. A very relevant feature of the code is its ability to handle complex-shaped rigid bodies and their full 6D motion. Examples of spacecraft close-proximity scenarios and their numerical simulations are shown.

Ferrari, Fabio↗

The effect of the solenoidal condition on the numerical magnetohydrodynamic simulation of coronal dynamics

Several MHD simulations of coronal dynamics performed more than 10 years ago violated the solenoidality condition for an 'open' magnetic field topology. Using an improved code, consideration is given to the effect on the physical validity of the numerical simulation (for a representative pulse disturbance) for the case when solenoidality is deliberately violated as compared to the case when it is preserved. It is found that the error incurred in this specific case ('open' topology) in the energy density and in the plasma density profiles is rather small and, hence, does not invalidate the earlier conclusions concerning mass and wave motion.

Wu, S. T.↗

Simulation of Dynamics of a Flexible Miniature Airplane

A short report discusses selected aspects of the development of the University of Florida micro-aerial vehicle (UFMAV) basically, a miniature airplane that has a flexible wing and is representative of a new class of airplanes that would operate autonomously or under remote control and be used for surveillance and/or scientific observation. The flexibility of the wing is to be optimized such that passive deformation of the wing in the presence of aerodynamic disturbances would reduce the overall response of the airplane to disturbances, thereby rendering the airplane more stable as an observation platform. The aspect of the development emphasized in the report is that of computational simulation of dynamics of the UFMAV in flight, for the purpose of generating mathematical models for use in designing control systems for the airplane. The simulations are performed by use of data from a wind-tunnel test of the airplane in combination with commercial software, in which are codified a standard set of equations of motion of an airplane, and a set of mathematical routines to compute trim conditions and extract linear state space models.

Waszak, Martin R.↗

Computational Fluid Dynamics (CFD) Simulation of Hypersonic Turbine-Based Combined-Cycle (TBCC) Inlet Mode Transition

Methods of computational fluid dynamics were applied to simulate the aerodynamics within the turbine flowpath of a turbine-based combined-cycle propulsion system during inlet mode transition at Mach 4. Inlet mode transition involved the rotation of a splitter cowl to close the turbine flowpath to allow the full operation of a parallel dual-mode ramjet/scramjet flowpath. Steady-state simulations were performed at splitter cowl positions of 0deg, -2deg, -4deg, and -5.7deg, at which the turbine flowpath was closed half way. The simulations satisfied one objective of providing a greater understanding of the flow during inlet mode transition. Comparisons of the simulation results with wind-tunnel test data addressed another objective of assessing the applicability of the simulation methods for simulating inlet mode transition. The simulations showed that inlet mode transition could occur in a stable manner and that accurate modeling of the interactions among the shock waves, boundary layers, and porous bleed regions was critical for evaluating the inlet static and total pressures, bleed flow rates, and bleed plenum pressures. The simulations compared well with some of the wind-tunnel data, but uncertainties in both the windtunnel data and simulations prevented a formal evaluation of the accuracy of the simulation methods.

Slater, John W.↗