Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Component sizing”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Development of a robust damage analysis scheme at the National Ignition Facility

The onset of laser-induced damage in optical materials is a limiting factor in the design and operation of most high-energy laser systems. As such, significant effort has been dedicated to developing laser damage testing protocols and procedures to inform laser system design and operating limits. These tests typically rely on multiple laser exposures for statistical validation. Historically, small beam single fluence (N / 1) or ramped fluence (R / 1) tests have been used to quantify the “laser-induced damage threshold” of a material. However, due to various distributions in damage precursor populations, the laser fluence at which damage occurs for a given sample varies with beam size and is therefore better described as the onset of observed damage for that testing geometry. In conclusion, we document the development of irradiation, measurement, and analysis methods for damage testing resulting in damage initiation density as a function of incident fluence or ρ(ϕ) which is suitable for both small and large beam area testing and more conducive to extrapolating from one to the other.

Beam diameter↗
DOE OSTI2020

Infrared Optical Constants from Pressed Pellets of Powders: I. Improved n and k Values of (NH4)(2)SO(4)from Single-Angle Reflectance

In combination with other parameters, the real, n(?), and imaginary, k(?), components of the complex refractive index, = n + ik, can be used to simulate the optical properties of a material in different forms, e.g. its infrared spectra. But obtaining reliable values of the optical constants n/k for solid materials is challenging due to the lack of optical quality specimens, usually crystals, large enough to measure. An alternative to crystals is to press the powder into a uniform disk. We have produced pellets from ammonium sulfate, (NH4)2SO4, powder and derived the pellets’ n and k values via single-angle reflectance using a specular reflectance device in combination with a Fourier transform infrared spectrometer. The single-angle technique measures amplitude of light reflected from the material as a function of wavelength over a wide spectral domain; the optical constants are determined from the reflectance data using the Kramers-Kronig relationship. We investigate several parameters associated with the pellets and pellet formation and their effects upon delivering the most reliable n/k values. Parameters studied include pellet diameter, mass and density (void space), drying, grinding, sieving and particle size in the pellet formation, as well as pressing pressure and duration. Of these parameters, using size-selected mixtures of dried, small (< 100 µm) particles and pressing at = 10 tons for at least 30 minutes were found key to forming highly reflective samples. Comparison to two sets of previous literature n(?) and k(?) values obtained from crystalline (NH4)2SO4 both as crystal reflectance as well as extinction spectra of aerosols measured in a flow tube show reasonable agreement, but suggests the present values, as confirmed from two independent techniques, represent a substantial improvement for n/k values for (NH4)2SO4, also demonstrating promise to measure the optical constants of other materials.

Optical constants, index of refraction, single-ang↗
DOE OSTI2023

Deep learning closure models for large-eddy simulation of flows around bluff bodies

Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications.

Mechanics↗
DOE OSTI2023

Behavior and properties of helium cavities in ductile-phase-toughened tungsten irradiated with He + ions at an elevated temperature

Ductile-phase toughened tungsten (DPT W) composites have demonstrated a great potential for plasma-facing components in fusion power systems. These materials can retain the outstanding thermomechanical properties of W while significantly enhancing the fracture toughness. The DPT W composite in this study consists of W particles (88 wt.%) embedded in a ductile NiFeW (12 wt.%) solution matrix. Irradiation of the composite was performed with 90 keV He + ions to 1.0 × 10 17 He + /cm 2 at 973 K. The He cavities formed in both W and NiFeW phases are examined using convergent-beam scanning transmission electron microscopy (CB-STEM) at an atomic-level resolution. Full-range depth profiles of the He cavity diameters and number densities are analyzed. The results show that the average cavity diameter in NiFeW is ~2.7 times that in W, while the number density is an order of magnitude lower. A quantitative analysis of the He density and pressure in the cavities is also achieved using STEM electron energy loss spectroscopy (STEM-EELS) mapping. The He atoms in the cavities in NiFeW are found to be under-pressurized, while those in similar-sized cavities in W are not detectable. Further, the He energy-loss shift exhibits a linear relationship with the He density in cavities in NiFeW, with the slope comparable to that for martensitic steel.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Underwater unexploded ordnance discrimination based on intrinsic target polarizabilities – A case study

Seabed unexploded ordnance that resulted partly from the high failure rate among munitions from more than 80 years ago and from decades of military training and testing of weapons systems poses an increasing concern all around the world. Although existing magnetic systems can detect clusters of debris, they are not able to tell whether a munition is still intact requiring special removal (e.g. in situ detonation) or is harmless scrap metal. The marine environment poses unique challenges, and transferring knowledge and approaches from land to a marine environment has not been easy and straightforward. On land, the background soil conductivity is much lower than the conductivity of the unexploded ordnance and the electromagnetic response of a target is essentially the same as that in free space. For those frequencies required for target characterization in the marine environment, the seawater response must be accounted for and removed from the measurements. The system developed for this study uses fields from three orthogonal transmitters to illuminate the target and four three-component receivers to measure the signal arranged in a configuration that inherently cancels the system's response due to the enclosing seawater, the sea–bottom interface and the air–sea interface for shallow deployments. The system was tested as a cued system on land and underwater in San Francisco Bay – it was mounted on a simple platform on top of a support structure that extended 1 m below and allowed the diver to place metal objects to a specific location even in low-visibility conditions. The measurements were stable and repeatable. Furthermore, target responses estimated from marine measurements matched those from land acquisition, confirming that the seawater and air–sea interface responses were removed successfully. Thirty-six channels of normalized induction responses were used for the classification, which was done by estimating the target principal dipole polarizabilities. Our results demonstrated that the system can resolve the intrinsic polarizabilities of the target, with clear distinctions between those of symmetric intact unexploded ordnance and irregular scrap metal. The prototype system was able to classify an object based on its size, shape and metal content and correctly estimate its location and orientation.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Computational Algorithms for Unit Commitment with AC Power Flows (Final Report)

Security-constrained unit commitment (SCUC) is a key component in power system operations. When AC power flow constraints are considered in the SCUC model (AC-SCUC), the problem becomes extremely difficult due to its discrete and non-convex nature, as described in “Grid Optimization Competition Challenge 3 Problem Formulation (GOCC)”. There are four main challenges: (i) Discrete decisions regarding unit online/offline status and start-up/shut-down procedures for every single unit. The number of discrete decision variables increases considerably when a system integrates multiple generators; (ii) Configuration-based combined-cycle formulations, and multi-commodity models that include ramping products, spin/non-spin products, and regulation up/down products. The combined-cycle units introduce additional discrete decision variables and auxiliary service products further complicate the model by connecting multi-commodity products’ continuous and discrete variables; (iii) SCUC models with AC power flow constraints are far more complex due to massive bilinear terms in the large-scale nonlinear power balance equations. The nonlinear power balance equations are further complicated by the discrete step control variables of shunts; (iv) N − 1 contingency analysis. The size of the model increases linearly with the number of contingencies considered, greatly increasing the size of the optimization model. Accordingly, there is an emergent need to develop a robust algorithm capable of deriving a high-quality solution in a short time and passing through contingency tests simultaneously. In this project, we explore innovative techniques to address this challenging problem by integrating advanced polyhedral theory, approximation methods, relaxation strategies, decomposition techniques, and parallel computing. Each technique approaches the problem from a different perspective, leveraging its specific strengths to tackle distinct challenges. Each individual method has demonstrated its effectiveness in the PI’s previous research. Their integration is expected to significantly reduce the computational time required to solve the proposed complex problem. Successful completion of this project has the potential to transform the industry by enhancing optimization solvers capable of handling large-scale day-ahead energy market clearing models within strict time constraints, while incorporating AC power flow constraints. This advancement will lead to reduced overall generation costs and, consequently, increased social welfare.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Predictive models enhance feedstock quality of corn stover via air classification

Feedstock heterogeneity is a fundamental obstacle to cost-competitive biobased products. Agricultural products like corn stover have anatomical components that vary in their chemical composition, mechanical properties, structure, and response to chemical and biological treatments. A technique that can enrich streams in select anatomical fractions would allow a tailored deconstruction approach to increase overall process efficiency. Air classification can be leveraged for such refining, however, fundamental characterization and understanding of the particle properties that underly the physics of air classification are only modestly documented. Here, we determine fundamental particle properties including mass-to-area ratio, drag coefficient, and partition velocity that describe how anatomical tissues of corn stover behave during air classification. In this work, mass-to-area ratios of anatomical tissues vary by nearly two orders of magnitude from 2.3 mg/mm 2 for cob to 0.04 mg/mm 2 for leaf. Drag coefficients of longer, fibrous materials (i.e., rind, husk, and sheath) are shown to correlate with particle area (p-value < 0.001) whereas granular tissues (i.e., cob, pith, and leaf) correlate better with mass-to-area ratio (p-values < 0.001). When compared to experimental observations, a simulated two-stage air classification and size reduction scenario predicts the overall partitioning of anatomical tissues within 15% for pith, husk, rind, and cob tissues. The model predicts an air-classified fraction preferentially enriched in cob (purity = 20%), rind (purity = 74%), and pith (purity = 4.5%) with a mass yield of 47%. Empirical relations for these properties can be used to predict the partitioning of corn stover during air classification based on anatomical type and size.

09 BIOMASS FUELS↗
DOE OSTI2023

Hardware Specification and Reference Design for the Low-Cost, Interoperable, User-Centric, Supervisory Controller Kit for Small and Medium Size Commercial Buildings

Commercial buildings are responsible for approximately 20 percent of the total United States energy consumption and greenhouse gas emissions. Over 85 percent of these buildings lack building automation systems to manage the various building systems they have. Many of these buildings are small (< 50,000 square feet), underserved, and use rooftop units (RTUs) for heating, ventilation, and air-conditioning needs. Because these buildings lack proper control systems, they have several operational deficiencies that lead to excess energy consumption. Studies have shown that managing the RTU’s heating and cooling setpoints, schedules, setbacks, and optimal start times can result in 20 to 25 percent reduction in electricity consumption in small and medium commercial buildings (SMBs). In addition, improving the demand flexibility of these buildings will result in additional cost savings for the building owner. To address the needs of the SMBs, the Department of Energy’s Building Technologies Office jointly funded Pacific Northwest National Laboratory (PNNL) and Oak Ridge National Laboratory (ORNL) to design, develop, test, and validate a low-cost, interoperable, user-centric, retrofit supervisory controller kit (SC-SMB) that can be used to continuously optimize energy consumption and deliver demand flexibility of SMBs, including all-electric buildings, and provide a means for maximizing decarbonization benefits. The team also includes industry partners Edo and Intellimation. Previously, the team drafted an SC-SMB system specification document (Goodman et al. 2023). This document describes a reference design of a low-cost supervisory control system for SMBs, including the various hardware components, software system, and an example implementation of the reference design. The final release of the reference design is planned for December 2025. Section 2.0 of the report documents the relevant building types that the SC-SMB system is suitable for. Section 3.0 documents the various hardware components, their function, cost, whether they are off-the-shelve, support standard communication, etc., and the software system. An example deployment in a 10,000-sf building with six rooftop units, a hot water heater, solar, and storage is described in Section 4.0. The planned next steps are described in Section 5.0, and references are listed in Section 6.0.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Optimized manufacturing process for multilayer two-dimensional focusing mirrors in laboratory X-ray applications

Recent advances in laboratory X-ray applications require high-performance optical components that achieve exceptional imaging resolution and beam uniformity within compact experimental setups. Montel mirrors have become a preferred solution due to their unique dual-reflection focusing mechanism and a space-efficient design. Here, in this study, we present an effective manufacturing process for producing Montel mirrors tailored to focus laboratory X-ray beams. The mirrors were fabricated from single-crystal silicon substrates, chosen for their high mechanical stability and compatibility with precision polishing techniques. Our approach begins with the integration of a deterministic chemo-mechanical polishing (CMP)-based pre-shaping step followed by ion beam figuring (IBF), significantly improving manufacturing efficiency. Subsequently, our custom-developed advanced metrology and IBF techniques were employed for fabricating an off-axis, elliptical cylinder Montel mirror system with a 6-mrad total slope, with stringent optical specifications. While post-IBF processes, including multilayer coating, dicing, and gluing, introduced minor surface errors, yet their impact on performance remained negligible. The Montel mirrors manufactured with the optimized process exhibited significantly improved beam uniformity and a reduced focal spot size. These findings validate our approach as a viable solution for high-precision Montel mirror fabrication and facilitate further advancements in laboratory X-ray applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Industrial Assessment Center

Established in 1990, San Diego State University’s (SDSU) Industrial Assessment Center (IAC) is proud of its years of service. During this period, it has served over 620 small and medium-sized manufacturing plants in Southern California. SDSU/IAC’s efforts to transfer state-of-the-art technologies to industry have increased revenues, cultivated creativity, improved energy efficiencies, and benefited the environment. The Center has contributed to the region's economic growth and stability by assisting small and medium size companies to better compete in the global market. It has helped mitigate climate change by reducing greenhouse gas emissions. IAC activities have fostered productive relationships between the University and local industry, assisted industrial sectors to improve their energy efficiency and enhance their manufacturing productivity, in turn impacting the material and working conditions of their employees. In addition to financial savings and environmental benefits, we have trained tens of students who became energy specialists in various companies. Thus, a substantial benefit of the IAC has been the ongoing training of engineering faculty and students. All IAC graduates were offered jobs before or within weeks of their graduation. The activities of the SDSU/IAC have expanded the institutional expertise of the College and improved the knowledge base of the faculties involved leading to several related publications, master’s theses, and senior student projects. Significant number of peer-reviewed publications of the IAC director at SDSU have greatly benefitted from the experience of the Center. As a result of this extensive exposure to manufacturing processes, the SDSU/IAC has grown to be an integral component of SDSU’s engineering research and training. We have successfully built upon these established achievements and academic excellence. IAC service to industry is particularly vital in Southern California, a region with one of the highest manufacturing concentrations in the country. SDSU/IAC has understood and implemented the overall objectives of DOE’s IAC program and guidelines except for the pandemic years when the country’s manufacturing sector was put in dire stress. In addition to student training and service to industry, IAC’s contribution to state and local governments as well as utility companies to assess energy policies and design rebate and incentive strategies cannot be undermined.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

WM26 Paper Multi-Robot Collaboration for Hazardous Environments

Hazardous nuclear and industrial facilities are rarely designed for robots. Work in these domains demand precise manipulation and robust mobility in cluttered, constrained spaces where off-the-shelf platforms struggle and “one-size-fits-all” machines become costly and complex. Idaho National Laboratory (INL) is developing an autonomous, multi-robot inspection system that coordinates task-specific platforms rather than relying on a single omni-tool robot. An electric truck serves as a power and compute hub for a custom manipulator co-developed with Florida International University (FIU), a commercial mini crawler, a pan–tilt–zoom camera, and a Nexxis Argus LiDAR mapping system. Working in concert, these robots generate spatial, radiation, and temperature maps of the pit environments at the Hanford Waste Tank Farms. These systems will capture visual records and environmental telemetry to allow for analysis post inspection. The system architecture uses Robot Operating System 2 (ROS 2) for publish/subscribe integration, NVIDIA Isaac Sim and Unity for simulation and visualization, and algorithms such as NVBlox to fuse data into unified 3D overlays. This robot-agnostic approach reduces operator burden by enabling autonomy across heterogeneous platforms and lets each robot be used where it is strongest. Having autonomous functions means operators don’t have to fully control multiple different components. The ease of use could allow for more widespread adoption of advanced robotics at waste management sites that see continued use. By coordinating simpler, purpose-built mechanisms, the approach lowers design and manufacturing complexity, reduces capital risk in contaminated settings, and improves controllability for complex inspection and manipulation tasks. We present the architecture, early results, and lessons learned from building and deploying this coordinated multi-robot system, with the goal of accelerating safe, cost-effective adoption of advanced robotics at waste-management sites.

42 - ENGINEERING↗
DOE OSTI2025

Flexible PCB Windings Size Optimization for Winding AC Resistance Minimization Under the Sinusoidal Voltage Regime

In this paper, optimization of the winding AC resistance of the flexible printed circuit board (FPCB) wire windings is performed. A one-dimensional model of the FPCB winding is introduced and an equation for the FPCB winding AC resistance and its low and medium frequency approximation are derived. The approximate FPCB winding AC resistance equation is used to derive optimum thickness of the FPCB winding conductor thickness at which minimum of winding AC resistance (global optimum) is achieved. The Finite Element Method analysis and experimental verification of derived equations (winding resistance, impedance, power loss, and temperature measurements) is performed in order validate derived model and winding resistance equation.

42 ENGINEERING↗
DOE OSTI2021

Structural Characterization of the Reaction and Substrate Specificity Mechanisms of Pathogenic Fungal Acetyl-CoA Synthetases

Acetyl CoA synthetases (ACSs) are Acyl-CoA/NRPS/Luciferase (ANL) superfamily enzymes that couple acetate with CoA to generate acetyl CoA, a key component of central carbon metabolism in eukaryotes and prokaryotes. Normal mammalian cells are not dependent on ACSs, while tumor cells, fungi, and parasites rely on acetate as a precursor for acetyl CoA. Consequently, ACSs have emerged as a potential drug target. As part of a program to develop antifungal ACS inhibitors, we characterized fungal ACSs from five diverse human fungal pathogens using biochemical and structural studies. ACSs catalyze a two-step reaction involving adenylation of acetate followed by thioesterification with CoA. Our structural studies captured each step of these two half-reactions including the acetyl-adenylate intermediate of the first half-reaction in both the adenylation conformation and the thioesterification conformation and thus provide a detailed picture of the reaction mechanism. We also used a systematic series of increasingly larger alkyl adenosine esters as chemical probes to characterize the structural basis of the exquisite ACS specificity for acetate over larger carboxylic acid substrates. Consistent with previous biochemical and genetic data for other enzymes, structures of fungal ACSs with these probes bound show that a key tryptophan residue limits the size of the alkyl binding site and forces larger alkyl chains to adopt high energy conformers, disfavoring their efficient binding. Together, our analysis provides highly detailed structural models for both the reaction mechanism and substrate specificity that should be useful in designing selective inhibitors of eukaryotic ACSs as potential anticancer, antifungal, and antiparasitic drugs.

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

Experimental Protocols for Studying Organic Non-aqueous Redox Flow Batteries

We report that Redox flow batteries (RFBs) are promising devices for grid-scale energy storage due to the decoupling of power and energy, which can be independently scaled by the electrode area and storage tank size, respectively. To date, only aqueous RFBs, such as the vanadium RFB, have been implemented commercially. Nevertheless, the limited energy densities and high-cost materials may preclude their wider market penetration. Organic non-aqueous redox-flow batteries (O-NRFBs), which utilize redox-active organic molecules (ROMs), have been offered as an attractive alternative. Many possible advantages include the use of earth-abundant elements (C, H, N, O, S, F), wherein the ROMs can be prepared from low-cost and sustainable materials. Additionally, a large variety of electroactive moieties are accessible as building blocks, providing a synthetic platform to tune the properties of ROMs through rational design. As a consequence, the development of novel ROMs has attracted researchers with diverse backgrounds, prompting remarkable innovations in the past decade. However, consistency in experimental protocols (e.g., electrochemical methods, cycling stability, experimental conditions) has not coincided with this uptick in research, leading to sometimes ambiguous and incomparable results. This is further convoluted by the complexity innate to battery development, such as cell design, detection, and characterization of reactions and active components. The O-NRFB application imposes stringent requirements on the physical or/and chemical processes involved in electrochemical cycling. Yet such considerations are often overlooked. Further, while quantum mechanical calculations provide a convenient tool for evaluating and selecting ROM candidates, this simulation is not always performed. Thus, in this Energy Focus, we detail a means to standardize experimental protocols for studies of O-NRFBs and suggest practices to facilitate fundamental understanding and development of ROMs.

25 ENERGY STORAGE↗
DOE OSTI2022

Using post-processing heat treatments to elucidate precipitate strengthening of additively manufactured superalloy 718

The poor machinability and extensive work hardening of Ni-based superalloys makes additive manufacturing an attractive option for producing geometrically complex components with distinct microstructures. Although previous studies show recovery of high strength at room temperature, very few studies demonstrate successful properties at elevated temperatures required for industrial applications. Here, the objectives of this study are to present a post-build heat treatment for high strength across a wide temperature range, determine the strength contribution of nanoscale precipitating phases to the overall mechanical properties of superalloy 718, and from these, provide a comprehensive microstructure-property relationship for wrought and AM 718 to guide efforts to simulate the properties of AM components. Laser powder bed fusion–produced superalloy 718 was characterized at multiple length scales using scanning electron microscopy and transmission electron microscopy in the as-built condition and with multiple heat treatments designed to form combinations of γ', γ", and δ precipitates. Uniaxial tensile tests performed from room temperature to 600 °C on subsize specimens determined the yield strength, elastic modulus, ultimate tensile strength, fracture stress, and uniform elongation. Precipitates in this work proved to be weak barriers to dislocation motion through a dispersed barrier model, but they provided strength to the alloy through their consistent high density. The relative contribution to the yield strength from γ" remained consistent between 48% and 57% of the total strength up to 600 °C, the primary influence on the high temperature strength of superalloy 718. The strength factors for γ" and δ precipitates were found to trend inversely with tensile test temperature and may be attributable to the differences in precipitate coherency. A post-build heat treatment is recommended to maintain high strength at elevated temperatures. A quantitative microstructure-property relationship, dependent on precipitate size, density, and morphology, was derived and can estimate the yield strength across a wide temperature range applicable to the operational regimes for superalloy 718.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Enabling machine learning-ready HPC ensembles with Merlin

With the growing complexity of computational and experimental facilities, many scientific researchers are turning to machine learning (ML) techniques to analyze large scale ensemble data. With complexities such as multi-component workflows, heterogeneous machine architectures, parallel file systems, and batch scheduling, care must be taken to facilitate this analysis in a high performance computing (HPC) environment. Here, we present Merlin, a workflow framework to enable large ML-friendly ensembles of scientific HPC simulations. By augmenting traditional HPC with distributed compute technologies, Merlin aims to lower the barrier for scientific subject matter experts to incorporate ML into their analysis. As a producer–consumer workflow model, Merlin enables multi-machine, cross-batch job, dynamically allocated yet persistent workflows capable of utilizing surge-compute resources. Key features of Merlin are a flexible HPC-centric interface, low per-task overhead, multi-tiered fault recovery, and a hierarchical sampling algorithm that allows for $\mathscr{O}$(N) task execution and $\mathscr{O}$(N ln N) task queuing to ensembles of millions of tasks. In addition to Merlin’s design, we test the algorithm’s performance in an HPC center and demonstrate the ability to enqueue 40 million simulations in 100 s, with a 30 millisecond per-task overhead that is independent of ensemble size. Finally, we describe some example applications that Merlin has enabled on leadership-class HPC resources, such as the ML-augmented optimization of nuclear fusion experiments and the calibration of infectious disease models to study the progression of and possible mitigation strategies for COVID-19.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Coarse woody debris decomposition assessment tool: Model development and sensitivity analysis

Coarse woody debris (CWD) is an important component in forests, hosting a variety of organisms that have critical roles in nutrient cycling and carbon (C) storage. We developed a process-based model using literature, field observations, and expert knowledge to assess woody debris decomposition in forests and the movement of wood C into the soil and atmosphere. The sensitivity analysis was conducted against the primary ecological drivers (wood properties and ambient conditions) used as model inputs. The analysis used eighty-nine climate datasets from North America, from tropical (14.2° N) to boreal (65.0° N) zones, with large ranges in annual mean temperature (26.5°C in tropical to -11.8°C in boreal), annual precipitation (6,143 to 181 mm), annual snowfall (0 to 612 kg m -2 ), and altitude (3 to 2,824 m above mean see level). The sensitivity analysis showed that CWD decomposition was strongly affected by climate, geographical location and altitude, which together regulate the activity of both microbial and invertebrate wood-decomposers. CWD decomposition rate increased with increments in temperature and precipitation, but decreased with increases in latitude and altitude. CWD decomposition was also sensitive to wood size, density, position (standing vs downed), and tree species. The sensitivity analysis showed that fungi are the most important decomposers of woody debris, accounting for over 50% mass loss in nearly all climatic zones in North America. The model includes invertebrate decomposers, focusing mostly on termites, which can have an important role in CWD decomposition in tropical and some subtropical regions. The role of termites in woody debris decomposition varied widely, between 0 and 40%, from temperate areas to tropical regions. Woody debris decomposition rates simulated for eighty-nine locations in North America were within the published range of woody debris decomposition rates for regions in northern hemisphere from 1.6° N to 68.3° N and in Australia.

54 ENVIRONMENTAL SCIENCES↗