Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamical System Scaling”

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 649 records · Page 36

Identification of the weak-to-strong transition in Alfvénic turbulence from space plasma

Abstract Plasma turbulence is a ubiquitous dynamical process that transfers energy across many spatial and temporal scales in astrophysical and space plasma systems. Although the theory of anisotropic magnetohydrodynamic (MHD) turbulence has successfully described natural phenomena, its core prediction of an Alfvénic transition from weak to strong MHD turbulence when energy cascades from large to small scales has not been observationally confirmed. Here we report evidence for the Alfvénic weak-to-strong transition in small-amplitude, turbulent MHD fluctuations in Earth’s magnetosheath using data from four Cluster spacecraft. Our observations demonstrate the universal existence of strong turbulence accompanied by weak turbulent fluctuations on large scales. Moreover, we find that the nonlinear interactions of MHD turbulence are crucial to the energy cascade, as they broaden the cascade directions and fluctuating frequencies. The observed connection between weak and strong MHD turbulence systems may be present in many astrophysical environments, such as star formation, energetic particle transport, turbulent dynamos, and solar corona or solar wind heating.

79 ASTRONOMY AND ASTROPHYSICS↗

A multisite decomposition of the tensor network path integrals

Tensor network decompositions of path integrals for simulating open quantum systems have recently been proven to be useful. However, these methods scale exponentially with the system size. This makes it challenging to simulate the non-equilibrium dynamics of extended quantum systems coupled with local dissipative environments. In this work, we extend the tensor network path integral (TNPI) framework to efficiently simulate such extended systems. The Feynman–Vernon influence functional is a popular approach used to account for the effect of environments on the dynamics of the system. In order to facilitate the incorporation of the influence functional into a multisite framework (MS-TNPI), we combine a matrix product state (MPS) decomposition of the reduced density tensor of the system along the sites with a corresponding tensor network representation of the time axis to construct an efficient 2D tensor network. The 2D MS-TNPI network, when contracted, yields the time-dependent reduced density tensor of the extended system as an MPS. The algorithm presented is independent of the system Hamiltonian. We outline an iteration scheme to take the simulation beyond the non-Markovian memory introduced by solvents. Applications to spin chains coupled to local harmonic baths are presented; we consider the Ising, XXZ, and Heisenberg models, demonstrating that the presence of local environments can often dissipate the entanglement between the sites. We discuss three factors causing the system to transition from a coherent oscillatory dynamics to a fully incoherent dynamics. The MS-TNPI method is useful for studying a variety of extended quantum systems coupled with solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Large-scale Multiphysics Simulations of Small Modular Reactors Operating in Natural Circulation

Thanks to the advancements in high-performance computing, advanced modeling and simulation have become crucial in driving the development and deployment of next-generation nuclear reactors, such as small modular reactors (SMRs). SMRs offer the promise of cost-effective baseload electricity production and improved safety, while addressing some of the challenges associated with large reactor designs, such as high capital costs and extended construction timelines. As part of the Exascale Computing Project, the large-scale multiphysics simulation of an entire SMR primary system has been achieved by combining computational fluid dynamics and neutronics. In addition to the successful demonstration of full-core SMR simulations, the current study integrated the impact of natural circulation into the system. Natural circulation is the primary mechanism driving coolant circulation in SMRs. The mass flow rate in the core depends on the core power, and a numerical model has been developed to predict it. The pressure drop caused by the helical coil steam generator was also accounted for by developing a pressure drop correlation based on high-fidelity large eddy simulation results, further improving prediction accuracy. In conclusion, the results of the study demonstrate that the implemented natural circulation model is effective in predicting the responses of SMR full-core multiphysics simulations.

ECP↗

Stochastic density functional theory combined with Langevin dynamics for warm dense matter

Here, this study overviews and extends a recently developed stochastic finite-temperature Kohn-Sham density functional theory to study warm dense matter using Langevin dynamics, specifically under periodic boundary conditions. The method's algorithmic complexity exhibits nearly linear scaling with system size and is inversely proportional to the temperature. Additionally, a linear-scaling stochastic approach is introduced to assess the Kubo-Greenwood conductivity, demonstrating exceptional stability for dc conductivity. Utilizing the developed tools, we investigate the equation of state, radial distribution, and electronic conductivity of hydrogen at a temperature of 30 000 K. As for the radial distribution functions, we reveal a transition of hydrogen from gaslike to liquidlike behavior as its density exceeds 4 g/cm 3 . As for the electronic conductivity as a function of the density, we identified a remarkable isosbestic point at frequencies around 7 eV, which may be an additional signature of a gas-liquid transition in hydrogen at 30 000 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identification of the weak-to-strong transition in Alfvénic turbulence from space plasma

Plasma turbulence is a ubiquitous dynamical process that transfers energy across many spatial and temporal scales in astrophysical and space plasma systems. Although the theory of anisotropic magnetohydrodynamic (MHD) turbulence has successfully described natural phenomena, its core prediction of an Alfvénic transition from weak to strong MHD turbulence when energy cascades from large to small scales has not been observationally confirmed. Here we report evidence for the Alfvénic weak-to-strong transition in small-amplitude, turbulent MHD fluctuations in Earth’s magnetosheath using data from four Cluster spacecraft. Our observations demonstrate the universal existence of strong turbulence accompanied by weak turbulent fluctuations on large scales. Moreover, we find that the nonlinear interactions of MHD turbulence are crucial to the energy cascade, as they broaden the cascade directions and fluctuating frequencies. The observed connection between weak and strong MHD turbulence systems may be present in many astrophysical environments, such as star formation, energetic particle transport, turbulent dynamos, and solar corona or solar wind heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling and Power-Hardware-in-the-Loop Validation of Synchronous Machine Governor

This paper introduces the development of a high-fidelity gas turbine governor model using a programmable logic controller for power-hardware-in-the-loop (PHIL) validation. The governor model is integrated with the National Renewable Energy Laboratory's (NREL) PHIL test bed, featuring a 2-MVA synchronous machine and a 2.5-MW variable-speed drive, to emulate NG-driven HRSGs and CTs under various operational scenarios. The primary objective of this research is to study the grid-connected and islanding operations of conventional generation sources, with representative startup sequences including turbine purge, ignition, speed ramp-up, synchronization, and breaker closure. Preliminary results of the generator governor model on NREL PHIL platform, particularly using the 2.5-MW dynamometer system, offered significant insights into the modeling techniques, hardware integration, scaling, and real-world simulation dynamics.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Visual/motion cue mismatch in a coordinated roll maneuver

The effects of bandwidth differences between visual and motion cueing systems on pilot performance for a coordinated roll task were investigated. Visual and motion cue configurations which were acceptable and the effects of reduced motion cue scaling on pilot performance were studied to determine the scale reduction threshold for which pilot performance was significantly different from full scale pilot performance. It is concluded that: (1) the presence or absence of high frequency error information in the visual and/or motion display systems significantly affects pilot performance; and (2) the attenuation of motion scaling while maintaining other display dynamic characteristics constant, affects pilot performance.

Shirachi, D. K.↗

Component Cost Analysis of Large Scale Systems

The ideas of cost decomposition is summarized to aid in the determination of the relative cost (or 'price') of each component of a linear dynamic system using quadratic performance criteria. In addition to the insights into system behavior that are afforded by such a component cost analysis CCA, these CCA ideas naturally lead to a theory for cost-equivalent realizations.

Skelton, R. E.↗

Modal techniques for analyzing airplane dynamics

A series of techniques are presented for analyzing airplane control and output characteristics. State-space matrix equations describing the linear perturbation dynamics are transformed from physical coordinates into scaled, modal coordinates. Techniques for analyzing the impacts of system inputs on the fundamental modes of motion and the appearance of these modes in the system outputs are explained. Scaled modal matrices are used to compute steady-state control inputs which optimize the steady-state response of selected system outputs. Graphics which promote quick understanding of the analysis are presented to display the resulting vectors and matrices. The defined analysis techniques are applied to an example airplane model, illustrating the insight which can be acquired using the described modal techniques.

Arbuckle, P. Douglas↗

Magnetic resonance imaging of convection in laser-polarized xenon

We demonstrate nuclear magnetic resonance (NMR) imaging of the flow and diffusion of laser-polarized xenon (129Xe) gas undergoing convection above evaporating laser-polarized liquid xenon. The large xenon NMR signal provided by the laser-polarization technique allows more rapid imaging than one can achieve with thermally polarized gas-liquid systems, permitting shorter time-scale events such as rapid gas flow and gas-liquid dynamics to be observed. Two-dimensional velocity-encoded imaging shows convective gas flow above the evaporating liquid xenon, and also permits the measurement of enhanced gas diffusion near regions of large velocity variation.

Non-NASA Center↗

Mesoscale Simulations of CRYSTAL-FACE 23 July 2002 Case

A key objective of the Cirrus Regional Study of Tropical Anvils and Cirrus Layers - Florida Area Cirrus Experiment (CRYSTAL-FACE) is to understand the relationships between properties of tropical convective cloud systems and the lifecycle of the extended cirrus anvils they produce. We report here on a case study of 23 July 2002 where a line of land-based convective storms was generated between Lake Okeechobee and the Florida east coast as a result of complex interactions between lake and sea breeze fronts and outflow boundaries. A central goal of this study is to develop a description of convective input to the anvil system and to quantify the ongoing dynamical forcing of anvil processes by mesoscale and large-scale dynamics. This information is then used to force high-resolution cloud simulations with a model that explicitly resolves cloud microphysical processes (bin model) for study of cirrus anvil microphysical development.

Starr, David↗

Satellite Soil Moisture Data Assimilation Impacts on Modeling Weather Variables and Ozone in the Southeastern US – Part 2: Sensitivity to Dry-Deposition Parameterizations

Ozone (O3) dry deposition is a major O3 sink. As a follow-up study of Huang et al. (2021), we quantify the impact of satellite soil moisture (SM) on model representations of this process when different dry-deposition parameterizations are implemented, based on which the implications for interpreting O3 air pollution levels and assessing the O3 impacts on human and ecosystem health are provided. The SM data from NASA's Soil Moisture Active Passive mission are assimilated into the Noah-Multiparameterization (Noah-MP) land surface model within the NASA Land Information System framework, semicoupled with Weather Research and Forecasting model with online Chemistry (WRF-Chem) regional-scale simulations covering the southeastern US. Major changes in the modeling system used include enabling the dynamic vegetation option, adding the irrigation process, and updating the scheme for the surface exchange coefficient. Two dry-deposition schemes are implemented, i.e., the Wesely scheme and a “dynamic” scheme, in the latter of which dry-deposition parameterization is coupled with photosynthesis and vegetation dynamics. It is demonstrated that, when the dynamic scheme is applied, the simulated O3 dry-deposition velocities vd and their stomatal and cuticular portions, as well as the total O3 fluxes Ft, are larger overall; vd and Ft are 2–3 times more sensitive to the SM changes due to the data assimilation (DA). Further, through case studies at two forested sites with different soil types and hydrological regimes, we highlight that, applying the Community Land Model type of SM factor controlling stomatal resistance (i.e., β factor) scheme in replacement of the Noah-type β factor scheme reduced the v(d) sensitivity to SM changes by ∼75 % at one site, while it doubled this sensitivity at the other site. Referring to multiple evaluation datasets, which may be associated with variable extents of uncertainty, the model performance of vegetation, surface fluxes, weather, and surface O3 concentrations shows mixed responses to the DA, some of which display land cover dependency. Finally, using model-derived concentration- and flux-based policy-relevant O3 metrics as well as their matching exposure–response functions, the relative biomass/crop yield losses for several types of vegetation/crops are estimated to be within a wide range of 1 %–17 %. Their sensitivities to the model's dry-deposition scheme and the implementation of SM DA are discussed.

Min Huang↗

New instrumentation and tools to study dynamic (non-equilibrium) physico-chemical interfacial processes across large length, rates, and time scales (Final Report)

The hydrophobic interaction dictates the behavior of many systems relevant to biological and industrial processes. However, the origin of the hydrophobic interaction has been extensively debated. In this article, we demonstrate the use of nanostructured surfaces to control the range of the hydrophobic interaction. We measure the strength and range of the hydrophobic interaction using bubble probe force microscopy, and interpret our results using the Stokes–Reynolds–Young–Laplace model. Our findings give insights into the mechanisms of the hydrophobic interaction and will lead to better understanding and control of diverse colloidal and interfacial phenomena.

Moskovits, Martin↗

Practical Hamiltonian learning with unitary dynamics and Gibbs states

We study the problem of learning the parameters for the Hamiltonian of a quantum many-body system, given limited access to the system. In this work, we build upon recent approaches to Hamiltonian learning via derivative estimation. We propose a protocol that improves the scaling dependence of prior works, particularly with respect to parameters relating to the structure of the Hamiltonian (e.g., its locality k). Furthermore, by deriving exact bounds on the performance of our protocol, we are able to provide a precise numerical prescription for theoretically optimal settings of hyperparameters in our learning protocol, such as the maximum evolution time (when learning with unitary dynamics) or minimum temperature (when learning with Gibbs states). Thanks to these improvements, our protocol has practical scaling for large problems: we demonstrate this with a numerical simulation of our protocol on an 80-qubit system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantitative comparison between sub-millisecond time resolution single-molecule FRET measurements and 10-second molecular simulations of a biosensor protein

Molecular Dynamics (MD) simulations seek to provide atomic-level insights into conformationally dynamic biological systems at experimentally relevant time resolutions, such as those afforded by single-molecule fluorescence measurements. However, limitations in the time scales of MD simulations and the time resolution of single-molecule measurements have challenged efforts to obtain overlapping temporal regimes required for close quantitative comparisons. Achieving such overlap has the potential to provide novel theories, hypotheses, and interpretations that can inform idealized experimental designs that maximize the detection of the desired reaction coordinate. Here, we report MD simulations at time scales overlapping with in vitro single-molecule Förster (fluorescence) resonance energy transfer (smFRET) measurements of the amino acid binding protein LIV-BP SS at sub-millisecond resolution. Computationally efficient all-atom structure-based simulations, calibrated against explicit solvent simulations, were employed for sampling multiple cycles of LIV-BP SS clamshell-like conformational changes on the time scale of seconds, examining the relationship between these events and those observed by smFRET. The MD simulations agree with the smFRET measurements and provide valuable information on local dynamics of fluorophores at their sites of attachment on LIV-BP SS and the correlations between fluorophore motions and large-scale conformational changes between LIV-BP SS domains. We further utilize the MD simulations to inform the interpretation of smFRET data, including Förster radius (R 0 ) and fluorophore orientation factor (κ 2 ) determinations. The approach we describe can be readily extended to distinct biochemical systems, allowing for the interpretation of any FRET system conjugated to protein or ribonucleoprotein complexes, including those with more conformational processes, as well as those implementing multi-color smFRET.

59 BASIC BIOLOGICAL SCIENCES↗

Data driven learning of robust nonlocal models

Nonlocal models use integral operators that embed length-scales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given physical system, such as laboratory mechanical property tests. In contrast, molecular dynamics (MD) does not require these integrands, but it suffers from computational limitations in the length and time scales it can address. To combine the strengths of both methods and to obtain a coarse-grained, homogenized continuum model that efficiently and accurately captures materials' behavior, we propose a learning framework to extract, from MD data, an optimal nonlocal model as a surrogate for MD displacements. Our framework guarantees that the resulting model is mathematically well-posed, physically consistent, and that it generalizes well to settings that are different from the ones used during training. The efficacy of this approach is demonstrated with several numerical tests for single layer graphene both in the case of perfect crystal and in the presence of thermal noise.

36 MATERIALS SCIENCE↗

Optimality and Noise Resilience of Critical Quantum Sensing

We compare critical quantum sensing to passive quantum strategies to perform frequency estimation, in the case of single-mode quadratic Hamiltonians. We show that, while in the unitary case both strategies achieve precision scaling quadratic with the number of photons, in the presence of dissipation this is true only for critical strategies. We also establish that working at the exceptional point or beyond threshold provides suboptimal performance. This critical enhancement is due to the emergence of a transient regime in the open critical dynamics, and is invariant to temperature changes. When considering both time and system size as resources, for both strategies the precision scales linearly with the product of the total time and the number of photons, in accordance with fundamental bounds. However, we show that critical protocols outperform optimal passive strategies if preparation and measurement times are not negligible. Our results are applicable to a broad variety of critical sensors whose phenomenology can be reduced to that of a single-mode quadratic Hamiltonian, including systems described by finite-component and fully connected models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗