2
Proceedings of the 12
thInternational Conference on Computational Fluid Dynamics in the Oil & Gas, Metallurgical and Process Industries
SINTEF
PROCEEDINGS
Progress in Applied CFD –
CFD2017
Editors:
Jan Erik Olsen and Stein Tore Johansen
Progress in Applied CFD – CFD2017
Proceedings of the 12
thInternational Conference on Computational Fluid Dynamics in the Oil & Gas, Metallurgical and Process Industries
SINTEF Proceedings
SINTEF Academic Press
SINTEF Proceedings no 2
Editors: Jan Erik Olsen and Stein Tore Johansen Progress in Applied CFD – CFD2017
Selected papers from 10th International Conference on Computational Fluid Dynamics in the Oil & Gas, Metal lurgical and Process Industries
Key words:
CFD, Flow, Modelling
Cover, illustration: Arun Kamath
ISSN 2387-4295 (online) ISBN 978-82-536-1544-8 (pdf)
© Copyright SINTEF Academic Press 2017
The material in this publication is covered by the provisions of the Norwegian Copyright Act. Without any special agreement with SINTEF Academic Press, any copying and making available of the material is only allowed to the extent that this is permitted by law or allowed through an agreement with Kopinor, the Reproduction Rights Organisation for Norway. Any use contrary to legislation or an agreement may lead to a liability for damages and confiscation, and may be punished by fines or imprisonment
SINTEF Academic Press
Address: Forskningsveien 3 B PO Box 124 Blindern
N-0314 OSLO
Tel: +47 73 59 30 00 Fax: +47 22 96 55 08 www.sintef.no/byggforsk www.sintefbok.no
SINTEF Proceedings
SINTEF Proceedings is a serial publication for peer-reviewed conference proceedings on a variety of scientific topics.
The processes of peer-reviewing of papers published in SINTEF Proceedings are administered by the conference organizers and proceedings editors. Detailed procedures will vary according to custom and practice in each scientific community.
PREFACE
This book contains all manuscripts approved by the reviewers and the organizing committee of the 12th International Conference on Computational Fluid Dynamics in the Oil & Gas, Metallurgical and Process Industries. The conference was hosted by SINTEF in Trondheim in May/June 2017 and is also known as CFD2017 for short. The conference series was initiated by CSIRO and Phil Schwarz in 1997.
So far the conference has been alternating between CSIRO in Melbourne and SINTEF in Trondheim.
The conferences focuses on the application of CFD in the oil and gas industries, metal production, mineral processing, power generation, chemicals and other process industries. In addition pragmatic modelling concepts and bio‐mechanical applications have become an important part of the conference. The papers in this book demonstrate the current progress in applied CFD.
The conference papers undergo a review process involving two experts. Only papers accepted by the reviewers are included in the proceedings. 108 contributions were presented at the conference together with six keynote presentations. A majority of these contributions are presented by their manuscript in this collection (a few were granted to present without an accompanying manuscript).
The organizing committee would like to thank everyone who has helped with review of manuscripts, all those who helped to promote the conference and all authors who have submitted scientific contributions. We are also grateful for the support from the conference sponsors: ANSYS, SFI Metal Production and NanoSim.
Stein Tore Johansen & Jan Erik Olsen
3
Organizing committee:
Conference chairman: Prof. Stein Tore Johansen Conference coordinator: Dr. Jan Erik Olsen Dr.
Bernhard Müller
Dr.Sigrid Karstad Dahl Dr.Shahriar Amini Dr.Ernst Meese Dr.Josip Zoric Dr.Jannike Solsvik Dr.Peter Witt
Scientific committee:
Stein Tore Johansen, SINTEF/NTNU Bernhard Müller, NTNU
Phil Schwarz, CSIRO
Akio Tomiyama, Kobe University
Hans Kuipers, Eindhoven University of Technology Jinghai Li, Chinese Academy of Science
Markus Braun, Ansys Simon Lo, CD‐adapco
Patrick Segers, Universiteit Gent Jiyuan Tu, RMIT
Jos Derksen, University of Aberdeen Dmitry Eskin, Schlumberger‐Doll Research Pär Jönsson, KTH
Stefan Pirker, Johannes Kepler University Josip Zoric, SINTEF
4
CONTENTS
PRAGMATIC MODELLING ... 9
On pragmatism in industrial modeling. Part III: Application to operational drilling ... 11
CFD modeling of dynamic emulsion stability ... 23
Modelling of interaction between turbines and terrain wakes using pragmatic approach ... 29
FLUIDIZED BED ... 37
Simulation of chemical looping combustion process in a double looping fluidized bed reactor with cu‐based oxygen carriers ... 39
Extremely fast simulations of heat transfer in fluidized beds ... 47
Mass transfer phenomena in fluidized beds with horizontally immersed membranes ... 53
A Two‐Fluid model study of hydrogen production via water gas shift in fluidized bed membrane reactors ... 63
Effect of lift force on dense gas‐fluidized beds of non‐spherical particles ... 71
Experimental and numerical investigation of a bubbling dense gas‐solid fluidized bed ... 81
Direct numerical simulation of the effective drag in gas‐liquid‐solid systems ... 89
A Lagrangian‐Eulerian hybrid model for the simulation of direct reduction of iron ore in fluidized beds... 97
High temperature fluidization ‐ influence of inter‐particle forces on fluidization behavior ... 107
Verification of filtered two fluid models for reactive gas‐solid flows ... 115
BIOMECHANICS ... 123
A computational framework involving CFD and data mining tools for analyzing disease in cartoid artery ... 125
Investigating the numerical parameter space for a stenosed patient‐specific internal carotid artery model ... 133
Velocity profiles in a 2D model of the left ventricular outflow tract, pathological case study using PIV and CFD modeling ... 139
Oscillatory flow and mass transport in a coronary artery ... 147
Patient specific numerical simulation of flow in the human upper airways for assessing the effect of nasal surgery ... 153
CFD simulations of turbulent flow in the human upper airways ... 163
OIL & GAS APPLICATIONS ... 169
Estimation of flow rates and parameters in two‐phase stratified and slug flow by an ensemble Kalman filter ... 171
Direct numerical simulation of proppant transport in a narrow channel for hydraulic fracturing application ... 179
Multiphase direct numerical simulations (DNS) of oil‐water flows through homogeneous porous rocks ... 185
CFD erosion modelling of blind tees ... 191
Shape factors inclusion in a one‐dimensional, transient two‐fluid model for stratified and slug flow simulations in pipes ... 201
Gas‐liquid two‐phase flow behavior in terrain‐inclined pipelines for wet natural gas transportation ... 207
NUMERICS, METHODS & CODE DEVELOPMENT ... 213
Innovative computing for industrially‐relevant multiphase flows ... 215
Development of GPU parallel multiphase flow solver for turbulent slurry flows in cyclone ... 223
Immersed boundary method for the compressible Navier–Stokes equations using high order summation‐by‐parts difference operators ... 233
Direct numerical simulation of coupled heat and mass transfer in fluid‐solid systems ... 243
A simulation concept for generic simulation of multi‐material flow, using staggered Cartesian grids ... 253
A cartesian cut‐cell method, based on formal volume averaging of mass, momentum equations ... 265
SOFT: a framework for semantic interoperability of scientific software ... 273
POPULATION BALANCE ... 279
Combined multifluid‐population balance method for polydisperse multiphase flows ... 281
A multifluid‐PBE model for a slurry bubble column with bubble size dependent velocity, weight fractions and temperature ... 285
CFD simulation of the droplet size distribution of liquid‐liquid emulsions in stirred tank reactors ... 295
Towards a CFD model for boiling flows: validation of QMOM predictions with TOPFLOW experiments ... 301
Numerical simulations of turbulent liquid‐liquid dispersions with quadrature‐based moment methods ... 309
Simulation of dispersion of immiscible fluids in a turbulent couette flow ... 317
Simulation of gas‐liquid flows in separators ‐ a Lagrangian approach ... 325
CFD modelling to predict mass transfer in pulsed sieve plate extraction columns ... 335
BREAKUP & COALESCENCE ... 343
Experimental and numerical study on single droplet breakage in turbulent flow ... 345
Improved collision modelling for liquid metal droplets in a copper slag cleaning process ... 355
Modelling of bubble dynamics in slag during its hot stage engineering ... 365
Controlled coalescence with local front reconstruction method ... 373
BUBBLY FLOWS ... 381
Modelling of fluid dynamics, mass transfer and chemical reaction in bubbly flows ... 383
Stochastic DSMC model for large scale dense bubbly flows ... 391
On the surfacing mechanism of bubble plumes from subsea gas release ... 399
Bubble generated turbulence in two fluid simulation of bubbly flow ... 405
HEAT TRANSFER ... 413
CFD‐simulation of boiling in a heated pipe including flow pattern transitions using a multi‐field concept ... 415
The pear‐shaped fate of an ice melting front ... 423
Flow dynamics studies for flexible operation of continuous casters (flow flex cc) ... 431
An Euler‐Euler model for gas‐liquid flows in a coil wound heat exchanger ... 441
NON‐NEWTONIAN FLOWS ... 449
Viscoelastic flow simulations in disordered porous media ... 451
Tire rubber extrudate swell simulation and verification with experiments ... 459
Front‐tracking simulations of bubbles rising in non‐Newtonian fluids ... 469
A 2D sediment bed morphodynamics model for turbulent, non‐Newtonian,
particle‐loaded flows ... 479
METALLURGICAL APPLICATIONS ... 491
Experimental modelling of metallurgical processes ... 493
State of the art: macroscopic modelling approaches for the description of multiphysics phenomena within the electroslag remelting process ... 499
LES‐VOF simulation of turbulent interfacial flow in the continuous casting mold ... 507
CFD‐DEM modelling of blast furnace tapping ... 515
Multiphase flow modelling of furnace tapholes ... 521
Numerical predictions of the shape and size of the raceway zone in a blast furnace ... 531
Modelling and measurements in the aluminium industry ‐ Where are the obstacles? ... 541
Modelling of chemical reactions in metallurgical processes ... 549
Using CFD analysis to optimise top submerged lance furnace geometries ... 555
Numerical analysis of the temperature distribution in a martensic stainless steel strip during hardening ... 565
Validation of a rapid slag viscosity measurement by CFD ... 575
Solidification modeling with user defined function in ANSYS Fluent ... 583
Cleaning of polycyclic aromatic hydrocarbons (PAH) obtained from ferroalloys plant ... 587
Granular flow described by fictitious fluids: a suitable methodology for process simulations ... 593
A multiscale numerical approach of the dripping slag in the coke bed zone of a pilot scale Si‐Mn furnace ... 599
INDUSTRIAL APPLICATIONS ... 605
Use of CFD as a design tool for a phospheric acid plant cooling pond ... 607
Numerical evaluation of co‐firing solid recovered fuel with petroleum coke in a cement rotary kiln: Influence of fuel moisture ... 613
Experimental and CFD investigation of fractal distributor on a novel plate and frame ion‐exchanger ... 621
COMBUSTION ... 631
CFD modeling of a commercial‐size circle‐draft biomass gasifier ... 633
Numerical study of coal particle gasification up to Reynolds numbers of 1000 ... 641
Modelling combustion of pulverized coal and alternative carbon materials in the blast furnace raceway ... 647
Combustion chamber scaling for energy recovery from furnace process gas: waste to value ... 657
PACKED BED ... 665
Comparison of particle‐resolved direct numerical simulation and 1D modelling of catalytic reactions in a packed bed ... 667
Numerical investigation of particle types influence on packed bed adsorber behaviour ... 675
CFD based study of dense medium drum separation processes ... 683
A multi‐domain 1D particle‐reactor model for packed bed reactor applications ... 689
SPECIES TRANSPORT & INTERFACES ... 699
Modelling and numerical simulation of surface active species transport ‐ reaction in welding processes ... 701
Multiscale approach to fully resolved boundary layers using adaptive grids ... 709
Implementation, demonstration and validation of a user‐defined wall function for direct precipitation fouling in Ansys Fluent ... 717
FREE SURFACE FLOW & WAVES ... 727 Unresolved CFD‐DEM in environmental engineering: submarine slope stability and
other applications... 729 Influence of the upstream cylinder and wave breaking point on the breaking wave
forces on the downstream cylinder ... 735 Recent developments for the computation of the necessary submergence of pump
intakes with free surfaces ... 743 Parallel multiphase flow software for solving the Navier‐Stokes equations ... 752
PARTICLE METHODS ... 759 A numerical approach to model aggregate restructuring in shear flow using DEM in
Lattice‐Boltzmann simulations ... 761 Adaptive coarse‐graining for large‐scale DEM simulations ... 773 Novel efficient hybrid‐DEM collision integration scheme ... 779 Implementing the kinetic theory of granular flows into the Lagrangian
dense discrete phase model ... 785 Importance of the different fluid forces on particle dispersion in fluid phase
resonance mixers ... 791 Large scale modelling of bubble formation and growth in a supersaturated liquid ... 798
FUNDAMENTAL FLUID DYNAMICS ... 807 Flow past a yawed cylinder of finite length using a fictitious domain method ... 809 A numerical evaluation of the effect of the electro‐magnetic force on bubble flow
in aluminium smelting process ... 819
A DNS study of droplet spreading and penetration on a porous medium ... 825
From linear to nonlinear: Transient growth in confined magnetohydrodynamic flows ... 831
12th International Conference on CFD in Oil & Gas, Metallurgical and Process Industries SINTEF, Trondheim, NORWAY
May 30th - June 1st, 2017
CFD 2017
AN EULER-EULER MODEL FOR GAS-LIQUID FLOWS IN A COIL WOUND HEAT EXCHANGER
Thomas ACHER1∗, Manuel KNAUP1, Konrad BRAUN1, Hans-Jörg ZANDER1
1The Linde Group - Engineering Division, Dr.-Carl-von-Linde-Straße 6-14, 82049 Pullach, GERMANY
∗E-mail: thomas.acher@linde.com
ABSTRACT
Coil-wound heat exchangers (CWHE) are commonly adapted in process engineering for the efficient transfer of heat between fluids which feature wide temperature and pressure ranges. The field of application for this apparatus ranges from heating or cooling of single-phase flows, over the evaporation or condensation of fluids, to the utilization as isothermal reactor. Due to their large specific heat transfer area accompanied by a compact design, coil-wound heat exchangers are widely used in various process plants (e.g., LNG plants). Depending on the application, two- phase flows may occur at both, the tube- as well as the shell- side of the apparatus. For the design of a CWHE, the fluid and thermodynamic processes in the unit are commonly represented by a system of one-dimensional correlations. This approach implies uniform thermohydraulic conditions on horizontal cutting planes of the exchanger. Fluid and thermodynamic effects in the apparatus which result in radial parameter variations are inaccessible to these conventional design tools. To this end, a multidimensional CFD model has been established to enhance the representation of fluid and thermodynamic phenomena in CWHE design. The shell- side of the CWHE and all tube-side sections are each numerically represented by separate domains which are coupled by source terms to account for the thermodynamic interaction between tube- and shell-side. In each flow region, the hydraulic effect of the tube bundle is modeled as a porous medium with corresponding fluid dynamic characteristics. The gas-liquid dynamics in each flow region is modeled based on an Euler-Euler approach. Unlike classical Euler-Euler models, local phase fractions and fluid properties are calculated from species relations as well as pressure and temperature fields. This model framework is augmented by locally evaluated correlations for pressure drop and heat transfer to account for apparatus internals and thermal coupling. The models for gas-liquid interaction forces are derived from standard correlations and augmented by findings from detailed CFD studies.
Remaining parameters are specified by a parameterization study based on experimental findings.
Keywords: process industry, heat exchanger, multiphase heat transfer, Euler-Euler approach.
NOMENCLATURE Greek Symbols
α Phase volume fraction,[−]
ε Porosity,[−]
ρ Mass density,[kg/m3] τ Stress tensor,[N/m3] Latin Symbols
a Velocity profile coefficient,[−].
A Area,[m2].
C Model constant,[−].
~f Volumetric force,[kg/ms2].
g Gravitational acceleration,[m/s2].
p Pressure,[Pa].
~u Velocity,[m/s].
Sub/superscripts ax Axial.
cent Centrifugal force.
g Gas.
h Homogeneous.
i Phase indicator.
j Phase indicator.
l Liquid.
li f t Lift force.
rad Radiali.
INTRODUCTION
In various process engineering applications coil-wound heat exchangers (CWHE) are adopted for the efficient transfer of heat between fluids. The potential operating conditions of this type of apparatus feature a wide range of temperature and pressure levels and the transfer of heat between multiple fluid streams can be realized in one heat exchanger. CWHEs are specifically robust to transients of flow properties (Pacio and Dorao, 2011) and are characterized by a large specific heat transfer area per volume in combination with a comparably compact design. In a CWHE, multiple layers of tubes are wound helically around a central pipe (see Fig. 1) where several tubes are grouped into sections with varying fluid streams. Apparatuses with up to 10,000 tubes enable complex heat exchange processes with heating surfaces of 30,000 m2 and more. Maximal bundle diameters of approx. 5,500 mm and unit heights up to 20 meters can be manufactured. (Walteret al., 2014)
In many large scale plants for the production of liquified natural gas (LNG) CWHEs represent a key equipment. In these applications, the gaseous natural gas is fed to tube sections at the bottom of the apparatus and is cooled, liquefied or subcooled while streaming upwards. In the 441
T. Acher, M. Knaup, K. Braun, H. Zander remaining tube sections the refrigerant is precooled by
streaming concurrent to the natural gas. Subsequently, the refrigerant is applied to the top of the shell-side as a dispersed two-phase stream which evaporates while streaming towards the bottom of the CWHE. The natural gas liquification is driven by the complex thermodynamic interaction across the tube bundle with the multiphase flow at the shell-side of the CWHE. The gas-liquid flow characteristics determine local heat transfer and affect the efficiency and functionality of the entire apparatus.
tube-side inlet
tube-side outlet
shell-side outlet shell-side inlet
Figure 1:Exemplary sketch of a coil-wound heat exchanger.
In general, one-dimensional correlations are applied to represent fluid and thermodynamic processes in order to design CWHE apparatuses (Wanget al., 2015). Merely axial variations of thermohydraulic parameters are captured by this approach while radial parameter variations, specifically fluiddynamic inhomogeneities, are neglected. In order to incorporate these potentially influential phenomena a two- dimensional CFD model has been developed. The following sections explicate this approach, specifically the gas-liquid flow representation as well as the effect of radial varying geometry parameters at the shell-side. Remaining model constants are specified by parameterization and CFD model results are presented.
CFD MODEL DESCRIPTION
A numerically efficient representation of the CWHE is realized by the consideration of a two-dimensional cross- section of the cylindrical form (see Fig. 2). While assuming rotational symmetry of the apparatus, radial variations of thermohydraulic parameters are represented by this modeling setup. The individual fluid dynamic regions on the shell- side and the tube-side of the CWHE are reflected by separate simulation domains, where multiple tube sections are represented individually. To model the heat transfer between the shell- and the tube-side the independent flow regions are thermodynamically coupled by respective source terms and the hydraulic effect of the tube bundle is represented as a porous medium with corresponding fluid dynamic characteristics. This setup enables a coherent decomposition of the numerical meshes of the individual domains to allow for an efficient computation on multiple processors.
shell-side tube-side
z }| {
section 1 section 2 section N
Q˙
Figure 2:Two-dimensional representation of the CWHE with a shell-side domain and multiple tube sections (Göllet al., 2013).
The fluid dynamics at both the tube- as well as the shell- side of the CWHE are modeled by an Euler-Euler approach (e.g., (Drew and Passman, 1999; Ishii and Hibiki, 2010)) since two-phase flows may occur at any flow region. The respective transport equations for momentum read
∂
∂tαiρi~ui+∇·αiρi~ui~ui=
−αi∇pi+∇·αiτi+αiρi~g+~frad,i+~fh,i (1) whereαi,ρi,~uiandτi denote the volume fraction, density, velocity and shear stress tensor of phase i. The fluid dynamic interaction between gas and liquid as well as the flow resistance in the porous media are modeled by the phase specific body forces in axial direction fh,iand radial direction
frad,i.
The Euler-Euler CFD solver twoPhaseEulerFoam of the program package OpenFOAM represents the basis for the model implementation which is interfaced to submodels for the evaluation of local void fractions and fluid parameters, pressure drop in the porous media and transferred heat between flow regions (see Fig. 3). In contrast to classical Euler-Euler models, local phase fractions are computed from the local composition of chemical species and thermodynamic conditions which requires the inclusion of transport equations of the comprised species instead of a respective formulation for local void fraction.
Local thermodynamic equilibrium is assumed to determine relevant fluid properties (i.e. viscosity, density, thermal 442
An Euler-Euler Model for Gas-Liquid Flows in a Coil Wound Heat Exchanger/ CFD 2017 conductivity and surface tension) and phase composition in
each numerical cell based on local species concentrations and thermodynamic relations (Jakobsen, 2008). Local thermodynamic conditions are computed from the energy transport equation which includes source terms for the heat transfer across the simulation domains.
As depicted in Fig. 3 the evaluation of fluid parameters at vapor/liquid-equilibrium in each numerical cell is realized by coupling the CFD solver to Linde’s inhouse software. This program is used for cell-wise conduction of flash calculations (pH-flash) in order to evaluate local phase fractions and fluid properties. Energy transport equation source terms which represent heat transfer effects are based on one-dimensional correlations used for CWHE design at Linde which are determined by flow parameters and single- or two-phase flow conditions. At the shell-side as well as at the tube- side, the flow resistance due to the tube bundle is integrated in the CFD model by a porous zone. The corresponding pressure drop properties are evaluated from CWHE design correlations analogous to the heat transfer calculations. For additional information on the CFD model framework the reader is referred to previous publications (Acher et al., 2016).
heat transfer
pressure drop fluid
properties
CFD solver (steady state) Q˙
p,~u,e
∆p|irr local thermo-
dynamic equilibrium
Figure 3:Depiction of numerical model components and their interaction.
As mentioned above, potential radial variations of fluid dynamic parameters might impact the heat transfer performance of the CWHE. These non-uniformities can originate from two-phase flow effects on the shell-side of the apparatus. The consideration of these phenomena necessitates the incorporation of submodels for an enhanced depiction of gas-liquid flow conditions.
Two-phase flow model development
The shell-side of the CWHE features a gas-liquid flow which is applied as a dispersed stream above the tube bundle and evaporates on its way to the bottom of the apparatus. The related thermohydraulic effects are commonly depicted by empirical correlations which resort to pseudo-homogenous parameters describing the multiphase flow. In contrast, in an Euler-Euler context the prevalent retaining and interfacial forces have to be modeled individually to yield phase- specific velocity fields (see Fig. 4).
In order to ensure consistency with the pseudo-homogeneous concept for a radially uniform flow field, the two- phase pressure drop characteristics along the bundle are
fh fh,l fh,g
liquid gaseous liquid gaseous
Figure 4:Schematic depiction of the volume force acting on fluid elements in a numerical cell (Knaup, 2015). Left:
pseudo-homogenous flow resistance. Right: phase- specific flow resistance.
retained as specified by the one-dimensional model and the homogeneous force field fh is divided into phase-specific contributionsfh,i:
fh,i=Cifh= fax,i
fax,j+fax,ifh (2) where the weighting coefficientCiis a function of the axial forces acting on the gas phase fax,gand the liquid phase fax,l inside the tube bundle.
The axial retaining and interfacial forces for the individual phases fax,gand fax,lare further subdivided according to the bilateral interaction between gas, liquid and tube bundle. To do so, a basic perception of the gas-liquid flow pattern in the tube bundle is required. Figure 5 depicts the assumed liquid flow structures on the shell-side of the tube bundle which are comprised of vertical films enclosing the coiled tubes and dispersed droplets between tube layers.
Figure 5:Sketch of the presumed gas-liquid flow conditions at the shell-side of the CWHE bundle
Due to their fundamentally different characteristics with regard to gas-liquid interaction the liquid film and droplets are modeled independently. Hence, the axial phase-specific forces at the shell-side are formulated as
fax,g= fgs+fgl,d+fgl,f, (3) 443
T. Acher, M. Knaup, K. Braun, H. Zander
fax,l=fls−fgl,d−fgl,f. (4) Table 1 gives a description on each of the phase-specific interfacial forces and states the modeling approach selected for their evaluation. As the gas-liquid interaction fgl,d is a function of the volume fraction of liquid droplets and drop sizes these parameters are derived from model correlations based on local fluid dynamic equilibrium considerations (Sirignano, 2010; Ishii and Mishima, 1989). Consistent to experimental observations, it is assumed that liquid velocities are dominated by the film around the tubes as the volumetric contribution of droplets is considered subordinate. This conception arguably justifies the adoption of simple approaches for the evaluation of droplet sizes and local volume ratios of droplet to film. Yet, this assumption is further validated by detailed 3D CFD simulations presented below.
Analogous to the representation of phase-specific axial forces the numerical model includes a formulation for radial volume forces frad,iin the Euler-Euler momentum equations Eqn. 1. Two prevalent mechanisms which induce radially acting volume forces were identified based on experimental findings, operational data analysis and detailed 3D CFD simulation of single- and multiphase flows. Table 2 gives a brief description of both the lift and centrifugal force. A more detailed explication with a focus on the origin of both effects with regard to the apparatus at hand is part of the following subsections.
The phase-specific cumulated effect of radial fluid dynamic processes is incorporated as
frad,l=fli f t,l+fcent,l, (5)
frad,g=−fli f t,g+fcent,g. (6) The formulations for both radially acting forces fli f t,i and fcent,i feature a model constant which has to be adjusted to experimental findings and is further explicated in the last section of this work.
The subdivision of retaining and interfacial forces fh,i, explicated above, is a prerequisite to the evaluation of radial fluid dynamic processes as the respective models for fli f t,i and fcent,idepend on individual phase velocitiesui. Locally determined phase-specific velocities which are inaccessible to common one-dimensional CWHE design correlations are required to investigate radial inhomogeneity of the two- phase flow on the shell-side of the apparatus. Radial parameter variations in the CWHE can be ascribed to radial inhomogeneity of the two-phase flow on the shell-side which impacts the efficiency of heat transfer processes.
Incorporation of the centrifugal force
The helical motion of gas in between the coiled tube rows as sketched in Fig. 6 effects a centrifugal force which is directed radially outward independent of the tube layer coiling direction. The gaseous flow pattern was confirmed by three-dimensional numerical simulations of single-phase flow through the resolved bundle geometry. It is assumed that liquid droplets follow this spiral motion (analogous to the flow in a cyclone used as a separator) and consequently are subject to the centrifugal force in radial direction.
Incorporation of the lift force
The lift force effect is initiated from radial variation of the axial gas flow velocity due to the geometrical inhomogeneity
Figure 6:Delineation of the helical gas flow alternating in orientation for each tube layer. The arrows indicate the direction of the gas motion.
of the tube bundle. To illustrate this connection, a section of a CWHE tube bundle geometry for an experimental setup is depicted in Fig. 7. It shows several layers of coil-wound tubes in between the central pipe (i.e. the mandrel) and the outer shroud of the apparatus. A constant distance between tube rows and layers is maintained by tube spacers which are depicted in Fig. 7 in an abstracted manner by vertical beams.
At different radii the spacers block varying portions of the free space between two layers of wound tubes. This results in a radial dependency of shell-side flow resistance which is to be adopted in the numerical modeling.
The non-uniform blocking of the free area between tube layers is represented in the CFD model via radially variant values of porosity defined as the ratio of actual void area to void area without any spacers:
ε=Avoid−Aspacer
Avoid (7)
Using this values of porosity a representative gas velocity for a heat exchanger without spacers can be derived from the local gas velocity:
ug,∞= 1
ε a
·ug (8)
This velocity is the input for the calculation of the homogenous flow resistance fhin equation 2 and leads to a radially varying flow resistance generating a velocity profile which can be adjusted by the coefficient a in equation 8.
The parameter a is deduced from single phase 3D CFD simulations of the shell side flow through a representative section with detailed resolution of the tube bundle geometry as explained in the following section.
Submodel validation and parametrization
The aforementioned setup for the calculation of the gas- liquid flow through a coil wound heat exchanger requires many different sub-models, whose accuracy is crucial for 444
An Euler-Euler Model for Gas-Liquid Flows in a Coil Wound Heat Exchanger/ CFD 2017
Table 1:Considered phase-specific axial forces on the two-phase flow at the shell-side of the CWHE bundle
Description Evaluation method
Gas-solid interaction fgs Retaining force on the gas flow by the partial blocking of the flow cross-section
Evaluation from single phase pressure drop correlation (Steinbauer and Hecht, 1996) Liquid-solid interaction fls Liquid film flow resistance
due to tube wall friction
Adaption of the Nusselt water skin theory to flow across tubes (Rogers, 1981) Gas-liquid droplet interactionfgl,d Drag force on droplets
induced by the gaseous flow
Calculation from conventional
drag force model (Schiller and Naumann, 1935) Gas-liquid film interactionfgl,f Shear forces at the fluidic interface
between gas and liquid film
Effect is considered as subordinate and therefore neglected
Table 2:Considered radial forces on the two-phase flow at the shell-side of the CWHE bundle
Description Evaluation method
Droplet lift force fli f t,i Radial force on the droplet due to inhomogeneous gas flow through the porosity
Determination from common lift force approach with an adapted coefficient (Ishii and Hibiki, 2010) Centrifugal force fcent,i Centrifugal effect on liquid droplets
due to rotational motion of the gas
Definition of the centrifugal force on the droplet with a model factor
Figure 7:CAD model of a section of an experimental CWHE tube bundle with highlighted tube spacers. (Kiewat, 2015)
the correct prediction of the heat transfer inside the bundle.
These models were derived under several assumptions and so have to be validated and parametrized with numerical and experimental data.
For the numerical investigation a CFD simulation resolving the gas-liquid and liquid-solid interface and most of the turbulence spectrum is required. In this case, the Volume- of-Fluid (VOF) method for the two-phase flow is coupled with an LES approach for the description of turbulence.
With this approach the phase interaction and the effects of turbulence are derived directly from the Navier-Stokes- Equations instead of physical models.
The VOF approach is accompanied by a high computational effort to resolve most of the relevant length scales. In this case the smallest cell size was chosen to have a dimension of about 10 micrometers to resolve the vast majority of expected droplets. The total number of cells is decreased by limiting the simulation domain to a representative cutout of a coil wound heat exchanger using periodic boundary conditions as depicted in Figure 8 for one tube layer.
Figure 8:Schematic depiction of the periodicities in one layer of a coil wound heat exchanger (Bassfeld, 2017). Left:
periodicity in axial direction; Right: periodicity in azimuthal direction.
To account for the interaction within the coil three layers were considered in the simulation setup. To this end, the non-overlapping boundaries resulting from the different orientated slope in the layers were associated with periodic boundary conditions as can be seen in Figure 9. As only fluiddynamical aspects were examined here, the flow was 445
T. Acher, M. Knaup, K. Braun, H. Zander regarded as incompressible without heat transfer from the
tube-side or vaporization of the shell-side flow.
Figure 9:Schematic depiction of the periodicities in radial direction. (Bassfeld, 2017). Same colour indicates periodic connection. For better comprehension also one section above and below is depicted in transparent.
Figure 10 shows an instantaneous depiction of the liquid phase distribution which is color-coded by the local velocity.
This confirms the current perception of gas-liquid pattern in the coil as shown in Figure 5. Especially the droplets’
volume fraction of the whole liquid could be determined to about 0.05, justifying the rather simple approach for the droplet interaction in the previous section.
velocity (m/s)
Figure 10:Isometric view of the liquid flow inside the tube bundle.
(Bassfeld, 2017).
Most of the liquid is present as a falling film while the reminder forms droplets which are generated by liquid separation from film waves. These droplets are transported away from the tube layer and cause the radial migration of liquid. The accumulated effect is captured by the radial fluid interaction submodels in the Euler-Euler context in an abstracted manner.
The helical gas motion as sketched in Figure 6 is confirmed by the averaged velocity field resulting from the VOF
simulation. Figure 11 depicts the respective azimuthal component of the time-averaged gas flow indicating the swirl flow pattern.
velocity (m/s)
Figure 11:Plot of the azimuthal velocity component of the gas flow. (Bassfeld, 2017).
To determine the coefficients for the radial profile of the axial gas velocity a single phase RANS simulation was performed.
The resulting profile is depicted in Figure 12. By adapting the coefficientain equation 8 to a value of 2.5 this velocity profile could be reproduced in the Euler-Euler model.
Normalizedaxialvelocity(-)
Normalized radius (-)
Figure 12:Plot of the axial velocity component of the gas flow over the radius. Local velocity values are normalized by average axial velocity.
To supplement the numerical data also an experimental campaign was launched using a pilot scale CWHE (see (Walteret al., 2014) for details). It completes the validation data set for the measurement of the radial liquid distribution for many different flow conditions with and without heat transfer. These results were used for a parametrization study in order to specify model constants in the definition of the radial forces fli f t,iand fcent,i
MODEL RESULTS
The final model setup was then used to simulate various experimental cases and the measured radial variation of the liquid distribution below the tube bundle was used to validate the CFD model setup.
The simulation results of the numerical model with a fixed parameter set show good accordance with the experimental findings for the entire spectrum of investigated fluid dynamic 446
An Euler-Euler Model for Gas-Liquid Flows in a Coil Wound Heat Exchanger/ CFD 2017 parameters. As an example, a comparison of numerical
and experimental results for two cases with strongly varying vapor quality is shown in Figures 13 and 14.
Experiment Simulation Uniform
Outer Middle Inner
Figure 13:Liquid flow distribution at the bottom of the tube bundle for an experimental case with high gas throughput.
Experiment Simulation Uniform
Outer Middle Inner
Figure 14:Liquid flow distribution at the bottom of the tube bundle for an experimental case with low gas throughput.
CONCLUSION
A two-dimensional CFD model was introduced in the present work to represent the thermo- and fluid-dynamic processes in a coil-wound heat exchanger. Specifically the representation of the gas-liquid hydrodynamics is explicated with a focus on the incorporation of retaining and interfacial forces.
The model framework represents a modified Euler-Euler approach adopting the concept of local thermodynamic equilibrium to determine volume fractions and phase properties. The local heat transfer across the exchanger fluid regions is represented by one-dimensional model formulations which are commonly used for coil-wound heat exchanger design. These thermohydraulic correlations also determine the characteristics of the porous media model to consider the impact of the tube bundle on the fluid dynamics through and within the tube bundle.
The representation of individual retaining and interfacial forces at the shell-side of the tube bundle enables the computation of phase-specific velocities. In association with the consideration of radial geometry variations the differences in the phase velocities effect a radial motion of the liquid. This effect is superimposed by a centrifugal effect which is induced by the helical motion of the gaseous phase phase following the coiled tubes. Both fluid- dynamic processes are represented in the numerical model by the incorporation of corresponding source terms in the momentum equations.
Detailed 3D CFD simulation and experimental investigations have been conducted in order to validate specifics of phase interaction models. The Volume-of-fluid simulations indicated local liquid flow conditions, specifically with regard to droplet sizes, quality and behavior. Additionally these investigations confirmed observations regarding helical gas flow for a multiphase system. The experimental results were used to obtain the coefficients for the radial forces.
A comparison with experimental findings showed that the present two-dimensional CFD model was able to simulate the inhomogeneity of cross-sectional liquid distribution with agreement to measured values over a broad range of flow conditions.
REFERENCES
ACHER, T., KNAUP, M., GÖLL, S., ZANDER, H.J., LENZ, S. and KERBER, C. (2016). “A two-dimensional CFD model for gas-liquid flows in a coil-wound heat exchanger”.Proceedings of the 9th International Conference on Multiphase Flow.
BASSFELD, J. (2017). Untersuchung der Hydrodynamik in gewickelten Wärmetauschern. Master’s thesis, Technische Hochschule Nürnberg Georg Simon Ohm.
DREW, D. and PASSMAN, S. (1999). Theory of multicomponent fluids. Springer.
GÖLL, S., ZANDER, H.J., LENZ, S. and KERBER, C. (2013). “Entwicklung eines 2D-CFD-Modells zur Berechnung der Wärmeübertragung in spiralgewickelten Wärmeaustauschern”. Jahrestreffen der ProcessNet- Fachgruppen CFD, Mischvorgänge und Rheologie Würzburg.
ISHII, M. and HIBIKI, T. (2010). Thermo-fluid dynamics of two-phase flow. Springer.
ISHII, M. and MISHIMA, K. (1989). “Droplet entrainment correlation in annular two-phase flow”. Int. J.
Heat Mass Transf.,32, 1835–1846.
JAKOBSEN, H. (2008). Chemical reactor modeling:
multiphase reactive flows. Springer.
KIEWAT, M. (2015). CFD SIMULATION OF COIL WOUND HEAT EXCHANGERS. Master’s thesis, Friedrich- Alexander-Universität Erlangen-Nürnberg.
KNAUP, M. (2015). Modellierung und numerische Simulation der Zweiphasenströmung in einem spiralgewickelten Wärmetauscher. Master’s thesis, Technische Universität München.
PACIO, J.C. and DORAO, C.A. (2011). “A review on heat exchanger thermal hydraulic models for cryogenic applications”.Cryogenics,51, 366–379.
ROGERS, J.T. (1981). “Laminar Falling Film Flow and Heat Transfer Characteristics on Horizontal Tubes”. Can. J.
Chem. Eng.,59, 213–222.
SCHILLER, L. and NAUMANN, K. (1935).
“Über die grundlegenden Berechnungen bei der Schwerkraftaufbereitung”. Zeitschrift des Vereins der deutschen Ingenieure,77, 316–320.
SIRIGNANO, W. (2010). Fluid Dynamics and Transport of Droplets and Sprays. Cambridge University Press.
STEINBAUER, M. and HECHT, T. (1996). “Optimised Calculation of Helical-Coiled Heat Exchangers in LNG Plants”.Eurogas 96 Conference.
WALTER, T., KERBER, C., BRAUN, K., RICHARDT, C. and STEINBAUER, M. (2014). “Coil wound heat exchangers for LNG - investigation of transport phenomena within the bundle”. Proceedings of the 3rd Trondheim Gas Technology Conference.
447
T. Acher, M. Knaup, K. Braun, H. Zander WANG, T., DING, G., DUAN, Z., REN, T., CHEN, J.
and PU, H. (2015). “A distributed-parameter model for LNG spiral wound heat exchanger based on graph theory”.Applied Thermal Engineering,81, 102 – 113.
448