Bayesian optimization is known to be able to search for the optimal conditions based on a small number of experiments. However, these experiments are insufficient to understand the experimental condition space. In contrast, we report the development of an algorithm that combines a low-confounding definitive screening design with Bayesian optimization, allowing for rapid optimization and ensuring sufficient experiments to understand the experimental condition space with a low confounding.

In process systems engineering (PSE), obtaining accurate process models for optimization can be expensive and time-consuming. Black-box Bayesian Optimization (BO) with Gaussian process (GP) surrogates offers a promising approach. However, full black-box optimization neglects valuable prior knowledge, which could otherwise improve the optimization process. This work explores methods of integrating prior knowledge in the form of low-fidelity data into BO by evaluating these methods on synthetic multi-fidelity test functions. Our results highlight possibilities for improved convergence of the BO optimization. However, our work further highlights potential pitfalls of these multi-fidelity models, such as bias, convergence to local optima, and overfitting on low-fidelity data. Hence, leveraging low-fidelity data in multi-fidelity models can improve BO convergence, but there are instances where the algorithms are more susceptible to failure.

This article aims to tutorial a few important extremum seeking control approaches that can be used for the model-free optimization of industrial processes in various fields. The application of several methods is illustrated with a simple case study related to the production of algal biomass in photobioreactors. Other methods and applications are briefly reviewed.

Hydrogen isotope separation is a critical component of the fusion fuel cycle, particularly for achieving the desired purity levels of deuterium and tritium while minimising tritium inventory. This study investigates the cryogenic distillation of hydrogen isotopes, with a focus on the effects of isotopic equilibrium reactions at reduced temperatures and different system configurations. A one-column architecture was analysed to evaluate the impact of feed and side stream equilibrator temperatures and flowrates on separation performance and tritium inventory. Additionally, a two-column architecture was studied, incorporating multiple isotopic equilibrators in interconnecting streams, to further reduce unwanted heteronuclear isotopologues and improve system efficiency. Comparative analysis of the proposed configurations highlights significant operational advantages of optimising equilibrator temperatures, including reduced tritium contamination and inventory. Results indicate that reducing... [more]

Vacuum (pressure) swing adsorption (V(P)SA) has received considerable attention in the past decades. Existing studies typically estimate vacuum pump energy consumption using an approximate constant energy efficiency or an empirical energy efficiency correlation, leading to inaccurate representation of realistic vacuum pump performance. In this paper an enhanced computational approach is proposed for simulation and optimisation of V(P)SA through simultaneous integration of realistic vacuum pump data and adsorption bed fluidisation limits. The computational results show that the developed prediction models accurately represent the actual performance curves of the vacuum pump. Incorporation of the vacuum pump prediction models and fluidisation constraints in V(P)SA optimisation leads to significantly different optimal solutions compared to when these factors are not considered.

By capturing carbon dioxide from biomass flue gases, energy processes with negative carbon footprint are achieved. Among carbon capture methods, the fluidized temperature swing adsorption (TSA) column is a promising low-pressure alternative, but it has been developed on small scales. This work aims to model, simulate and optimize a fluidized TSA multi-stage equilibrium system to obtain a cost estimate and a conceptual design for future process scale up. A mathematical model described adsorption in multiple stages, each with a heat exchanger, coupled to the desorption operation. The model was based on elementary macroscopic molar and energy balances, coupled to pressure drops in a fluidized bed designed to operate close to the minimum fluidization velocity, and coupled to thermodynamics of adsorption equilibrium of a mixture of carbon dioxide and nitrogen in solid sorbents (the Toth equilibrium isotherm was used). The complete fluidized TSA process has been optimized to minimize costs,... [more]

The dividing wall column (DWC) has significant energy saving potential compared to conventional column sequences. However, to reach these savings in practice, it is essential that the control structures can track the optimal operation point despite inevitable changes in feed properties, performance characteristics and other uncertainties. Otherwise, the energy consumption may rise significantly or, more commonly, the DWC becomes unable to produce pure products even at its maximum reboiler duty. Extremum seeking control (ESC) is a model-free optimisation technique that may mitigate off-optimal operation in this environment. By active perturbation of selected manipulative variables, the algorithm infers gradient properties of the measured cost function and, by that, enables tracking of a moving optimum. Extremum seeking control can be used also in combination with other approaches, e.g. self-optimising control. Applied to the DWC, the presented perturb-and-observe algorithm, which can be... [more]

This work explores Aotearoa-New Zealand’s potential to fully electrify and source industrial process heat demands from renewable energy for 286 industrial sites while exploring the feasibility of green methanol production using excess electricity. Most energy models rely on spatially aggregated supply and demand, which limits the accurate representation of energy value chains. To address this limitation, the model incorporates industrial sites with varied temperature profiles, enabling the use of diverse heating technologies such as heat pumps, electrode boilers, bubbling fluidised bed reactors and biomass boilers. The proposed Mixed-Integer Linear Programming energy model uses the Accelerated Branch-and-Bound (ABB) algorithm, which is implemented within the P-graph framework to optimise the system. The model considers different energy transportation modes, including road transport for biomass and grid infrastructure for electricity. The multi-period design determines optimal heating t... [more]

This work presents a plantwide model of a Haber-Bosch ammonia synthesis loop (HB-loop) in a PtA plant, consisting of heat exchangers, compressors, steam turbines, flash separators and catalytic reactor beds. The total electrical power utility of the HB-loop is a combination of compressor power, refrigeration power, and steam turbine power. We optimise the HB-loop operating parameters, subject to constraints for maximum reactor temperatures, compressor choke and stall, minimum steam temperature, and maximum loop pressure. The loop features six degrees of freedom (DOFs) for the optimisation: three reactor temperatures, reactor N2/H2-ratio, separator temperature, and loop pressure. The optimisation minimises the total loop power utility for a given hydrogen make-up feed flow, with the PtA load varied by ranging the hydrogen make-up feed flow from 10 % to 120 % of the nominal. Across this load range, different constraints become active, with the compressor surge limit being particularly cr... [more]

The optimization of the power conversion system, responsible for thermal-to-electrical energy conversion, for a pulsed fusion power plant is presented. A spherical tokamak is modelled as three heat sources, all pulsed, with different stream temperatures and available amounts of heat. A thermal energy storage system is considered in the design to compensate for the lack of thermal power during a dwell. Thermal storage enables continued power generation during a dwell and can avoid thermal transients in sensitive components like turbomachines. Multiple lower grade heat sources are integrated into the process through parallel preheating trains. The evaluation of a dynamic model of the power conversion system is used to define an objective function with multiple criteria. A bi-objective optimization problem is defined to investigate the trade-off between the size of the thermal energy storage system and the variability in turbine power output during a dwell. The set of non-dominated design... [more]

Solar Thermal Power Plants (STPPs) use solar energy for large-scale electricity production but face significant operational challenges. These include variations in solar radiation, cloud cover, electricity demand fluctuations, and the need for frequent shutdowns if energy storage is inadequate. Deciding an optimal STPP operating conditions is challenging due to these factors. While revenue maximization has been used as an objective in existing literature, current models are often static and fail to capture the dynamic nature of STPPs. In contrast, this work proposes a dynamic model-based revenue optimization approach that accounts for plant dynamics and operational constraints, such as solar radiation variability and changing electricity demand. The objective function is designed to maximize revenue while considering power generation and fluctuating electricity prices. A simulation model of 1 MWe hybrid solar thermal power plant in Gurgaon, India, featuring two solar fields—Parabolic T... [more]

Enhanced weathering (EW) of alkaline minerals in chemical reactors with a controlled environment is recognized as a promising approach for gigaton-level carbon dioxide removal. However, reactor configuration and operating conditions must be optimized to balance the interfacial areas between gas, liquid and solid phases prior to industrial application. We developed a physics-based and data-driven hybrid modelling approach, coupled with multi-objective optimization, to study and compare three typical chemical reactors, i.e., trickle bed, packed bubbling columns, and stirred slurry reactors, and the optimal design to improve CO2 capture rate and reduce energy and water consumptions. Then an adaptive optimization is proposed to dynamically adjust the operating of the reactors in response to intermittent CO2 emission and renewable energy supply. Results indicated that forced stirring enhances CO2 capture rates by accelerating mass transport but increases energy consumption. Trickle bed reac... [more]

Seawater desalination is considered as the only available solution that can cope with the increasing demand for freshwater around the world. Improving the desalination techniques may help to cut off the cost and increase sustainability. In this paper, a mathematical model describing the MED process is developed within gPROMs software. The model includes all the necessary mass and energy balance equations together with thermodynamic and physical properties equations. The model predictions are validated against the actual plant data before using the model for optimizing the process to achieve minimum heat transfer area. For two different operating conditions (summer and winter) and a fixed production demand, the heat transfer area is minimised while optimising different parameters of the MED process. The results showed that a 10.4% reduction in the heat transfer area can be achieved under summer operating conditions and around 26% decrease in the heat transfer area can be met under winte... [more]

Asserting the phase stability of multi-component mixtures is an important task in computer-aided mixture/blend design (CAMbD), but it is often hindered by the lack of reliable and tractable models. In this paper, we propose a group-contribution machine-learning (GC-ML) method to predict phase coexistence for a large set of ternary mixtures consisting of two solvents and one (fixed) solute. Each solvent is represented by a vector of functional group numbers, encoded by integer values. The solvent vectors are combined with mixture composition and temperature to form the input features to a GC-ML surrogate classifier, which distinguishes between four types of stable phase configurations as possible outputs: liquid (L), solid-liquid (SL), liquid-liquid (LL) or solid-liquid-liquid (SLL). To explore the performance of the trained GC-ML multi-classifier, it is embedded as a surrogate phase-stability constraint in the optimisation of an ibuprofen crystallisation process. A two-step solution s... [more]

Process systems are negatively impacted by manufacturing uncertainties, and increasingly by unknown-unknown disruptive events. To this effect, systems need to be designed with the inherent flexibility and resilience to overcome the impacts of uncertainties and disruptions respectively as it is more challenging to retrofit existing systems with such capabilities. To this end, we propose a methodology based on flexibility analysis to systematically explore the feasibility of design alternatives under parameter uncertainty and discrete disruption scenarios simultaneously. Multi-parametric programming is utilized to generate explicit relationships between design decisions and the resulting system’s ability to maintain feasible operations under uncertainty and disruptive events. We capture this ability by introducing the Combined Flexibility-Resilience Index (CFRI), which describes the likelihood that the system is feasible under the relevant uncertainty and disruption sets. With explicit f... [more]

Model-based control and optimisation strategies can play a key role in improving energy efficiency and reducing emissions into produced water re-injection facilities. However, building a model that adequately describes the plant is challenging and can also be used in online decision-making procedures. This work proposes a simple model based on a real water re-injection facility operating on the Norwegian continental shelf. The results demonstrate the model's flexibility, which could be fitted to different plant operating points while being fast to solve when embedded in optimisation problems. The developed model is expected to aid the implementation of strategies like self-optimising control and real-time optimisation on produced water re-injection facilities.

Sustainable industrial development depends on optimizing resource and energy integration within Eco-industrial parks (EIPs), combined with stringent carbon emissions reduction policies. The main challenge is ensuring transparency, accountability, and data privacy while optimizing the conversion of raw materials and energy into valuable products and controlling emissions within EIPs. This research introduces an innovative framework to design optimized EIPs and deploy a blockchain-enabled trading platform for resources and emissions management, tackling these key issues. The proposed framework integrates EIPs with emission control policies, supported by two distinct smart contracts: one dedicated to blockchain-based resource trading and another handling financial transactions related to emission control policies, including other regulations such as income tax. The resource trading platform fosters transparency, enabling accurate tracking of material and energy flows. Furthermore, the fra... [more]

Real-world flowshop problems which are very common in the chemical industry are often difficult to solve in a reasonable time with allocation, sequencing, and lot-sizing decisions. Although great progress has been made in the last 20 years regarding MILP model formulations and solution algorithms, realistically-sized flowshop problems with resource and buffer constraints are still difficult to solve. On the other hand, discrete-event simulation (DES) allows for very detailed modelling of process plants, but lacking of optimization capabilities. Simulation Optimization (SO) combines the high-detail DES with mathematical optimization. We show that is possible to integrate MILP and DES using Benders decomposition. We explain the Benders-DES (BDES) approach with a small motivation example with makespan minimization objective and apply it to a real-world case study of a formulation plant with seven formulation and filling lines with sequencing, allocation, and lot-sizing decisions. We show... [more]

In the last decades, the growing awareness of climate change and the high political sensitivity of critical resources such as energy have emphasized a need for local, renewable and optimized energy mixes at various scales. Several studies have therefore aimed to optimize renewable energy technologies and plant locations to develop more renewable and efficient Energy Mixes. Following this trend, this paper applies and compares Goal Programming, Branch-and-Cut and NSGA-II to a multi-objective combinatorial optimization problem focused on the energy mix of Oise, France. Results show more optimality for Goal Programming and Branch-and-Cut, accompanied by a high sensitivity to constraints, while NSGA-II provides more technological diversity in the computed solutions.

This work proposes a two-level model, combining process and supply chain models, and an optimization framework for an integrated biorefinery system to convert agricultural residues into renewable urea via gasification routes. The process model of the gasification, ammonia and urea synthesis was developed in Aspen Plus® to identify key performance indicators such as energy consumption and relative yields for urea for different biomasses and operating conditions; then, these key process data were used in a mixed-integer linear programming (MILP) model, designed to identify the optimal combination of energy source, technological route of urea production and plant location that maximizes the net present value of the system. The model was applied to the whole Brazilian territory, divided into 5569 cities and 558 micro-regions. Each region’s agricultural production was evaluated to estimate biomass supply and urea demand. The Assis microregion, in close proximity with sugarcane and soybean c... [more]

Perishable consumer products like food, cosmetics, and household chemicals face challenges in supply chain management due to limited shelf life and uncertainties in demand and transportation. To address some of these issues, this work proposes a robust optimization framework for jointly optimizing facility allocation and inventory management. The framework determines optimal locations for distribution centers and their assigned customers, as well as inventory policies that minimize the total costs related to transportation, distribution, and storage under uncertain demand in a robust setting. Specifically, we develop a two-stage mixed-integer linear programming (MILP) model is that incorporates First-In-First-Out (FIFO) inventory policy to reduce spoilage. The bilinear FIFO constraints are linearized to improve computational efficiency. Social equity is integrated by defining a fairness index and incorporating it in facility allocation. Demand uncertainty is tackled using a robust opti... [more]

The urgent need to mitigate climate change and reduce reliance on fossil fuels highlights green hydrogen as a key component of the global energy transition. This study assesses the feasibility of producing hydrogen offshore using wind energy, focusing on economic and environmental aspects. Offshore wind energy offers several advantages: access to water for electrolysis, potentially lower hydrogen export costs compared to electricity, and storage systems that stabilize wind energy output. However, significant challenges remain, including the high costs of storage solutions, capital expenditures (CAPEX), and operational costs (OPEX). A Mixed-Integer Linear Programming (MILP) model optimizes the production units, storage, and distribution processes. A case study in southern France examines hydrogen production from a 150 MW floating wind farm. While hydrogen produced from offshore wind ranks among the most environmentally friendly, its costs remain high, and production volumes fall short o... [more]

In response to global efforts to reduce carbon emissions, the oil refining sector, a major source of industrial emissions, has set ambitious decarbonization targets. This study analyzes trade-offs between minimizing CO2 emissions and costs through the use of Pareto optimal solutions. A superstructure optimization framework evaluates various technological pathways and timelines, employing a bi-criterion optimization approach using the ?-constraint method. Results show that cost-effective, higher-emission solutions often involve natural gas-based technologies with carbon capture, while expensive, low-emission solutions favor electric power-based technologies. The analysis incorporates detailed assumptions about grid carbon intensity of varying degrees and accounts for varying national policies. Comparative case studies across locations highlight how grid carbon profiles influence optimal strategies, providing insights to inform local policies and incentivize technologies.

Superstructure based approaches have long been employed for optimal process synthesis problems. Due to the difficulties of using rigorous process models and simultaneous solutions, shortcut calculations have been the preferred means of modeling unit operations within larger process network problems. However, even with the use of shortcut equations to model the behaviour of unit operations, the resulting mixed-integer programs can be challenging to solve. Furthermore, generating the problem superstructure has often been done manually, presenting issues for scaling to larger problems. We demonstrate the use of an algorithmic approach to generate network superstructures for synthesis problems coupled with equation reformulations to yield an MIQCP (Mixed-Integer Quadratically Constrained Program) for networks of thermally coupled distillation columns. The combination of rapid problem generation with the ability to leverage recent advances in the performance of QCP (Quadratically Constraine... [more]

The unit commitment problem determines the optimal strategy to meet the electricity demand at minimum cost by committing power generation units at each point of time. Solving the unit commitment problem gives rise to a challenging optimization problem due to its combinatorial complexity and potentially long solution time requirements. Our proposed solution approach utilizes a decomposition method in conjunction with alternative models from the EGRET library. Results of this decomposition approach tested against four benchmarking systems show that significant computational speed ups are achieved.