Simulation of Rough Electrodes Coupled with Structural Dynamics
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.
Inverters using phase-locked loops for control depend on voltages generated by synchronous machines to operate. This might be problematic if much of the conventional generation fleet is displaced by inverters. To solve this problem, grid-forming control for inverters has been proposed as being capable of autonomously regulating grid voltages and frequency. Presently, the performance of bulk power systems with massive penetration of grid-forming inverters has not been thoroughly studied as to elucidate benefits. Hence, this paper presents inverter models with two grid-forming strategies: virtual oscillator control and droop control. The two models are specifically developed to be used in positive-sequence simulation packages and have been implemented in PSLF. The implementations are used to study the performance of bulk power grids incorporating inverters with gridforming capability. Specifically, simulations are conducted on a modified IEEE 39-bus test system and the microWECC test system with varying levels of synchronous and inverter-based generation. The dynamic performance of the tested systems with gridforming inverters during contingency events is better than cases with only synchronous generation.
Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.
The recently developed real-time nuclear-electronic orbital (RT-NEO) approach provides an elegant framework for treating electrons and selected nuclei, typically pro- tons, quantum mechanically in nonequilibrium dynamical processes. However, the RT-NEO approach neglects the motion of the other nuclei, preventing a complete de- scription of the coupled nuclear-electronic dynamics and spectroscopy. In this work, we describe the dynamical interaction between the other nuclei and the electron-proton subsystem with the mixed quantum-classical Ehrenfest dynamics method. The NEO- Ehrenfest approach propagates the electrons and quantum protons in a time-dependent variational framework, while the remaining nuclei move classically on the correspond- ing average electron-proton vibronic surface. This approach includes the non-Born- Oppenheimer eects between the electrons and quantum protons with RT-NEO and between the classical nuclei and the electron-proton subsystem with Ehrenfest dynam- ics. Spectral features for vibrational modes involving both quantum and classical nuclei are resolved from the time-dependent dipole moments. Moreover, nuclear quantum ef- fects are included directly in nonadiabatic dynamics simulations. The excited state intramolecular proton transfer in the o-hydroxybenzaldehyde molecule is shown to be faster and to exhibit a larger kinetic isotope eect with NEO-Ehrenfest compared to classical Ehrenfest dynamics. This work shows that the NEO-Ehrenfest method is a powerful tool to study dynamical processes with coupled electronic and nuclear degrees of freedom.
Complexities of flow patterns in the azimuthal cross-section of a cylindrical magnetized helicon plasma and the corresponding plasma dynamics are investigated by means of a novel scheme for time delay estimation velocimetry. The advantage of this introduced method is the capability of calculating the time-averaged 2D velocity fields of propagating wave-like structures and patterns in complex spatiotemporal data. It is able to distinguish and visualize the details of simultaneously present superimposed entangled dynamics and it can be applied to fluid-like systems exhibiting frequently repeating patterns (e.g., waves in plasmas, waves in fluids, dynamics in planetary atmospheres, etc.). The velocity calculations are based on time delay estimation obtained from cross-phase analysis of time series. Each velocity vector is unambiguously calculated from three time series measured at three different non-collinear spatial points. This method, when applied to fast imaging, has been crucial to understand the rich plasma dynamics in the azimuthal cross-section of a cylindrical linear magnetized helicon plasma. The capabilities and the limitations of this velocimetry method are discussed and demonstrated for two completely different plasma regimes, i.e., for quasi-coherent wave dynamics and for complex broadband wave dynamics involving simultaneously present multiple instabilities.
Water’s extended hydrogen-bond network results in rich and complex dynamics on the sub-picosecond time scale. In this paper, we present a comprehensive analysis of the two-dimensional infrared (2D IR) spectrum of O–H stretching vibrations in liquid H{sub 2}O and their interactions with bending and intermolecular vibrations. By exploring the dependence of the spectrum on waiting time, temperature, and laser polarization, we refine our molecular picture of water’s complex ultrafast dynamics. The spectral evolution following excitation of the O–H stretching resonance reveals vibrational dynamics on the 50–300 fs time scale that are dominated by intermolecular delocalization. These O–H stretch excitons are a result of the anharmonicity of the nuclear potential energy surface that arises from the hydrogen-bonding interaction. The extent of O–H stretching excitons is characterized through 2D depolarization measurements that show spectrally dependent delocalization in agreement with theoretical predictions. Furthermore, we show that these dynamics are insensitive to temperature, indicating that the exciton dynamics alone set the important time scales in the system. Finally, we study the evolution of the O–H stretching mode, which shows highly non-adiabatic dynamics suggestive of vibrational conical intersections. We argue that the so-called heating, commonly observed within ∼1 ps in nonlinear IR spectroscopy of water, is a nonequilibrium state better described by a kinetic temperature rather than a Boltzmann distribution. Our conclusions imply that the collective nature of water vibrations should be considered in describing aqueous solvation.
Abstract Photorespiration is a dynamic process that is intimately linked to photosynthetic carbon assimilation. There is a growing interest in understanding carbon assimilation during dynamic conditions, but the role of photorespiration under such conditions is unclear. In this review, we discuss recent work relevant to the function of photorespiration under dynamic conditions, with a special focus on light transients. This work reveals that photorespiration is a fundamental component of the light induction of assimilation where variable diffusive processes limit CO2 exchange with the atmosphere. Additionally, metabolic interactions between photorespiration and the C3 cycle may help balance fluxes under dynamic light conditions. We further discuss how the energy demands of photorespiration present special challenges to energy balancing during dynamic conditions. We finish the review with an overview of why regulation of photorespiration may be important under dynamic conditions to maintain appropriate fluxes through metabolic pathways related to photorespiration such as nitrogen and one-carbon metabolism.
This paper presents a modeling approach for dynamic aeroelastic flight dynamic analysis of the Mach0.745 Transonic Truss-Braced Wing. The modeling approach is based on a transonic correction method to correct the Theodorsen’s theory for transonic flow. CFD unsteady Reynolds-averaged Navier-Stokes equations (RANS) simulations are conducted using FUN3D solver for a series of wing sections from the Mach 0.745Transonic Truss-Braced Wing in pitch and plunge oscillations. Unsteady lift and pitching moment coefficients are obtained and used to develop the correction terms in the Theodorsen’s theory to account for transonic aerodynamics. The unsteady lift and pitching moment derivatives with respect to the unsteady angle of attack are obtained as functions of the reduced frequency. These derivatives are used to compute the unsteady lift and pitching moment contributions by the angle of attack, pitch rate, and roll rate. They are then approximated using a frequency-domain regression to obtain the dynamic stability derivatives for the Mach 0.745 Transonic Truss-Braced Wing. The structural dynamic mode shapes of the Mach 0.745 TTBW are extracted from a NAS-TRAN finite-element model. These mode shapes are used to compute the generalized unsteady aerodynamic forces. The aerodynamic mass, damping, and stiffness and the aerodynamic lag states are constructed to couple the dynamic aeroelastic contribution to the flight dynamic model of the Mach 0.745 Transonic Truss-Braced Wing. The coupled dynamic aeroelastic flight dynamic equations of motion are formulated. The eigenvalues of the coupled system are computed. All the flight dynamic modes and structural dynamic modes are stable at Match 0.745. The effect of transonic aerodynamics generally causes all the dynamic modes to have lower damping values
Background: SPECT myocardial perfusion imaging (MPI) is a clinical mainstay that is typically performed with static imaging protocols and visually or semi-quantitatively assessed for perfusion defects based upon the relative intensity of myocardial regions. Dynamic cardiac SPECT presents a new imaging technique based on time-varying information of radiotracer distribution, which permits the evaluation of regional myocardial blood flow (MBF) and coronary flow reserve (CFR). In this work, a preliminary feasibility study was conducted in a small patient sample designed to implement a unique combined static-dynamic single-dose one-day visit imaging protocol to compare quantitative dynamic SPECT with static conventional SPECT for improving the diagnosis of coronary artery disease (CAD). Methods: Fifteen patients (11 males, four females, mean age 71 ± 9 years) were enrolled for a combined dynamic and static SPECT (Infinia Hawkeye 4, GE Healthcare) imaging protocol with a single dose of {sup 99m}Tc-tetrofosmin administered at rest and a single dose administered at stress in a one-day visit. Out of 15 patients, eleven had selective coronary angiography (SCA), 8 within 6 months and the rest within 24 months of SPECT imaging, without intervening symptoms or interventions. The extent and severity of perfusion defects in each myocardial region was graded visually. Dynamically acquired data were also used to estimate the MBF and CFR. Both visually graded images and estimated CFR were tested against SCA as a reference to evaluate the validity of the methods. Results: Overall, conventional static SPECT was normal in ten patients and abnormal in five patients, dynamic SPECT was normal in 12 patients and abnormal in three patients, and CFR from dynamic SPECT was normal in nine patients and abnormal in six patients. Among those 11 patients with SCA, conventional SPECT was normal in 5, 3 with documented CAD on SCA with an overall accuracy of 64%, sensitivity of 40% and specificity of 83%. Dynamic SPECT image analysis also produced a similar accuracy, sensitivity, and specificity. CFR was normal in 6, each with CAD on SCA with an overall accuracy of 91%, sensitivity of 80%, and specificity of 100%. The mean CFR was significantly lower for SCA detected abnormal than for normal patients (3.86±1.06 vs 1.94±0. 0.67, P < 0.001). Conclusions: The visually assessed image findings in static and dynamic SPECT are subjective, and may not reflect direct physiologic measures of coronary lesion based on SCA. The CFR measured with dynamic SPECT is fully objective, with better sensitivity and specificity, available only with the data generated from the dynamic SPECT method.
As a compelling case of sensitive structure–property relationship, an odd–even effect refers to the alternating trend of physical or chemical properties on odd/even number of repeating structural units. In crystalline or semicrystalline materials, such odd–even effects emerge as manifestations of differences in the periodic packing patterns of molecules. Therefore, due to the lack of long-range order, such an odd–even phenomenon is not expected for dynamic properties in amorphous state. Herein, we report the discovery of a remarkable odd–even effect of dynamical properties in the liquid phase. In a class of glass-forming diammonium citrate ionic liquids, using incoherent quasi-elastic neutron scattering measurements, we measured the dynamical properties including diffusion coefficient and rotational relaxation time. These directly measured molecular dynamics showed pronounced alternating trends with increased number of methylene (-CH2-) groups in the backbone. Meanwhile, the structure factor S(Q) showed no long-range periodic packing of molecules, while the pair distribution function G(r) revealed subtle differences in the local molecular morphology. The observed dynamical odd–even phenomenon in liquids showed that profound dynamical changes originate from subtle local structural differences.
Coupling of sloshing dynamics within a partially filled rotating dewar of superfluid helium 2 with spacecraft dynamics are investigated in response to the realistic environmental disturbance forces and torques acting on the spacecraft during normal operation. This study investigates: (1) the rotating bubble of superfluid helium 2 reacting to combined environmental disturbances, including gravity gradient, aerodynamic, and magnetic forces and torques; (2) characteristics of slosh reaction forces and torques coupling with spacecraft dynamics; (3) the contribution of slosh dynamics to over-all spacecraft dynamics; and (4) activating of attitude and translation control system. The numerical computation of sloshing dynamics is based on the rotational frame, while the spacecraft dynamics is associated with non-rotational frame. Results show that the contributions of spacecraft dynamics are driven by the environmental disturbances coupling with slosh dynamics. Without considering the effects of environmental disturbances-driven slosh dynamics acting on spacecraft coupling with the spacecraft dynamics may lead to the wrong results for the development of spacecraft system guidance and attitude control techniques.
We consider a problem of non-adiabatic dynamics of a 2D fermionic system with d + id-wave symmetry of pairing amplitude. Under the mean-field approximation, we determine the asymptotic behavior of the pairing amplitude following a sudden change of coupling strength. We also study an extended d + id pairing system for which the long-time asymptotic states of the pairing amplitude in the collisionless regime can be determined exactly. By using numerical methods, we have identified three non-equilibrium steady states described by different long-time asymptotes of the pairing amplitude for both the non-integrable and the integrable versions of d + id-wave models. We found that despite of its lack of integrability, long-time dynamics resulting from pairing quenches in the non-integrable d + id model are essentially similar to the ones found for its exactly integrable extended d + id model. We also obtain the long-time phase diagram of the extended d + id model through the Lax construction that exploits underlying integrability showing that the dynamic phases obtained by numerics are consistent with the dynamics of the exactly integrable approach. Both models describe a topological fermionic system with a topologically non-trivial BCS phase appearing at weak coupling strength. We show that the presence of oscillating order parameter region in the chiral d + id pairing dynamics differs from the d-wave (d{sub x{sup 2}−y{sup 2}}), which may be used to probe pairing symmetries of chiral superconductors.
Vitrimers are a class of polymer networks featuring dynamic covalent cross-links that can undergo associative bond exchange. There has been recent interest in these materials due to their promise as recyclable thermosets or self-healing polymers because of the ability of vitrimer networks to rearrange at the molecular level and undergo macroscopic flow. However, the practical use of these materials often occurs in the supercooled regime or glassy state, where the implications of dynamic bonds are complicated by the interplay between slow activated segmental dynamics, cross-link (i.e., bond-exchange) kinetics, and ultimately material properties. Here, in this paper, we combine coarse-grained molecular dynamics simulation and microscopic statistical mechanical theory to understand how cross-linking kinetics affect material dynamics and how this couples to segmental relaxation of the polymeric network strands across a spectrum of length and time scales, especially in the supercooled regime. We characterize the Kuhn segmental alpha relaxation time and bond exchange time for vitrimer systems across various cross-link densities, temperatures, and bond exchange rates. Simulation and theory both exhibit a bending-up behavior for bond exchange time upon cooling, suggesting a coupling between bond exchange dynamics and segmental relaxation that intensifies with faster bond exchange kinetics. We also found bond exchange dynamics have an impact on Kuhn segment alpha relaxation time, which is most significant at higher cross-link densities. Both these effects are most prominent when the bond exchange time is similar to the Kuhn segment alpha relaxation time, and the resulting coupling of these two relaxation processes is tied to both the probability of a free end to find a bonded pair and the time scale of the constraints imposed by the dynamic cross-links. This relationship is reflected by a cross-link dependence of a theoretical parameter which represents the quantitative degree of coupling between bond exchange and segmental dynamics. Overall, the combination of simulation and theory clarifies the intricate interaction between bond kinetics and segmental relaxation and demonstrates the ability to provide molecular-level insights into vitrimer dynamics over a wide temperature range.
We derive a coarse-grained description of the dynamics of a nanoparticle immersed in an isothermal simple fluid by performing a systematic coarse graining of the underlying microscopic dynamics. As coarse-grained or relevant variables, we select the position of the nanoparticle and the total mass and momentum density field of the fluid, which are locally conserved slow variables because they are defined to include the contribution of the nanoparticle. The theory of coarse graining based on the Zwanzing projection operator leads us to a system of stochastic ordinary differential equations that are closed in the relevant variables. We demonstrate that our discrete coarse-grained equations are consistent with a Petrov-Galerkin finite-element discretization of a system of formal stochastic partial differential equations which resemble previously used phenomenological models based on fluctuating hydrodynamics. Key to this connection between our “bottom-up” and previous “top-down” approaches is the use of the same dual orthogonal set of linear basis functions familiar from finite element methods (FEMs), both as a way to coarse-grain the microscopic degrees of freedom and as a way to discretize the equations of fluctuating hydrodynamics. Another key ingredient is the use of a “linear for spiky” weak approximation which replaces microscopic “fields” with a linear FE interpolant inside expectation values. For the irreversible or dissipative dynamics, we approximate the constrained Green-Kubo expressions for the dissipation coefficients with their equilibrium averages. Under suitable approximations, we obtain closed approximations of the coarse-grained dynamics in a manner which gives them a clear physical interpretation and provides explicit microscopic expressions for all of the coefficients appearing in the closure. Our work leads to a model for dilute nanocolloidal suspensions that can be simulated effectively using feasibly short molecular dynamics simulations as input to a FEM fluctuating hydrodynamic solver.
Ferroelectric materials such as BaTiO 3 show tremendous potential for emerging advances in memory devices, particular neuromorphic type devices. High density of memory can be obtained by stabilising polar domain walls at the nanoscale, regions of discontinuity between the well-defined polarization order parameter, but little is known about what controls their structure and dynamics in real nanoscale materials. Indeed, chiral polar domain walls have been observed in heterogeneous ferroelectrics, such as oxygen-deficient BaTiO 3 , but very little is known about how such polar-domains walls interact with defects. Indeed, a critical understanding of how dynamics of domain-walls depend on point-defects is crucial to create engineered ferroelectric memory devices. For this work, we perform large-scale simulations of nansocale domain-wall dynamics in pristine and defective BaTiO 3 using reactive force-field developed by us earlier (Phys. Chem. Chem. Phys., 2019, 21, 18240–18249), and capture their dynamical dependence on point defects using a graph dynamical neural-network approach, which we adapted to interrogate solids with well-defined order-parameters, and implemented using Pytorch based libraries. Our machine learning (ML) approach goes beyond the traditional post-processing methods to capture both spatial and temporal heterogeneities of large-scale molecular dynamics simulations of complex defective ferroelectric oxide materials. We crucially find that isolated oxygen vacancies introduce very localized spatial regions (~1–2 unit-cell in length) that show slow dipole relaxation due to formation of defect-dipoles, and that these defect-dipoles in turn slow the intrinsic dynamics of domain walls. Further, the roughness of domain walls, also influenced by vacancies, introduce dynamic heterogeneity along the domain-wall. As such we find a novel mechanism by which quenched disorder due to defects introduce dynamic heterogeneity thereby influencing response to external fields (particularly time varying fields) in a ferroelectric. Our study also emphasizes the need for creating digital twins of dynamical quantities to achieve autonomous in operando control of nanoscale switching.
The formation of stable radiation damage in crystalline solids often proceeds via complex dynamic annealing processes, involving migration and interaction of ballistically-generated point defects. The dominant dynamic annealing processes, however, remain unknown even for crystalline Si. Here, we use a pulsed ion beam method to study defect dynamics in Si bombarded in the temperature range from -20 to 140 °C with 500 keV Ar ions. Results reveal a defect relaxation time constant of ~10–0.2 ms, which decreases monotonically with increasing temperature. The dynamic annealing rate shows an Arrhenius dependence with two well-defined activation energies of 73 ± 5 meV and 420 ± 10 meV, below and above 60 °C, respectively. Rate theory modeling, bench-marked against this data, suggests a crucial role of both vacancy and interstitial diffusion, with the dynamic annealing rate limited by the migration and interaction of vacancies.
A six-degree-of-freedom nonlinear semi-definite model with time varying mesh stiffness has been developed for the dynamic analysis of spur gears. The model includes a spur gear pair, two shafts, two inertias representing load and prime mover, and bearings. As the shaft and bearing dynamics have also been considered in the model, the effect of lateral-torsional vibration coupling on the dynamics of gears can be studied. In the nonlinear model developed several factors such as time varying mesh stiffness and damping, separation of teeth, backlash, single- and double-sided impacts, various gear errors and profile modifications have been considered. The dynamic response to internal excitation has been calculated by using the 'static transmission error method' developed. The software prepared (DYTEM) employs the digital simulation technique for the solution, and is capable of calculating dynamic tooth and mesh forces, dynamic factors for pinion and gear, dynamic transmission error, dynamic bearing forces and torsions of shafts. Numerical examples are given in order to demonstrate the effect of shaft and bearing dynamics on gear dynamics.
Microgrids have repeatedly demonstrated the ability to provide uninterrupted service to critical end-use loads during normal outages, severe weather events, and natural disasters. While their ability to provide critical services is well documented, microgrids present a more dynamic operational environment than grid-connected distribution systems. The dynamics of a microgrid are commonly driven by the high inertia of rotating generators, which are common in many microgrids. In such high-inertia systems, the impact of end-use load dynamics are often not examined. However, with the increased penetration of inverter-based generation with little or no inertia, it is necessary to consider the impact that the dynamics of the end-use loads have on the operations of microgrids, particularly for a resiliency-based operation. These operations include, but are not limited to, switching operations, loss of generating units, and the starting of induction motors. This paper examines the importance of including dynamic models of the end-use load when evaluating the operations of low inertia microgrids, and shows that by properly representing their behavior, it is possible to cost effectively size equipment while supporting resilient operations of critical end-use loads.