Thermal separation processes, such as distillation, play a pivotal role in the chemical and petrochemical sectors, constituting a substantial portion of the industrial energy consumption. Consequently, owing to their huge application scales, these processes contribute significantly to greenhouse gas (GHG) emissions. Decarbonizing distillation units could mitigate carbon emissions substantially. Heat Pumps (HP), that recycle lower quality heat from the condenser to the reboiler by electric work present a unique opportunity to electrify distillation systems. In this research we try to answer the following question in the context of multi-component distillation – Do HPs actually reduce the effective fuel consumption or just merely shift the fuel demand from chemical industry to the power plant? If they do, what strategies consume minimum energy? To address these inquiries, we construct various simplified surrogate and shortcut models designed to effectively encapsulate the fundamental phy... [more]

Modern power grids coordinate electricity production and consumption via multi-scale wholesale energy markets. Historically, levelized cost metrics were the de facto standard for techno-economic analyses of energy systems and comparison of technology options. However, these metrics neglect the complexity of energy infrastructure including the time-varying value of electricity. An emerging alternative is multi-period optimization, which considers the locational marginal price of electricity as input data (parameters). In this work, we present a general interface for multi-period optimization with time-varying energy prices to facilitate rapid analysis and comparison of potential energy systems models. The PriceTakerModel class is written in the IDAES®-PSE platform and allows users to generate a multi-period, price-taker model instance, as well as automatically generate common operational constraints for their model, such as start-up and shutdown. We show this interface successfully gene... [more]

This study presents a comprehensive approach to optimizing hydrogen supply chain network (HSCN), focusing initially on Texas, with potential scalability to national and global regions. Utilizing mixed-integer nonlinear programming (MINLP), the research decomposes into two distinct modeling stages: broad supply chain modeling and detailed hub-specific analysis. The first stage identifies optimal hydrogen hub locations, considering county-level hydrogen demand, renewable energy availability, and grid capacity. It determines the number and placement of hubs, county participation within these hubs, and the optimal sites for hydrogen production plants. The second stage delves into each selected hub, analyzing energy mixes under variable solar, wind, and grid profiles, sizing specific production and storage facilities, and scheduling to match energy availability. Iterative refinement incorporates detailed insights back into the broader model, updating costs and configurations to converge upo... [more]

This study focuses on optimizing solid oxide electrolysis cell (SOEC) systems for efficient and durable long-term hydrogen (H2) production. While the elevated operating temperatures of SOECs offer advantages in terms of efficiency, they also lead to chemical degradation, which shortens cell lifespan. To address this challenge, dynamic degradation models are coupled with a steady-state, two-dimensional, non-isothermal SOEC model and steady-state auxiliary balance of plant equipment models, within the IDAES modeling and optimization framework. A quasi-steady state approach is presented to reduce model size and computational complexity. Long-term dynamic simulations at constant H2 production rate illustrate the thermal effects of chemical degradation. Dynamic optimization is used to minimize the lifetime cost of H2 production, accounting for SOEC replacement, operating, and energy expenses. Several optimized operating profiles are compared by calculating the Levelized Cost of Hydrogen (LC... [more]

Biogas is a promising energy source for sustainable hydrogen production due to its high concentration of CH4. However, determining the optimal process configuration is challenging due to the uncertainty of the fed biogas composition and the sensitivity of the operating conditions. This necessitates early-stage evaluation of the biomass-to-hydrogen process's performance, considering economics, energy efficiency, and environmental impacts. A data-driven model was introduced for early-stage assessment of hydrogen production from biogas without whole process simulation and optimization. The model was developed based on various biogas compositions and generated parameters for mass and energy balance. A database of unit processes was created using simulation models. Sensitivity analysis was performed under four techno-economic and environmental evaluation criteria: Unit Production Cost (UPC), Energy Efficiency (EEF), Net CO2 equivalent Emission (NCE), and Maximum H2 Production (MHP). The ea... [more]

Most integrated energy system (IES) optimization frameworks employ the price-taker approximation, which ignores important interactions with the market and can result in overestimated economic values. In this work, we propose a machine learning surrogate-assisted optimization framework to quantify IES/market interactions and thus go beyond price-taker. We use time series clustering to generate representative IES operation profiles for the optimization problem and use machine learning surrogate models to predict the IES/market interaction. We quantify the accuracy of the time series clustering and surrogate models in a case study to optimally retrofit a nuclear power plant with a polymer electrolyte membrane electrolyzer to co-produce electricity and hydrogen.

The chemical industry is actively pursuing energy transition and decarbonization through renewables and other decarbonization initiatives. However, navigating this transition is challenging due to uncertainties in capital investments, electricity costs, and carbon taxes. Adapting to decarbonization standards while preserving existing valuable infrastructure presents a dilemma. Early transitions may lead to inefficiencies, while delays increase the carbon footprint. This research proposes a framework to find an optimal retrofit decarbonization strategy for existing oil refineries. We start with a generic process flowsheet representing the refinery's current configuration and operations, and consider various decarbonization alternatives. Through superstructure optimization, we identify the most cost-effective retrofit strategy over the next three decades to achieve decarbonization goals. We develop a Mixed-Integer Linear Programming (MILP) model, integrating simplified process equations... [more]

In this work, we present a follow-up work of reinforcement learning (RL)-driven process design using the Institute for Design of Advanced Energy Systems Process Systems Engineering (IDAES-PSE) Framework. Herein, process designs are generated as stream inlet-outlet matrices and optimized using the IDAES platform, the objective function value of which is the reward to RL agent. Deep Q-Network is employed as the RL agent including a series of convolutional neural network layers and fully connected layers to compute the actions of adding or removing any stream connections, thus creating a new process design. The process design is then informed back to the RL agent to refine its learning. The iteration continues until the maximum number of steps is reached with feasible process designs generated. To further expedite the RL search of the design space which can comprise the selection of any candidate unit(s) with arbitrary stream connections, we investigate the role of RL reward function and... [more]

Continuous manufacturing in pharmaceutical industries has shown great promise to achieve process intensification. To better understand and justify such changes to the current status quo, a technoeconomic analysis of a continuous production must be conducted to serve as a predictive decision-making tool for manufacturers. This paper uses PharmaPy, a custom-made Python-based library developed for pharmaceutical flowsheet analysis, to simulate an annual production cycle for a given active pharmaceutical ingredient (API) of varying production volumes for a batch crystallization system and a continuous mixed suspension, mixed product removal (MSMPR) crystallizer. After each system is optimized, the generalized cost drivers, categorized as capital expenses (CAPEX) or operational expenses (OPEX), are compared. Then, a technoeconomic and sustainability cost analysis is done with the process mass intensity (PMI) as a green metric. The results indicate that while the batch system does have an ov... [more]

Computer-aided Molecular and Process Design (CAMPD) is an equation-oriented multi-scale decision making framework for designing both materials (molecules) and processes for separation, reaction, and reactive separation whenever material choice significantly impacts process performance. The inherent nonlinearity and nonconvexity in CAMPD optimization models, introduced through the property and process models, pose challenges to state-of-the-art solvers. Recently, quantum computing (QC) has shown promise for solving complex optimization problems, especially those involving discrete decisions. This motivates us to explore the potential usage of quantum optimization techniques for solving CAMPD problems. We have developed a technique for directly solving a class of mixed integer nonlinear programs using QC. Our approach represents both continuous and integer design decisions by a set of binary variables through encoding schemes. This transformation allows to reformulate certain types of CA... [more]

Antibody therapies are important in treating life-threatening ailments such as cancer and autoimmune diseases. Purity of the antibody is essential for successful applications and Protein A selective resin extraction is the standard step for antibody recovery. Unfortunately, such resins can cost up to 30% of the total cost of antibody production. Hence, the optimal design of this purification step becomes a critical factor in downstream processing to minimize the size of the column needed. An accurate predictive model, as a digital twin representing the purification process, is necessary where changes in the flow rates and the inlet concentrations are modeled via the Method of Moments. The system uncertainties are captured by including the stochastic Ito process model of Brownian motion with drift. Pontryagin’s Maximum Principle under uncertainty is then applied to predict the flowrate control strategy for optimized resin use, column design, and efficient capturing of the antibodies. In... [more]

In this paper, a hybrid heuristic rule-based and deterministic optimization-driven process decision framework is presented for the analysis and optimization of process flowsheets for end-to-end optimal (E2E0) pharmaceutical manufacturing. The framework accommodates various operating modes, such as batch, semi-batch and continuous, for the different unit operations that implement each manufacturing step. To address the challenges associated with solving process synthesis problems using a simulation-optimization approach, heuristic-based process synthesis rules are employed to facilitate the reduction of the superstructure into smaller sub-structures that can be more readily optimized. The practical application of the framework is demonstrated through a case study involving the end-to-end continuous manufacturing of an anti-cancer drug, lomustine. Alternative flowsheet structures are evaluated in terms of the sustainability metric, E-factor while ensuring compliance with the required pro... [more]

Explicitly incorporating the effects of chemical phenomena such as chemical pretreatment and mineral scaling during the design of treatment systems is critical; however, the complexity of these phenomena and limitations on data have historically hindered the incorporation of detailed water chemistry into the modeling and optimization of water desalination systems. Thus, while qualitative assessments and experimental studies on chemical pretreatment and scaling are abundant in the literature, very little has been done to assess the technoeconomic implications of different chemical pretreatment alternatives within the context of end-to-end water treatment train optimization. In this work, we begin to address this challenge by exploring the impact of pH control during pretreatment on the cost and operation of a high-recovery desalination train. We compare three pH control methods used in water treatment (H2SO4, HCl, and CO2) and assess their impact on the operation of a desalination plant... [more]

This paper proposes an AI-assisted approach aimed at accelerating chemical process design through causal incremental reinforcement learning (CIRL) where an intelligent agent is interacting iteratively with a process simulation environment (e.g., Aspen HYSYS, DWSIM, etc.). The proposed approach is based on an incremental learnable optimizer capable of guiding multi-objective optimization towards optimal design variable configurations, depending on several factors including the problem complexity, selected RL algorithm and hyperparameters tuning. One advantage of this approach is that the agent-simulator interaction significantly reduces the vast search space of design variables, leading to an accelerated and optimized design process. This is a generic causal approach that enables the exploration of new process configurations and provides actionable insights to designers to improve not only the process design but also the design process across various applications. The approach was valid... [more]

Design problems for process systems engineering applications often require multi-scale modeling integrating detailed process models. Consequently, black-box optimization and surrogate modeling have continued to play a fundamental role in mission-critical design applications. Inherent in surrogate modeling applications, particularly those constrained by “expensive” function evaluations, are the questions of how to properly balance “exploration” and “exploitation” and how to do so while harnessing parallel computing in techniques. We devise and investigate a one-step look-ahead GRASP heuristic for balancing exploration and exploitation in a parallel environment. Computational results reveal that our approach can yield equivalent or superior surrogate quality with near linear scaling in the number of parallel samples.

Chemical recycling of plastics is a promising technology to reduce carbon footprint and ease the pressure of waste treatment. Specifically, highly efficient conversion technologies for polyolefins will be the most effective solution to address the plastic waste crisis, given that polyolefins are the primary contributors to global plastic production. Significant challenges encountered by plastic waste valorization facilities include the uncertainty in the composition of the waste feedstock, process yield, and product price. These variabilities can lead to compromised performance or even render operations infeasible. To address these challenges, this work applied the robust optimization-based framework to design an integrated polyolefin chemical recycling plant. Data-driven surrogate model was built to capture the separation units’ behavior and reduce the computational complexity of the optimization problem. It was found that when process yield and price uncertainties were considered, wa... [more]

Solvent-based carbon capture processes typically suffer from the temperature rise of the solvent due to the heat of absorption of CO2. This increased temperature is not thermodynamically favorable and results in a significant reduction in performance in the absorber column. As opposed to interstage coolers, which only remove, cool, and return the solvent at discrete locations in the column, internal coolers that are integrated with the packing can cool the process inline, which can result in improved efficiency. This work presents the modeling of these internal coolers within an existing generic, equation-oriented absorber column model that can cool the process while allowing for simultaneous mass transfer. Optimization of this model is also performed, which is capable of optimally choosing the best locations to place these devices, such that heat removal and mass transfer area are balanced. Results of the optimization have shown that optimally placed cooling elements result in a signi... [more]

Cyclic adsorption processes attain a cyclic-steady state (CSS) condition by undergoing repeated cycles in time, owing to their transient and modular nature. Mathematically, solving a set of underlying nonlinear partial differential equations iteratively for different steps in a cycle until the CSS condition is attained presents a computational challenge, making the simulation and optimization of cyclic adsorption processes time-consuming. This paper focuses on expediting the CSS convergence in adsorption process simulations by implementing two vector-based acceleration methods that offer quadratic convergence akin to Newton’s methods. These methods are straightforward to implement, requiring no prior knowledge of the first derivatives (or Jacobian). The study demonstrates the efficacy of accelerated convergence by considering two adsorption processes that exhibit complex dynamics, namely, a four-step vacuum swing adsorption and a six-step temperature swing adsorption cycles for post-co... [more]

Developing methods for rapid, large-scale deployment of carbon capture systems is critical for meeting climate change goals. Optimization-based decisions can be employed at the design and manufacturing phases to minimize the costs of deployment and operation. Manufacturing standardization results in significant cost savings due to economies of numbers. Building on previous work, we present a process family design approach to design a set of carbon capture systems while explicitly including economies of numbers savings within the formulation. Our formulation optimizes both the number and characteristics of the common components in the platform and simultaneously designs the resulting set of carbon capture systems. Savings from economies of numbers are explicitly included in the formulation to determine the number of components in the platform. We show and discuss the savings we gain from economies of numbers.

In this work, we present recent algorithmic and implementation advances of the nonconvex two-stage robust optimization solver PyROS. Our advances include extensions of the scope of PyROS to models with uncertain variable bounds, improvements to the formulations and/or initializations of the various subproblems used by the underlying cutting set algorithm, and extensions to the pre-implemented uncertainty set interfaces. The effectiveness of PyROS is demonstrated through the results of an original benchmarking study on a library of over 8,500 small-scale instances, with variations in the nonlinearities, degree-of-freedom partitioning, uncertainty sets, and polynomial decision rule approximations. To demonstrate the utility of PyROS for large-scale process models, we present the results of a carbon capture case study. Overall, our results highlight the effectiveness of PyROS for obtaining robust solutions to optimization problems with uncertain equality constraints.

The determination of the design space (DS) in a pharmaceutical process is a crucial aspect of the quality-by-design (QbD) initiative which promotes quality built into the desired product. This is achieved through a deep understanding of how the critical quality attributes (CQAs) and process parameters (CPPs) interact that have been demonstrated to provide quality assurance. For computational inexpensive models, the original process model can be directly deployed to identify the design space. One such crucial process is the Tablet Press (TP), which directly compresses the powder blend into individual units of the final product or adds dry or wet granulation to meet specific formulation needs. In this work, we identify the design space of input variables in a TP such that there is a (probabilistic) guarantee that the tablets meet the quality constraints under a set of operating conditions. A reduced-order model of TP is assigned for this purpose where the effects of lubricants and glidan... [more]

Superstructure optimization for process synthesis is a challenging endeavour typically leading to large scale MINLP formulations. By the combination of phenomena-based building blocks, accurate thermodynamics, and structural screening we obtain a new framework for optimal process synthesis, which overcomes prior limitations regarding solution by deterministic MINLP solvers in combination with accurate thermodynamics. This is facilitated by MOSAICmodeling’s generic formulation of models in MathML / XML and subsequent decomposition and code export to GAMS and C++. A branch & bound algorithm is implemented to solve the overall MINLP problem, wherein the structural screening penalizes instances, which are deemed nonsensical and should not be further pursued. The general capabilities of this approach are shown for the distillation-based separation of a ternary system.

We present new modifications to superstructure optimization paradigms to i) enable their robust solution and ii) extend their applicability. Superstructure optimization of chemical process flowsheets on the basis of rigorous and detailed models of the various unit operations, such as in the state operator network (SON) paradigm, is prone to non-convergence. A key challenge in this optimization-based approach is that when process units are deselected from a superstructure flowsheet, the constraints that represent the deselected process unit can be numerically singular (e.g., divide by zero, logarithm of zero and rank-deficient Jacobian). In this paper, we build upon the recently-proposed modified state operator network (MSON) that systematically eliminates singularities due to unit deselection and is equally applicable to the context of both simulation-based and equation-oriented optimization. A key drawback of the MSON is that it is only applicable to the design of isobaric flowsheets... [more]

Maintaining product integrity in multi-product oil pipelines is crucial for efficiency and profit. This study presents a strategy combining design and process improvement to enhance flushing protocols, addressing the challenge of residual batch contamination. A pilot plant, mirroring industrial operations through dimensionless residence time distribution, was developed to identify and rectify bottlenecks during product transition. The pilot plant’s success in replicating industrial operations paves the way for targeted experiments and modelling to enhance optimized flushing, ensuring product quality and operational excellence.

This paper presents an algorithm for developing sparse surrogate models that satisfy mass/energy conservation even when the training data are noisy and violate the conservation laws. In the first step, we employ the Bayesian Identification of Dynamic Sparse Algebraic Model (BIDSAM) algorithm proposed in our previous work to obtain a set of hierarchically ranked sparse models which approximate system behaviors with linear combinations of a set of well-defined basis functions. Although the model building algorithm was shown to be robust to noisy data, conservation laws may not be satisfied by the surrogate models. In this work we propose an algorithm that augments a data reconciliation step with the BIDSAM model for satisfaction of conservation laws. This method relies only on known boundary conditions and hence is generic for any chemical system. Two case studies are considered-one focused on mass conservation and another on energy conservation. Results show that models with minimum bia... [more]