Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability and convergence”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Efficient CP Rounding Using Alternating Least Squares with QR Decomposition

The CANDECOMP/PARAFAC (CP) decomposition is widely used for analyzing multidimensional data, and the alternating least squares (CP-ALS) algorithm is a common method for its computation. CP rounding is the problem of computing a lower-rank CP decomposition of an input already in a higher-rank CP format. While the normal equations (NE) approach in CP-ALS is efficient for the CP rounding problem and frequently used, it becomes unstable in the presence of ill-conditioned subproblems. This paper presents a new QR-based CP-ALS method for CP rounding that preserves both numerical stability and computational efficiency. Here, our experiments show that the proposed method offers significant speedup over a previous QR-based approach and the Tensor Toolbox's NE-based implementation, particularly for higher-order tensors. Furthermore, our approach demonstrates a marked reduction in error for ill-conditioned problems, with error reductions several orders of magnitude smaller compared to the NE-based method, while achieving faster convergence and more accurate solutions. By using a more numerically stable approach, we can solve more problems in reduced working precision, which enables further reduction in time to solution.

CANDECOMP/PARAFAC↗
DOE OSTI2026

Assessment of the hydromechanical higher-order MPM for the simulation of geotechnical problems

The Material Point Method (MPM) has been increasingly used to simulate large strain deformations. Linear interpolation functions are commonly used to perform the spatial integration. It is well-known that the discontinuities in the interpolation function derivatives induce shock-like artifacts known as ‘cell-crossing’ error. These errors compound with volumetric locking errors when used with hydromechanical formulations for porous media, where different velocity fields are used for each phase. The capabilities of higher-order MPM frameworks have not been explored for real-scale geotechnical problems. As such, this paper aims to assess, validate, and further discuss a higher-order B-spline MPM (BS-MPM) framework. First, the BS-MPM framework is verified against the large-strain oedometer consolidation problem. Second, the framework is validated against a real-scale slope failure experiment triggered by pore water pressure recharge. Landslide features that are captured using the higher-order framework are specifically highlighted, and results (e.g., pore water pressure and deformation) are validated with field measurements. A generally convergent numerical solution is observed when using cubic interpolation functions. Third, a footing penetration problem is simulated using the multi-patch BS-MPM. Trends are examined with respect to penetration velocity and variation in hydraulic conductivity. The BS-MPM framework ultimately presents a stabilized numerical solution that captures plausible hydromechanical interaction trends important in geotechnical engineering applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Predicting and Synthesizing Interface Stabilized 2D Layers

The compound (Pb 2 MnSe 3 ) 0.6 VSe 2 was predicted to be kinetically stable based on density functional theory (DFT) calculations on an island of Pb 2 MnSe 3 between layers of VSe 2 . This approach provides a high degree of freedom by not forcing interlayer lattice match, making it ideal to investigate the likelihood of formation of new incommensurate layer misfit structures. The free space around the island is critical, as it allows atoms to diffuse and hence exploring the local energy landscape around the initial configuration. (Pb 2 MnSe 3 ) 0.6 VSe 2 was synthesized via a near diffusionless reaction from precursors where a repeating sequence of elemental layers matches the local composition and layer sequence of the predicted compound. The VSe 2 layer consists of a Se–V–Se trilayer with octahedral coordination of the V atoms. The Pb 2 MnSe 3 layer consists of three rock-salt-like planes, with a MnSe layer between the planes of PbSe. The center MnSe plane stabilizes the puckering of the outer PbSe layers. Electrical properties indicate that (Pb 2 Mn 1 Se 3 ) 0.6 VSe 2 undergoes a charge density wave transition at ~100 K and orders ferromagnetically at 35 K. Overall, the combination of theory and experiment enables a faster convergence to new heterostructures than either approach in isolation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Learning Provably Stable Local Volt/Var Controllers for Efficient Network Operation

Here this paper develops a data-driven framework to synthesize local Volt/Var control strategies for distributed energy resources (DERs) in power distribution grids (DGs). Aiming to improve DG operational efficiency, as quantified by a generic optimal reactive power flow (ORPF) problem, we propose a two-stage approach. The first stage involves learning the manifold of optimal operating points determined by an ORPF instance. To synthesize local Volt/Var controllers, the learning task is partitioned into learning local surrogates (one per DER) of the optimal manifold with voltage input and reactive power output. Since these surrogates characterize efficient DG operating points, in the second stage, we develop local control schemes that steer the DG to these operating points. We identify the conditions on the surrogates and control parameters to ensure that the locally acting controllers collectively converge, in a global asymptotic sense, to a DG operating point agreeing with the local surrogates. We use neural networks to model the surrogates and enforce the identified conditions in the training phase. AC power flow simulations on the IEEE 37-bus network empirically bolster the theoretical stability guarantees obtained under linearized power flow assumptions. The tests further highlight the optimality improvement compared to prevalent benchmark methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Deep-learning-aided extraction of optical constants in scanning near-field optical microscopy

Scanning near-field optical microscopy is one of the most effective techniques for spectroscopy of nanoscale systems. However, inferring optical constants from the measured near-field signal can be challenging because of a complicated and highly nonlinear interaction between the scanned probe and the sample. Conventional fitting methods applied to this problem often suffer from the lack of convergence or require human intervention. Here, we develop an alternative approach where the optical parameter extraction is automated by a deep learning network. The network provides an initial estimate that is subsequently refined by a traditional fitting algorithm. In conclusion, we show that this method demonstrates superior accuracy, stability against noise, and computational speed when applied to simulated near-field spectra.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A Class of Sparse Johnson–Lindenstrauss Transforms and Analysis of their Extreme Singular Values

The Johnson–Lindenstrauss (JL) lemma is a powerful tool for dimensionality reduction in modern algorithm design. The lemma states that any set of high-dimensional points in a Euclidean space can be projected into lower dimensions while approximately preserving pairwise Euclidean distances. Random matrices satisfying this lemma are called JL transforms (JLTs). Inspired by existing $s$-hashing JLTs with exactly $s$ nonzero elements on each column, the present work introduces an ensemble of sparse matrices encompassing so-called $s$-hashing-like matrices whose expected number of nonzero elements on each column is $s$. The independence of the sub-Gaussian entries of these matrices and the knowledge of their exact distribution play an important role in their analyses. Using properties of independent sub-Gaussian random variables, these matrices are demonstrated to be JLTs, and their smallest nontrivial singular values and largest singular values are estimated nonasymptotically using a technique from geometric functional analysis. As the dimensions of the matrix grow to infinity, these singular values are proved to converge almost surely to fixed quantities (by using the universal Bai–Yin law) and in distribution to the Gaussian orthogonal ensemble Tracy–Widom law after proper rescalings. Understanding the behaviors of extreme singular values is important in general because they are often used to define a measure of stability of matrix algorithms. For example, JLTs were recently used in derivative-free optimization algorithmic frameworks to select random subspaces in which are constructed random models or poll directions to achieve scalability, and hence estimating their smallest singular value in particular helps determine the dimension of these subspaces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Impact of reservoir parameters and wellbore permeability uncertainties on CO 2 and brine leakage potential at the Shenhua CO 2 Storage Site, China

Here, this paper describes a system-level risk assessment for the Shenhua CO 2 storage site, China, using the National Risk Assessment Partnership Integrated Assessment Model for Carbon Storage (NRAP-IAM-CS). Here, we begin by determining the optimal number of Monte Carlo (MC) simulations to achieve CO 2 and brine leakage result convergence. Then, we calculate mass CO 2 and brine leakage to the atmosphere and a hypothetical shallow aquifer. Finally, we assess the geochemical impacts in the event of leakage as if there were a shallow freshwater aquifer at the Shenhua site. Simulation results show that leakage results tend to stabilize after 300 Monte Carlo simulations. When the three wells on site are assigned a permeability of 10 -11 m 2 (representing significantly leaking wells), moderate CO 2 and brine leakages occur, and the percentage of CO 2 leakage exceeds the threshold value we set based on the Intergovernmental Panel on Climate Change (IPCC). This is, however, unlikely to be the case for the Shenhua site. For all the other scenarios, the CO 2 leakage is trivial although there is still the possibility of CO 2 leakage into the groundwater aquifer exceeding the 1% threshold over 1000 years, assuming constant legacy wellbore permeability. For the significantly leaking scenario, there is a 10% probability to have a moderate (2.7 × 10 7 m 3 ) leakage-affected volume in the shallow aquifer with the pH below 6.5, and a large (1.4 × 10 8 m 3 ) volume with the total dissolved solids (TDS) above 500 ppm, hence pH and TDS may be considered for site monitoring plans. Based on the simulation results, there is a very low probability of significant CO 2 /brine leakage through the existing wells at the Shenhua CO 2 storage site.

58 GEOSCIENCES↗
DOE OSTI2024

Advancing our understanding of instabilities in high-energy-density systems through the study of benchmark datasets with advanced modeling tools

Numerous types of pulsed power driven inertial confinement fusion (ICF) and high energy density (HED) systems rely on implosion stability to achieve desired temperatures, pressures, and densities. Sandia National Laboratories Pulsed Power Sciences Center’s main ICF platform, Magnetized Liner Inertial Fusion (MagLIF), suffers from implosion instabilities which limit attainable fuel conditions and can compromise fuel confinement. This Truman Fellowship research primarily focused on computationally exploring (a) methods for improving our understanding of hydrodynamic and magnetohydrodynamic instabilities that form during cylindrical liner implosions, (b) methods for mitigating implosion instabilities, particularly those that degrade performance of MagLIF targets, and (c) novel MagLIF target designs intended to improve target performance primarily via enhanced implosion stability. Several multi-dimensional computational tools were used, including the magnetohydrodynamics code ALEGRA, the radiation-magnetohydrodynamics code HYDRA, and the magnetohydrodynamics code KRAKEN. This research succeeded in executing and analyzing simulations of automagnetizing liner implosions, shockless MagLIF implosions, dynamic screw pinch driven cylindrical liner implosions, and cylindrically convergent HED instability studies. The methods and tools explored and developed in this Truman Fellowship research have been published in several peer-reviewed journal articles and will serve as useful contributions to the fields of pulsed power science and engineering, particularly pertaining to pulsed power ICF and HED science.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Scale-dependent enhancement of productivity and stability in xenic Nannochloropsis cultures

Although bacterial amendments can enhance algal productivity and stability, the development of algal-bacteria consortia for commercial-scale utilization is limited. Here, for this work, we used an established high-throughput approach to generate algal-bacteria consortia, and tested consortia performance at spatial scales from microplates to 320 L raceway ponds. We used both lab and field-reared strains of Nannochloropsis oceanica to build consortia. In some experiments, we imposed environmental perturbations to test the ability of bacteria to enhance algal culture stability. In repeated assays at the scale of well plates, flasks, and bioreactors, strong effects of bacterial amendments on N. oceanica were observed. These effects were most dramatic when cultures experienced stressors such as temperature perturbations or removal of CO 2 augmentation. Isolates that were advanced for field testing included species in the genera Algoriphagus, Oceanicaulis, and Marinobacter. When consortia were generated in the field, positive effects of bacterial amendments were not observed. The amended bacteria were outcompeted, and bacterial community composition across treatments converged after the first grow out. These results highlight the complexity of using consortia in open systems, where interactions between the existing bacterial community, inoculated bacteria, and changing environmental conditions are layered upon other differences in scale and cultivation regimes. Moreover, functionally redundant bacteria are likely present in the field. Following this work, we hypothesize that tight interactions (e.g., obligate relationships between partners) will scale more predictably to outdoor systems. We suggest relying on true synthetic ecological approaches in which the relationships between bacteria and algae partners are well understood, or synthetic ecological approaches coupled with high throughput approaches to design and test consortia. We also recommend future work to examine the effect of algae-bacteria inoculation ratios on productivity and stability, track dynamics of partners through time, and manage ponds to retain beneficial symbioses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Quasi-Lindblad pseudomode theory for open quantum systems

Here, we introduce a new framework to study the dynamics of open quantum systems with linearly coupled Gaussian baths. Our approach replaces the continuous bath with an auxiliary discrete set of pseudomodes with dissipative dynamics, but we further relax the complete positivity requirement in the Lindblad master equation and formulate a quasi-Lindblad pseudomode theory. We show that this quasi-Lindblad pseudomode formulation directly leads to a representation of the bath correlation function in terms of a complex weighted sum of complex exponentials, an expansion that is known to be rapidly convergent in practice and thus leads to a compact set of pseudomodes. The pseudomode representation is not unique and can differ by a gauge choice. When the global dynamics can be simulated exactly, the system dynamics is unique and independent of the specific pseudomode representation. However, the gauge choice may affect the stability of the global dynamics, and we provide an analysis of why and when the global dynamics can retain stability despite losing positivity. We showcase the performance of this formulation across various spectral densities in both bosonic and fermionic problems, finding significant improvements over conventional pseudomode formulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Asymptotic solutions and convergence studies of the resistive inner region equations

A recent publication presents a detailed prescription for constructing global eigenfunctions and complex growth rates for resistive instabilities in axisymmetric toroidal systems, using the method of matched asymptotic expansions. At each singular surface in the plasma, the small-x asymptotic solutions of the ideal MHD outer-region equations are matched to the large-x asymptotic solutions of the resistive inner-region equations, with x being the distance from the singular surface. The success of this method depends upon accurate asymptotic solutions of the resistive inner region equations. Extensive studies have shown that there are not only regimes of excellent agreement with straight-through methods, without asymptotic matching, but also regimes of strong disagreement. As such, in an effort to find the cause of this discrepancy, we present here a different method of constructing asymptotic solutions of the inner region equations and we study the convergence properties of the new and old solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Light-independent regulation of algal photoprotection by CO 2 availability

Photosynthetic algae have evolved mechanisms to cope with suboptimal light and CO 2 conditions. When light energy exceeds CO 2 fixation capacity, Chlamydomonas reinhardtii activates photoprotection, mediated by LHCSR1/3 and PSBS, and the CO 2 Concentrating Mechanism (CCM). How light and CO 2 signals converge to regulate these processes remains unclear. Here, we show that excess light activates photoprotection- and CCM-related genes by altering intracellular CO 2 concentrations and that depletion of CO 2 drives these responses, even in total darkness. High CO 2 levels, derived from respiration or impaired photosynthetic fixation, repress LHCSR3/CCM genes while stabilizing the LHCSR1 protein. Finally, we show that the CCM regulator CIA5 also regulates photoprotection, controlling LHCSR3 and PSBS transcript accumulation while inhibiting LHCSR1 protein accumulation. This work has allowed us to dissect the effect of CO 2 and light on CCM and photoprotection, demonstrating that light often indirectly affects these processes by impacting intracellular CO 2 levels.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Meshless discretization of the discrete-ordinates transport equation with integration based on Voronoi cells

The time-dependent, gray, linear radiation transport equation is discretized using the meshless local Petrov-Galerkin method with reproducing kernels. The integration is performed using a Voronoi tessellation, which creates a partition of unity that only depends on the position and extent of the kernels. The resolution of the integration automatically follows the particles and requires no manual adjustment. The discretization includes streamline-upwind Petrov-Galerkin stabilization to prevent oscillations and improve numerical conditioning. The angular quadrature is selectively refineable to increase angular resolution in chosen directions. The time discretization is done using backward Euler. The transport solve for each direction and the solve for the scattering source are both done using Krylov iterative methods. The results indicate first-order convergence in time and second-order convergence in space for linear reproducing kernels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Advances in 3D Gel Printing for Enzyme Immobilization

Incorporating enzymes with three-dimensional (3D) printing is an exciting new field of convergence research that holds infinite potential for creating highly customizable components with diverse and efficient biocatalytic properties. Enzymes, nature’s nanoscale protein-based catalysts, perform crucial functions in biological systems and play increasingly important roles in modern chemical processing methods, cascade reactions, and sensor technologies. Immobilizing enzymes on solid carriers facilitates their recovery and reuse, improves stability and longevity, broadens applicability, and reduces overall processing and chemical conversion costs. Three-dimensional printing offers extraordinary flexibility for creating high-resolution complex structures that enable completely new reactor designs with versatile sub-micron functional features in macroscale objects. Immobilizing enzymes on or in 3D printed structures makes it possible to precisely control their spatial location for the optimal catalytic reaction. Combining the rapid advances in these two technologies is leading to completely new levels of control and precision in fabricating immobilized enzyme catalysts. The goal of this review is to promote further research by providing a critical discussion of 3D printed enzyme immobilization methods encompassing both post-printing immobilization and immobilization by physical entrapment during 3D printing. Especially, 3D printed gel matrix techniques offer mild single-step entrapment mechanisms that produce ideal environments for enzymes with high retention of catalytic function and unparalleled fabrication control. Examples from the literature, comparisons of the benefits and challenges of different combinations of the two technologies, novel approaches employed to enhance printed hydrogel physical properties, and an outlook on future directions are included to provide inspiration and insights for pursuing work in this promising field.

3D printing↗
DOE OSTI2026

Mechanisms of Alkali Ionic Transport in Amorphous Oxyhalides Solid State Conductors

Amorphous oxyhalides have attracted significant attention due to their relatively high ionic conductivity (1 mS cm –1 ), excellent chemical stability, mechanical softness, and facile synthesis routes via standard solid‐state reactions. These materials exhibit an ionic conductivity that is almost independent of the underlying chemistry, in stark contrast to what occurs in crystalline conductors. In this work, we employ machine learning interatomic potentials to construct large‐scale molecular dynamics trajectories encompassing hundreds of nanoseconds to obtain statistically converged transport properties. We find that the amorphous state consists of chain fragments of metal‐anion tetrahedra of various lengths. By analyzing the residence time of alkali cations migrating around tetrahedrally‐coordinated metals, we find that oxygen anions limit alkali diffusion. By computing the full Einstein expression of the ionic conductivity, we demonstrate that the alkali transference number of these materials is strongly influenced by distinct‐particles correlations, while alkali transport is dictated by uncorrelated self‐diffusion. By extending this analysis to chemical compositions AMX 2.5 O 0.75 , spanning different alkaline (A = Li, Na, K), metallic (M = Al, Ga, In), and halogen (X = Cl, Br, I) species, we clarify why the diffusion properties of these materials remain largely insensitive to variations in atomic isovalent chemistry.

amorphous materials↗
DOE OSTI2020

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Sierra/SolidMechanics Verification Report

This document presents the Sierra/SolidMechanics (Sierra/SM) verification plan. This plan centers on the tests in the Sierra/SM verification test suite, a subset of which are documented in the Sierra/SM Verification Tests Manual. Most of these tests are run nightly with the Sierra/SM code suite, and the results of the tests are checked against analytic solutions. For each of the tests presented in the Verification Tests Manual, the test setup, a description of the analytic solution, and comparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergence is also checked on a nightly basis for several of these tests. This verification plan discusses these various types of tests and what they mean for Sierra/SM verification. Many other activities also contribute to Sierra/SM quality. These address code and solution quality and range from low-level unit tests, run nightly, up to full-fidelity acceptance tests, used to verify release stability. This acceptance test suite checks that new versions of Sierra/SM continue to yield the same answers for high-resolution analyst problems. Further code quality measures include an extensive suite of intermediate-size regression tests and automated nightly code quality checks. While these additional activities do not fall under a strict definition of verification, they greatly add to the quality, stability, and reliability of Sierra/SM, and are discussed here as well.

97 MATHEMATICS AND COMPUTING↗