Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear dynamical systems”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Systematic Characterization on Optical and Phonon Properties of Cs[Pb x Sn 1– x ]I 3 Alloyed Perovskites via First-Principles Modeling

Here, understanding the optical properties and phonon thermodynamic features of materials enables the determination of their viability for optoelectronic applications including optical-based gas sensors. Metal halide perovskites ABX 3 as well as their mixed alloys can be attractive in applications of this kind. In particular, the formation of alloyed perovskites A[B x B’ 1–x ]X 3 modulates the functional characteristics that can overtake the single-halide ABX 3 . To this end, we use the first-principles simulations to model the effect of B-site alloying on the optical and phonon properties of cubic Cs[Pb x Sn 1–x ]I 3 (x = 0.75, 0.50, and 0.25). We identified the systematic shifts of optical properties as a function of alloying concentration, such as dielectric constants, refraction indices, and absorption coefficients. It was also found that a linear correlation exists between the optical constants and the band gaps of these materials, which satisfies the Penn model. Modeling of phonons in Cs[Pb x Sn 1–x ]I 3 showed that dissipating the dynamically unstable phonons can be feasible at a carefully chosen alloying concentration, which are originally present in their single-halide Pb- and Sn-ends. Our findings pave a way to optimize the alloying conditions of the alloyed perovskite systems of interest for tuning their optical and phonon stability characteristics.

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

Additively Manufactured Compliant Hybrid Gas Thrust Bearing for Supercritical Carbon Dioxide Turbomachinery: Design and Proof of Concept Testing

The following paper presents a new type of gas lubricated thrust bearing that utilizes additive manufacturing or also known as direct metal laser melting (DMLM) to fabricate the bearing. The motivation for the new bearing concept is derived from the need for highly efficient supercritical carbon dioxide (sCO2) turbomachinery in the mega-watt power range. The paper provides a review of existing gas thrust-bearing technology, outlines the need for the new DMLM concept, and discusses proof-of-concept testing results. The new concept combines hydrostatic pressurization with individual tilting pads that are flexibly mounted using hermetic squeeze film dampers (HSFD) in the bearing-pad support. This paper describes the thrust-bearing concept and discusses the final design approach. Proof-of-concept testing in air for a 6.8 in. (173 mm) outer diameter thrust gas bearing was performed; with thrust loading, up to 1500 lbs (6.67 kN) and a thrust runner speed of 10krpm (91 m/s tip speed). The experiments were performed with a bent shaft resulting in thrust runner axial vibration magnitudes of 2.9 mils (74 μm) p-p and dynamic thrust loads of 270 lbs (1.2 kN) p-p. In addition, force deflection characteristics and stiffness coefficients of the bearing system are presented for an inlet hydrostatic pressure of 380 psi (2.62 MPa). Results at 10 krpm show that the pad support architecture was able to sustain high levels of dynamic misalignment equaling 6 times the nominal film clearance while demonstrating a unit load-carrying capacity of 55 psi (0.34 Mpa). Gas-film force deflection tests portrayed nonlinear behavior like a hardening spring, while the bearing pad support stiffness was measured to be linear and independent of gas film thickness.

Engineering↗
DOE OSTI2023

Elastic Flow Modeling for Hydropower Digital Twins

This report details the elastic, unsteady, one-dimensional flow equations that are used to model flow through a penstock in a hydropower facility. The elastic flow model is accurate even in cases of fast transients, long penstock length, and high gravitation head. The compressibility of the water and elasticity of the pipe walls are explicitly accounted for in the developed models so that the water hammer phenomena can be accurately captured. The elastic flow model is coupled to a mechanistic turbine model to resolve the dynamic feedbacks between the elastic water column and the turbine rotation rate. A finite volume method is employed to solve the governing equations, and the method is shown to have low numerical diffusion in sharp gradient phenomena, such as those encountered when simulating water hammer. The method is applied to single dimensional linear advection benchmark problem and numerical results are compared with analytic results. This is followed on by an application to a full hydropower system with a coupled turbine. The elastic flow model results are compared against an inelastic model.

13 HYDRO ENERGY↗
DOE OSTI2022

Time Evolution in Quantum Cosmology

The quantum description of time evolution in non-linear gravitational systems such as cosmological space-times is not well understood. We show, in the simplified setting of mini-superspace, that time evolution of this system can be obtained using a gauge fixed path integral, as long as one does not integrate over proper time. Using this gauge fixed action we can construct a Hamiltonian in the coherent - or classical - state basis. We show that by construction the coherent states satisfy the classical dynamical equations of General Relativity. They do not satisfy the Hamiltonian constraint. A consequence of this is that the Wheeler-DeWitt equation should not be satisfied in quantum gravity. Classical states have a natural non-trivial time evolution since they are not eigenstates of the Hamiltonian. A general feature of the unconstrained quantum theory of gravity is the prediction of a pressureless dark matter component of either sign energy density in the classical universe which may lead to novel phenomenology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Non-autoregressive time-series methods for stable parametric reduced-order models

Advection-dominated dynamical systems, characterized by partial differential equations, are found in applications ranging from weather forecasting to engineering design where accuracy and robustness are crucial. There has been significant interest in the use of techniques borrowed from machine learning to reduce the computational expense and/or improve the accuracy of predictions for these systems. These rely on the identification of a basis that reduces the dimensionality of the problem and the subsequent use of time series and sequential learning methods to forecast the evolution of the reduced state. Often, however, machine-learned predictions after reduced-basis projection are plagued by issues of stability stemming from incomplete capture of multiscale processes as well as due to error growth for long forecast durations. To address these issues, we have developed a non-autoregressive time series approach for predicting linear reduced-basis time histories of forward models. In particular, we demonstrate that non-autoregressive counterparts of sequential learning methods such as long short-term memory (LSTM) considerably improve the stability of machine-learned reduced-order models. Further, we evaluate our approach on the inviscid shallow water equations and show that a non-autoregressive variant of the standard LSTM approach that is bidirectional in the principal component directions obtains the best accuracy for recreating the nonlinear dynamics of partial observations. Moreover-and critical for many applications of these surrogates-inference times are reduced by three orders of magnitude using our approach, compared with both the equation-based Galerkin projection method and the standard LSTM approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Extended nonergodic regime and spin subdiffusion in disordered SU(2)-symmetric Floquet systems

In this work, we explore thermalization and quantum dynamics in a one-dimensional disordered SU(2)-symmetric Floquet model, where a many-body localized phase is prohibited by the non-Abelian symmetry. Despite the absence of localization, we find an extended nonergodic regime at strong disorder where the system exhibits nonthermal behaviors. In the strong disorder regime, the level spacing statistics exhibit neither a Wigner-Dyson nor a Poisson distribution, and the spectral form factor does not show a linear-in-time growth at early times characteristic of random matrix theory. The average entanglement entropy of the Floquet eigenstates is subthermal, although violating an area-law scaling with system sizes. We further compute the expectation value of local observables and find strong deviations from the eigenstate thermalization hypothesis. The infinite-temperature spin autocorrelation function decays at long times as $t^{–β}$ with $\textit{β}$ < 0.5, indicating subdiffusive transport at strong disorders.

36 MATERIALS SCIENCE↗
DOE OSTI2020

As a Matter of State: The Role of Thermodynamics in Magnetohydrodynamic Turbulence

Turbulence simulations play a key role in advancing the general understanding of the physical properties of turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence simulations are often conducted in the isothermal limit. Given that the majority of astrophysical systems are not governed by isothermal dynamics, we aim to quantify the impact of thermodynamics on the physics of turbulence, through varying adiabatic index, γ, combined with a range of optically thin cooling functions. Here, we present a suite of ideal magnetohydrodynamics simulations of thermally balanced stationary turbulence in the subsonic, super-Alfvénic, high ${\beta }_{{\rm{p}}}$ (ratio of thermal to magnetic pressure) regime, where turbulent dissipation is balanced by two idealized cooling functions (approximating linear cooling and free–free emission) and examine the impact of the equation of state by considering cases that correspond to isothermal, monatomic, and diatomic gases. We find a strong anticorrelation between thermal and magnetic pressure independent of thermodynamics, whereas the strong anticorrelation between density and magnetic field found in the isothermal case weakens with increasing γ. Similarly, the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing γ. This suggests that there exists a degeneracy in these relations with respect to thermodynamics and Mach number in this regime, which is dominated by slow magnetosonic modes. These results have implications for attempts to infer (e.g.,) Mach numbers from (e.g.,) Faraday rotation measurements, without additional information regarding the thermodynamics of the plasma. However, our results suggest that this degeneracy can be broken by utilizing higher-order moments of observable distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Site demonstration and performance evaluation of MPC for a large chiller plant with TES for renewable energy integration and grid decarbonization

Thermal energy storage (TES) for a cooling plant is a crucial resource for load flexibility. Traditionally, simple, heuristic control approaches, such as the storage priority control which charges TES during the nighttime and discharges during the daytime, have been widely used in practice, and shown reasonable performance in the past benefiting both the grid and the end-users such as buildings and district energy systems. However, the increasing penetration of renewables changes the situation, exposing the grid to a growing duck curve, which encourages the consumption of more energy in the daytime, and volatile renewable generation which requires dynamic planning. The growing pressure of diminishing greenhouse gas emissions also increases the complexity of cooling TES plant operations as different control strategies may apply to optimize operations for energy cost or carbon emissions. This paper presents a model predictive control (MPC), site demonstration and evaluation results of optimal operation of a chiller plant, TES and behind-meter photovoltaics for a campus-level district cooling system. The MPC was formulated as a mixed-integer linear program for better numerical and control properties. Compared with baseline rule-based controls, the MPC results show reductions of the excess PV power by around 25%, of the greenhouse gas emission by 10%, and of peak electricity demand by 10%.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Destabilization of binary mixing layer in supercritical conditions

Compressible mixing layer instabilities are of importance to a wide range of environmental and industrial applications. Past studies have focused on either ideal-gas or real-fluid thermodynamic regimes of single-species mixing layers. However, mixing layers of binary mixtures at supercritical conditions, commonly encountered in fuel injection systems, introduce additional complexities due to the added compositional degree-of-freedom. Moreover, the effect of strong variations in thermodynamic response functions across the Widom line on the binary mixing layer stability remains poorly understood. Thus, the objective of this study is to examine the coupling between the hydrodynamic instability and the real-fluid thermodynamics across the Widom line and its effects on the overall binary mixing layer dynamics. To this end, we develop a linear stability analysis of the full binary-species compressible transport equations coupled with the PC-SAFT equation of state. Analysis shows the existence of a novel instability mechanism that arises from juxtapositioning of the Widom-line transition and the hydrodynamic inflection point. This novel thermodynamically induced instability mechanism has the net effect of destabilizing the binary mixing layer at lowering supercritical conditions towards the critical pressure point. This is in contrast to previous stability analyses of supercritical single-species mixing layers, where increasing pressure destabilizes the flow due to its effect on reducing the density stratification. The discovered thermodynamically induced instability mechanism of binary mixing flows highlights the need for an extension of classical instability criteria to incorporate the effect of strong variations in the thermodynamic response functions across the Widom line on mixing layer instability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

When do molecular polaritons behave like optical filters?

This review outlines several linear optical effects featured by molecular polaritons arising in the collective strong light–matter coupling regime. Under weak laser irradiation and when the single-molecule light–matter coupling can be neglected (often in the limit when the number of molecules per photon mode is large), we show that the excited-state molecular dynamics under collective strong coupling can be exactly replicated without the cavity using a shaped (or “filtered”) laser, whose field amplitude is enhanced by the cavity quality factor, shining on the bare molecules. As a consequence, the absorption within a cavity can be understood as the overlap between the polariton transmission and the bare molecular absorption, suggesting that polaritons act in part as optical filters. This framework demystifies and provides a straightforward explanation for a large class of experiments and theoretical models in molecular polaritonics, highlighting that the same effects can be achieved without the cavity with shaped laser pulses. With a few modifications, this simple conceptual picture can also be adapted to understand the incoherent nonlinear response of polaritonic systems. This review establishes a clear distinction between polaritonic phenomena that can be fully explained through classical linear optics and those that require a quantum electrodynamics approach. It also highlights the need to differentiate between effects that necessitate polaritons (i.e., hybrid light–matter states) and those that can occur in the weak coupling regime. Here, we further discuss that certain quantum optical effects like fluorescence can be partially described as optical filtering, whereas some others like cavity-induced Raman scattering go beyond this. Further exploration in these areas is needed to uncover novel polaritonic phenomena beyond optical filtering.

Schwennicke, Kai [University of California, San Di↗
DOE OSTI2025

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

On applications of quantum computing to plasma simulations

Quantum computing is gaining increased attention as a potential way to speed up simulations of physical systems, and it is also of interest to apply it to simulations of classical plasmas. However, quantum information science is traditionally aimed at modeling linear Hamiltonian systems of a particular form that is found in quantum mechanics, so extending the existing results to plasma applications remains a challenge. Here, we report a preliminary exploration of the long-term opportunities and likely obstacles in this area. First, we show that many plasma-wave problems are naturally representable in a quantumlike form and thus are naturally fit for quantum computers. Second, we consider more general plasma problems that include non-Hermitian dynamics (instabilities, irreversible dissipation) and nonlinearities. We show that by extending the configuration space, such systems can also be represented in a quantumlike form and thus can be simulated with quantum computers too, albeit that requires more computational resources compared to the first case. Third, we outline potential applications of hybrid quantum–classical computers, which include analysis of global eigenmodes and also an alternative approach to nonlinear simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Electron-irradiation induced creep in amorphous alloys

Electron-irradiation induced creep rates in amorphous alloys, a-SiO2, Fe79B16Si5, Cu60Ta40, and Cu50Ti50, were measured at room temperature using a miniaturized beam-bending apparatus within a transmission electron microscope operated at 200 keV. The creep rates of these amorphous samples increased nearly linearly with both e-beam current density and applied stress, while a reference crystalline (c-)SiO2 sample failed to creep under the same conditions. The irradiation induced creep compliance of a-SiO2 was ~ 15 times larger than that of Fe79B16Si5 and over 1,000 times larger than that of the two Cu alloys. Molecular dynamics computer simulations were employed to simulate electron irradiation induced creep using interatomic potentials representing amorphous Cu75Zr25, Ni80P20, and SiO2 as model systems. The irradiation induced creep compliances calculated for Cu75Zr25 during 200 keV electron irradiation provided good quantitative agreement with the two Cu-based alloys, but that for a-SiO2 was ~ 180 times too small. These results indicate that unlike neutron or ion-beam induced creep in a-SiO2, creep under electron irradiation is dominated by the effects of ionization owing largely to the far higher ratio of electronic stopping to nuclear stopping for electrons than for ions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A molecular conveyor belt-associated protein controls the rotational direction of the bacterial type 9 secretion system

ABSTRACT Many bacteria utilize the type 9 secretion system (T9SS) for gliding motility, surface colonization, and pathogenesis. This dual-function motor supports both gliding motility and protein secretion, where rotation of the T9SS plays a central role. Fueled by the energy of the stored proton motive force and transmitted through the torque of membrane-anchored stator units, the rotary T9SS propels an adhesin-coated conveyor belt along the bacterial outer membrane like a molecular snowmobile, thereby enabling gliding motion. However, the mechanisms controlling the rotational direction and gliding motility of T9SS remain elusive. Shedding light on this mechanism, we find that in the gliding bacteriumFlavobacterium johnsoniae, deletion of the C-terminus of the conveyor belt-associated protein GldJ controls and, in fact, reverses the rotational direction of T9SS from counterclockwise (CCW) to clockwise (CW). This suggests that the interface between the conveyor belt-associated protein GldJ and the T9SS ring protein GldK plays an important role in controlling the directionality of T9SS, potentially by modulating its interaction with the stator complex GldLM, which drives motor rotation. Combined with MD simulation of the T9SS stator units GldLM, we suggest a “tri-component gearset” model where GldJ controls the rotational direction of its driver, the T9SS, thus providing adaptive sensory feedback to influence the motility of the gliding bacterium. IMPORTANCE The type 9 secretion system (T9SS) is fundamental to bacterial gliding motility, pathogenesis, and surface colonization. Our findings reveal that the C-terminal region of the conveyor belt-associated protein GldJ functions as a molecular switch which is capable of reversing the rotational direction of T9SS. Through the coordinated actions of the T9SS stator units (akin to a driving motor), the GldK ring (the gear that converts rotational energy into linear movement), and GldJ, this machinery forms a smart conveyor belt system reminiscent of flexible or cognitive mechanical conveyors. Such advanced conveyors can alter their direction to adapt to shifting demands. Here, we show that the bacterial T9SS similarly adjusts its rotational bias based on feedback from the conveyor belt-associated protein GldJ. This dual-role feedback mechanism underscores an evolved, controllable biological snowmobile, offering new avenues for studying how bacteria fine-tune motility in dynamic environments.

Microbiology↗
DOE OSTI2024

Radiation induced non-linear oscillations in ITER baseline scenario plasmas in DIII-D

Abstract This work shows how the radiation brought about by metals or metal-equivalent radiators such as Kr and Xe produces non-linear dynamics on otherwise stationary β N flattops of DIII-D ITER Baseline Scenario demonstration discharges. The Kr and Xe gases are used to reproduce the radiative loss rates of W in present machines that operate at core temperatures much lower than the expected ITER temperature. Experiments on DIII-D with injection of Kr and Xe, as well as with sources of intrinsic metals reach the range of radiated fraction values expected in the ITER core and experience slow oscillations in temperature and radiated power. In many cases of high radiated fraction, the core temperature decreases enough for the safety factor profile to rise above the 1/1 rational surface, naturally eliminating sawteeth and occasionally producing a persistent helical core. The oscillations can be reproduced by a modified Lotka–Volterra system for temperature and radiated fraction if diffusion and noise are included, which indicates that the interplay between temperature and radiation can be the main cause of the cyclic nature of the system. A new physics based model which includes equations for temperature, density and input power can also reproduce the oscillations observed in the experiments. The present results suggest that the non-linearity of the system can be increased by the inclusion of the inherently non-linear alpha heating term, which is proportional to ∼ n e 2 T i 2 , and obtains oscillations in the model when added to an otherwise more stationary system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Non-intrusive data-driven model reduction for differential–algebraic equations derived from lifting transformations

In this paper we present a non-intrusive data-driven approach for model reduction of nonlinear systems. The approach considers the particular case of nonlinear partial differential equations (PDEs) that form systems of partial differential–algebraic equations (PDAEs) when lifted to polynomial form. Such systems arise, for example, when the governing equations include Arrhenius reaction terms (e.g., in reacting flow models) and thermodynamic terms (e.g., the Helmholtz free energy terms in a phase-field solidification model). Using the known structured form of the lifted algebraic equations, the approach computes the reduced operators for the algebraic equations explicitly, using straightforward linear algebra operations on the basis matrices. The reduced operators for the differential equations are inferred from lifted snapshot data using operator inference, which solves a linear least squares regression problem. The approach is illustrated for the nonlinear model of solidification of a pure material. The lifting transformations reformulate the solidification PDEs as a system of PDAEs that have cubic structure. The operators of the lifted system for this solidification example have affine dependence on key process parameters, permitting us to learn a parametric reduced model with operator inference. Numerical experiments show the effectiveness of the resulting reduced models in capturing key aspects of the solidification dynamics.

42 ENGINEERING↗
DOE OSTI2024

FARM User Guidance and Instructions

The FARM (Feasible Actuator Range Modifier) module is a component of the RAVEN-based FORCE framework for the analysis of Integrated Energy Systems (IES). FARM aids HERON in the solution of the power dispatch problem by evaluating feasible set-point signals to be issued to the control systems of the different IES unit components. Set-points need to satisfy limits on both production variables (i.e., the variables to be optimized such as the electrical power, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). To enforce all these limits, a two-stage approach is adopted. First, the power dispatcher algorithm in HERON module estimates set-points that meet the constraints on the production variables, e.g., power levels and power ramp rate limits. These constraints are called explicit constraints. Then, if necessary, FARM adjusts these set-points to ensure the respect of the limits on the process variables of interest, given the knowledge of the system dynamics acquired through machine learning algorithms. These constraints are called implicit constraints. From this standpoint, FARM constitutes a bridge between the HERON power dispatcher that adopts a simplified description of the IES unit (low-resolution physics) and the HYBRID high-fidelity models (high-resolution physics). In this report, the general workflow and the software requirements of FARM module are summarized, and the detailed instructions for installing FARM software, running built-in example cases, deriving Linear Parameter-Varying (LPV) state-space models, and using FARM for user-defined power dispatch problems are provided.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗