Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source Term Analysis”

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

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

At least 397 records · Page 22

Coupled Thermal-Hydraulic Analysis and Species Mass Transport in a Versatile Experimental Salt Irradiation Loop (VESIL) using COBRA-TF

With the resurgence of interest in Molten Salt Reactor (MSR) technology among commercial and governmental agencies, there is a need for new experiments and new modelling capabilities to characterize the unique phenomena present in this fluid fuel system. A Versatile Experimental Salt Irradiation Loop (VESIL) is currently under investigation at Idaho National Laboratory. It is designed to be a dynamic fuel-bearing salt loop placed in the Advanced Test Reactor. One of the key phenomena this proposed experiment plans to elucidate is fission product speciation in the fuel-salt and the subsequent effects this can have on fuel-salt properties, fission product transport and interactions, source term generation, corrosion and redox potential control. Specifically, noble gases (Xe & Kr) will bubble out to a plenum or off-gas system and noble metals (Mo, Tc, Te etc.) will precipitate and deposit in specific zones in the loop. This work extends the mass transfer and species interaction models in the general species transport capability of COBRA-TF (CTF) and applies these models to give a preliminary estimation of fission product behavior in the proposed VESIL design.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modeling Radionuclide Inventories in MSR Off-Gas Systems with Radiochemical Transport Analysis

Molten salt reactors (MSR) contain unique characteristics that may require enhancements to modeling tools to accurately predict phenomena. One characteristic that may be advantageous to leverage during normal operation is on-line processing of the circulating fuel salt, such as an off-gas system (OGS) to remove volatile fission products. Therefore, new modeling tools must be developed to integrate spatial resolution and chemistry effects into fuel depletion tools to be able to account for these non-core sources of radioactivity. Such types of radiochemical transport analysis tools were used to estimate the removal rates for 12 elements within a flow model of the Molten Salt Reactor Experiment (MSRE) by optimizing against legacy experimental data of the gas-borne (GB) percentages of 12 nuclides. The removal rates were used in a depletion model to calculate the FP inventory that enters the MSRE OGS. Calculations are in good agreement with the empirical GB percentages reported for the 12 nuclides, which validates the approach and verifies each tool’s treatment of the radiochemical flow effects. The OGS inventory is discussed in terms of the largest nuclide contributors to activity, dose consequence, decay heat, and elemental composition. Finally, insights from the study allow recommendations to be made for future code development activities.

Shahbazi, Shayan↗

ZMPY3D: accelerating protein structure volume analysis through vectorized 3D Zernike moments and Python-based GPU integration

Abstract Motivation Volumetric 3D object analyses are being applied in research fields such as structural bioinformatics, biophysics, and structural biology, with potential integration of artificial intelligence/machine learning (AI/ML) techniques. One such method, 3D Zernike moments, has proven valuable in analyzing protein structures (e.g., protein fold classification, protein–protein interaction analysis, and molecular dynamics simulations). Their compactness and efficiency make them amenable to large-scale analyses. Established methods for deriving 3D Zernike moments, however, can be inefficient, particularly when higher order terms are required, hindering broader applications. As the volume of experimental and computationally-predicted protein structure information continues to increase, structural biology has become a “big data” science requiring more efficient analysis tools. Results This application note presents a Python-based software package, ZMPY3D, to accelerate computation of 3D Zernike moments by vectorizing the mathematical formulae and using graphical processing units (GPUs). The package offers popular GPU-supported libraries such as CuPy and TensorFlow together with NumPy implementations, aiming to improve computational efficiency, adaptability, and flexibility in future algorithm development. The ZMPY3D package can be installed via PyPI, and the source code is available from GitHub. Volumetric-based protein 3D structural similarity scores and transform matrix of superposition functionalities have both been implemented, creating a powerful computational tool that will allow the research community to amalgamate 3D Zernike moments with existing AI/ML tools, to advance research and education in protein structure bioinformatics. Availability and implementation ZMPY3D, implemented in Python, is available on GitHub (https://github.com/tawssie/ZMPY3D) and PyPI, released under the GPL License.

Lai, Jhih-Siang (ORCID:0000000156775890)↗

Scalable Pattern Matching in Metadata Graphs via Constraint Checking

Pattern matching is a fundamental tool for answering complex graph queries. Unfortunately, existing solutions have limited capabilities: They do not scale to process large graphs and/or support only a restricted set of search templates or usage scenarios. Moreover, the algorithms at the core of the existing techniques are not suitable for today’s graph processing infrastructures relying on horizontal scalability and shared-nothing clusters, as most of these algorithms are inherently sequential and difficult to parallelize. In this article we present an algorithmic pipeline that bases pattern matching on constraint checking. The key intuition is that each vertex and edge participating in a match has to meet a set of constraints implicitly specified by the search template. These constraints can be verified independently and typically are less expensive to compute than searching the full template. The pipeline we propose generates these constraints and iterates over them to eliminate all the vertices and edges that do not participate in any match, thus reducing the background graph to a subgraph that is the union of all template matches—the complete set of all vertices and edges that participate in at least one match. Additional analysis can be performed on this annotated, reduced graph, such as full match enumeration, match counting, or computing vertex/edge centrality. Furthermore, a vertex-centric formulation for constraint checking algorithms exists, and this makes it possible to harness existing high-performance, vertex-centric graph processing frameworks. This technique (i) enables highly scalable pattern matching in metadata (labeled) graphs; (ii) supports arbitrary patterns with 100% precision; (iii) enables tradeoffs between precision and time-to-solution, while always selects all vertices and edges that participate in matches, thus offering 100% recall; and (iv) supports a set of popular data analytics scenarios. We implement our approach on top of HavoqGT, an open-source asynchronous graph processing framework, and demonstrate its advantages through strong and weak scaling experiments on massive scale real-world (up to 257 billion edges) and synthetic (up to 4.4 trillion edges) labeled graphs, respectively, and at scales (1,024 nodes / 36,864 cores), orders of magnitude larger than used in the past for similar problems. This article serves two purposes: First, it synthesises the knowledge accumulated during a long-term project. Second, it presents new system features, usage scenarios, optimizations, and comparisons with related work that strengthen the confidence that pattern matching based on iterative pruning via constraint checking is an effective and scalable approach in practice. The new contributions include the following: (i) We demonstrate the ability of the constraint checking approach to efficiently support two additional search scenarios that often emerge in practice, interactive incremental search and exploratory search. (ii) We empirically compare our solution with two additional state-of-the-art systems, Arabsque and TriAD. (iii) We show the ability of our solution to accommodate a more diverse range of datasets with varying properties, e.g., scale, skewness, label distribution, and match frequency. (iv) We introduce or extend a number of system features (e.g., work aggregation, load balancing, and the ability to cap the generated traffic) and design optimizations and demonstrate their advantages with respect to improving performance and scalability. (v) We present bottleneck analysis and insights into artifacts that influence performance. (vi) We present a theoretical complexity argument that motivates the performance gains we observe.

97 MATHEMATICS AND COMPUTING↗

Advancements in Multiphysics Microdepletion Analysis of an eVinci TM -like Microreactor Leveraging OpenMC-CRAB Workflow

Nuclear microreactors (MRs) are a class of nuclear reactor technology, characterized by reduced dimensions, modular design, and reduced power output in contrast to conventional Light Water Reactors (LWRs). MRs are proposed for supplying electricity and eventual process heat to remote locations, such as military installations and disaster-affected areas. Current research work sponsored by the US Department of Energy Microreactor Program (MRP) is devoted to the development of novel modeling and simulation tools to better support MR vendors and regulatory bodies. Notably, the NRC is projected to utilize the CRAB multiphysics software driver for executing both design and beyond-design-basis accident analyses. Furthermore, the NRC has been utilizing the MELCOR code to calculate mechanistic source terms during accidents. Since MELCOR relies on isotopic inventory and reactor temperature/power profiles under accident conditions, which theoretically can be derived from CRAB, the goal is to establish a comprehensive CRAB-MELCOR computational framework. Past work was focused on testing and demonstrating CRAB's capability to generate results that can be used to inform mechanistic source term calculations in MELCOR. In particular, a computational workflow leveraging OpenMC-generated microscopic cross sections and CRAB was first applied to perform multiphysics microscopic depletion calculation followed by an accident scenario for a stylized microreactor problem. In fiscal year 2024, the research work has been focused on applying the OpenMC-CRAB workflow, which was first tested in fiscal year 2023, to a realistic 3D heat-pipe cooled MR problem representative of the eVinci TM design. The latter computational problem was developed with inputs from WEC to conserve selected neutronic and thermal characteristics of the eVinci TM design without releasing proprietary data. The results of this simulation, encompassing isotopic inventory, power density distribution, and kinetic parameters, will inform both MELCOR and the WEC-developed FATE code for mechanistic source terms calculations. The results from the two codes will then be compared for code verification purposes. This report contains the design characteristics of the realist heat pipe cooled microreactor developed as a use-case for the verification exercise, and the current results for the multiphysics microscopic depletion performed with the OpenMC-CRAB workflow. The results include eigenvalue as a function of time, power distribution at EOL, in addition to nuclides inventory's time evolution and spatial distribution. Finally, we report improvements to the workflow efficiency achieved through a collaboration with the NEAMS programs. Through this collaborative effort, we were able to strongly decrease the computational time for the multiphysics microdepletion calculation (i.e., from 17.4 hours to 5.7 hours on 280 processors) in addition to simplifying the interface to generate isotopics spatial distribution utilizable by FATE and MELCOR. Future work, including the improvement of the current microscopic cross-sections' library and the simulation of an accident scenario at EOL, is also discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Froggatt-Nielsen meets the SMEFT

We study the matching of Froggatt-Nielsen theories of flavour onto the Standard Model Effective Field Theory (SMEFT), upon integrating out a heavy Beyond-the-Standard-Model (BSM) scalar ‘flavon’ whose vacuum expectation value breaks an Abelian flavour symmetry at energies Λ FN well above the electroweak scale, Λ FN > Λ SM . We include matching contributions to the infrared d SM = 6 (Warsaw basis) SMEFT sourced from ultraviolet contact terms suppressed up to order 1/${\Lambda}_{\textrm{UV}}^2$ in the Froggatt-Nielsen Lagrangian, where Λ UV > Λ FN is an arbitrary ultraviolet scale where further unspecified BSM particles are dynamical. This includes tree-level (one-loop) ultraviolet diagrams with d FN = 6 (5) effective vertices. We first do so with a toy model, but then generalize our findings to arbitrary Frogatt-Nielsen charges. Our results indicate a rich and non-trivial signature of Froggatt-Nielsen theories on the (otherwise) model-independent operators of the SMEFT, and we briefly speculate on extending our analysis to broader classes of BSM flavour models, e.g. non-Abelian and/or gauged theories. We thus take an important step towards determining how to use rapidly developing theoretical and experimental SMEFT technologies to gain unambiguous insight into the SM’s longstanding fermion flavour puzzle.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling of Am-241 as an alternative fuel source in a radioisotope thermoelectric generator

Radioisotope thermoelectric generators (RTGs) have been utilized in the USA to power satellites and space exploration equipment for over half a century. In our preceding paper, Dustin and Borrelli (2021), we assessed the performance of 241 Am as an alternative to 238 Pu by modeling the specific decay-heat and gamma spectra of several potential radioisotopes in SCALE-ORIGEN, and comparing these with 238 Pu An analysis of heat output, heat stability, safety, availability and required mass showed that of the isotopes reviewed, 241 Am had the greatest potential for replacing 238 Pu, with additional design considerations. In our current work, we have assessed radiation transport for two configurations of an 241 Am powered Multi-Mission Radioisotope Thermoelectric Generator (MMRTG), due to the potential health effects for operators, fabricators, and others with the potential to be exposed to neutrons and lower energy gamma rays. We have modeled the use of americium as a fuel source replacement for plutonium in current MMRTG designs. The assessment has been performed using MCNP6, in terms of neutron and gamma emissions that might interfere with sensitive equipment and power system production. The radiation transport behavior in terms of dose rate of the 241 Am fueled models is compared with that of a 238 Pu fueled model. Finally, additional work was performed to validate modeling simplifications with regards to thermocouple material configurations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Mechanism of Increasing Summer Water Vapor Over the Tibetan Plateau

Rapid warming over the Tibetan Plateau (TP) has been associated with a wetting trend in atmospheric water vapor content, which is critical for recharging the Asian water tower. However, the mechanism associated with the wetting phenomenon remains unclear. Long-term changes in moisture balance and precipitation recycling processes are investigated using the ERA5 reanalysis from 1979 to 2019. The wetting trend over TP is mainly due to the summer water vapor trend with signi?cantly increases over the western TP. Based on the moisture balance analysis, it is found that precipitation, evaporation, and the convergence of moisture are all increasing over the western TP but decreasing over the eastern TP. Based on the dynamical precipitation recycling model, the results suggest that both internal and external cycles contributes to the wetting over TP, with 63.87% and 36.13% contributions respectively. Further analysis found that the atmospheric heating source is also increasing over the western TP, which could shift the moisture transportation from east to west at the southern boundary of TP. The increasing moisture convergence could enhance precipitation, and the enhanced latent heating in the mid-atmosphere can further induced moisture convergence, which forms a positive feedback. However, an opposite situation occurred over the eastern TP. The internal and external cycle of the recycling process can stimulate (suppress) each other through precipitation over the western (eastern) TP. This mechanism linked the changes in the precipitation recycling and atmospheric circulation, and induced the wetting trend over TP in summer.

54 ENVIRONMENTAL SCIENCES↗

Distribution Substation Planning Toolkit (dsp-toolkit) v1.0

The Distribution Substation Planning Toolkit (DSP Toolkit) is a software suite designed to streamline the planning and optimization of distribution substations. This toolkit offers a comprehensive set of tools and APIs for data curation, short-term electric load forecasting, and weather-sensitive load adjustment, making it an essential resource for utility companies, engineers, and researchers. Features • Data Preprocessing and Curation: Efficiently manage and preprocess large datasets to ensure high-quality input for analysis. • Short-Term Load Forecasting: Utilize data-driven models to predict short-term electric loads accurately. • Weather-Sensitive Modeling: Automatically adjust load forecasts based on weather data to predict future peak demands more precisely. Uses The DSP Toolkit is ideal for planning and optimizing distribution substations, providing a user-friendly interface and comprehensive documentation. It is suitable for both novice and experienced users, facilitating efficient and accurate planning processes. Advantages • Efficiency: Automates complex planning tasks, reducing manual effort and minimizing errors. • Scalability: Handles large datasets and complex models, making it suitable for large-scale projects. • Community and Support: Open-source with active community contributions, ensuring continuous improvement and support. • Extensibility: Easily extendable with custom modules and plugins, allowing users to tailor the toolkit to their specific needs. The DSP Toolkit stands out by offering a robust, flexible, and user-friendly solution for distribution substation planning. Public Abstract

Li, Han [Lawrence Berkeley National Laboratory (LB↗

Observed Relationships Between Cloud Droplet Effective Radius and Biogenic Gas Concentrations in Summertime Marine Stratocumulus Over the Eastern North Atlantic

Abstract Biogenic gases are a prominent component of the summertime marine boundary layer (MBL) over the Eastern North Atlantic. One of these gases, dimethyl sulfide (DMS), can produce sulfate cloud condensation nuclei (CCN) that, in theory, can brighten clouds through photolysis, and produces a reaction product, methane sulfonic acid (MSA). It is also possible that DMS can interact with sea‐salt or other marine aerosols changing their CCN activation spectrum, which could also modify cloud microphysical structure. Data collected aboard the G1 aircraft during the Aerosol Cloud Experiment Eastern North Atlantic (ACE‐ENA) in well‐mixed and decoupled marine boundary layers (MBLs) were used to examine relationships between the cloud droplet effective radii, , and the concentrations of DMS and MSA in constant cloud liquid water content (LWC) bins. A weak but statistically significant negative correlation was observed between CCN concentration and in most LWC bins, regardless of the source of the CCN, while a weak but statistically significant positive correlation between and DMS was observed. No correlation between the cloud droplet number concentration and DMS was found. The presence of MSA indicated that DMS‐to‐sulfate photolysis was likely occurring, but data sparsity prevented a statistically significant conclusion regarding the relationship between MSA and . Data sparsity in decoupled conditions also prevented statistically significant conclusions. To properly address biogenic gas impacts on cloud microphysics, it is recommended that aircraft data be supplemented by long‐term biogenic gas measurements at the surface in marine locations with appropriate remote and in‐situ cloud sensing capabilities, and the analysis limited to well‐mixed MBL's.

54 ENVIRONMENTAL SCIENCES↗

Off-axis representation of hyperbolic mirror shapes for X-ray beamlines

Mirror-centered, closed-form expressions for hyperbolic surfaces used in X-ray beamlines have been derived. Hyperbolic mirrors create a virtual focus or source point and can be used to lengthen or shorten the effective focal distance of a compound optical system. The derivations here express off-axis segments of a hyperbolic surface in terms of the real and virtual focal distances and the incident glancing angle at the center of the mirror. Conventional mathematical expressions of hyperbolic shapes describe the surfaces in Cartesian or polar coordinates centered on an axis of symmetry, necessitating cumbersome rotation and translation to mirror-centered coordinates. The representation presented here, with zero slope and the origin at the central point, is most convenient for modeling, metrology, aberration correction, and general surface analysis of off-axis configurations. The direct derivation avoids the need for nested coordinate transforms. A series expansion provides a helpful approximation; the coefficients of the implicit equation are also provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optimal dimensionality selection for independent component analysis of transcriptomic data

Independent component analysis is an unsupervised machine learning algorithm that separates a set of mixed signals into a set of statistically independent source signals. Applied to high-quality gene expression datasets, independent component analysis effectively reveals both the source signals of the transcriptome as co-regulated gene sets, and the activity levels of the underlying regulators across diverse experimental conditions. Two major variables that affect the final gene sets are the diversity of the expression profiles contained in the underlying data, and the user-defined number of independent components, or dimensionality, to compute. Availability of high-quality transcriptomic datasets has grown exponentially as high-throughput technologies have advanced; however, optimal dimensionality selection remains an open question. We computed independent components across a range of dimensionalities for four gene expression datasets with varying dimensions (both in terms of number of genes and number of samples). We computed the correlation between independent components across different dimensionalities to understand how the overall structure evolves as the number of user-defined components increases. We then measured how well the resulting gene clusters reflected known regulatory mechanisms, and developed a set of metrics to assess the accuracy of the decomposition at a given dimension. We found that over-decomposition results in many independent components dominated by a single gene, whereas under-decomposition results in independent components that poorly capture the known regulatory structure. From these results, we developed a new method, called OptICA, for finding the optimal dimensionality that controls for both over- and under-decomposition. Specifically, OptICA selects the highest dimension that produces a low number of components that are dominated by a single gene. We show that OptICA outperforms two previously proposed methods for selecting the number of independent components across four transcriptomic databases of varying sizes. OptICA avoids both over-decomposition and under-decomposition of transcriptomic datasets resulting in the best representation of the organism’s underlying transcriptional regulatory network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of methods for characterizing the fine particulate matter emissions from aircraft and other diffusion flame combustion aerosol sources

The U. S. Environmental Protection Agency in collaboration with the U. S. Air Force Arnold Engineering Development Complex conducted the VAriable Response In Aircraft nvPM Testing (VARIAnT) 3 and 4 test campaigns to compare nvPM emissions measurements from a variety of diffusion flame combustion aerosol sources (DFCASs), including a Cummins diesel engine, a diesel powered generator, two gas turbine start carts, a J85-GE-5 turbojet engine burning multiple fuels, and a Mini-CAST soot generator. The VARIAnT research program has conducted four test campaigns to date with the VARIAnT 3 and 4 test campaigns focusing on BC mass instrument performance, use of gas turbine engine sources for BC mass instrument calibration, and the characterization of the aerosols produced from the combustion sources. Results from both campaigns revealed agreement of about 20% between the AVL Micro-Soot Sensor, the Cavity Attenuated Phase Shift (CAPS PMSSA) monitor and the thermal-optical reference method for elemental carbon (EC), independent of the calibration source used. For the Artium LII-300, the measured mass concentrations in VARIAnT3 fall within 18% and in VARIAnT4 fall within 32% of the reference EC mass concentration when calibrated on a combustor rig in VARIAnT3 and on a LGT-60 start cart in VARIAnT4, respectively. It was also found that the three mass instrument types (MSS, CAPS PMSSA, and LII-300) can exhibit different BC to reference EC ratios depending on the emission source that appear to correlate to particle geometric mean mobility diameter, morphology, or some other parameter associated with particle geometric mean diameter (GMD) with the LII-300 showing a slightly stronger apparent trend with GMD. Systematic differences in LII-300 measured mass concentrations have been reduced by calibrating with a turbine combustion source (combustor or turbine engine). With respect to the particle size measurements made, the sizing instruments (TSI SMPS, TSI EEPS, and Cambustion DMS 500) were found to be in general agreement in terms of size distributions and concentrations with some exceptions. Teflon filter measurements of the total aerosol mass produced by the various DFCASs differed from the reference EC, BC, and integrated particle sized distribution measured aerosol masses. The Teflon filter measured mass was 30 to 40% greater than the reference EC and 50 to 60% greater than the BC mass for the J85 in both test campaigns. The measurements of particle size distributions and single particle analysis by miniSPLAT indicated the presence of large particles (>100 nm) having more compact morphologies, higher effective density, and a composition dominated by OC and ash. This increased large particle fraction is also associated with higher values of single scattering albedo measured by the CAPS PMSSA instrument and higher OC/EC measurements. These measurements indicate gas turbine engine emissions can be a more heterogeneous mix of particle types beyond the original E31 assumption that it is mainly composed of black carbon.

33 ADVANCED PROPULSION SYSTEMS↗

Finalizing Transition to the New Data Center at BNL

Computational science, data management and analysis have been key factors in the success of Brookhaven National Laboratory's scientific programs at the Relativistic Heavy Ion Collider (RHIC), the National Synchrotron Light Source (NSLS-II), the Center for Functional Nanomaterials (CFN), and in biological, atmospheric, and energy systems science, Lattice Quantum Chromodynamics (LQCD) and Materials Science, as well as our participation in international research collaborations, such as the ATLAS Experiment at Europe's Large Hadron Collider (LHC) at CERN (Switzerland) and the Belle II Experiment at KEK (Japan). The construction of a new data center is an acknowledgement of the increasing demand for computing and storage services at BNL in the near term and enable the Lab to address the needs of the future experiments at the High-Luminosity LHC at CERN and the Electron-Ion Collider (EIC) at BNL in the long term. The Computing Facility Revitalization (CFR) project is aimed at repurposing the former National Synchrotron Light Source (NSLS-I) building as the new data center for BNL. The construction of the new data center was finished in 2021Q3, and it was delivered for production in early FY2022 for all collaborations supported by the Scientific Data and Computing Center (SDCC), including STAR, PHENIX and sPHENIX experiments at RHIC collider at BNL, the Belle II Experiment at KEK (Japan), and the Computational Science Initiative at BNL (CSI). This paper highlights the key mechanical, electrical, and networking components of the new data center in its final configuration as used in production since 2021Q4 and gives an overview for the extension of the central network systems into the new data center and the migration of a significant portion of IT load and services from the old data center to the new data center carried out in 20212023, with expected completion of the main phase of the gradual IT equipment replacement and migration from the old data center into the new one set to the end of FY2023 (Sep 30, 2023).

99 GENERAL AND MISCELLANEOUS↗

Efficient Space–Time Reduced Order Model for Linear Dynamical Systems in Python Using Less than 120 Lines of Code

A classical reduced order model (ROM) for dynamical problems typically involves only the spatial reduction of a given problem. Recently, a novel space–time ROM for linear dynamical problems has been developed [Choi et al., Space–tume reduced order model for large-scale linear dynamical systems with application to Boltzmann transport problems, Journal of Computational Physics, 2020], which further reduces the problem size by introducing a temporal reduction in addition to a spatial reduction without much loss in accuracy. The authors show an order of a thousand speed-up with a relative error of less than 10−5 for a large-scale Boltzmann transport problem. In this work, we present for the first time the derivation of the space–time least-squares Petrov–Galerkin (LSPG) projection for linear dynamical systems and its corresponding block structures. Utilizing these block structures, we demonstrate the ease of construction of the space–time ROM method with two model problems: 2D diffusion and 2D convection diffusion, with and without a linear source term. For each problem, we demonstrate the entire process of generating the full order model (FOM) data, constructing the space–time ROM, and predicting the reduced-order solutions, all in less than 120 lines of Python code. We compare our LSPG method with the traditional Galerkin method and show that the space–time ROMs can achieve O(10−3) to O(10−4) relative errors for these problems. Depending on parameter–separability, online speed-ups may or may not be achieved. For the FOMs with parameter–separability, the space–time ROMs can achieve O(10) online speed-ups. Finally, we present an error analysis for the space–time LSPG projection and derive an error bound, which shows an improvement compared to traditional spatial Galerkin ROM methods.

97 MATHEMATICS AND COMPUTING↗

Simulating Catalysis with Realistic Pellet Geometries Using Mesoflow: A Case Study of Catalytic Propane Dehydrogenation

We present a case study of catalytic propane dehydrogenation with our open-source multiphysics solver, Mesoflow. The solver was developed to simulate reactive flow coupled to heterogeneous catalytic reactions and deactivation in the context of complex, mesoscale geometry. The method leverages cartesian block-structured adaptive mesh refinement to capture realistic catalyst microstructural features acquired directly from X-ray computed tomography data. A kinetic model for propane dehydrogenation and catalyst deactivation was developed based on temporal analysis of products (TAP) reactor experiments. The TAP reactor experiments allow for precise characterization of intrinsic kinetic reaction steps which are implemented into Mesoflow simulations to model the spatial and temporal evolution of reactants, products, and catalyst active sites. The short-term and long-term deactivation behavior is studied by using XCT data collected from fresh and aged catalyst pellets, which exhibit different microstructural features. This study employs time-splitting algorithms to connect disparate reaction and flow timescales, enabling the simulations to achieve realistic deactivation timescales on the order of minutes while the flow time-scales for small particles (100 microns) are several milliseconds. We also introduce a flexible automated python script that writes the necessary files to construct a Mesoflow simulation from user-created chemical mechanisms. We will also introduce a few new features that are added to Mesoflow such as higher order schemes, implicit chemistry integrators and the ability to run on AMD and NVIDIA graphics-processing-units.

AMReX↗

Molten salt electrolysis and room temperature ionic liquid electrochemical processes for refining rare earth metals: Environmental and economic performance comparison

Circularizing end of life products such as electronic waste can open a new source for rare earth elements (REEs) that will decrease the demand pressure on conventional virgin production. The last step to produce REEs consists of a refining process to convert rare earth oxide (REO) to REE metal. One promising technology is processing REO via a room temperature ionic liquid, which is characterized by low energy requirements in contrast to existing REEs extraction processes. In this work, life cycle and techno-economic assessments are performed and compared against molten salt electrolysis. Here, the results show that both processes have advantages and disadvantages in terms of environmental performance, and that they are similarly competitive in terms of economic performance. A breakeven analysis suggests that future research should focus on coupling the production of REOs with the refining processes to attempt to lower REO cost and make them economically feasible in the U.S.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbon Capture, Transport, And Storage (CTS) Cost Modeling Of The Onshore Gulf Coast

The onshore Gulf of Mexico region presents significant opportunities for CO2 capture, transport, and storage due to its numerous CO2 sources, such as power plants, refineries, and its substantial CO2 storage potential. However, operators face critical decisions in designing an efficient and cost-effective CO2 pipeline network. This study examines the economic implications of two primary strategies: constructing a trunkline with excess initial capacity versus developing dedicated pipelines incrementally as new CO2 sources come online. Building a trunkline first offers the advantage of future-proofing the network, allowing for the accommodation of increased CO2 volumes from various sources over time. However, this approach incurs higher upfront costs and risks underutilizing the transport capacity in the initial stages, potentially resulting in economic inefficiencies. Conversely, constructing dedicated pipelines for each new CO2 source as it becomes operational may avoid the initial overcapacity issue but fails to capitalize on the economies of scale. This could lead to higher overall costs due to the duplication of infrastructure and increased complexity in network management. This research employs a comprehensive cost-benefit analysis, integrating factors such as capital expenditure, operational costs, projected CO2 volumes, and potential economies of scale. Through this analysis, we aim to provide operators with insights into the most economically viable strategy for CO2 pipeline network design in the region. The findings underscore the importance of strategic planning and highlight the trade-offs between immediate capacity utilization and long-term cost savings, ultimately guiding stakeholders towards informed decision-making in the development of CO2 transport infrastructure. Presented at the 41st USAEE/IAEE North American Conference, 3-6 November 2024, Baton Rouge, LA, United States.

Shih, Chung Yan↗