Shell-side transfer in shell-and-tube heat exchangers

Local mass transfer coefficients were determined by using the electrochemical technique. A simple model of a heat exchanger with segmental nickel tube joined to p.v.c. rods replaced the exchanger tubes. Measurements were made for both no-Ieakage, semi-leakage and total leakage configurations. Baffle-spacings of 47.6 mm, 66.6 mm, 97 mm and 149.2 mm wer studied. Also studied were the overall exchanger pressure drops for each configuration. The comparison of the heat transfer data with this work showed good agreement at high flow rates for the no-leakage case, but the agreement became poor for lower flow rates and leakage configurations. This disagreement was explained by non-analogous driving forces existing in the two systems. The no-leakage data showed length-wise variation of transfer coefficients along the exchanger length. The end compartments showing transfer coefficients inferior by up to 26% compared to tbe internal compartments, depending on Reynolds number. With the introduction of leakage streams this variation however became smaller than the experimental accuracy. A model is outlined to show the characteristic behaviour of individual electrode segments within the compartment. This was able to discriminate between cross and window zones for the no- leakage case, but no such distinction could be made for the leakage case. A flow area was found which, when incorporated in the Reynolds number, enabled the correlation of baffle-cut and baffle-spacing parameters for the no-leakage case . This area is the free flow area determined at the baffle edge. Addition of the leakage area to this flow area resulted in correlation of all commercial leakage geometrical parameters. The procedures used to correlate the pressure drop data from a total of eighteen different configurations on a single curve are also outlined.

Computational fluid dynamical studies of structured distillation packings

In recent years structured packings have become more widely used in the process industries because of their improved volumetric efficiency. Most structured packings consist of corrugated sheets placed in the vertical plane The corrugations provide a regular network of channels for vapour liquid contact. Until recently it has been necessary to develop new packings by trial and error, testing new shapes in the laboratory. The orderly repetitive nature of the channel network produced by a structured packing suggests it may be possible to develop improved structured packings by the application of computational fluid dynamics (CFD) to calculate the packing performance and evaluate changes in shape so as to reduce the need for laboratory testing. In this work the CFD package PHOENICS has been used to predict the flow patterns produced in the vapour phase as it passes through the channel network. A particular novelty of the approach is to set up a method of solving the Navier Stokes equations for any particular intersection of channels. The flow pattern of the streams leaving the intersection is then made the input to the downstream intersection. In this way the flow pattern within a section of packing can be calculated. The resulting heat or mass transfer performance can be calculated by other standard CFD procedures. The CFD predictions revealed a circulation developing within the channels which produce a loss in mass transfer efficiency The calculations explained and predicted a change in mass transfer efficiency with depth of the sheets. This effect was also shown experimentally. New shapes of packing were proposed to remove the circulation and these were evaluated using CFD. A new shape was chosen and manufactured. This was tested experimentally and found to have a higher mass transfer efficiency than the standard packing.

BTX extraction by sulfolane in a mixer-settler cascade

The literature relating to the extraction of the aromatics, benzene, toluene and xylene (BTX) using different commercial solvents, and to mixer-settler design and performance, has been reviewed. Liquid-liquid equilibria of the ternary systems: hexane-benzene-sulfolane, n-heptane-toluene-sulfolane, and octane-xylene-sulfolane were determined experimentally at temperatures of 30oC, 35oC, and 40oC. The work was then extended to a multicomponent system. The data were correlated by using Hand's method and were found to be in a good agreement with theoretical predictions using the UNIFAC method. A study was made of the performance of a 10-stage laboratory mixer-settler cascade for the extraction of BTX from a synthetic reformate utilizing sulfolane as a solvent. Murphree stage efficiency decreased with stage number but 99% extraction was achievable within 4 stages. The effects of temperature, phase ratio, and agitator speed were investigated. The efficiency increased with agitator speed but > 1050 rpm resulted in secondary haze formation. An optimum temperature of 30oC was selected from the phase equilibria; the optimum solvent: feed ratio was 3:1 for 4 stages. The experimental overall mass transfer coefficients were compared with those predicted from single drop correlations and were in all cases greater, by a factor of 1.5 to 3, due to the surface renewal associated with drop break-up and coalescence promoted by agitation. A similar investigation was performed using real reformate from the Kuwait Oil Company. The phase ratios were in the range 0.5 to 1 to 3.25 to 1, the agitator speed 1050 rpm, and the operating temperature 30oC. A maximum recovery of 99% aromatics was achieved in 4 stages at a phase ratio of 3.25 to 1. A backflow model was extended to simulate conditions in the mixer-settler cascade with this multicomponent system. Overall mass transfer coefficients were estimated by obtaining the best fit between experimental and predicted concentration profiles. They were up to 10% greater than those with the synthetic feed but close agreement was not possible because the distribution coefficient and phase ratio varied with stage number. Sulfolane was demonstrated to be an excellent solvent for BTX recovery and a mixer-settler cascade was concluded to be a technically viable alternative to agitated columns for this process.

Transfer coefficients over inclined tube banks

The thesis presents experimental results for shell-side transfer coefficients and pressure drops across four different tube banks, using small-scale models, with yawed tubes, as found in many types of heat exchangers, boilers and nuclear reactors. The tube banks investigated have a staggered tube layout on a rotated square pitch, with a 1.25 pitch-to-diameter ratio. The angle of attack was varied between 45o and 90o. An extensive range of Reynolds number, i.e. 0.5. to 12,600, covering so-called laminar, transition and turbulent flows, was investigated. A diffusion-controlled electrochemical mass transfer technique has been employed to measure mass transfer coefficients. The heat transfer coefficients may be then readily obtained from the mass transfer values by applying the well-established Chilton-Colburn analogy. The results for the normal tube bank, which forms the base case for the study on inclined tube banks, show close agreement with previous work. The transfer coefficients and pressure drops of the inclined tube banks are compared with results from the ideal normal tube bank to examine the effect of inclination angle on heat transfer and pressure drop variations. The variation of the transfer coefficients row-by-row and the entrance and exit effects have also been investigated. An auxilary investigation has been carried out on the role of natural convection. A preliminary correlation of transfer coefficients and pressure drops against the variation in the yaw angle has been attempted. The results are discussed in the light of the few existing theoretical treatments and experimental data for these situations, and recommendations made for future work.

Flow patterns, performance and scale-up of distillation trays

This work is concerned with the nature of liquid flow across industrial sieve trays operating in the spray, mixed, and the emulsified flow regimes. In order to overcome the practical difficulties of removing many samples from a commercial tray, the mass transfer process was investigated in an air water simulator column by heat transfer analogy. The temperature of the warm water was measured by many thermocouples as the water flowed across the single pass 1.2 m diameter sieve tray. The thermocouples were linked to a mini computer for the storage of the data. The temperature data were then transferred to a main frame computer to generate temperature profiles - analogous to concentration profiles. A comprehensive study of the existing tray efficiency models was carried out using computerised numerical solutions. The calculated results were compared with experimental results published by the Fractionation Research Incorporation (FRl) and the existing models did not show any agreement with the experimental results. Only the Porter and Lockett model showed a reasonable agreement with the experimental results for cenain tray efficiency values. A rectangular active section tray was constructed and tested to establish the channelling effect and the result of its effect on circular tray designs. The developed flow patterns showed predominantly flat profiles and some indication of significant liquid flow through the central region of the tray. This comfirms that the rectangular tray configuration might not be a satisfactory solution for liquid maldistribution on sieve trays. For a typical industrial tray the flow of liquid as it crosses the tray from the inlet to the outlet weir could be affected by the mixing of liquid by the eddy, momentum and the weir shape in the axial or the transverse direction or both. Conventional U-shape profiles were developed when the operating conditions were such that the froth dispersion was in the mixed regime, with good liquid temperature distribution while in the spray regime. For the 12.5 mm hole diameter tray the constant temperature profiles were found to be in the axial direction while in the spray regime and in the transverse direction for the 4.5 mm hole tray. It was observed that the extent of the liquid stagnant zones at the sides of the tray depended on the tray hole diameter and was larger for the 4.5 mm hole tray. The liquid hold-up results show a high liquid hold-up at the areas of the tray with low liquid temperatures, this supports the doubts about the assumptions of constant point efficiency across an operating tray. Liquid flow over the outlet weir showed more liquid flow at the centre of the tray at high liquid loading with low liquid flow at both ends of the weir. The calculated results of the point and tray efficiency model showed a general increase in the calculated point and tray efficiencies with an increase in the weir loading, as the flow regime changed from the spray to the mixed regime the point and the tray efficiencies increased from approximately 30 to 80%.Through the mixed flow regime the efficiencies were found to remain fairly constant, and as the operating conditions were changed to maintain an emulsified flow regime there was a decrease in the resulting efficiencies. The results of the estimated coefficient of mixing for the small and large hole diameter trays show that the extent of liquid mixing on an operating tray generally increased with increasing capacity factor, but decreased with increasing weir loads. This demonstrates that above certain weir loads, the effect of eddy diffusion mechanism on the process of liquid mixing on an operating tray to be negligible.

Osmotic dehydration in plant tissues

The primary aim of the thesis is to provide a comprehensive investigation of the osmotic dehydration processes in plant tissue. Effort has been concentrated on the modelling for simulating the processes. Two mathematical models for simulating the mass transfer during osmotic dehydration processes in plant tissues are developed and verified using existing experimental data. Both models are based on the mechanism of diffusion and convection of any mobile material that can transport in plant tissues. The mass balance equation for the transport of each constituent is established separately for intracellular and extra-cellular volumes with taking into account the mass transfer across the cell membrane the intracellular and extra-cellular volumes and the shrinkage of the whole tissue. The contribution from turgor pressure is considered in both models. Model two uses Darcy’s law to build the relation between shrinkage velocity and hydrostatic pressure in each volume because the plant tissue can be considered as the porous medium. Moreover, it has been extended to solve the multi-dimensional problems. A lot of efforts have been made to the parameter study and the sensitivity analyses. The parameters investigated including the concentration of the osmotic solution, diffusion coefficient, permeability of the cell membrane, elastic modulus of the cell wall, critical cell volume etc. The models allow us to quantitatively simulate the time evolution of intracellular and extra-cellular volumes as well as the time evolution of concentrations in each cross-section.

New fuel cells system for portable application with low temperature cofired ceramic (LTCC) technology

Miniature direct methanol fuel cells (DMFCs) are promising micro power sources for portable appliction. Low temperature cofired ceramic (LTCC), a competitive technology for current MEMS based fabrication, provides cost-effective mass manufacturing route for miniature DMFCs. Porous silver tape is adapted as electrodes to replace the traditional porous carbon electrodes due to its compatibility to LTCC processing and other electrochemical advantages. Electrochemical evaluation of silver under DMFCs operating conditions demonstrated that silver is a good electrode for DMFCs because of its reasonable corrosion resistance, low passivating current, and enhanced catalytic effect. Two catalyst loading methods (cofiring and postfiring) of the platinum and ruthenium catalysts are evaluated for LTCC based processing. The electrochemical analysis exhibits that the cofired path out-performs the postfiring path both at the anode and cathode. The reason is the formation of high surface area precipitated whiskers. Self-constraint sintering is utilized to overcome the difficulties of the large difference of coefficient of thermal expansion (CTE) between silver and LTCC (Dupont 951) tape during cofiring. The graphite sheet employed as a cavity fugitive insert guarantees cavity dimension conservation. Finally, performance of the membrane electrode assembly (MEA) with the porous silver electrode in the regular graphite electrode based cell and the integrated cofired cell is measured under passive fuel feeding condition. The MEA of the regular cell performs better as the electrode porosity and temperature increased. The power density of 10 mWcm-2 was obtained at ambient conditions with 1M methanol and it increased to 16 mWcm -2 at 50°C from an open circuit voltage of 0.58V. For the integrated prototype cell, the best performance, which depends on the balance methanol crossover and mass transfer at different temperatures and methanol concentrations, reaches 1.13 mWcm-2 at 2M methanol solution at ambient pressure. The porous media pore structure increases the methanol crossover resistance. As temperature increased to 60°C, the device increases to 2.14 mWcm-2.

(Table T1) Density, porosity, void ratio, diffusivity and permeability of ODP Site 185-1149 sediments

The subduction of oceanic plates regulates crustal growth, influences arc volcanism, and refertilizes the mantle. Continental growth occurs by subduction of crustal material (seawater components, marine sediments, and basaltic crust). The geochemical and physical evolution of the Earth's crust depends, in large part, on the fate of subducted material at convergent margins (Armstrong, 1968, doi:10.1029/RG006i002p00175; Karig and Kay, 1981, 10.1098/rsta.1981.0108). The crustal material on the downgoing plate is recycled to various levels in the subduction zone. The recycling process that takes place in the "Subduction Factory" is difficult to observe directly but is clearly illuminated using chemical tracers. Von Huene and Scholl (1991, doi:10.1029/91RG00969) and Plank and Langmuir (1993, doi:10.1038/362739a0) preliminarily calculated a large flux of subducted materials. By mass balancing the chemical tracers and measuring the fractionations that occur between them, the Subduction Factory work and the effect on the Earth's evolution can be estimated. In order to elucidate this mass balance, Ocean Drilling Program Leg 185 drilled two deepwater shales into the oceanic crust situated in the Mariana-Izu Trenches and recovered core samples of incoming oceanic crust. The calculations of mass circulation in the subduction zone, however, did not take into account the mass transfer properties within subducted oceanic crust, although the dewatering fluid and diffused ions may play an important role in various activities such as seismogeneity, serpentine diapiring, and arc volcanism. Thus, this paper focuses on the quantitative measurements of the physical and mass transfer properties of subducted oceanic crust.

Laboratory Investigation of Diluted Bitumen Trapping and Dissolution in Gravel

The Canadian economy is largely dependent on the distribution of large volumes of oil to domestic and international markets by a long network of pipelines. Unfortunately, accidents occur, and oil can leak or spill from these pipelines before it reaches its destination. Of particular concern are the long-term consequences of oil spills in freshwater, which include sinking of oil in water and the contamination of sensitive areas, such as where fish (e.g., salmon) deposit their eggs in gravel-dominated river sediments. There is a knowledge gap regarding the fate and behaviour of oil in river sediment. To this end, this study aimed at finding the potential for diluted bitumen (dilbit) oil to become trapped in gravel and to transfer hydrocarbons into water by dissolution, which are harmful to aquatic life. Two sets of laboratory experiments were conducted to simulate conditions of an oil spill on an exposed bank or in shallow water. In the first set, by conducting capillary pressure-saturation (Pc-Sw) experiments it was found that dilbit can enter gravel pores without much resistance and approximately 14% of the pore volume can be occupied by discontinuous single or multipore blobs of dilbit following imbibition of water. Air-water Pc-Sw experiments done in laboratory 1-D columns required gravity correction and could be successfully scaled to predict dilbit-water Pc-Sw curves, except for the trapped amount of dilbit. Trapped dilbit constituents can be dissolved into the water flowing through gravel pores (hyporheic flow) at different velocities. In the second set, dissolution experiments suggested that for the duration of the test, hydrocarbons that cause acute toxicity dissolve rapidly, likely resulting in a decrease in their effective solubility. However, dilbit saturation changed only <2% within that time. Chronically toxic PAH compounds were also detected in the effluent water. The total concentration of all detected PAHs and alkylPAHs exceeded the threshold literature value to protect early-life stage fish. Observations of decreased concentrations with increased aqueous velocities as well as less than equilibrium concentrations indicated that the mass transfer was rate-limited. A correlation was developed for the mass transfer rate coefficient to understand the mass transfer behaviour beyond the conditions used in the experiments, which had a Reynolds number exponent similar to the studies of NAPL dissolution in groundwater.

Modeling the impact of large-scale solid mass transfers on the gravity field

The focus of this thesis is to explore and quantify the response of large-scale solid mass transfer events on satellite-based gravity observations. The gravity signature of large-scale solid mass transfers has not been deeply explored yet; mainly due to the lack of significant events during dedicated satellite gravity missions‘ lifespans. In light of the next generation of gravity missions, the feasibility of employing satellite gravity observations to detect submarine and surface mass transfers is of importance for geoscience (improves the understanding of geodynamic processes) and for geodesy (improves the understanding of the dynamic gravity field). The aim of this thesis is twofold and focuses on assessing the feasibility of using satellite gravity observations for detecting large-scale solid mass transfers and on modeling the impact on the gravity field caused by these events. A methodology that employs 3D forward modeling simulations and 2D wavelet multiresolution analysis is suggested to estimate the impact of solid mass transfers on satellite gravity observations. The gravity signature of various submarine and subaerial events that occurred in the past was estimated. Case studies were conducted to assess the sensitivity and resolvability required in order to observe gravity differences caused by solid mass transfers. Simulation studies were also employed in order to assess the expected contribution of the Next Generation of Gravity Missions for this application.

Constant wall heat flux boundary conditions in micro-channels filled with porous material with internal heat generation under local thermal non-equilibrium conditions

Forced convection heat transfer in a micro-channel filled with a porous material saturated with rarefied gas with internal heat generation is studied analytically in this work. The study is performed by analysing the boundary conditions for constant wall heat flux under local thermal non-equilibrium (LTNE) conditions. Invoking the velocity slip and temperature jump, the thermal behaviour of the porous-fluid system is studied by considering thermally and hydrodynamically fully-developed conditions. The flow inside the porous material is modelled by the Darcy–Brinkman equation. Exact solutions are obtained for both the fluid and solid temperature distributions for two primary approaches models A and B using constant wall heat flux boundary conditions. The temperature distributions and Nusselt numbers for models A and B are compared, and the limiting cases resulting in the convergence or divergence of the two models are also discussed. The effects of pertinent parameters such as fluid to solid effective thermal conductivity ratio, Biot number, Darcy number, velocity slip and temperature jump coefficients, and fluid and solid internal heat generations are also discussed. The results indicate that the Nusselt number decreases with the increase of thermal conductivity ratio for both models. This contrasts results from previous studies which for model A reported that the Nusselt number increases with the increase of thermal conductivity ratio. The Biot number and thermal conductivity ratio are found to have substantial effects on the role of temperature jump coefficient in controlling the Nusselt number for models A and B. The Nusselt numbers calculated using model A change drastically with the variation of solid internal heat generation. In contrast, the Nusselt numbers obtained for model B show a weak dependency on the variation of internal heat generation. The velocity slip coefficient has no noticeable effect on the Nusselt numbers for both models. The difference between the Nusselt numbers calculated using the two models decreases with an increase of the temperature jump coefficient.

DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

Computational modelling of non-equilibrium condensing steam flows in low-pressure steam turbines

The steam turbines play a significant role in global power generation. Especially, research on low pressure (LP) steam turbine stages is of special importance for steam turbine man- ufactures, vendors, power plant owners and the scientific community due to their lower efficiency than the high pressure steam turbine stages. Because of condensation, the last stages of LP turbine experience irreversible thermodynamic losses, aerodynamic losses and erosion in turbine blades. Additionally, an LP steam turbine requires maintenance due to moisture generation, and therefore, it is also affecting on the turbine reliability. Therefore, the design of energy efficient LP steam turbines requires a comprehensive analysis of condensation phenomena and corresponding losses occurring in the steam tur- bine either by experiments or with numerical simulations. The aim of the present work is to apply computational fluid dynamics (CFD) to enhance the existing knowledge and understanding of condensing steam flows and loss mechanisms that occur due to the irre- versible heat and mass transfer during the condensation process in an LP steam turbine. Throughout this work, two commercial CFD codes were used to model non-equilibrium condensing steam flows. The Eulerian-Eulerian approach was utilised in which the mix- ture of vapour and liquid phases was solved by Reynolds-averaged Navier-Stokes equa- tions. The nucleation process was modelled with the classical nucleation theory, and two different droplet growth models were used to predict the droplet growth rate. The flow turbulence was solved by employing the standard k-ε and the shear stress transport k-ω turbulence models. Further, both models were modified and implemented in the CFD codes. The thermodynamic properties of vapour and liquid phases were evaluated with real gas models. In this thesis, various topics, namely the influence of real gas properties, turbulence mod- elling, unsteadiness and the blade trailing edge shape on wet-steam flows, are studied with different convergent-divergent nozzles, turbine stator cascade and 3D turbine stator-rotor stage. The simulated results of this study were evaluated and discussed together with the available experimental data in the literature. The grid independence study revealed that an adequate grid size is required to capture correct trends of condensation phenomena in LP turbine flows. The study shows that accurate real gas properties are important for the precise modelling of non-equilibrium condensing steam flows. The turbulence modelling revealed that the flow expansion and subsequently the rate of formation of liquid droplet nuclei and its growth process were affected by the turbulence modelling. The losses were rather sensitive to turbulence modelling as well. Based on the presented results, it could be observed that the correct computational prediction of wet-steam flows in the LP turbine requires the turbulence to be modelled accurately. The trailing edge shape of the LP turbine blades influenced the liquid droplet formulation, distribution and sizes, and loss generation. The study shows that the semicircular trailing edge shape predicted the smallest droplet sizes. The square trailing edge shape estimated greater losses. The analysis of steady and unsteady calculations of wet-steam flow exhibited that in unsteady simulations, the interaction of wakes in the rotor blade row affected the flow field. The flow unsteadiness influenced the nucleation and droplet growth processes due to the fluctuation in the Wilson point.

FLOW REGIME DRIVEN THERMAL ENHANCEMENT IN INTERNALLY-GROOVED TUBES

Internally-grooved refrigeration tubes maximize tube-side evaporative heat transfer rates and have been identified as a most promising technology for integration into compact cold plates. Unfortunately, the absence of phenomenological insights and physical models hinders the extrapolation of grooved-tube performance to new applications. The success of regime-based heat transfer correlations for smooth tubes has motivated the current effort to explore the relationship between flow regimes and enhanced heat transfer in internally-grooved tubes. In this thesis, a detailed analysis of smooth and internally-grooved tube data reveals that performance improvement in internally-grooved tubes at low-to-intermediate mass flux is a result of early flow regime transition. Based on this analysis, a new flow regime map and corresponding heat transfer coefficient correlation, which account for the increased wetted angle, turbulence, and Gregorig effects unique to internally-grooved tubes, were developed. A two-phase test facility was designed and fabricated to validate the newly-developed flow regime map and regime-based heat transfer coefficient correlation. As part of this setup, a non-intrusive optical technique was developed to study the dynamic nature of two-phase flows. It was found that different flow regimes result in unique temporally varying film thickness profiles. Using these profiles, quantitative flow regime identification measures were developed, including the ability to explain and quantify the more subtle transitions that exist between dominant flow regimes. Flow regime data, based on the newly-developed method, and heat transfer coefficient data, using infrared thermography, were collected for two-phase HFE-7100 flow in horizontal 2.62mm - 8.84mm diameter smooth and internally-grooved tubes with mass fluxes from 25-300 kg/m²s, heat fluxes from 4-56 kW/m², and vapor qualities approaching 1. In total, over 6500 combined data points for the adiabatic and diabatic smooth and internally-grooved tubes were acquired. Based on results from the experiments and a reinterpretation of data from independent researchers, it was established that heat transfer enhancement in internally-grooved tubes at low-to-intermediate mass flux is primarily due to early flow regime transition to Annular flow. The regime-based heat transfer coefficient outperformed empirical correlations from the literature, with mean and absolute deviations of 4.0% and 32% for the full range of data collected.

AIR SIDE HEAT TRANSFER ENHANCEMENT IN HEAT EXCHANGERS UTILIZING INNOVATIVE DESIGNS AND THE ADDITIVE MANUFACTURING TECHNIQUE

Over the last decade, rapid development of additive manufacturing techniques has allowed the fabrication of innovative and complex designs. One field that can benefit from such technology is heat exchanger fabrication, as heat exchanger design has become more and more complex due to the demand for higher performance particularly on the air side of the heat exchanger. By employing the additive manufacturing, a heat exchanger design was successfully realized, which otherwise would have been very difficult to fabricate using conventional fabrication technologies. In this dissertation, additive manufacturing technique was implemented to fabricate an advanced design which focused on a combination of heat transfer surface and fluid distribution system. Although the application selected in this dissertation is focused on power plant dry cooling applications, the results of this study can directly and indirectly benefit other sectors as well, as the air-side is often the limiting side for in liquid or single phase cooling applications. Two heat exchanger designs were studied. One was an advanced metallic heat exchanger based on manifold-microchannel technology and the other was a polymer heat exchanger based on utilization of prime surface technology. Polymer heat exchangers offer several advantages over metals such as antifouling, anticorrosion, lightweight and often less expensive than comparable metallic heat exchangers. A numerical modeling and optimization were performed to calculate a design that yield an optimum performance. The optimization results show that significant performance enhancement is noted compared to the conventional heat exchangers like wavy fins and plain plate fins. Thereafter, both heat exchangers were scaled down and fabricated using additive manufacturing and experimentally tested. The manifold-micro channel design demonstrated that despite some fabrication inaccuracies, compared to a conventional wavy-fin surface, 15% - 50% increase in heat transfer coefficient was possible for the same pressure drop value. In addition, if the fabrication inaccuracy can be eliminated, an even larger performance enhancement is predicted. Since metal based additive manufacturing is still in the developmental stage, it is anticipated that with further refinement of the manufacturing process in future designs, the fabrication accuracy can be improved. For the polymer heat exchanger, by fabricating a very thin wall heat exchanger (150μm), the wall thermal resistance, which usually becomes the limiting side for polymer heat exchanger, was calculated to account for only up to 3% of the total thermal resistance. A comparison of air-side heat transfer coefficient of the polymer heat exchanger with some of the commercially available plain plate fin surface heat exchangers show that polymer heat exchanger performance is equal or superior to plain plate fin surfaces. This shows the promising potential for polymer heat exchangers to compete with conventional metallic heat exchangers when an additive manufacturing-enabled fabrication is utilized. Major contributions of this study are as follows: (1) For the first time demonstrated the potential of additive manufacturing in metal printing of heat exchangers that benefit from a sophisticated design to yield a performance substantially above the respective conventional systems. Such heat exchangers cannot be fabricated with the conventional fabrication techniques. (2) For the first time demonstrated the potential of additive manufacturing to produce polymer heat exchangers that by design minimize the role of thermal conductivity and deliver a thermal performance equal or better that their respective metallic heat exchangers. In addition of other advantages of polymer over metal like antifouling, anticorrosion, and lightweight. Details of the work are documented in respective chapters of this thesis.