Numerical simulation of randomly forced turbulent flows

Several authors have proposed studying randomly forced turbulent hows (e.g., E. A. Novikov, Soviet Physics JETP, 20(5), 1290 1965). More recently, theoretical investigations (e.g., renormalization group) have focused on whim-noise forced Navier-Stokes equations (V. Yakhot and S. A. Orszag, J.Sci.Comput. 1(1), 3 1986), The present article aims to provide an appropriate numerical method for the simulation of randomly forced turbulent systems. The spatial discretization is based on the classical Fourier spectral method. The time integration is performed by a second-order Runge-Kutta scheme. The consistency of an extension of this scheme to treat additive noise stochastic differential equations is proved. The random number generator is based on lagged Fibonacci series. Results are presented for two randomly forced problems: the Burgers and the incompressible Navier-Stokes equations with a white-noise in time forcing term characterized by a power-law correlation function in spectral space. A variety of statistics are computed for both problems, including the structure functions, The third-order structure functions are compared with their exact expressions in the inertial subrange. The influence of the dissipation mechanism (viscous or hyperviscous) on the inertial subrange is discussed. In particular, probability density functions of velocity increments are computed for the Navier-Stokes simulation, Finally, for both Burgers and Navier-Stokes problems, our results support the view that random sweeping is the dominant effect of the large-scale motion on the small-scales. (C) 1998 Academic Press.

Simulation of benchmark and industrial unsteady compressible turbulent fluid flows

Michele Mossi

Large-eddy simulation (LES) is a very promising technique for the numerical computation of unsteady turbulent flows, and seems to be the perfect tool to simulate the compressible air flow around a high-speed train in a tunnel, providing unsteady results for aerodynamic and aeroacoustic analysis. To look into this possible future application of LES, two major lines of investigation are pursued in this thesis: first, the study of the effective ability of shock-capturing schemes to predict fundamental turbulent phenomena; second, the analysis of the aerodynamic phenomena induced by a high-speed train in a tunnel. The numerical simulation of compressible flows requires the use of shock-capturing schemes. These schemes can be relatively dissipative and mask the subgrid-scale contribution introduced in a large-eddy simulation to account for the unresolved turbulence scales. To estimate their effective dissipation and their ability to resolve turbulence phenomena, shock-capturing schemes widely used for aeronautical applications, from second- to fifth-order space accuracy, are employed here for simulating well-known fundamental flows in subsonic and supersonic regimes. Direct and large-eddy numerical simulations are undertaken for the inviscid and viscous Taylor-Green vortex decay problem, the freely decaying homogeneous and isotropic turbulence, and the fully developed channel flow. For all the turbulent flows investigated, several turbulence statistics are computed and the numerical dissipation of the schemes tested is interpreted in terms of subgrid-scale dissipation in a LES context, yielding an equivalent Smagorinsky or dynamic coefficient. This coefficient is for all schemes of the same order of magnitude as the commonly accepted values in LES for the subgrid-scale term. On the grounds of this analysis and of the comparisons of the turbulence statistics with accurate data obtained in the literature, the addition of explicit subgrid-scale models to the shock-capturing schemes tested is not recommended. It is thus concluded that the use of the LES technique for compressible turbulent flows is not yet suitable for industrial applications. The aerodynamic phenomena generated by a high-speed train travelling in a tunnel are also discussed, their importance on the design of high-speed lies is pointed out, and the analysis tools commonly employed for their study are reviewed. To study numerically the three-dimensional, compressible and turbulent air flow around a high-speed train accelerating in a tunnel, by accounting for the unsteady effects at inlet and outlet boundaries due to the propagation of pressure waves generated at the train departure, new coupling conditions between one-dimensional and three-dimensional domains are developed. These conditions are applied successfully to the numerical simulation of the unsteady wake developing behind two- and three-dimensional vehicles, where the averaged Navier-Stokes equations are solved with the turbulence modelling approach. The influence on the wake of the length of the vehicle tail is also discussed and results of multi-dimensional simulations are compared with one-dimensional data.

EPFL1999

Analysis and numerical simulation of free surface flows

Alexandre Caboussat

Mathematical and numerical aspects of free surface flows are investigated. On one hand, the mathematical analysis of some free surface flows is considered. A model problem in one space dimension is first investigated. The Burgers equation with diffusion has to be solved on a space interval with one free extremity. This extremity is unknown and moves in time. An ordinary differential equation for the position of the free extremity of the interval is added in order to close the mathematical problem. Local existence in time and uniqueness results are proved for the problem with given domain, then for the free surface problem. A priori and a posteriori error estimates are obtained for the semi-discretization in space. The stability and the convergence of an Eulerian time splitting scheme are investigated. The same methodology is then used to study free surface flows in two space dimensions. The incompressible unsteady Navier-Stokes equations with Neumann boundary conditions on the whole boundary are considered. The whole boundary is assumed to be the free surface. An additional equation is used to describe the moving domain. Local existence in time and uniqueness results are obtained. On the other hand, a model for free surface flows in two and three space dimensions is investigated. The liquid is assumed to be surrounded by a compressible gas. The incompressible unsteady Navier-Stokes equations are assumed to hold in the liquid region. A volume-of-fluid method is used to describe the motion of the liquid domain. The velocity in the gas is disregarded and the pressure is computed by the ideal gas law in each gas bubble trapped by the liquid. A numbering algorithm is presented to recognize the bubbles of gas. Gas pressure is applied as a normal force on the liquid-gas interface. Surface tension effects are also taken into account for the simulation of bubbles or droplets flows. A method for the computation of the curvature is presented. Convergence and accuracy of the approximation of the curvature are discussed. A time splitting scheme is used to decouple the various physical phenomena. Numerical simulations are made in the frame of mould filling to show that the influence of gas on the free surface cannot be neglected. Curvature-driven flows are also considered.

EPFL2004

Parallel Algorithms for the Solution of Large-Scale Fluid-Structure Interaction Problems in Hemodynamics

Davide Forti

This thesis addresses the development and implementation of efficient and parallel algorithms for the numerical simulation of Fluid-Structure Interaction (FSI) problems in hemodynamics. Indeed, hemodynamic conditions in large arteries are significantly affected by the interaction of the pulsatile blood flow with the arterial wall. The simulation of fluid-structure interaction problems requires the approximation of a coupled system of Partial Differential Equations (PDEs) and the set up of efficient numerical solution strategies. Blood is modeled as an incompressible Newtonian fluid whose dynamics is governed by the Navier-Stokes equations. Different constituive models are used to describe the mechanical response of the arterial wall; specifically, we rely on hyperelastic isotropic and anistotropic material laws. The finite element method is used for the space discretization of both the fluid and structure problems. In particular, for the Navier-Stokes equations we consider a semi-discrete formulation based on the Variational Multiscale (VMS) method. Among a wide range of possible solution strategies for the FSI problem, here we focus on strongly coupled monolithic approaches wherein the nonlinearities are treated in a fully implicit mode. To cope with the high computational complexity of the three dimensional FSI problem, a parallel solution framework is often mandatory. To this end, we develop a new block parallel preconditioner for the coupled linearized FSI system obtained after space and time discretization. The proposed preconditioner, named FaCSI, exploits the factorized form of the FSI Jacobian matrix, the use of static condensation to formally eliminate the interface degrees of freedom of the fluid equations, and the use of a SIMPLE preconditioner for unsteady Navier-Stokes equations. In FSI problems, the different resolution requirements in the fluid and structure physical domains, as well as the presence of complex interface geometries make the use of matching fluid and structure meshes problematic. In such situations, it is much simpler to deal with discretizations that are nonconforming at the interface, provided however that the matching conditions at the interface are properly fulfilled. In this thesis we develop a novel interpolation-based method, named INTERNODES, for numerically solving partial differential equations by Galerkin methods on computational domains that are split into two (or several) subdomains featuring nonconforming interfaces. By this we mean that either a priori independent grids and/or local polynomial degrees are used to discretize each subdomain. INTERNODES can be regarded as an alternative to the mortar element method: it combines the accuracy of the latter with the easiness of implementation in a numerical code. The aforementioned techniques have been applied for the numerical simulation of large-scale fluid-structure interaction problems in the context of biomechanics. The parallel algorithms developed showed scalability up to thousands of cores utilized on high performance computing machines.

EPFL2016