Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Characterizing Porous and Nonporous Phenolic Resins from Molecular Dynamics Simulations

Phenolic resins are an important component of many ablative heat shield materials, which protect spacecrafts from the extreme temperatures reached during atmospheric entry. Examples include the high-density Heritage Carbon Phenolic (HCP) used in the Pioneer-Venus and Galileo missions, as well as the low-density Phenolic Impregnated Carbon Ablator (PICA) used in the Mars Science Laboratory and Mars 2020 missions. Additionally, recent developments within NASA have produced the mid-density Heatshield for Extreme Entry Environment Technology (HEEET) and its derivative 3D Woven Mid-Density Carbon Phenolic (3MDCP). Unlike the nonporous phenolic in HCP, PICA and HEEET/3MDCP are fabricated by infusing preforms with diluted phenolic formulations to obtain a lower density porous matrix. Despite the importance of the phenolic phase to the material response during entry, the variation in properties of porous and nonporous phenolic is not well understood. Here, we present an investigation of porous and nonporous phenolic resins using molecular dynamics (MD) simulations. Resin cure is mimicked in the simulations through the inclusion of representative reaction templates to generate accurate models of the complex crosslinked structures. To create porous models, explicit solvent molecules are included during the cure simulations. We observe nanoscale separation of the phenolic and solvent phases, which results in significant differences in the final structures of porous and nonporous models. In addition to a quantitative assessment of the network structure and porosity, we elucidate the effects of the phenolic formulation on the final material properties. These results are compared with experimental data as appropriate.

phenolic↗
NASA NTRS1993

Computational strategies in the dynamic simulation of constrained flexible MBS

This research focuses on the computational dynamics of flexible constrained multibody systems. At first a recursive mapping formulation of the kinematical expressions in a minimum dimension as well as the matrix representation of the equations of motion are presented. The method employs Kane's equation, FEM, and concepts of continuum mechanics. The generalized active forces are extended to include the effects of high temperature conditions, such as creep, thermal stress, and elastic-plastic deformation. The time variant constraint relations for rolling/contact conditions between two flexible bodies are also studied. The constraints for validation of MBS simulation of gear meshing contact using a modified Timoshenko beam theory are also presented. The last part deals with minimization of vibration/deformation of the elastic beam in multibody systems making use of time variant boundary conditions. The above methodologies and computational procedures developed are being implemented in a program called DYAMUS.

Amirouche, F. M. L.↗
NASA NTRS1998

Test of a flexible spacecraft dynamics simulator

There are a number of approaches one can take to modeling the dynamics of a flexible body. While one can attempt to capture the full dynamical behavior subject to disturbances from actuators and environmental torques, such a detailed description often is unnecessary. Simplification is possible either by limiting the amplitude of motion to permit linearization of the dynamics equations or by restricting the types of allowed motion. In this work, we study the nonlinear dynamics of bending deformations of wire booms on spinning spacecraft. The theory allows for large amplitude excursions from equilibrium while enforcing constraints on the dynamics to prohibit those modes that are physically less relevant or are expected to damp out fast. These constraints explicitly remove the acoustic modes (i.e., longitudinal sound waves and shear waves) while allowing for arbitrary bending and twisting, motions which typically are of lower frequency. As a test case, a spin axis reorientation maneuver by the Polar Plasma Laboratory (POLAR) spacecraft has been simulated. POLAR was chosen as a representative spacecraft because it has flexible wire antennas that extend to a length of 65 meters. Bending deformations in these antennas could be quite large and have a significant effect on the attitude dynamics of the spacecraft body. Summary results from the simulation are presented along, with a comparison with POLAR flight data.

Dichmann, Donald↗
NASA NTRS1989

Development of efficient computer program for dynamic simulation of telerobotic manipulation

Research in robot control has generated interest in computationally efficient forms of dynamic equations for multi-body systems. For a simply connected open-loop linkage, dynamic equations arranged in recursive form were found to be particularly efficient. A general computer program capable of simulating an open-loop manipulator with arbitrary number of links has been developed based on an efficient recursive form of Kane's dynamic equations. Also included in the program is some of the important dynamics of the joint drive system, i.e., the rotational effect of the motor rotors. Further efficiency is achieved by the use of symbolic manipulation program to generate the FORTRAN simulation program tailored for a specific manipulator based on the parameter values given. The formulations and the validation of the program are described, and some results are shown.

Chen, J.↗
NASA NTRS1988

Obstacles to high fidelity multibody dynamics simulation

The authors take a wide ranging look at issues involved in the dynamic modeling of complex, multibodied orbiting space systems. The nature of the multibody problem is reviewed and capabilities and limitations of two major codes (DISCOS, TREETOPS) are assessed. Problem areas that limit accuracy or contribute to uncertainty in the modeling, simulation process are identified. As a specific case, the significance of including, or of not including, second-degree geometric nonlinearity in the elastic displacement field (foreshortening) is evaluated. Conclusions and recommendations are made concerning the direction future development should take to achieve higher-fidelity and more computationally efficient multibody software solutions.

Lips, K. W.↗
NASA NTRS1974

The effect of simulator dynamics on pilot response

The effects of visual display dynamics on the altitude tracking performance of a subject in a fixed base flight simulator are considered. The subject, flying the linearized longitudinal equations of motion, attempts to maintain the same altitude as two airplanes positioned three hundred feet ahead, as in level formation flying. The horizon together with the two leading aircraft are represented symbolically on a CRT display. The subject's aircraft is disturbed by atmospheric turbulence. The data indicate a relationship between the bandwidth of the display dynamics and the short period characteristics of the simulated airplane. For an airplane with a relatively fast pitch response the presence of altitude display dynamics, with a bandwidth as high as five times the short period natural frequency, causes significant degradation of altitude tracking performance.

Weener, E. F.↗
NASA NTRSDate not supplied

Validation of Cryogenic Propellant Tank Filling using Computational Fluid Dynamics Simulation

Validation of cryogenic propellant tank filling was performed using the Computational Fluid Dynamics (CFD) solver Loci/STREAM-VoF. The validation effort helped identify modeling methodologies that enable NASA to best support its partners in both launch pad and on-orbit filling operations. Data from liquid hydrogen ground tests filled via jet injection were used for the validation effort which include a sensitivity to initial tank wall temperature. Initially hot walls are expected to yield rapid evaporation and possibly boiling. A Volume of Fluid (VoF) methodology was used to capture the gas-liquid interface. Rapid breakup of the liquid jet was observed in simulation results as liquid evaporated and expanded. Inflowing liquid transitioned to a contiguous jet as tank temperatures decreased, gas pressure increased, and saturation conditions at the incoming liquid temperature were approached. The final phase of filling was distinguished by rapid gas pressure rise due to a higher rate of gas volume compression than condensation at the liquid surface. Key physics of propellant tank filling were captured in the computational predictions, and opportunities for added simulation robustness and efficiency in future modeling efforts were identified.

CFM↗
NASA NTRS1997

Molecular Dynamics Simulation of a Multi-Walled Carbon Nanotube Based Gear

We used molecular dynamics to investigate the properties of a multi-walled carbon nanotube based gear. Previous work computationally suggested that molecular gears fashioned from (14,0) single-walled carbon nanotubes operate well at 50-100 gigahertz. The gears were formed from nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. The gear in this study was based on the smallest multi-walled nanotube supported by some experimental evidence. Each gear was a (52,0) nanotube surrounding a (37,10) nanotube with approximate 20.4 and 16,8 A radii respectively. These sizes were chosen to be consistent with inter-tube spacing observed by and were slightly larger than graphite inter-layer spacings. The benzyne teeth were attached via 2+4 cycloaddition to exterior of the (52,0) tube. 2+4 bonds were used rather than the 2+2 bonds observed by Hoke since 2+4 bonds are preferred by naphthalene and quantum calculations by Jaffe suggest that 2+4 bonds are preferred on carbon nanotubes of sufficient diameter. One gear was 'powered' by forcing the atoms near the end of the outside buckytube to rotate to simulate a motor. A second gear was allowed to rotate by keeping the atoms near the end of its outside buckytube on a cylinder. The ends of both gears were constrained to stay in an approximately constant position relative to each other, simulating a casing, to insure that the gear teeth meshed. The stiff meshing aromatic gear teeth transferred angular momentum from the powered gear to the driven gear. The simulation was performed in a vacuum and with a software thermostat. Preliminary results suggest that the powered gear had trouble turning the driven gear without slip. The larger radius and greater mass of these gears relative to the (14,0) gears previously studied requires a smaller rotation rate and multiple rows of teeth to avoid excessive force on the gear teeth resulting, in slip and failure of the driven gear to turn. We hope that studies such as these will eventually lead to synthesis of components that can be assembled into atomically precise fullerene machines. These machines, in turn, may someday be used in machine-phase fullerene materials with remarkable properties.

Han, Jie↗
NASA NTRS1990

Staggered solution procedures for multibody dynamics simulation

A computational procedure is presented for the direct integration of the multibody dynamical (MBD) equations with constraints. Numerical experiments conducted with this MBD procedure indicate that it yields robust solutions in cases where the step size furnishes more than 20 samples for the period of highest apparent response frequency in a given multibody system. The MBD solution procedure is implemented in two separate modules: the generalized coordinated solver, 'CINT', and the constraint Lagrange multiplier solver, 'LINT'.

Park, K. C.↗
DOE OSTI2025

First Principles Molecular Dynamics Simulations of Ammonia Adsorption onto MFI Zeolite Nanosheets

MFI zeolite nanosheets membranes are promising candidates for ammonia separation from nitrogen and hydrogen, yet questions remain on the origin of their high selectivity. Silanols, Si-OH, are present in high concentration at the surface of zeolite nanosheets, and force-field-based simulations indicate that surface adsorption at the silanols contributes to selectivity. Silanols can chemically react with ammonia, which may further contribute to the ability of zeolite nanosheet membranes to separate it from other gases. In this work, we use first-principles molecular dynamics techniques to simulate ammonia’s behavior within stacked MFI zeolite nanosheets. We find that at 523 K and a loading corresponding to 35 bar, conditions desired for the ammonia separation, about 30% of ammonia reacts with surface silanols. Our work explores H-bonding and proton transfer within this system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First-Principles Molecular Dynamics Simulations of Ammonia Adsorption onto MFI Zeolite Nanosheets

MFI zeolite nanosheet membranes are promising candidates for ammonia separation from nitrogen and hydrogen gases, yet questions remain on the origin of their high selectivity. Silanols, Si–OH, are present in high concentration at the surface of zeolite nanosheets, and force-field-based simulations indicate that surface adsorption at the silanols contributes to selectivity. Acidic silanol groups can chemically react with ammonia via transfer of a proton to form ammonium ions, which may further contribute to the ability of zeolite nanosheet membranes to separate ammonia from other gases. In this work, we used first-principles molecular dynamics techniques to simulate ammonia’s behavior within stacked MFI zeolite nanosheets. We found that at 523 K and a loading corresponding to 35 bar, conditions desired for the ammonia separation, about 30% of ammonia reacts with surface silanols. Our work explored H-bonding and proton transfer within this system.

ammonia↗
NASA NTRSDate not supplied

Validation of Cryogenic Propellant Tank Self-Pressurization by Leveraging Reduced Order Modeling within Computational Fluid Dynamics Simulation

Validation of cryogenic propellant tank self-pressurization was performed using a hybrid Computational Fluid Dynamics (CFD) and reduced order modeling methodology. Data from a liquid hydrogen ground test conducted at the K-site facility at the National Aeronautics and Space Administration (NASA) Glenn Research Center was used for the validation effort. Liquid phase dynamics were explicitly resolved with a CFD tool. Vapor phase dynamics were modeled as a point mass that communicated heat from the tank wall to the liquid phase via a boundary condition used at the gas-liquid interface. The method proved to be more accurate, robust, and efficient than explicit resolution of the dynamics using a standard Volume of Fluid (VOF) methodology. The subject pressurization process was found to be heavily dependent upon both the relatively high liquid temperature gradient near the gas-liquid interface and the natural convection flow path. Modeling the gas-liquid interface as an immovable surface eliminated temperature gradient destroying gas-liquid interface velocities observed in VOF simulations, and correspondingly enabled more rapid simulation since interface advection was not allowed. The single phase computational domain also facilitated the ability to demonstrate spatial resolution convergence of natural convection cells within the liquid which significantly impacted the tank pressurization rate. This work was used to demonstrate the critical physics for tank self-pressurization and numerical methodologies that may be used to best resolve those physics. The findings informed development and operation of production level CFD tools used in the Fluid Dynamics Branch at NASA Marshall Space Flight Center.

J. M. Brodnick↗
NASA NTRS2010

Thermal Protection System Cavity Heating for Simplified and Actual Geometries Using Computational Fluid Dynamics Simulations with Unstructured Grids

Thermal Protection System (TPS) Cavity Heating is predicted using Computational Fluid Dynamics (CFD) on unstructured grids for both simplified cavities and actual cavity geometries. Validation was performed using comparisons to wind tunnel experimental results and CFD predictions using structured grids. Full-scale predictions were made for simplified and actual geometry configurations on the Space Shuttle Orbiter in a mission support timeframe.

McCloud, Peter L.↗
NASA NTRS2016

Evaluation of Methods for Molecular Dynamics Simulation of Ionic Liquid Electric Double Layers

We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.

Haskins, Justin B.↗
NASA NTRS2019

Molecular Dynamics Simulations of Silicon Carbide, Boron Nitride and Silicon for Ceramic Matrix Composite Applications

A comprehensive computational molecular dynamics study is presented for crystalline α-SiC (6H, 4H, and 2H SiC), β-SiC (3C SiC), layered boron nitride, amorphous boron nitride and silicon, the constituent materials for high-temperature SiC/SiC compositions. Large-scale Atomic/Molecular Parallel Simulator software package was used. The Tersoff Potential force field was utilized to evaluate their mechanical characteristics of most of the materials, and the Reax force field was used to model silicon when the Tersoff Potential did not provide accurate results. Their mechanical behaviors were evaluated at a strain rate of 10(exp 7)/s and the results agree with the experimental data in the literature. The results are foundational for linking constituent behavior to composite performance, particularly when test data is unavailable or suspect.

Aluko, Olanrewaju↗
NASA NTRS2006

Dynamics Simulation Model for Space Tethers

This document describes the development of an accurate model for the dynamics of the Momentum Exchange Electrodynamic Reboost (MXER) system. The MXER is a rotating tether about 100-km long in elliptical Earth orbit designed to catch payloads in low Earth orbit and throw them to geosynchronous orbit or to Earth escape. To ensure successful rendezvous between the MXER tip catcher and a payload, a high-fidelity model of the system dynamics is required. The model developed here quantifies the major environmental perturbations, and can predict the MXER tip position to within meters over one orbit.

Levin, E. M.↗
NASA NTRS1986

Design and dynamic simulation of a fixed pitch 56 kW wind turbine drive train with a continuously variable transmission

The dynamic analysis of a horizontal axis fixed pitch wind turbine generator (WTG) rated at 56 kW is discussed. A mechanical Continuously Variable Transmission (CVT) was incorporated in the drive train to provide variable speed operation capability. One goal of the dynamic analysis was to determine if variable speed operation, by means of a mechanical CVT, is capable of capturing the transient power in the WTG/wind environment. Another goal was to determine the extent of power regulation possible with CVT operation.

Gallo, C.↗