On Two Flash Methods for Compositional Reservoir Simulations: Table Look-up and Reduced Variables
|
|
- Kerry Robinson
- 5 years ago
- Views:
Transcription
1 32 nd IEA EOR Annual Symposium & Workshop On Two Flash Methods for Compositional Reservoir Simulations: Table Look-up and Reduced Variables Wei Yan, Michael L. Michelsen, Erling H. Stenby, Abdelkrim Belkadi Center for Energy Resources Engineering (CERE) Technical University of Denmark October 18, 2011
2 Introduction Flash: for a mixture of composition z, will it split into two (or more) phases at specified T and P and what are the phase compositions and phase amounts? A summary of two recent studies: CSAT(table look-up): Belkadi et al., Comparison of two methods for speeding up flash calculations in compositional simulations, SPE compared with the shadow region method Reduced variables/reduction methods: Michelsen, M.L., Reduced variables - revisited, CERE Discussion Meeting 2011 compared with the conventional flash 2
3 Conventional flash Blind calculations without a priori information Two steps Stability analysis: whether the feed splits into two phases? The phase of composition z is stable at the specified (T,P) if and only if the tangent plane distance (TPD) ( ϕ ϕ ( )) tpd( w) = w ln w + ln ( w) ln z ln z 0 i i i i i i Phase split: calculate the equilibrium compositions using the initial estimates from the first step Old but robust, virtually no convergence problems More on safety than speed Michelsen, M. L. (1982a & b) Fluid Phase Equilibria 9: 1 19 & Michelsen and Mollerup (2007) Thermodynamic models: Fundamentals and Computational Aspects 3
4 Shadow region method Compositional simulations where information from previous calculations may be utilized (N F, x, tpd, ) Distinction between different regions by TPD Shadow region A. Unstable: one or two negative TPD B. Just stable: one trivial and one non trivial TPD=0 C. Single phase: one trivial and one non trivial TPD>0 D. Single phase: only trivial solutions Rasmussen et al. (2006) SPE Res Eval & Eng 9: Michelsen and Mollerup (2007) Thermodynamic models: Fundamentals and Computational Aspects. 4
5 Compositional Space Adaptive Tabulation (CSAT) method CSP/CST/CSAT 5 inspired by the 1D analytical solution of gas injection a few key tie-lines in the solution path. CSP based table look-up approach to replace stability test/phase split Procedure Tie-line tables constructed either in advance or adaptively For a new feed z Vapor fraction Criterion Voskov and Tchelepi (2007) SPE k z j x j β = k k y j x j k k 2 zi ( β yi + (1 β ) xi ) < ε i Voskov and Tchelepi (2008) Transport in Porous Media 75: k tieline index j component index arbitrarily chosen if one of the stored tie-lines satisfies these criteria for the composition of interest, the table is used to look up the flash results. Otherwise, a standard EOS based phase behavior procedure is employed.
6 Tie-line Table Look-up (TTL) our implementation of CSAT Only for phase split step to approximate flash results in twophase region A unique distance for each tie-line k k Shortest distance ( d = e ) from the feed z to tie-line k ( ( β i i ( β ) i )) 2 e k = ( d k ) = min z y k + 1 x k d k tie-line distance i The corresponding β is readily obtained as β = T ( z x k ) ( y k x k ) T ( y k x k ) ( y k x k ) If e k <ε, accept tie-line k as flash solution, and β calculated from Eq.(5) If e k >ε for all the M tie-lines in the table, flash the composition and update the tie-line table if it is two-phase. 2 Eq.(5) 6
7 Gas injection systems tested Tested with 1-D slimtube simulation with 500 cells System 1 System 2 System 3 System 4 Oil 13-component oil Zick Oil 1* Zick Oil 2* Zick Oil 3* Gas 0.8 CO C 1 Zick Gas 1* Zick Gas 2* Zick Gas 3* T (K) P (atm) EoS used SRK PR PR PR * 12-component fluid description from Jessen (2000) Ph.D. thesis or Orr (2007) Gas Injection Processes. 7
8 Analysis of CSAT using System 1 The influence of number of tie-lines M and the tolerance on simulation time and %skips of flash calculations Decreasing ε Increasing M ε =10-4 ε =10-5 ε =10-6 ε =10-7 M = 100 Time (sec) % skips 41% 0.1% 0.0% 0.0% M = 500 Time (sec) % skips 99.9% 10% 0.3% 0.2% M = 1000 Time (sec) % skips 99.9% 18% 0.9% 0.4% M = 5000 Time (sec) % skips 99.9% 64.7% 25.8% 7% Larger M increases simulation time and %skips Smaller ε increases simulation time but decreases %skips Sorting tie-lines gives limited help 8
9 Analysis of CSAT using System 1 1 Accurate solution 0.8 CSAT M=1000 eps=1e-4 CSAT M=1000 eps=1e-5 Gas saturation CSAT M=1000 eps=1e-6 CSAT M=1000 eps=1e Cell number ε =10-4 ε =10-5 ε =10-6 ε =10-7 M = 1000 Time (sec) % skips 99.9% 18% 0.9% 0.4% ε=10-4 not accurate; ε =10-6 or 10-7 too few skips. Higher M requires even smaller ε 9
10 TTL with pre-calculated tie-lines (TTL-PRE) The tie-line table can be calculated in advance to reduce simulation time Use M = and ε = 10-8 to find the most frequently used tie-lines during the simulation. 3 tie-lines are identified, accounting for 88% of hits 1 Gas saturation Accurate solution Recovery 0.8 CSAT-PRE eps=1e-4 80 CSAT-PRE eps=1e-5 70 Gas saturation CSAT-PRE eps=1e-6 CSAT-PRE eps=1e-7 Recovery (%) Accurate solution CSAT-PRE eps=1e-4 CSAT-PRE eps=1e-5 CSAT-PRE eps=1e-6 CSAT-PRE eps=1e Cell number PVI 10
11 System 1: simulation times 11 PVI=0.5 Time Direct Time (sec) approximation in (sec) two-phase * PVI=1.2 Direct approximation in two-phase * Conventional/ Full stability TTL M=100, ε = % % M=500, ε = % % M=1000, ε = % % M=5000, ε = % % TTL-PRE (three tielines) ε = % % ε = % % ε = % % ε = % % Shadow region * Reported numbers are percentages of total flashes in two-phase region
12 Tie-line Distance Based Approximation (TDBA) an alternative and simpler Just compare one tie-line in the same cell from a previous rigorous flash using tie-line distance. Procedure Calculate e k as before (only one) If e>ε, do new flash, and update the tie-line if it is two-phase If e<10-4 ε, use the previous results as a solution without any adjustment If ε >e>10-4 ε, use the previous results with adjustment Option 1 (TDBA1): use old K values to solve Rachford-Rice Eq. Option 2 (TDBA2): use vapor split factors θ i to adjust directly θ β y i i = β yi + ( 1 β ) xi old v z θ = l = z ( θ ) i i, new i i i, new 1 i 12
13 System 1: simulation times Conventional/ Full stability TTL-PRE (three tielines) Time (sec) PVI=0.5 Approx. with adjustment in two-phase * Direct approximation in two-phase * Time (sec) PVI=1.2 Approx. with adjustment in two-phase * Direct approximation in two-phase * ε = % % ε = % % TDBA1 ε = % 11% % 12% ε = % 11% % 11% ε = % 8% % 10% ε = % 7% % 10% TDBA2 ε = % 11% % 13% ε = % 10% % 11% ε = % 9% % 11% ε = % 6% % 9% Shadow region * Reported numbers are percentages of total flashes in two-phase region
14 TDBA1 results for System Accurate solution TDBA1 eps=1e-4 80 Gas saturation Gas saturation Accurate solution TDBA1 ε=10-4 TDBA1 ε=10-5 TDBA1 ε=10-6 TDBA1 ε=10-7 TDBA1 eps=1e-5 TDBA1 eps=1e-6 TDBA1 eps=1e-7 Recovery (%) Recovery (%) Accurate solution TDBA1 ε=10-4 TDBA1 ε=10-5 TDBA1 ε=10-6 TDBA1 ε=10-7 Accurate solution TDBA1 eps=1e-4 TDBA1 eps=1e-5 TDBA1 eps=1e-6 TDBA1 eps=1e Cell number Cell# PVI PVI 14
15 TDBA s potential : speeding up complicated EoS s 6-component gas injection simulated by PC-SAFT and SRK Simulation time (sec) SRK CPA PC-SAFT CPA new PC-SAFT new SRK+TDBA CPA+TDBA PC-SAFT+TDBA Simulation time ratio CPA PC-SAFT CPA new PC-SAFT new CPA+TDBA w.r.t. SRK+TDBA PC-SAFT+TDBA w.r.t. SRK+TDBA CPA+TDBA w.r.t. SRK PC-SAFT+TDBA w.r.t. SRK Number of components Number of components Speed-up 1: SPE (solid dashed ) Speed-up 2: TDBA1 (dashed dash-dot) 15
16 Summary on approximation methods CSAT/TTL saves the time for rigorous flash but managing the tieline table can be a significant overhead. The simulation time increases dramatically with the number of tie-lines used. Big tolerances lead to inaccurate results. TTL-PRE is better but gives limited speeding-up compared with the shadow region algorithm. TDBA is simpler and cuts the simulation time by 1/3 to 1/2. The approximation methods may have potential to speed up simulation with complicated EoS s. 16
17 Reduced variables methods basics Solution procedure for equilibrium calculations with a cubic EoS where the matrix of BIP s is of low rank B C = i b n i i P nrt A = V B V ( V + B ) C C = aijni n j aij = aiia jj (1 kij ) i j If all BIPs are zero, Ai = 2 aii aij n j and A ˆ ϕ = C + C A + C b ln i n a i b i C A A = = 2 n Consequently, the vector of ln ˆ ϕ i can be written as a linear combination of 3 pre-calculated vectors, with i th elements a1, ii and b i. Same applies to the lnk i. C a n i ij j i j j * ln ˆ ϕ i = C n + C a a ii + C b b i 17
18 A brief history First - to our knowledge - used 30 years ago by Michelsen and Heideman (1981) for critical point calculation Suggested for flash calculations by Michelsen (1986) Single nonzero BIP row/column, Jensen and Fredenslund (1987) Generalized for nonzero BIPs by Hendriks (1992) Extensively used in the generalized version for the last 20 years Its advantages first questioned in public by Haugen and Beckner in 2011 (SPE ) 18
19 Arguments against reduced variables Essentially restricted to the cubic EOS Difficult to be formulated as unconstrained minimization problems consequently, less safe. More cumbersome composition derivatives The simple algebraic operations required to evaluate A i are today very inexpensive (SIMD) Our experience: Effort of the conventional approach grows approximately linearly with C, not proportionally with C 2 or C 3. A fair comparison between minimization based reduced variables method and conventional flash requires substantial coding (perhaps modest potential for improvement) But recent development by Nichita and Graciaa (2010) enabled an adaption to Michelsen s existing code without extensive modifications! 19
20 Reduced variables by spectral decomposition Consider the matrix P with elements P ij =1-k ij. We calculate the spectral decomposition P C = k= 1 λ u u T k k k λ k where is the k th eigenvalue of P and u the corresponding eigenvector. The eigenvalues are numbered in decreasing magnitude. Assume now that the eigenvalues are negligible for k > M where M << C. For example: All BIP s equal to zero: M=1 Upper triangle of BIP s zero except for a few rows, i.e., k = 0 for i > m and j > i ij In this case, M = 2m + 1 and the match is exact 20
21 Cont d T We then get P λ u u and and A = 2 a n = 2 λ e e n = d e with Net results M = k= 1 k k k C M C M i ij j k ik jk j k ik k = 1 k = 1 j= 1 k = 1 ln ˆ ϕ i Vector of : Linear combination of 2m + 3 vectors Results identical to full approach Computational effort reduced from C 2 to 2CM plus overhead! Successive substitution Conventional implementation, where the reduction method is only implemented to calculate A i Acceleration as usual No effect on convergence behavior ij k ii ik jj jk k ik jk k= 1 k= 1 Used for stability analysis, as well as for phase split M a = λ a u a u = λ e e d M C = 2λ e n k k jk j j= 1 21
22 Second order Gibbs energy minimization ln K M + 2 i = cleil i, M 1 1 l= 1 e = + ei, M 2 + = b i Independent variables c 1, c 2,, and c M+2 Gradient C G G vi = or g c = Wg c v c v i j i= 1 i j where is vapor moles i and W = v / c (from Rachford-Rice equation) ij i j Hessian H c WHW T W ij looks complex to calculate, but simple algebraic expressions for the elements can be derived. 22
23 Procedure for the 2 nd order minimization Calculate the K-factors from c Solve the Rachford-Rice equation to get v i Calculate conventional gradient and Hessian Calculate transformation matrix W Calculate c-based gradient and Hessian Calculate corrected c using trust-region approach Similar procedure for Stability Analysis How does it compare to the classic approach? 23
24 Alternative simplification: sparse k C j= 1 a a (1 k ) n = a ( S S ) ii jj ij j ii i where S C = j= 1 a n jj j and S i C a k n i m jj ij j j= 1 = m jj ij j > j= 1 a k n i m Uses approximately 2mC multiplications. 24
25 Test examples Example 1 Modified SPE3 with 9 components. Modified such that all k ij = 0 for i > 3, j > 3. Non-zero interaction coefficients for N 2, CO 2 and CH 4. Example 2 Removal of N 2 and CO 2. Only CH 4 has nonzero BIPs. Phase diagram largely unmodified, as the content of the removed species was small. Only 5 variables in the reduced case! In both tests, the mixture was expanded to 27 or 25 components by subdividing the last species. One million flash calculations in an equidistant 1000 by 1000 grid in T and P. All calculations are blind. About 60% two-phase. 25
26 Test example 1 26
27 Test example 2 27
28 Summary on reduced variables Modest effect of increasing C: Linear, but less than proportional No advantage of reduction methods over alternatives Fastest result: utilize sparsity of BIP-matrix Computing times are in general very implementation dependent Other implementations of reduction methods might be more efficient than the one used here. To be convincing, they should be able to beat the current computing times. 28
29 Acknowledgment Danish Council for Technology and Production Sciences 29
Enhanced Oil Recovery with CO2 Injection
Enhanced Oil Recovery with CO2 Injection Wei Yan and Erling H. Stenby Department of Chemical Engineering Technical University of Denmark Contents Overview Mechanism of miscibility Experimental study of
More informationPVTpetro: A COMPUTATIONAL TOOL FOR ISOTHERM TWO- PHASE PT-FLASH CALCULATION IN OIL-GAS SYSTEMS
PVTpetro: A COMPUTATIONAL TOOL FOR ISOTHERM TWO- PHASE PT-FLASH CALCULATION IN OIL-GAS SYSTEMS A. M. BARBOSA NETO 1, A. C. BANNWART 1 1 University of Campinas, Mechanical Engineering Faculty, Energy Department
More informationTIE-SIMPLEX METHOD FOR THERMAL-COMPOSITIONAL SIMULATION
TIE-SIMPLEX METHOD FOR THERMAL-COMPOSITIONAL SIMULATION A DISSERTATION SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING AND THE COMMITTEE ON GRADUATE STUDIES OF STANFORD UNIVERSITY IN PARTIAL
More information18 a 21 de novembro de 2014, Caldas Novas - Goiás THERMODYNAMIC MODELING OF VAPOR-LIQUID EQUILIBRIUM FOR PETROLEUM FLUIDS
18 a 21 de novembro de 2014, Caldas Novas - Goiás THERMODYNAMIC MODELING OF VAPOR-LIQUID EQUILIBRIUM FOR PETROLEUM FLUIDS Antonio Marinho Barbosa Neto, aneto@dep.fem.unicamp.br 1 Jônatas Ribeiro, jonand@petrobras.com.br
More informationarxiv: v1 [physics.comp-ph] 27 Jan 2019
arxiv:1901.09380v1 [physics.comp-ph] 27 Jan 2019 Acceleration of the NVT-flash calculation for multicomponent mixtures using deep neural network models Yiteng Li 1, Tao Zhang 1 and Shuyu Sun 1,2 1 Physical
More informationMULTIPHASE EQUILIBRIUM CALCULATIONS WITH GAS SOLUBILITY IN WATER FOR ENHANCED OIL RECOVERY
MULTIPHASE EQUILIBRIUM CALCULATIONS WITH GAS SOLUBILITY IN WATER FOR ENHANCED OIL RECOVERY A THESIS SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT
More informationThe Rachford-Rice Contest
The Rachford-Rice Contest If you re a student, here s your chance to win $000! At least the first three who outperform the Rachford-Rice algorithm of phantom-engineer Aaron Zick. Professionals who beat
More informationTIE-LINE BASED PARAMETERIZATION FOR THERMAL COMPOSITIONAL RESERVOIR SIMULATION
TIE-LINE BASED PARAMETERIZATION FOR THERMAL COMPOSITIONAL RESERVOIR SIMULATION A REPORT SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT OF THE
More informationFluid Phase Equilibria
Fluid Phase Equilibria 267 (28) 6 17 Contents lists available at ScienceDirect Fluid Phase Equilibria journal homepage: www.elsevier.com/locate/fluid A robust negative flash based on a parameterization
More informationAn improved component retrieval method for cubic equations of state with non-zero binary interaction coefficients for natural oil and gas
J Petrol Explor Prod Technol (26) 6:243 25 DOI.7/s322-5-8-y ORIGINAL PAPER - PRODUCTION ENGINEERING An improved component retrieval method for cubic equations of state with non-zero binary interaction
More informationHamdi Tchelepi (Principal advisor)
THERMODYNAMIC EQUILIBRIUM COMPUTATION OF SYSTEMS WITH AN ARBITRARY NUMBER OF PHASES FLOW MODELING WITH LATTICE BOLTZMANN METHODS: APPLICATION FOR RESERVOIR SIMULATION A REPORT SUBMITTED TO THE DEPARTMENT
More informationCALCULATION OF THE COMPRESSIBILITY FACTOR AND FUGACITY IN OIL-GAS SYSTEMS USING CUBIC EQUATIONS OF STATE
CALCULATION OF THE COMPRESSIBILITY FACTOR AND FUGACITY IN OIL-GAS SYSTEMS USING CUBIC EQUATIONS OF STATE V. P. de MATOS MARTINS 1, A. M. BARBOSA NETO 1, A. C. BANNWART 1 1 University of Campinas, Mechanical
More informationCOMPOSITIONAL SPACE PARAMETERIZATION METHODS FOR THERMAL-COMPOSITIONAL SIMULATION
COMPOSITIONAL SPACE PARAMETERIZATION METHODS FOR THERMAL-COMPOSITIONAL SIMULATION A DISSERTATION SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING AND THE COMMITTEE ON GRADUATE STUDIES OF STANFORD
More informationUNIVERSITY OF CALGARY. Modelling Water-Hydrocarbon Mutual Solubility in Multiphase Equilibrium Calculations. Hongbo Yu A THESIS
UNIVERSITY OF CALGARY Modelling Water-Hydrocarbon Mutual Solubility in Multiphase Equilibrium Calculations by Hongbo Yu A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILMENT OF THE
More informationHydrocarbon Reservoirs and Production: Thermodynamics and Rheology
Hydrocarbon Reservoirs and Production: Thermodynamics and Rheology A comprehensive course by Prof. Abbas Firoozabadi RERI and Yale University and Prof. Gerald Fuller Stanford University Palo Alto, California
More informationP1: IML/FFX P2: IML/FFX QC: IML/FFX T1: IML AT029-FM AT029-Manual AT029-Manual-v8.cls December 11, :59. Contents
Contents Foreword Preface xvii xix Chapter 1 Introduction 1 Nomenclature 1 1.1 Nature of Petroleum Fluids 1 1.1.1 Hydrocarbons 3 1.1.2 Reservoir Fluids and Crude Oil 5 1.1.3 Petroleum Fractions and Products
More informationPETE 310 Lectures # 36 to 37
PETE 310 Lectures # 36 to 37 Cubic Equations of State Last Lectures Instructional Objectives Know the data needed in the EOS to evaluate fluid properties Know how to use the EOS for single and for multicomponent
More information"Energy Applications: Impact of Data and Models"
"Energy Applications: Impact of Data and Models" Energy Applications refers in this particular case to the wide application of equations of state upstream in the Production of Oil and Gas. The petroleum
More informationPreliminary Evaluation of the SPUNG Equation of State for Modelling the Thermodynamic Properties of CO 2 Water Mixtures
Available online at www.sciencedirect.com Energy Procedia 26 (2012 ) 90 97 2 nd Trondheim Gas Technology Conference Preliminary Evaluation of the SPUNG Equation of State for Modelling the Thermodynamic
More informationUNIVERSITY OF CALGARY. Development of a Four-Phase Compositional Simulator Using Equations of State. Yizheng Wei A THESIS
UNIVERSITY OF CALGARY Development of a Four-Phase Compositional Simulator Using Equations of State by Yizheng Wei A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILMENT OF THE REQUIREMENTS
More informationMachine Learning Applied to 3-D Reservoir Simulation
Machine Learning Applied to 3-D Reservoir Simulation Marco A. Cardoso 1 Introduction The optimization of subsurface flow processes is important for many applications including oil field operations and
More informationThe Pennsylvania State University. The Graduate School. Department of Energy and Mineral Engineering MATHEMATICS OF MULTIPHASE MULTIPHYSICS
The Pennsylvania State University The Graduate School Department of Energy and Mineral Engineering MATHEMATICS OF MULTIPHASE MULTIPHYSICS TRANSPORT IN POROUS MEDIA A Dissertation in Energy and Mineral
More information8.1 Concentration inequality for Gaussian random matrix (cont d)
MGMT 69: Topics in High-dimensional Data Analysis Falll 26 Lecture 8: Spectral clustering and Laplacian matrices Lecturer: Jiaming Xu Scribe: Hyun-Ju Oh and Taotao He, October 4, 26 Outline Concentration
More informationA Multi-Continuum Multi-Component Model for Simultaneous Enhanced Gas Recovery and CO 2 Storage in Stimulated Fractured Shale Gas Reservoirs Jiamin
A Multi-Continuum Multi-Component Model for Simultaneous Enhanced Gas Recovery and CO 2 Storage in Stimulated Fractured Shale Gas Reservoirs Jiamin Jiang M.S. Candidate Joined Fall 2013 1 Main Points Advanced
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 13
STAT 309: MATHEMATICAL COMPUTATIONS I FALL 208 LECTURE 3 need for pivoting we saw that under proper circumstances, we can write A LU where 0 0 0 u u 2 u n l 2 0 0 0 u 22 u 2n L l 3 l 32, U 0 0 0 l n l
More informationModeling and Computation of Phase Equilibrium. Using Interval Methods. Department of Chemical and Biomolecular Engineering, University of Notre Dame
Modeling and Computation of Phase Equilibrium Using Interval Methods Mark A. Stadtherr Department of Chemical and Biomolecular Engineering, University of Notre Dame Notre Dame, IN, USA 2nd Scandinavian
More informationEquations of State. Equations of State (EoS)
Equations of State (EoS) Equations of State From molecular considerations, identify which intermolecular interactions are significant (including estimating relative strengths of dipole moments, polarizability,
More informationChapter 7 PHASE EQUILIBRIUM IN A ONE-COMPONENT SYSTEM
Chapter 7 PHASE EQUILIBRIUM IN A ONE-COMPONENT SYSTEM 7.1 INTRODUCTION The intensive thermodynamic properties of a system are temperature, pressure, and the chemical potentials of the various species occurring
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More informationDISTILLATION. Keywords: Phase Equilibrium, Isothermal Flash, Adiabatic Flash, Batch Distillation
25 DISTILLATION Keywords: Phase Equilibrium, Isothermal Flash, Adiabatic Flash, Batch Distillation Distillation refers to the physical separation of a mixture into two or more fractions that have different
More informationPrediction of surface tension of binary mixtures with the parachor method
Prediction of surface tension of binary mixtures with the parachor method Tomáš Němec 1,a Institute of Thermomechanics ASCR, v.v.i., Dolejškova, 18 Praha 8, Czech Republic Abstract. The parachor method
More informationTo a large extent abandoned once gradient-based methods. One applicable class of methods is Generating set search
Derivative Free Optimization and Average Curvature Information Trond Steihaug Lennart Frimannslund Eighth US-Mexico Workshop on Optimization and its Applications, January th, 007 Compass search in R. Current
More informationELEMENT-BASED FORMULATIONS FOR COUPLED FLOW, TRANSPORT, AND CHEMICAL REACTIONS
ELEMENT-BASED FORMULATIONS FOR COUPLED FLOW, TRANSPORT, AND CHEMICAL REACTIONS A REPORT SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT OF THE
More informationBusiness. Final Exam Review. Competencies. Schedule Today. Most missed on Exam 3. Review Exam #3
Business Final Exam Review Online course evaluation (19/32 = 59%) Counts as a homework assignment (by Thurs) Professional program application Past due! Case study due today by 5 pm Leadership evaluation
More informationOn the estimation of water pure compound parameters in association theories
On the estimation of water pure compound parameters in association theories Andreas Grenner, Georgios M. Kontogeorgis, Michael L. Michelsen, Georgios K. Folas To cite this version: Andreas Grenner, Georgios
More informationChemical Equilibrium: A Convex Optimization Problem
Chemical Equilibrium: A Convex Optimization Problem Linyi Gao June 4, 2014 1 Introduction The equilibrium composition of a mixture of reacting molecules is essential to many physical and chemical systems,
More informationLinear Algebra Review. Vectors
Linear Algebra Review 9/4/7 Linear Algebra Review By Tim K. Marks UCSD Borrows heavily from: Jana Kosecka http://cs.gmu.edu/~kosecka/cs682.html Virginia de Sa (UCSD) Cogsci 8F Linear Algebra review Vectors
More informationA modification of Wong-Sandler mixing rule for the prediction of vapor-liquid equilibria in binary asymmetric systems
Korean J. Chem. Eng., 28(7), 16131618 (2011) DOI: 10.1007/s1181401005347 INVITED REVIEW PAPER A modification of WongSandler mixing rule for the prediction of vaporliquid equilibria in binary asymmetric
More informationThe Lanczos and conjugate gradient algorithms
The Lanczos and conjugate gradient algorithms Gérard MEURANT October, 2008 1 The Lanczos algorithm 2 The Lanczos algorithm in finite precision 3 The nonsymmetric Lanczos algorithm 4 The Golub Kahan bidiagonalization
More informationMath 471 (Numerical methods) Chapter 3 (second half). System of equations
Math 47 (Numerical methods) Chapter 3 (second half). System of equations Overlap 3.5 3.8 of Bradie 3.5 LU factorization w/o pivoting. Motivation: ( ) A I Gaussian Elimination (U L ) where U is upper triangular
More informationReliable Prediction of Phase Stability Using an Interval Newton Method
Reliable Prediction of Phase Stability Using an Interval Newton Method James Z. Hua a, Joan F. Brennecke b and Mark A. Stadtherr a a Department of Chemical Engineering, University of Illinois, 600 S. Mathews
More informationCS 246 Review of Linear Algebra 01/17/19
1 Linear algebra In this section we will discuss vectors and matrices. We denote the (i, j)th entry of a matrix A as A ij, and the ith entry of a vector as v i. 1.1 Vectors and vector operations A vector
More informationREV. CHIM. (Bucureºti) 58 Nr
1069 High-Pressure Vapour-Liquid Equilibria of Carbon Dioxide + 1-Pentanol System Experimental Measurements and Modelling CATINCA SECUIANU*, VIOREL FEROIU, DAN GEANÃ Politehnica University Bucharest, Department
More informationCH2351 Chemical Engineering Thermodynamics II Unit I, II Phase Equilibria. Dr. M. Subramanian
CH2351 Chemical Engineering Thermodynamics II Unit I, II Phase Equilibria Dr. M. Subramanian Associate Professor Department of Chemical Engineering Sri Sivasubramaniya Nadar College of Engineering Kalavakkam
More informationOctober 25, 2013 INNER PRODUCT SPACES
October 25, 2013 INNER PRODUCT SPACES RODICA D. COSTIN Contents 1. Inner product 2 1.1. Inner product 2 1.2. Inner product spaces 4 2. Orthogonal bases 5 2.1. Existence of an orthogonal basis 7 2.2. Orthogonal
More informationAzeotropic distillation Example 1
Azeotropic distillation Example 1 Separate water from iso-butanol. The phase behavior for this mixture is interesting. There is a minimum boiling azeotrope formed as well as a liquid-liquid phase separation
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
More informationAll Rights Reserved. Armando B. Corripio, PhD, P.E., Multicomponent Distillation Column Specifications... 2
Multicomponent Distillation All Rights Reserved. Armando B. Corripio, PhD, P.E., 2013 Contents Multicomponent Distillation... 1 1 Column Specifications... 2 1.1 Key Components and Sequencing Columns...
More informationContents and Concepts
Contents and Concepts 1. First Law of Thermodynamics Spontaneous Processes and Entropy A spontaneous process is one that occurs by itself. As we will see, the entropy of the system increases in a spontaneous
More informationContents and Concepts
Contents and Concepts 1. First Law of Thermodynamics Spontaneous Processes and Entropy A spontaneous process is one that occurs by itself. As we will see, the entropy of the system increases in a spontaneous
More informationCFD Simulation of Flashing and Boiling Flows Using FLUENT
CFD Simulation of Flashing and Boiling Flows Using FLUENT Hua Bai and Paul Gillis The Dow Chemical Company FLUENT UGM 2004 Liquid/Gas Phase Change found in many industrial chemical processes involves complex
More informationAccuracy of vapour ^ liquid critical points computed from cubic equations of state
High Temperatures ^ High Pressures 2000 volume 32 pages 449 ^ 459 15 ECTP Proceedings pages 433 ^ 443 DOI:10.1068/htwu303 Accuracy of vapour ^ liquid critical points computed from cubic equations of state
More informationNumerical Aspects of the SAFT Equation of State
University of Rhode Island DigitalCommons@URI Senior Honors Projects Honors Program at the University of Rhode Island 006 Numerical Aspects of the SAFT Equation of State Leah M. Octavio University of Rhode
More informationSimulating Condensation in a Supercritical Gas Jet
Simulating Condensation in a Supercritical Gas Jet L. Qiu, Y. Wang, H. Wang, Q. Jiao and R. D. Reitz * Engine Research Center University of Wisconsin Madison, WI 53706 USA Abstract The KIVA code was modified
More informationThermodynamic and Stochiometric Principles in Materials Balance
Thermodynamic and Stochiometric Principles in Materials Balance Typical metallurgical engineering problems based on materials and energy balance NiO is reduced in an open atmosphere furnace by excess carbon
More informationVapor-liquid Separation Process MULTICOMPONENT DISTILLATION
Vapor-liquid Separation Process MULTICOMPONENT DISTILLATION Outline: Introduction to multicomponent distillation Phase Equilibria in Multicomponent Distillation (Pg. 737) Bubble-point and dew-point calculation
More informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More informationModeling of Pressure Dependence of Interfacial Tension Behaviors of Supercritical CO 2. + Crude Oil Systems Using a Basic Parachor Expression
19 Modeling of Pressure Dependence of Interfacial Tension Behaviors of Supercritical + Crude Oil Systems Using a Basic Parachor Expression Saini Dayanand* California State University, Bakersfield, CA,
More informationS (13) Reference: FLUID 9791
Title: A Deduction of the MulticriticalityConditionsof, Mixturesfrom the Gibbs Tangent Plane Criterion Author: Nélio Henderson Wagner.F. Sacco Raimundo A. Rodrigues Jr PII: S0378-3812(13)00550-5 DOI: http://dx.doi.org/doi:10.1016/j.fluid.2013.09.043
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More informationA FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS
A FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS Ruben Juanes 1 and Knut-Andreas Lie 2 1 Stanford University, Dept. Petroleum Engineering, USA 2 SINTEF IKT, Dept., Norway ECMOR IX, August 30 September
More informationThe simultaneous prediction of vapor-liquid equilibrium and excess enthalpy. Kwon, Jung Hun. Thermodynamics and properties lab.
The simultaneous prediction of vapor-liquid equilibrium and excess enthalpy Kwon, Jung Hun. 2 Contents 1 A comparison of cubic EOS mixing rules for the simultaneous description of excess enthalpies and
More informationPrediction of methanol content in natural gas with the GC-PR-CPA model Hajiw, Martha; Chapoy, Antonin; Coquelet, Christophe; Lauermann, Gerhard
Heriot-Watt University Heriot-Watt University Research Gateway Prediction of methanol content in natural gas with the GC-PR-CPA model Hajiw, Martha; Chapoy, Antonin; Coquelet, Christophe; Lauermann, Gerhard
More informationChemical Potential. Combining the First and Second Laws for a closed system, Considering (extensive properties)
Chemical Potential Combining the First and Second Laws for a closed system, Considering (extensive properties) du = TdS pdv Hence For an open system, that is, one that can gain or lose mass, U will also
More informationrate of reaction forward conc. reverse time P time Chemical Equilibrium Introduction Dynamic Equilibrium Dynamic Equilibrium + RT ln f p
Chemical Equilibrium Chapter 9 of Atkins: Sections 9.1-9.2 Spontaneous Chemical Reactions The Gibbs Energy Minimum The reaction Gibbs energy Exergonic and endergonic reactions The Description of Equilibrium
More informationNumerical Analysis: Solving Systems of Linear Equations
Numerical Analysis: Solving Systems of Linear Equations Mirko Navara http://cmpfelkcvutcz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office
More informationRegularized Least Squares
Regularized Least Squares Ryan M. Rifkin Google, Inc. 2008 Basics: Data Data points S = {(X 1, Y 1 ),...,(X n, Y n )}. We let X simultaneously refer to the set {X 1,...,X n } and to the n by d matrix whose
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2018 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationReactor Design within Excel Enabled by Rigorous Physical Properties and an Advanced Numerical Computation Package
Reactor Design within Excel Enabled by Rigorous Physical Properties and an Advanced Numerical Computation Package Mordechai Shacham Department of Chemical Engineering Ben Gurion University of the Negev
More informationNonlinear FEM. Critical Points. NFEM Ch 5 Slide 1
5 Critical Points NFEM Ch 5 Slide Assumptions for this Chapter System is conservative: total residual is the gradient of a total potential energy function r(u,λ) = (u,λ) u Consequence: the tangent stiffness
More informationB024 RESERVOIR STREAMLINE SIMULATION ACCOUNTING
1 B024 RESERVOIR STREAMLINE SIMULATION ACCOUNTING FOR EFFECTS OF CAPILLARITY AND WETTABILITY R.A. BERENBLYUM, A.A. SHAPIRO, E.H. STENBY IVC-SEP, Department of Chemical Engineering, Technical University
More informationEffect of Sorption/Curved Interface Thermodynamics on Pressure transient
PROCEEDINGS, Twentieth Workshop on Geothermal Rey~volr Englneerlng Stanford Unhrenlty, Stanfoni, Callfornla, January 2426 1995 SGP-m-150 Effect of Sorption/Curved Interface Thermodynamics on Pressure transient
More informationImplementation of the Pertubed-Chain Statistical Association Fluid Theory Model
Implementation of the Pertubed-Chain Statistical Association Fluid Theory Model Bjørn Tore Løvfall Olaf Trygve Berglihn February 20, 2004 Abstract This report represents a workbook from the subject KP8108,
More informationUsing SVD to Recommend Movies
Michael Percy University of California, Santa Cruz Last update: December 12, 2009 Last update: December 12, 2009 1 / Outline 1 Introduction 2 Singular Value Decomposition 3 Experiments 4 Conclusion Last
More informationParameter Norm Penalties. Sargur N. Srihari
Parameter Norm Penalties Sargur N. srihari@cedar.buffalo.edu 1 Regularization Strategies 1. Parameter Norm Penalties 2. Norm Penalties as Constrained Optimization 3. Regularization and Underconstrained
More informationDerivation of the Kalman Filter
Derivation of the Kalman Filter Kai Borre Danish GPS Center, Denmark Block Matrix Identities The key formulas give the inverse of a 2 by 2 block matrix, assuming T is invertible: T U 1 L M. (1) V W N P
More informationOptimization. Benjamin Recht University of California, Berkeley Stephen Wright University of Wisconsin-Madison
Optimization Benjamin Recht University of California, Berkeley Stephen Wright University of Wisconsin-Madison optimization () cost constraints might be too much to cover in 3 hours optimization (for big
More informationPractical Linear Algebra: A Geometry Toolbox
Practical Linear Algebra: A Geometry Toolbox Third edition Chapter 12: Gauss for Linear Systems Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book www.farinhansford.com/books/pla
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationDECOUPLED AND MULTIPHYSICS MODELS FOR NON-ISOTHERMAL COMPOSITIONAL TWO-PHASE FLOW IN POROUS MEDIA
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 9, Number 1, Pages 17 28 c 2012 Institute for Scientific Computing and Information DECOUPLED AND MULTIPHYSICS MODELS FOR NON-ISOTHERMAL COMPOSITIONAL
More informationCMPSCI611: Three Divide-and-Conquer Examples Lecture 2
CMPSCI611: Three Divide-and-Conquer Examples Lecture 2 Last lecture we presented and analyzed Mergesort, a simple divide-and-conquer algorithm. We then stated and proved the Master Theorem, which gives
More informationARock: an algorithmic framework for asynchronous parallel coordinate updates
ARock: an algorithmic framework for asynchronous parallel coordinate updates Zhimin Peng, Yangyang Xu, Ming Yan, Wotao Yin ( UCLA Math, U.Waterloo DCO) UCLA CAM Report 15-37 ShanghaiTech SSDS 15 June 25,
More informationME 680- Spring Geometrical Analysis of 1-D Dynamical Systems
ME 680- Spring 2014 Geometrical Analysis of 1-D Dynamical Systems 1 Geometrical Analysis of 1-D Dynamical Systems Logistic equation: n = rn(1 n) velocity function Equilibria or fied points : initial conditions
More informationSGN Advanced Signal Processing Project bonus: Sparse model estimation
SGN 21006 Advanced Signal Processing Project bonus: Sparse model estimation Ioan Tabus Department of Signal Processing Tampere University of Technology Finland 1 / 12 Sparse models Initial problem: solve
More informationChemistry SEAS Q, K eq, and Equilibrium. K ; where G o = -RTlnK eq eq.
Chemistry 102 - SEAS Q, K eq, and Equilibrium At a given temperature and set of conditions (pressures or concentrations), we can tell if a reaction is already at equilibrium, or which way it will approach
More informationA generalized set of correlations for plus fraction characterization
370 Pet.Sci.(01)9:370-378 DOI 10.1007/s118-01-01-x A generalized set of correlations for plus fraction characterization JAMIALAHMADI Mohamad 1, ZANGENEH Hossein and HOSSEINI Seyed Sajad 1 Petroleum Engineering
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2016 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationQALGO workshop, Riga. 1 / 26. Quantum algorithms for linear algebra.
QALGO workshop, Riga. 1 / 26 Quantum algorithms for linear algebra., Center for Quantum Technologies and Nanyang Technological University, Singapore. September 22, 2015 QALGO workshop, Riga. 2 / 26 Overview
More information7.2 Linear equation systems. 7.3 Linear least square fit
72 Linear equation systems In the following sections, we will spend some time to solve linear systems of equations This is a tool that will come in handy in many di erent places during this course For
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationNew Fast Kalman filter method
New Fast Kalman filter method Hojat Ghorbanidehno, Hee Sun Lee 1. Introduction Data assimilation methods combine dynamical models of a system with typically noisy observations to obtain estimates of the
More informationRegularized Least Squares
Regularized Least Squares Charlie Frogner 1 MIT 2011 1 Slides mostly stolen from Ryan Rifkin (Google). Summary In RLS, the Tikhonov minimization problem boils down to solving a linear system (and this
More informationThe Missing-Index Problem in Process Control
The Missing-Index Problem in Process Control Joseph G. Voelkel CQAS, KGCOE, RIT October 2009 Voelkel (RIT) Missing-Index Problem 10/09 1 / 35 Topics 1 Examples Introduction to the Problem 2 Models, Inference,
More informationFundamentals of Engineering Analysis (650163)
Philadelphia University Faculty of Engineering Communications and Electronics Engineering Fundamentals of Engineering Analysis (6563) Part Dr. Omar R Daoud Matrices: Introduction DEFINITION A matrix is
More informationMonitoring CO 2 Injection at Weyburn Reservoir Using 3-D/3-C Seismic Datasets
Monitoring CO 2 Injection at Weyburn Reservoir Using 3-D/3-C Seismic Datasets Le Gao* and Igor Morozov, University of Saskatchewan, Saskatoon, Saskatchewan le.gao@usask.ca Summary In order to monitor and
More informationLinear Algebra V = T = ( 4 3 ).
Linear Algebra Vectors A column vector is a list of numbers stored vertically The dimension of a column vector is the number of values in the vector W is a -dimensional column vector and V is a 5-dimensional
More informationShortcut Design Method for Columns Separating Azeotropic Mixtures
3908 Ind. Eng. Chem. Res. 2004, 43, 3908-3923 Shortcut Design Method for Columns Separating Azeotropic Mixtures Guilian Liu, Megan Jobson,*, Robin Smith, and Oliver M. Wahnschafft Department of Process
More informationEvaluation of the SPUNG Equation of State for use in Carbon Capture and Storage Modelling
6 th Trondheim CCS-conference June 14-16 2011 Evaluation of the SPUNG Equation of State for use in Carbon Capture and Storage Modelling Geir Skaugen a Øivind Wilhelmsen a, Oddvar Jørstad b, Hailong Li
More informationAvailable online at Energy Procedia 00 (2011) TCCS-6
Available online at www.sciencedirect.com Energy Procedia 00 (2011) 000 000 Energy Procedia www.elsevier.com/locate/procedia TCCS-6 Evaluation of SPUNG # and other Equations of State for use in Carbon
More information