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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Machine learning method for enforcing variable independence in background estimation with LHC data: ABCDisCoTEC

A novel solution is presented for the problem of estimating the backgrounds of a signal search using observed data while simultaneously maximizing the sensitivity of the search to the signal. The 'ABCD method' provides a reliable framework for background estimation by partitioning events into one signal-enhanced region (A) and three background-enhanced control regions (B, C, and D) via two smoothly varying, statistically independent variables. In practice, even slight correlations between the two variables can significantly undermine the method's performance. Thus, choosing appropriate variables by hand can present a formidable challenge, especially when background and signal differ only subtly. To address this issue, the ABCD with distance correlation (ABCDisCo) method was developed to construct two learned variables via a neural network trained to provide strong signal-background discrimination with small values of the distance correlation (DisCo) measure between the two learned variables. However, relying solely on minimizing the DisCo can result in learned variables that may not have distributions of background events that are smoothly varying and localized at extreme values, as necessary for the validity of the background estimation. The ABCDisCo training enhanced with closure (ABCDisCoTEC) method is introduced to solve this issue by directly minimizing the nonclosure, expressed as a dedicated differentiable loss term. This extended method is applied to a data set of proton-proton collisions at a center-of-mass energy of 13 TeV recorded by the CMS detector at the CERN Large Hadron Collider. Additionally, given the complexity of the minimization problem with constraints on multiple loss terms, the modified differential method of multipliers is applied and shown to greatly improve the stability and robustness of the ABCDisCoTEC method, compared to grid search hyperparameter optimization procedures.

Hayrapetyan, Aram [Yerevan Phys. Inst.]↗

DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs

Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 x 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 x 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier--Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier--Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Direct implicit and explicit energy-conserving particle-in-cell methods for modeling of capacitively coupled plasma devices

Achieving large-scale kinetic modeling is a crucial task for the development and optimization of modern plasma devices. With the trend of decreasing pressure in applications, such as plasma etching, kinetic simulations are necessary to self-consistently capture the particle dynamics. The standard, explicit, electrostatic, momentum-conserving particle-in-cell method suffers from restrictive stability constraints on spatial cell size and temporal time step, requiring resolution of the electron Debye length and electron plasma period, respectively. This results in a very high computational cost, making the technique prohibitive for large volume device modeling. We investigate the direct implicit algorithm and the explicit energy conserving algorithm as alternatives to the standard approach, both of which can reduce computational cost with a minimal (or controllable) impact on results. These algorithms are implemented into the well-tested EDIPIC-2D and LTP-PIC codes, and their performance is evaluated via 2D capacitively coupled plasma discharge simulations. The investigation reveals that both approaches enable the utilization of cell sizes larger than the Debye length, resulting in a reduced runtime, while incurring only minor inaccuracies in plasma parameters. The direct implicit method also allows for time steps larger than the electron plasma period; however, care must be taken to avoid numerical heating or cooling. It is demonstrated that by appropriately adjusting the ratio of cell size to time step, it is possible to mitigate this effect to an acceptable level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Explicit Monotone Stable Super-Time-stepping Methods for Finite Time Singularities

We explore a novel way to numerically resolve the scaling behavior of finite-time singularities in solutions of nonlinear parabolic PDEs. The Runge–Kutta–Legendre (RKL) and Runge–Kutta–Gegenbauer (RKG) super-time-stepping methods were originally developed for nonlinear complex physics problems with diffusion. These are multistage single step second-order, forward-in-time methods with no implicit solves. The advantage is that the time-step size for stability scales with stage number 𝑠 as $\mathcal{O}$⁡(𝑠 2 ). Many interesting nonlinear PDEs have finite-time singularities, and the presence of diffusion often limits one to using implicit or semi-implicit time-step methods for stability constraints. Finite-time singularities are particularly challenging due to the large range of scales that one desires to resolve, often with adaptive spatial grids and adaptive time steps. Here, in this study, we show two examples of nonlinear PDEs for which the self-similar singularity structure has time and space scales that are resolvable using the RKL and RKG methods, without forcing even smaller time steps. Compared to commonly used implicit numerical methods, we achieve a significantly smaller run time while maintaining comparable accuracy. We also prove numerical monotonicity for both the RKL and RKG methods under their linear stability conditions for the constant coefficient heat equation, in the case of infinite domain and periodic boundary condition, leading to a theoretical guarantee of the superiority of the RKL and RKG methods over traditional super-time-stepping methods, such as the Runge-Kutta-Chebyshev and the orthogonal Runge-Kutta-Chebyshev methods. Code can be found at https://github.com/ZT220501/SRK-Singularity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Three-dimensional flow model development for thermal mixing and stratification modeling in reactor system transients analyses

Mixing, thermal-stratification, and mass transport phenomena in large pools or enclosures play major roles for the safety of reactor systems. Depending on the fidelity requirement and computational resources, various modeling methods, from the 0-D perfect mixing model to 3-D Computational Fluid Dynamics (CFD) models, are available. Each is associated with its own advantages and shortcomings. It is very desirable to develop an advanced and efficient thermal mixing and stratification modeling capability embedded in a modern system analysis code to improve the accuracy of reactor safety analyses and to reduce modeling uncertainties. An advanced system analysis tool, SAM, is being developed at Argonne National Laboratory for advanced non-LWR reactor safety analysis. While SAM is being developed as a system-level modeling and simulation tool, a reduced-order three-dimensional module is under development to model the mull-dimensional flow and thermal mixing and stratification in large enclosures of reactor systems. This paper provides an overview of the three-dimensional finite element flow model in SAM, including the governing equations, stabilization scheme, and solution methods. Additionally, several verification and validation tests are presented, including lid-driven cavity flow, natural convection inside a cavity, laminar flow in a channel of parallel plates. Here, based on the comparisons with the analytical solutions and experimental results, it is demonstrated that the developed 3-D fluid model can perform very well for a wide range of flow problems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1972

On the stability of numerical integration routines.

Numerical integration methods for the solution of initial value problems for ordinary vector differential equations may be modelled as discrete time feedback systems. The stability criteria discovered in modern control theory are applied to these systems and criteria involving the routine, the step size and the differential equation are derived. Linear multistep, Runge-Kutta, and predictor-corrector methods are all investigated.

Glover, K.↗

DOE OSTI2024

Dynamic flux surrogate-based partitioned methods for interface problems

Loosely coupled partitioned methods for multiphysics problems treat each subproblem as a separate entity and advance them independently in time. In so doing these methods enable code reuse, increase concurrency and provide a convenient framework for plug-and-play multiphysics simulations. However, mathematically loosely coupled schemes are equivalent to a single step of an iterative solution method, which can compromise their accuracy and stability. We present a new data-driven partitioned method for coupled parametric PDEs that can improve upon the accuracy of traditional loosely coupled methods without incurring a performance penalty. To that end, we replace conventional field transfers across the interface by a surrogate for the dynamics of the interface flux exchanged between the subdomains. To develop this surrogate we apply dynamic mode decomposition to a non-standard staggered-in-time state, comprising the interface flux and small solution patches near the interface. The new approach shifts the main computational burden to an offline training phase, whereas application of the surrogate in the online phase amounts to a single matrix–vector multiplication. In conclusion, we provide stability analysis of the surrogate-based partitioned scheme and include numerical results that demonstrate its potential.

Dynamic mode decomposition (DMD)↗

NASA NTRS1953

Wind-tunnel investigation of the effects of various dorsal-fin and vertical-tail configurations on the directional stability of a streamlined body at transonic speeds : transonic-bump method

Yawing-moment coefficients were obtained for several dorsal-fin and vertical-tail configurations in combination with a streamlined body through a large angle-of-sideslip range and a Mach number range of 0.59 to 1.11. The results indicated that dorsal fins improved the directional stability characteristics of the body alone and the body-vertical-tail configuration. A ring tail was more effective at small angles of sideslip and less effective at large angles of sideslip than a tapered low-aspect-ratio vertical tail.

TAIL-BODY COMBINATIONS - MISSILES↗

DOE OSTI2023

Development and Stability Analysis of Representative Future High Renewable Penetration Power Grid Models for California

The power grid in California is now experiencing a significant increase in renewable generation to achieve 100% clean energy by 2040. This paper discusses a method for developing the transient stability (TS) model of the future high renewable penetration in the power grid in California. The method consists of utilizing an existing future Western Electricity Coordinating Council (WECC) grid model in the database and locating the conventional plants in the base model that will potentially retire in the next few years in California. Next, the dynamic models of those identified conventional plants and controls in the existing model are modified to replicate the future power generation in the California region. Different scenarios of future power grid in California are developed for various renewable penetration levels and control strategies in renewable power plants. The effect of the control strategies on the stability of the power grid under events like loss of generation are evaluated. Based on these evaluations, certain control features in power plants are identified as necessary for stability of future power grids.

Samanta, Sayan↗

NASA NTRS1995

Euler/Navier-Stokes flow computations on flexible configurations for stability analysis

Longitudinal dynamic stability derivatives required for design of aircraft are computed by using the state-of-the-art numerical methods for wing-body configurations. The flow is modeled using the Euler/Navier-Stokes equations with turbulence models and solved using an efficient finite-difference scheme suitable for patched structured grids. Computations are made at a flow regime that is beyond the limits of the current linear methods mostly used for computing stability derivatives. Flow conditions include shockwaves and viscous dominated vortical flows. Effect of Mach number and angle-of-attack on stability derivatives are demonstrated for a typical wing-body configuration. For the same configuration the effects of wing flexibility on the magnitude and phase angles of stability derivatives are also demonstrated.

Guruswamy, G.↗

NASA NTRS1995

Euler/Navier-Stokes Flow Computations on Flexible Configurations for Stability Analysis

Abstract Longitudinal dynamic stability derivatives required for design of aircraft are computed by using the state-of-the-art numerical methods for wing-body configurations. The flow is modeled using the Euler/Navier-Stokes equations with turbulence models and solved using an efficient finite-difference scheme suitable for patched structured grids. Computations are made at a flow regime that is beyond the limits of the current linear methods mostly used for computing stability derivatives. Flow conditions include shock-waves and viscous dominated vortical flows. Effect of Mach number and angle-of-attack on stability derivatives are demonstrated for a typical wing-body configuration. For the same configuration the effects of wing flexibility on the magnitude and phase angles of stability derivatives are also demonstrated.

Guruswamy, Guru P.↗

NASA NTRSDate not supplied

TPSAS-NF1676L-10737-DND

The Clouds and the Earth's Radiant Energy System (CERES) project provides the community a quality climate TOA and surface fluxes and cloud properties and are key elements in evaluating the ability of climate models to predict both past and future climate states. One important aspect of successful climate monitoring is satellite calibration. The CERES instantaneous broadband calibration radiance stability is 0.1% per year, whereas the accuracy of the longwave (LW) and shortwave (SW) is better than 0.5% and 1% respectively. CERES uses the 3-hourly GEO fluxes and clouds to provide diurnally averaged daily and monthly TOA flux means from the CERES measurements onboard the Terra (10:30 AM local equator crossing time) and Aqua (13:30 PM) sun-synchronous satellites. The GEO radiances are calibrated against MODIS as a reference to provide consistent GEO retrieved cloud properties and derived broadband fluxes. These fluxes are then carefully normalized against the well-calibrated CERES fluxes to ensure climate quality fluxes, free of GEO artifacts. With the upcoming release of Edition3 products, CERES has just recalibrated all the CERES instruments to ensure that both Terra and Aqua fluxes are on the same radiometric scale and the removal of all known instrument trends. CERES will also recalibrate all GEO satellites using the entire GEO time record. CERES employs collocated, coincident ray-matched GEO/MODIS radiances as its primary technique, calibrating GEO radiances beginning in the year 2000. The advantage of the ray-matching technique is that it provides a full dynamic range of radiances to check the linearity of the GEO visible instrument. CERES uses deep convective clouds and invariant desert sites to validate the primary technique. All three methods can assess the stability of the GEO visible instrument. Examples of the CERES GEO calibration techniques and validation efforts will be shown at the conference. Special emphasis will be given to the newly incorporated desert site methodology.

R. Bhatt↗

DOE OSTI2026

Excipient screening by lyophilization provides insights into spray drying formulations for nanoparticle vaccines

Nanoparticles have shown great promise as delivery platforms in the development of tunable and safe vaccines. Nanolipoprotein particles (NLPs), also known as nanodiscs, are discoidal nanoparticles composed of a lipid bilayer stabilized at their periphery by apolipoproteins. Under the right conditions, the NLP self-assembly process is highly customizable in terms of lipids and apolipoprotein constituents, allowing for tunable physical and chemical characteristics. This flexibility allows a wide range of vaccine antigens and adjuvants to be incorporated onto the NLP platform for tailored vaccine design. The stability of NLPs during long term storage is a very important factor in developing a vaccine delivery platform suitable for widespread global use. When stored in a solution for extended periods of time, NLPs dissociate into their corresponding lipids and protein constituents, leading to particle degradation. Proper stabilization of NLPs can often be achieved by lyophilization (i.e. freeze-drying), a method widely used for various applications including pharmaceuticals. This process, however, can be damaging to particles without the presence of lyoprotectants or excipients that help maintain particle stability during lyophilization. Another method used to stabilize vaccines and pharmaceuticals is spray drying, a process that converts liquid formulations into dry powders through controlled heating and airflow. While spray drying is rapid, scalable, and cost-effective, lyophilization is typically a gentler process that better retains biomolecule structure and function. Both processes use excipients for particle stabilization, so lyophilization can be used as a surrogate to test stability of NLPs, to down-select formulations that may withstand the harsher conditions of spray drying. To screen different formulations, NLPs were synthesized and purified to homogeneity by size exclusion chromatography (SEC) and samples were prepared with a wide range of lyoprotectants and/or excipients. To assess the protective effects of excipients on NLPs upon spray drying, both pre- and post-lyophilized samples were analyzed by SEC. To assess the protective effects upon heating (encountered during the spray drying process), NLP samples were incubated at elevated temperatures prior to SEC analysis. The lyoprotectants and excipients evaluated in this study had different efficiencies in protecting NLPs during lyophilization and heating tests. Trehalose, for example, exhibits stabilization on NLPs both upon lyophilization and heating whereas leucine accelerated NLP dissociation. Although some lyoprotectants are effective by themselves, different combinations can decrease the stabilization of NLPs. Shelf-stable vaccines that do not require cold-chain storage are essential for global accessibility and our findings provide fundamental insight into how to advance NLP-based vaccines for these applications.

Serrano, Litzay J [Lawrence Livermore National Lab↗

DOE OSTI2021

Reactive Flows in Porous Media: Challenges in Theoretical and Numerical Methods

We review theoretical and computational research, primarily from the past 10 years, addressing the flow of reactive fluids in porous media. The focus is on systems where chemical reactions at the solid–fluid interface cause dissolution of the surrounding porous matrix, creating nonlinear feedback mechanisms that can often lead to greatly enhanced permeability. Here, we discuss insights into the evolution of geological forms that can be inferred from these feedback mechanisms, as well as some geotechnical applications such as enhanced oil recovery, hydraulic fracturing, and carbon sequestration. Until recently, most practical applications of reactive transport have been based on Darcy-scale modeling, where averaged equations for the flow and reactant transport are solved. We summarize the successes and limitations of volume averaging, which leads to Darcy-scale equations, as an introduction to pore-scale modeling. Pore-scale modeling is computationally intensive but offers new insights as well as tests of averaging theories and pore-network models. We include recent research devoted to validation of pore-scale simulations, particularly the use of visual observations from microfluidic experiments.

wormhole formation↗

DOE OSTI2022

Stabilization and control of persistent current magnets using variable inductance

Ultra-stable, tunable magnetic fields are desirable for a wide range of applications in medical imaging, electron microscopy, quantum science, and atomic physics. Superconducting magnets operated in persistent current mode, with device current flowing in a closed superconducting loop disconnected from a power source, are a common approach for applications with the most stringent requirements on temporal field stability. In this work, we present a method for active control of this persistent current by means of dynamic inductance change within the superconducting circuit. For a first realization of this general technique, we consider a variable superconducting inductor placed in series with the main magnet. The inductor acts as a dynamic flux storage device capable of transferring flux to or from the main magnet through inductance change. This allows for fine and fast adjustments of the persistent current without the use of thermal switches that limit the speed and accuracy of many present-day methods. With first experiments employing this technique, we demonstrate stabilization of a 1.95 T Nb–Ti round lens for electron microscopy against decay resulting from residual losses in the superconducting circuit, and more generally show flexibility for precise control over the magnitude and waveform of the persistent current.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1966

Calculated blade response at high tip-speed ratios

This paper presents the initial results of a digital computer study of helicopter rotor-blade-motion stability for a broad range of forward speeds encompassing proposed conventional and compound helicopter designs. The analysis treats the general case of nonlinear, coupled flapping and lagging motions of hinged rigid blades. The results of the study to date indicate that blade-motion stability boundaries cannot be specified for a given blade design in terms of a fixed value of rotor tip-speed ratio. The stability boundaries can shift significantly within the desired operational speed range, with the stability limits depending upon the rotor loading, the magnitude of the disturbance encountered by the blade, and the blade position at the instant the external disturbance is encountered. The sample methods suggested for improving blade-motion stability include increasing the effective hinge spring restraint by incorporating pitch-flap coupling or by using hingeless or teetering rotor systems. Reducing rotor loading also has a beneficial effect on the blade-motion stability. The results indicate that methods used to deal with the blade-motion stability problem can be expected to diminish the overall vibration levels of the helicopter rotor system across the entire design speed range.

TIP SPEED↗

NASA NTRS1988

Stability analysis and control of tethered satellites

Issues of stability and stabilization of the Shuttle/Tethered Satellite System (TSS) are considered. The basic, nonlinear equations of motion of the TSS are derived using the system Lagrangian. A tension control strategy in feedback form is introduced to stabilize the TSS during station-keeping. In addition, the stability properties of a constant angle control method for deployment and retrieval are investigated using Liapunov stability theory.

Liaw, D.-C.↗