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,297 records · Page 72

Quantum Algorithm for Simulating Hamiltonian Dynamics with an Off-diagonal Series Expansion

We propose an efficient quantum algorithm for simulating the dynamics of general Hamiltonian systems. Our technique is based on a power series expansion of the time-evolution operator in its off-diagonal terms. The expansion decouples the dynamics due to the diagonal component of the Hamiltonian from the dynamics generated by its off-diagonal part, which we encode using the linear combination of unitaries technique. Our method has an optimal dependence on the desired precision and, as we illustrate, generally requires considerably fewer resources than the current state-of-the-art. We provide an analysis of resource costs for several sample models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High frequency dynamic engine simulation

A digital computer simulation of a mixed flow, twin spool turbofan engine was assembled to evaluate and improve the dynamic characteristics of the engine simulation to disturbance frequencies of at least 100 Hz. One dimensional forms of the dynamic mass, momentum and energy equations were used to model the engine. A TF30 engine was simulated so that dynamic characteristics could be evaluated against results obtained from testing of the TF30 engine at the NASA Lewis Research Center. Dynamic characteristics of the engine simulation were improved by modifying the compression system model. Modifications to the compression system model were established by investigating the influence of size and number of finite dynamic elements. Based on the results of this program, high frequency engine simulations using finite dynamic elements can be assembled so that the engine dynamic configuration is optimum with respect to dynamic characteristics and computer execution time. Resizing of the compression systems finite elements improved the dynamic characteristics of the engine simulation but showed that additional refinements are required to obtain close agreement simulation and actual engine dynamic characteristics.

Schuerman, J. A.↗

Spectral solver for multi-scale plasma physics simulations with dynamically adaptive number of moments

A spectral method for kinetic plasma simulations based on the expansion of the velocity distribution function in a variable number of Hermite polynomials is presented. The method is based on a set of non-linear equations that is solved to determine the coefficients of the Hermite expansion satisfying the Vlasov and Poisson equations. In this paper, we first show that this technique combines the fluid and kinetic approaches into one framework. Second, we present an adaptive strategy to increase and decrease the number of Hermite functions dynamically during the simulation. The technique is applied to the Landau damping and two-stream instability test problems. Performance results show 21% and 47% saving of total simulation time in the Landau and two-stream instability test cases, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DYGABCD: A program for calculating linear A, B, C, and D matrices from a nonlinear dynamic engine simulation

A digital computer program, DYGABCD, was developed that generates linearized, dynamic models of simulated turbofan and turbojet engines. DYGABCD is based on an earlier computer program, DYNGEN, that is capable of calculating simulated nonlinear steady-state and transient performance of one- and two-spool turbojet engines or two- and three-spool turbofan engines. Most control design techniques require linear system descriptions. For multiple-input/multiple-output systems such as turbine engines, state space matrix descriptions of the system are often desirable. DYGABCD computes the state space matrices commonly referred to as the A, B, C, and D matrices required for a linear system description. The report discusses the analytical approach and provides a users manual, FORTRAN listings, and a sample case.

Geyser, L. C.↗

Efficient massively parallel simulation of dynamic channel assignment schemes for wireless cellular communications

Fast, efficient parallel algorithms are presented for discrete event simulations of dynamic channel assignment schemes for wireless cellular communication networks. The driving events are call arrivals and departures, in continuous time, to cells geographically distributed across the service area. A dynamic channel assignment scheme decides which call arrivals to accept, and which channels to allocate to the accepted calls, attempting to minimize call blocking while ensuring co-channel interference is tolerably low. Specifically, the scheme ensures that the same channel is used concurrently at different cells only if the pairwise distances between those cells are sufficiently large. Much of the complexity of the system comes from ensuring this separation. The network is modeled as a system of interacting continuous time automata, each corresponding to a cell. To simulate the model, conservative methods are used; i.e., methods in which no errors occur in the course of the simulation and so no rollback or relaxation is needed. Implemented on a 16K processor MasPar MP-1, an elegant and simple technique provides speedups of about 15 times over an optimized serial simulation running on a high speed workstation. A drawback of this technique, typical of conservative methods, is that processor utilization is rather low. To overcome this, new methods were developed that exploit slackness in event dependencies over short intervals of time, thereby raising the utilization to above 50 percent and the speedup over the optimized serial code to about 120 times.

Greenberg, Albert G.↗

On coarse projective integration for atomic deposition in amorphous systems

Direct molecular dynamics simulation of atomic deposition under realistic conditions is notoriously challenging because of the wide range of time scales that must be captured. Numerous simulation approaches have been proposed to address the problem, often requiring a compromise between model fidelity, algorithmic complexity, and computational efficiency. Coarse projective integration, an example application of the “equation-free” framework, offers an attractive balance between these constraints. Here, periodically applied, short atomistic simulations are employed to compute time derivatives of slowly evolving coarse variables that are then used to numerically integrate differential equations over relatively large time intervals. A key obstacle to the application of this technique in realistic settings is the “lifting” operation in which a valid atomistic configuration is recreated from knowledge of the coarse variables. Using Ge deposition on amorphous SiO{sub 2} substrates as an example application, we present a scheme for lifting realistic atomistic configurations comprised of collections of Ge islands on amorphous SiO{sub 2} using only a few measures of the island size distribution. The approach is shown to provide accurate initial configurations to restart molecular dynamics simulations at arbitrary points in time, enabling the application of coarse projective integration for this morphologically complex system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation Exploration through Immersive Parallel Planes

We present a visualization-driven simulation system that tightly couples systems dynamics simulations with an immersive virtual environment to allow analysts to rapidly develop and test hypotheses in a high-dimensional parameter space. To accomplish this, we generalize the two-dimensional parallel-coordinates statistical graphic as an immersive 'parallel-planes' visualization for multivariate time series emitted by simulations running in parallel with the visualization. In contrast to traditional parallel coordinate's mapping the multivariate dimensions onto coordinate axes represented by a series of parallel lines, we map pairs of the multivariate dimensions onto a series of parallel rectangles. As in the case of parallel coordinates, each individual observation in the dataset is mapped to a polyline whose vertices coincide with its coordinate values. Regions of the rectangles can be 'brushed' to highlight and select observations of interest: a 'slider' control allows the user to filter the observations by their time coordinate. In an immersive virtual environment, users interact with the parallel planes using a joystick that can select regions on the planes, manipulate selection, and filter time. The brushing and selection actions are used to both explore existing data as well as to launch additional simulations corresponding to the visually selected portions of the input parameter space. As soon as the new simulations complete, their resulting observations are displayed in the virtual environment. This tight feedback loop between simulation and immersive analytics accelerates users' realization of insights about the simulation and its output.

97 MATHEMATICS AND COMPUTING↗

Construction of dynamic stochastic simulation models using knowledge-based techniques

Over the past three decades, computer-based simulation models have proven themselves to be cost-effective alternatives to the more structured deterministic methods of systems analysis. During this time, many techniques, tools and languages for constructing computer-based simulation models have been developed. More recently, advances in knowledge-based system technology have led many researchers to note the similarities between knowledge-based programming and simulation technologies and to investigate the potential application of knowledge-based programming techniques to simulation modeling. The integration of conventional simulation techniques with knowledge-based programming techniques is discussed to provide a development environment for constructing knowledge-based simulation models. A comparison of the techniques used in the construction of dynamic stochastic simulation models and those used in the construction of knowledge-based systems provides the requirements for the environment. This leads to the design and implementation of a knowledge-based simulation development environment. These techniques were used in the construction of several knowledge-based simulation models including the Advanced Launch System Model (ALSYM).

Williams, M. Douglas↗

Model based co-simulation platform for integrated building system control and design optimization

Both steady-state and dynamic simulations have been widely used by HVAC&R industry to support product/equipment development for decades. Steady-state simulation focuses on the system mass, energy and momentum balance of an equilibrium state. It is based on high-fidelity components models, and thus is suitable for system and component design optimization. Dynamic simulation studies the system transient response and is generally used for controls development and verification. It usually does not require rigorous component models of high accuracy because 1) the commonly used PID control is feedback control whose control performance evaluation doesn’t require high fidelity system/plant model; 2) high-fidelity dynamic model significantly increases the number of equations and variables and creates tremendous challenge for math solver. For supervisory control, transactive control or optimization of an integrated building system, the HVAC&R equipment is often one of the sub-components to be controlled. High-fidelity equipment models are required for accurately evaluating control strategies. In addition, building equipment manufacturers have developed a lot of high-fidelity steady-state equipment/component models per their expertise. Thus, a platform that can integrate OEM high-fidelity steady-state model with dynamic building simulation and/or electric power system & grid simulation to support the development and verification of supervisory control for integrated building systems is necessary. In this study, ORNL’s heat pump design tool (HPDM) is utilized to develop a co-simulation platform for supervisory control and optimization in integrated building systems. It is based on a model that integrates high-fidelity steady-state simulation equipment models with dynamic building simulation. A practical case of using the proposed co-simulation platform to develop and evaluate the supervisory control and optimization is presented and discussed.

Sun, Jian↗

Generalized dynamic engine simulation techniques for the digital computer

Recently advanced simulation techniques have been developed for the digital computer and used as the basis for development of a generalized dynamic engine simulation computer program, called DYNGEN. This computer program can analyze the steady state and dynamic performance of many kinds of aircraft gas turbine engines. Without changes to the basic program DYNGEN can analyze one- or two-spool turbofan engines. The user must supply appropriate component performance maps and design-point information. Examples are presented to illustrate the capabilities of DYNGEN in the steady state and dynamic modes of operation. The analytical techniques used in DYNGEN are briefly discussed, and its accuracy is compared with a comparable simulation using the hybrid computer. The impact of DYNGEN and similar all-digital programs on future engine simulation philosophy is also discussed.

Sellers, J.↗

Generalized dynamic engine simulation techniques for the digital computer

Recently advanced simulation techniques have been developed for the digital computer and used as the basis for development of a generalized dynamic engine simulation computer program, called DYNGEN. This computer program can analyze the steady state and dynamic performance of many kinds of aircraft gas turbine engines. Without changes to the basic program, DYNGEN can analyze one- or two-spool turbofan engines. The user must supply appropriate component performance maps and design-point information. Examples are presented to illustrate the capabilities of DYNGEN in the steady state and dynamic modes of operation. The analytical techniques used in DYNGEN are briefly discussed, and its accuracy is compared with a comparable simulation using the hybrid computer. The impact of DYNGEN and similar all-digital programs on future engine simulation philosophy is also discussed.

Sellers, J.↗

Generalized dynamic engine simulation techniques for the digital computers

Recently advanced simulation techniques have been developed for the digital computer and used as the basis for development of a generalized dynamic engine simulation computer program, called DYNGEN. This computer program can analyze the steady state and dynamic performance of many kinds of aircraft gas turbine engines. Without changes to the basic program, DYNGEN can analyze one- or two-spool turbofan engines. The user must supply appropriate component performance maps and design point information. Examples are presented to illustrate the capabilities of DYNGEN in the steady state and dynamic modes of operation. The analytical techniques used in DYNGEN are briefly discussed, and its accuracy is compared with a comparable simulation using the hybrid computer. The impact of DYNGEN and similar digital programs on future engine simulation philosophy is also discussed.

Sellers, J.↗

Dynamic CFD Simulations of the MEADS II Ballistic Range Test Model

Dynamic, viscous, free-to-oscillate simulations of the Mars Entry Atmospheric Data System (MEADS) ballistic range model are performed using two different ow solvers, OVERFLOW and US3D. At the time of publication, data from the ballistic range experiment was not yet available, so the current work serves as a code-to-code exercise. Results from the analysis show good agreement between the predicted static aerodynamic coefficients for each solver. Both codes predict damped pitch oscillations for Mach 3:0 with initial amplitudes of 5⁰ and 30⁰, as well as for Mach 1:5 with initial amplitude of 30⁰. The two solvers predict undamped pitch oscillations for Mach 1:5 with initial amplitude of 5⁰. For most cases, US3D predicts less damping than OVERFLOW. The difference is attributed to higher pressures in the separated region of the wake, and the resultant effect on the backshell contribution to the pitching moment.

Stern, Eric C.↗

Auto_PDI

Protein-DNA Interaction Workflow (PDI Workflow), a pipeline that focuses on generating high-quality docking and molecular dynamics simulations for Protein-DNA complexes. This allows us to take DNA sequences with unknown tertiary structures, accurately predict their structure, dock them with the desired target protein, and then simulate their interactions using molecular dynamics simulations

Kumar, Neeraj↗

fp-tools

The output of a molecular dynamics simulation is a data file containing properties of the motion of the simulated particles as a function of time (e.g., positions, velocities). From this output, scientists can extract various properties which give insight into the physics of the system; for example, the radial distribution function can be computed from the particle positions, which gives insight into how close the system is to melting or freezing. Many publicly available codes to perform molecular dynamics simulations exist and have been well adopted by the scientific community (e.g., LAMMPS, VASP). However, it remains standard practice for scientists to write their own post-processing scripts to extract various properties of interest from the molecular dynamics simulation data. Since in many cases scientists are interested in the same set of “textbook” properties, there is a large duplication of effort in writing these post-processing scripts. Our proposed code, FP-Tools, is a C++ toolkit designed to help scientists extract and analyze quantities of interest from the output of a molecular dynamics simulation, eliminating the need for them to write these codes themselves and thus reducing this duplication of effort. The quantities which are computed by our code are well documented in the literature (either in publications or textbooks), and the algorithms we implement to compute these properties are also well known in the field. We are not introducing new science or methods here; rather, our primary goal is to provide a useful tool to the community in the form of a well-documented, easy to use package.

Hartman, Leah↗

Superhydrophilicity of $α$-alumina surfaces results from tight binding of interfacial waters to specific aluminols

Understanding the microscopic driving force of water wetting is challenging and important for design of materials. The relations between structure, dynamics and hydrogen bonds of interfacial water can be investigated using molecular dynamics simulations. Here, contact angles at the alumina (0001) and ($11\bar{2}0$) surfaces are studied using both classical molecular dynamics simulations and experiments. To test the superhydrophilicity, the free energy cost of removing waters near the interfaces are calculated using the density fluctuations method. The strength of hydrogen bonds is determined by their lifetime and geometry. Both surfaces are superhydrophilic and the (0001) surface is more hydrophilic. Interactions between surfaces and interfacial waters promote a templating effect whereby the latter are aligned in a pattern that follows the underlying lattice of the surfaces. Translational and rotational dynamics of interfacial water molecules are slower than in bulk water. Hydrogen bonds between water and both surfaces are asymmetric, water-to-aluminol ones are stronger than aluminol-to-water ones. Molecular dynamics simulations eliminate the impacts of surface contamination when measuring contact angles and the results reveal the microscopic origin of the macroscopic superhydrophilicity of alumina surfaces: strong water-to-aluminol hydrogen bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tipping Point for Expansion of Layered Aluminosilicates in Weakly Polar Solvents: Supercritical CO2

Layered aluminosilicates play a dominant role in the mechanical and gas storage properties of the subsurface, are used in diverse industrial applications, and serve as model materials for understanding solvent-ion-support systems. Although expansion in the presence of H2O is well known to be systematically correlated with the hydration free energy of the interlayer cation, in environments dominated by non-polar solvents (i.e. CO2), uptake into the interlayer is not well-understood. Using novel high pressure capabilities, we investigated the interaction of super-critical CO2 with Na+-, NH4+-, and Cs+-saturated montmorillonite, comparing results with predictions from molecular dynamics simulations. Despite the known trend in H2O, and that cation solvation energies in CO2 suggest a stronger interaction with Na+, both the NH4+- and Cs+-clays readily absorbed CO2 and expanded while the Na+-clay did not. The apparent inertness of the Na+-clay was not due to kinetics, as experiments seeking a stable expanded state showed that none exists. Molecular dynamics simulations revealed a large endothermicity to CO2 intercalation in the Na+-clay, but little or no energy barrier for the NH4+- and Cs+-clays. Consequently, we have shown for the first time that in the presence of a low dielectric constant gas swelling depends more on the strength of the interaction between interlayer cation and aluminosilicate sheets and less on that with solvent. The finding suggests a distinct regime in layered aluminosilicates swelling behavior triggered by low solvent polarizability, with important implications in geomechanics, storage and retention of volatile gases, and across industrial uses in gelling, decoloring, heterogeneous catalysis, and semi-permeable reactive barriers.

CO2, montmorillonite, molecular dynamics↗

Mapping 180° polar domains using electron backscatter diffraction and dynamical scattering simulations

A novel technique, which directly and nondestructively maps polar domains using electron backscatter diffraction (EBSD) is described and demonstrated. Through dynamical diffraction simulations and quantitative comparison to experimental EBSD patterns, the absolute orientation of a non-centrosymmetric crystal can be determined. With this information, the polar domains of a material can be mapped. The technique is demonstrated by mapping the non-ferroelastic, or 180°, ferroelectric domains in periodically poled LiNbO 3 single crystals. Furthermore, the authors demonstrate the possibility of mapping polarity using this technique in other polar materials system.

36 MATERIALS SCIENCE↗