Search NASA⌕ Search

SEARCH · Search NASA

Results for “formulation study”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

ENSO Diversity Simulated in a Revised Cane-Zebiak Model

The El Niño-Southern Oscillation (ENSO) phenomenon features rich sea surface temperature (SST) spatial pattern variations dominated by the Central Pacific (CP) and Eastern Pacific (EP) patterns during its warm phase. Understanding such ENSO pattern diversity has been a subject under extensive research activity. To provide a framework for unveiling the fundamental dynamics of ENSO diversity, an intermediate coupled model based on the Cane-Zebiak-type framework, named RCZ, is established in this study. Compared with the original Cane-Zebiak model, RCZ consists of revised model formulation and well-tuned parameterization schemes. All model components are carefully validated against the observations via the standalone mode, in which the observed anomalous SST (wind stress) forcing is prescribed to drive the atmospheric (oceanic) component. The superiority of RCZ’s model components over those in the original Cane-Zebiak model is evidenced by their better performance in simulating the observations. Coupled simulation with RCZ satisfactorily reproduces aspects of the observed ENSO characteristics, including the spatial pattern, phase-locking, amplitude asymmetry, and, particularly, ENSO diversity/bi-modality. RCZ serves as a promising tool for studying dynamics of ENSO diversity as it resolves most of the relevant processes proposed in the literature, including atmospheric nonlinear convective heating, oceanic nonlinear dynamical heating, and the ENSO/westerly wind burst interaction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Toward Predictive RANS and SRS Computations of Turbulent External Flows of Practical Interest

In this work, we investigate the main challenges to prediction of turbulent external flows of practical interest with Reynolds-Averaged Navier–Stokes equations (RANS) and Scale-Resolving Simulation (SRS) models. This represents a crucial step toward further developing and establishing these formulations so they can be confidently utilized in engineering problems without reference data. The study initiates by identifying the major challenges to prediction. A literature review is performed to illustrate their effects in RANS and SRS computations. Afterward, we evaluate the impact of the challenges to prediction by analyzing representative statistically steady and unsteady flows with prominent RANS and SRS methods. These include multiple turbulent viscosity and second-moment RANS closures, and hybrid and bridging SRS models. The results demonstrate the potential of the selected SRS models to predict engineering flows. Yet, they also show the importance of considering the challenges to prediction during the setup and conduction of numerical experiments. These can suppress the advantages of using SRS formulations. The data also indicate that only SRS models can confidently predict statistically unsteady flows. In contrast, the results demonstrate that mean-flow quantities of statistically steady flows can be efficiently calculated with RANS closures, especially second-moment closures. Among the selected SRS methods, bridging models reveal better suited for prediction due to their ability to prevent commutation errors and enable the robust evaluation of numerical and modeling errors. This last property allows the use of a new validation technique that does not require reference data.

42 ENGINEERING↗
DOE OSTI2020

Applying Improved Optical Recognition with Machine Learning on Sorting Cu Impurities in Steel Scrap

Annually, 20–55 million tons of electronic waste (e-waste) is produced worldwide (5% of all municipal solid waste). Although e-waste embodies only 2% of America’s municipal waste, it accounts for a significantly larger proportion of the heavy metals and flame retardants present in the waste stream. Currently, less than 20% of all e-waste is recycled in the United States because the heterogeneity of the feedstock limits the opportunity for reuse in high value products and processing of the waste itself is often too costly to justify handling. To address this concern, the aims of this study are: 1) to identify the major plastic and metal compositions within electronic shredder residue (ESR), 2) to formulate solvents and processing conditions to separate 90% of the plastics targeted from consumer shred ESR, and 3) to develop a process design model to estimate the cost and energy efficiency of the proposed solvent-based processing. In this study, a pre-sorted heterogeneous ESR feedstock (one where aluminum, magnetic components, and hazardous battery materials removed by an e-waste recycling facility) was used, with the major compositions of the ESR characterized. It was found that 25 wt.% of the feedstock was composed of plastics, 6 wt.% rubber, 27 wt.% printed circuit boards, 23% wire, and the remainder metals and capacitors. Within the plastic portion, polystyrene (PS, 40 wt.%), acrylonitrile butadiene styrene (ABS, 25 wt.%), and styrene-acrylonitrile (SAN, 9wt.%) were identified to compose the majority of the screened plastics using Fourier transform infrared spectroscopy (FTIR). Next, selective solvents were screened using Hansen Solubility Parameter Theory (HSP) for dissolving PS and ABS. The pre-screening results show that methylene chloride (dichloromethane, DCM) and tetrahydrofuran (THF) are capable of dissolving the most PS and ABS, while methanol (MeOH) and ethylene glycol (EG) are capable of precipitating the most PS and ABS. These solvents were subsequently used to recover polymers and remove. flame retardants within the ESR feedstock. By optimizing the dissolution time and the solvents used, the highest polymer dissolution yield (99 wt.%%) was achieved using DCM for 48 hr. Both pre-screened anti-solvents (MeOH and EG) showed the highest polymer precipitation yield (71 wt.%). In terms of flame retardant removal rate, EG was found to have a high phosphorus-containing flame retardant removal rate (up to 98%). Characterization shows that the proposed solvent-based processing can preserve a high molecular weight fraction of the polymers and effectively remove flame retardants. Cost analysis indicates that the amount of the solvent/anti-solvent recovered after the reaction would play a critical role in reducing the operating costs. The energy analysis shows that the proposed solvent-based processes can save up to 60% of the embodied energy used to manufacture plastics used in electronics (PS and ABS were used for calculations). The results from this project prove the potential of solvent-based processing to produce secondary materials (plastics and metals) from e-waste for cross-industry reuse.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Chemical Recycling of Mixed Plastics and Valuable Metals in the Electronic Waste Using Solvent-Based Processing

Annually, 20-55 million tons of electronic waste (e-waste) is produced worldwide (5% of all municipal solid waste). Although e-waste embodies only 2% of America’s municipal waste, it accounts for a significantly larger proportion of the heavy metals and flame retardants present in the waste stream. Currently, less than 20% of all e-waste is recycled in the United States because the heterogeneity of the feedstock limits the opportunity for reuse in high value products and processing of the waste itself is often too costly to justify handling. To address this concern, the aims of this study are: 1) to identify the major plastic and metal compositions within electronic shredder residue (ESR), 2) to formulate solvents and processing conditions to separate 90% of the plastics targeted from consumer shred ESR, and 3) to develop a process design model to estimate the cost and energy efficiency of the proposed solvent-based processing. In this study, a pre-sorted heterogeneous ESR feedstock (one where aluminum, magnetic components, and hazardous battery materials removed by an e-waste recycling facility) was used, with the major compositions of the ESR characterized. It was found that 25 wt.% of the feedstock was composed of plastics, 6 wt.% rubber, 27 wt.% printed circuit boards, 23% wire, and the remainder metals and capacitors. Within the plastic portion, polystyrene (PS, 40 wt.%), acrylonitrile butadiene styrene (ABS, 25 wt.%), and styrene-acrylonitrile (SAN, 9wt.%) were identified to compose the majority of the screened plastics using Fourier transform infrared spectroscopy (FTIR). Next, selective solvents were screened using Hansen Solubility Parameter Theory (HSP) for dissolving PS and ABS. The pre-screening results show that methylene chloride (dichloromethane, DCM) and tetrahydrofuran (THF) are capable of dissolving the most PS and ABS, while methanol (MeOH) and ethylene glycol (EG) are capable of precipitating the most PS and ABS. These solvents were subsequently used to recover polymers and remove. flame retardants within the ESR feedstock. By optimizing the dissolution time and the solvents used, the highest polymer dissolution yield (99 wt.%%) was achieved using DCM for 48 hr. Both pre-screened anti-solvents (MeOH and EG) showed the highest polymer precipitation yield (71 wt.%). In terms of flame retardant removal rate, EG was found to have a high phosphorus-containing flame retardant removal rate (up to 98%). Characterization shows that the proposed solvent-based processing can preserve a high molecular weight fraction of the polymers and effectively remove flame retardants. Cost analysis indicates that the amount of the solvent/anti-solvent recovered after the reaction would play a critical role in reducing the operating costs. The energy analysis shows that the proposed solvent-based processes can save up to 60% of the embodied energy used to manufacture plastics used in electronics (PS and ABS were used for calculations). The results from this project prove the potential of solvent-based processing to produce secondary materials (plastics and metals) from e-waste for cross-industry reuse.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Flash NanoPrecipitation as an Agrochemical Nanocarrier Formulation Platform: Phloem Uptake and Translocation after Foliar Administration

The increasing severity of pathogenic and environmental stressors that negatively affect plant health has led to interest in developing next-generation agrochemical delivery systems capable of precisely transporting active agents to specific sites within plants. In this work, we adapt Flash NanoPrecipitation (FNP), a scalable nanocarrier (NC) formulation technology used in the pharmaceutical industry, to prepare organic core–shell NCs and study their efficacy as foliar or root delivery vehicles. NCs ranging in diameter from 55 to 200 nm, with surface zeta potentials from -40 to +40 mV, and with seven different shell material properties were prepared and studied. Shell materials included synthetic polymers poly(acrylic acid), poly(ethylene glycol), and poly(2-(dimethylamino)ethyl methacrylate), naturally occurring compounds fish gelatin and soybean lecithin, and semisynthetic hydroxypropyl methylcellulose acetate succinate (HPMCAS). NC cores contained a gadolinium tracer for tracking by mass spectrometry, a fluorescent dye for tracking by confocal microscopy, and model hydrophobic compounds (alpha tocopherol acetate and polystyrene) that could be replaced by agrochemical payloads in subsequent applications. After foliar application onto tomato plants with Silwet L-77 surfactant, internalization efficiencies of up to 85% and NC translocation efficiencies of up to 32% were observed. Significant NC trafficking to the stem and roots suggests a high degree of phloem loading for some of these formulations. Results were corroborated by confocal microscopy and synchrotron X-ray fluorescence mapping. NCs stabilized by cellulosic HPMCAS exhibited the highest degree of translocation, followed by formulations with a significant surface charge. The results from this work indicate that biocompatible materials like HPMCAS are promising agrochemical delivery vehicles in an industrially viable pharmaceutical nanoformulation process (FNP) and shed light on the optimal properties of organic NCs for efficient foliar uptake, translocation, and delivery.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calculation of the First Moment of Energy Using D-T Reactivity Formalisms Under the Maxwell-Boltzmann Distribution (Pt. 1)

Nuclear fusion science is an example of a scientific field with a rich history of expert involvement and scientific publications, which together, form an expert-knowledge base. One example of this history is the utilization of published reaction rates from a variety of authors. Investigators for Deuterium-Tritium (D-T) ion fusion can choose from using frequently cited methods: the Bosch and Hale reactivity, thermonuclear reaction rates from Caughlan and Fowler, and the reactivity evaluation from Miley, Towner & Ivich which forms the basis of the Naval Research Lab (NRL) formulary. There are other choices available. Each of the reactivity formulations considered here, are based upon the Maxwell- Boltzmann velocity distribution for D-T fusion ion reactants. Numerical methods for computer codes simulating hot, energetic plasmas, include tabulations of the reactivity, and the first moment of energy. The purpose of this study is to answer the question: what is the first moment of energy, and how has it been formulated? This report is part of a series of reports by the authors on D-T fusion-reaction formalism. The present focus is on defining the mathematical relationship of the first moment of D-T fusion ion kinetic energy, $\langle$E$\rangle$, with the fusion crosssection, fusion reactivity and its derivative with respect to ion-temperature. Three variants of the first moment $\langle$E$\rangle$ are analytically developed and explored: 1) constant cross-section, 2) a normalized first moment, and 3) a particular function of the first moment. Much of the mathematical development is relegated to the Appendices, for a concise presentation of topics and results. A research component of this discussion is our comparison of some of the reactivity and its derivative results between evaluations (1) versus (2) and (3), which the authors have not found in any other publication to date. This investigation is useful for validation, verification (V&V) and uncertainty quantification (UQ) as it compares the work of several authors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Optimal Electrification Using Renewable Energies: Microgrid Installation Model with Combined Mixture k-Means Clustering Algorithm, Mixed Integer Linear Programming, and Onsset Method

Optimal planning and design of microgrids are priorities in the electrification of off-grid areas. Indeed, in one of the Sustainable Development Goals (SDG 7), the UN recommends universal access to electricity for all at the lowest cost. Several optimization methods with different strategies have been proposed in the literature as ways to achieve this goal. This paper proposes a microgrid installation and planning model based on a combination of several techniques. The programming language Python 3.10 was used in conjunction with machine learning techniques such as unsupervised learning based on K-means clustering and deterministic optimization methods based on mixed linear programming. These methods were complemented by the open-source spatial method for optimal electrification planning: onsset. Four levels of study were carried out. The first level consisted of simulating the model obtained with a cluster, which is considered based on the elbow and k-means clustering method as a case study. The second level involved sizing the microgrid with a capacity of 40 kW and optimizing all the resources available on site. The example of the different resources in the Togo case was considered. At the third level, the work consisted of proposing an optimal connection model for the microgrid based on voltage stability constraints and considering, above all, the capacity limit of the source substation. Finally, the fourth level involved a planning study of electrification strategies based mainly on microgrids according to the study scenario. The results of the first level of study enabled us to obtain an optimal location for the centroid of the cluster under consideration, according to the different load positions of this cluster. Then, the results of the second level of study were used to highlight the optimal resources obtained and proposed by the optimization model formulated based on the various technology costs, such as investment, maintenance, and operating costs, which were based on the technical limits of the various technologies. In these results, solar systems account for 80% of the maximum load considered, compared to 7.5% for wind systems and 12.5% for battery systems. Next, an optimal microgrid connection model was proposed based on the constraints of a voltage stability limit estimated to be 10% of the maximum voltage drop. The results obtained for the third level of study enabled us to present selective results for load nodes in relation to the source station node. Finally, the last results made it possible to plan electrification using different network technologies and systems in the short and long term. The case study of Togo was taken into account. The various results obtained from the different techniques provide the necessary leads for a feasibility study for optimal electrification of off-grid areas using microgrid systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Fast-, Light-Cured Scintillating Plastic for 3D-Printing Applications

Additive manufacturing techniques enable a wide range of possibilities for novel radiation detectors spanning simple to highly complex geometries, multi-material composites, and metamaterials that are either impossible or cost prohibitive to produce using conventional methods. The present work identifies a set of promising formulations of photocurable scintillator resins capable of neutron-gamma pulse shape discrimination (PSD) to support the additive manufacturing of fast neutron detectors. The development of these resins utilizes a step-by-step, trial-and-error approach to identify different monomer and cross-linker combinations that meet the requirements for 3D printing followed by a 2-level factorial parameter study to optimize the radiation detection performance, including light yield, PSD, optical clarity, and hardness. The formulations resulted in hard, clear, PSD-capable plastic scintillators that were cured solid within 10 s using 405 nm light. The best-performing scintillator produced a light yield 83% of EJ-276 and a PSD figure of merit equaling 1.28 at 450–550 keVee.

Frandsen, Brian G. (ORCID:0000000319964320)↗
DOE OSTI2025

Extrusion‐Spheronization of Mock Energetic Materials

The primary method for producing plastic bonded explosive (PBX) granules, or “prills”, has remained relatively unchanged for 70 years despite the complex nature of the process. In this work, we demonstrate the feasibility of using an extrusion‐spheronization technique to produce prills for PBX applications. We begin by detailing an inert formulation with similar properties of PBXs and then demonstrate the extrusion‐spheronization processing of these materials. A study is then performed where the spheronization process time of the extruded materials is varied and the resultant prills are morphologically characterized. Further, these prills are then pressed into high‐density articles and subject to compression testing to elucidate trends in process, properties, and performance. It was found that, for our formulation (95 wt.% melamine/5 wt.% polymer binder) and process, a spheronization time of 60 s yielded relatively uniform particles that exhibited improved poured, tapped, and pressed densities. Mechanical strength did not have a strong trend with process time as all spheronized materials had similar peak compression stress at failure. After further optimization, extrusion‐spheronization may be a promising path for future PBX formulation.

extrusion spheronization↗
DOE OSTI2022

In Situ Studies of the Formation of MoP Catalysts and Their Structure under Reaction Conditions for Higher Alcohol Synthesis: The Role of Promoters and Mesoporous Supports

Three formulations of a molybdenum phosphide (MoP) catalyst system were characterized for the higher alcohol synthesis (HAS) reaction using in situ X-ray absorption spectroscopy (XAS) and X-ray diffraction (XRD), monitoring the chemical phase evolution during activation and under reaction conditions. Here, the in situ study herein provides important insights into the effect of the support and of K that lead to high performance in HAS for K-promoted MoP supported on carbon, as evidenced by previous studies (ethanol selectivity: 29%; conversion: 5%). During the activation process, XAS shows that the carbon-supported samples reduce and reach a highly crystalline state at a lower temperature than the SBA-15-supported sample, indicating a substantial difference in catalyst activation. After activation, the samples are introduced to relevant reaction conditions resulting in spectra fairly similar to one another. XRD results corroborate the difference in the degree of crystallinity of these samples, in alignment with the XAS, and reveal the formation of crystalline potassium pyrophosphate (K 4 P 2 O 7 ) during the activation period of the K-promoted samples. This K 4 P 2 O 7 phase remains present under reaction conditions. Taken together, these results provide insights into the roles played by the carbon support and K promotion, connecting activity to electronic and crystal structure.

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

A Numerical Modeling Framework for Flocculation and Cohesive Sediment Transport in the Wave Bottom Boundary Layer

Flocculation, a critical process in coastal and estuarine systems, plays a significant role in sediment transport, nutrient cycling, and ecological health. This study develops a cohesive sediment transport modeling framework tailored to the wave bottom boundary layer under dilute and equilibrium conditions, explicitly incorporating flocculation effects via a Population Balance Equation (PBE). Using Direct Numerical Simulation, six baseline cases, each with a distinct sediment concentration profile resulting from a constant settling velocity and critical erosion shear stress, are generated to drive the PBE flocculation model for given floc yield strength and stickiness. Results reveal that flocculation significantly influences sediment concentration profiles promoting three distinct stages, well‐mixed, transition to lutocline, and well‐developed lutocline. At low concentrations with well‐mixed profiles, cohesive floc properties are less significant, and turbulence is a main flocculation driver. In contrast, as concentration increases, cohesive floc properties become crucial, facilitating lutocline formation. Here, the analysis also highlights limitations of depth‐averaged settling velocity as a parameterization. It is suitable for well‐mixed and transitional profiles but fails in well‐developed lutoclines, where empirical formulations that explicitly incorporate turbulent shear rate and sediment concentration better capture variability. This study underscores the necessity of incorporating flocculation effects into sediment transport models to enhance predictions of sediment dynamics in wave bottom boundary layers.

Penaloza‐Giraldo, Jorge A. [Oak Ridge National Lab↗
DOE OSTI2024

HybridSystemsSimulations.jl - Solving the Merchant Collocated Facilities with JuMP

The development of new clean-generation technologies also leads to new plant-level architectures that combine several generation and storage assets behind the point of connection. These co-located generation resources (Hybrid Systems) primarily operate as merchant assets that employ automated market bidding models and internal Energy Management Systems (EMS) to comply with the operator's signals. Formulating an optimal bidding model requires embedding the EMS control model into the bidding algorithm, resulting in a bi-level optimization problem. In this presentation, we first showcase using JuMP to formulate and solve this problem effectively for multiple merchant systems and the bidding outcomes considering different model formulations. Second, the bidding outcomes are later integrated into a PowerSimulations.jl (also built with JuMP) simulation to study the system-level effects of the various merchant bidding and the interactions between market-clearing models and the embedded EMS model. We will showcase simulations conducted in the RTS system considering different levels of merchant hybrid systems participation. The presentation provides the following specific insights on JuMP usage: 1) the Formulation of specialized bi-level problems with custom cuts to solve the merchant hybrid system bidding problem; 2) the integration of a modular model within a complex simulation workflow supported by JuMP in PowerSimulations.jl; 3) Accelerating the solution of power systems operations simulation that employ agent optimization problems using JuMP.

energy markets↗
DOE OSTI2022

Development, verification, and validation of comprehensive acoustic fluid-structure interaction capabilities in an open-source computational platform

The acoustic fluid-structure interaction (FSI) formulation is a practical numerical approach for the seismic analysis of fluid-filled tanks. However, there are no verification and validation studies reported in the literature that demonstrate the ability of an acoustic FSI numerical model to predict responses important to structural and mechanical design for intense translational and rotational earthquake inputs. Herein, an acoustic FSI formulation is implemented in the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE), and is formally verified and validated using analytical solutions and code-to-code verification, and experimental data, respectively. The analytical solutions are for small amplitude, unidirectional seismic inputs. The code-to-code verification utilizes a previously verified and validated Arbitrary Lagrangian-Eulerian (ALE) numerical model in the commercial finite element code LS-DYNA. The validation studies utilize a comprehensive data set assembled from results of 3D earthquake-simulator tests of a fluid-filled vessel. The acoustic numerical model in MOOSE is verified and validated for hydrodynamic pressures and support reactions except for cases that involve significant convective response. For small amplitude inputs, numerically predicted wave heights match those of the analytical solutions. The numerical model is not verified and validated for wave height calculations under intense 3D seismic inputs. The run times for the acoustic FSI simulations in MOOSE are an order of magnitude, or more, shorter than for the corresponding ALE simulations in LS-DYNA. The utility of the MOOSE acoustic FSI implementation is demonstrated by seismic analysis of a building equipped with a fluid-filled, advanced nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Joint Management and Optimization of Residential Natural Gas and Electricity Distribution Networks Coupled via Fuel Cells

The attractive features of natural gas as well as the growing electric power demand worldwide have created increasing interest in natural-gas-based distributed generation applications for electric distribution networks. Here, this paper investigates the interdependency between a residential natural gas network and an electric distribution network that are linked together via fuel cells. The modeling of the natural gas network is introduced first, and then the algorithm for gas flow study is presented. The optimal placement and sizing of fuel cell based distributed generation systems are formulated to minimize the losses in both the natural gas network and the electric distribution grid, subject to the constraints imposed by both networks. In addition, a probabilistic model for both gas and electricity demands is developed based on historical electricity and natural gas demand data. A K-means clustering method is used to determine the hourly load states to solve the joint probabilistic optimization problem. Simulation studies are carried out on an integrated system consisting of the IEEE 69-bus distribution network and a radial 27-node natural gas network to verify the developed optimization model and the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Intricacies of frustrated magnetism in the Kondo metal YbAgGe

The combination of localized magnetic moments, their frustration and interaction with itinerant electrons is a key challenge of condensed matter physics. Frustrated magnetic interactions promote degenerate ground states with enhanced fluctuations, a topic that is predominantly studied in magnetic insulators. The coupling between itinerant and localized electrons in metals add complexity to the problem, and is presently formulated only for extreme cases in which the itinerant electrons mediate exchange between localized spins (RKKY interaction) or suppress the formation of magnetic moments (Kondo screening). Here, we report an in-depth experimental study of the distorted Kagome metal YbAgGe, unravelling the open questions of how frustration, localized magnetism and itinerant electrons are intertwined in frustrated Kondo metals. We find that coupled itinerant and localized electrons give rise to dynamic magnetic correlations below T* ≈ 20 K. At lower temperature, frustrated magnetic interactions establish anisotropic magnetic short-range correlations that culminate into antiferromagnetic long-range order below T N = 0.68 K with a significantly reduced modulated magnetic moment. We show that local moment Hamiltonians can yield limited understanding of the microscopic behaviour in frustrated metals, and prompt the extension of more sophisticated model Hamiltonians incorporating itinerant effects.

Mazzone, Daniel G. [PSI Center for Neutron and Muo↗
DOE OSTI2025

Mixed-Integer Linear Programming Formulation with Embedded Machine Learning Surrogates for the Design of Chemical Process Families

In previous work, we introduced process family design. The main idea is to design a platform of common elements, and, allowing us to capture additional cost savings, simultaneously design a family of processes, and reducing both engineering and deployment timelines. We formulate this as an optimization problem, specifically a nonlinear generalized disjunctive program (GDP). We have proposed two approaches for reformulating and solving this problem: one based on full-discretization of the design space and one that uses Machine Learning (ML) surrogates to replace the nonlinear process models. Using ML surrogates to predict required system costs and performance indicators allows us to reformulate the nonlinearities in the GDP generate an efficient MILP formulation. In this work, we apply the ML surrogate approach to two case studies. One case study involves designing a family of carbon capture systems to cover a set of different flue gas flow rates and inlet CO 2 concentrations, where we consider the absorber and stripper as common unit module types. The second case study focuses on a water-desalination process, where we design a family of these processes for a variety of salt concentrations and flow rates. In both of these case studies, we demonstrate a scalable optimization approach that enables the design of multiple processes simultaneously, reducing the time-to-market and overall costs by maximizing the cost savings due to both economies of scale and economies of numbers.

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

Optimizing Desalination Operations for Energy Flexibility

Despite the value of energy optimization in desalination processes, modeling dynamic operations for monthly billing periods has remained a computational challenge. This work proposes a framework for energy flexibility optimization, which includes new modeling features for independent operation of parallel skids, start-up delays associated with chemical stabilization, the consideration of industrial energy tariff structures, and inclusion of hourly electrical carbon intensities. This is done using a modular and computationally efficient formulation that guarantees a globally optimal solution with standard optimization solvers. In this study, the approach is demonstrated in two distinct case studies: a seawater desalination plant in Santa Barbara, CA, and an indirect potable reuse facility in San Jose, CA. Trends predicted from the model are validated against operational facility measurements from a demand response shutdown event. Preliminary results show that optimizing energy flexibility can result in 18.51% monthly cost savings over energy efficiency-optimized operation. The value extracted from a facility-wide shutdown during peak electricity price hours is hampered by start-up delays in post-treatment chemical stabilization. In cases in which a facility does not have much excess capacity, using a flow equalization tank or operating over a wide recovery range may be cost-effective.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗