Distillation processes account for a substantial share of the industrial energy demand. Yet, these energy requirements can be reduced by a variety of energy integration methods, including various forms of direct heat integration, multi-effect distillation, thermal coupling and vapor recompression. Consequently, these intensification methods should be evaluated quantitatively in comparison to each other for individual separation tasks, instead of benchmarking single options with conventional sequences or relying on simplified heuristics. In order to overcome the computational burden of a broad assessment of a large number of process alternatives, a computationally-efficient framework for the energetic and economic evaluation of such energy integrated distillation processes is presented, which builds on thermodynamically-sound shortcut models that do not rely on constant relative volatility and constant molar overflow assumptions.

A new methodology for the assessment of the performance loss in catalytic reactors due to deactivation was developed and applied to fixed- and fluidized-bed CO methanation, with catalyst subject to coking. The methodology is based on the solution of heat and mass balances, by decoupling the reactor and deactivation dynamics. This is possible by using consecutive 1D, steady-state calculations for the characterization of the reactor performance. In this way, the progressively lower values of catalyst activity along the time on stream are computed with the integration of a dedicated dynamic model. This method has shown promising results in the characterization of the loss of performance of the reactor over time. The model correctly describes a progressive deactivation of the catalyst in fixed-bed reactors, while it shows that the decrease in activity is sudden for the whole reactor volume in fluidized bed reactors and occurs after a critical time-on-stream. Besides, it was observed that t... [more]

Hybrid modelling is widely employed in chemical engineering to generate highly accurate predictions. Such an approach merges first-principle modelling with machine learning techniques to identify and model the epistemic uncertainty from experimental data. Despite its advantages, this still requires cross-domain competencies that are difficult to find in the chemical industry and high human involvement. The possibility of automating the identification and training model would be significantly beneficial for the widespread adoption of hybrid modelling methodology within the chemical industry. This work presents a novel algorithm for the automatic identification of hybrid models (HMs) starting from the first-principle representation of the system, described by differential equation sets. The methodology formulates the problem as mixed-integer programming, identifying the equation running under uncertainty, identifying the machine learning model hyperparameters, and training the latter. Th... [more]

The introduction, peer review policy, and International Scientific Committee for Systems and Control Transactions volume 4 (ESCAPE 35 Proceedings)

This is the cover page and front matter for Systems and Control Transactions volume 4 (ESCAPE 35 Proceedings)

While process simulations often are either very rigid and accurate or very flexible and unprecise, informed decision making can only be maintained by establishing a detailed process model as early as possible within the project lifecycle while keeping relevant aspects of the process flexible enough. In this work, we present the development of a framework based on a dynamic interface between AspenPlus® process simulations and Python, enabling enhanced flexibility and automation for process modeling and optimization. This integration leverages the powerful simulation capabilities of AspenPlus® with the versatility of Py-thon for data analysis and optimization, delivering significant improvements in workflow efficiency and process control. By utilizing the dynamic simulation data exchange with Python, extensive parameter studies can be conducted.
In this provided dataset, the necessary input data, as well as the output files for each parameter run are provided. Furthermore, a .runtime an... [more]

This study proposes a model-based approach utilizing a hybrid population balance model (PBM) to optimize temperature profiles for minimizing agglomeration and enhancing crystal growth. The PBM incorporates key mechanisms—nucleation, growth, dissolution, agglomeration, and deagglomeration—and is applied to the crystallization of an industrial active pharmaceutical ingredient (API), Compound K. Parameters were estimated through prior design of experiments (DoE) and refined via additional thermocycle experiments. In-silico DoE simulations demonstrate that the hybrid PBM outperforms traditional methods in assessing process performance under agglomeration-prone conditions. Results confirm that thermocycles effectively reduce agglomeration and promote bulk crystal formation, though their efficiency plateaus beyond a certain cycle number. This model-based approach provides a more robust strategy for agglomeration control compared to conventional methods, offering valuable insights for industr... [more]

Mathematical modeling is essential for the effective control of many chemical engineering processes, including crystallization. However, most existing crystallization models used in industry and academia assume ideal mixing. As a result, the unclear effects of imperfect mixing on crystallization, reported in experimental studies, remain largely unexplained. In this work we aim to address this gap in understanding by examining antisolvent crystallization processes on a general theoretical level, using a novel dimensionless model. To address the impact of mixing on crystallization, we employ the Engulfment model coupled with a population balance, and we nondimensionalize the model equations. Using this model, we explore the dependence of the mean particle size on the homogenization rate, represented by the Damköhler number for crystallization. Moreover, we study the impact of mixing at various values of the model's kinetic parameters to simulate difference in properties of individual pro... [more]

As pharmaceutical manufacturing is transitioning from traditional batch to continuous manufacturing (CM), there is a lack of tools for CM design and development, which can integrate drug substance and drug product unit operations for overall evaluation. Recently, a Python-based PharmaPy framework was proposed to advance the design, simulation, and analysis of continuous pharmaceutical processes. However, the initial library of models only addressed upstream drug substance processing. In this work, new capabilities, including drug product unit operations such as feeder, blender, and tablet press, have been added to the PharmaPy framework, enabling end-to-end study and optimizing the effects of material properties and process conditions on solid oral dosage products. The platform supports computational efficiency and model accuracy by allowing the development of different mechanistic and semi-mechanistic models. Sensitivity analysis is performed on the integrated end-to-end simulator to... [more]

Mechanistic population balance modeling (PBM) has advanced the design of pharmaceutical crystallization processes, enabling the production of active pharmaceutical ingredient (API) crystals with desired critical quality attributes (CQAs), such as purity and crystal size distribution. However, PBM development can sometimes be resource-intensive, requiring extensive design of experiments (DoE) and high-quality process data, making it impractical under fast-paced industrial development timelines. This study proposes a machine learning (ML)-based workflow for developing ‘fit-for-purpose’ digital twins of crystallization processes, leveraging industrially available DoE data to link operating conditions with CQAs. Validated on industrial data for a commercial API with complex crystallization challenges, the workflow efficiently identifies optimal operating conditions, demonstrating the potential of data-driven digital twins to accelerate the development of pharmaceutical processes.

Tableting is a dry granulation process for compacting powder blends into tablets. In this process, a blend of active pharmaceutical ingredients (APIs) and excipients are fed into the hopper of a rotary tablet press via feeders. Inside the tablet press, rotating feed frame paddle wheels fill powder into dies, with tablet mass adjusted by the lower punch position during the die filling process. Pre-compression rolls press air out of the die, while main compression rolls apply the force necessary for compacting the powder into tablets. In this paper, process variables such as feeder screw speeds, feed frame impeller speed, lower punch position during die filling, and punch distance during main compression have been systematically varied. Corresponding responses, including pre-compression force, ejection force, and tablet porosity have been evaluated to optimize the tableting process. After implementing an open platform communications unified architecture (OPC UA) interface, process variab... [more]

Crystallization is a relevant process in the pharmaceutical industry for product purification and particle production. An efficient crystallization is characterized by crystals produced with the desired attributes. Therefore, modeling this process is a key point to achieve this goal. In this sense, the objective of this work is to propose a hybrid model to describe paracetamol dissolution in ethanol. The universal differential equations methodology is considered in the development of this model, using a neural network to predict the dissolution rate combined with the population balance equations to calculate the moments of the crystal size distribution (CSD) and the concentration. The model was developed using experimental batches. The dataset is composed of concentration measurements obtained using attenuated total reflectance-Fourier transform infrared (ATR-FTIR). The objective function of the optimization problem is to minimize the relative absolute difference between the experiment... [more]

In 2023, the International Council for Harmonisation (ICH) guideline for the development, implementation, and lifecycle management of pharmaceutical continuous manufacturing (PCM), was implemented in Europe. It promotes quality-by-design (QbD) and quality by control (QbC) strategies as well as the appropriate use of mathematical modelling. This development urges a harmonizing understanding across academia and industry for adoption of interpretable models instead of black-box models for advanced control strategies such as model predictive control (MPC), especially when applied in Good Manufacturing Practice (GMP) regulated areas. To this end, we first propose a comprehensive model development using Dynamic Mode Decomposition with Control (DMDc)to represent complex dynamics in a lower-dimensional space, disambiguating between underlying dynamics and actuation effects. Using data from a digital twin of PCM, our model demonstrates low computational complexity while effectively capturing no... [more]

The environmental impacts of pharmaceutical production underscore the need for comprehensive life cycle assessments (LCAs). Offshoring manufacturing, a common cost-saving strategy in the pharmaceutical industry, increases supply chain complexity and reliance on countries like India and China for active pharmaceutical ingredients (APIs). The COVID-19 pandemic exposed Europe’s vulnerability to global crises, prompting initiatives such as the French government’s re-industrialization plan to relocate the production of fifty critical drugs. Paracetamol production has been prioritized, with recent shortages highlighting the urgency to address supply chain risks while considering environmental impacts. This study uses process engineering to generate life cycle inventory (LCI) data for paracetamol production, offering an eco-design perspective. Aspen Plus was employed to model the API manufacturing process, integrating mass and energy balances to address the scarcity of LCI data. The results h... [more]

This study aims at developing a holistic approach to establish robust mathematical models of integrated and interactive multistep processes, while systematically identifying the corresponding design space and acceptable operating region (AOR). The overall objective is to reduce the experimentation costs, enhance accuracy of integrated metathetical models, and deliver built-in quality assurance based on a new Quality by Digital Design (QbDD) paradigm. This methodology starts with the construction of a set of model candidates for different unit operations, based on the prior knowledge and inherent assumptions. Several model candidates of the integrated multistep process are considered. A model discrimination based on model prediction performance reveals the best integrated model for the multistep process. In the next step, the estimability analysis and model-based design of experiment (MBDoE) are implemented to deliver information-rich data and systematically refine the integrated model.... [more]

This study presents a comprehensive comparison of the three alternative downstream manufacturing technologies for pharmaceuticals: i) Dry Granulation (DG) through roller compaction, ii) Direct Compaction (DC), and iii) Wet Granulation (WG) based on the economic, environmental and product quality performances. Firstly, the integrated dynamic mathematical models of the different downstream (drug product) processes were developed using gPROMS formulated products based on data from the literature or/and our recent experimental work. The process models were developed and simulated to reliably capture the impact of the different design options, process parameters, and material attributes. Uncertainty analysis was conducted using global sensitivity analysis to identify the set of critical process parameters (CPP) and critical material attributes (CMA) that mostly influence the quality and performance of the final pharmaceutical tablets in each case, captured by the critical quality attribute... [more]

Design space identification (DSId) and flexibility analysis are critical in process systems engineering, enabling efficient design of operating conditions. For bioprocess, these tasks are often hindered by the absence of reliable mechanistic models and limited experimental data. This paper presents an algorithm to address these challenges in bioprocesses. The methodology begins by constructing a Gaussian process (GP) model to predict key performance indicators (KPIs) from process inputs. Leveraging the probabilistic nature of GP predictions, we perform probabilistic design space identification (PDSId), characterizing each input point by its probability of feasibility which is the likelihood that constraints imposed on KPIs are satisfied. To visualize and analyse the feasibility space, contours at varying probability levels are identified using alpha shapes, which define deterministic boundaries corresponding to different confidence levels. This enables the quantification of volumetric... [more]

The pharmaceutical industry is a highly specialized field where strict quality control and accelerated time-to-market are essential for maintaining competitive advantage. Spherical crystallization has emerged as a promising approach in pharmaceutical manufacturing, offering significant potential to reduce equipment and operating costs, enhancing drug bioavailability, and facilitating compliance with product quality regulations. Emulsions, as an enabling technology for spherical crystallization, present unique advantages. However, the quality of spherical crystallization products derived from emulsions is significantly influenced by the intricate interactions between crystallization phenomena, formulation variables, and solution hydrodynamics. These complexities pose substantial challenges in determining optimal operational conditions to achieve the desired product characteristics. In this study, Bayesian optimization (BO) is employed to refine and optimize the operational conditions fo... [more]

The promising properties of lipid nanoparticles (LNPs) as drug carriers have been attracting significant attention in the field of drug delivery. However, further research is still required for a better understanding of their integration in the pharmaceutical industry. The Quality by Design (QbD) approach aims at ensuring the safety and efficiency in the development of new drugs, through an holistic, risk-based approach that gathers all sources of knowledge available about the system under analysis. One key resource of the QbD framework is the rich toolkit of Design of Experiments (DOE), to deepen the understanding of how the synthesis of LNPs by microfluidics can be effectively conducted and controlled. This study aimed to explore and understand the effectiveness of different DOE strategies, through an in silico study focused on the impact of factors related to the LNPs synthesis, namely the molar ratio of each lipid component in the lipidic mixture and the N/P ratio, while also consi... [more]

Understanding the complex dynamics of continuous processes in pharmaceutical manufacturing is essential to ensure product quality across the production line. This paper presents a data-driven modeling approach using Sparse Identification of Nonlinear Dynamics with Control (SINDYc) to capture the dynamics of a continuous direct compression (CDC) tableting line. By incorporating delayed control inputs into the candidate function library, the model effectively captures deviations from steady state in response to dynamic changes. The proposed model was developed by finding a balance between accuracy and sparsity, with focus on the ability to generalize to a wide range of operating conditions.

Reactive crystallization (RC) is a chemical process in which the reaction yields a crystalline product. It is used in various industries such as pharmaceutical manufacturing or water purification. In some cases, RC is the only feasible process pathway, such as the precipitation of certain ionic solids from solution. In other cases, a reaction can become a RC by changing the reaction environment to a solvent with low product-solubility. Despite the utility and prevalence of RC, it is not often emphasized in process design software. There are RC models that simulate the inner reactions and dynamics of a RC, but each has limiting assumptions, and are difficult to integrate with the rest of a process-line simulation. This modeling gap complicates RC process design and limits the exploration of the possible benefits to using RC as well as the ability to optimize a system that relies on it. To fill this gap, we built an open-source, customizable model that can be integrated with other unit o... [more]

This work introduces a Python-based interface that couples cradle-to-gate Life Cycle Assessment (LCA) with advanced process simulations in continuous biomanufacturing, resulting in dynamic process inventories and thus to dynamic LCA (dLCA). The open-source Brightway2.5 framework is used to dynamically track environmental inventories of the foreground process and LCA indicators (e.g. damage to ecosystems according to ReCiPE 2016) from the v3.10 cut-off ecoinvent database. The framework is applied to KTB1, a dynamic MATLAB–Simulink benchmark model of continuous Lovastatin production. 580 data points are computed across four different 24-hour scenarios. The difference between the hourly and the averaged foreground scenario is between 20-30%; a more pronounced deviation is observed when both background and foreground are averaged. The dLCA framework precisely identifies optimal periods for cleaner electricity usage, enabling future work on direct environmental feedback into process control... [more]

Emerging sectors within the biopharmaceutical industry are undergoing rapid scale-up due to the market boom of gene therapies and vaccine platform technologies. Manufacturers are pressured to orchestrate resources and plan investments under future demand uncertainty and, critically, an early-stage process uncertainty for platforms still under development. In this work, a multi-product multi-stage stochastic optimization problem integrating demand uncertainty is presented and augmented with a worst-case optimization approach with respect to process uncertainty. Results focus on a comparison between fixed equipment facilities and modular technologies, highlighting an inherent flexibility of the latter option due to shorter recourse actions for capacity scale-out. The impact of process uncertainty integration is quantified. With more conservative decisions taken in first-stages of the time horizon, expected costs result lower for modular single-use equipment. This suggests that capacity a... [more]

The increased use of pharmaceuticals globally over the past two decades has contributed to an increase in unused pharmaceuticals and a corresponding surge in pharmaceutical waste. Thus, there is an impetus for the development of processes for the recovery of the active pharmaceutical ingredients (APIs) from these unused drugs. This study introduces a decision framework for solvent selection to enable the recovery of APIs using a general separation train where cooling crystallization is the final step. The framework is designed to base solvent selection not just on the solubilities of the formulation contents but also considers the overall recovery that can be achieved in the process. In addition, the environmental sustainability of the framework is analyzed using the process mass intensity metric (PMI). The effectiveness of this framework is demonstrated by using paracetamol (PA) as a model API in a formulation consisting of five of the excipients commonly found in PA formulations. The... [more]

Choosing suitable solvents for crystallisation processes can be a challenging task when developing new pharmaceuticals, given the vast number of candidates available. To streamline this task, however, process modelling tools can be used to systematically probe the behaviour of different crystallisation setups entirely in-silico. In fact, it is possible to couple thermodynamic models with basic solid-liquid equilibria (SLE) principles to determine the impact of key process variables (e.g., temperature and solvent choice) on process performance, prior to conducting lab-scale experiments. In light of this, in this study we have used thermodynamic computational modelling tools (implemented within MATLAB®) to evaluate the cost and environmental impact of different batch crystallisation processes that may be used to manufacture flurbiprofen – a non-steroidal anti-inflammatory drug (NSAID) that can be used to treat various forms of arthritis. To complete this work, we have used the Apelblat e... [more]