Abstract:The UV flash problem is ubiquitous in process engineering for both dynamic simulation and dynamic optimization of vapor-liquid equilibrium processes. Such vapor-liquid equilibrium governed processes include flash separators, distillation columns, as well as enhanced oil recovery strategies for production of oil from an oil reservoir. A core modeling element in all these applications is the UV flash problem which describes equilibrium between the phases of a mixture. Coupling the flash equilibrium conditions with conservation of mass and energy results in a differential-algebraic model. We present an adjoint gradient-based single shooting algorithm for dynamic optimization of processes that involve UV flash constraints. This algorithm is relevant for nonlinear model predictive control of vapor-liquid equilibrium governed processes. We test the algorithm using two different approaches for numerical integration of the differential-algebraic equations. In one approach, the equilibrium equations obtained from the flash problem are solved separately from the conservation equations resulting in a set of nested loops. In the alternative approach, the equilibrium and conservation equations are solved simultaneously, thus in general reducing the total number of loops. Numerical integration of the differential-algebraic system is a key computational part of the dynamic optimization algorithm and the efficiency relies heavily on the underlying linear solver and thermodynamic library. Depending on the size of the system, evaluating thermodynamic properties may be significantly more time-consuming than solving the linear systems and thus becomes the bottleneck with regard to performance. We consider a numerical example involving dynamical UV flash constraints. The UV flash is considered particularly difficult as none of the intensive variables, such as pressure or temperature, are specified. This example illustrates that simultaneous numerical integration is faster than nested numerical integration and that the computation time of the underlying thermodynamic library is more than an order of magnitude higher than that of solving the linear systems arising in the numerical integration. The algorithm is implemented in Matlab for both types of numerical integration and in C using simultaneous numerical integration together with a thermodynamic library implemented in Fortran. The single shooting algorithm transforms the infinitedimensional dynamic optimization problem into a nonlinear program which may be solved using appropriate software libraries. The efficiency of the dynamic optimization algorithm is directly related to the efficiency of the optimization algorithm with respect to convergence, since every iteration will require up to several integrations of the differential-algebraic equations and computation of corresponding gradients. There exists a large number of commercial and noncommercial optimization libraries for nonlinear programs. We compare the efficiency of Matlabs fmincon, the open-source IPOPT library as well as the commercial libraries, KNITRO and NPSOL. We consider various combinations of using Intel and GNU compilers for the C and Fortran code as well as using Netlibs LAPACK library and Intels MKL for solving the linear systems. The performance tests show that the commercial software yields the best performance. More specifically a combination of using NPSOL together with Intel compilers and Intel MKL gives an improvement in performance of around sixty compared to a pure Matlab implementation.

GENERAL FLASH CALCULATION BY THE NEWTON- RAPHSON METHOD

Flash Calculation Using Successive Substitution Accelerated by the General Dominant Eigenvalue Method in Reduced-Variable Space: Comparison and New Insights

A design for flash on ray-tracing

A numerical method for two-dimensional nonlinear transient inverse heat conduction problems for orthotropic material

Artificial neural network assisted two-phase flash calculations in isothermal and thermal compositional simulations

A Line-Search-Based Algorithm for Multiphase Flash Calculations with CO 2 -Hydrocarbon System

An Arbitrary Lagrangian–Eulerian Method Based on the Adaptive Riemann Solvers for General Equations of State

Analytical Method of Nonlinear Coupled Constitutive Relations for Rarefied Non-Equilibrium Flows

Power Flow Calculation Based Newton-Raphson Calculation

A Novel Method to Accelerate the Three-phase Flash Calculation of the Hydrocarbon-CO 2 System

Quick Calculation and Applications of Adiabatic Flash

Accelerated Newton–Raphson Power Flow

Computational Aspects of the UV Flash Problem for Dynamic Optimization and Nonlinear Model Predictive Control

Rapid multi-component phase-split calculations using volume functions and reduction methods

CFD-Based Numerical Simulation of Water Film Flash Evaporation with a New Flash Evaporation Model

An Efficient and Robust Partial Differential Equation Solver by Flash-Based Computing in Memory

HIGH RESOLUTION KFVS FINITE VOLUME METHODS FOR MULTICOMPONENT FLOW CALCULATIONS

Investigation of 3D Internal Flow Using New Flux-Reconstruction High Order Method

Construction of a Minimum Energy Path for the VT Flash Model by an Exponential Time Differencing Scheme with the String Method

Construction of a Minimum Energy Path for the VT Flash Model by the String Method Coupled with the Exponential Time Differencing Scheme

Artificial Neural Network Accelerated Flash Calculation for Compositional Simulations