Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,567 records · Page 87

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↗

A model for simulating adaptive, dynamic flows on networks: Application to petroleum infrastructure

Simulation models can greatly improve decisions meant to control the consequences of disruptions to critical infrastructures. We describe a dynamic flow model on networks purposed to inform analyses by those concerned about consequences of disruptions to infrastructures and to help policy makers design robust mitigations. We conceptualize the adaptive responses of infrastructure networks to perturbations as market transactions and business decisions of operators. We approximate commodity flows in these networks by a diffusion equation, with nonlinearities introduced to model capacity limits. To illustrate the behavior and scalability of the model, we show its application first on two simple networks, then on petroleum infrastructure in the United States, where we analyze the effects of a hypothesized earthquake.

02 PETROLEUM↗

Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil" in the Journal of Geophysical Research-Biogeosciences. We simulated biogeochemical cycling in arctic soils using the PFLOTRAN geochemical model combined with measurements from previous NGEE Arctic incubations of polygonal permafrost soils in northern Alaska (Zheng et al., 2018). Simulated iron cycling, carbon dioxide production, and methane production were compared with incubation measurements and the parameterized model was then used to simulate coupled iron and carbon cycling over repeated oxic-anoxic cycles at different levels of carbon substrate availability and pH. The most recent data version (2.0) in the archive incorporates changes to the model and simulations as suggested by reviewers during the manuscript review process. These changes include an updated parameterization of the model; a new set of simulations omitting the iron cycle for direct evaluation of how iron cycle processes affect modeled outcomes; and a set of simulations testing different scenarios of carbon substrate availability in addition to scenarios of initial soil pH. This archive contains simulation code, model output, and analysis code for PFLOTRAN simulations. All scripts are python except the batch script for submitting multiprocessor jobs. Note that the model also requires compiled versions of the Alquimia interface and the NGEE Arctic fork of the PFLOTRAN geochemical simulator (see the README_INSTALL document for basic instructions). The Output directory contains eight data files in netCDF format generated by the model. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 10-year research effort (2012-2022) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Structural properties of Y 2 O 3 –Al 2 O 3 liquids and glasses: An overview

Liquids in the system Y 2 O 3 –Al 2 O 3 have been the subject of considerable study because of the reported occurrence of a first-order density and entropy-driven liquid-liquid phase transition (LLPT) in the supercooled liquid state. The observations have become controversial because of the presence of crystalline material that can be formed simultaneously and that can mask the nucleation and growth of the lower density liquid. The previous work is summarized here along with arguments for and against the different viewpoints. Also two studies have been undertaken to investigate the LLPT in this refractory system with emphasis on determining the structure of unequivocally amorphous materials. These include the in situ high energy X-ray diffraction (HEXRD) of supercooled Y 2 O 3 –Al 2 O 3 liquids and the low frequency vibrational dynamics of recovered glasses. Manybody molecular dynamics simulations are also used to interpret the results of both studies. The HEXRD measurements, combined with aerodynamic levitation and rapid data acquisition techniques, show that for the 20 mol% Y 2 O 3 (i.e. AlY20) liquid there is a shift in the position of the first peak in the diffraction pattern over a narrow temperature range (2100-1800 K) prior to crystallization. Microbeam Raman spectroscopy measurements made on AlY20 glasses clearly show contrasting spectra in the low frequency part of the spectrum for low(LDA) and high-density (HDA) glassy regions. The molecular dynamics simulations identify contrasting coordination environments around oxygen anions for the high- (HDL) and low-density (LDL) liquids. (C) 2014 Elsevier B.V. All rights reserved.

36 MATERIALS SCIENCE↗

Coupling LAMMPS and the vl3 Framework for Co-Visualization of Atomistic Simulations

In this work we perform in-situ visualization of molecular dynamics simulations, which can help scientists to visualize simulation output on-the-fly, without incurring storage overheads. We present a case study to couple LAMMPS, the large-scale molecular dynamics simulation code with vl3, our parallel framework for large-scale visualization and analysis. Our motivation is to identify effective approaches for covisualization and exploration of large-scale atomistic simulations at interactive frame rates.We propose a system of coupled libraries and describe its architecture, with an implementation that runs on GPU-based clusters. We present the results of strong and weak scalability experiments, as well as future research avenues based on our results.

co-visualization↗

Data-assisted combustion simulations with dynamic submodel assignment using random forests

This investigation outlines a data-assisted approach that employs random forest classifiers for local and dynamic submodel assignment in turbulent-combustion simulations. This method is demonstrated in simulations of a single-element GOX/GCH4 rocket combustor; a priori as well as a posteriori assessments are conducted to (i) evaluate the accuracy and adjustability of the classifier for targeting different quantities of interest (QoIs), and (ii) assess improvements, resulting from the data-assisted combustion model assignment, in predicting target QoIs during simulation runtime. Results from the a priori study show that random forests, trained with local flow properties as input variables and combustion model errors as training labels, assign three different combustion models – finite-rate chemistry (FRC), flamelet progress variable (FPV) model, and inert mixing (IM) – with reasonable classification performance even when targeting multiple QoIs. Applications in a posteriori studies demonstrate improved predictions from data-assisted simulations, in temperature and CO mass fraction, when compared with monolithic FPV calculations. An additional a posteriori data-assisted simulation of a modified configuration demonstrates that the present approach can be successfully applied to different configurations, as long as thermophysical behavior can be represented by the training data. Furthermore, these results demonstrate that this data-driven framework holds promise for dynamic combustion submodel assignments in reacting flow simulations.

42 ENGINEERING↗

Discussion on molecular dynamics (MD) simulations of the asphalt materials

The application of asphalt materials in pavement engineering has been increasingly widespread and sophisticated over the past several decades. Variations in the properties of asphalt binder during mixing, transportation, and paving can affect the performance of asphalt pavement. However, the asphalt material is a non-homogeneous and complex organic substance, consisting of various molecules with widely various molecular weights, elemental compositions, and structures. This complexity leads to difficulties for researchers to clearly and immediately understand the properties of asphalt materials and their variations. The multi-scale research approach combines macroscopic experimental data and microscopic simulation results from a practical engineering perspective. It helps to improve the understanding of asphalt materials. The molecular dynamics (MD) simulation proposes a corresponding molecular model of asphalt material based on experimental data, and the simulation algorithm is able to derive properties similar to those of real asphalt. Here, this paper provides a comprehensive review of the current studies on MD simulation of asphalt materials, including modeling, properties, and multi-scale analysis. As a key part of the computational simulation, this paper discusses the typical asphalt binder and asphalt-aggregate interface models constructed by different groups, and also presents their differences from real samples and their feasibility based on fundamental properties. After the introduction of molecular models, the extensive work made by researchers based on molecular models is categorically reviewed and discussed. The strengths and weaknesses of MD simulation methods in the study of asphalt materials are also summarized in order to provide the reader with a more comprehensive understanding of the relevant contents and to guide subsequent research.

36 MATERIALS SCIENCE↗

Off-fault deformations and shallow slip deficit from dynamic rupture simulations with fault zone plasticity

Kinematic source inversions of major (M≥7) strike-slip earthquakes show that the slip at depth exceeds surface displacements measured in the field, and it has been suggested that this shallow slip deficit (SSD) is caused by distributed plastic deformation near the surface. We perform dynamic rupture simulations of M 7.2–7.4 earthquakes in elastoplastic media and analyze the sensitivity of SSD and off-fault deformation (OFD) to rock quality parameters. While linear simulations clearly underpredict observed SSD and OFDs, nonlinear simulations for a moderately fractured fault damage zone predict a SSD of 44–53% and OFDs of 39–48%, consistent with the 30–60% SSD and 46 ± 10% (1σ) OFD reported for the 1992 M 7.3 Landers earthquake. Finally, both SSD and OFDs are sensitive to the quality of the fractured rock mass inside the fault damage zone, and surface rupture is almost entirely suppressed in poor quality material.

58 GEOSCIENCES↗

Data-driven material models for atomistic simulation

The central approximation made in classical molecular dynamics simulation of materials is the interatomic potential used to calculate the forces on the atoms. Great effort and ingenuity is required to construct viable functional forms and find accurate parametrizations for potentials using traditional approaches. Machine learning has emerged as an effective alternative approach to develop accurate and robust interatomic potentials. Beginning with a very general model form, the potential is learned directly from a database of electronic structure calculations and therefore can be viewed as a multiscale link between quantum and classical atomistic simulations. Risk of inaccurate extrapolation exists outside the narrow range of time and length scales where the two methods can be directly compared. Here, we use the spectral neighbor analysis potential (SNAP) and show how a fit can be produced with minimal interpolation errors which is also robust in extrapolating beyond training. To demonstrate the approach, we have developed a tungsten-beryllium potential suitable for the full range of binary compositions. Subsequently, large-scale molecular dynamics simulations were performed of high energy Be atom implantation onto the (001) surface of solid tungsten. The machine learned W-Be potential generates a population of implantation structures consistent with quantum calculations of defect formation energies. A very shallow ( < 2 nm ) average Be implantation depth is predicted which may explain ITER diverter degradation in the presence of beryllium.

36 MATERIALS SCIENCE↗

A Comparative Study of Interface Techniques for Transmission and Distribution Dynamic Co-Simulation

Transmission and distribution (T&D) dynamic co-simulation is a practical and effective approach to leverage existing simulation tools for transmission and distribution systems to simulate dynamic stability and performance of T&D systems in a systematic manner. Given that these tools are developed as stand-alone programs and there are inherent differences between them, interface techniques become critical to “bridge” them. Two important unsolved questions are: 1) which interface technique is better and should be used, and 2) how the modeling and simulation capabilities in these tools that are available and can be exploited for co-simulation should be considered when selecting an interface technique. This paper presents a comprehensive study for different interface techniques that can be employed for T&D dynamic co-simulation. The study provided insights to the pros and cons of each interface technique, and helped researchers make informed decisions on choosing a proper interface technique.

Huang, Qiuhua↗