Vincenzo Citro

Global stability and sensitivity analysis of boundary-layer flows past a hemispherical roughness element

We study the full three-dimensional instability mechanism past a hemispherical roughness element immersed in a laminar Blasius boundary layer. The inherent three-dimensional flow pattern beyond the Hopf bifurcation is characterized by coherent vortical structures usually called hairpin vortices. Direct numerical simulation results are used to analyze the formation and the shedding of hairpin vortices inside the shear layer. The first bifurcation is investigated by global-stability tools. We show the spatial structure of the linear direct and adjoint global eigenmodes of the linearized Navier-Stokes equations and use the structural-sensitivity field to locate the region where the instability mechanism acts. The core of this instability is found to be symmetric and spatially localized in the region immediately downstream of the roughness element. The effect of the variation of the ratio between the obstacle height k and the boundary layer thickness δk∗ is also considered. The resulting bifurcation scenario ...

Three-dimensional stability, receptivity and sensitivity of non-Newtonian flows inside open cavities

We investigate the stability properties of flows over an open square cavity for fluids with shear-dependent viscosity. Analysis is carried out in context of the linear theory using a normal-mode decomposition. The incompressible Cauchy equations, with a Carreau viscosity model, are discretized with a finite-element method. The characteristics of direct and adjoint eigenmodes are analyzed and discussed in order to understand the receptivity features of the flow. Furthermore, we identify the regions of the flow that are more sensitive to spatially localized feedback by building a spatial map obtained from the product between the direct and adjoint eigenfunctions. Analysis shows that the first global linear instability of the steady flow is a steady or unsteady three-dimensionl bifurcation depending on the value of the power-law index n. The instability mechanism is always located inside the cavity and the linear stability results suggest a strong connection with the classical lid-driven cavity problem.

Efficient stabilization and acceleration of numerical simulation of fluid flows by residual recombination

Abstract The study of the stability of a dynamical system described by a set of partial differential equations (PDEs) requires the computation of unstable states as the control parameter exceeds its critical threshold. Unfortunately, the discretization of the governing equations, especially for fluid dynamic applications, often leads to very large discrete systems. As a consequence, matrix based methods, like for example the Newton–Raphson algorithm coupled with a direct inversion of the Jacobian matrix, lead to computational costs too large in terms of both memory and execution time. We present a novel iterative algorithm, inspired by Krylov-subspace methods, which is able to compute unstable steady states and/or accelerate the convergence to stable configurations. Our new algorithm is based on the minimization of the residual norm at each iteration step with a projection basis updated at each iteration rather than at periodic restarts like in the classical GMRES method. The algorithm is able to stabilize any dynamical system without increasing the computational time of the original numerical procedure used to solve the governing equations. Moreover, it can be easily inserted into a pre-existing relaxation (integration) procedure with a call to a single black-box subroutine. The procedure is discussed for problems of different sizes, ranging from a small two-dimensional system to a large three-dimensional problem involving the Navier–Stokes equations. We show that the proposed algorithm is able to improve the convergence of existing iterative schemes. In particular, the procedure is applied to the subcritical flow inside a lid-driven cavity. We also discuss the application of Boostconv to compute the unstable steady flow past a fixed circular cylinder (2D) and boundary-layer flow over a hemispherical roughness element (3D) for supercritical values of the Reynolds number. We show that Boostconv can be used effectively with any spatial discretization, be it a finite-difference, finite-volume, finite-element or spectral method.

Stability and Sensitivity Analysis of the Secondary Instability in the Sphere Wake

The three-dimensional flow past a fixed sphere placed within a uniform stream is investigated. This paper focuses on the second bifurcation, which is responsible for the onset of the unsteadiness. ...

Computation of the bluff-body sound generation by a self-consistent mean flow formulation

The sound generated by the flow around a circular cylinder is numerically investigated by using a finite-element method. In particular, we study the acoustic emissions generated by the flow past the bluff body at low Mach and Reynolds numbers. We perform a global stability analysis by using the compressible linearized Navier-Stokes equations. The resulting direct global mode provides detailed information related to the underlying hydrodynamic instability and data on the acoustic field generated. In order to recover the intensity of the produced sound, we apply the self-consistent model for non-linear saturation proposed by Mantic-Lugo, Arratia, and Gallaire [“Self-consistent mean flow description of the nonlinear saturation of the vortex shedding in the cylinder wake,” Phys. Rev. Lett. 113, 084501 (2014)]. The application of this model allows us to compute the amplitude of the resulting linear mode and the effects of saturation on the mode structure and acoustic field. Our results show excellent agreement...