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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Data-Driven Security Assessment of Power Grids Based on Machine Learning Approach: Preprint

Data-driven security assessment provides key indicators on power system stability using simulations on scheduling models, as opposed to dynamic simulations that are more time-consuming. This paper investigates data-driven security assessment of power grids based on machine learning. Multivariate random forest regression is used as the machine learning algorithm due to its high robustness to the input data. Three stability issues are analyzed using the proposed machine learning tool, including transient stability, frequency stability and small signal stability. The estimation values from machine learning tool are compared with those from dynamic simulations. Results show that the proposed machine learning tool can effectively predict the stability margins for the three stability metrics.

14 SOLAR ENERGY↗
DOE OSTI2018

Molecular dynamics studies of fundamental bulk properties of palladium hydrides for hydrogen storage

Solid-state hydrogen storage materials undergo complex phase transformations whose behavior are collectively determined by thermodynamic (e.g., Gibbs free energy), mechanical (e.g., lattice and elastic constants), and mass transport (e.g., diffusivity) properties. These properties depend on the reaction conditions and evolve continuously during (de)hydrogenation. Thus, they are difficult to measure in experiments. Because of this, past progress to improve solid-state hydrogen storage materials has been prolonged. Using PdH x as a representative example for interstitial metal hydride, we have recently applied molecular dynamics simulations to quantify hydrogen diffusion in the entire reaction space of temperature and composition. In this paper, we have further applied molecular dynamics simulations to obtain well-converged expressions for lattice constants, Gibbs free energies, and elastic constants of PdH x at various stages of the reaction. Finally, our studies confirm significant dependence of elastic constants on temperature and composition. Specifically, a new dynamic effect of hydrogen diffusion on elastic constants is discovered and discussed.

08 HYDROGEN↗
DOE OSTI2023

Simulation-Based Recovery Action Analysis Using the EMRALD Dynamic Risk Assessment Tool

Recovery human action is defined as the action that prevents deviant conditions from producing unwanted effects. Analyzing recovery actions has been a critical part in human reliability analysis (HRA). However, there are a couple of limitations to treating recovery actions using only the current HRA methods available. Representatively, the existing recovery analysis does not specifically consider recovery actions as they have occurred in actual nuclear power plants (NPPs). To handle the challenges in the existing recovery analyses, this study suggests a way to analyze recovery actions under a dynamic HRA method, the Procedure-based Risk Investigation MEthod-Human Reliability Analysis (PRIME-HRA) method. The PRIME-HRA method suggests a way on how to develop dynamic simulation models using dynamic risk assessment tools such as the Event Modeling Risk Assessment Using Linked Diagram (EMRALD) [1] and the Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER) [2]. EMRALD and HUNTER are the dynamic probabilistic risk assessment and HRA tools developed at Idaho National Laboratory. In this paper, differences on analyzing recovery actions in the Technique for Human Error-Rate Prediction (THERP), the Cause-Based Decision Tree (CBDT) and the Korean Standard HRA (K-HRA) and challenges of these approaches are introduced. How we have developed the PRIME-HRA is also introduced in this paper. Then, the proposed approach to analyzing recovery human actions in dynamic context is partially discussed with an example.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS2016

Evolution of Flexible Multibody Dynamics for Simulation Applications Supporting Human Spaceflight

During the course of transition from the Space Shuttle and International Space Station programs to the Orion and Journey to Mars exploration programs, a generic flexible multibody dynamics formulation and associated software implementation has evolved to meet an ever changing set of requirements at the NASA Johnson Space Center (JSC). Challenging problems related to large transitional topologies and robotic free-flyer vehicle capture/ release, contact dynamics, and exploration missions concept evaluation through simulation (e.g., asteroid surface operations) have driven this continued development. Coupled with this need is the requirement to oftentimes support human spaceflight operations in real-time. Moreover, it has been desirable to allow even more rapid prototyping of on-orbit manipulator and spacecraft systems, to support less complex infrastructure software for massively integrated simulations, to yield further computational efficiencies, and to take advantage of recent advances and availability of multi-core computing platforms. Since engineering analysis, procedures development, and crew familiarity/training for human spaceflight is fundamental to JSC's charter, there is also a strong desire to share and reuse models in both the non-realtime and real-time domains, with the goal of retaining as much multibody dynamics fidelity as possible. Three specific enhancements are reviewed here: (1) linked list organization to address large transitional topologies, (2) body level model order reduction, and (3) parallel formulation/implementation. This paper provides a detailed overview of these primary updates to JSC's flexible multibody dynamics algorithms as well as a comparison of numerical results to previous formulations and associated software.

Multibody dynamics↗
DOE OSTI2017

Quantifying Economic and Environmental Impacts of Transportation Network Disruptions with Dynamic Traffic Simulation

Several transportation network vulnerability models have been proposed. However, most only consider disruptions as a static snapshot in time and the impact on total travel time. These approaches cannot consider the time-varying nature of travel demand nor other undesirable outcomes that follow from transportation network disruptions. This paper proposes an algorithmic approach to assess the vulnerability of a transportation network that considers the time-varying demand with an open source dynamic transportation simulation tool. The open source nature of the tool allows us to systematically consider many disruption scenarios and quantitatively compare their relative criticality. This is far more efficient than traditional approaches which would require days or weeks of a transportation engineers time to manually set up, run, and assess these simulations. In addition to travel time, we also collect statistics on additional fuel consumed and the corresponding carbon dioxide emissions. Our approach, thus provides a more systematic approach that is both time-varying and can consider additional negative consequences of disruptions for decision makers to evaluate.

Shekar, Venkateswaran↗
DOE OSTI2016

Nanocarbon synthesis by high-temperature oxidation of nanoparticles

High-temperature oxidation of silicon-carbide nanoparticles (nSiC) underlies a wide range of technologies from high-power electronic switches for efficient electrical grid and thermal protection of space vehicles to self-healing ceramic nanocomposites. Here, multimillion-atom reactive molecular dynamics simulations validated by ab initio quantum molecular dynamics simulations predict unexpected condensation of large graphene flakes during high-temperature oxidation of nSiC. Initial oxidation produces a molten silica shell that acts as an autocatalytic ‘nanoreactor’ by actively transporting oxygen reactants while protecting the nanocarbon product from harsh oxidizing environment. Percolation transition produces porous nanocarbon with fractal geometry, which consists of mostly sp 2 carbons with pentagonal and heptagonal defects. Furthermore, this work suggests a simple synthetic pathway to high surface-area, low-density nanocarbon with numerous energy, biomedical and mechanical-metamaterial applications, including the reinforcement of self-healing composites.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1999

Simulation of Dynamics of PVT Growth: ZnSe

An unsteady model for the simulation of PVT growth process has been developed and applied to growth of ZnSe. The model is capable of tracking the unsteady aspects of the growth process caused by, for example, changes in the system temperature associated with growth, source-depletion, translation of the charge, out-gassing at high temperatures, as well as changes in the source composition during growth. The model predicts the evolution of growth interface morphology as well as the shape of source-gas interface. Thermo-fluid transport in the gas phase is handled by a Low-Mach number unsteady compressible flow formulation. Heat transfer in the crystal, source, and the quartz housing is calculated by the Monte-Carlo technique to capture the variations in the spectral transmittance and index of refraction of the participating media. Simulation results are presented for the growth interface morphology, thermal stresses in the crystal, and growth rate at various stages of growth. The influence of impurity gases on the growth process is also discussed, The simulation results are compared with the experimental observations. The capabilities of the developed tool in addressing other PVT growth processes such as SiC are discussed.

Worlikar, A.↗
NASA NTRS2018

ASPIRE Flight Mechanics Modeling and Post Flight Analysis

The Advanced Supersonic Parachute Inflation Research and Experiment (ASPIRE) is a series of sounding rocket flights aimed at understanding the dynamics of supersonic parachutes that are used for Mars robotic applications. SR01 was the first sounding rocket flight of ASPIRE that occurred off the coast of Wallops Island, VA on Oct. 4, 2017 and showed the successful deployment and inflation of a Mars Science Laboratory built-to- print parachute in flight conditions similar to the 2012 Mars Science Laboratory (MSL) mission. SR02 was the second sounding rocket flight that also occurred off the coast of Wallops Island on March 31, 2018 and showcased the successful deployment and inflation of a new strengthened parachute being considered for the Mars 2020 mission at fifty percent higher dynamic pressure than observed on MSL. Prior to both flights, a multi-body flight dynamics simulation was developed to predict the parachute dynamics and was used, in conjunction with other tools, to target Mars-relevant flight conditions. After each flight, the reconstructed trajectory was used to validate the pre-flight dynamics simulation and recommend changes to improve predictions for future flights planned for the ASPIRE pro- gram. This paper describes the flight mechanics simulation and the post flight reconciliation process used to validate the flight models.

Dutta, Soumyo↗
NASA NTRS2020

ASPIRE Parachute Modeling and Comparison to Post-Flight Reconstruction

The Advanced Supersonic Parachute Inflation Research and Experiment (ASPIRE) was a series of sounding rocket flights aimed at understanding the dynamics of supersonic parachutes that are used for Mars robotic applications. Three flights for ASPIRE occurred off the coast of Wallops Island, VA in Oct. 2017, Mar. 2018, and Sept. 2018 and successfully demonstrated deployment and inflation of the Mars Science Laboratory and Mars 2020 mission parachute. Prior to all three flights, a multi-body flight dynamics simulation was developed to predict the parachute dynamics and was used, in conjunction with other tools, to target Mars-relevant flight conditions. After each flight, the reconstructed trajectory was used to validate the preflight dynamics simulation and recommend changes to improve predictions for future flights in the ASPIRE program. This paper describes the parachute models and flight mechanics simulation used to target conditions for the three flights and the post-flight comparison of the tools.

Soumyo Dutta↗
NASA NTRS2020

Advanced Supersonic Parachute Inflation Research Experiment Preflight Trajectory Modeling and Postflight Reconstruction

The Advanced Supersonic Parachute Inflation Research and Experiments (ASPIRE) was a series of sounding rocket flights aimed at understanding the dynamics of supersonic parachutes that are used for Mars robotic applications. The 2012 Mars Science Laboratory (MSL) had a successful deployment of a supersonic parachute, but based on post-flight analysis of parachute margins, the ASPIRE project was created as a risk-reduction program to improve quantification of these margins and qualify a supersonic parachute for Mars 2020, the follow-on mission to MSL. The first sounding rocket (SR01) flight of ASPIRE occurred near Wallops Island, Virginia on Oct. 4, 2017 and demonstrated the successful deployment and inflation of a MSL build-to-print parachute in flight conditions similar to the 2012 MSL mission. The ASPIRE SR02 and SR03 were successful follow-on flights on Mar. 31, 2018 and Sep. 7, 2018 that demonstrated the new, strengthened supersonic parachute designed for the Mars 2020 project. The SR02 and SR03 parachuteswere targeted to 100% and 140% of the expected flight limit load for Mars 2020 to confirm new margins expected from the strengthened parachute. Prior to all flights, a multi-body flight dynamics simulation was developed to predict the parachute dynamics and was used, in conjunction with other tools, to target Mars-relevant flight conditions. After each flight, the on-board data were used to reconstruct the flight trajectory and to validate the pre-flight dynamics simulation. Post-flight analysis showed that all three tests achieved their targeted conditions and pre-flight modeling bounded the key performance metrics for the parachute. This paper describes the flight mechanics simulation, post-flight reconstruction, and the reconciliation process used to validate the flight models.

Soumyo Dutta↗
DOE OSTI2018

Quantitative Prediction of Uranium Speciation and Amidoxime Binding in Seawater from Advanced Simulation Techniques

The goal of this project was to quantitatively predict stability of uranium species and thermochemistry of uranium binding with amidoxime ligands in realistic seawater conditions. The specific aims of the project were: (i) determine the stable species of uranium and its solvation structure in realistic seawater conditions; also speciation and solvation of uranium's major competing species in seawater, vanadium; (ii) accurately predict the free energies (and thereby stability constants) of binding between uranium and amidoxime; (iii) accurately predict the reaction enthalpies and entropies of binding between uranium and amidoxime as well as between vanadium and amidoxime. Main results and conclusions are as follows. First-principles molecular dynamic simulations found that the structure of the Ca 2 UO 2 (CO 3 ) 3 complex is very stable and that one Ca ion binds to the center UO 2 (CO 3 ) 3 4- anion stronger than the other Ca ion. This finding suggests that using time-resolved EXAFS spectra may confirm the asymmetry in binding of the two Ca ions in the aqueous Ca 2 UO 2 (CO 3 ) 3 complex. To consider the common ions, solvation of the Ca 2 UO 2 (CO 3 ) 3 complex in seawater was simulated by classical molecular dynamics simulation. It was found that the structure of the Ca 2 UO 2 (CO 3 ) 3 complex is very stable in the model seawater. A Na+ ion was found to be closely associated with the Ca 2 UO 2 (CO 3 ) 3 complex by indirectly interacting with one axial oxygen atom of the UO 2 group bridged by a water molecule. In addition, the Na + ion interacts closely with one Ca 2+ ion than the other. The present simulations revealed the key role of common ions such as Na + in impacting the solvation, structure, and apparent charge of the Ca 2 UO 2 (CO 3 ) 3 complex in seawater. The sequential processes of UO 2 2+ binding with the three CO 3 2- groups and the two Ca 2+ ions in pure water and in 0.1 M NaCl were simulated by classical molecular dynamics with both the non-constrained brute-force approach and umbrella sampling. The simulated free energies show excellent agreement with the experiment. The coordination structures of UO 2 (CO 3 ) 3 4- and CaUO 2 (CO 3 ) 3 2- were found to be significantly affected by the presence of Na + ions, leading to a monodentate binding of a carbonate group to U and a water molecule entering the first coordination shell of U. To shed light on the binding between uranyl and glutardiamidoxime (H 2 B), an important model and state-of-the-art ligand for seawater uranium extraction, a suite of computational methods was employed. From molecular dynamics (MD) simulations, it was found that the binding configurations of B 2- with uranyl favor the twofold distorted η 2 binding between the oximate ends (C=N-O-) and U, while HB - prefers the chelating mode for the oximate end with the neutral end being solvated by water. The free energies of sequential ligand binding to form UO 2 B, [UO 2 B 2 ] 2- , and [UO 2 (HB)B] - were simulated with umbrella sampling and very good agreement with the experimental values was achieved, which corroborates the structural insights into the binding mode. State-of-the-art polymeric sorbents employ both amidoximate and carboxylate groups on the side chains to achieve optimal U uptake and selectivity, so we simulated the binding of a model amidoximate–carboxylate bifunctional ligand with uranyl. Classical MD and free-energy simulations in 0.5 M NaCl showed that the carboxylate group binds first to uranyl, leading to a loose intermediate state, and then, the amidoximate group binds, resulting in a more stable and tight chelate state. Binding of the second bifunctional ligand follows a similar process, and the two ligands prefer a trans configuration around the uranyl group. The simulated free energies indicate that the two bifunctional ligands bind with uranyl 14 kcal/mol stronger than the two ligands with only amidoximate groups, confirming an important synergy between amidoximate and carboxylate groups in binding uranyl.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Long timescale solvation dynamics and confinement: The case of non-ionic deep eutectic solvents of lauric acid and N-methylacetamide

Microscopic segregation and molecular heterogeneities in complex liquids are the result of the interplay between different intermolecular forces, all of which contribute to the energy landscape of the system. A consequence of the intricate energy landscape is the nontrivial effect on the solvation dynamics. Here, the impact of molecular heterogeneities on the solvation dynamics is studied using infrared spectroscopies and molecular dynamics simulations. In particular, this study focuses on the dynamical effect of nanoscopic heterogeneities present in deep eutectic solvents (DESs) composed of lauric acid (LA) and N-methylacetamide (NMA). To this end, a molecular probe containing a carbon triple bond is used as an infrared reporter. The results show that the vibrational probe is likely to be located in the NMA polar domains. Furthermore, the probe solvation dynamics derived from the 2DIR spectra presents a slowdown of its timescale with increasing LA concentration in the DES. Kubo modeling of the probe solvation dynamics shows a correlation between the amplitude of its long time component and the presence of molecular heterogeneities in the sample. Semiclassical modeling of the vibrational line shape of the triple bond stretch demonstrates that the heterogeneities affect the whole solvation dynamics of the system through the amplitudes of the frequency fluctuations. Molecular dynamics simulations confirm the experimental results and their interpretation by showing a slowdown of the solvation dynamics when the LA heterogeneities are present. Overall, the study presents a molecular framework to explain the effect of confinement created by nanoscopic LA heterogeneities on the solvation dynamics of the system.

Chemical compounds↗
DOE OSTI2017

Hybrid Cascading Outage Analysis of Extreme Events with Optimized Corrective Actions

Power system are vulnerable to extreme contingencies (like an outage of a major generating substation) that can cause significant generation and load loss and can lead to further cascading outages of other transmission facilities and generators in the system. Some cascading outages are seen within minutes following a major contingency, which may not be captured exclusively using the dynamic simulation of the power system. The utilities plan for contingencies either based on dynamic or steady state analysis separately which may not accurately capture the impact of one process on the other. We address this gap in cascading outage analysis by developing Dynamic Contingency Analysis Tool (DCAT) that can analyze hybrid dynamic and steady state behavior of the power system, including protection system models in dynamic simulations, and simulating corrective actions in post-transient steady state conditions. One of the important implemented steady state processes is to mimic operator corrective actions to mitigate aggravated states caused by dynamic cascading. This paper presents an Optimal Power Flow (OPF) based formulation for selecting corrective actions that utility operators can take during major contingency and thus automate the hybrid dynamic-steady state cascading outage process. The improved DCAT framework with OPF based corrective actions is demonstrated on IEEE 300 bus test system.

Cascading outage↗
DOE OSTI2019

Direct measurement of nanostructural change during in situ deformation of a bulk metallic glass

Abstract To date, there has not yet been a direct observation of the initiation and propagation of individual defects in metallic glasses during deformation at the nanoscale. Here, we show through a combination of in situ nanobeam electron diffraction and large-scale molecular dynamics simulations that we can directly observe changes to the local short to medium range atomic ordering during the formation of a shear band. We observe experimentally a spatially resolved reduction of order prior to shear banding due to increased strain. We compare this to molecular dynamics simulations, in which a similar reduction in local order is seen, and caused by shear transformation zone activation, providing direct experimental evidence for this proposed nucleation mechanism for shear bands in amorphous solids. Our observation serves as a link between the atomistic molecular dynamics simulation and the bulk mechanical properties, providing insight into how one could increase ductility in glassy materials.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗
NASA NTRS2018

Simulating Forest Dynamics of Lowland Rainforests in Eastern Madagascar

Ecological modeling and forecasting are essential tools for the understanding of complex vegetation dynamics. The parametric demands of some of these models are often lacking or scant for threatened ecosystems, particularly in diverse tropical ecosystems. One such ecosystem and also one of the world's biodiversity hotspots, Madagascar's lowland rainforests, have disappeared at an alarming rate. The processes that drive tree species growth and distribution remain as poorly understood as the species themselves. We investigated the application of the process-based individual-based FORMIND model to successfully simulate a Madagascar lowland rainforest using previously collected multi-year forest inventory plot data. We inspected the model's ability to characterize growth and species abundance distributions over the study site, and then validated the model with an independently collected forest-inventory dataset from another lowland rainforest in eastern Madagascar. Following a comparative analysis using inventory data from the two study sites, we found that FORMIND accurately captures the structure and biomass of the study forest, with r(squared) values of 0.976, 0.895, and 0.995 for 1:1 lines comparing observed and simulated values across all plant functional types for aboveground biomass (tonnes/ha), stem numbers, and basal area (m(squared)/ha), respectively. Further, in validation with a second study forest site, FORMIND also compared well, only slightly over-estimating shade-intermediate species as compared to the study site, and slightly under-representing shade-tolerant species in percentage of total aboveground biomass. As an important application of the FORMIND model, we measured the net ecosystem exchange (NEE, in tons of carbon per hectare per year) for 50 ha of simulated forest over a 1000-year run from bare ground. We found that NEE values ranged between 1 and -1 t Cha(exp -1)year(exp -1), consequently the study forest can be considered as a net neutral or a very slight carbon sink ecosystem, after the initial 130 years of growth. Our study found that FORMIND represents a valuable tool toward simulating forest dynamics in the immensely diverse Madagascar rainforests.

Forest Productivity↗
DOE OSTI2025

Assessing the Effect of Explicit Polarizability on Models of Carbon Dioxide Solvation in Ionic Liquids

Ionic liquids are an important possible carbon capture material because of their anomalously high sorption selectivity for carbon dioxide over other gases common in air. Many research groups have investigated the molecular origins of this property and provided important insights, including using 1D and 2D-IR spectroscopy. Molecular dynamics simulations have been indispensable to the interpretation of these experiments. In prior molecular dynamics simulation work, charge-scaled force fields have typically been used to provide a mean-field treatment of effects vital to ionic liquid systems such as charge transfer and polarization. Here, we compare models of carbon dioxide solvated in ionic liquids with explicit polarization to models of the same with implicit polarizability through charge-scaling. We calculate structural, dynamical, and spectroscopic properties, and make comparisons to the same items measured in experiment. In this study, we focus on two ionic liquids: 1-butyl-3- methylimidazolium (BMIM + ) paired with bis(trifluoromethane sulfonyl imide) (Tf 2 N − ) and 1-butyl-3-methylimidazolium (BMIM+) paired with hexafluorophosphate (PF 6 − ). We find that many structural, dynamical, and spectroscopic properties are changed when polarization is modeled explicitly. We also find that explicit polarizability softens local ion cages around the carbon dioxide and that the long-time diffusion of the carbon dioxide is gated by the reorganization of the ionic liquid molecules. Comparisons to experiment show modest improvement of many observables compared with experiment for the explicitly polarizable model over the charge-scaled model. Overall, our results show that charge-scaled force fields are likely sufficient to compute spectroscopic properties of carbon dioxide in ionic liquids and suggest some interpretive rules for understanding their structural and dynamical properties. Those using charge-scaled force fields should generally assume that the ion cages around solutes such as carbon dioxide are too stiff and cation-rich in their models and adjust their interpretations and predictions accordingly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

System Identification Applied to Dynamic CFD Simulation and Wind Tunnel Data

Demanding aerodynamic modeling requirements for military and civilian aircraft have provided impetus for researchers to improve computational and experimental techniques. Model validation is a key component for these research endeavors so this study is an initial effort to extend conventional time history comparisons by comparing model parameter estimates and their standard errors using system identification methods. An aerodynamic model of an aircraft performing one-degree-of-freedom roll oscillatory motion about its body axes is developed. The model includes linear aerodynamics and deficiency function parameters characterizing an unsteady effect. For estimation of unknown parameters two techniques, harmonic analysis and two-step linear regression, were applied to roll-oscillatory wind tunnel data and to computational fluid dynamics (CFD) simulated data. The model used for this study is a highly swept wing unmanned aerial combat vehicle. Differences in response prediction, parameters estimates, and standard errors are compared and discussed

Murphy, Patrick C.↗