Matching Items (6)
Filtering by
- All Subjects: Computational Fluid Dynamics
- Creators: Calhoun, Ronald
Description
ABSTRACT A vortex tube is a device of a simple structure with no moving parts that can be used to separate a compressed gas into a hot stream and a cold stream. Many studies have been carried out to find the mechanisms of the energy separation in the vortex tube. Recent rapid development in computational fluid dynamics is providing a powerful tool to investigate the complex flow in the vortex tube. However various issues in these numerical simulations remain, such as choosing the most suitable turbulent model, as well as the lack of systematic comparative analysis. LES model for the vortex tube simulation is hardly used in the present literatures, and the influence of parameters on the performance of the vortex tube has scarcely been studied. This study is aimed to find the influence of various parameters on the performance of the vortex tube, the best geometric value of vortex tube and the realizable method to reach the required cold out flow rate 40 kg/s . First of all, setting up an original 3-D simulation vortex tube model. By comparing experiment results reported in the literature and our simulation results, a most suitable model for the simulation of the vortex tube is obtained. Secondly, we perform simulations to optimize parameters that can deliver a set of desired output, such as cold stream pressure, temperature and flow-rate. We also discuss the use of the cold air flow for petroleum engineering applications.
ContributorsCang, Ruijin (Author) / Chen, Kangping (Thesis advisor) / Huang, Hueiping (Committee member) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2013
Description
Modern day gas turbine designers face the problem of hot mainstream gas ingestion into rotor-stator disk cavities. To counter this ingestion, seals are installed on the rotor and stator disk rims and purge air, bled off from the compressor, is injected into the cavities. It is desirable to reduce the supply of purge air as this decreases the net power output as well as efficiency of the gas turbine. Since the purge air influences the disk cavity flow field and effectively the amount of ingestion, the aim of this work was to study the cavity velocity field experimentally using Particle Image Velocimetry (PIV). Experiments were carried out in a model single-stage axial flow turbine set-up that featured blades as well as vanes, with purge air supplied at the hub of the rotor-stator disk cavity. Along with the rotor and stator rim seals, an inner labyrinth seal was provided which split the disk cavity into a rim cavity and an inner cavity. First, static gage pressure distribution was measured to ensure that nominally steady flow conditions had been achieved. The PIV experiments were then performed to map the velocity field on the radial-tangential plane within the rim cavity at four axial locations. Instantaneous velocity maps obtained by PIV were analyzed sector-by-sector to understand the rim cavity flow field. It was observed that the tangential velocity dominated the cavity flow at low purge air flow rate, its dominance decreasing with increase in the purge air flow rate. Radially inboard of the rim cavity, negative radial velocity near the stator surface and positive radial velocity near the rotor surface indicated the presence of a recirculation region in the cavity whose radial extent increased with increase in the purge air flow rate. Qualitative flow streamline patterns are plotted within the rim cavity for different experimental conditions by combining the PIV map information with ingestion measurements within the cavity as reported in Thiagarajan (2013).
ContributorsPathak, Parag (Author) / Roy, Ramendra P (Thesis advisor) / Calhoun, Ronald (Committee member) / Lee, Taewoo (Committee member) / Arizona State University (Publisher)
Created2013
Description
The flow around a golf ball is studied using direct numerical simulation (DNS). An immersed boundary approach is adopted in which the incompressible Navier-Stokes equations are solved using a fractional step method on a structured, staggered grid in cylindrical coordinates. The boundary conditions on the surface are imposed using momentum forcing in the vicinity of the boundary. The flow solver is parallelized using a domain decomposition strategy and message passing interface (MPI), and exhibits linear scaling on as many as 500 processors. A laminar flow case is presented to verify the formal accuracy of the method. The immersed boundary approach is validated by comparison with computations of the flow over a smooth sphere. Simulations are performed at Reynolds numbers of 2.5 × 104 and 1.1 × 105 based on the diameter of the ball and the freestream speed and using grids comprised of more than 1.14 × 109 points. Flow visualizations reveal the location of separation, as well as the delay of complete detachment. Predictions of the aerodynamic forces at both Reynolds numbers are in reasonable agreement with measurements. Energy spectra of the velocity quantify the dominant frequencies of the flow near separation and in the wake. Time-averaged statistics reveal characteristic physical patterns in the flow as well as local trends within dimples. A mechanism of drag reduction due to the dimples is confirmed, and metrics for dimple optimization are proposed.
ContributorsSmith, Clinton E (Author) / Squires, Kyle D (Thesis advisor) / Balaras, Elias (Committee member) / Herrmann, Marcus (Committee member) / Adrian, Ronald (Committee member) / Stanzione, Daniel C (Committee member) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2011
Description
A moving overlapping mesh methodology that achieves spectral accuracy in space and up to second-order accuracy in time is developed for solution of unsteady incompressible flow equations in three-dimensional domains. The targeted applications are in aerospace and mechanical engineering domains and involve problems in turbomachinery, rotary aircrafts, wind turbines and others. The methodology is built within the dual-session communication framework initially developed for stationary overlapping meshes. The methodology employs semi-implicit spectral element discretization of equations in each subdomain and explicit treatment of subdomain interfaces with spectrally-accurate spatial interpolation and high-order accurate temporal extrapolation, and requires few, if any, iterations, yet maintains the global accuracy and stability of the underlying flow solver. Mesh movement is enabled through the Arbitrary Lagrangian-Eulerian formulation of the governing equations, which allows for prescription of arbitrary velocity values at discrete mesh points.
The stationary and moving overlapping mesh methodologies are thoroughly validated using two- and three-dimensional benchmark problems in laminar and turbulent flows. The spatial and temporal global convergence, for both methods, is documented and is in agreement with the nominal order of accuracy of the underlying solver.
Stationary overlapping mesh methodology was validated to assess the influence of long integration times and inflow-outflow global boundary conditions on the performance. In a turbulent benchmark of fully-developed turbulent pipe flow, the turbulent statistics are validated against the available data.
Moving overlapping mesh simulations are validated on the problems of two-dimensional oscillating cylinder and a three-dimensional rotating sphere. The aerodynamic forces acting on these moving rigid bodies are determined, and all results are compared with published data. Scaling tests, with both methodologies, show near linear strong scaling, even for moderately large processor counts.
The moving overlapping mesh methodology is utilized to investigate the effect of an upstream turbulent wake on a three-dimensional oscillating NACA0012 extruded airfoil. A direct numerical simulation (DNS) at Reynolds Number 44,000 is performed for steady inflow incident upon the airfoil oscillating between angle of attack 5.6 and 25 degrees with reduced frequency k=0.16. Results are contrasted with subsequent DNS of the same oscillating airfoil in a turbulent wake generated by a stationary upstream cylinder.
The stationary and moving overlapping mesh methodologies are thoroughly validated using two- and three-dimensional benchmark problems in laminar and turbulent flows. The spatial and temporal global convergence, for both methods, is documented and is in agreement with the nominal order of accuracy of the underlying solver.
Stationary overlapping mesh methodology was validated to assess the influence of long integration times and inflow-outflow global boundary conditions on the performance. In a turbulent benchmark of fully-developed turbulent pipe flow, the turbulent statistics are validated against the available data.
Moving overlapping mesh simulations are validated on the problems of two-dimensional oscillating cylinder and a three-dimensional rotating sphere. The aerodynamic forces acting on these moving rigid bodies are determined, and all results are compared with published data. Scaling tests, with both methodologies, show near linear strong scaling, even for moderately large processor counts.
The moving overlapping mesh methodology is utilized to investigate the effect of an upstream turbulent wake on a three-dimensional oscillating NACA0012 extruded airfoil. A direct numerical simulation (DNS) at Reynolds Number 44,000 is performed for steady inflow incident upon the airfoil oscillating between angle of attack 5.6 and 25 degrees with reduced frequency k=0.16. Results are contrasted with subsequent DNS of the same oscillating airfoil in a turbulent wake generated by a stationary upstream cylinder.
ContributorsMerrill, Brandon Earl (Author) / Peet, Yulia (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Kostelich, Eric (Committee member) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2016
Description
A numerical study of chemotaxis in 3D turbulence is presented here. Direct Numerical
Simulation were used to calculate the nutrient uptake for both motile and non-motile bacterial
species and by applying the dynamical systems theory the effect of flow topology on the
variability of chemotaxis is analyzed. It is done by injecting a highly localized patch of nutrient
in the turbulent flow, and analyzing the evolution of reaction associated with the observed
high and low stretching regions. The Gaussian nutrient patch is released at different locations
and the corresponding nutrient uptake is obtained. The variable stretching characteristics of
the flow is depicted by Lagrangian Coherent Structures and the roles they play in affecting the
uptake are analyzed. The Lagrangian Coherent Structures are quantified by the Finite Time
Lyapunov Exponents which is a measure of the average stretching experienced by the flow in
finite time. It is found that in high stretching regions, the motile bacteria are attracted to the
nutrient patch very quickly, but also dispersed quickly; whereas in low stretching regions the
bacteria respond slower towards the nutrient patch. However the total uptake is intricately
determined by stretching history. These reaction characteristics are reflected in the several
realizations of simulations. This helps in understanding turbulence intensity and how it affects
the uptake of the nutrient.
Simulation were used to calculate the nutrient uptake for both motile and non-motile bacterial
species and by applying the dynamical systems theory the effect of flow topology on the
variability of chemotaxis is analyzed. It is done by injecting a highly localized patch of nutrient
in the turbulent flow, and analyzing the evolution of reaction associated with the observed
high and low stretching regions. The Gaussian nutrient patch is released at different locations
and the corresponding nutrient uptake is obtained. The variable stretching characteristics of
the flow is depicted by Lagrangian Coherent Structures and the roles they play in affecting the
uptake are analyzed. The Lagrangian Coherent Structures are quantified by the Finite Time
Lyapunov Exponents which is a measure of the average stretching experienced by the flow in
finite time. It is found that in high stretching regions, the motile bacteria are attracted to the
nutrient patch very quickly, but also dispersed quickly; whereas in low stretching regions the
bacteria respond slower towards the nutrient patch. However the total uptake is intricately
determined by stretching history. These reaction characteristics are reflected in the several
realizations of simulations. This helps in understanding turbulence intensity and how it affects
the uptake of the nutrient.
ContributorsGeorge, Jino (Author) / Tang, Wenbo (Thesis advisor) / Peet, Yulia (Thesis advisor) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2017
Description
Autonomic closure is a recently-proposed subgrid closure methodology for large eddy simulation (LES) that replaces the prescribed subgrid models used in traditional LES closure with highly generalized representations of subgrid terms and solution of a local system identification problem that allows the simulation itself to determine the local relation between each subgrid term and the resolved variables at every point and time. The present study demonstrates, for the first time, practical LES based on fully dynamic implementation of autonomic closure for the subgrid stress and the subgrid scalar flux. It leverages the inherent computational efficiency of tensorally-correct generalized representations in terms of parametric quantities, and uses the fundamental representation theory of Smith (1971) to develop complete and minimal tensorally-correct representations for the subgrid stress and scalar flux. It then assesses the accuracy of these representations via a priori tests, and compares with the corresponding accuracy from nonparametric representations and from traditional prescribed subgrid models. It then assesses the computational stability of autonomic closure with these tensorally-correct parametric representations, via forward simulations with a high-order pseudo-spectral code, including the extent to which any added stabilization is needed to ensure computational stability, and compares with the added stabilization needed in traditional closure with prescribed subgrid models. Further, it conducts a posteriori tests based on forward simulations of turbulent conserved scalar mixing with the same pseudo-spectral code, in which velocity and scalar statistics from autonomic closure with these representations are compared with corresponding statistics from traditional closure using prescribed models, and with corresponding statistics of filtered fields from direct numerical simulation (DNS). These comparisons show substantially greater accuracy from autonomic closure than from traditional closure. This study demonstrates that fully dynamic autonomic closure is a practical approach for LES that requires accuracy even at the smallest resolved scales.
ContributorsStallcup, Eric Warren (Author) / Dahm, Werner J.A. (Thesis advisor) / Herrmann, Marcus (Committee member) / Calhoun, Ronald (Committee member) / Kim, Jeonglae (Committee member) / Kostelich, Eric J. (Committee member) / Arizona State University (Publisher)
Created2020