SEARCH · Search NASA

Results for “Iterative methods”

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 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

DOE OSTI2021

A persistent adjoint method with dynamic time-scaling and an application to mass action kinetics

In this article, we consider an optimization problem where the objective function is evaluated at the fixed-point of a contraction mapping parameterized by a control variable, and optimization takes place over this control variable. Since the derivative of the fixed-point with respect to the parameter can usually not be evaluated exactly, an adjoint dynamical system can be used to estimate gradients. Using this estimation procedure, the optimization algorithm alternates between derivative estimation and an approximate gradient descent step. We analyze a variant of this approach involving dynamic time-scaling, where after each parameter update the adjoint system is iterated until a convergence threshold is passed. Here, we prove that, under certain conditions, the algorithm can find approximate stationary points of the objective function. We demonstrate the approach in the settings of an inverse problem in chemical kinetics, and learning in attractor networks.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Direct Estimation of Parameters in ODE Models Using WENDy: Weak-Form Estimation of Nonlinear Dynamics

Abstract We introduce the Weak-form Estimation of Nonlinear Dynamics (WENDy) method for estimating model parameters for non-linear systems of ODEs. Without relying on any numerical differential equation solvers, WENDy computes accurate estimates and is robust to large (biologically relevant) levels of measurement noise. For low dimensional systems with modest amounts of data, WENDy is competitive with conventional forward solver-based nonlinear least squares methods in terms of speed and accuracy. For both higher dimensional systems and stiff systems, WENDy is typically both faster (often by orders of magnitude) and more accurate than forward solver-based approaches. The core mathematical idea involves an efficient conversion of the strong form representation of a model to its weak form, and then solving a regression problem to perform parameter inference. The core statistical idea rests on the Errors-In-Variables framework, which necessitates the use of the iteratively reweighted least squares algorithm. Further improvements are obtained by using orthonormal test functions, created from a set of $$C^{\infty }$$ C ∞ bump functions of varying support sizes.We demonstrate the high robustness and computational efficiency by applying WENDy to estimate parameters in some common models from population biology, neuroscience, and biochemistry, including logistic growth, Lotka-Volterra, FitzHugh-Nagumo, Hindmarsh-Rose, and a Protein Transduction Benchmark model. Software and code for reproducing the examples is available at https://github.com/MathBioCU/WENDy .

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Experimental confirmation of first-principles thermal conductivity in Zirconium-doped ThO 2

The degradation of thermal conductivity in advanced nuclear fuels due to the accumulation of fission products and irradiation-induced defects is inevitable, and must be considered as part of safety and efficiency analyses of nuclear reactors. Here, this study examines the thermal conductivity of a zirconium-doped ThO 2 crystal, synthesized via the hydrothermal method using a spatial domain thermoreflectance technique. Zirconium is one of the soluble fission products in oxide fuels that can effectively scatter heat-carrying phonons in the crystalline lattice of fuel. Thus, thermal property measurements of zirconium-doped ThO 2 single crystals provide insights into the effects of substitutional zirconium doping, isolated from extrinsic factors such as grain boundary scattering. The experimental results are compared with first-principles calculations of the lattice thermal conductivity of ThO 2 , employing an iterative solution of the Peierls-Boltzmann transport equation. Additionally, the non-perturbative Green's function methodology is utilized to compute phonon-point defect scattering rates, accounting for local distortions around point defects, including mass difference changes, interatomic force constants, and structural relaxation. The congruence between the predicted results from first-principles calculations and the measured temperature-dependent thermal conductivity validates the computational methodology. Furthermore, the methodologies employed in this study enable systematic investigations of thermal conductivity reduction by fission products, potentially leading to the development of more accurate fuel performance codes.

36 - MATERIALS SCIENCE↗

DOE OSTI2020

Global adjoint tomography—model GLAD-M25

SUMMARY Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 ‘Titan’ and the IBM Power 9 ‘Summit’ at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.

58 GEOSCIENCES↗

DOE OSTI2025

Learning with Adaptive Conservativeness for Distributionally Robust Optimization: Incentive Design for Voltage Regulation

Information asymmetry between the Distribution System Operator (DSO) and Distributed Energy Resource Aggregators (DERAs) obstructs designing effective incentives for voltage regulation. To capture this effect, we employ a Stackelberg game-theoretic framework, where the DSO seeks to overcome the information asymmetry and refine its incentive strategies by learning from DERA behavior over multiple iterations. We introduce a model-based online learning algorithm for the DSO, aimed at inferring the relationship between incentives and DERA responses. Given the uncertain nature of these responses, we also propose a distributionally robust incentive design model to control the probability of voltage regulation failure and then reformulate it into a convex problem. This model allows the DSO to periodically revise distribution assumptions on uncertain parameters in the decision model of the DERA. Finally, we present a gradient-based method that permits the DSO to adaptively modify its conservativeness level, measured by the size of a Wasserstein metric-based ambiguity set, according to historical voltage regulation performance. The effectiveness of our proposed method is demonstrated through numerical experiments.

adaptation models↗

DOE OSTI2022

Application-Driven Creation of Building Metadata Models with Semantic Sufficiency

Semantic metadata models such as Brick, RealEstateCore, Project Haystack, and BOT promise to simplify and lower the cost of developing software for smart buildings, enabling the widespread deployment of energy efficiency applications. However, creating these models remains a challenge. Despite recent advances in creating models from existing digital representations like point labels and architectural models, there is still no feedback mechanism to ensure that the human input to these methods results in a model that can actually support the desired software. In this paper, we introduce the notion of semantic sufficiency, a practical principle for semantic metadata model creation that asserts that a model is "finished" when it contains the metadata necessary to support a given set of applications. To support semantic sufficiency, we design a standard representation for capturing application metadata requirements and a templating system for generating common metadata model components with limited user input. We then construct an iterative model creation workflow that integrates metadata requirements to direct the model creation effort, and present several novel optimizations that increase the model utility while minimizing the effort by a human operator. These new abstractions for model creation and validation lower model development costs and ensure the utility of the resulting model, thus facilitating the adoption of intelligent building applications.

applications↗

DOE OSTI2024

Learning with Adaptive Conservativeness for Distributionally Robust Optimization: Incentive Design for Voltage Regulation: Preprint

distribution system operator↗

DOE OSTI2020

A Novel Low-Cost Method of Manufacturing Nb3Sn Superconductors with Multiple-Tin-Tube Sources (CRADA FRA-2008-0001 Final Report)

In order to successfully sustain a fusion reaction, peak magnetic fields on the order of 12-13 Tesla will be required. Magnetic fields of this magnitude can only be accomplished by advanced superconductors such as Nb3Sn. However, the economic success of a fusion machine will depend on further improvements in the cost-performance of the Nb3Sn conductor. This project will develop a novel, low-cost, multiple-tin-tube process as a new manufacturing approach for large-scale Nb3Sn-conductor production. The process will be suitable for the efficient production of larger strands (0.83 mm), consistent with ITER (International Thermonuclear Experimental Reactor) specifications. In Phase I, Cu/Nb composites were manufactured and shaped into sub-elements. These sub-elements were tin-coated by electroplating, assembled into a precise-fit restack billet tube, and drawn to 0.83 mm diameter. In Phase II, the process will be scaled-up to full production levels. The improved Nb3Sn conductor should have an immediate benefit for high-field magnet applications. A prototype fusion machine, based on a cost effective Nb3Sn conductor, would have enormous economic and social benefits. In addition, the conductor should be applicable to nuclear magnetic resonance (NMR), which has requirements for use in chemical research, biochemistry, pharmaceutical chemistry, polymer science, petroleum research, agricultural chemistry, and medicine. Leszek Motowidlo, Principal Investigator, will be responsible for the overall coordination of the effort. He and others of SupraMagnetics staff will design and fabricate protoype PIT Nb3Sn conductors for evaluation at Fermilab. Emanuela Barzi will be responsible for the Fermilab subcontract and will coordinate and supervise critical current testing Nb3Sn strands and cable fabrication and evaluation.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Fuzzified PaCcET for Economic-Emission Scheduling of Microgrids

In this paper, a new approach is proposed to solve a multi-objective economic-emission scheduling problem in microgrids (MGs) by simultaneously minimizing the energy and emission costs of the MG with various distributed energy resources (DERs). The proposed approach is an extension of a computationally effective multiobjective optimization technique, Pareto concavity elimination transformation (PaCcET). The proposed approach, referred to as Fuzzified-PaCcET, employs a fuzzy logic controller to dynamically revise crossover and mutation rates in the original PaCcET leading to the faster convergence of the solution. The proposed approach finds the best Pareto front, also referred to as a Non-dominated set (NDS) of solutions, instead of finding a single optimal solution. In order to find the solutions on concave areas of the Pareto front, an iterative objective space transformation is performed in the PaCcET algorithm to allow a linear combination of objective functions (in the transformed objective space). The proposed Fuzzified-PaCcET-based scheduling is implemented on a MG with various dispatchable and non-dispatchable DERs to find the set of optimal solutions according to the total fuel cost of DERs, as well as the most optimum environmental cost. In order to extract the best compromise solution (BCS) among NDS of solutions, a fuzzy-based method is implemented. The comparison of the simulation results of the Fuzzified-PaCcET with that of PaCcET shows that Fuzzified-PaCcET can generate better solution with less computational burden.

Gautam, Mukesh↗

DOE OSTI2023

Developing a SARS-CoV-2 main protease binding prediction random forest model for drug repurposing for COVID-19 treatment

The coronavirus disease 2019 (COVID-19) global pandemic resulted in millions of people becoming infected with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus and close to seven million deaths worldwide. It is essential to further explore and design effective COVID-19 treatment drugs that target the main protease of SARS-CoV-2, a major target for COVID-19 drugs. In this study, machine learning was applied for predicting the SARS-CoV-2 main protease binding of Food and Drug Administration (FDA)-approved drugs to assist in the identification of potential repurposing candidates for COVID-19 treatment. Ligands bound to the SARS-CoV-2 main protease in the Protein Data Bank and compounds experimentally tested in SARS-CoV-2 main protease binding assays in the literature were curated. These chemicals were divided into training (516 chemicals) and testing (360 chemicals) data sets. To identify SARS-CoV-2 main protease binders as potential candidates for repurposing to treat COVID-19, 1188 FDA-approved drugs from the Liver Toxicity Knowledge Base were obtained. A random forest algorithm was used for constructing predictive models based on molecular descriptors calculated using Mold2 software. Model performance was evaluated using 100 iterations of fivefold cross-validations which resulted in 78.8% balanced accuracy. The random forest model that was constructed from the whole training dataset was used to predict SARS-CoV-2 main protease binding on the testing set and the FDA-approved drugs. Model applicability domain and prediction confidence on drugs predicted as the main protease binders discovered 10 FDA-approved drugs as potential candidates for repurposing to treat COVID-19. Our results demonstrate that machine learning is an efficient method for drug repurposing and, thus, may accelerate drug development targeting SARS-CoV-2.

Research & Experimental Medicine↗

DOE OSTI2021

Fuzzified PaCcET for Economic-Emission Scheduling of Microgrids

42 ENGINEERING↗

DOE OSTI2022

Joint iterative reconstruction and 3D rigid alignment for X-ray tomography

X-ray tomography is widely used for three-dimensional structure determination in many areas of science, from the millimeter to the nanometer scale. The resolution and quality of the 3D reconstruction is limited by the availability of alignment parameters that correct for the mechanical shifts of the sample or sample stage for the images that constitute a scan. In this paper we describe an algorithm for marker-free, fully automated and accurately aligned and reconstructed X-ray tomography data. Our approach solves the tomographic reconstruction jointly with projection data alignment based on a rigid-body deformation model. We demonstrate the robustness of our method on both synthetic phantom and experimental data and show that our method is highly efficient in recovering relatively large alignment errors without prior knowledge of a low resolution approximation of the 3D structure or a reasonable estimate of alignment parameters.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Energetic Particles - microturbulence interaction thrust (EP SIWG white paper)

This white paper outlines a strategic approach to Energetic Particle (EP) research, aimed at applying the results to the future burning devices such as ITER. Our goal is to deepen the understanding of the interplay between plasma background microturbulence and diverse species of energetic ions in burning plasma devices. This strategy is poised to advance our knowledge in this crucial area of fusion science, paving the way for significant breakthroughs in plasma physics and fusion energy research. We consider the confinement of energetic ions such as auxiliary heating beams, minority ions from ICRH, and fusion product alpha particles. Recent studies have underscored the potential of EP/microturbulence interaction to enhance fusion plasma performance by stabilizing microturbulence. However, it also poses risks, such as exacerbating Alfvénic eigenmode (AE) instabilities, which could compromise the sustainability of plasma discharge by causing EP losses. The interplay between EPs and background microturbulence is evident in the phenomenon of effective pitch angle scattering, a crucial aspect of the quasilinear (QL) theory. This theory is integral to developing numerically-efficient yet comprehensive and self-consistent approaches, recently employed to investigate the relaxation of energetic particle populations in the holistic modeling of fusion-grade plasmas. The primary objectives of this white paper thrust encompass a dual focus. First is that we need to delve into the mechanisms through which microturbulence engenders effective pitch angle scattering. This investigation will entail the examination of the presence of known or self-consistently predicted spectra of modes accountable for microturbulence. In these explorations, electrostatic microturbulence serves as a logical initial stride towards achieving the trust's objectives. This endeavor is anticipated to yield formulations expressing the parametric dependencies of the effective pitch angle scattering frequency on variables such as thermal plasma electron and/or ion temperatures, as well as their respective thermal conductivities. In our second objective, we aim to delve into the intricate formation of zonal flow (ZF) structures amidst the complex interplay of microturbulence and Alfvenic eigenmodes (AEs). This endeavor poses greater challenges as we seek to unravel the macroscopic manifestations influenced by microturbulence, termed as zonal structures (ZS), stemming from microturbulence-induced ZF. Drawing from this understanding, we anticipate employing a QL approach to yield comprehensive insights into the distribution function of energetic particles (EP) within phase space. This method entails resolving the dominant multidimensional phase space diffusion processes while effectively averaging over the rapid ballistic responses. Nevertheless, substantial strides remain imperative to realize a comprehensive whole-device modeling framework. This entails meticulous verification and validation exercises against experimental observations, as well as rigorous benchmarking against theoretical frameworks and numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Quarter 4 Report: Report on Final Findings and Opportunities for Future Work in the Use of Mixed Precision in Iterative Solvers

The fourth quarter of the project was spent developing an error analysis of the s-step Lanczos and CG algorithms. Our theoretical bounds and numerical experiments show that the numerical behavior of the algorithm can be significantly improved by using extra precision in a small part of the computation related to the computation and application of the Gram matrix. We have published a technical report which includes all steps of the analysis [8]; a shortened version for journal submission is in preparation. We plan to submit this paper in the following weeks. Activities related to this also include a collaboration with Ichitaro Yamazaki on gathering performance results for these new mixed precision s-step Krylov subspace methods using single/double precision on GPUs. Namely, we would like to obtain performance results that show that the performance overhead of using double the working precision in these select computations is minimal. Other activities include attending biweekly xSDK meetings and presenting a pitch talk on this work to the group on February 25, 2021. In the remainder of the document, we summarize our findings on the potential for mixed precision in classical Krylov subspace methods and s-step Krylov subspace methods, as well as key opportunities for future work.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Feature Engineering and Ensemble Methods for Imbalanced ICS Intrusion Detection: Pipeline Audit and Constrained Evaluation

Industries are becoming increasingly connected and are more vulnerable to cyberattacks due to the widened attack surface. Industrial Control Systems (ICS) are among the most critical sectors that malicious actors can target, as such attacks can cause significant operational disruption and physical damage. It is imperative to detect such attacks as early as possible. This paper evaluates constraint-conditioned optimistic performance estimates for traditional ML models in ICS intrusion detection (i.e., estimates obtained under contiguous, non-shuffled temporal evaluation without test-set alteration, but with pre-split feature engineering that may introduce temporal leakage, due to dataset constraints). Our findings are threefold. First, we quantify how iterative feature engineering affects tree-based ensemble performance and examine how pipeline decisions (split strategy, sampling scope, and cleaning policy) can inflate or reduce reported IDS results under constraint-bound evaluation. Second, we compare intrinsic class-imbalance handling across ensemble models. Third, under our current pipeline constraints (including pre-split feature engineering), CatBoost achieves the best performance on Water Storage Tank (accuracy: 0.9831, class-1 F1: 0.9682), while Light- GBM achieves the best performance on Gas Pipeline (accuracy: 0.9618, class-1 F1: 0.9086).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Use of a segmented cell for the combinatorial development of platinum group metal-free electrodes for polymer electrolyte fuel cells

The development of novel platinum group metal (PGM)-free catalysts is a difficult task, particularly when coupled with the integration into a cathode catalyst layer (CCL). The optimization of such PGM-free electrode structures is often non-linear and iterative, making it a demanding task due to the high number of parameters that affect performance. To accelerate both materials discovery and electrode development, this work demonstrates the application of a high-resolution segmented cell (SC) with 121 segments of 0.413 cm 2 size for combinatorial high-throughput PGM-free catalyst screening and CCL optimization. The approach utilizes three flow-field strips with active areas of 4.45 cm 2 each, distributed over 11 segments. Electrodes with identical catalyst material and fabrication method result in reproducible data with typically less than 10% variation between each segment and match in performance with differential single cell data. High throughput testing of combinatorial sample sets, with varied PGM-free electrocatalyst materials or electrode composition, demonstrated the ability to rapidly discern high performing outliers, and establish or confirm trends in electrode optimization. The results indicated significant performance benefits of electrodes with ionomer content of 45 wt% over lower values, and of water rich inks with 82% water content over those with 50%.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2023

DOE Bioenergy Technologies Office (BETO) 2023 Project Peer Review: WBS 2.4.1.100 Bench Scale Research & Development

Bench Scale Integration develops and optimizes fermentation processes to produce bio-based fuels and chemicals for commercial scale-up. The project uses fermentation science to achieve high titers and production rates by, for example, manipulating how the microorganisms are fed biomass sugars and nutrients, modifying fermentation conditions (pH, temperature, aeration) or developing online control strategies for better fermentation operations and high titer, rates, and yield (TRY). For this period of performance, we continued our development of a commercial-ready 2,3-butanediol (BDO) fermentation from biomass sugars utilizing NREL's proprietary Zymomonas mobilis microorganism. The engineered Z. mobilis can use all the main sugars in corn stover biomass, which are glucose, xylose, and arabinose. BDO is a versatile, low-carbon chemical which can be catalytically upgraded to a variety of hydrocarbon fuels and chemicals. The project had three goals during this review period; evaluate the technical feasibility of using whole slurry pretreated corn stover to achieve the techno-economic analysis (TEA) performance goals, continue optimizing a liquor-based fed-batch fermentation process for high titer, and develop strategies to enable scale-up. After evaluating different iterations of a whole slurry fermentation that did not meet the TEA goals, a Go/No-Go decision was made to pivot to liquor-only with new TEA performance targets, the main one being 140 g/L titer. We successfully met this goal, producing 141 g/L BDO at 1 g/L-hr productivity and 84% process yield. This titer and productivity attracted industrial interest to scale the fermentation resulting in a Technology Commercialization Fund project award in 2022. We used an NREL developed near-infra red (NIR) spectroscopy method for rapid analysis which allowed for changes to aeration levels and sugar feeding during the fermentation to maximize BDO production. The NIR analysis can be done using a hand-held spectrometer, essentially taking the analysis on to the plant floor during scale-up and preliminary work shows the feasibility of using an online probe for continuous monitoring and control. The other scale-up tool is mapping oxygen transfer coefficient (kLa) and oxygen transfer rate (OTR) in various vessels to find conditions that match the optimized 500 mL vessels. Showing a correlation to the mapping work, which is done with water and a dissolved oxygen probe, can reduce the risk of failed fermentations during scale-up. Our end-of-project goal is to meet the design target BDO titer (140-150 g/L) from DMR corn stover liquor at 1000L or larger scale to demonstrate BDO process design case and transfer the technology to industrial fermentation stakeholders for commercialization.

2 3-butanediol↗

DOE OSTI2022

Large-Scale Optimization with Linear Equality Constraints Using Reduced Compact Representation

For optimization problems with linear equality constraints, we prove that the (1,1) block of the inverse KKT matrix remains unchanged when projected onto the nullspace of the constraint matrix. In this work, we develop reduced compact representations of the limited-memory inverse BFGS Hessian to compute search directions efficiently when the constraint Jacobian is sparse. Orthogonal projections are implemented by a sparse QR factorization or a preconditioned LSQR iteration. In numerical experiments two proposed trust-region algorithms improve in computation times, often significantly, compared to previous implementations of related algorithms and compared to IPOPT.

97 MATHEMATICS AND COMPUTING↗