The continued research and development of high-order methods in Computational Fluid Dynamics (CFD) is primarily motivated by their potential to significantly reduce the computational cost and memory usage required to obtain a solution to a desired level of accuracy. In this work, a high-order Central Essentially Non-Oscillatory (CENO) finite-volume scheme is developed for the Euler and Navier-Stokes equations in three dimensions. The proposed CENO scheme is based on a hybrid solution reconstruction procedure using a fixed central stencil. A solution smoothness indicator facilitates the hybrid switching between a high-order k-exact reconstruction technique, and a monotonicity preserving limited piecewise linear reconstruction algorithm. The resulting scheme is applied to the compressible forms of the Euler and Navier-Stokes equations in three dimensions. The latter of which includes the application of this high-order work to the Large Eddy Simulation (LES) of turbulent non-reacting flows.

A 3d parallel CFD code is written to investigate the characteristics of and differences between Large Eddy Simulation (LES) models in the context of simulating a thermal buoyant plume. An efficient multigrid scheme is incorporated to solve the Poisson equation, resulting from the fractional step, projection method used to solve the Low Mach Number (LMN) Navier-Stokes equations. A wide range of LES models are implemented, including a variety of eddy models, structure models, mixed models and dynamic models, for both the momentum stresses and the temperature fluxes. Generalised gradient flux models are adapted from their RANS counterparts, and also tested. A number of characteristics are observed in the LES models relating to the thermal plume simulation in particular and turbulence in general. Effects on transition, dissipation, backscatter, equation balances, intermittency and energy spectra are all considered, as are the impact of the governing equations, the discretisation scheme, and the effect of grid coarsening. Also characteristics to particular models are considered, including the subgrid kinetic energy for the one-equation models, and constant histories for dynamic models. The argument that choice of LES model is unimportant is shown to be incorrect as a general statement, and a recommendation for when the models are best used is given.

In this work we try to analyse the dynamics of the Navier-Stokes equations in a problem without domain complexities as is the case of the plane Poiseuille flow. The Poiseuille problem is described as the flow of a viscous incompressible fluid, in a channel between two infinite parallel plates. We have considered it in two dimensions for the most common boundary conditions used to drive the fluid: mean constant pressure gradient or constant flux through the channel. We also specify the relation between this two formulations.We give the details of the direct numerical solution of the full two-dimensional, time-dependent, incompressible Navier-Stokes equations, formulated by means of spectral methods on the spatial variables and finite differences for time. Unlike other authors we have considered the classical formulation in terms of primitive variables for velocity and pressure. We also describe the approach adopted to eliminate the pressure and the cross-stream component of the velocity, obtaining thus a reduced system of ordinary differential equations from an original system of differential-algebraic equations. This is translated to a reduction of two thirds in the dimension of the original system and, in addition, it allows us to study the stability of fixed points by means of the analytical Jacobian matrix.We reproduce previous calculations on travelling waves (which are time-periodic orbits) and its stability to superharmonic disturbances. These solutions are observed as stationary in a Galilean reference in the streamwise direction. We begin by reviewing some results of the Orr-Sommerfeld equation which serve as a starting point to obtain the bifurcating solutions of time-periodic flows for several values of the periodic length in the streamwise direction. In turn, we also calculate several Hopf bifurcations that appear on the branch of periodic flows, for both cases of imposed constant flux and pressure.Likewise, for each unstable periodic flow, we study the connection of its unstable manifold to other attracting solutions.Starting at the Hopf bifurcations found for periodic flows, we analyse the bifurcating branches of quasi-periodic solutions at the two first Hopf bifurcations for the case of imposed constant pressure and the first one for constant flux. Those solutions are found as fixed points of an appropriate Poincaré map since, by the symmetry of the channel, they may be viewed as periodic flows in an appropriate moving frame of reference. We also study their stability by analysing the linear part of the Poincaré map. In the case of constant flux we have found a branch of quasi-periodic solutions which, on increasing the Reynolds number, changes from stable to unstable, giving rise to an attracting family of quasi-periodic flows with 3 frequencies. The results referring to the first Hopf bifurcation for constant pressure, are not in qualitative agreement with those of Soibelman & Meiron (1991),which yield a different bifurcation picture and stability properties for the obtained quasi-periodic flows. From the computed unstable flows we follow their unstable invariant manifold and describe what new attracting solution they are conducted to.

En la presente tesis se han estudiado métodos de paso fraccionado para la resolución numérica de la ecuación de Navier-Stokes incompresible mediante el método de los elementos finitos; dicha ecuación rige el movimiento de un fluido incompresible viscoso. Partiendo del análisis del método de proyección clásico, se desarrolla un método para el problema de Stokes (lineal y estacionario) con iguales propiedades en cuanto a discretizacion espacial que aquel, explicando así sus propiedades de estabilización de la presión. Se da también una extensión del nuevo método a la ecuación de Navier-Stokes incompresible estacionaria (no lineal).En la segunda parte de la tesis, se desarrolla un método de paso fraccionado para el problema de evolución que supera un inconveniente del método de proyección relativo a la imposición de las condiciones de contorno.Para todos los métodos desarrollados, se demuestran teoremas de convergencia y estimaciones de error, se proponen implementaciones eficientes y se proporcionan numerosos resultados numéricos.

Hydropower is a more environmental friendly way of producing electric power than many other alternatives today. Though, the effects of constructing mega dams are much tangible for the local eco systems in addition to changing many people’s lives forever. In order to prevent floods, riverbank erosions or landslides, proper investigations of the environmental impact from dam constructions must be performed. One of the key parameters in such investigations is the flow discharge velocity. This master thesis treats experimental measurements and numerical simulations of the velocity downstream a model of Xiluodu dam. The Xiluodu dam is a mega dam under construction in China and will have a total capacity of 12 600 MW when completed. The model is in scale 1:100 and the experiments have been performed at Department of Hydraulic Engineering, Tsinghua University, Beijing, China. The velocity profile shows that the velocity in the middle of the river is larger than the velocity at the surface and near the riverbank. The comparison between the measured and the simulated velocities shows a difference of less than 20 percent in almost all points which can be considered as a good result. In those points where the difference is more than 20 percent, this is believed to be due to the position of these points. Some of them were located near a vortex and others very close to the bottom. This is a problem when sparsely measured topography in combination with linear interpolation makes the boundaries of the simulations incorrect. In order to perform better simulations, more densely topography data and better flow boundary conditions should be used. More measuring points of the velocity could also improve the result.

This work develops a nonlinear, three-dimensional spectral collocation method for the simulation of the incompressible Navier-Stokes equations for geophysical and environmental flows. These flows are often driven by the interaction of stratified fluid with topography, which is accurately accounted for in this model using a mapped coordinate system. The spectral collocation method used here evaluates derivatives with a Fourier trigonometric or Chebyshev polynomial expansion as appropriate, and it evaluates the nonlinear terms directly on a collocated grid. The coordinate mapping renders ineffective fast solution methods that rely on separation of variables, so to avoid prohibitively expensive matrix solves this work develops a low-order finite-difference preconditioner for the implicit solution steps. This finite-difference preconditioner is itself too expensive to apply directly, so it is solved pproximately with a geometric multigrid method, using semicoarsening and line relaxation to ensure convergence with locally anisotropic grids. The model is discretized in time with a third-order method developed to allow variable timesteps. This multi-step method explicitly evaluates advective terms and implicitly evaluates pressure and viscous terms. The model’s accuracy is demonstrated with several test cases: growth rates of Kelvin-Helmholtz billows, the interaction of a translating dipole with no-slip boundaries, and the generation of internal waves via topographic interaction. These test cases also illustrate the model’s use from a high-level programming perspective. Additionally, the results of several large-scale simulations are discussed: the three-dimensional dipole/wall interaction, the evolution of internal waves with shear instabilities, and the stability of the bottom boundary layer beneath internal waves. Finally, possible future developments are discussed to extend the model’s capabilities and optimize its performance within the limits of the underlying numerical algorithms.

<p>A new method of stabilizing low-order, proper orthogonal decomposition based reduced-order models of the Navier Stokes equations is proposed. Unlike traditional approaches, this method does not rely on empirical turbulence modeling or modification of the Navier-Stokes equations. It provides spatial basis functions different from the usual proper orthogonal decomposition basis function in that, in addition to optimally representing the solution, the new proposed basis functions also provide stable reduced-order models. The proposed approach is illustrated with two test cases: two-dimensional flow inside a square lid-driven cavity and a two-dimensional mixing layer.</p> / Dissertation

In the present study, a fifth-order weighted essentially non-oscillatory (WENO) scheme was built for solving the surface-capturing level-set equation. Combined with the level-set equation, the three-dimensional Reynolds averaged Navier-Stokes (RANS) equations were employed for the prediction of nonlinear wave-interaction and wave-breaking phenomena over sloping beaches. In the level-set finite-analytic Navier-Stokes (FANS) method, the free surface is represented by the zero level-set function, and the flows are modeled as immiscible air-water two phase flows. The Navier-Stokes equations for air-water two phase flows are formulated in a moving curvilinear coordinate system and discretized by a 12-point finite-analytical scheme using the finite-analytic method on a multi-block over-set grid system. The Pressure Implicit with Splitting of Operators / Semi-Implicit Method for Pressure-Linked Equation Revised (PISO/SIMPLER) algorithm was used to determine the coupled velocity and pressure fields. The evolution of the level-set method was solved using the third-order total variation diminishing (TVD) Runge-Kutta method and fifth-order WENO scheme. The accuracy was confirmed by solving the Zalesak's problem. Two major subjects are discussed in the present study. First, to identify the WENO scheme as a more accurate scheme than the essentially non-oscillatory scheme (ENO), the characteristics of a nonlinear monochromatic wave were studied systematically and comparisons of wave profiles using the two schemes were conducted. To eliminate other factors that might produce wave profile fluctuation, different damping functions and grid densities were studied. To damp the reflection waves efficiently, we compared five damping functions. The free-surface elevation data collected from gauges distributed evenly in a numerical wave tank are analyzed to demonstrate the damping effect of the beach. Second, as a surface-tracking numerical method built on curvilinear coordinates, the level-set RANS model was tested for nonlinear bichromatic wave trains and breaking waves on a sloping beach with a complex free surface. As the wave breaks, the velocity of the fluid flow surface became more complex. Numerical modeling was performed to simulate the two-phase flow velocity and its corresponding surface and evolution when the wave passed over different sloping beaches. The breaking wave test showed that it is an efficient technique for accurately capturing the breaking wave free surface. To predict the breaking points, different wave heights and beach slopes are simulated. The results show that the dependency of wave shape and breaking characteristics to wave height and beach slope match the results provided by experiments.

In this thesis, the two-dimensional, transient, laminar flow of viscous and incompressible fluids is solved by using the dual reciprocity boundary element method (DRBEM). Natural convection and mixed convection flows are also solved with the addition of energy equation. Solutions of natural convection flow of nanofluids and micropolar fluids in enclosures are obtained for highly large values of Rayleigh number. The fundamental solution of Laplace equation is used for obtaining boundary element method (BEM) matrices whereas all the other terms in the differential equations governing the flows are considered as nonhomogeneity. This is the main advantage of DRBEM to tackle the nonlinearities in the equations with considerably small computational cost. All the convective terms are evaluated by using the DRBEM coordinate matrix which is already computed in the formulation of nonlinear terms. The resulting systems of initial value problems with respect to time are solved with forward and central differences using relaxation parameters, and the fourth-order Runge-Kutta method. The numerical stability analysis is developed for the flow problems considered with respect to the choice of the time step, relaxation parameters and problem constants. The stability analysis is made through an eigenvalue decomposition of the final coefficient matrix in the DRBEM discretized system. It is found that the implicit central difference time integration scheme with relaxation parameter value close to one, and quite large time steps gives numerically stable solutions for all flow problems solved in the thesis. One-and-two-sided lid-driven cavity flow, natural and mixed convection flows in cavities, natural convection flow of nanofluids and micropolar fluids in enclosures are solved with several geometric configurations. The solutions are visualized in terms of streamlines, vorticity, microrotation, pressure contours, isotherms and flow vectors to simulate the flow behaviour.

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.