No description available.

No description available.

No description available.

No description available.

Fuel cells are promising technology to meet the energy need of the future. This alternative energy source is clean and efficient, and with the continuous decrease in fossil fuel resources, one of the best bets towards sustaining our power needs. Fuel cells are being used in automobiles as well as to fulfill portable power needs. In this work a computational model has been developed for fuel cells which can be used to simulate traditional as well as passive proton exchange membrane fuel cell behavior. The model is unsteady, two phase, nonisothermal in nature, and also capable of handling natural convection or buoyancy driven flows. The model also takes into account electrochemical reactions at catalyst sites. The model has been implemented and validated against experiments. It is used to carry out unsteady simulations to study start-up characteristic of proton exchange membrane fuel cells and to follow the behavior of liquid water as well as heat transfer within the cell. The buoyancy model is used to simulate a natural convection region and a passive fuel cell (used for portable applications). Design of passive fuel cells is driven by high temperature regimes and that issue has been further explored.

Unsteady flow fields over a NACA 0012 at an angle of attack are investigated. The first is the classical pitching motion about the airfoil's quarter chord. The second is the flow over a fixed airfoil immersed in the wake of the pitching airfoil. Large reduced frequencies are considered. Measurements were obtained in a water tunnel by Laser-Doppler velocimetry. Ensemble-averaged velocity measurements were obtained in the vicinity of the trailing edges of both the pitching and the fixed airfoils. The flowfields in the wake and at the trailing edges of both airfoils were studied visually. The validity of the quasi-steady and an extension to an unsteady Kutta condition are examined. A new dynamic similarity parameter is proposed. An analytical method based on the dynamics of discrete vortices is employed. Numerical calculations of the flow over the fixed airfoil are compared with experimental results. / Ph. D.

Unsteady flowfields are integral to high-speed applications, demanding precise modelling to characterize their unsteady features accurately. The simulation of unsteady supersonic and hypersonic flows is inherently computationally expensive, requiring a highly refined mesh to capture these unsteady effects. While anisotropic metric-based adaptive mesh refinement has proven effective in achieving accuracy with much less complexity, current algorithms are primarily tailored for steady flow fields. This thesis presents a novel approach to address the challenges of anisotropic grid adaptation of unsteady flows by leveraging a data-driven technique called Dynamic Mode Decomposition (DMD). DMD has proven to be a powerful tool to model complex nonlinear flows, given its links to the Koopman operator, and also its easy mathematical implementation. This research proposes the integration of DMD into the process of anisotropic grid adaptation to dynamically adjust the mesh in response to evolving flow features. The effectiveness of the proposed approach is demonstrated through numerical experiments on representative unsteady flow configurations, such as a cylinder in a subsonic flow and a cylinder in a supersonic channel flow. Results indicate that the incorporation of DMD enables an accurate representation of unsteady flow dynamics. Overall, this thesis contributes to making advances in the adaptation of unsteady flows. The novel framework proposed makes it computationally tractable to track the evolution of the main coherent features of the flowfield without losing out on accuracy by using a data-driven method. / Master of Science / Simulating unsteady, high-speed fluid flows around objects like aircraft and rockets poses a significant computational challenge. These flows exhibit rapidly evolving, intricate pattern structures that demand highly refined computational meshes to capture accurately. However, using a statically refined mesh for the entire simulation is computationally prohibitive. This research proposes a novel data-driven approach to enable efficient anisotropic mesh adaptation for such unsteady flow simulations. It leverages a technique called Dynamic Mode Decomposition (DMD) to model the dominant coherent structures and their evolution from snapshot flow field data. DMD has shown powerful capabilities in identifying the most energetic flow features and their time dynamics from numerical or experimental data. By integrating DMD into the anisotropic mesh adaptation process, the computational mesh can be dynamically refined anisotropically just in regions containing critical time-varying flow structures. The efficacy of this DMD-driven anisotropic adaptation framework is demonstrated in representative test cases - an unsteady subsonic flow over a circular cylinder and a supersonic channel flow over a cylinder. Results indicate that it enables accurate tracking and resolution of the key unsteady flow phenomena like vortex shedding using far fewer computational cells compared to static mesh simulations. In summary, this work makes anisotropic mesh adaptation computationally tractable for unsteady flow simulations by leveraging data-driven DMD modelling of the evolving coherent structures. The developed techniques pave the way for more accurate yet efficient unsteady CFD simulations.

The unsteady behavior of a tip leakage flow downstream of a simulated axial compressor rotor has been studied. The Virginia Tech low speed linear cascade wind tunnel was adapted to model the unsteady tip leakage flow produced by a rotor operating in the vortical wakes of a set of stator vanes. The cascade, consisting of 8 GE rotor B blades, has adjustable tip gap, inlet angle of 65.1 degrees, turning angle of 11.8 degrees and solidity of 1.076. The cascade Reynolds number, based on blade chord, was 393,000. A moving end wall was used to simulate the relative motion between rotor and casing, and vortex generators attached to the moving end wall were used to produce an idealized periodic unsteady vortical inflow similar to that shed by the junction of a row of inlet guide vanes. Measurements of the vortical inflow to the cascade produced by the generators and of the mean blade loading at the mid span are presented. The periodic and aperiodic behavior of the tip leakage flow downstream of the cascade, produced by this vortical disturbance, is also presented using phase and time averaged 3-component turbulence and pressure fluctuation measurements. These measurements are made for tip gap from 0.83% to 3.3% chord and streamwise locations from 0.772% to 1.117% blade spacing axially downstream of the cascade. The phase averaged inflow measurements reveal that the inflow produced by the vortex generators consists of a pair asymmetric counter-rotating vortices embedded in a thin (4.6% chord) endwall boundary layer. The vortices extend some 7.4% chord from the end wall. Their strength is about two orders smaller than the typical circulation of the tip leakage vortices produced by the cascade. Phase averaged single point three component hot-wire measurements downstream of the cascade reveal that the vortical inflow is, however, capable of producing significant large scale fluctuations in the size, strength, structure and position of the tip leakage vortex. These effects increase in magnitude with increase of tip gap. For small tip gaps these effects appear to be due to simple superposition between the inflow vortices and the tip leakage vortex. However for larger tip gaps these effects appear primarily a consequence of the inflow vortices interfering with the shedding of circulation from the blade tip. The fact that the circulation fluctuation is consistent with the inviscid unsteady loading prediction suggests that the inviscid response may be a major mechanism for generating the tip leakage unsteadiness. Although there is large periodic fluctuation in the tip leakage flow disturbed by the inflow, there is a larger aperiodic component. Two point correlation measurements and linear stochastic estimation are used to reveal the structure of this aperiodic part for a tip gap of 3.3% chord. The aperiodic fluctuation, containing most of the turbulence energy, is found appearing to be organized structures in large scale, and making the estimated instantaneous velocity field significantly different from the phase averaged periodic velocity field. Phase averaged pressure fluctuation measurements made using a microphone in the tip leakage vortex downstream of the cascade reveal that there are significant periodic fluctuating pressure waves and intense mean square fluctuation of the aperiodic fluctuating pressure. They are consistent with the measured periodic flow and aperiodic flow field respectively. These microphone measurements are validated using fluctuating pressure gradient estimates determined from the hot-wire measurements. / Ph. D.

A series of experiments were conducted in a F109 turbofan engine to investigate the unsteady wake profiles of an Inlet Guide Vane (IGV) at a typical spacing to the downstream fan at subsonic and transonic relative blade velocities. The sharp trailing-edge vanes were designed to produce a wake profile consistent with modern IGV. Time averaged baseline measurements were first performed with the IGV located upstream of the aerodynamic influence of the fan. Unsteady experiments were performed with an IGV-fan spacing of 0.43 fan chords. High-frequency on-vane pressure measurements showed strong peak-to-peak amplitudes at the blade passing frequency (BPF) of 4.7 psi at the transonic fan speeds. High-frequency total pressure measurements of the IGV wake were taken between the IGV and fan. Results showed that the total pressure loss coefficient of the time averaged IGV wake is reduced by 30% for the subsonic fan, and increased by a factor of 2 for the transonic fan compared to the baseline. Time resolved wake profiles for subsonic fan speeds show constructive and destructive interactions over each blade pass generated by the fan potential flow field. Time resolved wake profiles for the transonic fan speeds show that shock interactions with the IGV surface result in the wake shedding off of the vane at the BPF. Furthermore, the effectiveness of trailing edge blowing (TEB) flow control was investigated. TEB is the method of injecting air aft of the IGV to reduce the low pressure regions (deficits) in the viscous wakes shed by the vanes. Minimizing the IGV wakes reduces the forcing function on the downstream fan blades, thereby reducing high cycle fatigue. The TE span of the vane contains discrete holes at the axial centerline for TEB. Baseline results showed that TEB eliminates the IGV wake, while using only 0.03% of the total engine mass flow per IGV. TEB for the subsonic fan at the close spacing shows complete wake filling using the same mass flow as the baseline. TEB for the transonic fan shows a reduction of 68% in the total pressure loss coefficient, while requiring 2.5 times the mass flow as the baseline. / Ph. D.

The interaction between rotor blade tip leakage vortex and inflow disturbances, such as encountered in shrouded marine propulsors, was simulated in the Virginia Tech Linear Cascade Wind Tunnel equipped with a moving endwall system. Upstream of the blade row, idealized periodic inflow unsteadiness was generated using vortex generator pairs attached to the endwall at the same spacing as the blade spacing. At three tip gap settings, 1.7%c, 3.3%c and 5.7%c, the flow near the lower endwall of the center blade passage was investigated through three-component mean velocity and turbulence distributions measured by four-sensor hotwires. Besides time-averaged data, the measurements were processed for phase-locked analysis, with respect to pitchwise locations of the vortex generators relative to the blade passage. Moreover, surface pressure distributions at the blade tip were acquired at eight tip gaps from 0.87%c to 12.9%c. Measurements of pressure-velocity correlation were also performed with wall motion but without inflow disturbances. Achieved in this study is an understanding of the characteristics and structures of the tip leakage vortex at its initial formation. The mechanism of the tip leakage vortex formation seems to be independent of the tip gap setting. The tip leakage vortex consists of a vortical structure and a region of low streamwise-momentum fluid next to the endwall. The vortical structure is initially attached to the blade tip that creates it. This structure picks up circulation shed from that blade tip, as well as those from the endwall boundary layer, and becomes stronger with downstream distance. Partially induced by the mirror images in the endwall, the vortical structure starts to move across the passage resulting in a reduction in its rotational strength as the cross sectional area of the vortex increases but little circulation is added. The larger the tip gap, the longer the vortical structure stays attached to the blade tip, and the stronger the structure when it reaches downstream of the passage. Phased-averaged data show that the inflow disturbances cause small-scale responses and large-scale responses upstream and downstream of the vortex shedding location, respectively. This difference in scale is possibly dictated by a variation in the shedding location since the amount of circulation in the vortex is dependent on this location. The inflow disturbances possibly cause a variation in the shedding location by manipulating the separation of the tip leakage flow from the endwall and consequently the flow's roll-up process. Even though this manipulation only perturbs the leakage flow in a small scale, the shedding mechanism of the tip leakage vortex amplifies the outcome. / Ph. D.