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.

Pavement surface temperature is calculated using a fundamental energy balance model developed previously. It can be studied using a one-dimensional mathematical model. The input to the model is changed, to study the effect of different properties of pavement on its diurnal surface temperatures. It is observed that the pavement surface temperature has a microclimatic effect on the air temperature above it. A major increase in local air temperature is caused by heating of solid surfaces in that locality. A case study was done and correlations have been established to calculate the air temperature above a paved surface. Validation with in-situ pavement surface and air temperatures were made. Experimental measurement for the city of Phoenix shows the difference between the ambient air temperature of the city and the microclimatic air temperature above the pavement is approximately 10 degrees Fahrenheit. One mitigation strategy that has been explored is increasing the albedo of the paved surface. Although it will reduce the pavement surface temperature, leading to a reduction in air temperature close to the surface, the increased pavement albedo will also result in greater reflected solar radiation directed towards the building, thus increasing the building solar load. The first effect will imply a reduction in the building energy consumption, while the second effect will imply an increase in the building energy consumption. Simulation is done using the EnergyPlus tool, to find the microclimatic effect of pavement on the building energy performance. The results indicate the cooling energy savings of an office building for different types of pavements can be variable as much as 30%.

This study uses the Weather Research and Forecasting (WRF) model to simulate and predict the changes in local climate attributed to the urbanization for five desert cities. The simulations are performed in the fashion of climate downscaling, constrained by the surface boundary conditions generated from high resolution land-use maps. For each city, the land-use maps of 1985 and 2010 from Landsat satellite observation, and a projected land-use map for 2030, are used to represent the past, present, and future. An additional set of simulations for Las Vegas, the largest of the five cities, uses the NLCD 1992 and 2006 land-use maps and an idealized historical land-use map with no urban coverage for 1900.

The study finds that urbanization in Las Vegas produces a classic urban heat island (UHI) at night but a minor cooling during the day. A further analysis of the surface energy balance shows that the decrease in surface Albedo and increase effective emissivity play an important role in shaping the local climate change over urban areas. The emerging urban structures slow down the diurnal wind circulation over the city due to an increased effective surface roughness. This leads to a secondary modification of temperature due to the interaction between the mechanical and thermodynamic effects of urbanization.

The simulations for the five desert cities for 1985 and 2010 further confirm a common pattern of the climatic effect of urbanization with significant nighttime warming and moderate daytime cooling. This effect is confined to the urban area and is not sensitive to the size of the city or the detail of land cover in the surrounding areas. The pattern of nighttime warming and daytime cooling remains robust in the simulations for the future climate of the five cities using the projected 2030 land-use maps. Inter-city differences among the five urban areas are discussed.

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.

With a ground-based Doppler lidar on the upwind side of a wind farm in the Tehachapi Pass of California, radial wind velocity measurements were collected for repeating sector sweeps, scanning up to 10 kilometers away. This region consisted of complex terrain, with the scans made between mountains. The dataset was utilized for techniques being studied for short-term forecasting of wind power by correlating changes in energy content and of turbulence intensity by tracking spatial variance, in the wind ahead of a wind farm. A ramp event was also captured and its propagation was tracked.

Orthogonal horizontal wind vectors were retrieved from the radial velocity using a sector Velocity Azimuth Display method. Streamlines were plotted to determine the potential sites for a correlation of upstream wind speed with wind speed at downstream locations near the wind farm. A "virtual wind turbine" was "placed" in locations along the streamline by using the time-series velocity data at the location as the input to a modeled wind turbine, to determine the extractable energy content at that location. The relationship between this time-dependent energy content upstream and near the wind farm was studied. By correlating the energy content with each upstream location based on a time shift estimated according to advection at the mean wind speed, several fits were evaluated. A prediction of the downstream energy content was produced by shifting the power output in time and applying the best-fit function. This method made predictions of the power near the wind farm several minutes in advance. Predictions were also made up to an hour in advance for a large ramp event. The Magnitude Absolute Error and Standard Deviation are presented for the predictions based on each selected upstream location.

The role of environmental factors that influence atmospheric propagation of sound originating from freeway noise sources is studied with a combination of field experiments and numerical simulations. Acoustic propagation models are developed and adapted for refractive index depending upon meteorological conditions. A high-resolution multi-nested environmental forecasting model forced by coarse global analysis is applied to predict real meteorological profiles at fine scales. These profiles are then used as input for the acoustic models. Numerical methods for producing higher resolution acoustic refractive index fields are proposed. These include spatial and temporal nested meteorological simulations with vertical grid refinement. It is shown that vertical nesting can improve the prediction of finer structures in near-ground temperature and velocity profiles, such as morning temperature inversions and low level jet-like features. Accurate representation of these features is shown to be important for modeling sound refraction phenomena and for enabling accurate noise assessment. Comparisons are made using the acoustic model for predictions with profiles derived from meteorological simulations and from field experiment observations in Phoenix, Arizona. The challenges faced in simulating accurate meteorological profiles at high resolution for sound propagation applications are highlighted and areas for possible improvement are discussed.



A detailed evaluation of the environmental forecast is conducted by investigating the Surface Energy Balance (SEB) obtained from observations made with an eddy-covariance flux tower compared with SEB from simulations using several physical parameterizations of urban effects and planetary boundary layer schemes. Diurnal variation in SEB constituent fluxes are examined in relation to surface layer stability and modeled diagnostic variables. Improvement is found when adapting parameterizations for Phoenix with reduced errors in the SEB components. Finer model resolution (to 333 m) is seen to have insignificant ($<1\sigma$) influence on mean absolute percent difference of 30-minute diurnal mean SEB terms. A new method of representing inhomogeneous urban development density derived from observations of impervious surfaces with sub-grid scale resolution is then proposed for mesoscale applications. This method was implemented and evaluated within the environmental modeling framework. Finally, a new semi-implicit scheme based on Leapfrog and a fourth-order implicit time-filter is developed.

A municipal electric utility in Mesa, Arizona with a peak load of approximately 85 megawatts (MW) was analyzed to determine how the implementation of renewable resources (both wind and solar) would affect the overall cost of energy purchased by the utility. The utility currently purchases all of its energy through long term energy supply contracts and does not own any generation assets and so optimization was achieved by minimizing the overall cost of energy while adhering to specific constraints on how much energy the utility could purchase from the short term energy market. Scenarios were analyzed for a five percent and a ten percent penetration of renewable energy in the years 2015 and 2025. Demand Side Management measures (through thermal storage in the City's district cooling system, electric vehicles, and customers' air conditioning improvements) were evaluated to determine if they would mitigate some of the cost increases that resulted from the addition of renewable resources.

In the 2015 simulation, wind energy was less expensive than solar to integrate to the supply mix. When five percent of the utility's energy requirements in 2015 are met by wind, this caused a 3.59% increase in the overall cost of energy. When that five percent is met by solar in 2015, it is estimated to cause a 3.62% increase in the overall cost of energy. A mix of wind and solar in 2015 caused a lower increase in the overall cost of energy of 3.57%. At the ten percent implementation level in 2015, solar, wind, and a mix of solar and wind caused increases of 7.28%, 7.51% and 7.27% respectively in the overall cost of energy.

In 2025, at the five percent implementation level, wind and solar caused increases in the overall cost of energy of 3.07% and 2.22% respectively. In 2025, at the ten percent implementation level, wind and solar caused increases in the overall cost of energy of 6.23% and 4.67% respectively.

Demand Side Management reduced the overall cost of energy by approximately 0.6%, mitigating some of the cost increase from adding renewable resources.

Climate change has been one of the major issues of global economic and social concerns in the past decade. To quantitatively predict global climate change, the Intergovernmental Panel on Climate Change (IPCC) of the United Nations have organized a multi-national effort to use global atmosphere-ocean models to project anthropogenically induced climate changes in the 21st century. The computer simulations performed with those models and archived by the Coupled Model Intercomparison Project - Phase 5 (CMIP5) form the most comprehensive quantitative basis for the prediction of global environmental changes on decadal-to-centennial time scales. While the CMIP5 archives have been widely used for policy making, the inherent biases in the models have not been systematically examined. The main objective of this study is to validate the CMIP5 simulations of the 20th century climate with observations to quantify the biases and uncertainties in state-of-the-art climate models. Specifically, this work focuses on three major features in the atmosphere: the jet streams over the North Pacific and Atlantic Oceans and the low level jet (LLJ) stream over central North America which affects the weather in the United States, and the near-surface wind field over North America which is relevant to energy applications. The errors in the model simulations of those features are systematically quantified and the uncertainties in future predictions are assessed for stakeholders to use in climate applications. Additional atmospheric model simulations are performed to determine the sources of the errors in climate models. The results reject a popular idea that the errors in the sea surface temperature due to an inaccurate ocean circulation contributes to the errors in major atmospheric jet streams.

The implications of a changing climate have a profound impact on human life, society, and policy making. The need for accurate climate prediction becomes increasingly important as we better understand these implications. Currently, the most widely used climate prediction relies on the synthesis of climate model simulations organized by the Coupled Model Intercomparison Project (CMIP); these simulations are ensemble-averaged to construct projections for the 21st century climate. However, a significant degree of bias and variability in the model simulations for the 20th century climate is well-known at both global and regional scales. Based on that insight, this study provides an alternative approach for constructing climate projections that incorporates knowledge of model bias. This approach is demonstrated to be a viable alternative which can be easily implemented by water resource managers for potentially more accurate projections. Tests of the new approach are provided on a global scale with an emphasis on semiarid regional studies for their particular vulnerability to water resource changes, using both the former CMIP Phase 3 (CMIP3) and current Phase 5 (CMIP5) model archives. This investigation is accompanied by a detailed analysis of the dynamical processes and water budget to understand the behaviors and sources of model biases. Sensitivity studies of selected CMIP5 models are also performed with an atmospheric component model by testing the relationship between climate change forcings and model simulated response. The information derived from each study is used to determine the progressive quality of coupled climate models in simulating the global water cycle by rigorously investigating sources of model bias related to the moisture budget. As such, the conclusions of this project are highly relevant to model development and potentially may be used to further improve climate projections.

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.