Earlier investigations of Combustion Powered Actuation (COMPACT) have demonstrated its utility for high-speed aerodynamic flow control. In this actuation approach, momentary (pulsed) actuation jets are produced by the ignition of a mixture of gaseous fuel and oxidizer within a cubic-centimeter scale chamber. The combustion process yields a high pressure burst and the ejection of a high-speed exhaust jet. The present thesis focuses on characterization of the effects of the internal flow (which is altered through the fuel and oxidizer inlet streams) on mixing and flame propagation within the actuator's combustion chamber, and thereby on actuator operation and performance. A test chamber with a grid of interchangeable air and fuel inlets was used for parametric investigations of the effects of inlet size and location. Actuator performance is characterized using dynamic pressure measurements and phase-locked Particle Image Velocimetry (PIV) of the combustor's internal flow field in the presence and absence of the active combustion process. Over the range tested, increased momentum of the air inlet jet for a given flow rate improves the actuator performance by increasing bulk velocities and small-scale motions within the chamber, thus yielding net higher flame propagation speed and subsequently faster pressure rise and higher pressure peak. Variation in inlet location that results in swirling flow within the chamber yields higher internal pressures while air flow over the spark ignition site yields lower internal pressures and erratic combustion. Improved refill and combustion processes will lead to enhanced performance combustor designs.

This study was concerned with the most common reverse flow type of cyclones where the flow enters the cyclone through a tangential inlet and leaves via an axial outlet pipe at the top of the cyclone. Numerical computations of two different cyclones were based on the so-called Stairmand cyclone. The difference in geometry between these two cyclones was basically characterized by the geometrical swirl number Sg of 3.5 and 4. Turbulent secondary flows inside a straight square channel have been studied numerically by using Large Eddy Simulation (LES) in order to verify the implementation process. Prandtl’s secondary motion calculated by LES shows satisfying agreement with both, Direct Numerical Simulation (DNS) and experimental results. Numerical calculations were carried out at various axial positions and at the apex cone of a gas cyclone separator. Two different NS-solvers (a commercial one, and a research code), based on a pressure correction algorithm of the SIMPLE method have been applied to predict the flow behavior. The flow was assumed as unsteady, incompressible and isothermal. A k − epsilon turbulence model has been applied first using the commercial code to investigate the gas flow. Due to the nature of cyclone flows, which exhibit highly curved streamlines and anisotropic turbulence, advanced turbulence models such as RSM (Reynolds Stress Model) and LES (Large Eddy Simulation) have been used as well. The RSM simulation was performed using the commercial package CFX4.4, while for the LES calculations the research code MISTRAL/PartFlow-3D code developed in our multiphase research group has been applied utilizing the Smagorinsky model. It was found that the k − epsilon model cannot predict flow phenomena inside the cyclone properly due to the strong curvature of the streamlines. The RSM results are comparable with LES results in the area of the apex cone plane. However, the application of the LES reveals qualitative agreement with the experimental data, but requires higher computer capacity and longer running times than RSM. These calculations of the continuous phase flow were the basis for modeling the behavior of the solid particles in the cyclone separator. Particle trajectories, pressure drop and the cyclone separation efficiency have been studied in some detail. This thesis is organized into five chapters. After an introduction and overview, chapter 2 deals with continuous phase flow turbulence modeling including the governing equations. The emphasis will be based on LES modelling. Furthermore, the disperse phase motion is treated in chapter 3. In chapter 4, the validation process of LES implementation with channel flow is presented. Moreover, prediction profiles of the gas flow are presented and discussed. In addition, disperse phase flow results are presented and discussed such as particle trajectories; pressure drop and cyclone separation efficiency are also discussed. Chapter 5 summarizes and concludes the thesis.

Zugl.: Berlin, Techn. Univ., Diss., 2009

A number of clinical, in vitro and computational studies have shown the potential for thromboembolic complications in bileaflet mechanical heart valves (BMHV), primarily due to the complex and unsteady flows in the valve hinges. These studies have focused on quantitative and qualitative parameters such as velocity magnitude, turbulent shear stresses, vortex formation and platelet activation to identify potential for blood damage. However, experimental characterization of the whole flow fields within the valve hinges has not yet been conducted. This information can be utilized to investigate instantaneous damage to blood elements and also to validate numerical studies focusing on the hinge’s complex fluid dynamics. The objective of this study was therefore to develop a high-resolution imaging system to characterize the flow fields and global velocity maps in a BMHV hinge. Subsequently, the present study investigated the effect of hinge gap width on flow fields in a St. Jude Medical BMHV. The results from this study suggest that the BMHV hinge design is a delicate balance between reduction of fluid shear stresses and areas of flow stasis during leakage flow, and needs to be optimized to ensure minimal thromboembolic complications. Overall, the current study demonstrates the ability of high-resolution Micro Particle Image Velocimetry to assess the fluid flow fields within the hinges of bileaflet mechanical heart valves, which can be extended to investigate micro-scale flow domains in critical regions of other cardiovascular devices to assess their blood damage potential.

The stability and turbulence characteristics of a vortex filament emanating from a single-bladed rotor in hover are investigated using proper orthogonal decomposition. The rotor is operated at a tip chord Reynolds number and a tip Mach number of 218,000 and 0.22, respectively, and with a blade loading of CT /σ = 0.066. In-plane components of the velocity field (normal to the axis of the vortex filament) are captured by way of 2D particle image velocimetry with corrections for vortex wander being performed using the Γ1 method. Using the classical form of POD, the first POD mode alone is found to encompass nearly 75% of the energy for all vortex ages studied and is determined using a grid of sufficient resolution as to avoid numerical integration errors in the decomposition. The findings reveal an equal balance between the axisymmetric and helical modes during vortex roll-up which immediately transitions to helical mode dominance at all other vortex ages. This helical mode is one of the modes of the elliptic instability. While the snapshot POD is shown to reveal similar features of the first few energetic modes, the classical POD is employed here owing to the easier interpretation of the Fourier-azimuthal modes. The spatial eigenfunctions of the first few Fourier-azimuthal modes associated with the most energetic POD mode are shown to be sensitive to the choice of the wander correction technique used. Higher Fourier-azimuthal modes are observed in the outer portions of the vortex and appeared not to be affected by the choice of the wander correction technique used. / text

Simulations are becoming increasingly popular in science and engineering. One type of simulation is Computation Fluid Dynamics (CFD) that is used when closed forms solutions are impractical. The field of Verification & Validation emerged from the need to assess simulation accuracy as they often contain approximations and calibrations. Validation involves the comparison of experimental data with simulation outputs and is the focus of this work. Errors in simulation predictions may be assessed in this way. Validation requires highly-detailed data and description to accompany these data, and uncertainties are very important. The purpose of this work is to provide highly complete validation data to assess the accuracy of CFD simulations. This aim is fundamentally different from the typical discovery experiments common in research. The measurement of these physics was not necessarily original but performed with modern, high fidelity methods. Data were tabulated through an online database for direct use in Reynolds-Averaged Navier Stokes simulations. Detailed instrumentation and documentation were used to make the data more useful for validation. This work fills the validation data gap for steady and transient mixed convection. The physics in this study included mixed convection on a vertical flat plate. Mixed convection is a condition where both forced and natural convection influence fluid momentum and heat transfer phenomena. Flow was forced over a vertical flat plate in a facility built for validation experiments. Thermal and velocity data were acquired for steady and transient flow conditions. The steady case included both buoyancy-aided and buoyancy-opposed mixed convection while the transient case was for buoyancy-opposed flow. The transient was a ramp-down flow transient, and results were ensemble-averaged for improved statistics. Uncertainty quantification was performed on all results with bias and random sources. An independent method of measuring heat flux was devised to assess the accuracy of commercial heat flux sensors used in the heated wall. It measured the convective heat flux by the temperature gradient in air very near the plate surface. Its accuracy was assessed by error estimations and uncertainty quantification.

No description available.

The self-similar evolution of a turbulent Rayleigh-Taylor (R-T) mix is investigated through experiments and numerical simulations. The experiments consisted of velocity and density measurements using thermocouples and Particle Image Velocimetry techniques. A novel experimental technique, termed PIV-S, to simultaneously measure both velocity and density fields was developed. These measurements provided data for turbulent correlations, power spectra, and energy balance analyses. The self-similarity of the flow is demonstrated through velocity profiles that collapse when normalized by an appropriate similarity variable and power spectra that evolve in a shape-preserving form. In the self-similar regime, vertical r.m.s. velocities dominate over the horizontal r.m.s. velocities with a ratio of 2:1. This anisotropy, also observed in the velocity spectra, extends to the Taylor scales. Buoyancy forcing does not alter the structure of the density spectra, which are seen to have an inertial range with a -5/3 slope. A scaling analysis was performed to explain this behavior. Centerline velocity fluctuations drive the growth of the flow, and can hence be used to deduce the growth constant. The question of universality of this flow was addressed through 3D numerical simulations with carefully designed initial conditions. With long wavelengths present in the initial conditions, the growth constant was found to depend logarithmically on the initial amplitudes. In the opposite limit, where long wavelengths are generated purely by the nonlinear interaction of shorter wavelengths, the growth constant assumed a universal lower bound value of

The tendency of aortic bileaflet mechanical heart valves (BiMHVs) to promote thrombosis has been well documented in the literature. The relationship of thrombosis to valve fluid dynamics has prompted numerous studies of aortic BiMHV flow. In this study, steady flow was investigated downstream of a model Carbomedics No. 25 BiMHV in an axisymmetric aortic sinus using stereoscopic particle image velocimetry (SPIV). The Reynolds number based on inlet diameter was 7600, and the measurement plane was perpendicular to the leaflet axes at the centerline of the aortic sinus. The typical formation of three jets was observed: the upper and lower lateral orifice jets, and the central jet. Flow separation from the valve ring was seen, and large scale vortices were identified in both the upper and lower sinus regions. An asymmetry in the reverse flow was found, and possible causes were discussed. All three jets were seen to decay similarly to free rectangular jets, with zero decay initially, followed by a 'linear' decay rate in which Umax^2~X. The central jet was also seen to be self similar in the linear decay region. Analysis of the out-of-plane velocity yielded two alternate explanations of streamwise vortex (i.e. Wx) structure, with either a four-cell or an eight-cell streamwise vortex structure being present in the mean velocity field. Organization of large scale three dimensional flow structures was thus apparent. Calculation of in-plane Reynolds stresses showed that values were highest in the outer shear layers of the lateral orifice jets. Elevated Reynolds shear stress values were also found in the leaflet wake regions, and the shear layers of the central jet.

The tendency of aortic bileaflet mechanical heart valves (BiMHVs) to promote thrombosis has been well documented in the literature. The relationship of thrombosis to valve fluid dynamics has prompted numerous studies of aortic BiMHV flow. In this study, steady flow was investigated downstream of a model Carbomedics No. 25 BiMHV in an axisymmetric aortic sinus using stereoscopic particle image velocimetry (SPIV). The Reynolds number based on inlet diameter was 7600, and the measurement plane was perpendicular to the leaflet axes at the centerline of the aortic sinus. The typical formation of three jets was observed: the upper and lower lateral orifice jets, and the central jet. Flow separation from the valve ring was seen, and large scale vortices were identified in both the upper and lower sinus regions. An asymmetry in the reverse flow was found, and possible causes were discussed. All three jets were seen to decay similarly to free rectangular jets, with zero decay initially, followed by a 'linear' decay rate in which Umax^2~X. The central jet was also seen to be self similar in the linear decay region. Analysis of the out-of-plane velocity yielded two alternate explanations of streamwise vortex (i.e. Wx) structure, with either a four-cell or an eight-cell streamwise vortex structure being present in the mean velocity field. Organization of large scale three dimensional flow structures was thus apparent. Calculation of in-plane Reynolds stresses showed that values were highest in the outer shear layers of the lateral orifice jets. Elevated Reynolds shear stress values were also found in the leaflet wake regions, and the shear layers of the central jet.