Electrohydrodynamics in convective boiling heat exchangers has been studied since the early 1990’s and has been shown to result in a large variation in the average performance enhancement of these systems. The behaviour of EHD assisted convective boiling heat exchangers, is still largely unpredictable owing to a number of conflicting parameters which are rarely kept constant in empirical studies, i.e. flow pattern and heat flux. In this thesis, it is hypothesised that by reducing the number of confounding variables in the experimental test conditions, and understanding the behaviour of EHD in convective boiling systems from a flow pattern dependent point of view, this can allow for the development of flow pattern dependent experimental correlations & numerical models to develop a methodology for performance prediction, control strategies and system integration for an EHD assisted convective boiling heat exchange device. A 30 cm long, smooth, concentric, annular test section is used to analyse the effect of EHD on convective boiling performance under constant flow pattern, constant, low heat flux, and negligible free charge conditions. Saturated boiling conditions for flow-rates between 60 kg/m2s and 180 kg/m2s and thermodynamic quality range of 0.25 - 0.55 were tested. Heat transfer enhancement ranged from 0.95 to 2.3 fold and pressure drop penalty varied from 1.4-3 fold over these test conditions. The local EHD behaviour was found to be more consistent along the axial length of the test section compared to empirical data in the literature, which uses much longer test section lengths, where flow pattern can vary. An experimental database of EHD convective boiling data for horizontal annular electrode geometries was compiled to be used for analysis purposes. The performance of the heat exchanger in both free-field and high voltage conditions could be explained by looking at the flow patterns in each case. Electrostatic modelling was used to determine electric field strength distributions and interfacial stress due to the dielectrophoretic and electrostriction forces on the liquid vapour interface, which induce liquid extraction based flow pattern re-distribution in two phase dielectric flows. A fully coupled 2D, adiabatic numerical model for the effect of the electric body force on two phase flow pattern distribution was developed. Charge was neglected in this model. Two different models for the interface migration were used and compared; a moving mesh (MM) interface tracking model and a volume of fluid (VOF) interface tracking mode. Both were verified against published experimental data. For the liquid extraction verification case, the VOF model suffers interface stretching up to 300% resulting in a 42% slower extraction time and underestimated forces. However, it is useful to use the VOF model when simulating complex flow patterns which are subject to topological changes like bubble detachment or droplet coalescence as these cannot be simulated with the moving mesh model. The moving mesh model can be used to determine the error in forces and phase velocities when using the VOF model. A methodology for generating two-phase EHD flow pattern maps was developed by incorporating the electric Froude number into each of the flow pattern transition equations. A semi-analytical model was developed to determine the maximum interfacial stress due to EHD for stratified flows to reduce the requirement of numerical modeling, and thus the flow pattern map generation methodology is fully equation based. Although transition equations developed by multiple researchers were used and compared, it is recommended that the Steiner transitions equations be used for EHD two-phase flow pattern mapping, until more fundamental experimental data can be gathered to modify the semi-empirical transition equations used in more state-of-the-art maps. EHD was found to significantly affect the “stratified-stratified wavy (SSW)” and “stratified wavy – intermittent/annular (SWIA)” transitions for concentric horizontal geometries, with minimal effect on the transition to dryout and no effect on the “intermittent dispersed bubbly (IB)” transition. The EHD flow pattern maps were generated and compared against data from the present study and a database of experimental EHD convective boiling studies. The regions where maximum enhancement were seen in the literature correlate well with those regions predicted by the maps. Performance correlations for the EHD convective boiling heat transfer and pressure drop were developed. They are based on the free-field Kandlikar correlation [1] for two-phase heat transfer and the Chisholm-Laird [2] correlations for two-phase pressure drop, respectively. The EHD flow pattern map is used to determine what the flow pattern for a given applied voltage will be, and flow pattern based enhancement linear multipliers are then used to determine the EHD performance above the free-field case. EHD is a form of active enhancement, i.e. it requires power. Thus, it would be used in systems that require performance control or regulation, in addition to some niche applications like space where it can be used instead of gravity. A method for EHD controller design was established and an EHD control algorithm was designed and implemented on the test section for the flow pattern and applied waveforms that were determined to be optimal to maximize enhancement in this geometry. System identification was performed empirically to determine the transfer function between EHD voltage and heat load to be controlled for. This resulted in a 1st order plus dead-time model to which proportional-integral controller constants were tuned. Two controllers were developed; a PID control system and a Smith model predictive control system and these were compared based on their ability to regulate the output quality of the heat exchanger when subject to dynamic heat loading. Regulation was achieved for a dynamic heat load within ±25% bound from the designed steady state load. These controllers operate on one flow pattern as the test section is 30 cm long. Flow pattern dependent controller design would be required for a full length convective boiling heat exchanger. / Thesis / Doctor of Philosophy (PhD) / Control of boiling heat transfer using electric fields is hard to predict. This thesis presents a set a design guidelines based on how the electric field enhances the flow pattern.

This research addresses the critical need for precise two-phase flow data in the development of computer simulation models, with a specific focus on the annular flow regime's droplet behavior. The study aims to contribute to the evaluation of safety and efficiency in nuclear reactors that handle fluids transitioning between liquid and gas states for thermal energy transport. Central to the investigation is the collection and analysis of droplet size and velocity distribution data, particularly to help with developing models for the water-cooled nuclear power plants. The experimental setup employs advanced tools, including a high-speed camera, lens, teleconverter, and a selected light source, to capture high-resolution images of droplets. Calibration procedures, incorporating depth of field testing, are implemented to ensure accurate droplet size measurements. A critical component of the research is the introduction of a droplet identification program, developed using Matlab, which facilitates efficient processing of experimental data. Preliminary results from the Virginia Tech test facility demonstrate the system's capability to eliminate out-of-focus droplets and obtain precise droplet data in a reasonable amount of time. Experimental results from the Rensselaer Polytechnic Institute test facility provide droplet size and velocity distributions for a variety of annular flow conditions. This facility has a concurrent two-flow system that pumps air and water at different rates through a 9.525 mm inner diameter tube. The conditions tested include gas superficial velocities ranging from 22 to 40 m/s and liquid superficial velocities ranging from 0.09 to 0.44 m/s. The measured flow has a temperature of 21°C and a pressure of 1 atm. / Master of Science / This research explores the behavior of small droplets as fluids transition between liquid and gas states, particularly within the context of the cooling water in nuclear power plants. The overarching goal is to collect data on these droplets to improve computer simulations that help design nuclear reactors and assess their safety. This is important because it is often infeasible due to safety, monetary, or time restrictions to physically test some nuclear reactor equipment. The study employs state-of-the-art technology, including high-speed cameras and specialized imaging tools, to capture and analyze droplet size distribution data. This investigation is pivotal in ensuring the fuel in nuclear reactors remain adequately cooled during part of the boiling process. The research methodology includes the development of a droplet identification program using Matlab, ensuring efficient processing of experimental data. Preliminary findings from experimental tests at Virginia Tech showcase the program's capability to filter out irrelevant data and provide accurate droplet information. Experimental results from the Rensselaer Polytechnic Institute annular flow test facility provide droplet size and velocity data for a range of conditions that cooling water may face. Beyond its contributions to nuclear engineering, this research holds promise for influencing advancements in various applications that involve liquid droplets, opening avenues for innovation in the broader scientific and engineering communities.

External loop airlift reactors are widely used for biochemical applications such as syngas fermentation and wastewater treatment. To further understand the inherent gas-liquid flow physics within the reactors, computational modeling and simulations of hydrodynamics for air-water external loop airlift reactors were investigated. The gas-liquid flow dynamics in a bubble column were simulated using a FORTRAN code developed by Los Alamos National Laboratory, CFDLib, which employs an Eulerian-Eulerian ensemble averaged method. A two-dimensional Cartesian coordinate system was used to conduct an extensive grid resolution study; it was found that grid cells smaller than the bubble diameter produced unstable solutions. Next, closure models for drag force and turbulent viscosity were investigated for a simple bubble column geometry. The effects of using a bubble pressure model and two drag coefficient models, the White model and the Schiller-Naumann model, were investigated. The bubble pressure model performed best for homogeneous (low velocity) flows and the Schiller-Naumann model was best for all flow regimes. Based on the studies for bubble column flows, an external loop airlift reactor was simulated using both two- and three-dimensional coordinates and results for gas holdup and riser velocity agreed better with experimental data for the 3D simulations. It was concluded that when performing 2D and 3D simulations, care must be taken when specifying the effective bubble diameter size, especially at high flow rates. Population balance models (PBM) for bubble break-up and coalescence were implemented into CFDLib, validated with experiments, and simulated for the external loop airlift reactor at high inlet superficial gas velocities. The PBM predictions for multiple bubble sizes were comparable with the single bubble size simulations; however, the PBM simulations better predicted the formation of the gas bubble in the downcomer. The 3D PBM simulations also gave better predictions for the average bubble diameter size in the riser. It was concluded that a two-dimensional domain is adequate for gas-liquid flow simulations of a simple bubble column geometry, whereas three-dimensional simulations are required for the complex airlift reactor geometry. To conclude, a two-fluid Eulerian-Eulerian model coupled with a PBM is needed for quantitative as well as physical predictions of gas-liquid external loop airlift reactor flows at high inlet superficial gas velocities. / Ph. D.

Geological carbon dioxide (CO2) sequestration (GCS) in deep saline aquifers is viewed as a viable solution to dealing with the impact of anthropogenic CO2 emissions on global warming. The trapping mechanisms that control GCS include capillary trapping, structural trapping, dissolution trapping, and mineral trapping. Wettability and density-driven convection play an important role in GCS, because wettability significantly affects the efficiency of capillary trapping, and density-driven convection greatly decreases the time scale of dissolution trapping. This work focuses on the role of wettability on multiphase flow in porous media, density-driven convection in porous media, and their implications for GCS in deep saline aquifers. Wettability is a critical control over multiphase fluid flow in porous media. However, our understanding on the wettability heterogeneity of a natural rock and its effect on multiphase fluid flow in a natural rock is limited. This work innovatively models the heterogeneous wettability of a rock as a correlated random field. The realistic wetting condition of a natural rock can be reconstructed with in-situ measurements of wettability on the internal surfaces of the rock. A Bentheimer sandstone was used to demonstrate the workflow to model and reconstruct a wettability field. Relative permeability, capillary pressure-water saturation relation are important continuum-scale properties controlling multiphase flow in porous media. This work employed lattice Boltzmann method to simulate the displacement process. We found that pore-scale surface wettability heterogeneity caused noticeable local scCO2 and water redistributions under less water-wet conditions at the pore scale. At the continuum scale, the capillary pressure-water saturation curve under the heterogeneous wetting condition was overall similar to that under the homogeneous wetting condition. This suggested that the impact of local wettability heterogeneity on the capillary pressure-water saturation curve was averaged out at the entire-sample scale. The only difference was that heterogeneous wettability led to a negative entry pressure at the primary drainage stage under the intermediate-wet condition. The impact of pore-scale wettability heterogeneity was more noticeable on the relative permeability curves. Particularly, the variation of the scCO2 relative permeability curve in the heterogeneous wettability scenario was more significant than that in the homogenous wettability scenario. Results showed that higher wettability heterogeneity (i.e., higher standard deviation and higher correlation length) increased the variations in the CO2/brine relative permeability curves. Dissolution of CO2 into brine is a primary mechanism to ensure the long-term security of GCS. CO2 dissolved in brine increases the CO2-brine solution density and thus can cause downward convection. Onset of density-driven instability and onset of convective dissolution are two critical events in the transition process from a diffusion-dominated regime to a convection-dominated regime. In the laboratory, we developed an empirical correlation between light intensity and in-situ solute concentration. Based on the novel and well-controlled experimental methods, we measured the critical Rayleigh-Darcy number and critical times for the onset of density-driven instability and convective dissolution. To further investigate the impact of permeability heterogeneity on density-driven convection, a three-dimensional (3D) fluidics method was proposed to advance the investigation on density-driven convection in porous media. Heterogeneous porous media with desired spatial correlations were efficiently built with 3D-printed elementary porous blocks. In the experiments, methanol-ethylene-glycol (MEG), was used as surrogate fluid to CO2. The heterogeneous porous media were placed in a transparent tank allowing visual observations. Results showed that permeability structure controlled the migration of MEG-rich water. Permeability heterogeneity caused noticeable uncertainty in dissolution rates and uncertainty in dissolution rates increases with correlation length. To sum up, this work comprehensively employed novel experimental methods and large-scale direct simulations to investigate the sequestration of CO2 in saline aquifers at a pore scale and a continuum scale. The findings advanced our understanding on the role of wettability heterogeneity and permeability heterogeneity on GCS in deep saline aquifers. / Doctor of Philosophy / Global warming caused by anthropogenic CO2 emissions is a pressing issue to address of our time. The storage of CO2 in deep saline aquifers is a promising solution because of saline aquifers' vast storage capacity. Property heterogeneity exists extensively in saline aquifers from a continuum scale to a pore scale. The implications of pore-scale wettability heterogeneity and continuum-scale permeability heterogeneity for the storage of CO2 in saline aquifers are not clear. This work is to employ novel experimental methods and powerful simulation tools to investigate the role of wettability heterogeneity and permeability heterogeneity on the storage of CO2 in saline aquifers. This work measured contact angles on the scanned micro-CT images of a Bentheimer sandstone after a CO2 flooding. A correlated lognormal wettability model was put forward with the statistical information of the contact angle measurements. Simulations on the CO2/brine flow in the Bentheimer sandstone were performed. Results showed that the wettability heterogeneity caused noticeable redistributions of CO2/brine compared to scenarios under homogeneous wettability. Impact of wettability on capillary pressure-water saturation curve was not noticeable because the effects were averaged out through the entire rock sample. The standard deviation and correlation length caused variations on the relative permeabilities. This means that we need to take them into consideration in simulating the migration of CO2 in saline aquifers at a reservoir scale. After CO2 pools beneath the impermeable cap rock, dissolution of CO2 into brine dominates the trapping process. Convection caused by CO2 dissolution can greatly accelerate the dissolution rate. The onset of convection is a critical issue and lack of experimental evidence. This work firstly determined the onset time of instability. To further investigate the heterogeneity on the convection, this work proposed a 3D-print-based method to efficiently build heterogeneous porous media with a designed permeability distribution. The experiments were conducted, and results showed that heterogeneity structure of porous media can cause great variations on the dissolution rate of CO2. The findings of this work advanced our understanding on the migration of CO2 in saline aquifers, provided solid basis for assessment and decision on the storage of CO2 into saline aquifers.

Void fraction is an essential parameter for understanding the interfacial structure, and heat and mass transfer mechanisms in various gas-liquid flow systems. It becomes critically important to accurately measure void fraction as advanced high fidelity two-phase flow models require high-quality validation data. However, void fraction measurement remains a challenging task to date due to the complexity and rapid-changing characteristic of the gas-liquid boundary flow structure. This study aims to develop an advanced void fraction measurement system based on x-ray and fast line detector technologies. The dissertation has covered the major components necessary to develop a complete measurement system. Spectral analysis of x-ray attenuation in two-phase flow has been performed, and a new void fraction model is developed based on the analysis. The newly developed pixel-to-radial conversion algorithm is capable of converting measured void fraction along with the detector array to the radial distribution in a circular pipe for a wide range of void fraction conditions. The x-ray system attains the radial distributions of key measurable factors such as void fraction and gas velocity. The data are compared with the double-sensor conductivity probe and gas flowmeter for various flow conditions. The results show reasonable agreements between the x-ray and the other measurement techniques. Finally, various 2-D tomography algorithms are implemented for the non-axisymmetric two-phase flow reconstruction. A comprehensive summary of classical absorption tomography for the two-phase flow study is provided. An in-depth sensitivity study is carried out using synthetic bubbles, aiming to investigate the effect of various uncertainty factors such as background noise, off-center shift, void profile effect, etc. The sensitivity study provides a general guideline for the performance of existing 2-D reconstruction algorithms. / Doctor of Philosophy / Gas-liquid flow phenomenon exists in an extensive range of natural and engineering systems, for example, hydraulic pipelines in a nuclear reactor, heat exchanger, pump cavitation, and boilers in the gas-fired power stations. Accurate measurement of the void fraction is essential to understand the behaviors of the two-phase flow phenomenon. However, measuring void fraction distribution in two-phase flow is a difficult task due to its complex and fast-changing interfacial structure. This study developed a comprehensive suite of the non-intrusive x-ray measurement techniques, and a pixel-to-radial conversion algorithm to process the line- and time-averaged void fraction information. The newly developed algorithm, called the Area-based Onion-Peeling (ABOP) method, can convert the pixel measurement to the radial void fraction distribution, which is more useful for studying and modeling axisymmetric flows. Various flow conditions are measured and evaluated for the benchmarking of the algorithm. Finally, classical 2-D reconstruction algorithms are investigated for the void fraction measurement in non-axisymmetric flows. A comprehensive summary of the performance of these algorithms for a two-phase flow study is provided. An in-depth sensitivity study using synthetic bubbles has been performed to examine the effect of uncertainty factors and to benchmark the algorithms for the non-axisymmetric flows.

This thesis presents an investigation into the motions of the horseshoe vortices and the passage vortex, within a plane turbine blade cascade. Fluid motion was determined using a tracer gas technique. Ethylene was injected into the pressure-side and suction-side legs of the horseshoe vortex, near the leading edge of the cascade. Ethylene concentrations were determined at two downstream locations using a flame ionization detector. It was found that the pressure-side leg of the horseshoe vortex moved toward the suction side of the passage, starting the formation of the passage vortex, and was distributed throughout the passage vortex. The suction-side leg of the horseshoe vortex convected once around the periphery of the passage vortex before passing the cascade trailing edge. Downstream of the trailing edge, most of the fluid from the suction-side leg diffused into the passage vortex. However, twice as much fluid from the suction-side leg, as opposed to the pressure-side leg, mixed within the blade wake. At a location 40% of the axial chord downstream of the trailing edge, the passage vortex (shown previously to account for 60% of the overall total pressure losses) contained over 65% of the fluids from both legs. / M.S.

Biomass is considered a biorenewable alternative energy resource that can potentially reduce the use of natural gas and provide low cost power production or process heating needs. Biomass hydrodynamics in a fluidized bed are extremely important to industries that are using biomass material in gasfication processes to yield high quality producer gas. However, biomass particles are typically difficult to fluidize due to their peculiar shape and a second inert material, such as sand, is typically added to the bed. The large differences in size and density between the biomass and inert particles lead to nonuniform distribution of the biomass within the fluidized bed, and particle interactions and mixing become major issues. The main goal of this research was to use CFD as a tool for modeling and analyzing the hydrodynamic behavior of biomassas a single material or as part of a mixture in a fluidized bed. The first part of this research focused on the characterization of biomass particles in a fluidized bed and validation of a numerical model with experimental results obtained from pressure measurements and CT and X-ray radiograph images. For a 2D fluidized bed of glass beads, the pressure drop, void fraction and mean bed height expansion were in quantitative agreement between the experiments and simulations using Syamlal-O'Brien and Gidaspow drag models. It was encouraging that the Gidaspow model predictions were in close agreement because the model does not require knowing the minimum fluidization as an input. Ground walnut shells were used to represent biomass because the material fluidizes uniformly and is classified as a Geldart type B particle. Two-dimensional simulations of ground walnut shells were analyzed to determine parameters that cannot easily be measured experimentally. The parametric study for ground walnut shell indicated that the material can be characterized with a medium sphericity (~0.6) and a relatively large coefficient of restitution (~0.85). In the second part of this work numerical simulations of a ground walnut shell fluidizing bed with side air injection were compared to CT data for the gas-solid distribution to demonstrate the quantitative agreement for bed fluidization. The findings showed that 2D simulations overpredicted the fluidized bed expansion and the results did not demonstrate a uniformly fluidizing bed. The 3D simulations compared well for all cases. This study demonstrates the importance of using a 3D model for a truly 3D flow in order to capture the hydrodynamics of the fluidized bed for a complicated flow and geometry. Finally, CFD modeling of pressure fluctuations was performed on sand and cotton-sand fluidized beds operating at inlet velocities ranging from 1.0-9.0Umf with the objective of predicting characteristic features of bubbling, slugging, and turbulent fluidization regimes. It was determined that the fluidized bed can be modeled using MUSCL discretization and the Ahmadi turbulence model. Three-dimensional sand fluidized beds were simulated for different fluidization regimes. Fluidized beds for all the regimes behaved as second-order dynamic systems. Bubbling fluidized beds showed one broad peak with a maximum at 2.6 Hz while slugging and turbulent showed two distinct peaks. It was observed that the peak at low frequency increased in magnitude as the flow transitioned from a slugging to a turbulent fluidization regime. CFD simulations of fluidized beds with the purpose of studying pressure fluctuations have demonstrated to be a useful tool to obtain hydrodynamic information that will help determine the fluidization regime. Prediction of slugging and turbulent fluidization regimes using CFD have not been reported to date. The work presented here is the first of its kind and can be an important advantage when designing a reactor and evaluating different operation conditions without the need to test them in a pilot plant or a prototype. / Ph. D.

We present 2-D numerical simulations of two phase flow in seafloor hydrothermal systems using the finite control volume numerical scheme FISHES. The FISHES code solves the coupled non-linear equations for mass, momentum, energy, and salt conservation in a NaCl-H2O fluid to model the seafloor hydrothermal processes. These simulations use homogeneous box geometries at a fixed seafloor pressure of 25 MPa with constant bottom temperature boundary conditions that represent a sub-axial magma chamber to explore the effects of permeability, maximum bottom temperature and system depth on the evolution of vent fluid temperature and salinity, and heat output. We also study the temporal and spatial variability in hydrothermal circulation. The two-phase simulation results show that permeability plays an important role in plume structure and heat output of hydrothermal systems, but it has little effect on vent fluid temperature and salinity, given the same bottom temperature. For some permeability values, multiple plumes can vent at the seafloor above the simulated magma chamber. Temporal variability of vent fluid temperature and salinity and the complexity of phase separation suggest that pressure and temperature conditions at the top of the axial magma chamber cannot be easily inferred from vent fluid temperature and salinity alone. Vapor and brine derived fluids can vent at the seafloor simultaneously, even from neighboring locations that are fed by the same plume. / Master of Science

Biopharmaceutical manufacturing is a rigorous and expensive process. Due to the medicinal nature of the product, a high purity level is required and several expensive purification steps must be utilized. Cost-effective production and purification is essential for any biopharmaceutical product to be successful and development of the fastest, most economical, and highest-yielding purification scheme is a constant engineering challenge. Commercial-scale purification schemes currently revolve around the use of multiple chromatography steps for the purification of biopharmaceutical products. Chromatography has many shortcomings including high cost, limited throughput, and complex scale up. The goal of this research was to develop an alternative, non-chromatography purification step for the separation of an acidic model protein, recombinant β-glucuronidase (rGUS), from transgenic tobacco with high yield and purity. Aqueous two-phase extraction (ATPE) is a powerful technique for separation and purification of proteins, and has the potential to replace an expensive chromatography step for the initial purification of recombinant proteins. ATPE enables high levels of target protein recovery and concentration while removing large amounts of impurities from the initial extract. Fractional factorial designs and response surface methodology were used to determine an optimized aqueous two-phase system for the purification of rGUS from transgenic tobacco. In a 13.4 % (w/w) PEG/18% (w/w) potassium phosphate system, 74% of the rGUS was recovered in the top PEG-rich phase while 90% of the native tobacco proteins were removed in the interphase and the bottom phase. A purification factor of about 20 was achieved in this process. / Master of Science

Fluidized bed reactor technology has been widely adopted within the industry as vital component for numerous manufacturing, power generation and gasification processes due to its enhanced mixing characteristics. Computational modeling of fluidized bed hydrodynamics is a significant challenge that has to be tackled for increasing predictive accuracy. The distributor plate of a fluidized bed is typically modeled using a uniform inlet condition, when in reality the inlet is non-uniform inlet. The regions of bed mass that do not fluidize because of the non-uniform inlet conditions form deadzones and remain static between the jets. A new model based on the mass that contributes to the pressure drop is proposed to model a fluidized bed, and has been investigated for a cylindrical reactor for glass beads, ceramic single solids particles, and glass-ceramic, and ceramic-ceramic binary mixtures. The adjusted mass model was shown to accurately predict fluidization characteristics such as pressure drop and minimum fluidization velocity. The effectiveness of the adjusted mass model was further illustrated by applying it to fluidized beds containing coal, switchgrass, poplar wood, and cornstover biomass particles and coal-biomass binary mixtures. The adjusted mass model was further analyzed for bed expansion heights of different mixtures, and for solids distribution by analyzing the solids volume fraction. The effect of increasing the percent biomass in the mixture was also investigated. To further model the non-uniform inlet condition, two different distributor configurations with 5 and 9 jets was considered for a quasi-2D bed, and simulations were performed in both 2D and 3D. Fluidization characteristics and mixing of the bed were analyzed for the simulation. Furthermore, the deadzones formed due to multiple jet configurations of the distributor are quantified and their distributions over the plate were analyzed. / Ph. D.