Development of a one-dimensional boiling model: Part I A two-phase flow pattern map for a heavy hydrocarbon feedstock Pieter Verhees, Abdul Akhras Rahman, Kevin M. Van Geem, Geraldine J. Heynderickx Laboratory for Chemical Technology, Ghent University http://www.lct.ugent.be CHISA, Prague, 28/08/2016 1
Fouling of heat exchangers Liquid Heavy oil (Partly) Evaporated oil Tube wall temperature = Fouling Asphaltene precipitation Shutdown Cleaning needed Heat transfer Pressure drop Huge cost 2
Fouling of heat exchangers Liquid Heavy oil (Partly) Evaporated oil Tube wall temperature = Fouling Asphaltene precipitation Shutdown Cleaning needed Heat transfer Pressure drop Huge cost Prediction fouling = 1D modeling 3
Fouling Flow regime Evaporating flow 4
Fouling Flow regime Evaporating flow Stratified flow Annular flow Mist flow Liquid Vapor High tube temperature Low tube temperature Very high tube temperature 5
Fouling Flow regime Evaporating flow Stratified flow Annular flow Mist flow Liquid Vapor High tube temperature Low tube temperature Very high tube temperature Tube temperature Fouling 6
Flow regime! 1D modeling evaporation Fouling prediction Calculation heat transfer coefficient Step 1 Determination flow regime Step 2 HTC calculation Step 1 Fouling model (~T) Step 2 Determination flow regime Step 3 Estimation wall temperature 7
Flow regime! 1D modeling evaporation Fouling prediction Calculation heat transfer coefficient Step 1 Determination flow regime Step 2 HTC calculation Accurate flow regime prediction is important Step 1 Fouling model (~T) Flow regime map Step 2 Determination flow regime Step 3 Estimation wall temperature 8
This work Regime maps in literature are developed based on experiments for Water Cooling fluids In literature correction is proposed to other fluids Valididation of the flow regime maps for heavy hydrocarbons CFD simulations 9
Adiabatic flow regime map Baker (1954) Based on experimental databases Water + air Generalization to other fluids by correction parameters λ = ρ g ρ air ρ l ρ water 0.5 ψ = σ water σ μ l μ water ρ water ρ l 2 1/3 Valid for heavy hydrocarbons? Baker, O., Design of pipe lines for simultaneous flow of oil and gas. Oil and Gas Journal 1954, 26 10
ANSYS FLUENT 13.0 CFD Model Multiphase model: Volume of Fluid Interface tracking: Geo-reconstruct Interface reconstruction Actual interface Piece wise linear interface Accurate when mesh size is an order of magnitude smaller than radius of curvature Validated with Water + Air simulations 11
CFD simulations Feed Gasoil is complex hydrocarbon mixture Represented by 1 pseudocomponent Simulation conditions 7 simulations are performed Operating conditions All seven regimes are observed W A M Sl B St P 12
CFD results W A M B Sl Wavy St P Plug Stratified Slug Liquid Vapor De Schepper, S. C. K.; Heynderickx, G. J.; Marin, G. B., CFD modeling of all gas-liquid and vapor-liquid flow regimes predicted by the Baker chart. Chemical Engineering Journal 2008, 138, (1-3), 349-357 13
CFD results Mist Annular Bubble W A M B Sl St P Liquid Vapor 14
CFD results Mist Annular Bubble Evaporation has an effect on flow regimes Diabatic flow regime map W A M B Sl St P Liquid Vapor 15
Kattan, Thome and Favrat (1998) Diabatic flow regime map Based on experimental databases cooling fluid R134a Correlations are developed to construct the boundaries Boundaries dependent on Geometry Fluid property Wall heat flux Correlations valid for heavy hydrocarbons? Kattan, N.; Thome, J. R.; Favrat, D., Flow boiling in horizontal tubes: Part 1 - Development of a diabatic two-phase flow pattern map. Journal of Heat Transfer-Transactions of the Asme 1998, 120, (1), 140-147 16
CFD simulations CFD model adaptation Phase change model Evaporation (T > T sat ) Condensation (T < T sat ) coeff = 0.1 s -1 Simulation conditions 370 kg/(m 2 s) 17
CFD results 370 kg/(m 2 s) 18
CFD results Annular flow is not observed Grid adaption is applied Vapor/liquid Interface Refine mesh 370 kg/(m 2 s) Liquid Vapor 19
CFD results Grid adaption 370 kg/(m 2 s) Liquid Vapor 20
CFD results Grid adaption No grid adaptation Grid adaptation Liquid Vapor 21
CFD results Grid adaption No grid adaptation Annular flow still not stable Nucleation at the wall not modeled Lack of nucleation models Grid adaptation Accurate regime map? Liquid Vapor 22
CFD results Grid adaption No grid adaptation Annular flow not stable Wall is less wetted as expected Heat transfer coefficient Wall temperature Grid adaptation Liquid Vapor 23
Conclusions & future work Adiabatic CFD model was validated by Baker chart Diabatic CFD simulations did not completely agree with the constructed regime map Proper regime map? Nucleation at the wall 24
acknowledgement Long Term Structural Methusalem Funding by the Flemish Government UGent High-Performance Computing Stevin 25