VTechWorks Repository :: Browsing by Author "Coutier-Delgosha, Olivier"

Browsing by Author "Coutier-Delgosha, Olivier"

Now showing 1 - 20 of 23

Aerodynamic Interactions in Vortex Tube Separator Arrays
Acharya, Aditya Sudhindra (Virginia Tech, 2023-06-22)
Helicopter turboshaft engines may ingest large amounts of foreign particles (most commonly sand/dust), which can cause significant compressor blade damage and even engine failure. In many helicopters, this issue is mitigated by separating the particles from the intake airstream. An effective device for engine air-particle separation is the vortex tube separator (VTS), which uses centrifugal forces in a vortical flow to radially filter foreign particles from a duct with an annular exit. Dozens or hundreds of these devices are linked together on a shared manifold known as a VTS array. There is a distinct lack of scientific literature regarding these arrays, which likely feature significantly more complex flowfields than singular VTSs due to aerodynamic interactions between the devices. The research presented in this dissertation identifies and explains flow features unique to arrays by means of an experimental investigation downstream of various VTS configurations in a wind tunnel. Mean PIV flowfields reveal that the VTS array rapidly generates a strong central recirculation zone while a single VTS does not, implying the existence of axial flow gradients within associated separators that could affect filtration efficiency. The key factor here is the global swirl intensity, which is increased in array flows due to high angular momentum contributions from separators that are radially distant from the duct center. A preliminary momentum integral model is constructed to predict the onset of recirculation in VTS flows. Analysis is then extended to the unsteady flowfield, where it is shown that VTS-generated turbulence contains only low levels of anisotropy. Spectral proper orthogonal decomposition is conducted on the array flow; it reveals the existence of low-frequency harmonic behavior composed of back-and-forth pumping motions downstream of the central VTS. Additionally, a unique precession motion is found in the same region at a slightly higher frequency. Similar precessing vortex cores have been shown to reduce separation efficiency in other cyclone separators. Both of these coherent structures may be associated with the central recirculation zone and may interfere with VTS array filtration given their timescales relative to potential particle relaxation timescales. This dissertation opens the door for future experimental and computational studies of fluid and particle dynamics in VTS flows with the goal of improving VTS array-specific design philosophies.
An Analysis on Hydrodynamic Loads for Surface-Piercing Propellers Using Computational Fluid Dynamics
Brookshire, Kaleb (Virginia Tech, 2022-07-18)
A surface piercing propeller (SPP) is a propeller that is partially submerged in water and is considered a possible solution to high-speed vessels (greater than 50 knots) where cavitation plays a vital role due to its ever-increasing detrimental effects. Computational Fluid Dynamics (CFD) has become a more prevalent solution in recent years due to lower costs and the ability to evaluate varying setups. However, Computational Fluid Dynamics has had problems accurately solving the hydrodynamic loads for an SPP as recently as a few years ago. Accurately predicting these loads is of great importance because it will allow future simulations to add more effects such as cavitation, shaft inclination effects, multiple propellers, and fluid-structure interaction. Using FINE/Marine, a CFD software specifically designed for marine applications simulations with the 841-B SPP model and changing the Froude number (Fn) and advance coefficient (J), an in-depth validation process and extending upon previous results found when combining CFD and surface-piercing propellers was performed. Several cases between J = 0.6 to J = 1 and Fn = 2 to Fn = 6 are first performed to validate the models against experiments, then more complex features such as multiple propellers and shaft inclination angles were included to extend upon previous work of CFD for surface-piercing propellers. This analysis of the results suggests that CFD models could genuinely be validated against current experimental setups, and therefore more complex additions could also be made and with stronger accuracy than in previous years.
Computational Modeling of Droplet Impact Dynamics on Solid Substrates
Saravanan Manikkam, Pratulya Rajan (Virginia Tech, 2023-01-31)
A computational model is developed to simulate the impact dynamics of a droplet on solid substrates with the purpose of predicting the droplet spreading characteristics over time. Previous studies focused on finding relations between the impact parameters and outcome dynamics. A modified approach like the one used in this project revolves around modeling the moving contact lines at the interface in a multiphase flow environment. Focusing on research from an aircraft de-icing point of view, this study is considered a prerequisite in understanding the physics of droplet impact. The primary focus is on extending the application to incorporate super-cooled environments. Development of the model involved the use of the Volume-of-Fluid function coupled with the High-Resolution Interface Capturing scheme to model the moving contact line. The evolution of the moving contact line is modeled with contact angles as their inputs to understand the effect of the surface tension forces. Contact angle modeling is based on the Blended-Kistler method, which captures the contact angle evolution based on the surface tension and capillary number. Preliminary validation performed on the model proves its effectiveness in accurately simulating the impact behavior when compared to the literature, where the spread diameter and height agree well with experiments. The validated model is also compared to the in-house experiments performed at the Cavitation and Multiphase flow laboratory using different substrate materials. The substrates each show unique behavior - Impact on Glass results in the droplet depositing on the surface. Aluminum results in a full rebound and PET-G, results in a drop ejection. Based on inputs from the experiments - contact angles, spread diameter, and the maximum spread $beta$, show good agreement in comparison to the literature.
Data Analysis of an Unsteady Cavitating Flow on a Venturi-type Profile
Nemati Kourabbasloo, Navid (Virginia Tech, 2021-12-01)
The instability modes and non-linear behavior of a cavitating flow have been studied based on the experimental data obtained from planar Particle Image Velocimetry (PIV). Three data-driven techniques, Proper Orthogonal Decomposition (POD), Dynamic Mode Decomposition (DMD), and Clustered-based Reduced Order Modeling (CROM), are applied to the snapshots of the fluctuating component of velocity to investigate instability modes of the cavitating flow. DMD and POD analysis yield multiple modes are corresponding to slow-varying drift flow, cloud-shedding, and Kelvin-Helmholtz (KH) instability for a fixed inlet flow condition. The high coherence measure obtained from the instabilities suggests a transfer of energy from the largest scales, fluctuating mean flow, to the smaller scales such as cloud cavitation and Kelvin-Helmholtz (KH) instability. It is demonstrated that the POD decorrelation of length scales yields inherently quasi-periodic time dynamics, e.g., incommensurate frequencies. Moreover, the eigenvalue obtained from DMD revealed multiple harmonic with different decay rates associated with the cloud cavitation. The above-mentioned intermittent transition between distinct cloud shedding regimes is investigated via Clustered-based Reduced Order Modeling (CROM). Four aperiodic shedding regimes are identified. 68% of the time, triplets of vortices are formed, while 28% of the time, a pair of vortices are formed in the near wake of the throat. Dominant mechanisms governing the momentum transport and the turbulence kinetic energy production, destruction, and redistribution in distinct regions of the flow field have been identified using Gaussian Mixture Models (GMMs). The preceding data-driven techniques and in-depth analysis of the results facilitated modeling of the cavitation inception and break-up. Accordingly, a phase transition field model is developed using the ultra-fast Time-Resolved Particle Image Velocimetry (TR-PIV) and vapor void fraction spatial and temporal data. The data acquisition is implemented in a Venturi-type test section. The approximate Reynolds number based upon the throat height is 10,000, and the approximate cavitation number is 1.95. The non-equilibrium cavitation model assumes that the phase production and destruction are correlated to the static pressure field, pressure spatial derivatives, void fraction, and divergence of the velocity field. Finally, the dependence of the model on the empirical constants has been investigated.
Droplet Rebound and Atomization Characteristics of Vibrating Surfaces
Kendurkar, Chinmay (Virginia Tech, 2023-02-28)
Icing on aircraft wings is one of the leading causes of aircraft crashes. It is mainly caused due to accumulation of ice or snow on the wing surface due to impact with supercooled droplets when passing through clouds at high altitudes, causing loss of lift obtained by the wings. It was found that droplet impact characteristics are dependent on droplet size, surface roughness, surface material hydrophobicity, and droplet impact velocity. As a continuation of the study of droplet impact contact characteristics by varying surface roughness and impact velocity, this study focuses on droplets impacting the vibrating surface at frequencies between 2-7 kHz. Atomization (water drop splitting into smaller droplets and ejecting from the surface) has been observed at different rates for all frequencies. The first set of data is collected by keeping roughness constant and increasing the amplitude of the vibration to observe the critical amplitude at which atomization is initiated. The surface roughness is varied for the second set of experiments. The data is quantified using image processing of the high-speed videos to obtain the rate of ejection for each case.
The Dynamics of Single and Double Cavitation Bubbles and Interaction Between Bubbles and Different Materials
Zhao, Ben (Virginia Tech, 2022-09-06)
We present two distinct projects in this article. In the first project, an experiment aiming to quantify the impacts of material acoustic impedance and thickness on single laser-induced cavitation bubble dynamics with measurements of exerted pressure on a specific material in order to identify the primary sources most responsible for material damages is presented in this article. Two types of major pressure sources have been identified. For bubble collapsing near a rigid wall, when standoff ratio γ < 0.6, the ring collapse is the most prominent pressure source. The jet takes the strongest effects at γ = 1.12. The pressure is minimal at γ = 0.913. After the first jet impingement, a second ring collapse will follow and input the maximum pressure to the wall. By further increasing γ, a similar pressure profile of the second collapse to the first collapse is achieved, during which the pressure for the second collapse is minimal at γ = 1.41 and for the jet is maximum at γ = 1.79. Compared with the maximum pressure dealt by the first jet, the second ring collapse and jet are increasing much faster with the bubble size and eventually overwhelm the first jet. However, the first ring collapse is still the most dominant pressure source responsible for material damages. For wall featuring smaller acoustic impedance or thickness that cannot be approximated to a rigid body, the ring collapse and jet occur at smaller standoff ratios. The cavity shrinking rate suggests the maximum pressure exerted on the wall at applicable standoff ratios should be smaller than that on a rigid wall. In the second project, a comprehensive collection of dynamics of one and two laser-induced cavitation bubbles collapsing near different boundaries is presented in this article by measuring the velocity fields using particle image velocimetry (PIV) techniques. Cases include a single bubble collapsing near the hard, medium, and soft walls characterized by acoustic impedance, free collapse of two bubbles, and two bubbles collapsing near the hard and soft walls. We implemented the most significant velocity and top velocity regions derived from each velocity field to analyze the features of these cases. Before converging to free collapse, the bubble near the hard wall experienced a significant velocity decrease before collapse, the bubble near the medium wall was severely damped at a specific standoff distance, and the bubble near the soft wall collapsed much earlier and preserved a linear velocity region at low speed. Free collapse of two same bubbles underwent a decrease of acceleration before collapse. Decreasing the size of one bubble caused a jet in the other. With the presence of a hard wall near two bubbles, the bubble closer to it may be stretched to a cavity with a high aspect ratio, leading to very mild collapse. With a bigger bubble between a smaller one and the soft wall, the merging cavity may suppress the tendency of jet formation, making the velocity stay at low levels throughout the lifetime. For configurations regarding single bubbles collapsing near a wall and free collapse of two same bubbles, we performed data scaling to study the velocity variations for different bubble sizes by controlling the standoff ratios and assessed the data quality aided by curving fitting and statistics. Results indicated measured velocity regarding a single bubble collapsing near the wall over its diameter remained the same given a standoff ratio, while measured velocity did not change given a standoff ratio for free collapse of two same bubbles within the scope of the experiment. In addition, we detailed the experimental setup and water treatment for better signal-to-noise ratios as well as validated the system from both the PIV and high speed imaging approaches using free collapse of a single bubble to ensure the reliability of this experiment.
Experimental and Numerical Investigation of Forward and Aft Swept Stepped Planing Hulls in Calm Water and Regular Waves
Husser, Nicholas Alexander (Virginia Tech, 2023-02-22)
Stepped hull forms are hulls with a vertical step in the hull bottom to improve performance at top speed. Stepped hulls are well documented anecdotally and scientifically to improve calm water performance at high speeds, but commonly demonstrate dangerous and unexpected dynamic instabilities during initial trials. These hulls also operate practically in waves, but their performance characteristics in waves are not well understood and rarely evaluated prior to full scale trials. To expand the scientific understanding of stepped hull performance, a systematic set of experiments and Reynolds Averaged Navier Stokes (RANS) computational fluid dynamics (CFD) simulations are used to evaluate the calm water performance, dynamic stability, and regular head wave response of two stepped hull models. Calm water experiments on two stepped hull models at varying displacement, longitudinal center of gravity location and forward speed offer data which can be used in the design to interpolate hull performance throughout expected operating conditions. CFD simulations in calm water are validated using the experimental results and numerical modeling approaches for stepped hull simulations are recommended. The calm water dynamic stability of both stepped hulls is investigated experimentally and numerically and procedures to evaluate the dynamic stability using both approaches are recommended. The performance of both stepped hulls in regular head waves is studied through experiments, which are used to validate CFD simulations of the hull in regular waves. System identification is used on five calm water CFD simulations to identify a reduced order model for the prediction of stepped hull response in waves.
An Experimental Investigation into the Passive Reconfiguration of Flexible Plates Near a Free Surface
Scianna, Nicholas Alexander (Virginia Tech, 2022-05-26)
Reconfiguration refers to the ability of a flexible structure to change its shape, allowing it to reduce its area perpendicular to the flow, to reduce drag. Decreasing the flexural rigidity of human-made structures can lead to improved designs that operate at higher propulsive efficiencies. The work presented in this thesis examines the physics surrounding a flexible plate under prescribed oscillatory heaving motions. White light movies were recorded at constant frequency and varying proximity to the free surface to investigate the change in reconfiguration as the plate approaches the free surface. Results, analyzed in terms of deformed plate shape, deflection, and plate tip kinematics, found that free surface effects increase the deflection of the plate as the plate approaches the free surface. Expanding on the initial experiments, a variety of frequencies were tested. The results show that each heaving frequency has a different critical height to the free surface in which deep water behavior is distinguished from shallow water behavior. At the critical depth, the plate deflection becomes asymmetric due to free surface effects. The second stage of experiments focused on measuring the fluid loading and fluid flow surrounding the flexible plate. The fluid loading, or drag force, acting on the plate was estimated by using a strain gauge load cell. Results of these experiments found that the drag force is equivalent on plates with lower heaving frequencies when compared to the highest heaving frequency tested due to increased reconfiguration at the higher frequency. The fluid moved from the keel to the edge of the plate as seen in the particle image velocimetry experiments. Higher heaving frequencies created faster fluid flow off the plate and stronger tip vortices being shed from the plate. When the flexible plate operated at large distances from the free surface, the fluid dynamics showed the same behavior for the upstroke and downstroke of the plate. Whereas, when the plate operated close to the free surface, a vortex only forms on the upstroke, leading to asymmetric loading and deformations.
An Experimental Investigation on the Performance of a Shape Changing, Bio-inspired F2MC Panel
Johansson, Oscar (Virginia Tech, 2024-05-23)
The purpose of this thesis is to explore the performance of a bio-inspired plate undergoing oscillatory heave motions and active shape change. The shape change will be achieved using a panel embedded with Fluidic Flexible Matrix Composite (F2MC) tubes for actuation. A beam, or plate strip, model is presented as a means of verifying that F2MC tubes can effectively serve as a means of actuation. This model was actuated in air and water at several internal tube pressures. The static experimental deflections were compared to two beam models relying on Euler-Bernoulli and Timoshenko beam theories with concentrated tip moments and a distributed moment. It was found that the Euler-Bernoulli model with a concentrated tip moment best approximated the static experimental deflections. Following the success of the plate strip, and panel with 10 embedded F2MC tubes was manufactured. The plate panel was constructed with Dragon Skin Silicone and embedded with two rows of five F2MC tubes which provide the means of shape actuation. Experimental results from actuating the panel in static conditions showed that F2MC tubes are an effective means of prescribing a repeatable shape change to a silicone panel. Then, Classical Plate Theory and First-Order Shear Deformation Plate Theory were used with a concentrated tip moment at the free edge to provide a means of modeling the full panel. When comparing the static experimental results to the numerical models, it was found that the deflected plate shape could be most accurately predicted at lower pressures for upward deflection and higher pressures for downward deflections. When tested in unsteady conditions in a heaving experiment (0.5 Hz to 2.3 Hz), the force measured at frequencies above 1.5 Hz were up to 3.6 times greater than those measured for frequencies below 1.5 Hz. Additionally, the phase difference between the tip deflection and force with respect to the keel position decreased for force as frequency increased, while the opposite was true for the tip deflection. At 1.5 Hz, the tip deflection and force were equally out of phase with the keel. When the panel was subjected to an oscillatory heaving motion while asymmetrically actuated, it was found that faster heaving frequencies resulted in higher maximum force values for all actuation pressures, actuation directions, and depths below the free surface. However, when subjected to dual actuation by pressurizing the top and bottom tubes at the same pressure, the tip amplitude was highly dependent on specific combinations of heaving frequency, actuation pressure, and depth below the free surface. This indicates that the actuation pressure must be tuned to the depth and frequency of operation to obtain the desired tip amplitude for a given application. These findings further the knowledge of shape-changing F2MC panels operating near a free surface and lay a groundwork for developing flapping propulsors that mimic marine animals.
Experimental Study of Pressure and Velocity Fluctuations Induced by Cavitation in a Small Venturi Channel
Zhang, Linrong; Zhang, Guangjian; Ge, Mingming; Coutier-Delgosha, Olivier (MDPI, 2020-12-08)
The purpose of this paper is to investigate experimentally the influence of the cavitation extent on the pressure and velocity fluctuations in a small convergent–divergent channel. The mean cavity length is determined from high-speed photography images. The mean pressure and the intensity of the pressure fluctuations are obtained from the transient pressure signals recorded by two pressure transducers at the inlet and outlet of the test section. The statistical turbulence quantities are derived from the instantaneous velocity fields measured by the laser-induced fluorescent particle image velocimetry (PIV-LIF) technique. The experimental results show that the decrease of the cavitation number (the increase in the extent of cavitation) leads to a rise in the turbulent fluctuations in the wake region due to the impact of vapour clouds collapsing, while the presence of a vapour phase is found to reduce the streamwise and cross-stream velocity fluctuations in the attached cavity. It might be attributed to two mechanisms: the presence of a vapour phase modifies the vortex-stretching process, and the cavitation compressibility damps out the turbulent fluctuations. Similar effects of cavitation are also observed in the pressure fluctuations.
Experimental Study of Two-Phase Cavitating Flows and Data Analysis
Ge, Mingming (Virginia Tech, 2022-05-25)
Cavitation can be defined as the breakdown of a liquid (either static or in motion) medium under very low pressure. The hydrodynamic happened in high-speed flow, where local pressure in liquid falls under the saturating pressure thus the liquid vaporizes to form the cavity. During the evolution and collapsing of cavitation bubbles, extreme physical conditions like high-temperature, high-pressure, shock-wave, and high-speed micro-jets can be generated. Such a phenomenon shall be prevented in hydraulic or astronautical machinery due to the induced erosion and noise, while it can be utilized to intensify some treatment processes of chemical, food, and pharmaceutical industries, to shorten sterilization times and lower energy consumption. Advances in the understanding of the physical processes of cavitating flows are challenging, mainly due to the lack of quantitative experimental data on the two-phase structures and dynamics inside the opaque cavitation areas. This dissertation is aimed at finding out the physical mechanisms governing the cavitation instabilities and making contributions in controlling hydraulic cavitation for engineering applications. In this thesis, cavitation developed in various convergent-divergent (Venturi) channels was studied experimentally using the ultra-fast synchrotron X-ray imaging, LIF Particle Image Velocimetry, and high-speed photography techniques, to (1) investigate the internal structures and evolution of bubble dynamics in cavitating flows, with velocity information obtained for two phases; (2) measure the slip velocity between the liquid and the vapor to provide the validation data for the numerical cavitation models; (3) consider the thermodynamic effects of cavitation to establish the relation between the cavitation extent and the fluid temperature, then and optimize the cavitation working condition in water; (4) seek the coherent structures of the complicated high-turbulent cavitating flow to reduce its randomness using data-driven methods.
Hydrodynamic Study of Pisciform Locomotion with a Towed Biolocomotion Emulator
Nguyen, Khanh Quoc (Virginia Tech, 2021-06-04)
The ability of fish to deform their bodies in steady swimming action is gaining from robotic designers. While bound by the same physical laws, fish have evolved to move in ways that often outperform artificial systems in critical measures such as efficiency, agility, and stealth through thousands of years of natural selection. As we expand our presence in the ocean with deep-sea exploration or offshore drilling for petroleum and natural gas, the demand for prolonging underwater operations is growing significantly. Therefore, it is critical for robotic designers to understand the physics of pisciform (fish-like) locomotion and learn how to effectively implement the propulsive mechanisms into their designs to create the next generation of aquatic robots. Aiming to assist this process, this thesis presents an experimental apparatus called Towed Biolocomotion Emulator (TBE), which is capable of imitating the undulating action of different fish species in steady swimming and can be quickly adapted to different configurations with modular modules. Using the TBE device, an experiment is performed to test its hydrodynamic performance and evaluate the effectiveness of the bio-inspired locomotion implemented on such a mechanical system. The analysis of hydrodynamic data collected from the experiment shows that there exists a small range of kinematic parameters where the undulating motion of the device produces the optimal performance. This result confirms the benefits of utilizing pisciform locomotion for small-scale underwater vehicles. In addition, this thesis also proposes a reduced-order flow model using the unsteady vortex lattice method (UVLM) to predict the hydrodynamic performance of such a system. The proposed model is then validated with the experimental data collected earlier. The tool developed can be employed to quickly explore the possible design space early in the conceptual design stage for such a bio-mimetic vehicle.
Impact Dynamics of Water Droplet on Solid Surfaces: Effect of Impact Reynolds Number, Hydrophobicity, Surface Roughness and Temperature
Naveed, Ahsan (Virginia Tech, 2023-06-23)
One of the most complicated issues the aerospace and aviation industries are dealing with is aircraft icing. The impact and freezing process of a water droplet on a cold surface has been investigated over time in order to develop preventative methods for avoiding icing. In the present study, we examined the behavior of a water droplet impacting on an aluminum plate with a surface roughness of 0.01µm and surface temperature variation from room temperature to 0oC, −5oC, −10oC and −15oC. The effect of droplet impact Reynolds number along with surface temperature variation on non-dimensional parameters like spread factor, retraction rate, and spread velocity is analyzed. The increase in impact Reynolds number and droplet spread factor is observed with a rise in the initial height of the droplet. At a higher Reynolds number, inertial forces are dominant over viscous and capillary forces, while at a lower Reynolds number, surface temperature shows a significant effect. The graphical representation of droplet retraction rate indicates a decrease with lower surface temperature and a rise with higher Reynolds numbers. Moreover, the spread velocity of the droplet is higher with an increased Reynolds number, and surface temperature does not have a notable effect on it. A rapid transition of momentum from vertical to horizontal direction occurs, and droplet dissipates energy in overcoming the viscous effects. The effect of surface roughness variation coupled with surface temperature is investigated in detail for three different surface roughness of aluminum and glass. The increase in surface roughness and temperature enhance hydrophobic behavior by repelling the droplet, while reduced surface temperatures show hydrophilic behavior by causing adhesion of the droplet on surface.
Left Ventricle Assist Device
Patel, Rudra Dharmeshbhai (Virginia Tech, 2023-07-10)
This research paper highlights the significance of mock circulatory loops in aiding the development and testing of left ventricle assist devices (LVADs). Heart disease is the leading cause of death in America, and LVAD offers an effective solution to patients with severe coronary artery disease who are not eligible for heart transplants. However, before testing LVADs in-vivo, extensive performance and reliability testing is required as per Food and Drug Administration guidelines. Mock circulatory loops are useful for simulating the cardiac cycle and capturing pressure and flow meter readings. The research focused on developing and testing a mock circulatory loop that accurately captures pressure and flow meter readings. The team experimented with various silicone elastomers and a urethane-based material, ClearFlex 30, to create a transparent phantom. The report highlights the importance of surface roughness for optical clarity, with an average surface roughness of 0.186 μ being ideal for optimal clarity achieved with 800-grit sandpaper. The impact of pressure differences between the aortic and mitral inlet/outlets was also studied, and it was found that the loop's resistance can be modified to achieve elevated pressure in the aortic outlet. The report further emphasized the importance of refractive index matching to perform particle image velocity. Matching the refractive index of the phantom with the medium is critical to avoid distortion and refraction of the light. Glycerol water was found to be an effective medium for refractive index matching. Overall, the successful testing of the mock circulatory loop provides valuable insights into the flow structure within the heart, aiding in the development of future LVADs. This research is significant in advancing the cardiology field and will potentially benefit many patients suffering from heart disease.
Long-Pulsed Laser-Induced Cavitation: Laser-Fluid Coupling, Phase Transition, and Bubble Dynamics
Zhao, Xuning (Virginia Tech, 2024-02-29)
This dissertation develops a computational method for simulating laser-induced cavitation and investigates the mechanism behind the formation of non-spherical bubbles induced by long-pulsed lasers. The proposed computational method accounts for the laser emission and absorption, phase transition, and the dynamics and thermodynamics of a two-phase fluid flow. In this new method, the model combines the Navier-Stokes (NS) equations for a compressible inviscid two-phase fluid flow, a new laser radiation equation, and a novel local thermodynamic model of phase transition. The Navier-Stokes equations are solved using the FInite Volume method with Exact two-phase Riemann solvers (FIVER). Following this method, numerical fluxes across phase boundaries are computed by constructing and solving one-dimensional bi-material Riemann problems. The new laser radiation equation is derived by customizing the radiative transfer equation (RTE) using the special properties of laser, including monochromaticity, directionality, high intensity, and a measurable focusing or diverging angle. An embedded boundary finite volume method is developed to solve the laser radiation equation on the same mesh created for the NS equations. The fluid mesh usually does not resolve the boundary and propagation directions of the laser beam, leading to the challenges of imposing the boundary conditions on the laser domain. To overcome this challenge, ghost nodes outside the laser domain are populated by mirroring and interpolation techniques. The existence and uniqueness of the solution are proved for the two-dimensional case, leveraging the special geometry of the laser domain. The method is up to second-order accuracy, which is also proved, and verified using numerical tests. A method of latent heat reservoir is developed to predict the onset of vaporization, which accounts for the accumulation and release of latent heat. In this work, the localized level set method is employed to track the bubble surface. Furthermore, the continuation of phase transition is possible in laser-induced cavitation problems, especially for long-pulsed lasers. A method of local correction and reinitialization is developed to account for continuous phase transitions. Several numerical tests are presented to verify the convergence of these methods. This multiphase laser-fluid coupled computational model is employed to simulate the formation and expansion of bubbles with different shapes induced by different long-pulsed lasers. The simulation results show that the computational method can capture the key phenomena in the laser-induced cavitation problems, including non-spherical bubble expansion, shock waves, and the ``Moses effect''. Additionally, the observed complex non-spherical shapes of vapor bubbles generated by long-pulsed laser reflect some characteristics (e.g., direction, width) of the laser beam. The dissertation also investigates the relation between bubble shapes and laser parameters and explores the transition between two commonly observed shapes -- namely, a rounded pear-like shape and an elongated conical shape -- using the proposed computational model. Two laboratory experiments are simulated, in which Holmium:YAG and Thulium fiber lasers are used respectively to generate bubbles of different shapes. In both cases, the predicted bubble nucleation and morphology agree reasonably well with the experimental observation. The full-field results of laser radiance, temperature, velocity, and pressure are analyzed to explain bubble dynamics and energy transmission. It is found that due to the lasting energy input, the vapor bubble's dynamics is driven not only by advection, but also by the continued vaporization at its surface. Vaporization lasts less than 1 microsecond in the case of the pear-shaped bubble, compared to over 50 microseconds for the elongated bubble. It is thus hypothesized that the bubble's morphology is determined by a competition between the speed of bubble growth due to advection and continuous vaporization. When the speed of advection is higher than that of vaporization, the bubble tends to grow spherically. Otherwise, it elongates along the laser beam direction. To test this hypothesis, the two speeds are defined analytically using a model problem and then estimated for the experiments using simulation results. The results support the hypothesis and also suggest that when the laser's power is fixed, a higher laser absorption coefficient and a narrower beam facilitate bubble elongation.
Numerical Investigation on Shape Impact of Deformable Droplets on Evaporation and Combustion: Method Development and Characterization
Setiya, Meha (Virginia Tech, 2023-08-21)
Inspired by the dilute spray regime in spray combustion, this dissertation explores the evaporation and combustion of an isolated droplet. Under a highly convective environment inside a gas combustor, due to imbalance of inertial and surface tension forces, the droplets of larger size in sprays exhibit notable deformations from spherical to non-spherical shapes. Such shape changes are generally observed but not quantified in experimental studies. Therefore, the effect of this deformation on droplet combustion dynamics is unknown yet. To bridge this gap, a comprehensive investigation of an isolated freely deforming droplet can be insightful as it can reveal more about the interaction of droplet shape with its evaporation and combustion. This work attempts to analyze and quantify the impact of such deformations on evaporation and combustion using interface-capturing Direct Numerical Simulation approach. With the focus on small-scale processes involved in evaporation as it is a pre-step for combustion, this dissertation first covers a thorough examination on evaporation of a deformable droplet under both natural and forced convection. A single component jet-fuel surrogate n-decane is chosen. To ensure that the droplet remains stationary throughout its lifetime, a novel numerical method called "gravity update method" is developed and implemented. The results obtained from these two separate studies are validated against experimental results and analytical correlations respectively. The findings from the investigation of droplet evaporation under forced convective flow at moderate Reynolds numbers are noteworthy. The droplet shape under such flow conditions is governed by Weber number (We) which is a ratio of inertial force to surface tension force. The results demonstrated upto 20% en- hancement in total evaporation rate for highly deformed droplets. This improvement is a net results of increased droplet surface area and alteration in the distribution of local evaporation flux ( m'' ). It is found that m'' is proportional to its curvature up to the point of flow separation which agrees with low Re theories on droplet evaporation by Tonini and Cossalli (International Journal of Heat and Mass Transfer 2013), Palmore (Journal of Heat Transfer 2022). Beyond the flow separation point, evaporation flux distribution depends on the boundary layer development and flow evolution downstream of the droplet. For highly deformed droplets, a larger wake region creates favorable fuel vapor gradients and promotes mixing in droplet wake, hence higher evaporation flux. Such positive impact of droplet deformation on total evaporation rate motivated further investigation on droplet combustion under a low Reynolds number convective flow. High pressure and temperature gas flow leads to Damköhler number is higher than 1. This fa- vors the generation of envelope type flame. The results show overall little sensitivity to combustion related parameters despite the droplet shape change and significant (upto 9%) enhancement in total evaporation rate. It is also noted that while burning, droplets do not reach critical deformation conditions and break-up even beyond the critical Weber number, suggesting the suppression of deformation due to faster evaporation rate. The findings presented in these studies provide substantial evidence for the interaction between droplet shape and flow dynamics. Therefore, it demonstrates the potential for enhancing the existing numerical models and analytical correlations by accounting the influence of droplet shape.
OpenFOAM Implementation of Microbubble Models for Ocean Applications
Harris, David Benjamin (Virginia Tech, 2021-07-27)
An investigation was carried out on the current state of the art in bubble modelling for computational fluid dynamics, and comparisons made between the different methods for both polydisperse and monodisperse multiphase flows. A multigroup method for polydisperse bubbly flows with the bubbles binned in terms of mass was selected from the various alternatives, which included other multigroup models and moment methods. The latter of these involve the integration of moments of the bubble number density function and transport of these quantities. The equations from this multigroup solver were then changed to more accurately and efficiently model cases involving extremely small bubbles over significant amounts of time, as the original model which was subsequently adapted had, as its primary purpose, simulation of larger bubbles over shorter periods of time. This was done by decoupling the gas and liquid momentum equations and adding an empirical rise velocity term for the bubbles. This new model was then partially implemented into OpenFOAM. The functioning of this new solver was confirmed by comparisons between the results and basic analytical solutions to the problems, as well as by means of comparison with another similar multiphase CFD solver (pbeTransportFoam). Following this confirmation of its functionality, the bubble model was implemented into another solver specifically designed for modelling wakes. Finally, the newly created solver was used to run some cases of interest involving a submerged wake.
Optimization of the critical overlap length of hudraulic distributors regarding cavitation inception
Maiga, Mahamadou Adama; Coutier-Delgosha, Olivier; Buisine, Daniel (2021-10)
The present paper is focused on the optimization of the overlap length in various types of distributors in hydraulic systems. It has been shown that this length has a strong influence on cavitation in rotating distributors, so its optimization would be beneficial to reduce the detrimental effects of cavitation, such as noise, vibrations, and instabilities. The objective of the present study is to generalize the modelling of the overlap length to different types of distributors and to determine the critical length at cavitation inception. The results show that the pressure losses in rotating distributors and in other types of distributors are similar, including the slide-valves that are the vast majority of distributors in hydraulic systems. It is confirmed that the overlap length has a stronger influence on the cavitation development than the distributor rotation speed. A critical overlap length of 0.071 mm is obtained, below which no significant development of cavitation is obtained, at all rotation speeds.
Pressure Shielding Mechanisms in Bio-Inspired Unidirectional Canopy Surface Treatments
Nurani Hari, Nandita (Virginia Tech, 2022-06-27)
Reduction of surface pressure fluctuations is desirable in various aerodynamic and hydrodynamic applications. Over the past few years, studies on canopy surface treatments have been conducted to investigate the fundamental mechanisms of surface pressure attenuation termed as pressure shielding. This work talks about the design, development and experimental testing of unidirectional canopy surface treatments which are evenly spaced arrays of streamwise rods placed parallel to the wall without an entrance condition. The canopy designs are based on surface treatments tested by Clark et al. (2014) inspired by the downy coating on owl wings. The main objective of the work is to establish fundamental physical and mathematical basis for treatments that shield aerodynamic surfaces from turbulent pressure fluctuations, while maintaining the wall-normal transport of momentum and low aerodynamic drag. Experimental testing of these canopy treatments are performed in the Anechoic Wall-Jet facility at Virginia Tech. Different canopy configurations are designed to understand the effect of various geometric parameters on the surface pressure attenuation. The treatment is found to exhibit broadband reduction in the surface pressure spectrum. Attenuation develops in two frequency regions which scale differently depending on two different mechanisms. Canopies seems to reduce the large-scale turbulent fluctuations up to nearly twice the height. Semi-analytical model is developed to predict surface pressure spectra in a wall-jet and canopy flow. The rapid term model shows that the inflection in the streamwise mean velocity profile is the most dominant source of surface pressure fluctuations. Synchronized pressure and velocity measurements elucidate significant features of the sources that could be affecting surface pressure fluctuations. Overall, this study explores the qualitative and quantitative physics behind pressure shielding mechanism which find application particularly in trailing edge noise reduction.
A qualitative assessment and optimization of URANS modelling for unsteady cavitating flows
Apte, Dhruv Girish (Virginia Tech, 2024-06-07)
Cavitation is characterized by the formation of vapor bubbles when the pressure in a working fluid drops sharply below the vapor pressure. These bubbles, upon exiting the low-pressure region burst emanating tremendous amounts of energy. Unsteady cavitating flows have been influential in several aspects from being responsible for erosion damage and vibrations in hydraulic engineering devices to being used for non-invasive medical surgeries and drilling for geothermal energy. While the phenomenon has been investigated using both experimental and numerical methods, it continues to pose a challenge for numerical modelling techniques due to its flow unsteadiness and the cavitation-turbulence interaction. One of the principal aspects to modelling cavitation requires the coupling of a cavitation and a turbulence model. While, scale-resolving turbulence modelling techniques like Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) upto a certain extent may seem an intuitive solution, the physical complexities involved with cavitation result in extremely high computational costs. Thus, Unsteady Reynolds-Averaged Navier-Stokes (URANS) models have been widely utilized as a workhorse for cavitating simulations. However, URANS models are unable to reproduce the periodic vapor shedding observed in experiments and thus, are often corrected by empirical correction. Recently, some models termed as hybrid RANS-LES models that behave as RANS or LES depending on location of flow have been introduced and employed to model cavitating flows. In addition, there has also been a rise in defining some frameworks that use data from high-fidelity simulations or experiments to drive numerical algorithms and aid standard turbulence modelling procedures for accurately simulating turbulent flows. This dissertation is aimed at (1) evaluating the abilities of these corrections, traditional URANS and hybrid RANS-LES models to model cavitation and (2) optimizing the URANS modelling strategy by designing a methodology driven by experimental data to augment the turbulence modelling to simulate cavitating flow in a converging-diverging nozzle.

Browsing by Author "Coutier-Delgosha, Olivier"

Results Per Page

Sort Options