In this study, atmospheric and turbulent flow solutions are obtained using meteorological flowfield and topographical terrain data in high resolution. The terrain topology of interest, which may be obtained in various resolution levels, is accurately modeled using structured or unstructured grids depending on whether high-rise building models are present or not. Meteorological weather prediction software MM5, is used to provide accurate and unsteady boundary conditions for the solution domain. Unsteady turbulent flow solutions are carried out via FLUENT with the help of several User Defined Functions developed. Unsteady flow solutions over topographical terrain of METU campus are computed with 25m x 25m x 15m resolution using structured grids. These FLUENT solutions are compared with the MM5 solutions. Also, the accuracy of the boundary layer velocity profiles is assessed. Finally, effects of surface roughness model extracted from MM5 for the region of interest is investigated. In addition, unsteady flow solutions over METU campus are repeated in presence of high-rise building models using unstructured grids with resolution varying from 5 meters around buildings to 80 meters further away. The study shows that unsteady, turbulent flow solutions can be accurately obtained using low resolution atmospheric weather prediction models and high resolution Navier-Stokes solutions over topographical terrains.

A mesh-less method for solving incompressible, multi-phase flow problems has been developed and is discussed along with the presentation of benchmark results showing good agreement with theoretical and experimental results. Results of a systematic, parametric study of the single phase flow around a 2D circular cylinder at Reynolds numbers up to 1000 are presented and discussed. Simulation results show good agreement with experimental results. Extension of the method to deal with multiphase flow including liquid-to-vapor phase transition along with applications to cavitating flow are discussed. Insight gleaned from numerical experiments of the cavity closure problem are discussed along with recommendations for additional research. Several conclusions regarding the use of the method are made.

A computer code is developed using finite volume technique for solving steady twodimensional and axisymmetric compressible Euler and Navier-Stokes equations for internal flows by &ldquo / multi-block&rdquo / technique. For viscous flows, both laminar and turbulent flow properties can be used. Explicit one step second order accurate Lax-Wendroff scheme is used for time integration. Inviscid solutions are verified by comparing the results of test cases of a support project which was supported by ONERA/France for Turkey T-108, named &ldquo / 2-D Internal Flow Applications for Solid Propellant Rocket Motors&rdquo / . For laminar solutions, analytical flat plate solution is used for planar case and theoretical pipe flow solution is used for axisymmetric case for verification. Prandtl turbulent flow analogy is used in a flat plate solution to verify the turbulent viscosity calculation. The test cases solved with single-block code are compared with the ones solved with multi-block technique to verify the multi-block algorithm and good similarity is observed between single-block solutions and multi-block solutions. For the burning simulation of propellant of Solid Propellant Rocket Motors, injecting boundary is used. Finally, a segmented solid propellant rocket motor case is solved to show the multi-block algorithm&rsquo / s flexibility in solving complex geometries.

A high resolution spectrally accurate three-dimensional flow solver is developed in order to simulate convection dominated fluid flows. The governing incompressible Navier Stokes equations along with the energy equation for temperature are discretized using a second-order accurate projection method which utilizes Adams Bashforth and Backward Differentiation formula for temporal discretization of the non-linear convective and linear viscous terms, respectively. Spatial discretization is performed using a Fourier/Chebyshev spectral method. Extensive tests on three-dimensional Taylor Couette flow are performed and it is shown that the method successfully captures the different states ranging from formation of Taylor vortices to wavy vortex regime. Next, the code is validated for convection dominated flows through a comprehensive comparison of the results for two dimensional Rayleigh Benard convection with the theoretical and experimental results from the literature. Finally, fully parallel simulations, with efficient utilization of computational resources and memory, are performed on a model three-dimensional axially homogeneous Rayleigh Benard convection problem in order to explore the high Rayleigh number flows and to test the scaling of global properties.