Conforming Chebyshev spectral collocation methods for the solution of laminar flow in a constricted channel

The numerical simulation of steady planar two-dimensional, laminar flow of an incompressible fluid through an abruptly contracting channel using spectral domain decomposition methods is described. The key features of the method are the decomposition of the flow region into a number of rectangular subregions and spectral approximations which are pointwise C(1) continuous across subregion interfaces. Spectral approximations to the solution are obtained for Reynolds numbers in the range 0 to 500. The size of the salient corner vortex decreases as the Reynolds number increases from 0 to around 45. As the Reynolds number is increased further the vortex grows slowly. A vortex is detected downstream of the contraction at a Reynolds number of around 175 that continues to grow as the Reynolds number is increased further

Fourier spectral computation of geometrically confined two-dimensional flows

Large-scale flow phenomena in the atmosphere and the oceans are predominantly two-dimensional (2D) due to the large aspect ratio of the typical horizontal and vertical length scales in the flow. The 2D nature of large-scale geophysical flows motivates the use of a conceptual approach known as "2D turbulence". It usually involves the (forced/damped) Navier-Stokes equations on a square domain with periodic boundaries or on a spherical surface. This setup may be useful for numerical studies of atmospheric flow. For the oceans, on the other hand, geometrical confinement due to the continental shelves is of crucial importance. The physically most relevant boundary condition for oceanographic flow is probably the no-slip condition. Previous numerical and experimental studies have shown that confinement by no-slip boundaries dramatically affects the dynamics of (quasi-)2D turbulence due to its role as vorticity source. An important process is the detachment of high-amplitude vorticity filaments from the no-slip sidewalls that subsequently affect the internal flow. The first part of the thesis concerns the development and extensive testing of a Fourier spectral scheme for 2D Navier-Stokes flow in domains bounded by rigid noslip walls. An advantage of Fourier methods is that higher-order accuracy can, in principle, be achieved. Moreover, these methods are fast, relatively easy to implement even for performing parallel computations. The no-slip boundary condition is enforced by using an immersed boundary technique called "volume-penalization". In this method an obstacle with no-slip boundaries is modelled as a porous medium with a small permeability. It has recently been shown that in the limit of infinitely small permeability the solution of the penalized Navier-Stokes equations converges towards the solution of the Navier-Stokes equations with no-slip boundaries. Therefore the penalization error can be controlled with an arbitrary parameter. A possible drawback is that the sharp transition between the fluid and the porous medium can trigger Gibbs oscillations that might deteriorate the stability and accuracy of the scheme. Using a very challenging dipole-wall collision as a benchmark problem, it is, however, shown that higher-order accuracy is retrieved by using a novel 159 post-processing procedure to remove the Gibbs effect. The second topic of the thesis is the dynamics of geometrically confined 2D turbulent flows. The role of the geometry on the flow development has been studied extensively. For this purpose high resolution Fourier spectral simulations have been conducted where different geometries are implemented by using the volumepenalization method. A quantity that is of particular importance on a bounded domain is the angular momentum. On a circular domain production of angular momentum is virtually absent. Therefore the amount of angular momentum carried by the initial flow has important consequences for the evolution of the flow. The results of the simulations are consistent with previous numerical and experimental work on this topic performed in a lower Reynolds number regime. The typical vortex structures of the late time evolution of the flow are explained by means of a minimum enstrophy principle and the presence of weak viscous dissipation. For an elliptic geometry it is shown that strong spin-up events of the flow occur even for small eccentricities. The spin-up phenomenon can be related to the role of the pressure along the boundary of the domain. It is found that the magnitude of the torque exerted on the internal fluid can be scaled with the eccentricity. Furthermore, it is observed that angular momentum production in a non circular geometry is not restricted to moderate Reynolds numbers. Significantly higher Reynolds number flow computations in a square geometry clearly reveal strong and rapid spin-up of the flow. Finally the scale-dependence of the vorticity and velocity statistics in forced 2D turbulence on a bounded domain has been studied. A challenging aspect is that a statistically steady state can be achieved by a balance between the injection of kinetic energy by the external forcing and energy dissipation at the no-slip sidewalls. It is important to note that on a double periodic domain a steady state is usually achieved by introducing volumetric drag forces. Several studies reported that this strongly affects the spatial scaling behaviour of the flow. Therefore it is very interesting to quantify the small-scale statistics in the bulk of statistically steady flow on a domain with no-slip boundaries in the absence of bottom drag. It is observed that the internal flow shows extended self-similar, locally homogeneous and isotropic scaling behaviour at small scales. It is further demonstrated that a direct enstrophy cascade develops in the interior of the flow domain. Some deviations from the classical scaling theory of 2D turbulence developed independently by Kraichnan, Batchelor and Leith may be associated to the presence of coherent structures in the flow. It is, however, anticipated that higher-resolution simulations are required in order to draw more decisive conclusions. The parallel Fourier spectral scheme with volume-penalization is very suitable for pursuing such simulations on high performance machines in the near future. In summary the thesis contributes to both the development of numerical techniques and understanding of wall-bounded two-dimensional flows. The Fourier spectral scheme with volume-penalization is found very suitable for pursuing direct numerical simulations in complex geometries. The high-resolution simulations considered in the thesis clearly reveal that spontaneous production of angular momentum due to interaction with non-circular domain boundaries is present for significantly higher Reynolds numbers than considered previously

Effect of high order interpolation in the stability and efficiency of the time-integration process in vorticity-velocity CFD algorithms

The numerical solution of the incompressible Navier-Stokes Equations offers an effective alternative to the experimental analysis of Fluid-Structure interaction i.e. dynamical coupling between a fluid and a solid which otherwise is very complex, time consuming and very expensive. To have a method which can accurately model these types of mechanical systems by numerical solutions becomes a great option, since these advantages are even more obvious when considering huge structures like bridges, high rise buildings, or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the KLE to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ODE time integration schemes and thus allows us to tackle the various multiphysics problems as separate modules. The current algorithm for KLE employs a structured or unstructured mesh for spatial discretization and it allows the use of a self-adaptive or fixed time step ODE solver while dealing with unsteady problems. This research deals with the analysis of the effects of the Courant-Friedrichs-Lewy (CFL) condition for KLE when applied to unsteady Stoke’s problem. The objective is to conduct a numerical analysis for stability and, hence, for convergence. Our results confirmthat the time step ∆t is constrained by the CFL-like condition ∆t ≤ const. hα, where h denotes the variable that represents spatial discretization

Effect of mesh distortion on the accuracy of high order vorticity-velocity CFD approaches

The numerical solution of the incompressible Navier-Stokes equations offers an alternative to experimental analysis of fluid-structure interaction (FSI). We would save a lot of time and effort and help cut back on costs, if we are able to accurately model systems by these numerical solutions. These advantages are even more obvious when considering huge structures like bridges, high rise buildings or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the Kinematic Laplacian Equation (KLE) to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ordinary differential equations (ODE) time integration schemes, allowing us to tackle each problem as a separate module. The current algortihm for the KLE uses an unstructured quadrilateral mesh, formed by dividing each triangle of an unstructured triangular mesh into three quadrilaterals for spatial discretization. This research deals with determining a suitable measure of mesh quality based on the physics of the problems being tackled. This is followed by exploring methods to improve the quality of quadrilateral elements obtained from the triangles and thereby improving the overall mesh quality. A series of numerical experiments were designed and conducted for this purpose and the results obtained were tested on different geometries with varying degrees of mesh density