Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “process loads”

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 415 records · Page 23

Cellulose molecular weight of corn stover with cotreatment

Here we present a dataset of cellulose molecular weights from corn stover after a cotreatment process. The corn stover sample was subjected to a combination of mechanical disruption (ball milling) and consolidated bioprocessing (CBP) process using the bacterium C. thermocellum at 60 grams/L solids loadings. Cellulose was isolated from control, no cotreatment control, and cotreatment corn stover. Cellulose was then derivatized using anhydrous pyridine and phenyl isocyanate over 48 h at 70 ºC. The reaction was quenched by anhydrous methanol. Methanol and water mixture (7:3, v/v) was added dropwise to promote precipitation of the cellulose derivative. The solids were collected by filtration and thoroughly washed with the methanol and water mixture followed by water. The cellulose derivative was dried overnight under vacuum at 40°C. The obtained cellulose tricarbanilates was then dissolved in tetrahydrofuran (THF) at 1.0 mg/mL and the solution was filtered through a 0.45 µm PTFE filter and placed in an auto-sampler vial. The molecular weight of cellulose was measured by Agilent GPC SECurity 1200 system equipped with four Waters Styragel columns (i.e., HR0.5, HR2, HR4, and HR6) and a UV detector (270 nm). THF was used as the mobile phase. Data collection and processing was performed by Polymer Standards Service WinGPC Unity software (Build 6807). The molecular weight was calculated relative to the polystyrene calibration curve. The data on molecular weights provided information about the effects of cotreatment on corn stover cellulose molecular changes.

Cellulose, molecular weights, corn stover, cotreat↗

Identifying Bounds of Inorganic Content in Solventless Processing of Hybrid Solid Electrolytes

Solid-state lithium batteries require safe, robust electrolytes to enable higher energy densities and improved safety over conventional cells. Hybrid polymer–ceramic electrolytes are a promising solution, combining the processability of polymers with the high ionic conductivity and mechanical strength of inorganic fillers. In this work, we demonstrate a solventless, UV-curing method to produce hybrid solid electrolytes using a poly(ethylene glycol) dimethyl ether (PEGDME)-based photocurable matrix incorporating Li 1.5 Al 0.5 Ge 1.5 (PO 4 ) 3 (LAGP) or Li 7 La 3 Zr 2 O 12 (LLZO) ceramic electrolyte. Inorganic filler loadings up to ∼55 wt.% could be successfully incorporated via this process which was the highest inorganic content at which the slurry remains processable and cured into a uniform film. The resulting UV-cured composite electrolytes remain flexible and exhibit room-temperature ionic conductivities on the order of 10 −4 S·cm −1 , along with notably improved lithium-ion transference numbers compared to conventional polymer electrolytes. Similar performance and processing limits were observed for both LAGP and LLZO, indicating that ceramic filler chemistry does not significantly affect the UV-curing process or the electrolyte's ion transport properties in this regime. Eliminating solvents from fabrication not only simplifies processing and mitigates environmental concerns but also enables higher solid contents that enhance mechanical strength and help suppress lithium dendrite formation. In conclusion, this scalable approach thus paves the way for manufacturing robust composite solid electrolytes for next-generation solid-state batteries (SSBs).

Batteries↗

Full-Field Strain Measurement Integrated with Two Dimension Regression Analysis to Evaluate the Bi-Modulus Elastic Properties of Isotropic and Transversely Isotropic Materials

Background: Measuring the physical properties of shale is critical for optimizing engineering activities such as geothermal energy generation and hydraulic fracturing. Shale is a transversely isotropic material. Furthermore, this material can also include micro and macro cracks at different locations and orientations that cause it to behave differently under tensile or compressive loading. Objective: In this work, a combined experimental–numerical approach is proposed to evaluate the bi-modulus elastic properties of isotropic and transversely isotropic materials. Methods: Full-field strain measurements for a circular disk under diametral compression are integrated with a regression analysis technique to evaluate the elastic properties of bi-modulus materials subjected to tensile and compressive loads using two loading configurations on the same specimen. Digital Image Correlation (DIC) is used to measure the full-field strains. Subsequently, in the case of an isotropic material, a linear least-squares approach is utilized to process the experimentally determined strains in conjunction with analytical expressions of the stress fields (in terms of far-field loading) to determine the elastic modulus E, the shear modulus G, and the Poisson’s ratio $v$. In the case of a transversely isotopic material, such as shale, a finite element model is implemented to determine the stress fields (again in terms of far-field loading), which is followed by repeating the previous regression analysis in an iterative process to estimate the elastic parameters. Results: The results show that the proposed technique successfully provides a complete set of elastic properties as a function of both the loading condition and the principal material directions. The technique is validated by measurements on a known isotropic material and then applied to determine the properties of shale. Conclusion: In this work, the proposed approach is successfully used to calculate the bi-modulus elastic response of poly(methyl meth- acrylate) (PMMA) and shale. As expected, PMMA exhibits an isotropic response with no bi-modulus effect, however, shale exhibits both transverse isotropy and a bi-modulus effect. Therefore, this approach holds promise for investigating the elastic properties of materials like rocks and fiber-reinforced composite laminates as functions of the principal material directions and the loading conditions.

42 ENGINEERING↗

Impacts of constraints and uncertainties on projected amount of Hanford low-activity waste glasses

A systematic evaluation of Hanford low-activity waste (LAW) loading in alkali-aluminoborosilicate glasses was performed. A set of 22 waste compositions were selected to represent the range of wastes anticipated across the Hanford mission. Glass composition was formulated for each waste to simultaneously satisfy a full cadre of property and composition constraints while maximizing waste loading. Glass properties were estimated using property models recently developed for Hanford LAW Facility operation. The impacts from the property and composition constraints on the glass mass to be produced were calculated and compared. Electrical conductivity, Vapor Hydration Test and SO 3 melter tolerance, and canister-centerline cooling crystal constraints have a large effect on reducing total glass mass, while the impact from Product Consistency Test, viscosity, and isothermal crystal constraints is relatively small. Additionally, effects of model prediction and process/composition uncertainties were quantified. The results identified areas of future research that could be useful to improve the efficiency of plant operation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Small to Full Height Scale Comparisons of Cesium Ion Exchange Performance with Crystalline Silicotitanate

The U.S. Department of Energy’s (DOE) Hanford Site houses 56 million gallons of high-level radioactive waste generated from plutonium production from 1944 to 1988. The supernatant waste, currently stored in underground tanks, is intended to be vitrified following filtration and 137 Cs removal at the Hanford Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Vitrification Facility. The WTP Pretreatment Facility will not be operational for several years. The Tank Side Cesium Removal (TSCR) system is a technology demonstration that will remove cesium from tank waste supernate to support directly feeding LAW to the vitrification facility. The 137 Cs removal is important to meet the WTP LAW contract specification and ultimately for creating a contact-handled waste form. The waste acceptance criteria (WAC) limit for the WTP LAW Facility is <3.18E-5 Ci 137 Cs per mole Na. The TSCR system is skid mounted and employs two key technologies: dead end filtration for solids removal and ion exchange (IX) for cesium removal. Filtration is necessary to protect the functionality of the IX columns. The IX process utilizes sodium form crystalline silicotitanate (CST) IX media produced by UOP Honeywell, LLC (Des Plaines, IL) under the product name IONSIV R9140-B, 18 x 50 mesh, in a lead-lag-polish column configuration. Each column contains a CST bed height of 92 inches and a volume of approximately 157 gallons of CST IX media. The full-scale TSCR operation will run at a residence time of 1.9 bed volumes (BV)/h, which results in a superficial velocity of 7.3 cm/min and a flowrate of ~5 gal/min. Column testing at small (2.5% of the full bed height), medium (12% of the full bed height), and full-height scales has been previously conducted to evaluate process variables and scale up performance of Cs exchange onto the CST. Cesium load performances from various sieve cuts at the small scale indicated that a 30-mesh sieve cut be tested to determine if it better reflects the 12% and full-height column performances at the small scale. Two process flowrates were tested in the small-column configuration with <30-mesh CST and simulant solutions. These tests were compared to full-height column tests at the same residence times to assess CST particle size effect on column scaling. Table S.1 summarizes the observed column performance determined for the two flowrates juxtaposed to the previous work with 5.6 M Na simulant at higher scales (used as benchmarks). The WAC breakthroughs between the small and full-height tests at 1.3 BV/h differed by ~59 BVs. The WAC breakthroughs for small, medium, and full-height tests at 1.8 BV/h were consistent at 240 BVs. The 50% Cs breakthroughs were nominally equivalent for all column tests. The common transition zones and onsets of Cs breakthrough at the 1.8 BV/h tests indicated that the Cs mass transfers were equivalent and thus the <30-mesh CST Cs load performance at the small scale successfully modeled that of the full-height system. It is recommended that the <30-mesh CST be used in subsequent 10-mL CST bed tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Forming Complex Nuclear Fuel Shapes in High-Loaded Silicide Surrogates

This work provides proof of the concept that high silicide loading nuclear fuel meat surrogates with complicated geometries can be produced with uniform density through an application of cold isostatic pressing (CIP). Dispersion fuels with high volumetric loading of U_3 Si_2 have challenges in fabrication. Fabrication involves a series of processes including powder compaction, rolling pack assembly, and roll forming. PNNL has undertaken extensive experimental work using MoSi2 and WSi2 as surrogates for U3Si2 to explore the feasibility of reducing or eliminating the issues through the application of CIP in the powder pressing step prior to rolling pack assembly. The composites were prepared at >40 vol% silicide loading, which was representative of 4.8 gU/cm^3 and formed via CIP at 50 KSI (Kilopound per Square Inch) pressure. The CIP mold design was taken through a design process which aimed to reduce defects and increase precision. The application of CIP here provides a broad strategy for producing highly loaded dispersion fuels with complex geometries and uniform density.

Clelland, Dustin T.↗

Automated Lane Centering: An Off-the-Shelf Computer Vision Product vs. Infrastructure-Based Chip-Enabled Raised Pavement Markers

Safe autonomous vehicle (AV) operations depend on an accurate perception of the driving environment, which necessitates the use of a variety of sensors. Computational algorithms must then process all of this sensor data, which typically results in a high on-vehicle computational load. For example, existing lane markings are designed for human drivers, can fade over time, and can be contradictory in construction zones, which require specialized sensing and computational processing in an AV. But, this standard process can be avoided if the lane information is simply transmitted directly to the AV. High definition maps and road side units (RSUs) can be used for direct data transmission to the AV, but can be prohibitively expensive to establish and maintain. Additionally, to ensure robust and safe AV operations, more redundancy is beneficial. A cost-effective and passive solution is essential to address this need effectively. In this research, we propose a new infrastructure information source (IIS), chip-enabled raised pavement markers (CERPMs), which provide environmental data to the AV while also decreasing the AV compute load and the associated increase in vehicle energy use. CERPMs are installed in place of traditional ubiquitous raised pavement markers along road lane lines to transmit geospatial information along with the speed limit using long range wide area network (LoRaWAN) protocol directly to nearby vehicles. This information is then compared to the Mobileye commercial off-the-shelf traditional system that uses computer vision processing of lane markings. Our perception subsystem processes the raw data from both CEPRMs and Mobileye to generate a viable path required for a lane centering (LC) application. To evaluate the detection performance of both systems, we consider three test routes with varying conditions. Our results show that the Mobileye system failed to detect lane markings when the road curvature exceeded ±0.016 m -1 . For the steep curvature test scenario, it could only detect lane markings on both sides of the road for just 6.7% of the given test route. On the other hand, the CERPMs transmit the programmed geospatial information to the perception subsystem on the vehicle to generate a reference trajectory required for vehicle control. The CERPMs successfully generated the reference trajectory for vehicle control in all test scenarios. Moreover, the CERPMs can be detected up to 340 m from the vehicle’s position. Our overall conclusion is that CERPM technology is viable and that it has the potential to address the operational robustness and energy efficiency concerns plaguing the current generation of AVs.

33 ADVANCED PROPULSION SYSTEMS↗

An iterative bidirectional gradient boosting approach for CVR baseline estimation

Here this paper presents a novel Iterative Bidirectional Gradient Boosting Model (IBi-GBM) for estimating the baseline of Conservation Voltage Reduction (CVR) programs. In contrast to many existing methods, we treat CVR baseline estimation as a missing data retrieval problem. The approach involves dividing the load and its corresponding temperature profiles into three periods: pre-CVR, CVR, and post-CVR. To restore the missing load profile during the CVR period, the method employs a three-step process. First, a forward-pass GBM is executed using data from the pre-CVR period as inputs. Subsequently, a backward-pass GBM is applied using data from the post-CVR period. The two restored load profiles are reconciled, considering pre-calculated weights derived from forecasting accuracy, and only the leftmost and rightmost points are retained. The newly restored points are then included as inputs for the subsequent iteration. This iterative procedure continues until the original load data in the CVR period is fully restored. We develop IBi-GBM using actual smart meter and Supervisory Control and Data Acquisition (SCADA) data. Our results demonstrate that IBi-GBM exhibits robust performance across various data resolutions and in different seasons and outperforms existing methods by achieving a 1-2% reduction in normalized Root Mean Square Error (nRMSE).

42 ENGINEERING↗

Optimizing high energy density sulfur cathodes: A multivariate approach to electrode formulation and processing

Lithium-sulfur (Li-S) batteries involve complex solid-liquid-solid phase transformations during both discharging and charging processes, where cathode materials, formulation, and structure play a crucial role. Here, a design of experiments (DoE) methodology and an empirical model are developed to systematically explore the interactions and trade-offs among cathode factors and process variables, and to obtain generalizable effects estimates for the multivariate system. Compared to the conventional one-factor-at-a-time (OFAT) approach, this work demonstrates advantages in both efficiency and accuracy by allowing the data to guide future research and decisions. Further, an optimized cathode formulation and processing parameters are predicted and validated experimentally, achieving over 1000 mAh g -1 in discharge capacity and improved cycling under practical lean electrolyte (4 µL mg -1 S) and high S-loading cathodes (>4 mg cm -2 ) conditions. The optimized cathode was scaled up and assembled into Li-S pouch cells, achieving 316 Wh kg -1 in cell-level energy, proving that the comprehensive and rigorous framework for optimizing complex systems with DoE leads to improved performance in a practical pouch cell system.

25 ENERGY STORAGE↗

Coupled heat and mass transfer in internally cooled silica gel-coated aluminum foam heat exchangers for dehumidification

Solid desiccant dehumidifiers offer significant potential to improve humidity control in buildings due to their ability to independently regulate latent (humidity) and sensible (temperature) loads. However, adsorption of water vapor by desiccants is an exothermic process; the resulting heat of adsorption increases both the desiccant bed temperature and the temperature of the air flowing over the bed, thereby limiting the performance of desiccant dehumidifiers. To address this problem, desiccant-coated heat exchanger prototypes were developed by coating silica gel onto aluminum foam substrates and integrated with chilled water circulation for effective removal of the heat of adsorption. Prototype heat exchangers with two different pore densities (10 and 20 pores per inch (PPI)) and two different desiccant loadings (∼10.36% and ∼17.80%) were fabricated. Experiments were conducted by varying the inlet air face velocity, chilled water inlet temperature, and chilled water flow rate to investigate the coupled heat and mass transfer characteristics of the prototype heat exchangers over a wide range of operating conditions. The underlying transport mechanisms were analyzed using transient desiccant bed temperature, outlet air temperature, and outlet humidity profiles. The results demonstrate that internal cooling significantly enhances dehumidification performance by maintaining more favorable thermodynamic conditions for water vapor adsorption. Under identical operating conditions, a 20 PPI heat exchanger coated with 25.05 g of silica gel achieved 98.98% higher cumulative water uptake and 97.75% higher moisture removal capacity than the corresponding non-cooled case after 5 min of adsorption. This substantial enhancement was further supported by a validated one-dimensional coupled heat and mass transfer model developed in the present study. To date, this is one of the first experimental studies to systematically investigate the combined effects of internal cooling, foam pore density, and desiccant loading on the performance of silica-gel-coated metal foam heat exchangers under HVAC-relevant operating conditions.

42 ENGINEERING↗

Fast and accurate reduced-order modeling of a MOOSE-based additive manufacturing model with operator learning

One predominant challenge in additive manufacturing (AM) is to achieve specific material properties by manipulating manufacturing process parameters during the runtime. Such manipulation tends to increase the computational load imposed on existing simulation tools employed in AM. The goal of the present work is to construct a fast and accurate reduced-order model (ROM) for an AM model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, ultimately reducing the time/cost of AM control and optimization processes. Our adoption of the operator learning (OL) approach enabled us to learn a family of differential equations produced by altering process variables in the laser’s Gaussian point heat source. More specifically, we used the Fourier neural operator (FNO) and deep operator network (DeepONet) to develop ROMs for time-dependent responses. Furthermore, we benchmarked the performance of these OL methods against a conventional deep neural network (DNN)-based ROM. Ultimately, we found that OL methods offer comparable performance and, in terms of accuracy and generalizability, even outperform DNN at predicting scalar model responses. The DNN-based ROM afforded the fastest training time. Furthermore, all the ROMs were faster than the original MOOSE model yet still provided accurate predictions. FNO had a smaller mean prediction error than DeepONet, with a larger variance for time-dependent responses. Unlike DNN, both FNO and DeepONet were able to simulate time series data without the need for dimensionality reduction techniques. Finally, the present work can help facilitate the AM optimization process by enabling faster execution of simulation tools while still preserving evaluation accuracy.

36 MATERIALS SCIENCE↗

Final Report DE-FE0031785 Mohsen Ahmadian, Ph.D. Demonstration of Proof of Concept of a Multiphysics Approach for Real-Time Remote Monitoring of Dynamic Changes in Pressure and Salinity in Hydraulically Fractured Networks

Hydraulic fracturing has evolved into a multistep process with varying flow rates, carrier fluids (e.g., gel or slickwater), proppant loadings, and proppant grain sizes. As a result, primary recovery from a hydraulically fractured tight-oil reservoir is often a tiny fraction of the original oil in place, ranging between 5 and 10%. As stated in the FOA1990, “part of this problem is due to the inability of current well completion processes to effectively stimulate the entire reservoir volume in contact with the wellbore. Innovative technologies are needed that can help improve the effectiveness of reservoir completion methods, maximize stimulated reservoir volumes, and optimize recovery over the entire producing life span of a well”. We first need to enhance the current fracture diagnostic techniques to improve a well-completion design. However, detecting and delineating a subsurface hydraulic fracture is extremely difficult because the induced fracture network is only fractionally propped, and these propped fractures are generally very thin. Microseismic and tiltmeter monitoring techniques can provide information on the fracture extent but provide little or no information on the movement and final distribution of proppant or production fluids. On the other hand, electromagnetic (EM) imaging has shown the capability to monitor proppant distribution throughout the fracture area, especially in the presence of Electrically Active Proppants (EAPs). A previous EM survey of hydraulic fracturing at the Devine Fracture Pilot Site (DFPS) and subsequent EM code developments demonstrated this survey as a robust technique to remotely interrogate the extent of the EAP-filled hydraulic fracture during its propagation. The objectives of the project were threefold: (1) to capitalize on the material properties of an EAP to demonstrate remote monitoring of relative changes in pressure, pressure, and flow that are commonly encountered during production from a hydraulically fractured reservoir; (2) to evaluate EM imaging tools, to achieve Objective 1 in near real-time; and (3) to develop a multi-physics joint inversion approach to precisely predict flow patterns and physiochemical changes within an EAP-filled fracture network. This research project was built upon our previous work at the Devine Test Site managed by the Bureau of Economic Geology (BEG) at The University of Texas at Austin (UT-Austin). It also leveraged a significant investment from the Advanced Energy Consortium (AEC) to address the DOE's interest in subsurface flow, containment, and characterization by multiple signals. This three-year and three-month project succeeded in demonstrating the feasibility of a real-time dynamic fluid flow mapping technique at Technology Readiness Level 5 (TRL5) by utilizing a commercially available surface-based Controlled-Source Electromagnetic (CSEM) method (Objectives 1, 2). We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, the observed electric field in our study is influenced by bottomhole pressure, flow rate, and salinity, which is demonstrated by comparing these parameters with the electrical field potentials. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce the whole measured electric field magnitudes. Preliminary estimates showed that including Streaming Potential (SP) in our geophysical model is likely needed to reduce the simulation misfit.

02 PETROLEUM↗

High-Capacity Enzymatic Degradation of Postconsumer Poly(ethylene terephthalate): Modeling and Experimental Investigations

Enzymatic degradation is considered as one of the key steps for biorecycling of polyethylene terephthalate (PET), a widely used plastic. Recently, we have developed an optimized leaf-branch compost cutinase PelB-LCC ICCG and demonstrated ∼80% degradation of 200 g L −1 recycled PET (RPET) within two days in bioreactors. In this study, continued research efforts were made to achieve a complete degradation of up to 300 g L −1 RPET in bioreactors. First, low-capacity reaction conditions with ≤5 g L −1 PET were examined to identify the key variables for PET degradation with PelB-LCC ICCG , including temperature, enzyme loading, and product inhibition. To better understand and further optimize the enzymatic degradation process, a mechanism-based model was established to describe the kinetics of PET degradation and formation of the main product terephthalic acid (TPA) and byproducts MHET and BHET in high-capacity bioreactors. Model simulation suggested that a minimal enzyme loading of ∼1.4 mg PelB-LCC ICCG g −1 RPET is required to achieve a nearly complete degradation of RPET within 48 h, which was used to guide more high-capacity experiments with 100−300 g L −1 RPET in fully controlled 1 L bioreactors. A higher temperature (≥65 °C) was found not only to enable rapid degradation in the beginning but also to induce a gradual increase in RPET’s crystallinity and significantly slow down the degradation after 48 h. A high loading of RPET solids and the accumulation of the produced insoluble TPA pose a big challenge on mixing and mass transfer in the stirred bioreactor, which can be addressed by increasing the stirring speed. The results pave the way toward biorecycling of PET at a large scale.

enzymatic degradation↗

Ammonia-salt solvent promotes cellulosic biomass deconstruction under ambient pretreatment conditions to enable rapid soluble sugar production at ultra-low enzyme loadings

In this paper, we report a novel ammonia : ammonium salt solvent based pretreatment process that can rapidly dissolve crystalline cellulose into solution and eventually produce highly amorphous cellulose under near-ambient conditions. Pre-activating the cellulose I allomorph to its ammonia–cellulose swollen complex (or cellulose III allomorph) at ambient temperatures facilitated rapid dissolution of the pre-activated cellulose in the ammonia-salt solvent ( i.e. , ammonium thiocyanate salt dissolved in liquid ammonia) at ambient pressures. For the first time in reported literature, we used time-resolved in situ neutron scattering methods to characterize the cellulose polymorphs structural modification and understand the mechanism of crystalline cellulose dissolution into a ‘molecular’ solution in real-time using ammonia-salt solvents. We also used molecular dynamics simulations to provide insight into solvent interactions that non-covalently disrupted the cellulose hydrogen-bonding network and understand how such solvents are able to rapidly and fully dissolve pre-activated cellulose III. Importantly, the regenerated amorphous cellulose recovered after pretreatment was shown to require nearly ~50-fold lesser cellulolytic enzyme usage compared to native crystalline cellulose I allomorph for achieving near-complete hydrolytic conversion into soluble sugars. Lastly, we provide proof-of-concept results to further showcase how such ammonia-salt solvents can pretreat and fractionate lignocellulosic biomass like corn stover under ambient processing conditions, while selectively co-extracting ~80–85% of total lignin, to produce a highly digestible polysaccharide-enriched feedstock for biorefinery applications. Unlike conventional ammonia-based pretreatment processes ( e.g. , Ammonia Fiber Expansion or Extractive Ammonia pretreatments), the proposed ammonia-salt process can operate at near-ambient conditions to greatly reduce the pressure/temperature severity necessary for conducting effective ammonia-based pretreatments on lignocellulose.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Distributed Multi-GPU Community Detection on Exascale Computing Platforms

Community detection is a fundamental operation in graph mining, and by uncovering hidden structures and patterns within complex systems it helps solve fundamental problems pertaining to social networks, such as information diffusion, epidemics, and recommender systems. Scaling graph algorithms for massive networks becomes challenging on modern distributed-memory multi-GPU (Graphics Processing Unit) systems due to limitations such as irregular memory access patterns, load imbalances, higher communication-computation ratios, and cross-platform support. We present a novel algorithm HiPDPL-GPU (distributed parallel Louvain) to address these challenges. We conduct experiments involving different partitioning techniques to achieve optimized performance of HiPDPL-GPU on the two largest supercomputers: Frontier and Summit. Remarkably, HiPDPL-GPU processes a graph with 4.2 billion edges in less than 3 minutes using 1024 GPUs. Qualitatively performance of HiPDPL-GPU is similar or better compared to other state-of-the-art CPU- and GPU-based implementations. While prior GPU implementations have predominantly employed CUDA, our first-of-its-kind implementation for community detection is cross-platform, accommodating both AMD and NVIDIA GPUs.

graph algorithms, high performance comptuing↗

Distributed Multi-GPU Community Detection on Exascale Computing Platforms

Community detection is a fundamental operation in graph mining, and by uncovering hidden structures and patterns within complex systems it helps solve fundamental problems pertaining to social networks, such as information diffusion, epidemics, and recommender systems. Scaling graph algorithms for massive networks becomes challenging on modern distributed-memory multi-GPU (Graphics Processing Unit) systems due to limitations such as irregular memory access patterns, load imbalances, higher communication-computation ratios, and cross-platform support. We present a novel algorithm HiPDPL-GPU (Distributed Parallel Louvain) to address these challenges. We conduct experiments involving different partitioning techniques to achieve an optimized performance of HiPDPL-GPU on the two largest supercomputers: Frontier and Summit. Remarkably, HiPDPL-GPU processes a graph with 4.2 billion edges in less than 3 minutes using 1024 GPUs. Qualitatively, the performance of HiPDPL-GPU is similar or better compared to other state-of-the-art CPU- and GPU-based implementations. While prior GPU implementations have predominantly employed CUDA, our first-of-its-kind implementation for community detection is cross-platform, accommodating both AMD and NVIDIA GPUs.

Sattar, Naw Safrin↗

Lignin content of Populus trichocarpa residues after CELF pretreatment and CBP fermentation

Here we present a dataset of lignin content from a woody energy crop (Populus trichocarpa) residues after a series of co-solvent enhanced lignocellulosic fractionation (CELF) pretreatment and consolidated bioprocessing (CBP) process. The natural poplar variant GW-9947 from the Center for Bioenergy Innovation (CBI) was used. The poplar was knife milled and passed through a 1 mm sieve and CELF pretreatment was performed in a Parr autoclave reactor with 7.5 wt % solids loading, 0.5 wt% H2SO4 as catalyst at 150°C with 5, 15, 25 and 30 minutes, respectively. Tetrahydrofuran was added in a 1:1 mass ratio with water as the pretreatment solvent. The residues from CELF pretreatment were then subjected to CBP using the bacterium C. thermocellum DSM 1313. CBP fermentations were performed at 60 °C in a shaker at 50, 75, and 100 grams/L solids loadings. The Klason lignin was measured using a two-step acid hydrolysis process. In brief, the poplar samples were first hydrolyzed by 72 wt% sulfuric acid at 30 oC for an hour. The hydrolysates were then diluted to 4 wt% sulfuric acid using deionized water and subsequently autoclaved at 121 oC for 1 h. Upon the completion of the two-step hydrolysis, the resulting solution was cooled to room temperature and the precipitate was then filtered through a G8 glass fiber filter through a crucible, dried at oven for overnight, and weighed to get the Klason lignin content. The lignin data provides information about lignin content changes after CELF and C. thermocellum CBP process.

Lignin Populus trichocarpa CELF pretreatment CBP f↗

Toward Guided Mutagenesis: Gaussian Process Regression Predicts MHC Class II Antigen Mutant Binding

Antigen-specific immunotherapies (ASI) require successful loading and presentation of antigen peptides into the major histocompatibility complex (MHC) binding cleft. One route of ASI design is to mutate native antigens for either stronger or weaker binding interaction to MHC. Exploring all possible mutations is costly both experimentally and computationally. To reduce experimental and computational expense, here we investigate the minimal amount of prior data required to accurately predict the relative binding affinity of point mutations for peptide-MHC class II (pMHCII) binding. Using data from different residue subsets, we interpolate pMHCII mutant binding affinities by Gaussian process (GP) regression of residue volume and hydrophobicity. We apply GP regression to an experimental data set from the Immune Epitope Database, and theoretical data sets from NetMHCIIpan and Free Energy Perturbation calculations. We find that GP regression can predict binding affinities of nine neutral residues from a six-residue subset with an average R 2 coefficient of determination value of 0.62 ± 0.04 (±95% CI), average error of 0.09 ± 0.01 kcal/mol (±95% CI), and with an receiver operating characteristic (ROC) AUC value of 0.92 for binary classification of enhanced or diminished binding affinity. Similarly, metrics increase to an R2 value of 0.69 ± 0.04, average error of 0.07 ± 0.01 kcal/mol, and an ROC AUC value of 0.94 for predicting seven neutral residues from an eight-residue subset. Our work finds that prediction is most accurate for neutral residues at anchor residue sites without register shift. This work holds relevance to predicting pMHCII binding and accelerating ASI design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗