Multiphase flow provides a solution to the high heat flux and precision required by modern-day gadgets and heat transfer devices as phase change processes make high heat transfer rates achievable at moderate temperature differences. An application of multiphase flow commonly used in industry is the condensation of refrigerants in inclined tubes. The identification of two-phase flow patterns, or flow regimes, is fundamental to the successful design and subsequent optimisation given that the heat transfer efficiency and pressure gradient are dependent on the flow structure of the working fluid. This study showed that with visualisation data and artificial neural networks (ANN), a machine could learn, and subsequently classify the separate flow patterns of condensation of R-134a refrigerant in inclined smooth tubes with more than 98% accuracy. The study considered 10 classes of flow pattern images acquired from previous experimental works that cover a wide range of flow conditions and the full range of tube inclination angles. Two types of classifiers were considered, namely multilayer perceptron (MLP) and convolutional neural networks (CNN). Although not the focus of this study, the use of a principal component analysis (PCA) allowed feature dimensionality reduction, dataset visualisation, and decreased associated computational cost when used together with multilayer perceptron neural networks. The superior two-dimensional spatial learning capability of convolutional neural networks allowed improved image classification and generalisation performance across all 10 flow pattern classes. In both cases, the classification was done sufficiently fast to enable real-time implementation in two-phase flow systems. The analysis sequence led to the development of a predictive tool for the classification of multiphase flow patterns in inclined tubes, with the goal that the features learnt through visualisation would apply to a broad range of flow conditions, fluids, tube geometries and orientations, and would even generalise well to identify adiabatic and boiling two-phase flow patterns. The method was validated by the prediction of flow pattern images found in the existing literature. / Dissertation (MEng)--University of Pretoria, 2021. / NRF / Mechanical and Aeronautical Engineering / MEng / Restricted

The rapid depletion of hydrocarbon fields around the world has led the industry to search for these resources in ever increasing water depths. In this context, the large diameter (D > 100mm) vertical riser has become a subject of great interest. In this research work, a major investigation was undertaken to determine the two phase flow hydrodynamics in a 254mm vertical riser. Two types of experiments were performed for range of air-water superficial velocities. The first experimental campaign addresses the issue of the two gas injector’s performances (conventional vs. novel design gas injector) in the large diameter vertical riser. The experimental results show that the novel design gas injector should be the preferential choice. The second set of the experimental work investigates the two phase flow hydrodynamics in the vertical riser in detail. The two phase flow patterns and their transitions were identified by combination of visual observations and statistical features. Based on the results, the experimental flow regime map was developed and compared with the existing vertical upflow regime maps/models. None of the flow regime transition models adequately predicted the flow regimes transitions in large diameter vertical risers as a whole. In this regard, the Taitel et al. (1980) bubble to slug flow transition model has been modified for large diameter vertical upflow conditions, based on the physical mechanism observed. The general trends of modified criteria agreed well with the current and other large diameter experimental results. The effect of upstream conditions on the vertical riser flow behaviour was also investigated in detail by two different inlet configurations (i) near riser base injection and (ii) upstream flowline injection. It was found that no significant differences exist in flow behaviour at low air-water superficial velocities for both the inlet configuration, at high air-water superficial velocities, the intermittent flow behavior in flowline influences the riser flow pattern characteristics and thereby controls the riser dynamics. It is found that liquid slugs from the flowline naturally dissipate to some extent in the riser as a consequence of compression of succeeding bubble that rapidly expands and break through the liquid slug preceding it when it enters the riser. The experimental work corroborates the general consensus that slug flow does not exist in large diameter vertical upflow condition. Experimental data has been further compared to increase the confidence on the existing two phase flow knowledge on large diameter vertical riser: (a) by comparing with other experimental studies on large diameter vertical upflow in which generally, a good agreement was found, (b) by assessing the predictive capability of void fraction correlations/pressure gradient methods. The important implication of this assessment is that the mechanistic approach based on specific flow regime in determining the void fraction and pressure gradient is more successful than conventional empirical based approaches. The assessment also proposes a proposed set a of flow regime specific correlations that recommends void fraction correlations based on their performances in the individual flow regimes. Finally, a numerical model to study the hydrodynamic behaviour in the large diameter horizontal flowline-vertical riser system is developed using multiphase flow simulator OLGA. The simulated results show satisfactory agreement for the stable flows while discrepancies were noted for highly intermittent flows. The real time boundary application was partially successful in qualitatively reproducing the trends. The discrepancies between the predicted results and experimental data are likely to be related to the incorrect closure relations used based on incorrect flow regimes predictions. The existence of the multiple roots in the OLGA code is also reported for the first time.

The rapid depletion of hydrocarbon fields around the world has led the industry to search for these resources in ever increasing water depths. In this context, the large diameter (D > 100mm) vertical riser has become a subject of great interest. In this research work, a major investigation was undertaken to determine the two phase flow hydrodynamics in a 254mm vertical riser. Two types of experiments were performed for range of air-water superficial velocities. The first experimental campaign addresses the issue of the two gas injector’s performances (conventional vs. novel design gas injector) in the large diameter vertical riser. The experimental results show that the novel design gas injector should be the preferential choice. The second set of the experimental work investigates the two phase flow hydrodynamics in the vertical riser in detail. The two phase flow patterns and their transitions were identified by combination of visual observations and statistical features. Based on the results, the experimental flow regime map was developed and compared with the existing vertical upflow regime maps/models. None of the flow regime transition models adequately predicted the flow regimes transitions in large diameter vertical risers as a whole. In this regard, the Taitel et al. (1980) bubble to slug flow transition model has been modified for large diameter vertical upflow conditions, based on the physical mechanism observed. The general trends of modified criteria agreed well with the current and other large diameter experimental results. The effect of upstream conditions on the vertical riser flow behaviour was also investigated in detail by two different inlet configurations (i) near riser base injection and (ii) upstream flowline injection. It was found that no significant differences exist in flow behaviour at low air-water superficial velocities for both the inlet configuration, at high air-water superficial velocities, the intermittent flow behavior in flowline influences the riser flow pattern characteristics and thereby controls the riser dynamics. It is found that liquid slugs from the flowline naturally dissipate to some extent in the riser as a consequence of compression of succeeding bubble that rapidly expands and break through the liquid slug preceding it when it enters the riser. The experimental work corroborates the general consensus that slug flow does not exist in large diameter vertical upflow condition. Experimental data has been further compared to increase the confidence on the existing two phase flow knowledge on large diameter vertical riser: (a) by comparing with other experimental studies on large diameter vertical upflow in which generally, a good agreement was found, (b) by assessing the predictive capability of void fraction correlations/pressure gradient methods. The important implication of this assessment is that the mechanistic approach based on specific flow regime in determining the void fraction and pressure gradient is more successful than conventional empirical based approaches. The assessment also proposes a proposed set a of flow regime specific correlations that recommends void fraction correlations based on their performances in the individual flow regimes. Finally, a numerical model to study the hydrodynamic behaviour in the large diameter horizontal flowline-vertical riser system is developed using multiphase flow simulator OLGA. The simulated results show satisfactory agreement for the stable flows while discrepancies were noted for highly intermittent flows. The real time boundary application was partially successful in qualitatively reproducing the trends. The discrepancies between the predicted results and experimental data are likely to be related to the incorrect closure relations used based on incorrect flow regimes predictions. The existence of the multiple roots in the OLGA code is also reported for the first time.

The use of a swirling flow can provide a more uniform velocity distribution and a calmer filling condition according to previous studies of both ingot and continuous casting processes of steel. However, the existing swirling flow generation methods developed in last decades all have some limitations. Recently, a new swirling flow generator, the TurboSwirl device, was proposed. In this work, the convergent nozzle was studied with different angles. The maximum wall shear stress can be reduced by changing the convergent angle between 40º and 60º to obtain a higher swirl intensity. Also, a lower maximum axial velocity can be obtained with a smaller convergent angle. Furthermore, the maximum axial velocity and wall shear stress can also be affected by moving the location of the vertical runner. A water model experiment was carried out to verify the simulation results of the effect of the convergent angle on the swirling flow pattern. The shape of the air-core vortex in the water model experiment could only be accurately simulated by using the Reynolds Stress Model (RSM). The simulation results were also validated by the measured radial velocity in the vertical runner by the ultrasonic velocity profiler (UVP). The TurboSwirl was reversed and connected to a traditional SEN to generate the swirling flow. The periodic characteristic of the swirling flow and asymmetry flow pattern were observed in both the simulated and measured results. The detached eddy simulation (DES) turbulence model was used to catch the time-dependent flow pattern and the predicted results agree well with measured axial and tangential velocities. This new design of the SEN with the reverse TurboSwirl could provide an almost equivalent strength of the swirling flow generated by an electromagnetic swirling flow generator. It can also reduce the downward axial velocities in the center of the SEN outlet and obtain a calmer meniscus and internal flow in the mold. / Tidigare studier visar att ett roterande flöde kan ge en mer likformig hastighetsfördelning och en lugnare fyllning i både göt- och stränggjutning av stål. De befintliga metoderna för att generera ett roterande flöde har vissa begränsningar. En ny metod för att generera det roterande flödet, en så kallad TurboSwirl, föreslogs nyligen. I detta arbete undersöktes ett konvergent munstycke med olika vinklar för att se hur detta påverkade det roterande flödet som genererades i anordningen. Resultaten visar att skjuvspänningen i systemet kan reduceras genom att ändra munstyckets vinkel mellan 40º till 60º. En lägre maximal axiell hastighet kan också uppnås med en mindre konvergent vinkel på munstycket. Det är även möjligt att påverka den maximala axiella hastigheten och skjuvspänningen i systemet genom att förflytta den vertikala kanalen i anordningen. Vattenmodellexperiment har utförts för att validera simuleringsresultaten. Det kraftigt roterande flödet kunde endast beskrivas väl av Reynolds Stress Model (RSM). Validering utfördes också genom att mäta den radiella hastigheten i den vertikala kanalen med en Ultrasonic Velocity Profiler (UVP). TurboSwirl-anordningen vändes och kopplades till gjutröret för att generera det roterande flödet. Detta studerades både med numeriska modeller och med vattenmodellering. Ett periodiskt asymmetriskt roterande flöde observerades både i numeriska modellerna och i vattenmodellerna. För att modellera detta periodiska flöde så användes detached eddy simulation (DES) modellen. Resultaten då denna modell användes stämmer väl med de experimentella mätningarna. Denna nya design med TurboSwirl kan uppnå liknande styrka på det roterande flödet som när elektromagnetisk omrörning användes. Det resulterande roterande flödet leder till en lägre axiell hastighet i gjutröret samt en lugnare yta och ett lugnare flöde i kokillen. / <p>QC 20161123</p>

Em linhas gerais, esta dissertação de mestrado envolve o estudo de padrões de escoaomento durante a ebulição convectiva em microcanais. Resultados experimentais foram levantados para um tubo com diâmetro de 2,32 mm durante a evaporação convectiva dos refrigerantes R134a e R245fa. Para a investigação, técnicas experimentais e de análise foram desenvolvidas. A caracterização dos padrões de escoamento envolveu o tratamento simultâneo de sinais provenientes dos seguintes dispositivos: um par diodo/sensor-laser tendo um tubo transparente entre eles no interior do qual ocorre o escoamento bifásico; um transdutor de pressão piezo-elétrico de tamanho reduzido com o objetivo de determinar a variação local da pressão do escoamento; e de um micro-termopar em contato com o fluido refrigerante. A técnica de tratamento de dados utilizada envolve a aglomeração progressiva de dados que apresentem características médias similares através do algoritmo k-means. Os sinais de pressão, intensidade de radiação e temperatura foram adquiridos simultaneamente a uma freqüência de 25 kHz. Imagens simultâneas do escoamento bifásico a uma velocidade de captura em torno de 10.000 imagens/s foram levantadas através de uma câmera de filmagem rápida (até 100.000 imagens/s), e os padrões de escoamentos observados contrastados aos resultados fornecidos pelo método proposto. Baseado nesta análise, mapas de escoamento foram propostos, os quais incorporaram não apenas critérios subjetivos como a visualização, mas também objetivos como as variações transientes da pressão local do escoamento e da morfologia do escoamento através do seu efeito na dispersão da radiação emitida pelo foto-diodo. Os resultados previstos pelo método objetivo apresentam concordância razoável com os dados caracterizados com base em visualizações. Adicionalmente, características de bolhas alongadas foram determinadas. / The present research has been focused on the study of flow patterns inside channels of diameter less than 3 mm during the convective evaporation of refrigerants such as R134a and R245fa. For the investigation of such topics, experimental techniques and methods of analysis of results were developed. A broad database was gathered in an experimental test facility. The characterization of flow patterns involved the simultaneous processing of signals from the following devices: a pair diode / laser-sensor having a transparent tube between them, within which occurs the two-phase flow; a micro piezoelectric pressure transducer to determine the local variation of pressure of the flow and a micro-thermocouple fixed within the fluid. The technique used in data processing involves the gradual agglomeration of data having similar average characteristics; this method was developed based on the k-means clustering algorithm. The signals from the transducers were acquired simultaneously at a frequency of 25 kHz. The program for the acquisition and for processing of the signals was developed using LabView. Simultaneous images of two-phase flow at a speed of capture around 10,000 images / s were obtained through a high speed camera and the observed flow patterns were contrasted to the results provided by the objective method. Based on this analysis, flow maps were proposed, which incorporate not only subjective criteria such as visualization, but also objective criteria like the transient variations of local pressure of the flow, temperature of the fluid and the effect of the flow morphology based on the dispersion of light which effect was captured by the photo-diode. The maps obtained by the objective method were compared against flow pattern segregated based on visualization and a reasonable agreement was obtained. Besides the elongated bubble characteristics were determined.

Transportation of unprocessed multiphase reservoir fluids from deep/ultra deep offshore through a long subsea tieback/pipeline is inevitable. This form of transportation is complex and requires accurate knowledge of critical transport velocity, flow pattern changes, phase velocity, pressure drop, particle drag & lift forces, sand/liquid/gas holdup, flow rate requirement and tieback sizing etc at the early design phase and during operation for process optimisation. This research investigated sand transport characteristics in multiphase, water‐oil‐gas‐sand flows in horizontal, inclined and vertical pipes. Two critical factors that influence the solid particle transport in the case of multiphase flow in pipes were identified; these are the transient phenomena of flow patterns and the characteristic drag & lift coefficients ( D C , L C ). Therefore, the equations for velocity profile were developed for key flow patterns such as dispersed bubble flow, stratified flow, slug flow and annular flow using a combination of analytical equations and numerical simulation tool (CFD). The existing correlations for D C & L C were modified with data acquired from multiphase experiment in order to account for different flow patterns. Minimum Transport Velocity (MTV) models for suspension and rolling were developed by combining the numerically developed particle velocity profile models with semi‐empirical models for solid particle transport. The models took into account the critical parameters that influence particle transport in pipe flow such as flow patterns and particle drag & lift coefficients, thus eliminate inaccuracies currently experienced with similar models in public domain. The predictions of the proposed MTV models for suspension and rolling in dispersed bubble, slug flow and annular flow show maximum average error margin of 12% when compared with experimental data. The improved models were validated using previously reported experimental data and were shown to have better predictions when compared with existing models in public domain. These models have the potential to solve the problems of pipe and equipment sizing, the risk of sand deposition and bed formation, elimination of costs of sand unloading, downtime and generally improve sand management strategies.

Fluidized beds have been applied in many industrial processes (e.g. coal combustion, gasification and granulation) as an effective means for providing excellent gas and solids contact and mixing, as well as good heat transfer. Although research on the fluidized bed has been carried out for more than 70 years, uncertainties and difficulties still remain. These challenges exist primarily due to the complex and dynamic flow structure within fluidized beds and the lack of reliable measurement techniques. The positron emission particle tracking (PEPT) technique, developed at the University of Birmingham, enables individual particles to be tracked non-invasively in opaque three-dimensional (3-D) fluidized beds and offers favourable temporal and spatial resolutions. PEPT is considered to be a powerful tool for fluidized bed studies and was utilized in the current study to investigate the dynamic behaviour of solid and gas in fluidized beds. The experiments in this study were conducted in a 150-mm inner diameter (I.D.) column and operated in the bubbling fluidization regime at ambient conditions. The effects of various factors on the solid flow structure were examined: solid properties, superficial gas velocity, bed height-to-diameter aspect ratio (H/D) and pore size of the air distributor. The solid flow structure was classified into four patterns, namely patterns A, B, C and D, in which pattern C was newly observed in this thesis. The solid motion, bubble behaviour (i.e., bubble spatial distribution, bubble size and bubble rise velocity) and solid mixing were assessed for each flow pattern to understand their unique fluidization behaviours. This assessment was achieved by the development of three methods: a method to reconstruct bubble behaviours based on solid motion, and two methods for estimating the solid mixing profile in this thesis. The results were discussed and compared with the published literature. The bubble rise velocity and bubble size calculated in this research from the PEPT-measured data was in agreement with other research, particularly that of Kunii and Levenspiel, Yasui and Johanson, and Mori and Wen. Finally, a parameter was developed to predict and control flow patterns based on particle kinetic energy and various factors. The outcomes of this study advance the understanding of the complicated dynamics of bubbling fluidized beds and may benefit several industries in the enhancement of fluidized bed design and control to achieve desirable qualities and efficiencies.

Ultrasound, applied in combination with microbubbles, has potential as a means to enhance the uptake of therapeutic agents, which could include drugs and nucleic acids, into biological cells. This process is commonly referred to as 'sonoporation', and the enhanced uptake can be caused through the incident ultrasonic pressure fi eld causing radial oscillations (cavitation) in the microbubbles, amongst other possibilities. However, the mechanisms responsible for any resultant increase in cell membrane permeability are not yet fully understood. This project focussed on achieving a more fundamental understanding of these salient processes by building on a platform of previous work within the group. One strand of the project involved a complete characterisation of the performance of a rotating mirror high speed camera (Cordin 550-62) that was previously used by our group [and others] to investigate microbubble cavitation phenomena and interactions with proximal cell membranes. Speci cally, I present herein an investigation into the image formation process with this type of camera, the essence of which stymied previous data interpretations. I demonstrate that an inherent asynchrony in the exposure of pixels within individual image frames leads to a temporal anomaly. This was achieved using low cost, flashing LED lights and resulted in the extraction of an algorithm to correct for the temporal anomaly. In a slightly diff erent context, the delivery of suitable ultrasonic fields is necessary to achieve a uniform treatment across a therapeutic target. This thesis also reports on a study on the design of ultrasonic lenses to alter the focal region of a focussed ultrasound transducer with the aim of producing focal regions that can enable sonoporation of tumours of varying sizes. We show that the use of lenses can be an inexpensive alternative to more complex systems such as phased array transducers. Design modelling and experimental testing of lens prototypes are presented along with preliminary results with tissue mimicking polyacrylamide gel phantoms. The target environment in which the process of sonoporation will be clinically useful (i.e. in the physiological circulation) can be simpli ed as a microfluidic system. One strategy for bubble mediated therapy involves the use of a pro-drug approach, that is, when two otherwise benign ingredients are loaded onto separate microbubble populations, but can become mixed at the anatomical target site by the action of focussed ultrasound whereupon a potent drug is produced. The required mixing can be achieved by the violent coalescence of nearby cavitating bubbles, their reaction product then being released and di ffused into the interiour of nearby cells through sonoporation. A study related to this field is presented here where laser induced thermocapillary flows are shown to cause mixing of the content of a drop in a microfluidic channel in a bid to understand the mixing process at a level that may assist future microbubble engineering strategy. To summarise then, the work presented in this thesis has consolidated earlier unpublished data sets achieved by the group, providing new and exacting experimental evidence and an accurate algorithm that will facilitate post-processing of that earlier data (Chapters 2-3). Moreover, group aspirations to translate earlier in-vitro work on sonoporation towards next phase medical-phantom exposures have been boosted through the provision of a new direction involving acoustic lensing, the experimental data from which was used to completely validate existing models for our own design scenarios (Chapter 4). Finally, previous unpublished observations on microbubble coalescence undertaken by the group suggested a means to implement pro-drug delivery with direct in-situ mixing. Such suggestions were explored within microfluidic contexts using lasers to control and visualise the mixing processes that might arise in such situations (Chapter 5). All of these new insights have served to consolidate the group's previous and as yet unpublished data, opening the way for dissemination with confidence in the integrity of that data, and have also extended group capability and expertise in the areas of MHz-rate high speed framing cameras, the fabrication of acoustic lenses, and with microfluidic mixing.

O presente trabalho envolve um estudo teórico e experimental do escoamento bifásico externo a banco de tubos. Inicialmente, apresenta-se uma ampla revisão da literatura sobre padrões de escoamento, fração de vazio e perda de pressão, durante escoamentos monofásicos e bifásicos externos a banco de tubos. Nesta análise são também descritos os métodos de previsão destes parâmetros. Verificam-se diferenças significativas entre as estimativas proporcionadas por eles, fato que indica a inexistência de métodos generalizados. Posteriormente é apresentada uma descrição detalhada da bancada experimental projetada e construída durante o doutoramento. O aparato completo compõe-se da seção de testes, circuito de água, sistema de compressão e condicionamento de ar, e seções de injeção dos fluxos e condicionamento do escoamento. A seção de testes consiste em um banco de tubos distribuídos segundo configuração triangular normal, com os tubos apresentando diâmetro externo de 19,1 mm, comprimento de 381 mm, e espaçamento transversal de 24 mm. Os experimentos foram realizados para escoamento vertical ascendente de misturas água-ar e velocidades superficiais da fase líquida e gás de 0,020 a 1,500 m/s e de 0,10 a 10,00 m/s, respectivamente. Neste estudo foram desenvolvidas técnicas inéditas para determinação experimental da fração de vazio superficial no interior do banco de tubos baseadas em sistemas óptico e de sensoriamento capacitivo. Os padrões de escoamento foram identificados subjetivamente através de visualização de imagens e vídeos do escoamento, e objetivamente com o auxílio do método de agrupamento de dados k-means utilizando parâmetros baseados no sinal de perda de pressão e do sensoriamento capacitivo. Identificou-se subjetivamente os padrões de escoamento bolhas, bolhas dispersas, bolhas grandes, agitante, intermitente e anular. Constatou-se equivalência entre os padrões de escoamento identificados através dos métodos objetivo e subjetivo. Resultados experimentais para fração de vazio foram obtidos através de técnicas óptica e capacitiva. Constatou-se que o traçador rodamina B utilizado no método óptico altera as condições do escoamento, ainda que em concentrações reduzidas. A partir dos resultados obtidos com o sensoriamento capacitivo estimou-se a fração de vazio para o padrão bolhas. Resultados para a parcela friccional da perda de pressão também foram levantados. Constata-se o incremento da fração de vazio e da parcela friccional da perda de pressão com as velocidades superficiais das fases líquida e gás. Os resultados para fração de vazio foram comparados com métodos de previsão da literatura, e de maneira geral os métodos preveem as tendências dos resultados experimentais apenas para vazões de líquido reduzidas. Analogamente, os resultados para perda de pressão foram comparados com estimativas segundo métodos da literatura, concluindo que os métodos não preveem satisfatoriamente os resultados obtidos. Desta forma, foram propostos novos métodos de previsão para padrões de escoamento, fração de vazio e parcela friccional da perda de pressão, desenvolvidos a partir de análises dos mecanismos dominantes do escoamento, e adotando parâmetros adimensionais para correlacionar os dados. Os métodos propostos preveem satisfatoriamente os resultados experimentais deste estudo e da literatura para escoamentos bifásicos água-ar. / The present thesis concerns a theoretical and experimental study of external two-phase flows across tube bundles. Initially, a comprehensive literature review covering flow patterns, void fraction and pressure drop for single and two-phase flows across tubes bundle is presented. The review also describes predictive methods for these parameters. A comparison of these methods reveals reasonable disagreement among their predictions, indicating the absence of generalized methods. Subsequently, the apparatus and instrumentation designed and built to obtain the experimental data are described. The experimental apparatus comprises the test section, a water loop, air compression and conditioning systems, and sets for fluid flow injections and conditioning. The test section is a normal triangular tube bundle, with 19.1 mm OD tubes, 381 mm long and transversal pitch of 24 mm. The experiments were performed for air-water upward vertical flow, for superficial liquid and gas velocities ranging from 0.020 to 1.500 m/s and 0.10 to 10.00 m/s, respectively. Innovative techniques to evaluate the void fraction within the bundle were developed based on capacitive and optical methods. The flow patterns were identified subjectively and objectively by k-means clustering method, using as clustering parameters the pressure drop and the capacitive signals. Bubbles, dispersed bubbles, large bubbles, churn, intermittent and annular flow patterns were identified subjectively. The data groups identified by the objective method are representative of the flow patterns. Void fraction measurements were obtained for bubbly flow using both techniques (optical and capacitive). The void fraction data based on the optical method had its experimental range limited due to changes in the flow characteristics caused by the addition of the fluorescent dye Rhodamine B. The experimental results indicate that the void fraction increases with increasing the superficial velocities of both phases. In general, the void fraction predictive methods available in the literature capture the trends of the experimental results only for reduced liquid flow rates. According to the experimental results, the frictional pressure drop increases asymptotically with increasing the flow rates of both phases. None of the predictive methods from literature evaluated in the present study predicted satisfactorily the experimental results. Methods for prediction of flow patterns, void fraction and frictional pressure drop parcel were also developed in the present study. These methods provided reasonable predictions of the experimental results obtained in the present study, and also from the literature for air and water flows across tube bundles.

Drilled shafts are cylindrical, cast-in-place concrete deep foundation elements. Their construction involves drilled excavation of soil or rock using large diameter augers, and placement of the necessary reinforcing steel in the excavation followed by concreting. Where a high water table is encountered, drilling slurry is used to support the excavation walls and concreting is tremie-placed. Even though the history of drilled shaft construction goes back to the 1950s, the occurrence of anomalies persists in the form of soil inclusions, reduction in shaft cross-sectional area and exposure of reinforcement. One of the main reasons for the anomalies is attributed to the kinematics of concrete flowing radially from within the reinforcing cage to the surrounding annulus/concrete cover region. In view of this radial component of concrete flow and thus radially flowing interfaces between the concrete and slurry, the region outside the cage is more likely to contain veins of poorly cemented or high water-cement ratio material. These veins contain trapped slurry, which oftentimes consists of bentonite, jeapordizing the integrity of the shafts. This research program focuses on the numerical evaluation of self-consolidating concrete (SCC) for drilled shaft application by taking into account realistic non-Newtonian concrete flow properties and the shaft structural blockages. For this objective, a 3-D computational fluid dynamics (CFD) model of the concrete flow in the shaft excavation is developed in ANSYS-Fluent. As a precursor to 3-D modeling, 2-D CFD modeling is carried out using COMSOL Multiphysics. In both 2-D and 3-D models, the Volume of Fluid method is used for computing the motion of the interface between the concrete and the drilling slurry. The models predict the flow patterns and volume fraction of concrete and slurry. The results are encouraging as the flow pattern from the simulation shows both horizontal and vertical creases in the concrete cover region. Moreover the flow pattern shows the concrete head differential developed between the inside and the outside the reinforcement cage. Further, the 3-D model is evaluated by studying the influence of the size of drilled shaft and arrangement of the bars and the results obtained are realistic. With this 3-D model developed as a tool, the simulation of SCC and the normal standard concrete (NC) flow in drilled shaft concreting are studied in terms of creases and concrete head differential encountered in the flow. From the simulation, it is observed that in the flow pattern of SCC, the creases are very few compared to the one obtained from the flow pattern of NC. Moreover, the concrete head differential in the flow pattern of SCC is much less, than the head differential obtained from the flow pattern of NC flow. In the case of SCC, the head differential encountered about one inch. In the case of NC, the concrete head differential is 4-inch when the vertical rebars are spaced at 7-inch apart and 10-inch when the rebars are placed at 3.5-inch apart. Based on this numerical evaluation of SCC flow in the drilled shaft excavation, it is concluded that the performance of SCC is better than the performance of NC in filling the cover annular region of drilled shafts.