SEARCH · Engineering Papers

Results for “model scaling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Computer vision models enable mixed linear modeling to predict arbuscular mycorrhizal fungal colonization using fungal morphology

Abstract The presence of Arbuscular Mycorrhizal Fungi (AMF) in vascular land plant roots is one of the most ancient of symbioses supporting nitrogen and phosphorus exchange for photosynthetically derived carbon. Here we provide a multi-scale modeling approach to predict AMF colonization of a worldwide crop from a Recombinant Inbred Line (RIL) population derived from Sorghum bicolor and S. propinquum . The high-throughput phenotyping methods of fungal structures here rely on a Mask Region-based Convolutional Neural Network (Mask R-CNN) in computer vision for pixel-wise fungal structure segmentations and mixed linear models to explore the relations of AMF colonization, root niche, and fungal structure allocation. Models proposed capture over 95% of the variation in AMF colonization as a function of root niche and relative abundance of fungal structures in each plant. Arbuscule allocation is a significant predictor of AMF colonization among sibling plants. Arbuscules and extraradical hyphae implicated in nutrient exchange predict highest AMF colonization in the top root section. Our work demonstrates that deep learning can be used by the community for the high-throughput phenotyping of AMF in plant roots. Mixed linear modeling provides a framework for testing hypotheses about AMF colonization phenotypes as a function of root niche and fungal structure allocations.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

A data-driven multiscale model for reactive wetting simulations

Here, we describe a data-driven, multiscale technique to model reactive wetting of a silver–aluminum alloy on a Kovar™ (Fe-Ni-Co alloy) surface. We employ molecular dynamics simulations to elucidate the dependence of surface tension and wetting angle on the drop’s composition and temperature. A design of computational experiments is used to efficiently generate training data of surface tension and wetting angle from a limited number of molecular dynamics simulations. The simulation results are used to parameterize models of the material’s wetting properties and compute the uncertainty in the models due to limited data. The data-driven models are incorporated into an engineering-scale (continuum) model of a silver–aluminum sessile drop on a Kovar™ substrate. Model predictions of the wetting angle are compared with experiments of pure silver spreading on Kovar™ to quantify the model-form errors introduced by the limited training data versus the simplifications inherent in the molecular dynamics simulations. The paper presents innovations in the determination of “convergence” of noisy MD simulations before they are used to extract the wetting angle and surface tension, and the construction of their models which approximate physio-chemical processes that are left unresolved by the engineering-scale model. Together, these constitute a multiscale approach that integrates molecular-scale information into continuum scale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

PathTrace and MPVEASI: A Path Analysis Comparative Validation Study

Developed in 2018, PathTrace is a software package built with the intention of making path analysis simple and intuitive. PathTrace is a top-down pathway analysis software where a user is able to explore vulnerable pathways into a facility. The intention of utilizing a software tool like PathTrace is to characterize an existing physical protection system (PPS) and to upgrade the system to achieve a high level of response interruption, or probability of interruption (P I ) of the adversary. There are four steps for conducting path analysis using PathTrace. The first step is to identify an image to use to build the model and scale the model within PathTrace using a section of known distance (wall or fence perimeter, for example). The scaling process will produce a grid of cells through which the user is able to build a model. The second step is to fill out the grid of cells with four categories of materials: Barriers, Detection Areas, Jumps, and Targets. These materials apply associated delay and detection values to the cells in which they are applied. The third step is to represent the adversary and response forces. The adversaries are represented by their capabilities in interacting with the materials identified in step two, and the response is represented by how quickly they will be able to respond to an adversary attack. Finally, the user is able to take all of the information from the previous three steps and perform a Most Vulnerable Path (MVP) analysis. In this stage, the user is able to visualize vulnerable adversary pathways and reason about how to upgrade these pathways to provide a high level of P I .

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Numerical wave modeling for operational and survival analyses of wave energy converters at the US Navy Wave Energy Test Site in Hawaii

The US Navy Wave Energy Test Site (WETS) is currently operating three grid-connected berths off Marine Corps Base Hawaii in support of technological development through scaled-model testing and pre-commercial prototyping. The location on the east shore of Oahu is open to year-round wind waves and seasonal north swells as well as severe seas from subtropical systems and hurricanes. We have assembled a spectral wave model system comprising WAVEWATCH III and SWAN on a hierarchy of global, regional, and nearshore computational grids. With wind forcing from global forecast and rean-alysis datasets as well as their regional downscaling, the system produces operational 7-day wave forecasts and a long-term hindcast. The daily forecasts, validated with real-time buoy measurements, facilitate safe deployment, operation, and retrieval of wave energy converters (WECs). The wave hindcast from 1979 to 2017 defines the intra- and inter-annual variations for statistical analysis of sea states and energy resources. Due to limited occurrences of hurricanes in Hawaii waters, their wind forcing is determined parametrically with simulated tracks from climate model downscaling. Finally, systematic analysis of the hindcast and simulation datasets provides significant wave heights from trade winds, swells, subtropical systems, and hurricanes with return periods up to 100 years for WECs survival analysis.

16 TIDAL AND WAVE POWER↗

DOE OSTI2024

Ground Heat Flux Reconstruction Using Bayesian Uncertainty Quantification Machinery and Surrogate Modeling

Ground heat flux (G 0 ) is a key component of the land-surface energy balance of high-latitude regions. Despite its crucial role in controlling permafrost degradation due to global warming, G 0 is sparsely measured and not well represented in the outputs of global scale model simulation. In this study, an analytical heat transfer model is tested to reconstruct G 0 across seasons using soil temperature series from field measurements, Global Climate Model, and climate reanalysis outputs. The probability density functions of ground heat flux and of model parameters are inferred using available G 0 data (measured or modeled) for snow-free period as a reference. When observed G 0 is not available, a numerical model is applied using estimates of surface heat flux (dependent on parameters) as the top boundary condition. These estimates (and thus the corresponding parameters) are verified by comparing the distributions of simulated and measured soil temperature at several depths. Aided by state-of-the-art uncertainty quantification methods, the developed G 0 reconstruction approach provides novel means for assessing the probabilistic structure of the ground heat flux for regional permafrost change studies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Investigating the Evolution of Ice Particle Distributions in Mixed-Phase Clouds

The goal of this project is to conduct modeling studies that focus on the processes that control the macrophysical and microphysical properties of mixed phase clouds, such as the partitioning of water phase and the processes that promote precipitation and govern cloud evolution. The process specific to this project is the collection of vapor-grown non-spherical ice crystals and snow hydrometeors through aggregation, as uncertainty remains in the evolution of particle properties as they aggregate. To aid in the understanding of ice particle evolution from pristine monomers to aggregates, the aggregation process is depicted using an offline simulator for ice crystal collection, namely the Ice Particle and Aggregate Simulator (IPAS). Of particular interest in the context of this work is the influence of monomer habit on aggregation and subsequently on cloud microphysical structure. This overarching goal has been completed by conducting modeling studies using the Adaptive Habit Model (AHM) in mixed-phase cloud systems, which predicts and evolves ice shape. While the AHM is designed with physical detail in mind, those details may be lost on larger spatial or temporal scales. Hence, how the contributions of particle growth and collection control the evolution of ice particle size distributions (PSDs) and redistribute mass within the cloud and at the surface has been investigated. Analyses include microphysical sensitivity investigations that inform on future improvement capabilities of microphysical parameterizations in larger scale models. Further, improvements of processes within ice clouds have resulted from this work as the environmental processes controlling ice crystal growth, microphysical processes, and precipitation are inherently integrated within this investigation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Design Principles to Govern Electrode Fabrication for the Lithium Trivanadate Cathode

A full depth of discharge mathematical model for the lithium trivanadate cathode, considering lithiation of the layered α-phase, phase change, and lithiation of the rock-salt like β-phase at lower potentials, is developed. The coupled electrode-scale and crystal-scale model is fit to electrochemical data, and additionally validated with operando EDXRD. There is good agreement between the simulated and measured spatial variation of the volume fraction of the β-phase. This mathematical model is used to guide electrode fabrication, accounting for both ionic and electronic transport effects. Values of three design parameters—electrode thickness, porosity, and volume fraction of conductor—are identified, and the sensitivity of the energy density to these design parameters is quantified. The model is also used to investigate electrode design to create electrodes that deliver the maximum achievable energy density under the constraint that the α to β-phase transition is avoided, since phase change has been demonstrated to reduce cycle life. The energy density sacrificed to avoid phase change decreases at higher discharge rates, but the target values for electrode fabrication remain the same as those when optimizing the electrode for the full depth of discharge.

25 ENERGY STORAGE↗

DOE OSTI2021

Networked HIL Simulation System for Modeling Large-scale Power Systems

This paper presents a network hardware-in-the-loop (HIL) simulation system for modeling large-scale power systems. Researchers have developed many HIL test systems for power systems in recent years. Those test systems can model both microsecond-level dynamic responses of power electronic systems and millisecond-level transients of transmission and distribution grids. By integrating individual HIL test systems into a network of HIL test systems, we can create large-scale power grid digital twins with flexible structures at required modeling resolution that fits for a wide range of system operating conditions. This will not only significantly reduce the need for field tests when developing new technologies but also greatly shorten the model development cycle. In this paper, we present a networked OPAL-RT based HIL test system for developing transmission-distribution coordinative Volt-VAR regulation technologies as an example to illustrate system setups, communication requirements among different HIL test systems, and connection mechanisms. Impacts of communication delays, information exchange cycles, and computing delays are illustrated. Simulation results show that the performance of a networked HIL test system is satisfactory.

digital twin, distribution system, co-simulation, ↗

DOE OSTI2022

Parked aeroelastic field rotor response for a 20% scaled demonstrator of a 13‐MW downwind turbine

Abstract Aeroelastic parked testing of a unique downwind two‐bladed subscale rotor was completed to characterize the response of an extreme‐scale 13‐MW turbine in high‐wind parked conditions. A 20% geometric scaling was used resulting in scaled 20‐m‐long blades, whose structural and stiffness properties were designed using aeroelastic scaling to replicate the nondimensional structural aeroelastic deflections and dynamics that would occur for a lightweight, downwind 13‐MW rotor. The subscale rotor was mounted and field tested on the two‐bladed Controls Advanced Research Turbine (CART2) at the National Renewable Energy Laboratory's Flatiron Campus (NREL FC). The parked testing of these highly flexible blades included both pitch‐to‐run and pitch‐to‐feather configurations with the blades in the horizontal braked orientation. The collected experimental data includes the unsteady flapwise root bending moments and tip deflections as a function of inflow wind conditions. The bending moments are based on strain gauges located in the root section, whereas the tip deflections are captured by a video camera on the hub of the turbine pointed toward the tip of the blade. The experimental results are compared against computational predictions generated by FAST, a wind turbine simulation software, for the subscale and full‐scale models with consistent unsteady wind fields. FAST reasonably predicted the bending moments and deflections of the experimental data in terms of both the mean and standard deviations. These results demonstrate the efficacy of the first such aeroelastically scaled turbine test and demonstrate that a highly flexible lightweight downwind coned rotor can be designed to withstand extreme loads in parked conditions.

17 WIND ENERGY↗

DOE OSTI2023

Dynamics of Current-Sharing Within a REBCO Tape-Stack Cable

Current redistribution in a non-insulating REBCO tape-stack cable may enable the use of such a cable in a high field dipole without provision for transposition of the constituent tapes. As current is increased in the REBCO cable a dynamic resistance arises in the superconducting layer near operation at critical current. The critical current of each tape in the stack is determined by the local magnetic field variation within the cable. The superconducting-normal transition in REBCO occurs over a working range of current and, within that range, the longitudinal Ohmic electric field produces a transverse electric field between adjacent tapes with different critical currents. Circuit models indicate that, as a given tape approaches critical operation, current will naturally redistribute to neighboring tapes with higher current carrying capacity due to the dynamic rise in resistivity thus preventing premature quench. Here, a multi-scale model is being developed to study the dynamics of current-sharing, the limits of stability, and the impact of fluctuations in critical-current density along each tape.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Cloud response to co-condensation of water and organic vapors over the boreal forest

Abstract. Accounting for the condensation of organic vapors along with water vapor (co-condensation) has been shown in adiabatic cloud parcel model (CPM) simulations to enhance the number of aerosol particles that activate to form cloud droplets. The boreal forest is an important source of biogenic organic vapors, but the role of these vapors in co-condensation has not been systematically investigated. In this work, the environmental conditions under which strong co-condensation-driven cloud droplet number enhancements would be expected over the boreal biome are identified. Recent measurement technology, specifically the Filter Inlet for Gases and AEROsols (FIGAERO) coupled to an iodide-adduct chemical ionization mass spectrometer (I-CIMS), is utilized to construct volatility distributions of the boreal atmospheric organics. Then, a suite of CPM simulations initialized with a comprehensive set of concurrent aerosol observations collected in the boreal forest of Finland during spring 2014 is performed. The degree to which co-condensation impacts droplet formation in the model is shown to be dependent on the initialization of temperature, relative humidity, updraft velocity, aerosol size distribution, organic vapor concentration, and the volatility distribution. The predicted median enhancements in cloud droplet number concentration (CDNC) due to accounting for the co-condensation of water and organics fall on average between 16 % and 22 %. This corresponds to activating particles 10–16 nm smaller in dry diameter that would otherwise remain as interstitial aerosol. The highest CDNC enhancements (ΔCDNC) are predicted in the presence of a nascent ultrafine aerosol mode with a geometric mean diameter of ∼ 40 nm and no clear Hoppel minimum, indicative of pristine environments with a source of ultrafine particles (e.g., via new particle formation processes). Such aerosol size distributions are observed 30 %–40 % of the time in the studied boreal forest environment in spring and fall when new particle formation frequency is the highest. To evaluate the frequencies with which such distributions are experienced by an Earth system model over the whole boreal biome, 5 years of UK Earth System Model (UKESM1) simulations are further used. The frequencies are substantially lower than those observed at the boreal forest measurement site (< 6 % of the time), and the positive values, peaking in spring, are modeled only over Fennoscandia and the western parts of Siberia. Overall, the similarities in the size distributions between observed and modeled (UKESM1) are limited, which would limit the ability of this model, or any model with a similar aerosol representation, to project the climate relevance of co-condensation over the boreal forest. For the critical aerosol size distribution regime, ΔCDNC is shown to be sensitive to the concentrations of semi-volatile and some intermediate-volatility organic compounds (SVOCs and IVOCs), especially when the overall particle surface area is low. The magnitudes of ΔCDNC remain less affected by the more volatile vapors such as formic acid and extremely low- and low-volatility organic compounds (ELVOCs and LVOCs). The reasons for this are that most volatile organic vapors condense inefficiently due to their high volatility below the cloud base, and the concentrations of LVOCs and ELVOCs are too low to gain significant concentrations of soluble mass to reduce the critical supersaturations enough for droplet activation to occur. A reduction in the critical supersaturation caused by organic condensation emerges as the main driver of the modeled ΔCDNC. The results highlight the potential significance of co-condensation in pristine boreal environments close to sources of fresh ultrafine particles. For accurate predictions of co-condensation effects on CDNC, also in larger-scale models, an accurate representation of the aerosol size distribution is critical. Further studies targeted at finding observational evidence and constraints for co-condensation in the field are encouraged.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Predictive Modeling of Electrodeposition in a Single Pore Flow-Through Electrode: From Electronucleation to Coating Thickness Uniformity

The emergence of advanced manufacturing methods capable of producing porous three-dimensional structures has expanded the design space for next-generation functional components. The ability to fabricate ordered 3D foams for use in electrocatalysis reactors has increased the need for controlled deposition of catalytic metals onto porous support materials, such as carbon. However, there is a lack of clear design guidelines for electrodeposition onto 3D substrates, due to the geometric complexity and multi-scale nature of the problem. Furthermore, electro-nucleation phenomena are often overlooked in macro-scale models of current distribution during deposition. Here, a graphite flow-through electrode (FTE) is used as a model system for copper deposition within a single pore. Potential distributions and electro-nucleation phenomena are coupled in a continuum level model by incorporating nucleation size and density as a function of overpotential, determined experimentally using in situ atomic force microscopy (AFM). The model predictions are validated by measuring the coating uniformity in the pore using micro-computed X-ray tomography (μCT). A scaling analysis comprising dimensionless parameters such as the Wagner number is presented. The simplified scaling relationship framework can guide the electrodeposition process and electrode design to optimize plating of porous substrates under fluid flow conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Advanced PGM-free Cathode Engineering for High Power Density and Durability

Polymer electrolyte fuel cells (PEFCs) are among the most promising technologies for future electric vehicles by using clean H2 with much-improved energy conversion efficiency, longer range, and rapid refueling. However, due to a large amount of platinum group metal (PGM) catalyst used in their electrodes, their prohibitively high cost hinders broad commercialization of PEFCs for transportation. Therefore, there is a critical need to develop low-cost, high-performance PGM-free cathode catalysts that have the potential to dramatically transform the economics of PEFC commercialization by reducing catalyst costs by one to two orders of magnitude. However, before PGM-free cathodes become viable, several technical challenges associated with PGM-free cathodes must be addressed, including insufficient activity and stability of the catalysts, as well as severe water flooding and large transport losses in the electrodes. Overcoming those barriers and ultimately meeting the challenging automotive PEFC performance targets was the focus of this comprehensive research and development effort on new PGM-free cathodes. To this end, we assembled a team including leading researchers from universities and industry with different but complementary expertise and capabilities. The project combined three novel and promising approaches: Advanced metal-organic framework (MOF)-derived M-N-C catalysts with a high activity and impressive durability, Novel PGM-free specific cathode architectures and fabrication strategies capable of addressing the substantial flooding and transport resistances in thicker cathodes by introducing engineered hydrophobicity through additives and support layers, and Advanced electrode ionomers with high proton conductivity for low ohmic losses across the electrode and more uniform catalyst utilization. The implementation of these new materials and electrode designs was supported by a suite of advanced experimental and simulation tools that allows us to identify performance and durability bottlenecks, devise solutions, and establish rational material design and synthesis targets. These methods include advanced electrochemical characterization, high-resolution imaging, and multi-scale modeling. In addition, the project team leveraged a broad cross-section of the ElectroCat consortium’s national laboratory facilities and expertise in advancing these materials and design strategies. Finally, the industry partners on the project facilitated the evaluation of scaled-up synthesis and manufacturing in the United States. Over its four-year period, the project made significant year-over-year advances in PGM-free cathode performance and viability. A combination of high activity and highly durable catalysts were developed through novel catalyst synthesis strategies, which met several performance and durability targets. More specifically, a catalyst prepared from MOFs and Fe2O3 nanoparticles with ammonium chloride and chemical vapor deposition treatments yielded a significant advancement in PGM-free cathode durability. Several novel strategies for fabricating cathodes were demonstrated, including those designed to reduce flooding and thickness of the cells for significantly increased volumetric power density. An optimized cathode with high conductivity ionomer and tuned ink processing for hydrophobicity yielded high fuel cell performance with new levels power density and maximum current. The scientific studies and modeling assessment also provided an outlook for future efforts, including a focus on catalysts with an increased density of the highly stable active sites developed in this project.

08 HYDROGEN↗

DOE OSTI2022

Cavitation in a soft porous material

Abstract We study the collapse and expansion of a cavitation bubble in a deformable porous medium. We develop a continuum-scale model that couples compressible fluid flow in the pore network with the elastic response of a solid skeleton. Under the assumption of spherical symmetry, our model can be reduced to an ordinary differential equation that extends the Rayleigh–Plesset equation to bubbles in soft porous media. The extended Rayleigh–Plesset equation reveals that finite-size effects lead to the breakdown of the universal scaling relation between bubble radius and time that holds in the infinite-size limit. Our data indicate that the deformability of the porous medium slows down the collapse and expansion processes, a result with important consequences for wide-ranging phenomena, from drug delivery to spore dispersion.

Leng, Yu↗

DOE OSTI2024

The role of ion-scale micro-turbulence in pedestal width of the DIII-D wide-pedestal QH mode

The low-edge rotation, intrinsically ELM-free, and improved confinement wide-pedestal quiescent H-mode (QH-mode), discovered in DIII-D tokamak, has pedestal widths exceeding the EPED-kinetic-ballooning mode (KBM) model scaling typically by at least 25%. Ion-scale ($k_yp_s$ < 1) microturbulence and its role in setting the pedestal structure is investigated using the radially local δ$f$ gyrokinetic code CGYRO. The electromagnetic trapped electron mode (TEM) is unstable at the pedestal top, while plasma beta (β$_e$) is ∼60% below the KBM onset threshold and the electron temperature gradient mode is found to be unstable in the peak gradient region. Nonlinear simulation reveals that the ion-scale turbulence could produce electron energy flux consistent with the flux inferred from power balance at the pedestal top, with a reasonable variation of the local shearing rate; and the local neoclassical transport from NEO is dominant over the simulated turbulent transport in the ion energy flux channel. The simulated ion-scale turbulence produces much lower electron energy flux than inferred from experiment in the pedestal peak gradient region. A correction to the EPED-KBM pedestal width scaling is obtained based on the two-dimensional scan of pedestal top plasma beta (β$_e$) and normalized electron density and temperature scale lengths,${a}$/${L_n}_e$, ${a}$/${L_T}_n$using CGYRO linear simulations. Mode transitions among TEM, micro-tearing mode, ion-temperature gradient mode and KBM, are observed in the 2D scan at the pedestal top. A fixed normalized growth rate for these drift-type modes is taken to determine the pedestal width scaling, which shows good consistency with the QH experimental database on pedestal heights and widths. The onset of KBM instabilities and the local E × B shear suppression criterion set the lower and upper limit for the pedestal width of standard QH-mode, wide-pedestal QH-mode and type-I ELMy H mode. A potentially higher and wider pedestal is expected from the new scaling of pedestal width. This work presents an improved understanding of the ion-scale micro-turbulence of wide-pedestal QH-mode and sheds light on a promising scenario for future reactors, including ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Membrane Degradation in PEM Fuel Cells: Part I. Modeling Gas Crossover and the Pt Band

Understanding chemical degradation of the proton-exchange membrane in fuel cells is crucial for extending their lifetimes. Herein, various degradation reactions reported in literature are organized and analyzed, including direct radical generation and an indirect (Fenton) pathway. To understand the transport of dissolved H 2 and O 2 crossover gases as they relate to membrane degradation, an agglomerate-scale model is introduced, treating gas, ionomer, and catalyst as discrete phases. The model reveals a key phenomenon: at working potentials, dissolved gases are mostly consumed at the interface between the catalyst layer and the membrane, leaving little gas to cross the membrane. Under open-circuit conditions, dissolved gases are not consumed and can then cross the membrane. This explains high H 2 O 2 concentrations and degradation rates seen in experiments but not captured in previous models. Following mixed-potential theory, crossover gases supply the hydrogen-oxidation and oxygen-reduction reactions, which occur simultaneously on individual Pt particles comprising the Pt band in the membrane, forming reactive species (H 2 O 2 , OH·). Results show crossover gas almost entirely reacts on the Pt band, allowing little to reach the opposite electrode. Furthermore, the micro-scale geometry of the catalyst-layer/membrane interface impacts the gas crossover at working potentials, indicating that cell construction affects membrane durability.

Johnson, Evan F. [Lawrence Berkeley National Labor↗

DOE OSTI2024

Scaling neural simulations in STACS

Abstract As modern neuroscience tools acquire more details about the brain, the need to move towards biological-scale neural simulations continues to grow. However, effective simulations at scale remain a challenge. Beyond just the tooling required to enable parallel execution, there is also the unique structure of the synaptic interconnectivity, which is globally sparse but has relatively high connection density and non-local interactions per neuron. There are also various practicalities to consider in high performance computing applications, such as the need for serializing neural networks to support potentially long-running simulations that require checkpoint-restart. Although acceleration on neuromorphic hardware is also a possibility, development in this space can be difficult as hardware support tends to vary between platforms and software support for larger scale models also tends to be limited. In this paper, we focus our attention on Simulation Tool for Asynchronous Cortical Streams (STACS), a spiking neural network simulator that leverages the Charm++ parallel programming framework, with the goal of supporting biological-scale simulations as well as interoperability between platforms. Central to these goals is the implementation of scalable data structures suitable for efficiently distributing a network across parallel partitions. Here, we discuss a straightforward extension of a parallel data format with a history of use in graph partitioners, which also serves as a portable intermediate representation for different neuromorphic backends. We perform scaling studies on the Summit supercomputer, examining the capabilities of STACS in terms of network build and storage, partitioning, and execution. We highlight how a suitably partitioned, spatially dependent synaptic structure introduces a communication workload well-suited to the multicast communication supported by Charm++. We evaluate the strong and weak scaling behavior for networks on the order of millions of neurons and billions of synapses, and show that STACS achieves competitive levels of parallel efficiency.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Screening green solvents for multilayer plastic film recycling processes

Multilayer (ML) plastic films are essential packaging materials that help protect products from diverse external factors; however, only 5% of all ML films are recycled in the United States. Solvent-based technologies are a promising alternative for recycling ML films because they enable recovery of constituent polymer resins. For example, the Solvent Targeted Recovery and Precipitation (STRAPTM) process sequentially dissolves and separates polymer components using a series of targeted solvent washes. A crucial design aspect of this process is the impact of selected solvents on human health and on the environment. Here, this work introduces a computational framework that integrates molecular modeling, process modeling, techno-economic analysis (TEA), and life-cycle analysis (LCA) to quickly screen green solvents for solvent-based ML recycling processes. Initial screening for solvents based on selectivity is performed by estimating temperature-dependent solubilities using molecular-scale models. Subsequent screening uses basic estimates of energy use and octanol-water partition coefficients (logP) as key measures of health, safety, and environmental hazards. Detailed process modeling, TEA, and LCA are used on a reduced set of promising solvents identified in early screening steps to more accurately determine how solvent selection and associated operating conditions impact overall economics and environmental impacts. The framework is used for the identification of green solvents (from a database of 1,000 solvents) that separate an industrial ML film composed of polyethylene (PE), ethylene vinyl alcohol (EVOH), and polyethylene terephthalate (PET). Our analysis shows the effectiveness of the framework and reveals fundamental trade-offs between solvent greenness, solubility, and economics. Our work emphasizes the importance of taking a holistic systems view during solvent design and aims to inform the development of new processes for ML film recycling and the identification of new ML films that are easier to recycle.

economics↗