Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Robust stability”

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 541 records · Page 30

IPL-Annealed Mixed-Cation Perovskites with Robust Coating Window toward Scalable Manufacturing of Commercial Perovskite Solar Cells

Perovskite solar cells (PSCs) are a promising alternative solar technology, but the technical challenges of (1) stability/durability, (2) efficiency at scale, and (3) manufacturability must be overcome to achieve widespread PSC commercialization. The challenge of balancing solution ink formulation and scalable manufacturing is often overlooked in the literature, where focus is on adopting inks from processes that will not scale. In this study, we apply a classical roll-to-roll manufacturing perspective, utilizing both compositional engineering and intense pulsed light (IPL) annealing, to develop a mixed-cation perovskite ink with a robust coating window that simultaneously solve issues of stability and manufacturability for PSCs. Our method resulted in blade-coated, flexible, mixed-cation PSCs on ITO-PET substrates with a champion power conversion efficiency (PCE) efficiency of 16.7% using IPL annealing of the absorber layer and, to our knowledge, is one of the fastest processing methods for the perovskite layer. This overall reduction in processing time with a stable ink represents an advance toward the scaled production of perovskite solar cells on flexible substrates.

flexible perovskite solar cells↗

Digital twin framework for PIP-II linac: AI-driven multi-scale modeling from ion source to 800 MeV

The PIP-II superconducting linac at Fermilab is designed to deliver multi-megawatt proton beams for neutrino physics and other high-intensity applications. To expedite commissioning and enhance operational reliability, we have developed an EPICS-based data flow framework that seamlessly integrates digital twins (DT) with physical twins (PT). These digital twins comprise high-fidelity beam dynamics models or data-driven surrogate models connected to their physical counterparts through real-time diagnostics and advanced machine-learning algorithms.Central to this framework is Linac_Gen, an accelerated simulation tool that incorporates convolutional neural networks, random forests, and genetic algorithms to provide up to a tenfold speedup in optimizing the accelerator geometry model. An EPICS translator layer ensures interoperability by efficiently mapping lattice parameters across diverse simulation platforms.Our EPICS-based framework supports multiple operational modes—monitoring, passive learning, closed-loop control, and online learning—covering the entire machine lifecycle. By leveraging HPC resources and multi-objective optimization techniques, the digital twin enables adaptive trajectory correction, real-time fault detection, and predictive modeling of beam stability. This comprehensive approach paves the way for robust, high-intensity operation and data-driven accelerator R&D at Fermilab.

Pathak, Abhishek [Fermilab]↗

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advanced controls for stability assessment of solar dynamics space power generation

In support of the power requirements for the Space Station Alpha (SSA), a joint program by the U.S. and Russia for a permanently manned space station to be launched into orbit by 1998, a robust control scheme is needed to assure the stability of the rotating machines that will be integrated into the power subsystem. A framework design and systems studies for modeling and analysis is presented. It employs classical d-q axes machine model with voltage/frequency dependent loads. To guarantee that design requirements and necessary trade studies are done, a functional analysis tool CORE is used for the study. This provides us with different control options for stability assessment. Initial studies and recommendations using advanced simulation tools are also presented. The benefits of the stability/control scheme for evaluating future designs and power management are discussed.

Momoh, James A.↗

Decentralized control of large-scale systems: Fixed modes, sensitivity and parametric robustness

The problem of stabilization/pole placement under structural constraints of large scale linear systems is discussed. The existence of a solution to this problem is expressed in terms of fixed modes. The aim is to provide a bibliographic survey of the available results concerning the fixed modes (characterization, elimination, control structure selection to avoid them, control design in their absence) and to present the author's contribution to this problem which can be summarized by the use of the mode sensitivity concept to detect or to avoid them, the use of vibrational control to stabilize them, and the addition of parametric robustness considerations to design an optimal decentralized robust control.

Tarras, A.↗

Asymptotic Linearity of Optimal Control Modification Adaptive Law with Analytical Stability Margins

Optimal control modification has been developed to improve robustness to model-reference adaptive control. For systems with linear matched uncertainty, optimal control modification adaptive law can be shown by a singular perturbation argument to possess an outer solution that exhibits a linear asymptotic property. Analytical expressions of phase and time delay margins for the outer solution can be obtained. Using the gradient projection operator, a free design parameter of the adaptive law can be selected to satisfy stability margins.

Nguyen, Nhan T.↗

A generalized class of strongly stable and dimension-free T-RPMD integrators

Here, recent work shows that strong stability and dimensionality freedom are essential for robust numerical integration of thermostatted ringpolymer molecular dynamics (T-RPMD) and path-integral molecular dynamics, without which standard integrators exhibit non-ergodicity and other pathologies [R. Korol et al., J. Chem. Phys. 151, 124103 (2019) and R. Korol et al., J. Chem. Phys. 152, 104102 (2020)]. In particular, the BCOCB scheme, obtained via Cayley modification of the standard BAOAB scheme, features a simple reparametrization of the free ring-polymer sub-step that confers strong stability and dimensionality freedom and has been shown to yield excellent numerical accuracy in condensed-phase systems with large time steps. Here, we introduce a broader class of T-RPMD numerical integrators that exhibit strong stability and dimensionality freedom, irrespective of the Ornstein–Uhlenbeck friction schedule. In addition to considering equilibrium accuracy and time step stability as in previous work, we evaluate the integrators on the basis of their rates of convergence to equilibrium and their efficiency at evaluating equilibrium expectation values. Within the generalized class, we find BCOCB to be superior with respect to accuracy and efficiency for various configuration-dependent observables, although other integrators within the generalized class perform better for velocity-dependent quantities. Extensive numerical evidence indicates that the stated performance guarantees hold for the strongly anharmonic case of liquid water. Both analytical and numerical results indicate that BCOCB excels over other known integrators in terms of accuracy, efficiency, and stability with respect to time step for practical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

State-Constrained Grid-Forming Inverter Control for Robust Operation of AC Microgrids

The control of grid-forming inverters has become an important problem to achieve high efficiency, reliability, and resilience of the future power grid. The distributed control that utilizes a sparse communication network has been proposed for AC microgrids, in which each DG only has access to the information of a small number of its neighboring DGs, reducing the computational complexity and improving the resiliency to faults and unknown parameters. However, these approaches have not provided a clear analysis of the coupling between active/reactive power injections and the voltage magnitude/frequency or a distributed updating law with a solid theoretical foundation. Additionally, in the existing approaches, tracking the references from the secondary control may violate operational constraints.

Xu, Ying↗

Sliding Mode Control of the X-33 with an Engine Failure

Ascent flight control of the X-3 is performed using two XRS-2200 linear aerospike engines. in addition to aerosurfaces. The baseline control algorithms are PID with gain scheduling. Flight control using an innovative method. Sliding Mode Control. is presented for nominal and engine failed modes of flight. An easy to implement, robust controller. requiring no reconfiguration or gain scheduling is demonstrated through high fidelity flight simulations. The proposed sliding mode controller utilizes a two-loop structure and provides robust. de-coupled tracking of both orientation angle command profiles and angular rate command profiles in the presence of engine failure, bounded external disturbances (wind gusts) and uncertain matrix of inertia. Sliding mode control causes the angular rate and orientation angle tracking error dynamics to be constrained to linear, de-coupled, homogeneous, and vector valued differential equations with desired eigenvalues. Conditions that restrict engine failures to robustness domain of the sliding mode controller are derived. Overall stability of a two-loop flight control system is assessed. Simulation results show that the designed controller provides robust, accurate, de-coupled tracking of the orientation angle command profiles in the presence of external disturbances and vehicle inertia uncertainties, as well as the single engine failed case. The designed robust controller will significantly reduce the time and cost associated with flying new trajectory profiles or orbits, with new payloads, and with modified vehicles

Shtessel, Yuri B.↗

Data efficiency and extrapolation trends in neural network interatomic potentials

Abstract Recently, key architectural advances have been proposed for neural network interatomic potentials (NNIPs), such as incorporating message-passing networks, equivariance, or many-body expansion terms. Although modern NNIP models exhibit small differences in test accuracy, this metric is still considered the main target when developing new NNIP architectures. In this work, we show how architectural and optimization choices influence the generalization of NNIPs, revealing trends in molecular dynamics (MD) stability, data efficiency, and loss landscapes. Using the 3BPA dataset, we uncover trends in NNIP errors and robustness to noise, showing these metrics are insufficient to predict MD stability in the high-accuracy regime. With a large-scale study on NequIP, MACE, and their optimizers, we show that our metric of loss entropy predicts out-of-distribution error and data efficiency despite being computed only on the training set. This work provides a deep learning justification for probing extrapolation and can inform the development of next-generation NNIPs.

36 MATERIALS SCIENCE↗

High Energy Aqueous Lithium-Ion Batteries (Final Report)

Lithium-ion batteries (LIBs) are a widely used energy storage device. For decades, researchers have been studying aqueous electrolytes due to their non-flammability and environmental friendless. However, the hydrogen evolution reaction (HER) and oxygen evolution reaction (OER) of water greatly limit the voltage window of the aqueous electrolyte. The recent concept of ‘solvent-in-salt’ electrolytes, namely ‘water-in-salt’ electrolytes (WISE) for aqueous batteries and ‘super-concentrated’ electrolytes for non-aqueous batteries, greatly expands the electrochemical window of these battery systems and consequently their energy density. This expansion occurs because the large number of salt aggregates in ‘solvent-in-salt’ electrolytes enables preferential salt decomposition over solvent decomposition to form robust interphase layers on electrodes, thus enhancing their stability in a wide range of operating voltages. However, even when the salt concentration reached 63m, regardless of cost, the cathodic potential of the WISE expanded only to 1.75V with an overall electrochemical stability window of <3.25V. To design the next generation of high energy aqueous LIBs, electrolytes that possess a wider electrochemical stability window and lower cost but still maintain the merit of non-flammability are much needed. In this project, we designed the aqueous LIBs with widened voltage window. We reported a cheap and non-flammable ternary eutectic electrolyte that can form a robust solid electrolyte interface (SEI)so that the cathodic limiting potential can be decreased to <1.5 V for Li 4 Ti 5 O 12 anodes, while maintaining both a high ionic conductivity of >10 -3 S cm –1 and a low viscosity. We demonstrated the super-stable2.5V LiMn 2 O 4 ||Li 4 Ti 5 O 12 pouch cell with mass loading of 2.5 mAh/cm. Further, LiCl-LiBr-Graphite/Li full cells were designed with high energy density. Finally, we use a bi-layer design to enlarge the voltage window of aqueous electrolyte, which can endow the stable cycling of aqueous Lithium metal batteries.

25 ENERGY STORAGE↗

Upper Limit of Weights in TAI Computation

The international reference time scale International Atomic Time (TAI) computed by the Bureau International des Poids et Mesures (BIPM) relies on a weighted average of data from a large number of atomic clocks. In it, the weight attributed to a given clock depends on its long-term stability. In this paper the TAI algorithm is used as the basis for a discussion of how to implement an upper limit of weight for clocks contributing to the ensemble time. This problem is approached through the comparison of two different techniques. In one case, a maximum relative weight is fixed: no individual clock can contribute more than a given fraction to the resulting time scale. The weight of each clock is then adjusted according to the qualities of the whole set of contributing elements. In the other case, a parameter characteristic of frequency stability is chosen: no individual clock can appear more stable than the stated limit. This is equivalent to choosing an absolute limit of weight and attributing this to to the most stable clocks independently of the other elements of the ensemble. The first technique is more robust than the second and automatically optimizes the stability of the resulting time scale, but leads to a more complicated computatio. The second technique has been used in the TAI algorithm since the very beginning. Careful analysis of tests on real clock data shows that improvement of the stability of the time scale requires revision from time to time of the fixed value chosen for the upper limit of absolute weight. In particular, we present results which confirm the decision of the CCDS Working Group on TAI to increase the absolute upper limit by a factor of 2.5. We also show that the use of an upper relative contribution further helps to improve the stability and may be a useful step towards better use of the massive ensemble of HP 507IA clocks now contributing to TAI.

Thomas, Claudine↗

Stabilization of Photocatalytically Active Uranyl Species in a Uranyl–Organic Framework for Heterogeneous Alkane Fluorination Driven by Visible Light

When photoactivated, the uranyl ion is a robust oxidant capable of abstracting hydrogen atoms from nonactivated C–H bonds. Yet, the highly reactive singly reduced [U V O 2 ] + intermediate is unstable with respect to disproportionation to the uranyl dication and insoluble tetravalent uranium phases, which limits the usage of uranyl ions as robust photocatalysts. Herein, we demonstrate that photoactivated uranyl ions can be stabilized by immobilizing and separating them spatially in a uranyl–organic framework heterogeneous catalyst, NU-1301. The visible-light-photoactivated uranyl ions in NU-1301 demonstrated longer-lived U(V) and radicals than those in homogeneous counterparts, as evidenced by X-ray photoelectron spectroscopy and time-dependent electron paramagnetic resonance, leading to higher turnovers and enhanced stability for the fluorination of nonactivated alkanes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗

Virtual Deformation Control of the X-56A Model with Simulated Fiber Optic Sensors

A robust control law design methodology is presented to stabilize the X-56A model and command its wing shape. The X-56A was purposely designed to experience flutter modes in its flight envelope. The methodology introduces three phases: the controller design phase, the modal filter design phase, and the reference signal design phase. A mu-optimal controller is designed and made robust to speed and parameter variations. A conversion technique is presented for generating sensor strain modes from sensor deformation mode shapes. The sensor modes are utilized for modal filtering and simulating fiber optic sensors for feedback to the controller. To generate appropriate virtual deformation reference signals, rigid-body corrections are introduced to the deformation mode shapes. After successful completion of the phases, virtual deformation control is demonstrated. The wing is deformed and it is shown that angle-ofattack changes occur which could potentially be used to an advantage. The X-56A program must demonstrate active flutter suppression. It is shown that the virtual deformation controller can achieve active flutter suppression on the X-56A simulation model.

Suh, Peter M.↗

Virtual Deformation Control of the X-56A Model with Simulated Fiber Optic Sensors

A robust control law design methodology is presented to stabilize the X-56A model and command its wing shape. The X-56A was purposely designed to experience flutter modes in its flight envelope. The methodology introduces three phases: the controller design phase, the modal filter design phase, and the reference signal design phase. A mu-optimal controller is designed and made robust to speed and parameter variations. A conversion technique is presented for generating sensor strain modes from sensor deformation mode shapes. The sensor modes are utilized for modal filtering and simulating fiber optic sensors for feedback to the controller. To generate appropriate virtual deformation reference signals, rigid-body corrections are introduced to the deformation mode shapes. After successful completion of the phases, virtual deformation control is demonstrated. The wing is deformed and it is shown that angle-of-attack changes occur which could potentially be used to an advantage. The X-56A program must demonstrate active flutter suppression. It is shown that the virtual deformation controller can achieve active flutter suppression on the X-56A simulation model.

Modal Filter↗

Experimental and analytical study of the hydrodynamic and single and two-phase convective heat transfer performance of flexible PDMS microchannels with micropillar arrays

Various copper and silicon based thermal management systems are used in the cooling of electronics. However, the rigid nature of these materials along with their high thermal and electrical conductivity pose a difficulty in developing direct contact embedded flexible cooling systems that can offer robust cooling performance. The low density, thermal stability, chemical inertness, and electrical insulation of Polydimethylsiloxane (PDMS) make it an ideal material to develop lightweight direct contact thermal management systems for electronics. Its ease of fabrication with tunable flexibility provides the opportunity to go beyond traditional electronics and develop advanced active and passive thermal management systems for a wide–range of applications in foldable and wearable electronics, liquid cooling garments, microgravity, and electric motors. In this study, a flexible PDMS based microchannel with micropillar arrays, which enhance the thermal performance of the device through capillary-assisted flow, has been developed. The hydrodynamic and convection heat transfer performance of three PDMS wick pillar geometries, ranging from a porosity of 0.8–0.91, are investigated and compared under single-phase and two-phase conditions. Dielectric coolant FC-3283 is employed and permeability measurements are made for mass fluxes ranging from 53 kg/m 2 s to 369 kg/m 2 s. Given its conformability, the device demonstrates a deviation from Darcy’s Law, within the laminar regime, with an increasing permeability with mass flux at the rate of ~0.5–0.8 Darcy/(kg/m 2 s). A semi-analytical model has been developed and reported to quantify the conformability of the device. The heat transfer performance is experimentally evaluated using the same dielectric fluid for mass flux ranging from 105 kg/m 2 s to 420 kg/m 2 s with heat fluxes ranging from 1.5 W/cm 2 to 16 W/cm 2 . Heat transfer coefficients of up to 7000 W/m 2 K are observed, which are comparable to copper and silicon microchannels. The effect of porosity on the single phase thermal performance has been evaluated against the pumping power to provide a basis for thermal management system design. Finally, high-speed imaging is performed to study the two-phase flow characteristics to provide insight into the vapor formation and removal.

42 ENGINEERING↗

Novel PtNi single-atom–nanocluster (SA–NC) ensembles promote Tafel kinetics and ampere-class AEM hydrogen evolution

The development of efficient, durable, and low-PGM electrocatalysts for the hydrogen evolution reaction (HER) in alkaline media is critical for next-generation electrolysis technologies. We report a facile two-step synthesis of highly dispersed PtNi and PtNi-nitride nanoclusters (NCs) (~2.3 nm) with ultralow Pt content (0.5 at.%) anchored on N-doped Vulcan carbon. Structural and compositional characterization via XAS, XPS, HAADF-STEM, HRTEM, and EDS mapping established key structure–activity relationships across varying Pt/Ni ratios and pyrolysis temperatures. The Pt0.5Ni0.5/C-750 catalyst, an ensemble of PtNi M-N-C type single-atom (SA) moieties with neighboring PtNi nanoclusters (NC), exhibited superior HER performance in alkaline media, achieving overpotentials of 30, 115, and 210 mV at 10, 100, and 500 mA cm−2, respectively. Despite at a lower Pt content, this novel SA–NC ensemble outperformed commercial Pt/C by ~36%. A standardized literature comparison with contemporary Pt- and Ru-doped analogues reveals the as-prepared Pt0.5Ni0.5/C-750 to sit at the apex of Tafel-limited kinetics and low overpotential at 100 mA cm−2. Tafel-limited Tafel slopes in both alkaline and acidic regimes confirm favorable proton recombination kinetics. Mass activities at 200 mV reached 13.8 and 18.84 A mgPt−1 in alkaline and acidic media, respectively. However, excessive nitridation (e.g., at 650 °C) adversely altered Pt electronic structure and HER kinetics. While Ni enhanced alkaline HER, acidic HER favored Ni-free analogues. Pt0.5Ni0.5/C-750 also demonstrated robust temperature responsiveness and 300-h operational stability at high current densities (0.5–1.0 A cm−2) in MEA tests. This work presents a scalable strategy for designing thermally responsive, durable, and compositionally tunable NC catalysts with neighboring SA moieties for alkaline electrolysis.

AEMWE↗