Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence rate”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Construction of a Code Verification Matrix for Heat Conduction With Finite Element Code Applications

When establishing the pedigree of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model. During this process, numerical results on various meshes are systematically compared to a reference analytic solution. The selection of analytic solutions can be a laborious process, as it is difficult to establish adequate code confidence without performing redundant work. In this report we address this issue by applying a physics-based process that establishes a set of reference problems. In this process, code simulation options are categorized and systematically tested, which ensures that gaps in testing are easily identified and addressed. The resulting problems are primarily intended for code verification analysis but may also be useful for comparison to other simulation codes, troubleshooting activities, or training exercises. The process is used to select fifteen code verification problems relevant for the one-dimensional steady-state heat conduction equation. These problems are applicable to a wide variety of simulation tools, but, in this work, a demonstration is performed using the finite element-based nuclear fuel performance code BISON. Convergence to the analytic solution at the theoretical rate is quantified for a selection of the problems, which establishes a baseline pedigree for the code. Not only can this standard set of conduction solutions be used for verification of other codes, but also the physics-based process for selecting problems can be utilized to quantify and expand testing for any simulation tool.

42 ENGINEERING↗
DOE OSTI2024

Algorithm 1049: The Delaunay Density Diagnostic

Accurate approximation of a real-valued function depends on two aspects of the available data: the density of inputs within the domain of interest and the variation of the outputs over that domain. There are few methods for assessing whether the density of inputs is sufficient to identify the relevant variations in outputs—i.e., the “geometric scale” of the function—despite the fact that sampling density is closely tied to the success or failure of an approximation method. In this article, we introduce a general purpose, computational approach to detecting the geometric scale of real-valued functions over a fixed domain using a deterministic interpolation technique from computational geometry. The algorithm is intended to work on scalar data in moderate dimensions (2–10). Our algorithm is based on the observation that a sequence of piecewise linear interpolants will converge to a continuous function at a quadratic rate (in L 2 norm) if and only if the data are sampled densely enough to distinguish the feature from noise (assuming sufficiently regular sampling). We present numerical experiments demonstrating how our method can identify feature scale, estimate uncertainty in feature scale, and assess the sampling density for fixed (i.e., static) datasets of input–output pairs. Finally, we include analytical results in support of our numerical findings and have released lightweight code that can be adapted for use in a variety of data science settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Novel transition dynamics of topological solitons

Continuous phase transitions can be classified into ones characterized by local-order parameters and others that need additional topological constraints. The critical dynamics near the former transitions have been extensively studied, but the latter is less understood. We fill this gap in knowledge by studying the transition dynamics to a parity-breaking topological ground state called the chiral soliton lattice in quantum chromodynamics at finite temperature, baryon chemical potential, and an external magnetic field. We find a slowing down of the soliton’s translational motion as the critical magnetic field approaches, while the local relaxation rate to the solitonic state remains finite. The time required to converge to the stationary state strongly depends on the symmetry of the initial configuration, determining whether translational motion occurs during the dynamic process. Published by the American Physical Society 2024

Nishimura, Kentaro↗
DOE OSTI2025

Pulse-by-pulse transient thermal deformation in crystal optics under high-repetition-rate FEL

Time-domain modeling of the thermal deformation of crystal optics can help define acceptable operational ranges across the pulse-energy repetition-rate phase space. In this paper, we have studied the transient thermal deformation of a water-cooled diamond crystal for a cavity-based X-ray free-electron laser (CBXFEL), either an X-ray free-electron laser oscillator (XFELO) or a regenerative amplifier X-ray free-electron laser (RAFEL), by numerical simulations including finite-element analysis and advanced data processing. Pulse-by-pulse transient thermal deformation of a 50 µm-thick diamond crystal has been performed with X-ray pulse repetition rates between 50 kHz and 1 MHz. Results for temperature and thermal deformation have been compared with the results of transient analysis using a continuous wave (CW) power loading. Temperature and thermal deformation results from pulse-by-pulse transient analysis vary with time about the results for the CW case for the same average power. The variation amplitude increases with pulse energy and decreases with repetition rate. When the repetition rate increases to infinity, both temperature and thermal deformation converge to the results for the CW case. Two critical time scales for the operation of crystal optics in a CBXFEL are (1) first-turn time, i.e. the time for the XFEL pulse to complete the first turn around the cavity so that the crystal sees the recirculated XFEL pulse, and (2) period-end time, i.e. the time that the next electron bunch arrives for the amplification, so that the crystal outcouples the amplified FEL power. For the same average power, simulation results show that the crystal thermal deformation seen by the XFEL beam decreases with repetition rate at the first-turn time of a 300 m-long cavity and increases with repetition rate at the period-end time. For the wavefront preservation requirement of the crystal optics, a pulse-energy versus repetition-rate phase space has been established. The upper bounds of the pulse energy at both first-turn and period-end times decreases with repetition rate, especially at the period-end time. The upper bound of the thermal deformation of the crystal at the period-end time for any repetition frequency can be estimated from the CW case. For a water-cooled diamond crystal of dimension 5 mm × 5 mm × 0.05 mm, the time to reach a quasi steady-state is about 50 ms for temperature and 50 µs for thermal deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Locally adaptive activation functions with slope recovery for deep and physics-informed neural networks

Here we propose two approaches of locally adaptive activation functions namely, layer-wise and neuron-wise locally adaptive activation functions, which improve the performance of deep and physics-informed neural networks. The local adaptation of activation function is achieved by introducing a scalable parameter in each layer (layer-wise) and for every neuron (neuron-wise) separately, and then optimizing it using a variant of stochastic gradient descent algorithm. In order to further increase the training speed, an activation slope-based slope recovery term is added in the loss function, which further accelerates convergence, thereby reducing the training cost. On the theoretical side, we prove that in the proposed method, the gradient descent algorithms are not attracted to sub-optimal critical points or local minima under practical conditions on the initialization and learning rate, and that the gradient dynamics of the proposed method is not achievable by base methods with any (adaptive) learning rates. We further show that the adaptive activation methods accelerate the convergence by implicitly multiplying conditioning matrices to the gradient of the base method without any explicit computation of the conditioning matrix and the matrix–vector product. The different adaptive activation functions are shown to induce different implicit conditioning matrices. Furthermore, the proposed methods with the slope recovery are shown to accelerate the training process.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Microstructure Scale Lithium-Ion Battery Modeling, Part IV: The Representativity of Microstructure Parameters and Electrochemical Response

Lithium-ion battery electrochemical models require an accurate description of the electrodes microstructures to be predictive, that can be achieved through nanoscale imaging. Such observations are however limited by their field of view (FOV), as they provide only a subset of the whole electrode volume that does not necessarily represent the whole electrode microstructure heterogeneity, and therefore can bias the analysis. A representativity analysis has been performed on the microstructure parameters and, in a novel way, on the full cell electrochemical response to evaluate the predictions representativeness, and thus relevance, of a microstructure scale electrochemical model. The microstructure parameter deviation propagations to the electrochemical response have been quantified for different charge rates. This defines a threshold for the microstructure parameters FOV for a desired maximum deviation of the electrochemical response. Electrochemical model shows cell representative section areas are increasing with C-rate, due to higher in-plane heterogeneities, indicating larger FOVs are required specifically for fast charge modeling. Representativity analysis determines a cell FOV of 144.4 × 154.4 μm 2 is large enough to establish a convergence on the representative section areas for low-intermediate C-rate (≤2.5 C), therefore positively concludes on the model representativeness for these rates, but is not large enough to conclude for higher rates.

25 ENERGY STORAGE↗
DOE OSTI2022

Surrogate modeling for efficiently, accurately and conservatively estimating measures of risk

We present a surrogate modeling framework for conservatively estimating measures of risk from limited realizations of an expensive physical experiment or computational simulation. Risk measures combine objective probabilities with the subjective values of a decision maker to quantify anticipated outcomes. Given a set of samples, we construct a surrogate model that produces estimates of risk measures that are always greater than their empirical approximations obtained from the training data. These surrogate models limit over-confidence in reliability and safety assessments and produce estimates of risk measures that converge much faster to the true value than purely sample-based estimates. We first detail the construction of conservative surrogate models that can be tailored to a stakeholder’s risk preferences and then present an approach, based on stochastic orders, for constructing surrogate models that are conservative with respect to families of risk measures. Our surrogate models include biases that permit them to conservatively estimate the target risk measures. We provide theoretical results that show that these biases decay at the same rate as the L 2 error in the surrogate model. Numerical demonstrations confirm that risk-adapted surrogate models do indeed overestimate the target risk measures while converging at the expected rate.

42 ENGINEERING↗
DOE OSTI2020

Lignin lags, leads, or limits the decomposition of litter and soil organic carbon

Lignin’s role in litter and soil organic carbon (SOC) decomposition remains contentious. Lignin decomposition was traditionally thought to increase during midstage litter decomposition, when cellulose occlusion by lignin began to limit mass loss. Alternatively, lignin decomposition could be greatest in fresh litter as a consequence of co-metabolism, and lignin might decompose faster than bulk SOC. To test these competing hypotheses, we incubated 10 forest soils with C 4 grass litter (amended with 13 C-labeled or unlabeled lignin) over 2 yr and measured soil respiration and its isotope composition. Early lignin decomposition was greatest in 5 of 10 soils, consistent with the co-metabolism hypothesis. However, lignin decomposition peaked 6–24 months later in the other five soils, consistent with the substrate-limitation hypothesis; these soils were highly acidic. Rates of lignin, litter, and SOC decomposition tended to converge over time. Cumulative lignin decomposition was never greater than SOC decomposition; lignin decomposition was significantly lower than SOC decomposition in six soils. Net nitrogen mineralization predicted lignin decomposition ratios relative to litter and SOC. Although the onset of lignin decomposition can indeed be rapid, lignin still presents a likely bottleneck in litter and SOC decomposition, meriting a reconsideration of lignin’s role in modern decomposition paradigms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Receptance coupling substructure analysis and chatter frequency-informed machine learning for milling stability

This paper describes a milling stability identification approach that simultaneously considers: physics-based models for the tool tip frequency response functions and stability predictions; the binary result from a milling test (automatically labeled as stable or unstable based on frequency content); chatter frequency when an unstable result is obtained; and user risk tolerance. The algorithm applies probabilistic Bayesian machine learning with adaptive, parallelized Markov Chain Monte Carlo sampling to update the probability of stability with each milling test. Furthermore, the result is a robust solution for rapid convergence to optimized milling parameters for maximum metal removal rate using all available information.

42 ENGINEERING↗
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 OSTI2025

Depth-Resolved Lithiated Gradients in Pristine and Laser-Ablated Anodes During Fast Charging

Laser ablating 3D electrode microstructures is a technique to improve Li-ion battery fast-charge performance. This technique has been theoretically proposed and electrochemically validated previously in the literature. The fundamental principle underlying laser ablation is that the ablated features reduce Li-ion transport pathways, improving access to the electrode active material near the current collector. This, in turn, promotes more homogeneous electrode utilization. The present study seeks to directly affirm the physics attributed to laser ablation using operando high-speed synchrotron X-ray diffraction. In this study, depth-resolved graphite lithiation gradients are measured operando during high-rate (15 min) charging. The depth-resolved lithiation dynamics of both ablated and non-ablated anodes are compared. The results highlight that the laser-ablated graphite electrode has notably more homogeneous utilization as compared to the non-ablated electrode. Additionally, the ablated electrode has a significant delay in reaching the maximum graphite lithiation at the separator, indicating less propensity for lithium plating. During low rate delithiation/discharge (2 hr), the two cells' lithiation gradients converge. Notably, a calibrated physics-based electrochemical model accurately reflects experimental findings, suggesting the potential to use pseudo-4D models not only to optimize laser ablation parameters in fast-charge capable electrodes but also to guide fast charging protocols that avoid lithium plating.

25 ENERGY STORAGE↗
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.

42 ENGINEERING↗
DOE OSTI2021

Surrogate Modeling For Efficiently Accurately and Conservatively Estimating Measures of Risk

We present a surrogate modeling framework for conservatively estimating measures of risk from limited realizations of an expensive physical experiment or computational simulation. We adopt a probabilistic description of risk that assigns probabilities to consequences associated with an event and use risk measures, which combine objective evidence with the subjective values of decision makers, to quantify anticipated outcomes. Given a set of samples, we construct a surrogate model that produces estimates of risk measures that are always greater than their empirical estimates obtained from the training data. These surrogate models not only limit over-confidence in reliability and safety assessments, but produce estimates of risk measures that converge much faster to the true value than purely sample-based estimates. We first detail the construction of conservative surrogate models that can be tailored to the specific risk preferences of the stakeholder and then present an approach, based upon stochastic orders, for constructing surrogate models that are conservative with respect to families of risk measures. The surrogate models introduce a bias that allows them to conservatively estimate the target risk measures. We provide theoretical results that show that this bias decays at the same rate as the L 2 error in the surrogate model. Our numerical examples confirm that risk-aware surrogate models do indeed over-estimate the target risk measures while converging at the expected rate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Multi-variance replica exchange SGMCMC for inverse and forward problems via Bayesian PINN

Physics-informed neural network (PINN) has been successfully applied in solving a variety of nonlinear non-convex forward and inverse problems. However, the training is challenging because of the non-convex loss functions and the multiple optima in the Bayesian inverse problem. In this work, we propose a multi-variance replica exchange stochastic gradient Langevin dynamics method to tackle the challenge of the multiple local optima in the optimization and the challenge of the multiple modal posterior distribution in the inverse problem. Replica exchange methods are capable of escaping from the local traps and accelerating the convergence; two chains with different temperatures are designed where the low temperature chain aims for the local convergence, and the target of the high temperature chain is to travel globally and explore the whole loss function entropy landscape. However, it may not be efficient to solve mathematical inversion problems by using the vanilla replica method directly since the method doubles the computational cost in evaluating the forward solvers (likelihood functions) in the two chains. To address this issue, we propose to make different assumptions on the energy function estimation and this facilities one to use solvers of different fidelities in the likelihood function evaluation. More precisely, one can use a solver with low fidelity in the high temperature chain while using a solver with high fidelity in the low temperature chain. Our proposed method significantly lowers the computational cost in the high temperature chain, meanwhile preserving the accuracy and converging very fast. Here we give an unbiased estimate of the swapping rate and give an estimation of the discretization error of the scheme. To verify our idea, we design and solve four inverse problems which have multiple modes. The proposed method is also employed to train the Bayesian PINN to solve the forward and inverse problems; faster and more accurate convergence has been observed when compared to the stochastic gradient Langevin dynamics (SGLD) method and vanilla replica exchange methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Reliability Improvement and Effective Switching Layer Model of Thin-Film MoS 2 Memristors

2D memristors have demonstrated attractive resistive switching characteristics recently but also suffer from the reliability issue, which limits practical applications. Previous efforts on 2D memristors have primarily focused on exploring new material systems, while damage from the metallization step remains a practical concern for the reliability of 2D memristors. Here, the impact of metallization conditions and the thickness of MoS 2 films on the reliability and other device metrics of MoS 2 -based memristors is carefully studied. The statistical electrical measurements show that the reliability can be improved to 92% for yield and improved by ≈16× for average DC cycling endurance in the devices by reducing the top electrode (TE) deposition rate and increasing the thickness of MoS 2 films. Intriguing convergence of switching voltages and resistance ratio is revealed by the statistical analysis of experimental switching cycles. An “effective switching layer” model compatible with both monolayer and few-layer MoS 2 , is proposed to understand the reliability improvement related to the optimization of fabrication configuration and the convergence of switching metrics. In conclusion, the Monte Carlo simulations help illustrate the underlying physics of endurance failure associated with cluster formation and provide additional insight into endurance improvement with device fabrication optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

End-gas autoignition fraction and flame propagation rate in laser-ignited primary reference fuel mixtures at elevated temperature and pressure

Knock in spark-ignited (SI) engines is initiated by autoignition of the unburned gasses upstream of spark-ignited, propagating, turbulent premixed flames. Knock propensity of fuel/air mixtures is typically quantified using research octane number (RON), motor octane number (MON), or methane number (MN; for gaseous fuels), which are measured using single-cylinder, variable compression ratio engines. In this study, knock propensity of SI fuels was quantified via observations of end-gas autoignition (EGAI) in unburned gasses upstream of laser-ignited, premixed flames at elevated pressures and temperatures in a rapid compression machine. Stoichiometric primary reference fuel (PRF; n-heptane/isooctane) blends of varying reactivity (50 ≤ PRF ≤ 100) were ignited using an Nd:YAG laser over a range of temperatures and pressures, all in excess of 545 K and 16.1 bar. Laser ignition produced outwardly-propagating premixed flames. High-speed pressure measurements and schlieren images indicated the presence of EGAI. The fraction of the total heat release attributed to EGAI (i.e., EGAI fraction) varied with fuel reactivity (i.e., octane number) and the time-integrated temperature of the end-gas prior to ignition. Flame propagation rates, which were measured using schlieren images, were only weakly correlated with octane number but were affected by turbulence caused by variation in piston timing. Under conditions of low turbulence, measured flame propagation rates approached one-dimensional premixed laminar flame speed computations performed at the same conditions. Experiments were simulated with a three-dimensional CONVERGE™ model using reduced chemical kinetics (121 species, 538 reactions). The simulations accurately captured the measured flame propagation rates, as well as the variation in EGAI fraction with fuel reactivity and time-integrated end-gas temperature. The simulations also revealed low-temperature heat release as well as formaldehyde and hydrogen peroxide formation in the end-gas upstream of the propagating flame, which increased the temperature and degree of chain branching in the end-gas, ultimately leading to EGAI.

42 ENGINEERING↗
DOE OSTI2026

Convergent evolution of aerobic fermentation through divergent mechanisms acting on the same key glycolytic genes

Expression profiling by high throughput sequencing The model yeast Saccharomyces cerevisiae has evolved aerobic fermentation in part through a whole genome duplication and subsequent retention of glycolytic genes. To identify independent yeast lineages with a high glycolytic rate phenotype we developed an assay measuring extracellular acidification rates (ECAR) to evaluate the glycolytic rates across diverse yeast species. This assay identified a novel group of yeasts in the genus Saturnispora with rapid glycolytic rates. Through comparative approaches, we found that several glycolytic genes encoding hexokinases and enzymes in lower glycolysis had higher expression and modifications in promoter sequences in rapid ECAR species. Intriguingly, many of the upregulated genes are the same genes that are duplicated in S. cerevisiae. When the transcription factor required for their activation was deleted, the mutants had a slow glycolytic rate and more respiratory phenotype. This work shows how the aerobic fermentation phenotype has convergently evolved in Saturnispora through transcriptional rewiring. This divergent mechanism nevertheless impacts the same glycolytic genes, which suggests that there are evolutionary constraints on how aerobic fermentation can arise.

Horianopoulos, Linda C. [Department of Food Scienc↗
DOE OSTI2023

Delayed leaf greening involves a major shift in the expression of cytosolic and mitochondrial ribosomes to plastid ribosomes in the highly phosphorus-use-efficient Hakea prostrata (Proteaceae)

Abstract Background and aims Hakea prostrata (Proteaceae) is a highly phosphorus-use-efficient plant native to southwest Australia. It maintains a high photosynthetic rate at low leaf phosphorus (P) and exhibits delayed leaf greening, a convergent adaptation that increases nutrient-use efficiency. This study aimed to provide broad physiological and gene expression profiles across leaf development, uncovering pathways leading from young leaves as nutrient sinks to mature leaves as low-nutrient, energy-transducing sources. Methods To explore gene expression underlying delayed greening, we analysed a de novo transcriptome for H. prostrata across five stages of leaf development. Photosynthesis and respiration rates, and foliar pigment, P and nitrogen (N) concentrations were determined, including the division of P into five biochemical fractions. Key results Transcripts encoding functions associated with leaf structure generally decreased in abundance across leaf development, concomitant with decreases in foliar concentrations of 85% for anthocyanins, 90% for P and 70% for N. The expression of genes associated with photosynthetic function increased during or after leaf expansion, in parallel with increases in photosynthetic pigments and activity, much later in leaf development than in species that do not have delayed greening. As leaves developed, transcript abundance for cytosolic and mitochondrial ribosomal proteins generally declined, whilst transcripts for chloroplast ribosomal proteins increased. Conclusions There was a much longer temporal separation of leaf cell growth from chloroplast development in H. prostrata than is found in species that lack delayed greening. Transcriptome-guided analysis of leaf development in H. prostrata provided insight into delayed greening as a nutrient-saving strategy in severely phosphorus-impoverished landscapes.

59 BASIC BIOLOGICAL SCIENCES↗