SEARCH · Search NASA

Results for “core optimization”

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 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Improved Safety Injection Flow Map Associated with Target RCS Depressurization to Maintain Core Coolability of OPR1000

In Korea, Severe Accident Management Guidance (SAMG) for Optimized Power Reactor (OPR1000) has been developed through the research results related to the WOG SAMG, in which the operation phase shifts from the Emergency Operation Procedure (EOP) to the SAMG when the Core Exit Temperature (CET) reaches 923 K. Among 7 mitigation strategies envisioned, the third management strategy, injecting coolant into Reactor Coolant System (RCS) by depressurizing the RCS is one of the key strategies to fortify the in-vessel retention capability. As an aid to such mitigation strategy, a safety injection (SI) flow map has been developed by Lee et al. for predicting required flow rate associated with target RCS depressurization to recover in-vessel coolability. The mechanistic SI flow map used the CET increase rate and core water level decrease rate obtained from MELCOR code simulation. It should be noted that during the severe accident, the total heat generated is sum of the decay heat and the oxidation heat. Thus, it is necessary to consider additional heat sources by oxidation reaction between cladding, structure, control rod, and hot steam prior to the injection of coolant. It was believed that the SI model estimated required flow rate rather conservatively by considering the expected amount decay heat and the accumulated heat in the core based on the RCS pressure. However, it was confirmed that the SI model developed by Lee et al. could not represent the total heat removal especially since start of oxidation. Under transient conditions, the heat sources accumulated in the heat structures (fuel, cladding, supporting structures) should be embarked in the model for accurate contribution of the total heat source to the required flow rate. Therefore, the objective of this study is to improve the previous SI flow map predicting the injection flow rate associated with target RCS depressurization to satisfy the in-vessel coolability during severe accidents. Additional terms were developed to reflect accurate heat accumulating and included in the existing model. Resulting SI flow map includes the target RCS depressurization, which enables discharge of the required flow rate by utilizing the flow rate-pressure curve with operation of two high pressure safety injections (HPSIs). Also recent MELCOR simulation result of a postulated steam generator tube rupture (SGTR) accident verified the accuracy of the improved map.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Optimization Methodology of Pebble Bed HTGR Start-Up and Running-in Strategy

In recent years interest in advanced reactor technologies has increased significantly. However, the methods used for analysis of traditional nuclear reactors are insufficient to consider all the different and varied advanced reactor designs without further development. One promising advanced reactor design is the pebble bed reactor (PBR). PBRs possess unique operational and fuel cycle features that require the development of specific analysis methodologies to adequately design and analyze the systems. There is a need in PBR research for a capability to analyze and optimize the process of transitioning from the start-up reactor core to the equilibrium core (known as the “running-in” of the reactor). The start-up of a PBR and the transition to the equilibrium core is a complex, multi-physics challenge that has not yet been well researched and has many opportunities for design, analysis, and optimization of the process. In this research, a methodology is defined to consider the potential strategies in PBR start-up and run-in to the equilibrium core. Multiple candidate software are considered, and their pros and cons are discussed for the PBR-specific application in the methodology. A preliminary software selection for the physics engine is made, and initial verification of key modules is performed. Software to support optimization of the reactor run-in strategies through reduced order modeling (ROM) and machine learning are considered. Challenges for a full implementation of the methodology are discussed. Additional code selections and verification of models relevant to the application are needed before full demonstration of the methodology can be achieved and an optimal strategy determined.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2006

Constraining Numerical Geodynamo Modeling with Surface Observations

Numerical dynamo solutions have traditionally been generated entirely by a set of self-consistent differential equations that govern the spatial-temporal variation of the magnetic field, velocity field and other fields related to dynamo processes. In particular, those solutions are obtained with parameters very different from those appropriate for the Earth s core. Geophysical application of the numerical results therefore depends on correct understanding of the differences (errors) between the model outputs and the true states (truth) in the outer core. Part of the truth can be observed at the surface in the form of poloidal magnetic field. To understand these differences, or errors, we generate new initial model state (analysis) by assimilating sequentially the model outputs with the surface geomagnetic observations using an optimal interpolation scheme. The time evolution of the core state is then controlled by our MoSST core dynamics model. The final outputs (forecasts) are then compared with the surface observations as a means to test the success of the assimilation. We use the surface geomagnetic data back to year 1900 for our studies, with 5-year forecast and 20-year analysis periods. We intend to use the result; to understand time variation of the errors with the assimilation sequences, and the impact of the assimilation on other unobservable quantities, such as the toroidal field and the fluid velocity in the core.

Kuang, Weijia↗

DOE OSTI2020

Feedback on Forty-year Long Clean-up Operations of a Contaminated Soil for Environmental Purpose - 20026

One of the CEA facilities in France is the place of miscellaneous mapping and clean-up operations since the 1980's. The final purpose of the conducted work deals with the environmental remediation of the building. In France, the absence of a regulatory framework for the management of sites polluted by radioactive substances, especially as the lack of release thresholds, has led to develop different approaches of soil radiological characterization. The geostatistical approach is part of them. Preliminary investigations are an unavoidable step prior to clearance and remediation. First, historical analysis enables to define site perimeters to investigate as well as to precisely locate expected contaminated parts. For radioactive pollutants, the media (water, soil, air) and transfer routes participate in the definition of the investigation perimeter and thus also need to be identified through a geological study. At the same time, functional analysis leads to choose the best measurement means which could provide useful information when combined with visual inspection. An equipment well adapted to easy-to-measure radionuclides as gamma emitters, such as gamma probe or on-site gamma spectrometry, is usually selected as a non-destructive assay. Coupled with spatial positions, processed on-site radiological data can be mapped. Data processing consists in mapping and kriging before a geostatistical analysis in order to identify the zones of interest to be targeted. Selected zones are not only linked to high level count rates but also to low level ones so that a pollution-free reference is known in the area. This data processing also enables to categorize waste types according to their origin and contamination levels. An expected volume of different waste then follows as well as the waste characterization equipments. In-depth investigations, e.g. core drillings and samplings, naturally ensue to aim at remediation optimization through the assessment of environmental impact. After clean-up operations, non-destructive assays as well as destructive samplings enable to check the effective clearance. Following this methodology in the present remediation work, results of the first non-destructive assay campaign highlighted a long-lived transuranic radionuclides contamination in the concrete flagstone. Successive non-destructive assays were then performed in the defined zones of interest. Count rates measured with a surface probe were first mapped. The best localization of destructive assays, i.e. core drillings at a 1 meter depth, came from this surface mapping followed by a geostatistical analysis. The concrete flagstone was totally removed. Then, the depth of soil to remove under the flagstone was optimized from the results of radiological activities measured in core drillings samples. A new surface mapping was then drawn after this partial soil excavation. The sand and demolition rubble samples, homogeneously constituted from the flagstone fragments, were analyzed by gamma spectrometry. Meanwhile, waste drums and other containers were also measured by gamma spectrometry and passive neutron measurement devices. All these results have brought miscellaneous pieces of information. In this context, major issues appear in terms of physical constraints such as the premises tininess as well as of soil sampling techniques and in terms of measurement performances assessment and the representativeness of homogeneous samples. The choice of radiological soil and waste characterization devices is described in the light of enhanced performances by keeping in mind radiation protection requirements and a willing to always optimize waste categories (very low level waste) and volume. The improvement of the characterization methods appears through this forty-year long work while sharing the feedback of encountered difficulties. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Microgrid Energy Management System Integration with Advanced Distribution Management System

The Integrated Distribution Management System (IDMS) project was initiated to demonstrate the interactive operation of microgrid systems and the distribution systems with which they interconnect. The key technologies for this are the microgrid management system and the utility distribution management system. The IDMS project successfully demonstrated that a utility’s advanced distribution management system/distributed energy resources management system (ADMS/DERMS) could effectively manage microgrids to provide visibility and control functionalities as well as use the microgrid as a dispatchable resource to support the utility grid. The DERMS accomplishes this by determining active and reactive power needs at the point of common coupling (PCC) using advanced applications like volt/VAR watt optimization (VVWO) and load relief (LR) with the underlying core applications state estimation (SE) and load flow (LF). The ADMS/DERMS can use the microgrid as a resource to resolve and prevent violations in the grid and to optimize the operational working state of the grid. The IDMS project demonstrated that a utility-operated ADMS with embedded DERMS functionality can flexibly manage a variety of microgrids and other aggregated distributed energy resources (DER) in concert with the wider distribution grid. Microgrids can provide grid services in any number of different ways to meet the operational needs of the distribution utility. The manner of aggregation — microgrid or virtual power plant — is not necessarily relevant to the utility as long as the grid services from the aggregated DER are available and can be managed by its ADMS/DERMS for the stability and reliability of the grid. The IDMS project demonstrated this integrated ADMS concept by combining hardware/software-in-the-loop testing with commercial products from different vendors and a utility’s network model. The project integrated the Schneider Electric EcoStruxure™ ADMS with DERMS with the Schweitzer Engineering Laboratories (SEL) POWERMAX ® Microgrid Control System, and the simulated resources with energy company PECO’s model of a utility-owned microgrid, establishing an operational relationship in which the utility manages the operational functionality of a microgrid at the PCC and that provides the utility with the capability to control the comprehensive power system, inclusive of the macrogrid and microgrid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2017

Efficient Machine Learning Approach for Optimizing Scientific Computing Applications on Emerging HPC Architectures

Efficient parallel implementations of scientific applications on multi-core CPUs with accelerators such as GPUs and Xeon Phis is challenging. This requires - exploiting the data parallel architecture of the accelerator along with the vector pipelines of modern x86 CPU architectures, load balancing, and efficient memory transfer between different devices. It is relatively easy to meet these requirements for highly structured scientific applications. In contrast, a number of scientific and engineering applications are unstructured. Getting performance on accelerators for these applications is extremely challenging because many of these applications employ irregular algorithms which exhibit data-dependent control-ow and irregular memory accesses. Furthermore, these applications are often iterative with dependency between steps, and thus making it hard to parallelize across steps. As a result, parallelism in these applications is often limited to a single step. Numerical simulation of charged particles beam dynamics is one such application where the distribution of work and memory access pattern at each time step is irregular. Applications with these properties tend to present significant branch and memory divergence, load imbalance between different processor cores, and poor compute and memory utilization. Prior research on parallelizing such irregular applications have been focused around optimizing the irregular, data-dependent memory accesses and control-ow during a single step of the application independent of the other steps, with the assumption that these patterns are completely unpredictable. We observed that the structure of computation leading to control-ow divergence and irregular memory accesses in one step is similar to that in the next step. It is possible to predict this structure in the current step by observing the computation structure of previous steps. In this dissertation, we present novel machine learning based optimization techniques to address the parallel implementation challenges of such irregular applications on different HPC architectures. In particular, we use supervised learning to predict the computation structure and use it to address the control-ow and memory access irregularities in the parallel implementation of such applications on GPUs, Xeon Phis, and heterogeneous architectures composed of multi-core CPUs with GPUs or Xeon Phis. We use numerical simulation of charged particles beam dynamics simulation as a motivating example throughout the dissertation to present our new approach, though they should be equally applicable to a wide range of irregular applications. The machine learning approach presented here use predictive analytics and forecasting techniques to adaptively model and track the irregular memory access pattern at each time step of the simulation to anticipate the future memory access pattern. Access pattern forecasts can then be used to formulate optimization decisions during application execution which improves the performance of the application at a future time step based on the observations from earlier time steps. In heterogeneous architectures, forecasts can also be used to improve the memory performance and resource utilization of all the processing units to deliver a good aggregate performance. We used these optimization techniques and anticipation strategy to design a cache-aware, memory efficient parallel algorithm to address the irregularities in the parallel implementation of charged particles beam dynamics simulation on different HPC architectures. Experimental result using a diverse mix of HPC architectures shows that our approach in using anticipation strategy is effective in maximizing data reuse, ensuring workload balance, minimizing branch and memory divergence, and in improving resource utilization.

Arumugam, Kamesh↗

DOE OSTI2020

Federal Best Practices: Core Principles of Sustainable Fleet Management

This document is designed to help agency fleet managers understand and implement optimal petroleum reduction strategies for each fleet location by evaluating the most appropriate combination of the four core principles of sustainable fleet management: Right-sizing the fleet to agency mission by implementing a vehicle allocation methodology study; Minimizing vehicle miles traveled; Increasing fleet fuel efficiency by replacing inefficient vehicles with more fuel-efficient vehicles, maintaining vehicles, driving more efficiently, and avoiding excessive idling; Optimizing cost-effective alternative fuel use, including maximizing utilization of existing alternative fuel infrastructure, installing alternative fuel infrastructure where practical (including electric vehicle charging stations), and aligning deployment of alternative fuel vehicles with fueling infrastructure.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2024

CORE DESIGN AND NEUTRONIC ANALYSIS OF THE EUROPEAN SODIUM FAST REACTOR WITH METALLIC FUEL

The current ESFR (European Sodium Fast Reactor) design was proposed and in-depth evaluated in the frame of the past ESFR-SMART project. As a follow-up project, the ESFR-SIMPLE has been launched with the aim of challenging the current commercial-size ESFR design in terms of safety features and economic performance. Among the new safety measures to be developed and assessed in ESFR-SIMPLE, the current oxide fuel ESFR design will be challenged by a modified version of the core with metallic fuel. This intends to conclude on what types of benefits can be obtained with high-density fuel, under similar safety and design constraints. In this paper, the designing approach for enabling the use of metallic fuel in the current ESFR core is described and a preliminary neutronic evaluation is carried out. The optimal configuration is established through the optimization of key neutronic parameters aiming at the potential reduction of the plutonium inventory. The resulting core configuration serves as a basis for further safety assessment analyses, which will provide insight into the advantages and drawbacks of the two types of fuels.

Jiménez-Carrascosa, Antonio↗

NASA NTRS2018

Performance Enhancements for the Lattice-Boltzmann Solver in the LAVA Framework

Performance enhancements in NASA's recently developed Lattice Boltzmann solver within the Launch Ascent and Vehicle Aerodynamics (LAVA) framework are presented. Two key algorithmic developments are highlighted. A coarse-fine interface treatment that discretely conserves mass and momentum has been implemented and successfully verified and validated. Code optimizations targeting improved serial and parallel performance were presented. For a simple turbulent Taylor-Green Vortex problem, we were able to demonstrate a 2.3 times speedup over the baseline code for a single Skylake-SP node containing 40 physical cores, and a 2.14 times speedup for 64 nodes containing 2560 physical cores. In addition, we were able to show that the optimizations enabled us to scale the code almost perfectly to 20480 physical cores where, including ghost cells, the problem size was 10 billion cells.

Barad, Michael↗

NASA NTRS2018

Performance Enhancements for the Lattice-Boltzmann Solver in the LAVA Framework

LAVA↗

DOE OSTI2025

Advanced Transmission Technologies – GETs and HPCs Session 2: Advanced Power Flow Control and Transmission Topology Optimization

The INL TADA GETs Cohort Session 2, held on November 7, 2025, conducted in collaboration with ScottMadden, focused on two core Advanced Transmission Technologies (ATTs): Advanced Power Flow Control (APFC) and Transmission Topology Optimization (TTO). These technologies are pivotal in enhancing grid flexibility, reliability, and cybersecurity resilience. APFC, particularly through modular FACTS devices like Modular Static Synchronous Series Compensators (M-SSSCs), enables dynamic voltage injection to reroute power flows. The session highlighted the deployment benefits of APFC, such as rapid installation, minimal civil works, and re-deployability. Regulatory drivers like FERC Order 2023 mandate the inclusion of Grid-Enhancing Technologies (GETs) in interconnection studies. Case studies from Central Hudson, CAISO, and National Grid (UK) demonstrated APFC’s effectiveness in congestion relief and cost savings. The session also addressed cybersecurity concerns, including firmware vulnerabilities, SCADA integration risks, and supply chain dependencies. Participants engaged in interactive exercises to rank cybersecurity and supply chain risks, emphasizing the need for robust digital assurance strategies. TTO involves software-based reconfiguration of transmission networks to optimize power flow without new infrastructure. The session showcased its operational value, with examples from SPP, PJM, and MISO showing significant congestion cost reductions. Cybersecurity vulnerabilities were discussed, particularly in API security and software supply chains, referencing incidents like SolarWinds and attacks on Danish utilities. Digital assurance exercises explored worst-case scenarios, attack paths, and mitigation responsibilities between vendors and utilities. Reliability challenges such as algorithm stability, vendor dependency, and operator trust were also examined. Cross-cutting themes emphasized the importance of digital assurance tools, including Software Bills of Materials (SBOMs) and hardware-in-loop testing. Human performance, training, and operational confidence were identified as critical enablers of technology adoption. The session concluded with a preview of Session 3, which will focus on High Performance Conductors (HPCs) and risk-based cybersecurity tools. Session 2 of 3.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Example Evaluation of a Representative Heat Pipe Test Article Design for Structural Acceptability using ASME Design Rules

This report describes a design analysis of a test article mimicking a core block of a heat-pipe mi-croreactor using the ASME Boiler & Pressure Vessel Code Section III, Division 5, Subsection HB, Subpart B rules. These rules cover the design and construction of elevated temperature nu-clear reactor structural components, like the core block of a heat-pipe microreactor emulated by the test article. The purposes of the report are to: (1) provide a step-by-step example of applying the ASME rules to a reasonable core block geometry and corresponding loading conditions, (2) to provide design feedback on the test article, and (3) to evaluate the ASME rules for use in the design of complex heat-pipe microreactor components. The test article design has a relatively short ASME design life, limited by high thermal stresses. The current ASME rules are adequate for evaluating the core block component designs, though they could be optimized to increase the design efficiency and reduce the amount of designer/analyst time required to execute the design process. The report covers both the current, 2019 edition ASME rules and modified, currently unapproved, rules designed to improve the design process for components with complex geometries and high secondary stresses, like the core block test article.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Direct optimization of neoclassical ion transport in stellarator reactors

Abstract We directly optimize stellarator neoclassical ion transport while holding neoclassical electron transport at a moderate level, creating a scenario favorable for impurity expulsion and retaining good ion confinement. Traditional neoclassical stellarator optimization has focused on minimizing ϵ eff , the geometric factor that characterizes the amount of radial transport due to particles in the 1 / ν regime. Under expected reactor-relevant conditions, core electrons will be in the 1 / ν regime and core fuel ions will be in the ν regime. Traditional optimizations thus minimize electron transport and rely on the radial electric field ( E r ) that develops to confine the ions. This often results in an inward-pointing E r that drives high- Z impurities into the core, which may be troublesome in future reactors. In this work, we increase the ratio of the thermal transport coefficients L 11 e / L 11 i , which previous research has shown can create an outward-pointing E r . This effect is very beneficial for impurity expulsion. We obtain self-consistent density, temperature, and E r profiles at reactor-relevant conditions for an optimized equilibrium. This equilibrium is expected to enjoy significantly improved impurity transport properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Near–Band–Edge Enhancement in Perovskite Solar Cells via Tunable Surface Plasmons

Plasmonic perovskite solar cells (PSCs) using core-shell type plasmonic particles are designed, which possess the plasmon resonance in the near-infrared range. This can selectively strengthen the interaction of the perovskite layer with low-energy photons. The mesoporous PSCs employing the plasmonic particles have delivered a 10%-15% enhancement of external quantum efficiency in the plasmonic resonance range. This surface-plasmonic effect has been analyzed using complementary techniques, including selective wavelength excitation and time-dependent photoluminescence. Here it is shown that the metal-based core-shell-type plasmonic structures in PSCs optimize the scattering and absorption of incident light and the dynamics of photogenerated carriers. Furthermore, both optical and electronic effects increase the power conversion efficiency of PSCs from 17.49% to 19.88%, paving a way toward controlling the thickness of the photoactive layer for advanced devices such as tandem solar cells.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Correlation consistent basis sets for actinides. II. The atoms Ac and Np–Lr

New correlation consistent basis sets optimized using the all-electron third-order Douglas-Kroll-Hess (DKH3) scalar relativistic Hamiltonian are reported for the actinide elements Ac and Np through Lr. These complete the series of sets reported previously for Th–U [K. A. Peterson, J. Chem. Phys. 142, 074105 (2015); M. Vasiliu et al., J. Phys. Chem. A 119, 11422 (2015)]. The new sets range in size from double- to quadruple-zeta and encompass both those optimized for valence (6s6p5f7s6d) and outer-core electron correlations (valence + 5s5p5d). The final sets have been contracted for both the DKH3 and eXact 2-component (X2C) Hamiltonians, yielding cc-pVnZ-DK3/cc-pVnZ-X2C sets for valence correlation and cc-pwCVnZ-DK3/cc-pwCVnZ-X2C sets for outer-core correlation (n = D, T, Q in each case). In order to test the effectiveness of the new basis sets, both atomic and molecular benchmark calculations have been carried out. In the first case, the first three atomic ionization potentials (IPs) of all the actinide elements Ac–Lr have been calculated using the Feller-Peterson-Dixon (FPD) composite approach, primarily with the multireference configuration interaction (MRCI) method. Excellent convergence towards the respective complete basis set (CBS) limits is achieved with the new sets, leading to good agreement with experiment, where these exist, after accurately accounting for spin-orbit effects using the 4-component Dirac-Hartree-Fock method. For a molecular test, the IP and atomization energy (AE) of PuO 2 have been calculated also using the FPD method but using a coupled cluster approach with spin-orbit coupling accounted for using the 4-component MRCI. The present calculations yield an IP 0 for PuO 2 of 159.8 kcal/mol, which is in excellent agreement with the experimental electron transfer bracketing value of 162 ± 3 kcal/mol. Likewise, the calculated 0 K AE of 305.6 kcal/mol is in very good agreement with the currently accepted experimental value of 303.1 ± 5 kcal/mol. The ground state of PuO 2 is predicted to be the 5 Σ$+\atop{0g}$ state.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

High-temperature superconducting CORC ® wires with record-breaking axial tensile strain tolerance present a breakthrough for high-field magnets

Cuprate high-temperature superconductors (HTS), such as RE-Ba 2 Cu 3 O 7–δ (REBCO, RE = rare earth), (Bi,Pb) 2 Sr 2 Ca 2 Cu 3 O 10–x and Bi 2 Sr 2 CaCu 2 O 8–x , have enabled the development of high-field superconducting magnets capable of generating magnetic fields far exceeding 20 T. The brittle nature of HTS requires elaborate means to protect them against the high stresses and strains associated with high-field magnet operation, and so far, has prevented reliable high-field HTS magnets from becoming a reality. Here we report a more than tenfold increase in the irreversible strain limit under axial tension (epsilonirr) to over 7% in optimized high-current conductor on round core (CORC®) conductors, compared to the REBCO tapes from which the CORC® conductor is wound. Minimizing the tape winding pitch of the helical wind mechanically decouples the brittle REBCO film from the overall conductor. The REBCO tapes behave as springs, limiting the rate at which applied strain is transferred to the ceramic film. In addition, high-strength alloy cores allow the critical stress (ϵ crit ) under axial tension at which initial degradation of CORC® conductors occurs to exceed 600 MPa, making them one of the strongest superconductors available. Mechanically decoupling the ceramic REBCO films from the overall CORC® conductor allows effective protection against the high operating stresses in high-field magnets. Furthermore, this breakthrough presents a monumental shift for HTS magnet technology, bringing reliable high-field superconducting magnets for compact fusion machines, the next generation of particle accelerators, and 40–60 T research solenoids within reach.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Graph-Based Representations and Applications to Process Simulation

Rapid and robust convergence of a process flowsheet is critical to enable large-scale simulations that address core scientific questions related to process design, optimization, and sustainability. However, due to the highly coupled and nonlinear nature of chemical processes, efficiently solving a flowsheet remains a challenge. In this work, we show that graph representations of the underlying physical phenomena in unit operations may help identify potential avenues to systematically reformulate the network of equations and enable more robust topology-based convergence of flowsheets. To this end, we developed graph abstractions of the governing equations of vapor-liquid and liquid-liquid equilibrium separation equipment. These graph abstractions consist of a mesh of interconnected variable nodes and equation nodes that are systematically generated through PhenomeNode, a new open-source library in Python developed in this study. We show that partitioning the graph into separate mass, energy, and equilibrium subgraphs can help decouple nonlinearities and guide decomposition algorithms. By employing the graph abstraction on an industrial separation process for separating glacial acetic acid from water, we implemented a new block decomposition scheme in BioSTEAM and demonstrated that this can accelerate convergence over a traditional sequential modular approach.

Distillation↗

NASA NTRS2016

Parallel Adjective High-Order CFD Simulations Characterizing SOFIA Cavity Acoustics

This paper presents large-scale MPI-parallel computational uid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft fuselage of a Boeing 747SP. These simulations focus on how the unsteady ow eld inside and over the cavity interferes with the optical path and mounting structure of the telescope. A temporally fourth-order accurate Runge-Kutta, and spatially fth-order accurate WENO- 5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh re nement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32k CPU cores and 4 billion compu- tational cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregular numerical cost associated with blocks con- taining boundaries. Limits to scaling beyond 32k cores are identi ed, and targeted code optimizations are discussed.

CFD↗