Characterization of Load Uncertainty on Simulated Dynamic Responses
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Coupled cluster theory is one of the most accurate electronic structure methods for predicting ground and excited state chemistry. However, the presence of numerical artifacts at electronic degeneracies, such as complex energies, has made it difficult to apply the method in nonadiabatic dynamics simulations. While it has already been shown that such numerical artifacts can be fully removed by using similarity constrained coupled cluster (SCC) theory [J. Phys. Chem. Lett. 2017, 8(19), 4801–4807], simulating dynamics requires efficient implementations of gradients and nonadiabatic couplings. Here, we present an implementation of nuclear gradients and nonadiabatic derivative couplings at the similarity constrained coupled cluster singles and doubles (SCCSD) level of theory, thereby making possible nonadiabatic dynamics simulations using a coupled cluster theory that provides a correct description of conical intersections between excited states. We present a few numerical examples that show good agreement with literature values and discuss some limitations of the method.
In this paper, time-dependent electronic structure methods are growing in popularity as tools for modeling ultrafast and/or nonlinear processes, for computing spectra, and as the electronic structure component of mean-field molecular dynamics simulations. Time-dependent configuration interaction (TD-CI) offers several advantages over the widely used real-time time-dependent density functional theory: namely, that it correctly models Rabi oscillations; it offers a spin-pure description of open-shell systems; and a hierarchy of TD-CI methods can be defined that systematically approach the exact solution of the time-dependent Schrodinger equation (TDSE). In this work, we present a novel TD-CI approach that extends TD-CI to large complete active-space configuration expansions. Such extension is enabled by use of a direct configuration interaction approach that eliminates the need to explicitly build, store, or diagonalize the Hamiltonian matrix. Graphics processing unit (GPU) acceleration enables fast solution of the TDSE even for large active spaces—up to 12 electrons in 12 orbitals (853776 determinants) in this work. A symplectic split operator propagator yields long-time norm conservation. We demonstrate the applicability of our approach by computing the response of a large molecule with a strongly correlated ground state, decacene (C 42 H 24 ), to various pulses (δ-function, transform limited, chirped). Our simulations predict that chirped pulses can be used to induce dipole-forbidden transitions. Simulations of decacene using the 6-31G(d) basis set and a 12 electrons/12 orbitals active space took 20.1 h to propagate for 100 fs with a 1 attosecond time step on a single NVIDIA K40 GPU. Convergence with respect to time step is found to depend on the property being computed and the chosen active space.
The objective of this multi-year project between the LLNL and Kyoto University, is the investigation of effects of rupture dynamics on near-fault ground motion amplitude for strike slip earthquakes through rupture dynamics simulations. The work presented in this report was in part performed by Jikai Sun, a graduate student from Kyoto university, during his visit to the LLNL as a summer student. He was trained by Dr. Pitarka in using computer programs for generating rupture models and performing rupture dynamics simulations.
Dynamic disturbances of many orders of magnitude greater than those of conventional spacecraft will be routine for the Space Station. Accurate knowledge of inflight structural dynamics and disturbances will be lacking. System identification will reduce uncertainties. To deal with the remaining model errors and time-varying elements, adaptive control may be required. Dynamic models for two Space Station configurations are derived. A direct model reference adaptive control algorithm is synthesized and evaluated with respect to the station model parameter errors and plant dynamics truncations. Both the rigid body and the flexible modes are treated. Simulation results show that convergence of the adaptive algorithm can be achieved in 100 to about 300 seconds with reasonable performance even during Shuttle hard docking operations in which station mass and inertia are instantaneously changed by more than 100 percent.
Power grid planning and operation decisions are made based on simulation of the dynamic behavior of the system. Enabling substantial energy savings while increasing the reliability of the aging North American power grid through improved utilization of existing transmission assets hinges on the adoption of wide-area measurement systems (WAMS) for power system stabilization. However, adoption of WAMS alone will not suffice if the power system is to reach its full entitlement in stability and reliability. It is necessary to enhance predictability with "faster than real-time" dynamic simulations that will enable the dynamic stability margins, proactive real-time control, and improve grid resiliency to fast time-scale phenomena such as cascading network failures. Present-day dynamic simulations are performed only during offline planning studies, considering only worst case conditions such as summer peak, winter peak days, etc. With widespread deployment of renewable generation, controllable loads, energy storage devices and plug-in hybrid electric vehicles expected in the near future and greater integration of cyber infrastructure (communications, computation and control), monitoring and controlling the dynamic performance of the grid in real-time would become increasingly important. The state-of-the-art dynamic simulation tools have limited computational speed and are not suitable for real-time applications, given the large set of contingency conditions to be evaluated. These tools are optimized for best performance of single-processor computers, but the simulation is still several times slower than real-time due to its computational complexity. With recent significant advances in numerical methods and computational hardware, the expectations have been rising towards more efficient and faster techniques to be implemented in power system simulators. This is a natural expectation, given that the core solution algorithms of most commercial simulators were developed decades ago, when High Performance Computing (HPC) resources were not commonly available.
Dynamic airfoil stall simulation in wind tunnels, considering pitch rate, Reynolds number, oscillation and test equipment effects
2-Propanol was investigated, in both the liquid and supercooled states, as a model system to study how hydrogen bonds affect the structural relaxation and the dynamics of mesoscale structures, of approximately several Ångstroms, employing static and quasi-elastic neutron scattering and molecular dynamics simulation. Dynamic neutron scattering measurements were performed over an exchanged wave-vector range encompassing the pre-peak, indicative of the presence of H-bonding associates, and the main peak. The dynamics observed at the pre-peak is associated with the formation and disaggregation of the H-bonded associates and is measured to be at least one order of magnitude slower than the dynamics at the main peak, which is identified as the structural relaxation. Additionally, the measurements indicate that the macroscopic shear viscosity has a similar temperature dependence as the dynamics of the H-bonded associates, which highlights the important role played by these structures, together with the structural relaxation, in defining the macroscopic rheological properties of the system. Importantly, the characteristic relaxation time at the pre-peak follows an Arrhenius temperature dependence whereas at the main peak it exhibits a non-Arrhenius behavior on approaching the supercooled state. The origin of this differing behavior is attributed to an increased structuring of the hydrophobic domains of 2-propanol accommodating a more and more encompassing H-bond network, and a consequent set in of dynamic cooperativity.
Subcutaneous injection is the preferred route of administration for many antibody therapeutics for reasons that include its speed and convenience. However, the small volume limit (typically ≤2 mL) for subcutaneous delivery often necessitates antibody formulations at high concentrations (commonly ≥100 mg/mL), which may lead to physicochemical problems. For example, antibodies with large hydrophobic or charged patches can be prone to self-interaction giving rise to high viscosity. Here, we combined X-ray crystallography with computational modeling to predict regions of an anti-glucagon receptor (GCGR) IgG 1 antibody prone to self-interaction. An extensive mutational analysis was undertaken of the complementarity-determining region residues residing in hydrophobic surface patches predicted by spatial aggregation propensity, in conjunction with residue-level solvent accessibility, averaged over conformational ensembles from molecular dynamics simulations. Dynamic light scattering (DLS) was used as a medium throughput screen for self-interaction of ~200 anti-GCGR IgG 1 variants. A negative correlation was found between the viscosity determined at high concentration (180 mg/mL) and the DLS interaction parameter measured at low concentration (2–10 mg/mL). Additionally, anti-GCGR variants were readily identified with reduced viscosity and antigen-binding affinity within a few fold of the parent antibody, with no identified impact on overall developability. The methods described here may be useful in the optimization of other antibodies to facilitate their therapeutic administration at high concentration.
Large-scale atomistic models for complex polycyclic aromatic hydrocarbon systems help understand the chemical properties and behaviors of complex feedstocks such as coal or petroleum. However, the development and utilization of large-scale models remain limited due to the difficulty in achieving the varied structural characteristics necessary to capture stochastic nature of these feedstocks. Here we demonstrate a systematic workflow to construct stochastic molecular systems from a broad analytical suite: high-resolution transmission electron microscopy (HRTEM), carbon-13 nuclear magnetic resonance spectroscopy ( 13 C NMR), laser desorption ionization mass spectroscopy (LDI-MS), and elemental analysis. We present a model construction and analysis utility of a new Python-based module. We selected one subbituminous and three high-volatile bituminous coals to construct large-scale models (~40,000 atoms). The constructed models were utilized to examine the affinity for solvent extraction (naphthalene or tetralin) and the effect of structural properties (e.g., aromatic cluster size, functional groups, and cross-linking) in reactive molecular dynamics simulations. Complex chemical reactions were monitored with bond order transitions, intermediates formation, and mass distributions. Reactive molecular dynamics simulations suggest a plausible chemical extraction process and products for the complex fossil feedstocks. The results indicated that radical formations with bond breaking of bridging oxygens and carbons were required at high temperatures to facilitate hydrogeneration and extraction of gas molecules from radical-free molecules. We observed that aliphatic chains of tetralin were easily decomposed and combined with radicals to form small size of molecules with aryl bonding, mainly increasing molecules in the 500–1000 Da, while naphthalene had little impact on chemical extraction process.
Carbon fiber-reinforced polymer (CFRP) composites have shown promise as a material for structures designed to absorb energy in crush-style impact loading. In complementing the lightweight and tailorable characteristics of CFRP structures, the benefits of impact energy absorption are particularly interesting in aerospace vehicles that require excellent crashworthiness performance. However, simulating the behavior of CFRP structures in dynamic, crushing-style impact loading is challenging because of the many damage and failure modes that are essential to capture in the model. For the present investigation, CFRP stanchions were tested using a crash sled experimental test rig. The stanchions were designed by the Composite Materials Handbook-17 (CMH-17) Crashworthiness Working Group for the purpose of comparing experimental crash sled tests to corresponding simulations of the tests. The stanchions are C-channel shaped and represent a geometry common in the interior of aerospace vehicle structures. Explicit simulations in LSDYNA were performed using a well-established composite material model (MAT58) and a next-generation material model (MAT213). Fully integrated shell elements were utilized instead of 3D solid elements to limit computation time. In all models, each of the sixteen plies were represented with individual layers of shell elements with tiebreak contact between each adjacent ply pair to simulate interlaminar fracture. Simulating the crushing of the stanchions occurred in two phases. First, the material models were calibrated using flat specimens that were manufactured with the same layup as the stanchions. While the original goal was calibration of material-related properties, meshdependent behavior was observed in simulations with either material model, and an unstructured mesh was selected to remediate undesirable mesh-dependent failure modes. Additionally, for both the MAT58 and MAT213 models, it was found that either the crush force or the failure mode could be modeled accurately, but no set of parameters could be identified to attain both results in the same model. Once satisfactory calibration was achieved, the same material parameters were applied to the stanchion crush simulations. The stanchion simulations showed that MAT213 more accurately predicted the experimentally determined crush force, and both material models predicted key aspects of the experimentally observed failure modes.
A software process is developed to convert distribution network models from a quasi-static time-series tool (OpenDSS) to a real-time dynamic phasor simulator (ePHASORSIM). The description of this process in this paper would be helpful for researchers who intend to perform similar conversions. The converter could be utilized directly by users of real-time simulators who intend to perform software-in-the-loop or hardware-in-the-loop tests on large distribution test feeders for a range of use cases, including testing functions of advanced distribution management systems against a simulated distribution system. In the future, the developers intend to release the conversion tool as open source to enable use by others.
Multizone models are widely used in building airflow and energy performance simulations due to their fast computing speed. However, multizone models assume that the air in a room is well mixed, consequently limiting their application. In specific rooms where this assumption fails, the use of computational fluid dynamics (CFD) models may be an alternative option. Previous research has mainly focused on coupling CFD models and multizone models to study airflow in large spaces. While significant, most of these analyses did not consider the coupled simulation of the building airflow with the building’s Heating, Ventilation, and Air-Conditioning (HVAC) systems. This paper tries to fill the gap by integrating the models for HVAC systems with coupled multizone and CFD simulations for airflows, using the Modelica simulation platform. To improve the computational efficiency, we incorporated a simplified CFD model named fast fluid dynamics (FFD). We first introduce the data synchronization strategy and implementation in Modelica. Then, we verify the implementation using two case studies involving an isothermal and a non-isothermal flow by comparing model simulations to experiment data. Afterward, we study another three cases that are deemed more realistic. Furthermore, this is done by attaching a variable air volume (VAV) terminal box and a VAV system to previous flows to assess the capability of the models in studying the dynamic control of HVAC systems. Finally, we discuss further research needs on the coupled simulation using the models.
Abstract not provided
The efficacy of vegetation dynamics simulations in offline land surface models (LSMs) largely depends on the quality and spatial resolution of meteorological forcing data. In this study, the Princeton Global Meteorological Forcing Data (PMFD) and the high spatial resolution and upscaled China Meteorological Forcing Data (CMFD) were used to drive the Simplified Simple Biosphere model version 4/Top-down Representation of Interactive Foliage and Flora Including Dynamics (SSiB4/TRIFFID) and investigate how meteorological forcing datasets with different spatial resolutions affect simulations over the Tibetan Plateau (TP), a region with complex topography and sparse observations. By comparing the monthly Leaf Area Index (LAI) and Gross Primary Production (GPP) against observations, we found that SSiB4/TRIFFID driven by upscaled CMFD improved the performance in simulating the spatial distributions of LAI and GPP over the TP, reducing RMSEs by 24.3% and 20.5%, respectively. The multi-year averaged GPP decreased from 364.68 gC m –2 yr –1 to 241.21 gC m –2 yr –1 with the percentage bias dropping from 50.2% to –1.7%. When using the high spatial resolution CMFD, the RMSEs of the spatial distributions of LAI and GPP simulations were further reduced by 7.5% and 9.5%, respectively. This study highlights the importance of more realistic and high-resolution forcing data in simulating vegetation growth and carbon exchange between the atmosphere and biosphere over the TP.