Matching Items (7)
Filtering by
- All Subjects: Computational Fluid Dynamics
- All Subjects: Wind power plants--Environmental aspects.
- Creators: Calhoun, Ronald
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
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
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
This research work uses the Weather Research and Forecasting Model to study the effect of large wind farms with an area of 900 square kilometers and a high power density of 7.58 W/m2 on regional climate. Simulations were performed with a wind farm parameterization scheme turned on in south Oregon. Control cases were also run with the parameterization scheme turned off. The primary emphasis was on offshore wind farms. Some analysis on onshore wind farms was also performed. The effects of these wind farms were studied on the vertical profiles of temperature, wind speed, and moisture as well as on temperature and on wind speed near the surface and at hub height. The effects during the day and at night were compared. Seasonal variations were also studied by performing simulations in January and in July. It was seen that wind farms produce a reduction in wind speed at hub height and that the downward propagation of this reduction in wind speed lessens as the atmosphere becomes more stable. In all the cases studied, the wind farms produced a warming effect near the surface, with greater atmospheric stability leading to higher near-surface temperatures. It was also observed that wind farms caused a drying effect below the hub height and a moistening effect above it, because they had facilitated vertical transport of moisture in the air from the lower layers of the atmosphere to the layers of the atmosphere above the wind farm.
ContributorsGeorge, Sushant (Author) / Huang, Huei-Ping (Thesis advisor) / Wang, Zhihua (Committee member) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2016
Description
Due to decrease in fossil fuel levels, the world is shifting focus towards renewable sources of energy. With an annual average growth rate of 25%, wind is one of the foremost source of harnessing cleaner energy for production of electricity. Wind turbines have been developed to tap power from wind. As a single wind turbine is insufficient, multiple turbines are installed forming a wind farm. Generally, wind farms can have hundreds to thousands of turbines concentrated in a small region. There have been multiple studies centering the influence of weather on such wind farms, but no substantial research focused on how wind farms effect local climate. Technological advances have allowed development of commercial wind turbines with a power output greater than 7.58 MW. This has led to a reduction in required number of turbines and has optimized land usage. Hence, current research considers higher power density compared to previous works that relied on wind farm density of 2 to 4 W/m 2 . Simulations were performed using Weather Research and Forecasting software provided by NCAR. The region of simulation is Southern Oregon, with domains including both onshore and offshore wind farms. Unlike most previous works, where wind farms were considered to be on a flat ground, effects of topography have also been considered here. Study of seasonal effects over wind farms has provided better insight into changes in local wind direction. Analysis of mean velocity difference across wind farms at a height of 10m and 150m gives an understanding of wind velocity profiles. Results presented in this research tends to contradict earlier belief that velocity reduces throughout the farm. Large scale simulations have shown that sometimes, more than 50% of the farm can have an increased wind velocity of up to 1m/s
at an altitude of 10m.
at an altitude of 10m.
ContributorsKadiyala, Yogesh Rao (Author) / Huang, Huei-Ping (Thesis advisor) / Rajagopalan, Jagannathan (Committee member) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2015
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