Soybean is a highly valuable global commodity due to its versatility and numerous derivative products. During harvest, all non-seed materials become “straw”. Currently, this waste is primarily used for low-value purposes such as animal feed, landfilling, and incineration. To address this, the present work proposes a conceptual biorefinery aimed at converting soybean straw into higher-value products. The study began with data collection to identify potential conversion routes. Based on this information, a superstructure was developed, comprising seven conversion routes: four thermochemical routes (pyrolysis, combustion, hydrothermal gasification, and liquefaction), two biological routes (fermentation and anaerobic fermentation), and one chemical route (alkaline extraction). Each process was evaluated based on product yields, conversion times, and associated capital and operating costs. Using this data, an MILP (Mixed-Integer Linear Programming) optimization model was built in Pyomo usin... [more]

Renewable energy and energy efficiency are increasingly recognised as crucial for creating new economic opportunities and mitigating environmental impacts. Anaerobic digestion (AD) transforms organic materials into a clean, renewable energy source. Co-digestion of various organic wastes and energy crops addresses the disadvantages of single-substrate digestion, increasing production flexibility yet adding process complexity and sensitivity. This study employs a two-pronged approach to optimise biogas production while considering global warming potential: a nonlinear programming (NLP) model for dynamic system economic optimisation with a model predictive control (MPC) strategy for precise temperature regulation within the digester. The NLP model integrates a combined heat and power (CHP) system to leverage dynamic electricity, heat, and gas prices, accounting for physical and economic parameters such as biomethane potential, chemical oxygen demand, and substrate density. A cardinal temp... [more]

Process control can be one of the most exciting and engaging chemical engineering undergraduate courses! This paper describes our experience transforming Chemical Process Control into Data Analytics, Optimization, and Control at the University of Notre Dame (second semester required course in the junior year). Our modern course is built around six hands-on experiments in which students practice data-centric modeling and analysis using the Arduino-based Temperature Control Lab (TCLab) hardware. We argue that state-space dynamic modeling and optimization are more critical for educating modern chemical engineers than topics such as frequency domain analysis and controller synthesis emphasized in many classical undergraduate control courses. All the course material is available online at https://ndcbe.github.io/controls.

Robust optimisation is a powerful approach for addressing uncertainty ensuring constraint satisfaction for all uncertain parameter realisations. While convex robust optimisation problems are effectively tackled using robust reformulations and cutting plane methods, extending these techniques to non-convex problems remains largely unexplored. In this work we propose a method that is based on a parallel robustness and optimality search. We introduce a novel spatial Branch-and-Bound algorithm integrated with robust cutting-planes for solving non-convex robust optimisation problems. The algorithm systematically incorporates global and robust optimisation techniques, leveraging McCormick relaxations. The proposed algorithm is evaluated on benchmark pooling problems with uncertain feed quality, demonstrating algorithm stability and solution robustness. The computational time for the examined case studies is within the same order of magnitude as state-of-the-art. The findings of this work hig... [more]

Process Systems Engineering (PSE) provides the advanced conceptual framework and software tools to formulate and optimise well-considered integrated solutions that could accelerate the sustainability transition within the industrial sector. The landscape of advanced PSE is poised to undertake a considerable transformation with the rise in popularity of open-source and script-based software platforms with predictive modelling capabilities based on modern mathematical optimization techniques. This paper highlights three leading equation-based platforms—IDAES, Modelica, and GEKKO-that are increasingly utilised for the modelling, simulation, and optimisation of complex systems within the advanced PSE domain, alongside the strengths and limitations of each approach. Following this, we present a framework through which emerging techniques within the domain of Software Engineering could be leveraged to address these limitations, with a vision of improving the accessibility and flexibility of... [more]

Dynamic or Programmable Catalysis is an innovative strategy to improve heterogeneous catalysis processes by modulating the binding energies (BE) of adsorbates on a catalytic surface. The technique enables the periodic favoring of different reaction steps, overcoming limitations imposed by the Sabatier Principle and allowing for higher overall reaction rates, otherwise unattainable. Previously, we implemented a simultaneous simulation approach using the algebraic modeling language Pyomo and the solver IPOPT to obtain cyclic steady state results for a unimolecular reactive system with up to four-order of magnitude increases in computational performance compared to the previously reported sequential approach. The flexibility of the method allowed for the investigation of the influence of forcing signal parameters on system behavior and provided a framework for waveform design. In this study, we use a hybrid framework that combines the sequential and the simultaneous simulation approaches... [more]

Dynamic modelling of wind turbines and their simulation is a very useful tool for studying their behaviour. One of the key elements concerning the physical models of wind turbines is the power coefficient Cp, which acts as an efficiency in the extraction of power from the wind. Unfortunately, this coefficient is often unknown a priori, as it does not usually appear in the information provided by manufacturers. This paper first describes a methodology for obtaining the power coefficient parameters of a commercial wind turbine model using the power curve provided by the manufacturer, which indicates the theoretical power that the wind turbine can produce at each wind speed. To achieve this, a parameter estimation problem is formulated and solved to determine the power coefficient parameters. Nevertheless, this information is often insufficient, requiring additional knowledge, such as operational data, to improve the fit. Finally, a new parameter estimation is performed using only real da... [more]

Soybean is a highly valuable global commodity due to its versatility and numerous derivative products. During harvest, all non-seed materials become “straw”. Currently, this waste is pri-marily used for low-value purposes such as animal feed, landfilling, and incineration. To address this, the present work proposes a conceptual biorefinery aimed at converting soybean straw into higher-value products. The study began with data collection to identify potential conversion routes. Based on this information, a superstructure was developed, comprising seven conversion routes: four thermochemical routes (pyrolysis, combustion, hydrothermal gasification, and lique-faction), two biological routes (fermentation and anaerobic fermentation), and one chemical route (alkaline extraction). Each process was evaluated based on product yields, conversion times, and associated capital and operating costs. Using this data, an MILP (Mixed-Integer Linear Program-ming) optimization model was built in Pyomo u... [more]

Equation-Oriented (EO) modeling techniques have been gaining popularity as an alternative for simulating and optimizing process systems due to their flexibility and ability to leverage state-of-the-art solvers inaccessible to many procedural modeling approaches. Despite these advantages, adopting EO modeling tools remains challenging due to the significant learning curve and effort required to build and solve models. Many techniques are available to help diagnose problems with EO process models and reduce the effort required to create and use them. However, these techniques still need to be integrated into EO modeling environments, and many modelers are unaware of sophisticated EO diagnostic tools. To survey the availability of model diagnostic tools and common workflows, the U.S. Department of Energy’s Institute for the Design of Advanced Energy Systems (IDAES) has conducted user experience interviews of users of the IDAES Integrated Platform (IDAES-IP) for process modeling. The inter... [more]

Equation-oriented (EO) modeling has the potential to enable the effective design and optimization of the operation of advanced energy systems. However, advanced modeling of energy systems results in a large number of variables and non-linear equations, and it can be difficult to search through these to identify the culprit(s) responsible for convergence issues. The Institute for the Design of Advanced Energy Systems Integrated Platform (IDAES-IP) contains a tool to identify poorly scaled constraints and variables by searching for rows and columns of the Jacobian matrix with small L2-norms so they can be rescaled. A further singular value decomposition can be performed to identify degenerate sets of equations and remaining scaling issues. This work presents an EO model of a flowsheet developed for post-combustion carbon capture using a monoethanolamine (MEA) solvent system as a case study. The IDAES diagnostics tools were successfully applied to this flowsheet to identify problems to im... [more]

Reverse electrodialysis (RED) is a nascent renewable technology that generates clean, baseload electricity from salinity differences between two water streams, a renewable source known as salinity gradient energy (SGE). Full-scale RED progress calls for robust techno-economic and environmental assessments. Using generalized disjunctive programming (GDP) and life cycle assessment (LCA) principles, this work proposes cost-optimal and sustainable RED process designs involving different RED stack sizes and width-over-length ratios to guide the design and operation from the demonstration to full-scale phases. Results indicate that RED units will benefit from larger aspect ratios with a relative increase in net power of over 30% with 6 m2 membrane size. Commercial RED unit sizes (0.25–3 m2) require larger aspect ratios to reach an equal relative increase in net power but exhibit higher power densities. The GDP model devises profitable RED process designs for all the assessed aspect ratios in... [more]

The oil and gas energy sector uses billions of gallons of water for hydraulic fracturing each year to extract oil and gas. The water injected into the ground for fracturing along with naturally occurring formation water from the oil wells surfaces back in the form of produced water. Produced water can contain high concentrations of total dissolved solids and is unfit for reuse outside the oil and gas industry without desalination. In semi-arid shale plays, produced water desalination for beneficial reuse could play a crucial role in alleviating water shortages and addressing extreme drought conditions. In this paper we co-optimize the design and operation of desalination technologies along with operational decisions across produced water networks. A multi-period produced water network model with simplified split-fraction-based desalination nodes is developed. Rigorous steady-state desalination mathematical models based on mechanical vapor recompression are developed and embedded at the... [more]

Post-combustion carbon capture technologies have the potential to contribute significantly to achieving the environmental goals of reducing CO2 emissions in the short term. However, these technologies are energy and cost-intensive, and the variability of flue gas represents important challenges. The optimal design and optimization of such systems are critical to reaching the net zero and net negative goals, in this context, the use of computer-aided process design can be very effective in overcoming these issues. In this study, we explore the implementation of carbon capture technologies within an industrial complex, by considering the pooling of CO2 streams. We present an optimization formulation to design carbon capture plants with the goal of enhancing efficiency and minimizing the capture costs. Capital and operating costs are represented via surrogate models (SMs) that are trained using rigorous process models in Aspen Plus, each data point is obtained by solving an optimization p... [more]

Adaptive and self-optimizing intelligent systems such as digital twins are increasingly important in science and engineering. Digital twins utilize mathematical models to provide added precision to decision-making. However, physics-informed models are challenging to build, calibrate, and validate with existing data science methods. Model-based design of experiments (MBDoE) is a popular framework for optimizing data collection to maximize parameter precision in mathematical models and digital twins. In this work, we apply MBDoE, facilitated by the open-source package Pyomo.DoE, to train and validate mathematical models for batch crystallization. We quantitatively examined the estimability of the model parameters for experiments with different cooling rates. This analysis provides a quantitative explanation for the heuristic of using multiple experiments at different cooling rates.

Water chemistry plays a critical role in the design and operation of water treatment processes. Detailed chemistry modeling tools use a combination of advanced thermodynamic models and extensive databases to predict phase equilibria and reaction phenomena. The complexity and formulation of these models preclude their direct integration in equation-oriented modeling platforms, making it difficult to use their capabilities for rigorous water treatment process optimization. Neural networks (NN) can provide a pathway for integrating the predictive capability of chemistry software into equation-oriented models and enable optimization of complex water treatment processes across a broad range of conditions and process designs. Herein, we assess how NN architecture and training data impact their accuracy and use in equation-oriented water treatment models. We generate training data using PhreeqC software and determine how data generation and sample size impact the accuracy of trained NNs. The... [more]

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 main purpose of the work was examining various methods of decomposition of a network optimization problem of simultaneous routing and bandwidth allocation based on Lagrangian relaxation. The problem studied is an NP-hard mixed-integer nonlinear optimization problem. Multiple formulations of the optimization problem are proposed for the problem decomposition. The decomposition methods used several problem formulations and different choices of the dualized constraints. A simple gradient coordination algorithm, cutting-plane coordination algorithm, and their more sophisticated variants were used to solve dual problems. The performance of the proposed decomposition methods was compared to the commercial solver CPLEX and a heuristic algorithm.