Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stabilization method”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Resolving phase stability in the Ti-O binary with first-principles statistical mechanics methods

The Ti-O system consists of a multitude of stable and metastable oxides that are used in wide ranging applications. In this work we investigate phase stability in the Ti-O binary from first principles. We perform a systematic search for ground state structures as a function of oxygen concentration by considering oxygen-vacancy and/or titanium-vacancy orderings over four parent crystal structures: (i) hcp Ti, (ii) ω-Ti, (iii) rocksalt, and (iv) hcp oxygen containing interstitial titanium. We explore phase stability at finite temperature using cluster expansion Hamiltonians and Monte Carlo simulations. The calculations predict a high oxygen solubility in hcp Ti and the stability of suboxide phases that undergo order-disorder transitions upon heating. Vacancy ordered rocksalt phases are also predicted at low temperature that disorder to form an extended solid solution at high temperatures. Predicted stable and metastable phase diagrams are qualitatively consistent with experimental observations, however, important discrepancies are revealed between first-principles density functional theory predictions of phase stability and the current understanding of phase stability in this system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Low-iridium stabilized ruthenium oxide anode catalyst for durable proton-exchange membrane water electrolysis

While mixing iridium (Ir) with ruthenium oxide (RuO 2 ) has proven to be an effective strategy for reducing Ir loading in anode catalysts for proton-exchange membrane (PEM) water electrolysers, achieving industrially relevant long-term stability typically requires an Ir-rich, Ru-lean combination. Here, in this study, by combining density functional theory with Metropolis Monte Carlo methods, we discovered that sufficient stabilization in the RuO 2 lattice could be achieved with less than 50 at.% of Ir, and that Ir in the first subsurface layer plays a critical role. By effectively dispersing Ir dopants within the RuO 2 lattice, we demonstrated an Ir:Ru atomic ratio of only 1:6 that exhibited exceptional stability for over 1,500 h of continuous water electrolysis at 2 A cm −2 . Our Ru 6 IrO x catalyst has the potential to reduce Ir loading by 80% compared with current commercial PEM water electrolysers, and its stability was further validated under industrial testing conditions in a 25-cm 2 PEM electrolyser.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Semi-Implicit Computation of Fast Modes in a Scheme Integrating Slow Modes by a Leapfrog Method Based on a Selective Implicit Time Filter

Abstract A scheme for integration of atmospheric equations containing terms with differing time scales is developed. The method employs a filtered leapfrog scheme utilizing a fourth-order implicit time filter with one function evaluation per time step to compute slow-propagating phenomena such as advection and rotation. The terms involving fast-propagating modes are handled implicitly with an unconditionally stable method that permits application of larger time steps and faster computations compared to fully explicit treatment. Implementation using explicit and recurrent formulation is provided. Stability analysis demonstrates that the method is conditionally stable for any combination of frequencies involved in the slow and fast terms as they approach the origin. The implicit filter used in the method damps the computational modes without noticeably sacrificing the accuracy of the physical mode. TheO[(Δt 4 )] accuracy for amplitude errors achieved by the implicitly filtered leapfrog is preserved in applications where terms responsible for fast propagation are integrated with a semi-implicit method. Detailed formulation of the method for soundproof nonhydrostatic anelastic equations is provided. Procedures for implementation in global spectral shallow-water models are also given. Examples comparing numerical and analytical solutions for linear gravity waves demonstrate the accuracy of the scheme. The performance is also shown in more practical nonlinear applications, where numerical solutions accomplished by the method are evaluated against those computed from a scheme where the slow terms are handled by the third-order Runge–Kutta scheme. It demonstrates that the method is able to accurately resolve fine-scale dynamics of Kelvin–Helmholtz shear instabilities, the evolution of density current, and nonlinear drifts of twin tropical cyclones.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part I: the Linear Setting

Fully implicit Runge--Kutta (IRK) methods have many desirable properties as time integration schemes in terms of accuracy and stability, but high-order IRK methods are not commonly used in practice with numerical PDEs due to the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic preconditioning framework for solving the systems of equations that arise from IRK methods applied to linear numerical PDEs (without algebraic constraints). Additionally, this framework also naturally applies to discontinuous Galerkin discretizations in time. Under quite general assumptions on the spatial discretization that yield stable time integration, the preconditioned operator is proven to have condition number bounded by a small, order-one constant, independent of the spatial mesh and time-step size, and with only weak dependence on number of stages/polynomial order; for example, the preconditioned operator for 10th-order Gauss IRK has condition number less than two, independent of the spatial discretization and time step. The new method can be used with arbitrary existing preconditioners for backward Euler-type time-stepping schemes and is amenable to the use of three-term recursion Krylov methods when the underlying spatial discretization is symmetric. The new method is demonstrated to be effective on various high-order finite-difference and finite element discretizations of linear parabolic and hyperbolic problems, demonstrating fast, scalable solution of up to 10th-order accuracy. The new method consistently outperforms existing block preconditioning approaches, and in several cases, the new method can achieve 4th-order accuracy using Gauss integration with roughly half the number of preconditioner applications and wallclock time as required using standard diagonally IRK methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Efficient Synchronization Stability Metrics for Fault Clearing

Direct methods can provide rapid screening of the dynamical security of large numbers fault and contingency scenarios by avoiding extensive time simulation. We introduce a computationally-efficient direct method based on optimization that leverages efficient cutting plane techniques. The method considers both unstable equilibrium points and the effects of additional relay tripping on dynamical security[1]. Similar to other direct methods, our approach yields conservative results for dynamical security, however, the optimization formulation potentially lends itself to the inclusion of additional constraints to reduce this conservatism.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

Electronic representation of wave equation

The Taylor series method for solving differential equations represents a non-traditional way of a numerical solution. Even though this method is not much preferred in the literature, experimental calculations done at the Department of Intelligent Systems of the Faculty of Information Technology of TU Brno have verified that the accuracy and stability of the Taylor series method exceeds the currently used algorithms for numerically solving differential equations. This paper deals with solution of Telegraph equation using modelling of a series small pieces of the wire. Corresponding differential equations are solved by the Modern Taylor Series Method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

SU-E-T-410: Fringe Stability and Phase Shift Measurements in a Michelson Interferometer for Optical Calorimetry

Purpose: To identify the variables limiting the resolution of a Michelson interferometer used to measure phase shifts (PS) in water as part of a radiometric calorimeter. Methods: We investigated the output stability of a He-Ne laser and a laser diode. The short and long term stability of the fringe pattern in a Michelson interferometer was tested with different types of lasers, thermal insulation arrangements, damping systems and optical mounts to optimize system performance. PS were induced by electrically heating water in a 1 cm quartz cuvette located in one of the interferometer arms. The PS was calculated from fringe intensity changes and compared to a calculated PS using thermocouple-measured temperature changes in the water. Results: The intensity of the laser diode is more stable, but the gas laser’s profile is more suitable for fringe analysis and has better temporal coherence. The laser requires a warm-up time of 4 hours before its output is stabilized (SNR>95). The fringe’s stability strongly depends on the thermal insulation. When the interferometer is exposed to ambient temperature swings of 0.7 K, it is not possible to stabilize the fringe pattern. Enclosing the system in a 2.5 cm-thick Styrofoam box improves the SNR, but further insulation will be needed to increase the SNR above 50. High frequency noise is significantly reduced by damping the system.Inducing a temperature rise in water, starting at 299 K, the average temperature increase for a 2π PS is 0.29 ± 0.02 K and the proportionality constant is -21.1 ± 0.8 radians/K. This is 5.8% lower than the calculated value using the thermocouple. Conclusion: Interferometric PS measurements of temperature may provide an alternative to thermistors for water calorimetry. The resolution of the current prototype is limited by ambient temperature stability. Calculated and measured thermally-induced PS in water agreed to within 5.8%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Estimation of turbulence dissipation rate and its variability from sonic anemometer and wind Doppler lidar during the XPIA field campaign

Despite turbulence being a fundamental transport process in the boundary layer, the capability of current numerical models to represent it is undermined by the limits of the adopted assumptions, notably that of local equilibrium. Here we leverage the potential of extensive observations in determining the variability in turbulence dissipation rate (ϵ). These observations can provide insights towards the understanding of the scales at which the major assumption of local equilibrium between generation and dissipation of turbulence is invalid. Typically, observations of ϵ require time- and labor-intensive measurements from sonic and/or hot-wire anemometers. We explore the capability of wind Doppler lidars to provide measurements of ϵ. We refine and extend an existing method to accommodate different atmospheric stability conditions. To validate our approach, we estimate ϵ from four wind Doppler lidars during the 3-month XPIA campaign at the Boulder Atmospheric Observatory (Colorado), and we assess the uncertainty of the proposed method by data intercomparison with sonic anemometer measurements of ϵ. Our analysis of this extensive dataset provides understanding of the climatology of turbulence dissipation over the course of the campaign. Further, the variability in ϵ with atmospheric stability, height, and wind speed is also assessed. Finally, we present how ϵ increases as nocturnal turbulence is generated during low-level jet events.

17 WIND ENERGY↗
DOE OSTI2016

Computation of resistive instabilities by matched asymptotic expansions

Here, we present a method for determining the linear resistive magnetohydrodynamic (MHD) stability of an axisymmetric toroidal plasma, based on the method of matched asymptotic expansions. The plasma is partitioned into a set of ideal MHD outer regions, connected through resistive MHD inner regions about singular layers where q = m/n, with m and n toroidal mode numbers, respectively, and q the safety factor. The outer regions satisfy the ideal MHD equations with zero-frequency, which are identical to the Euler-Lagrange equations for minimizing the potential energy delta W. The solutions to these equations go to infinity at the singular surfaces. The inner regions satisfy the equations of motion of resistive MHD with a finite eigenvalue, resolving the singularity. Both outer and inner regions are solved numerically by newly developed singular Galerkin methods, using specialized basis functions. These solutions are matched asymptotically, providing a complex dispersion relation which is solved for global eigenvalues and eigenfunctions in full toroidal geometry. The dispersion relation may have multiple complex unstable roots, which are found by advanced root-finding methods. These methods are much faster and more robust than the previous numerical methods. The new methods are applicable to more challenging high-pressure and strongly shaped plasma equilibria and generalizable to more realistic inner region dynamics. In the thermonuclear regime, where the outer and inner regions overlap, they are also much faster and more accurate than the straight-through methods, which treat the resistive MHD equations in the whole plasma volume.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2016

Stability of rigidly rotating supermassive stars against gravitational collapse

We revisit secular stability against quasi-radial collapse for rigidly rotating supermassive stars (SMSs) in general relativity. We suppose that the SMSs are in a nuclear-burning phase and can be modeled by polytropic equations of state with the polytropic index n{sub p} slightly smaller than 3. The stability is determined in terms of the turning point method. We find a fitting formula of the stability condition for the plausible range of n{sub p} (2.95≲n{sub p}≲3) for SMSs. This condition reconfirms that while non-rotating SMSs with a mass of ∼10{sup 5}M{sub ⊙}–10{sup 6}M{sub ⊙} may undergo a general relativistically induced quasi-radial collapse, rigidly rotating SMSs with a ratio of rotational to gravitational potential energy (β) of ∼10{sup −2} are likely to be stable against collapse unless they are able to accrete ∼5 times more mass during the (relatively brief) hydrogen-burning phase of their evolution. We discuss the implications of our results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

A Comparative Study of Two Widely Used Grid-Forming Droop Controls on Microgrid Small Signal Stability

Throughout history, two similar grid-forming droop controls are widely reported in literature—the single-and multi-loop droop controls. Although being very similar, the authors find that the dynamic performance and stability characteristics of each control method are very different in a microgrid. Compared with single-loop droop control, the multi-loop droop control is prone to be less damped and loses stability more easily under some circumstances. This paper provides a novel insight into the different dynamic responses of the two basic controls. It points out that the two similar controls adjust the angular frequency and voltage magnitude at different locations within the inverter, resulting in different coupling reactances that impact the dynamic response and stability of microgrids differently. The use of single-loop droop control results in a larger coupling reactance, which helps to improve dynamic response and stability. This novel insight is verified through full-order small signal analysis, offline electromagnetic transient simulation, and real-time hardware-in-the-loop simulation experiments. Results show that the microgrid has a larger small signal stability boundary when using single-loop droop control, and this difference increases as the value of an inverter’s inner filter inductance increases.

42 ENGINEERING↗
DOE OSTI2023

Optimizing cloud motion estimation on the edge with phase correlation and optical flow

Abstract. Phase correlation (PC) is a well-known method for estimating cloud motion vectors (CMVs) from infrared and visible spectrum images. Commonly, phase shift is computed in the small blocks of the images using the fast Fourier transform. In this study, we investigate the performance and the stability of the blockwise PC method by changing the block size, the frame interval, and combinations of red, green, and blue (RGB) channels from the total sky imager (TSI) at the United States Atmospheric Radiation Measurement user facility's Southern Great Plains site. We find that shorter frame intervals, followed by larger block sizes, are responsible for stable estimates of the CMV, as suggested by the higher autocorrelations. The choice of RGB channels has a limited effect on the quality of CMVs, and the red and the grayscale images are marginally more reliable than the other combinations during rapidly evolving low-level clouds. The stability of CMVs was tested at different image resolutions with an implementation of the optimized algorithm on the Sage cyberinfrastructure test bed. We find that doubling the frame rate outperforms quadrupling the image resolution in achieving CMV stability. The correlations of CMVs with the wind data are significant in the range of 0.38–0.59 with a 95 % confidence interval, despite the uncertainties and limitations of both datasets. A comparison of the PC method with constructed data and the optical flow method suggests that the post-processing of the vector field has a significant effect on the quality of the CMV. The raindrop-contaminated images can be identified by the rotation of the TSI mirror in the motion field. The results of this study are critical to optimizing algorithms for edge-computing sensor systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Stability-Constrained Learning for Frequency Regulation in Power Grids With Variable Inertia

The increasing penetration of converter-based renewable generation has resulted in faster frequency dynamics, and low and variable inertia. As a result, there is a need for frequency control methods that are able to stabilize a disturbance in the power system at timescales comparable to the fast converter dynamics. This paper proposes a combined linear and neural network controller for inverter-based primary frequency control that is stable at time-varying levels of inertia. We model the time-variance in inertia via a switched affine hybrid system model. We derive stability certificates for the proposed controller via a quadratic candidate Lyapunov function. We test the proposed control on a 12-bus 3-area test network, and compare its performance with a base case linear controller, optimized linear controller, and finite-horizon Linear Quadratic Regulator (LQR). Our proposed controller achieves faster mean settling time and over 50% reduction in average control cost across 100 inertia scenarios compared to the optimized linear controller. Unlike LQR which requires complete knowledge of the inertia trajectories and system dynamics over the entire control time horizon, our proposed controller is real-time tractable, and achieves comparable performance to LQR.

data-driven control↗
DOE OSTI2021

Implicit extensions of an explicit multirate Runge–Kutta scheme

In this work, we propose a new method that extends conservative explicit multirate methods to implicit explicit-multirate methods. We develop extensions of order one and two with different stability properties on the implicit side. The method is suitable for time-stepping adaptive mesh refinement PDE discretizations with different degrees of stiffness. A numerical example with an advection-diffusion problem illustrates the new method's properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Humidity-Induced Photoluminescence Hysteresis in Variable Cs/Br Ratio Hybrid Perovskites

Hybrid organic–inorganic perovskites containing Cs are a promising new material for light-absorbing and light-emitting optoelectronics. However, the impact of environmental conditions on their optical properties is not fully understood. In this paper, we elucidate and quantify the influence of distinct humidity levels on the charge carrier recombination in Cs x FA 1–x Pb(I y Br 1–y ) 3 perovskites. Using in situ environmental photoluminescence (PL), we temporally and spectrally resolve light emission within a loop of critical relative humidity (rH) levels. Our measurements show that exposure up to 35% rH increases the PL emission for all Cs (10–17%) and Br (17–38%) concentrations investigated here. Spectrally, samples with larger Br concentrations exhibit PL redshift at higher humidity levels, revealing water-driven halide segregation. The compositions considered present hysteresis in their PL intensity upon returning to a low-moisture environment due to partially reversible hydration of the perovskites. Our findings demonstrate that the Cs/Br ratio strongly influences both the spectral stability and extent of light emission hysteresis. Finally, we expect our method to become standard when testing the stability of emerging perovskites, including lead-free options, and to be combined with other parameters known for affecting material degradation, e.g., oxygen and temperature.

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

Templating Routes to Supported Oxide Catalysts by Design

The rational design and understanding of supported oxide catalysts requires at least three advancements, in order of increasing complexity: the ability to quantify the number and nature of active sites in a catalytic material, the ability to place external controls on the number and structure of these active sites, and the ability to assemble these active sites so as to carry out more complex functions in tandem. As part of an individual investigator research program that is integrated with the Northwestern University Institute for Catalysis in Energy Processes (ICEP) as of 2015, significant advances were achieved in these three areas. First, phosphonic acids were utilized in the quantitative assessment of the number of active and geometrically-available sites in MO x -SiO 2 catalysts, including nanocrystalline composites, co-condensed materials, and grafted structures, for M=Ti, Zr, Hf, Nb, and Ta. That work built off progress in understanding supported Fe, Cu, and Co oxide catalysts from chelating and/or multinuclear precursors to maximize surface reactivity. Secondly, significant progress was made in the new area of using thin oxide overcoats containing ‘nanocavities’ from organic templates as a method to control the dispersion and thermal stability of subsequently deposited metal nanoparticles or other catalytic domains. Similar methods were used to control surface reactivity in SiO 2 -Al 2 O 3 acid catalysts and to control reactant selectivity in Al 2 O 3 -TiO 2 photocatalysts. Finally, knowledge gained from the first two areas has been combined to synthesize a tandem catalyst for hydrotreating reactions and an orthogonal tandem catalyst system where two subsequent reactions in a reaction network are independently controlled by light and heat. Overall, work carried out under this project significantly advanced the knowledge of synthesis-structure-function relationships in supported oxide catalysts for energy applications.

02 PETROLEUM↗
DOE OSTI2019

Hybrid Fixed-Point Fixed-Stress Splitting Method for Linear Poroelasticity

Efficient and accurate poroelasticity models are critical in modeling geophysical problems such as oil exploration, gas-hydrate detection, and hydrogeology. We propose an efficient operator splitting method for Biot’s model of linear poroelasticity based on fixed-point iteration and constrained stress. In this method, we eliminate the constraint on time step via combining our fixed-point approach with a physics-based restraint between iterations. Three different cases are considered to demonstrate the stability and consistency of the method for constant and variable parameters. The results are validated against the results from the fully coupled approach. In case I, a single iteration is used for continuous coefficients. The relative error decreases with an increase in time. In case II, material coefficients are assumed to be linear. In the single iteration approach, the relative error grows significantly to 40% before rapidly decaying to zero. This is an artifact of the approximate solutions approaching the asymptotic solution. The error in the multiple iterations oscillates within 10 − 6 before decaying to the asymptotic solution. Nine iterations per time step are enough to achieve the relative error close to 10 − 7 . In the last case, the hybrid method with multiple iterations requires approximately 16 iterations to make the relative error 5 × 10 − 6 .

Delgado, Paul M.↗
DOE OSTI2018

STABILIZATION OF BIO-OIL TO ENABLE ITS HYDROTREATING TO PRODUCE BIOFUELS

Fast pyrolysis is considered to be the simplest and most cost-effective approach to produce liquid oil (bio-oil) from biomass. Bio-oil is not suitable to substitute for petroleum as high-quality fuels and significant upgrading such as hydrotreating is required to remove oxygen, add hydrogen, and rearrange carbon backbone of bio-oil. However, the grand challenge in bio-oil hydrotreating technology is bio-oil instability, which limits the lifetime of catalyst and operation. To enable a sustainable and economically viable process for bio-oil hydrotreating, it is vital to develop effective technologies for stabilizing bio-oils. This chapter will be devoted to bio-oil stabilization. The current understating of the major cause of bio-oil instability, condensation of reactive species such as sugar, aldehydes, ketones, and phenolics, is elucidated. The reported physical and chemical methods for bio-oil stabilization are summarized in detail, with a specific focus on bio-oil catalytic hydrogenation for stabilization. The impact of stabilization on bio-oil hydrotreating is discussed as well.

Wang, Huamin↗