Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence analysis”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Metrics for Intercomparison of Remapping Algorithms (MIRA) protocol applied to Earth system models

Abstract. Strongly coupled nonlinear phenomena such as those described by Earth system models (ESMs) are composed of multiple component models with independent mesh topologies and scalable numerical solvers. A common operation in ESMs is to remap or interpolate component solution fields defined on their computational mesh to another mesh with a different combinatorial structure and decomposition, e.g., from the atmosphere to the ocean, during the temporal integration of the coupled system. Several remapping schemes are currently in use or available for ESMs. However, a unified approach to compare the properties of these different schemes has not been attempted previously. We present a rigorous methodology for the evaluation and intercomparison of remapping methods through an independently implemented suite of metrics that measure the ability of a method to adhere to constraints such as grid independence, monotonicity, global conservation, and local extrema or feature preservation. A comprehensive set of numerical evaluations is conducted based on a progression of scalar fields from idealized and smooth to more general climate data with strong discontinuities and strict bounds. We examine four remapping algorithms with distinct design approaches, namely ESMF Regrid (Hill et al., 2004), TempestRemap (Ullrich and Taylor, 2015), generalized moving least squares (GMLS) (Trask and Kuberry, 2020) with post-processing filters, and WLS-ENOR (Li et al., 2020). By repeated iterative application of the high-order remapping methods to the test fields, we verify the accuracy of each scheme in terms of their observed convergence order for smooth data and determine the bounded error propagation using challenging, realistic field data on both uniform and regionally refined mesh cases. In addition to retaining high-order accuracy under idealized conditions, the methods also demonstrate robust remapping performance when dealing with non-smooth data. There is a failure to maintain monotonicity in the traditional L2-minimization approaches used in ESMF and TempestRemap, in contrast to stable recovery through nonlinear filters used in both meshless GMLS and hybrid mesh-based WLS-ENOR schemes. Local feature preservation analysis indicates that high-order methods perform better than low-order dissipative schemes for all test cases. The behavior of these remappers remains consistent when applied on regionally refined meshes, indicating mesh-invariant implementations. The MIRA intercomparison protocol proposed in this paper and the detailed comparison of the four algorithms demonstrate that the new schemes, namely GMLS and WLS-ENOR, are competitive compared to standard conservative minimization methods requiring computation of mesh intersections. The work presented in this paper provides a foundation that can be extended to include complex field definitions, realistic mesh topologies, and spectral element discretizations, thereby allowing for a more complete analysis of production-ready remapping packages.

58 GEOSCIENCES↗
DOE OSTI2025

Cosmological constraints using Minkowski functionals from the first year data of the Hyper Suprime-Cam

We use Minkowski functionals to analyse weak lensing convergence maps from the first-year data release of the Subaru Hyper Suprime-Cam (HSC-Y1) survey. Minkowski functionals provide a description of the morphological properties of a field, capturing the non-Gaussian features of the Universe matter-density distribution. Using simulated catalogues that reproduce survey conditions and encode cosmological information, we emulate Minkowski functionals predictions across a range of cosmological parameters to derive the best-fit from the data. By applying multiple scales cuts, we rigorously mitigate systematic effects, including baryonic feedback and intrinsic alignments. From the analysis, combining constraints of the angular power spectrum and Minkowski functionals, we obtain S8≡σ8Ωm/0.3=0.808−0.046+0.033 and Ωm=0.293−0.043+0.157⁠. These results represent a 40 per cent improvement on the S8 constraints compared to using power spectrum only. Minkowski functionals results are consistent with other two-point, and higher order statistics constraints using the same data, being in agreement with CMB results from the Planck S8 measurements. Our study demonstrates the power of Minkowski functionals beyond two-point statistics to constrain and break the degeneracy between Ωm and σ8⁠.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A Critical Review of Heliostat Design for Concentrating Solar Thermal Technologies

Heliostat designs have undergone a widespread and eclectic development process, with many unique designs demonstrated. However, recent developments do not show the overall cohort converging toward a globally accepted universal design. Here, this study characterizes heliostats by breaking down and evaluating design traits based on emergent patterns from a comprehensive compilation of known heliostat designs spanning several decades. Four main categories for evaluation emerged: heliostat base, heliostat primary axis, heliostat drive, and facet support. Each of these four categories is further defined by four subtypes so that all heliostats fall into a single subtype within each category. The classified heliostats are ranked, yielding several view slices into the heliostat compilation or a breakdown of heliostats by type. An analysis of the breakdown shows several trends: a scale-up of established designs, new approaches at small to medium scales, and a movement toward greater adoption of linear drives. These trends reflect the most meaningful contributing factors to a proposed trajectory for a new era of heliostat designs striving to meet widely considered cost targets of $\$$50/m 2 or $\$$75/m 2 .

14 SOLAR ENERGY↗
DOE OSTI2026

Analysis of a Heat Balance Experiment in the TREAT Reactor with SPARTA and OpenMC

This article presents the interpretation of a heat balance experiment in TREAT using the SPARTA methodology and an OpenMC model. The study examines the effects of control rods and temperature variations on the effective multiplication factor ($k_{\rm eff}$) and correction factors. Results indicate that both core temperature and rod positions have a linear effect on $k_{\rm eff}$ and correction factors in the upper-west channel. A critical hypothesis of SPARTA was validated: reactivities and correction factors can be obtained by linearly combining individual effects. The previously considered negligible effect of temperature on correction factors was found to be significant, necessitating updates to the SPARTA methodology. By fitting an exponential curve to mean temperature data, SPARTA's accuracy and convergence have been improved. Comparisons between SPARTA and OpenMC showed consistent reactivity removal due to temperature increase. While SPARTA effectively interprets complex transients, simpler methods may suffice when correction factors have minimal influence on reactivity calculations.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Parameter Sensitivity Analysis of the SparTen High Performance Sparse Tensor Decomposition Software (Extended Analysis)

Tensor decomposition models play an increasingly important role in modern data science applications. One problem of particular interest is fitting a low-rank Canonical Polyadic (CP) tensor decomposition model when the tensor has sparse structure and the tensor elements are nonnegative count data. SparTen is a high-performance C++ library which computes a low-rank decomposition using different solvers: a first-order quasi-Newton or a second-order damped Newton method, along with the appropriate choice of runtime parameters. Since default parameters in SparTen are tuned to experimental results in prior published work on a single real-world dataset conducted using MATLAB implementations of these methods, it remains unclear if the parameter defaults in SparTen are appropriate for general tensor data. Furthermore, it is unknown how sensitive algorithm convergence is to changes in the input parameter values. This report addresses these unresolved issues with large-scale experimentation on three benchmark tensor data sets. Experiments were conducted on several different CPU architectures and replicated with many initial states to establish generalized profiles of algorithm convergence behavior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

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

Trajectory sampling and finite-size effects in first-principles stopping power calculations

Abstract Real-time time-dependent density functional theory (TDDFT) is presently the most accurate available method for computing electronic stopping powers from first principles. However, obtaining application-relevant results often involves either costly averages over multiple calculations or ad hoc selection of a representative ion trajectory. We consider a broadly applicable, quantitative metric for evaluating and optimizing trajectories in this context. This methodology enables rigorous analysis of the failure modes of various common trajectory choices in crystalline materials. Although randomly selecting trajectories is common practice in stopping power calculations in solids, we show that nearly 30% of random trajectories in an FCC aluminum crystal will not representatively sample the material over the time and length scales feasibly simulated with TDDFT, and unrepresentative choices incur errors of up to 60%. We also show that finite-size effects depend on ion trajectory via “ouroboros” effects beyond the prevailing plasmon-based interpretation, and we propose a cost-reducing scheme to obtain converged results even when expensive core-electron contributions preclude large supercells. This work helps to mitigate poorly controlled approximations in first-principles stopping power calculations, allowing 1–2 order of magnitude cost reductions for obtaining representatively averaged and converged results.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Elevated nocturnal respiratory rates in the mitochondria of CAM plants: current knowledge and unanswered questions

Abstract Crassulacean acid metabolism (CAM) is a metabolic adaptation that has evolved convergently in 38 plant families to aid survival in water-limited niches. Whilst primarily considered a photosynthetic adaptation, CAM also has substantial consequences for nocturnal respiratory metabolism. Here, we outline the history, current state and future of nocturnal respiration research in CAM plants, with a particular focus on the energetics of nocturnal respiratory oxygen consumption. Throughout the 20th century, research interest in nocturnal respiration occurred alongside initial discoveries of CAM, although the energetic and mechanistic implications of nocturnal oxygen consumption and links to the operation of the CAM cycle were not fully understood. Recent flux balance analysis (FBA) models have provided new insights into the role that mitochondria play in the CAM cycle. Several FBA models have predicted that CAM requires elevated nocturnal respiratory rates, compared to C3 species, to power vacuolar malic acid accumulation. We provide physiological data, from the genus Clusia, to corroborate these modelling predictions, thereby reinforcing the importance of elevated nocturnal respiratory rates for CAM. Finally, we outline five unanswered questions pertaining to nocturnal respiration which must be addressed if we are to fully understand and utilize CAM plants in a hotter, drier world.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Comparative analysis of sublimation and thermal decomposition of TATB

This experimental study investigated the effects of confinement, starting mass, and heating rate on TATB thermal decomposition and sublimation using a combined Thermo-Gravimetric Analyzer and Differential Scanning Calorimetry (TGA/DSC) instrument. The confinement of volatile products was varied using different pinhole sizes with TGA/DSC pans. The measurements showed the open pan experiments without lids/pinholes resulted in complete sublimation of TATB between 320 °C and 360 °C. The heat of sublimation was determined to be 176 kJ/mol (42 kcal /mol), consistent with literature data obtained from other experimental techniques. The use of pinholes suppressed the sublimation of TATB such that the decrease in pinhole size resulted in 1) an increase in the enthalpy of reaction and an increase in the amount of carbonaceous material remaining at the end of decomposition, and 2) convergence of the two peak temperatures corresponding to maximum heat flow and maximum weight loss. Also, a transition from a two-exotherm thermal decomposition behavior towards a single-exotherm occurred as the pinhole size was decreased for a given starting mass or as the starting mass was increased for a given pinhole size. These results indicate the kinetics of TATB sublimation, TATB thermal decomposition, and gas diffusion out of a TGA/DSC pan can all compete and result in significantly different enthalpies, amounts of remaining materials, and peak temperatures depending on the pinhole size and starting mass used in the measurements. Furthermore, the results also indicate precise control of process variables (pinhole size, starting mass, and heating rate) in TGA/DSC measurements is required for thermal safety assessment of explosives.

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

Numerical discreteness errors in multispecies cosmological N -body simulations

ABSTRACT We present a detailed analysis of numerical discreteness errors in two-species, gravity-only, cosmological simulations using the density power spectrum as a diagnostic probe. In a simple set-up where both species are initialized with the same total matter transfer function, biased growth of power forms on small scales when the solver force resolution is finer than the mean interparticle separation. The artificial bias is more severe when individual density and velocity transfer functions are applied. In particular, significant large-scale offsets in power are measured between simulations with conventional offset grid initial conditions when compared against converged high-resolution results where the force resolution scale is matched to the interparticle separation. These offsets persist even when the cosmology is chosen so that the two particle species have the same mass, indicating that the error is sourced from discreteness in the total matter field as opposed to unequal particle mass. We further investigate two mitigation strategies to address discreteness errors: the frozen potential method and softened interspecies short-range forces. The former evolves particles under the approximately ‘frozen’ total matter potential in linear theory at early times, while the latter filters cross-species gravitational interactions on small scales in low-density regions. By modelling closer to the continuum limit, both mitigation strategies demonstrate considerable reductions in large-scale power spectrum offsets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Discrete modeling of a transformer with ALEGRA

We report progress on a task to model transformers in ALEGRA using the “Transient Magnetics” option. We specifically evaluate limits of the approach resolving individual coil wires. There are practical limits to the number of turns in a coil that can be numerically modeled, but calculated inductance can be scaled to the correct number of turns in a simple way. Our testing essentially confirmed this “turns scaling” hypothesis. We developed a conceptual transformer design, representative of practical designs of interest, and that focused our analysis. That design includes three coils wrapped around a rectangular ferromagnetic core. The secondary and tertiary coils have multiple layers. The tertiary has three layers of 13 turns each; the secondary has five layers of 44 turns; the primary has one layer of 20 turns. We validated the turns scaling of inductance for simple (one-layer) coils in air (no core) by comparison to available independent calculations for simple rectangular coils. These comparisons quantified the errors versus reduced number of turns modeled. For more than 3 turns, the errors are <5%. The magnetic field solver failed to converge (within 5000 iterations) for >10 turns. Including the core introduced some complications. It was necessary to capture the core surfaces in thin grid sheaths to minimize errors in computed magnetic energy. We do not yet have quantitative benchmarks with which to compare, but calculated results are qualitatively reasonable.

42 ENGINEERING↗
DOE OSTI2022

Peer-to-Peer Energy Trading under Network Constraints Based on Generalized Fast Dual Ascent

We report the wide deployment of renewable energy resources, combined with a more proactive demand-side management, is inducing a new paradigm in both power system operation and electricity market trading, which especially boosts the emergence of the peer-to-peer (P2P) market. A more flexible local market mechanism is highly desirable in response to fast changes in renewable power generation at the distribution network level. Moreover, large-scale implementation of P2P energy trading inevitably affects the secure and economic operation of the distribution network. This paper presents a new P2P electricity trading framework with distribution network security constraints considered using the generalized fast dual ascent method. First, an event-driven local P2P market framework is presented to facilitate short-term or immediate local energy transactions. Then, the sensitivity analysis of nodal voltage and network loss with respect to nodal power injections is used to evaluate the impacts of P2P transactions on the distribution network, which ensures the secure operation of the distribution system. Thereby, the external operational constraints are internalized, and the cost of P2P energy trading can be appropriately allocated in an endogenous way. Moreover, a generalized fast dual ascent method is employed to implement distributed market-clearing efficiently. Finally, numerical results indicate that the proposed model could guarantee secure operation of the distribution system with P2P energy trading, and the solution method enjoys good convergence performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Self-consistent analysis for the η c → γ γ process

The next-to-next-to-leading-order (NNLO) perturbative QCD (pQCD) predictions for both the decay width and the transition form factor in the η c → γ γ process, based on nonrelativistic QCD (NRQCD), deviate from precise experimental measurements. These significant discrepancies have cast doubt on the applicability of NRQCD to charmonium processes. In this paper, we analyze the η c → γ γ process by applying the principle of maximum conformality (PMC), a systematic method for eliminating renormalization scheme-and-scale ambiguities. The PMC renormalization scales are determined by absorbing the nonconformal β -terms that govern the behavior of the QCD running coupling via the renormalization group equation. We obtain the PMC scale Q ⋆ = 4.49 m c for the η c → γ γ decay width. Even after using the PMC method, the convergence of the pQCD series is still poor, which indicates the importance of uncalculated next-to-next-to-next-to-leading-order and higher-order terms. The resulting value for Γ η c → γ γ is in agreement with the Particle Data Group’s reported value of Γ η c → γ γ = 5.1 ± 0.4 keV within the bounds of uncertainties. Moreover, the transition form factor obtained using the PMC is also in good agreement with precise experimental measurements. The application of the PMC suggests a potential resolution to η c → γ γ puzzle and supports the applicability of NRQCD to charmonium processes. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

siRNA high throughput screening identifies regulators of chloropicrin and hydrogen fluoride injury in human corneal epithelial cell models

Corneal injuries induced by various toxicants result in similar clinical presentations such as corneal opacity and neovascularization. Many studies suggest that several weeks post-exposure a convergence of the molecular mechanisms drives these progressive pathologies. However, chemical agents vary in toxicological properties, and early molecular responses are anticipated to be somewhat dissimilar for different toxicants. In this report we chose 3120 targets from the Dharmacon Human Druggable genome to screen for chloropicrin (CP) and hydrogen fluoride (HF) corneal injury as we hypothesized that targets identified in vitro may be effective as therapeutic targets in future studies. Human immortalized corneal epithelial cells (SV40-HCEC) were used for screening. Cell viability and IL-8 were analyzed to down-select hits into validation studies, where multiplex cytokine analysis and high content analysis were performed to understand toxicant effect and target function. Some endpoints were also evaluated in a second human immortalized corneal epithelial cell line, TCEpi. Over 20 targets entered validation studies for CP and HF; of these, only three targets were shared: NR3C1, RELA, and KMT5A. These findings suggest that early molecular responses to different toxicants may be somewhat distinctive and present dissimilar targets for possible early intervention.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

PolyDODT: a macrocyclic elastomer with unusual properties

The effect of reaction conditions on the structure of poly(3,6-dioxa-1,8-octanedithiol) (polyDODT) made by Reversible Radical Recombination Polymerization (R3P) using triethylamine (TEA), H 2 O 2 and air was investigated. 800 MHz (1 and 2D) NMR was used to investigate the polymer structures. Sensitivity analysis provided direct evidence for high purity cyclic polyDODT up to M n ~ 100 000 g mol –1 . Here, comparative analysis by High Resolution Multidetector Size Exclusion Chromatography (SEC) using integrated data showed that the cycles had lower viscosity and were more compact (both for R g and R h ) than linear samples of similar molecular weight. However, differential data revealed unusual behavior. While lower molecular weight cyclic polymers eluted later and had lower intrinsic viscosity than their linear counterparts at the same molecular weight, at higher molecular weights the polymers showed strange behavior: both the diffusion coefficient measured by Quasielastic Light Scattering (QELS) and Mark–Houwink–Sakurada plots of intrinsic viscosity for linear and cyclic polyDODT were found to converge. R3P, an aqueous based “green” method is capable of producing polymers at the 10–100 g scale in the lab, which will allow more detailed studies of this new class of biodegradable elastomers so further experimentation can be performed to elucidate the reasons for the unusual findings.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Optimization using pathwise algorithmic derivatives of electromagnetic shower simulations

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Effects of initial conditions and transport on ram pressure, Mach number, and uniformity for plasma liner formation and implosion

Plasma jet driven magneto-inertial fusion involves the merging of high Mach number plasma jets in a spherically converging geometry to form an imploding plasma liner. A three-dimensional smooth particle hydrodynamic simulation tool was used to study the number of plasma jets, jet density, and implosion velocity on the effects of formation and subsequent implosion of a plasma liner. The effects of different transport physics are presented. Thermal radiation and two-temperature modeling were found to have the most significant impact on the ram pressure. Solid angle averaged ram pressure and Mach number degradation were of particular interest in this study. Higher jet numbers were found to result in higher peak ram pressure. A spherical harmonic analysis was conducted for each case, providing a quantitative measure of the non-uniformity of the liner during the implosion process and how non-uniformities are affected by the different parameters explored in this study.

Schillo, Kevin↗
DOE OSTI2020

Large-scale genome sequencing of mycorrhizal fungi provides insights into the early evolution of symbiotic traits

Mycorrhizal fungi are mutualists that play crucial roles in nutrient acquisition in terrestrial ecosystems. Mycorrhizal symbioses arose repeatedly across multiple lineages of Mucoromycotina, Ascomycota, and Basidiomycota. Considerable variation exists in the capacity of mycorrhizal fungi to acquire carbon from soil organic matter. Here, we present a combined analysis of 135 fungal genomes from 73 saprotrophic, endophytic and pathogenic species, and 62 mycorrhizal species, including 29 new mycorrhizal genomes. This study samples ecologically dominant fungal guilds for which there were previously no symbiotic genomes available, including ectomycorrhizal Russulales, Thelephorales and Cantharellales. Our analyses show that transitions from saprotrophy to symbiosis involve (1) widespread losses of degrading enzymes acting on lignin and cellulose, (2) co-option of genes present in saprotrophic ancestors to fulfill new symbiotic functions, (3) diversification of novel, lineage-specific symbiosis-induced genes, (4) proliferation of transposable elements and (5) divergent genetic innovations underlying the convergent origins of the ectomycorrhizal guild.

59 BASIC BIOLOGICAL SCIENCES↗