Response of a plate in turbulent channel flow: Analysis of fluid–solid coupling
Introduction
The coupling between a turbulent flow and the resulting structural excitation is a problem of interest in marine, civil and aerospace engineering. In this paper, we investigate this coupling in a canonical setting — linear one-way coupled (fluid affects solid, but not vice versa) response of an elastic plate in turbulent channel flow (Pope, 2001) due to wall-pressure fluctuations alone. Specifically, we address the question — how much do the fluid sources at different wall-normal locations contribute to the plate excitation for different frequencies, and what are the salient features of these fluid sources? We answer this question with a novel formulation that combines Direct Numerical Simulation (DNS) data, Green’s function formulation and spectral Proper Orthogonal Decomposition (POD). For brevity, we will sometimes refer to wall-pressure fluctuations as just wall-pressure.
The one-way coupling between the fluid sources and plate excitation can be broken into two parts: (i) fluid source–wall-pressure fluctuation coupling, and (ii) wall-pressure fluctuation–plate excitation coupling. Note that we neglect the wall-shear stress contribution to the plate forcing. We further classify the techniques to investigate the fluid source–wall-pressure fluctuation coupling into — scaling variables-based, Green’s function-based and conditional averaging-based techniques. We discuss some features of the wall-pressure fluctuation sources identified by each of these techniques.
Identification of the scaling variables for the power-spectral density (PSD)/wavenumber spectrum of wall-pressure fluctuation yields qualitative information of the wall-normal region of the fluid sources. The wall-pressure PSD in the low (), mid () and high frequency ranges () scale with the potential flow variables (), outer flow variables (), and inner flow variables (), respectively (Farabee and Casarella, 1991), where is the angular frequency, is the fluid density, is the boundary layer thickness, is the displacement thickness of the boundary layer, is the centerline velocity, is the wall-shear stress, is the friction velocity, is the kinematic viscosity of the fluid, and the friction Reynolds number is defined as . Thus, the sources responsible for the low, mid and high frequency wall-pressure fluctuations are predominantly in the potential, outer and inner region of the turbulent boundary layer, respectively.
The Green’s function-based techniques (Chang III et al., 1999, Anantharamu and Mahesh, 2020) yield quantitative information of the sources of wall-pressure fluctuation. The premultiplied streamwise wavenumber spectrum and the PSD of the wall-pressure fluctuations in a turbulent channel show peaks at (Panton et al., 2017) and (Hu et al., 2006) for , respectively, where is the streamwise wavelength, and + indicates normalization with viscous units ( and ). The dominant contributors to this inner peak are in the buffer region of the channel (Anantharamu and Mahesh, 2020). The approach of Anantharamu and Mahesh (2020) that identified this dominant contribution (i) combines DNS data with the Green’s function formulation to express the wall-pressure PSD () as integrated contribution () from all wall-parallel plane pairs, , (ii) accounts for the relative phase difference between the contributions from different wall-parallel planes neglected in the previous Green’s function approach of Chang III et al. (1999), and (iii) yields a distribution of sources in the wall-normal direction instead of a wall-normal region as indicated by the scaling variables. Further, the methodology identified decorrelated features of wall-pressure fluctuation sources using spectral Proper Orthogonal Decomposition (POD). The identified dominant wall-pressure source at the linear and premultiplied wall-pressure PSD peak frequency resembled tall and inclined patterns, respectively.
The conditional averaging-based technique (Ghaemi and Scarano, 2013) yields patterns of the flow structure that are correlated to a particular wall-pressure fluctuation event. The time history of the wall-pressure fluctuation signal at a point on the wall shows occasional positive and negative high amplitude wall-pressure peaks. The conditionally averaged flow fields show coupling between a hairpin vortex and the high amplitude peaks (Ghaemi and Scarano, 2013). The flow structure responsible for the positive and negative high amplitude wall-pressure peak at a point are the sweep and ejection event occurring above it, respectively. The ejection event responsible for the negative peak occurs upstream of the hairpin head in between the quasi-streamwise vortices. The sweep event that leads to the positive peak occurs downstream of the hairpin head.
The dynamic linear elasticity equations describe the wall-pressure fluctuation–plate excitation coupling. This one-way coupled FSI approach is valid for small linear deformation () of the plate, where is the wall-normal displacement. The approach generally uses (i) plate theories (e.g. Poisson–Kirchhoff) to describe the deformation, and modal superposition to obtain the response, (ii) frequency domain since steady state response is usually the quantity of interest, and (iii) a model wavenumber-frequency spectrum (Corcos, 1964, Chase, 1980, Hwang, 1998) for the spatially homogeneous wall-pressure fluctuations as input. Note that the model wavenumber-frequency spectrum usually requires a model PSD (Bull, 1967, Smol’Iakov and Tkachenko, 1991, Goody, 2004). The mode shapes and natural frequencies of the plate required to perform modal superposition can be obtained analytically for simple boundary conditions and geometry. For complicated boundary conditions and geometry, Finite Element Method (FEM) is used to compute the modal decomposition.
The wall-pressure fluctuation–plate excitation coupling has been previously investigated in wavenumber space (Hwang and Maidanik, 1990, Blake, 2017). The modal force PSD of the plate can be expressed as the wavenumber integral (Blake, 2017) where and are the origins of the plate in the streamwise and spanwise directions, and are the lengths of the plate in the streamwise and spanwise directions, is the modal force PSD of the th mode shape, is the wall-pressure wavenumber-frequency spectrum and is the modal shape function. From the above equation, we observe that the modal shape function couples the wall-pressure wavenumber-frequency spectrum to the modal force. The relative contribution of different wavenumber regions to the modal force spectra depends on the mode order (), boundary conditions, and the ratio of the streamwise modal wavenumber () to the convective wavenumber at the natural frequency of the mode (Hwang and Maidanik, 1990). The high streamwise wavenumber () portion of decays as , and for clamped, simply supported and free boundary conditions on all edges (Blake, 2017). Thus, plates with free boundary conditions accept more of the high streamwise wavenumber component of the wall-pressure fluctuations. Further, special wall-pressure fluctuation models that separately approximate the high and low wavenumber portion of the wall-pressure fluctuation wavenumber-frequency spectrum can be derived and used to obtain the response of plates (Hambric et al., 2004). Hambric et al. (2004) showed good agreement between FEM response of a plate excited by the modified Corcos model of Hwang (1998) and an equivalent edge forcing model which only models the convective component in the modified Corcos model for a plate with three edges clamped and one edge free. This shows the importance of the convective region of wall-pressure fluctuation spectrum for plates with free boundary conditions. For a plate with all four edges clamped, FEM response from a low wavenumber excitation model showed good agreement with the modified Corcos model, thus highlighting the dominance of low wavenumber contribution for clamped boundary condition.
Experiments by Zhang et al. (2017) have shown coupling between flow structures and the response of a compliant wall in a turbulent channel flow. The large positive and negative deformation of the compliant wall is coupled to the ejection and sweep events, respectively, occurring above it (Zhang et al., 2017). Conditionally averaged flow fields show that these events are related to the high amplitude pressure peaks and hairpin vortices that surround the local deformation of the compliant wall. For large deformation of the compliant wall, the plate deflection affects the near-wall turbulence. The compliant wall deflection into the buffer layer breaks the near-wall streaks and the associated quasi-streamwise vortices, and induces more spanwise coherence (Rosti and Brandt, 2017).
In this paper, we develop a formulation to obtain the wall-normal distribution of intensity and relative phase of the fluid sources responsible for the plate excitation. Previous research works do not yield such quantitative information of the fluid sources. The main idea is to express the plate averaged displacement PSD as a double wall-normal integral of the ‘net displacement source’ cross-spectral density (CSD) across the height of the channel. The analysis framework combines the volumetric DNS data, Green’s function solution of the pressure fluctuation and modal superposition, and builds on the previous work of Anantharamu and Mahesh (2020). We then apply the framework to explain the one-way coupled FSI simulation results of an elastic plate in turbulent channel flow at and . The fluid and solid simulations make use of finite volume DNS and time-domain FEM, respectively. Further, the decorrelated fluid sources that contribute the most to plate response are obtained using spectral POD of the net displacement source CSD.
The organization of the paper is as follows: In Section 2, we describe the computational domain, mesh resolution, and the FSI simulation details. Section 3 discusses the novel one-way coupling analysis framework. In Section 4.1, we discuss the obtained one-way coupled FSI results. We compare the fluid and solid length scales for the first three modes in Section 4.2. Section 4.3 discusses the spectral features of the net displacement source CSD and in Section 4.4, we identify the decorrelated features of the fluid source using spectral POD. Finally, we summarize the results in Section 5.
Note that , and denote the streamwise, wall-normal and spanwise coordinates, respectively. Superscripts/subscripts and denote fluid and solid quantities, respectively.
Section snippets
Computational domain
Fig. 1 shows a schematic of the fluid and solid computational domain and Table 1 shows the domain extents. The fluid computational domain is a Cartesian box of size . We choose , and , where is the half channel height. Long streamwise and spanwise domains include the contribution of large scale structures to pressure fluctuations. The solid computational domain is a rectangular plate clamped on all sides placed at the bottom wall of the channel. The plate is
Theory
The goal is to express the plate averaged displacement PSD as a double integral over all the wall-parallel plane pairs. To accomplish this, we first express the steady-state component of the bottom wall displacement as a wall-normal integral, Here, is called the ‘net displacement source’ (exact form is derived later). It gives the contribution of each wall-parallel plane to the surface displacement of the plate. We define the plate averaged
FSI simulation results
Fig. 2 shows an instantaneous visualization of the FSI simulation. The vertical and horizontal slices show the fluid streamwise velocity and wall-pressure fluctuations, respectively. The center patch denotes the deformed plate. The isosurfaces are of Q-criterion at non-dimensional values of 500 and 1000. The colored overlayed on the isosurface denotes the streamwise component of vorticity. We use different colormaps for each quantity. The instantaneous field clearly shows the fine scales
Summary
In summary, we present a novel framework to investigate the fluid–solid coupling in a canonical setting — linear one-way coupled excitation of an elastic plate in turbulent channel flow. We apply the framework to explain the response of a clamped plate obtained using the in-house FSI solver — MPCUGLES-SOLID at and .
The obtained plate response at and with fixed non-dimensional Young’s modulus have overlapping plate averaged low frequency spectrum in outer units.
CRediT authorship contribution statement
Sreevatsa Anantharamu: Formal analysis, Writing - original draft. Krishnan Mahesh: Writing -review & editing, Supervision.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgments
This work is supported by the United States Office of Naval Research (ONR) under grant N00014-17-1-2939 with Dr. Ki-Han Kim as the technical monitor. The computations were made possible through computing resources provided by the US Army Engineer Research and Development Center (ERDC) in Vicksburg, Mississippi on the Cray machines, Copper and Onyx of the High Performance Computing Modernization Program . We also thank for the computing resources provided by the US Air Force Research Laboratory
References (27)
- et al.
Turbulent channel flow simulations in convecting reference frames
J. Comput. Phys.
(2013) Modeling the wavevector-frequency spectrum of turbulent boundary layer wall pressure
J. Sound Vib.
(1980)- et al.
Vibrations of plates with clamped and free edges excited by low-speed turbulent boundary layer flow
J. Fluids Struct.
(2004) - et al.
A wavenumber analysis of the coupling of a structural mode and flow turbulence
J. Sound Vib.
(1990) - et al.
A scaled thickness conditioning for solid-and solid-shell discretizations of thin-walled structures
Comput. Methods Appl. Mech. Engrg.
(2011) - et al.
A numerical method for large-eddy simulation in complex geometries
J. Comput. Phys.
(2004) - et al.
Analysis of wall-pressure fluctuation sources from direct numerical simulation of turbulent channel flow
J. Fluid Mech.
(2020) Finite Element Procedures
(2006)- et al.
Random Data: Analysis and Measurement Procedures, Vol. 729
(2011) Mechanics of Flow-Induced Sound and Vibration, Volume 1 and 2
(2017)
Wall-pressure fluctuations associated with subsonic turbulent boundary layer flow
J. Fluid Mech.
Relationship between wall pressure and velocity-field sources
Phys. Fluids
The structure of the turbulent pressure field in boundary-layer flows
J. Fluid Mech.
Cited by (7)
Direct numerical simulation-based vibroacoustic response of plates excited by turbulent wall-pressure fluctuations
2023, Journal of Fluid MechanicsNonlinear Vibration and Dynamic Bifurcation of Axially Moving Plates Under Subsonic Airflow in a Narrow Space
2024, Iranian Journal of Science and Technology - Transactions of Mechanical EngineeringGlobal stability analysis and direct numerical simulation of boundary layers with an isolated roughness element
2022, Journal of Fluid MechanicsLarge-eddy simulation for the aero-vibro-acoustic analysis: plate-cavity system excited by turbulent channel flow
2022, Acta Mechanica Sinica/Lixue Xuebao