Next Article in Journal
Disinfection of Escherichia coli by Mixing with Bulk Ultrafine Bubble Solutions
Next Article in Special Issue
Surrogate Models for Heat Transfer in Oscillating Flow with a Local Heat Source
Previous Article in Journal
Six-Field Theory for a Polyatomic Gas Mixture: Extended Thermodynamics and Kinetic Models
Previous Article in Special Issue
CFD Modeling of Wind Turbine Blades with Eroded Leading Edge
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coupled Simulation of Fluid Flow and Vibro-Acoustic Processes in the Channel with a Circular Cylinder

by
Konstantin Volkov
Institute of Hydrodynamics and Control Processes, St. Petersburg State Marine Technical University, 190121 St. Petersburg, Russia
Fluids 2022, 7(12), 382; https://doi.org/10.3390/fluids7120382
Submission received: 8 November 2022 / Revised: 5 December 2022 / Accepted: 7 December 2022 / Published: 11 December 2022
(This article belongs to the Special Issue Industrial CFD and Fluid Modelling in Engineering)

Abstract

:
Vibro-acoustic processes are an interacting set of pulsations of the working fluid and vibrations of mechanical structural elements. The simulation of vibro-acoustic processes in a long pipe with an elastic round cylinder is considered. The mathematical model is developed in a coupled formulation, when not only pressure pulsations cause pipe vibrations, but also vibrations of the mechanical subsystem affect sound wave propagation in the working fluid. The influence of vortex formation processes in the channel on the system dynamics is taken into account. The fluid flow is found using delayed detached eddy simulation. The flow regimes around a single round cylinder corresponding to various Reynolds numbers are investigated to validate the computational algorithm. The distributions of the flow quantities and vibro-acoustic behavior of the system are discussed.

1. Introduction

Turbulent pressure fluctuations on a streamlined surface play an important role in acoustic measurements in moving media. In particular, in flexible extended towed antennas, which are used in seismic exploration to search for minerals under the seabed, sensitive receiving elements are enclosed in an acoustically transparent cylindrical shell and are located close enough to the streamlined surface. Sources of hydrodynamic noise and vibrations are the phenomena associated with the flow around various structural elements [1] (boundary layer formed on a streamlined surface, separation regions near bluff bodies, and jet streams). An important role in excitation of vibrations is played by pressure fluctuations in the boundary layer on the walls and acoustic oscillations excited by the separation of vortices. Pseudo-sonic pulsations of near-wall pressure on the streamlined surface of carriers of hydro-acoustic equipment lead to the appearance of hydrodynamic noise.
The generation of coherent vortices and their interaction with each other and with the streamlined surface lead to the emission of powerful pressure pulsations into the environment. The variety of scales of vortex structures in a turbulent boundary layer and their transfer rates lead to a broadband spectrum of velocity and pressure fluctuations. The appearance of dominant vortex formations, generated in one way or another, determines the presence of discrete or tonal peaks in the spectral dependencies of pressure fluctuations measured both above the streamlined surface and on its surface.
The hydrodynamic noise of the turbulent boundary layer, the sources of which are coherent vortex structures interacting with each other and with the streamlined surface, has sonic and pseudo-sonic components [2]. Unlike sound, which has a wave nature and propagates into the environment at the speed of sound, pseudo-sonic pressure pulsations are transferred at a speed close to the flow speed. Pseudo-sound generated by the nonlinear interaction of the vortex structures of the boundary layer does not obey the superposition principle. The pseudo-sound pressure does not depend on the average ambient pressure and decreases inversely with the square of the distance from the source. When gas flows around elongated bodies (strings, cables, long rods, high pipes), acoustic vibrations arise due to the separation of vortices (Aeolian tones). The appearance of Aeolian tones is associated with fluctuations in the lift force acting on a streamlined body and its drag [3].
Pseudo-sonic pressure fluctuations are localized within the turbulent flow, while sonic pressure fluctuations propagate over large distances outside the region occupied by the turbulent flow. The study of the physics of turbulent pressure fluctuations is important not only because they are sources of hydrodynamic noise. They affect the surface in a flow, generating vibrational oscillations, but they are also present in the correlation dependencies in the Reynolds stress transfer equations and the dissipation tensor energy.
One of the test problems that make it possible to determine the characteristic features of hydrodynamic acoustic processes is the problem of flow past a round cylinder. The fluid flow around a circular cylinder is considered in a large number of experimental and computational studies [4,5,6,7]. A round cylinder in a transverse flow of fluid is a typical element of many structures and is the object of intensive theoretical, numerical, and experimental research. The formation and shedding of vortices behind the cylinder and other bluff bodies leads to unwanted vibrations and destruction of the structure. In this regard, it is necessary to be able to control the process of vortex shedding, reduce the drag force of the streamlined body and the amplitude of fluctuations of the forces applied to it [8,9]. The results of calculations of the loads experienced by an obstacle during its flow are used to optimize the methods of attaching the cylinder to the channel walls.
Modeling and simulation of the flowfield around cylinders have been tackled by many researchers at various Reynolds and Mach numbers [10,11,12,13,14,15,16,17,18,19,20,21,22,23]. The results obtained form a basis for validation of new models and computational algorithms.
Studies of the flow around a cylinder near the wall include data on the effect of the gap between the cylinder and the wall on the frequency of vortex formation (von Karman road), as well as information on the interaction of shear layers on the surface of the cylinder and on the wall depending on the gap and the Reynolds number [5,24,25,26]. A system of large-scale horseshoe-shaped vortices is formed at the junction of the shedder body and the walls parallel to the flow. Such vortex formations have a significant impact on the flow, determining local shear, acoustic noise, drag and lift coefficients [27,28].
Pseudo-sonic pulsations of near-wall pressure behind an annular obstacle on a longitudinally streamlined flexible extended cylinder are studied in [29]. The integral and spectral statistical characteristics of the pressure fluctuation field behind the obstacle are obtained, and its influence on the structure of the turbulent boundary layer is studied. An increase in the obstacle diameter and flow velocity leads to an increase in the low-frequency spectral components of pressure fluctuations and a weakening of high-frequency pressure fluctuations in comparison with the boundary layer on a hydraulically smooth cylinder.
The results of theoretical, numerical, and experimental studies of the characteristics of pseudo-sonic pressure oscillations caused by the interaction of the flow and the sound field inside a three-dimensional spherical dimple are presented in [30]. Symmetric and asymmetric large-scale vortex systems inside the dimple are found, the existence of which depends on the flow regime, their location and periodicity of ejection to the outside are indicated.
An important task in assessing the vibration strength of technical systems is to determine their vibration characteristics. Many of the developed approaches are empirical in nature and are not universal. The available methods for calculating the vibration parameters of a system under its force loading by a pulsating flow of the working fluid are characterized by significant idealization. Mathematical modeling of vibro-acoustic processes in technical systems is an important task, the solution of which will make it possible to predict vibro-acoustic loads at the design stage and ensure the operability of the system by rationally choosing its parameters, locations, and characteristics of supports, devices for correcting dynamic characteristics (vibration dampers, dynamic vibration dampers, vibration dampers of the working liquids).
In [31], simulation of vibro-acoustic interaction processes is implemented for a rectangular channel with one compliant and three rigid bounding surfaces. The form of vibration of the compliant plate, excited by the pulsations of the pressure of the working medium, is a superposition of at least two natural modes of vibration. The possibilities of eddy-resolving approaches to modeling turbulent flow near a system of ordered circular cylinders and predicting their flow characteristics, including pressure fluctuations, are shown in [32]. One of the main requirements for the methods of solving vibro-acoustics problems is the high resolution of the difference scheme and its low dissipation in the calculation of convective flows, which is important for the correct transfer of acoustic disturbances with a small amplitude.
The solution of coupled problems is characterized by large velocity and pressure gradients. In such problems, there is also a wide range of flow quantity variations, which creates challenges for numerical methods. Many previous studies, for example, Refs. [29,30,31,32] are limited to the assumption that it is possible to describe acoustic processes in a working medium by a one-dimensional wave equation. In addition, in many previous works, an assumption is made about the harmonic law of changing the parameters of the working environment. Insufficient attention is also paid to the issues of determining variable stresses from the impact of vibro-acoustic loads when calculating the strength of systems.
In this study, large vortices of vibro-acoustic processes are simulated in a pipe with a transversely located round cylinder. The generalized mathematical model is considered in a coupled formulation, when not only pressure pulsations cause pipe vibration, but also vibrations of the mechanical subsystem affect wave processes in the working fluid. Based on the results of numerical simulation, the vibro-acoustic behavior of the system is discussed.

2. Coupled Problems

Vibro-acoustic problems (Acoustic–Structure Interaction, ASI) imply joint modeling of elastic waves in solids, acoustic pressure waves in a liquid, and taking into account the interaction between them. Acoustic waves in different types of materials are described by different equations.
The loads on the surfaces found from the CFD problem are applied to the boundaries of the computational domain. These loads define boundary conditions for the solution of the stress–strain state. With a single transfer of data on the impact of a flow on a certain structural element from a CFD package to an FEA package, the interaction is unidirectional and may be used when the reverse effect of the deformation process on the flow is weak. In the general case, a bidirectional data exchange is required (the CFD package transmits loads, but receives the displacement values of the nodes of the interface boundary).
With a single data exchange during a time step, an explicit scheme for pairing software systems is implemented. Multiple data exchange at the time step of solving the adjoint problem leads to an implicit conjugation scheme. The need for multiple data exchange is associated with the problem of stability of the calculation process, which is explained by the high inertia of the liquid due to its density and leads to large reaction forces acting from the flow on the moving interface. The iterative procedure [33,34] is used to obtain a convergent solution of the adjoint problem at a time step.
The need to modify the grid structure is due to the displacement of the boundaries that form the object and are the boundaries of the computational domain of the hydrodynamic problem. Modification of the grid consists in its local rebuilding or deformation (stretching/compression) of the set of grid lines while maintaining the relationships of nodes, edges, and faces of control volumes.
For situations that do not express a high degree of interaction nonlinearity, the explicit method seems to be more efficient and flexible, allowing flow calculations and finite element analysis to be performed independently and then iteratively linked. Pairing on each time step is implemented by a recursive procedure based on iterations between different calculated modules until the specified level of convergence is reached.

3. Governing Equations

Modeling of vibro-acoustic processes consists of description of the non-stationary movement of the working medium, the interaction of the working medium with moving or stationary surrounding or surrounding elements of the mechanical subsystem, their interaction with the environment, and other elements that, for some reason, were not included in this system.

3.1. Sound Propagation in Fluid

Flows are considered in whichever period of change of flow parameters is much greater than the frequency of sound propagating in the flowfield. The flow quantities are presented as sum of average and fluctuating components ( u = u ¯ + u ). The prime corresponds to small fluctuations of flow quantities.
The governing equations describing small fluctuations of the flow parameters have the form (the line above the average values is omitted)
v t + [ rot v , v ] + [ rot v , v ] + ( v , v ) = 1 ρ p + 1 ρ 2 ( p · ρ ) ; ρ t + ( v , ρ ) + ( v , ρ ) + ρ div v + ρ div v = 0 ; s t + ( v , s ) + ( v , s ) = 0 .
The equation of state is written as
p = c 2 ρ + h s , c 2 = p ρ s , h = p s ρ .
The system of Equation (1) is valid for a flow with non-uniform entropy ( s 0 ) and in the presence of vorticity ( rot v 0 ). Approximations consist only in the fact that only linear terms are taken into account, and irreversible processes in the sound wave are not taken into account; therefore, the propagation of sound waves is considered as an adiabatic process d ( s + s ) / d t = 0 . Auxiliary relations, which are necessary for the transition from entropy to temperature, have the form
s = c p T T β p ρ p , p s ρ = ρ c 2 c p β p T ,
where c p is the specific heat capacity, β p = ( 1 / V ) ( V / T ) p is the volume expansion coefficient.
Assuming that the flow is isentropic ( s = 0 ) and irrotational ( rot v = 0 ), for a small adiabatic compression or expansion of the gas, the equation of state takes the form
p = c 2 ρ .
The pressure potential, Π , and the velocity potential of sound vibrations, φ , are introduced using the relations
Π = p ρ , v = φ .
The equation for the velocity potential is
d 2 φ d t 2 = c 2 Δ φ + Π 0 , φ + d φ d t ( v , l n c 2 ) .
The potential of the initial flow is
Π 0 = d p ρ .
The convective derivative is defined as
d φ d t = φ t + ( v , φ ) .
It follows from Equation (2) that in the absence of vorticity and entropy gradient, no sound is generated in the medium.
In a medium at rest v = 0 and rot z , v = 0 , the wave equation for the velocity potential (2) takes the form
2 φ t 2 c 2 Δ φ = 0 ,
where c 2 = γ p / ρ is the speed of sound.
The boundary condition on the surface of an elastic body has the following form
φ n = v 0 n ,
where v 0 n is the normal velocity of the body surface. For an absolute rigid body,
φ n = 0 .
For the correct statement of the boundary value problem, it is necessary to set the initial conditions for φ and φ / t .

3.2. Oscillations of an Elastic Body

The equation of motion for an elastic body is written in the following form
ρ a = ρ F + div σ ,
where a is acceleration, F is external force, σ is stress tensor. Acceleration is represented in terms of derivatives of speed
a = d V d t = V t + V i V x i .
Slow deformations are considered (velocities V i are small) and it is assumed that a V / t . The equations in displacements u for an elastic body take the form
ρ 2 u t 2 = ρ F + div σ .
The stress tensor is determined by Hooke’s law for an isotropic elastic medium in the Lame form
σ = λ I trace ( ε ) + 2 μ ε .
The function trace ( ε ) denotes the trace of the strain tensor. The strain tensor has the following form
ε = 1 2 grad u + ( grad u ) .
The stress tensor satisfies the relation
σ = λ I div u + μ grad u + ( grad u ) .
The Lame coefficients are determined by the relations
μ = E 2 ( 1 + ν ) , λ = ν E ( 1 + ν ) ( 1 2 ν ) ,
where E is Young’s modulus, ν is Poisson’s ratio.
Substituting the stress tensor from Hooke’s law into the equations of motion, we obtain the equation in small displacements
ρ 2 u t 2 = ρ F + ( μ + λ ) grad div u + μ Δ u .
Equation (3) is used to describe small oscillations in an elastic body.

4. Coupling Scheme

The pressure fluctuations obtained from CFD calculations under the assumption of stationary pipe walls are not sound, since the signal drifts with the flow velocity. The fluid acts on the pipe walls and barriers, which are made of elastic material and begin to deform and oscillate. In a solid material, not only body waves propagate, as in a fluid, but also shear waves. Interaction of wall vibrations with the fluid leads to the generation of acoustic waves and their propagation. Acoustic waves also propagate in the surrounding air. Simulation of vibro-acoustic characteristics is a coupled problem, in which waves in a solid material and acoustic waves in the fluid interact. To simplify the model, it is assumed that acoustic waves have a small effect on the flow. In this case, the vibro-acoustic problem is considered as three independent problems: wave propagation in a solid material, acoustic wave propagation in water and ambient air.
At stage 1, a water flow in a channel with rigid walls and unsteady pressure and friction forces on the pipe walls are found. At stage 2, the problem of the propagation of vibro-acoustic waves is solved under the influence of non-stationary forces on the wall from the inside of the pipe.
The internal diameter of the channel is d = 50 mm. Wall thickness is h = 5 mm, and length of the channel is L = 500 mm. The density of water is ρ = 998.5 kg/m3. Kinematic viscosity of water is ν = 1.003 · 10−3 m/s2. The pipe ends are fixed (Figure 1). The diameter of cylinder is dc = 8 mm. The cylinder is located at l1 = 0.1 m from the inlet. The inlet water velocity is V = 5 m/s. Relative pressure at the outlet is zero. In the solution of the CFD problem, the pipe walls are stationary. No-slip and no-penetration boundary conditions are applied to the walls. Air pressure and temperature are p = 105 Pa and T = 300 K. Material of the pipe is steel. When calculating the vibro-acoustic problem on the inner walls of the pipe, the pressure fluctuations (pressure amplitudes for each frequency after the Fourier transform) are interpolated from the hydrodynamic problem.

5. Estimations of Flow Quantities

The flow regime is defined by the Reynolds number Re h = ρ V d h / μ , where ρ is density, V is flow velocity, d h is hydraulic diameter, μ is dynamic viscosity. If the flow velocity is V = 5 m/s, the Reynolds number equals Re h = 248 , 804 (turbulent regime of flow).
In fluid flows past a cylinder, the vortex shedding frequency depends on the Reynolds number [10,11,12]. The Strouhal number is found as St = f d c / V , where f is characteristic frequency, d c is cylinder diameter, V is flow velocity. The Reynolds number calculated from the cylinder diameter is Re c = 39,809 (turbulent regime of flow). For Reynolds numbers larger than about 1000, the Strouhal number remains almost constant ( St = 0.21 ). The vortex street during the flow in a pipe with a cylinder is the dominant non-stationary process; therefore, the fundamental frequency is approximately equal to f m = St V / d c = 131 Hz, and the oscillation period is T m = 7.62 × 10 3 s.
To describe a periodic process, it is necessary to measure at least two points per period (Nyquist’s theorem). The calculations use approximately 100 points per period to resolve the dominant frequency. The time step is chosen equal to Δ t = T m / 100 = 7.62 × 10 5 s. In this case, processes with frequency f max = 1 / ( 10 Δ t ) = 1312 Hz are calculated with 10 points per period.
To perform a Fourier transform on the fundamental tone of a signal, it is necessary to have data for a certain number of periods of this tone. A total of 100 periods of the pitch signal is stored in the calculations. After the oscillations reach a stable mode, the data are saved for t f = 100 T m = 0.76 s. In this case, it is necessary to perform N = t f / Δ t = 10 4 steps in time. During the time t f , there are 10 periods with frequency f min = 10 / t f = 12.5 Hz. Proper resolution is observed in the frequency range from f min to f max (greater errors occur above or below these frequencies).

6. Scheme of Simulation

The delayed detached-eddy simulation (DDES) method with a shear stress transport (SST) k ω turbulence model is used to resolve vortex structures. DDES is an alternative to Reynolds-averaged Navier–Stokes (RANS) equations and large-eddy simulation (LES). RANS models are not capable of providing acceptable accuracy for calculating flows with large separation zones [7], and LES imposes high requirements on computational resources when modelling separated flows, which makes their practical use difficult [19].
The following approach is used to simulate the unsteady flowfield. A preliminary stationary calculation of the flow in the pipe is performed using the RANS equations with the SST-model. The semi-implicit method for pressure-linked equations-consistent (SIMPLEC) method is used to integrate the equations. Then, unsteady calculations of the flow in the pipe are performed with DDES and the SST-model. The calculations are carried out until a steady state flow is realized. To integrate the equations of motion, a modified pressure implicit with a splitting of operator (PISO) method is used. The pressure correction step is performed more times by the SIMPLEC method. The time step is found from the stability conditions. An unsteady flowfield in a pipe with a fixed time step is found. The pressure values are preserved in each cell bordering the pipe and cylinder walls at every time step.
A block-structured mesh is used. Nodes are clustered near the cylinder surface, wake region, and pipe walls to resolve the boundary layer (Figure 2). In the pipe, DDES operates in the LES mode, and the mesh is uniform. In the DDES approach, inviscid fluxes are discretized with a hybrid scheme and the Roe method to maintain accuracy in the LES region and stability of the scheme in the RANS region with a non-uniform and coarse mesh. The distributions of flow quantities found from the RANS solution with the SST-model are used as initial conditions.
The mesh nodes are distributed uniformly along the circumferential coordinate. In the inner region ( 0.5 r 5 ), which is a circle, an O-type mesh is used. In the near-wall region ( 0.5 r 1 ), the minimum step along the radial coordinate is r min = 0.003 , the growth factor is q = 1.032 (about 16 nodes along the radial coordinate are used to resolve the boundary layer). In the outer region, the minimum step in radial direction is r min = 0.2 , and the growth factor is q = 1.04 . In the buffer region, which is located at a considerable distance from the cylinder and allows avoiding the reflection of sound waves from the boundaries of the computational domain, the maximum step in radial direction is r min = 20 , and the growth factor is q = 1.09 .
The total number of mesh nodes is about 1.5 million. In this case, the condition y + < 1.2 is satisfied in the near-wall region. The dimensionless time step is 0.1056. The calculations are carried out on a time interval equal to 162.
Calculations are based on a fully compressible solver [34]. Convergence to a steady state is accelerated by the use of multigrid techniques, and by the application of block-Jacobi preconditioning for high-speed flows, with a separate low Mach number preconditioning method for use with low-speed flows.
With a known flowfield, the calculations of pressure fluctuations that occur on the channel walls due to the influence of an unsteady fluid flow are performed. The boundary condition is set in the form of a known amplitude of pressure fluctuations for each frequency.
The fast Fourier transformation (FFT) allows to obtain dependencies of pressure on time over the pipe walls and cylinder. The data obtained are input for solving the vibro-acoustic problem. The solution of the fluid problem, which describes the distribution of amplitudes of pressure fluctuations along the walls of the pipe and cylinder, contains data with a fine frequency step. For vibro-acoustic calculation, such a set of frequencies is redundant, and therefore, the data are averaged over a coarser frequency grid ( 1 / n octave spectrum filters). The most common are 1/3-octave and 1/8-octave filters.

7. Flow around Cylinder

A numerical simulation of an unsteady flow of a viscous compressible fluid around a cylinder is considered. The calculation results are given for two Reynolds numbers equal to Re = 150 (regime 1) and Re = 3900 (regime 2), which are widely used for testing numerical methods in the literature. In both cases, the Mach number is set equal to M = 0.2 . The Reynolds numbers correspond to the flow regime near the cylinder with flow separation from its surface. Using the fast Fourier transform, the values of the Strouhal numbers, the amplitudes of oscillations of the aerodynamic coefficients are determined, and the phases of the formation of vortices and separation of the flow are found from the values of the maximum amplitudes.

7.1. Computational Domain and Boundary Conditions

A transverse flow around a circular cylinder located in a viscous compressible gas, which is at rest for t 0 , is simulated. For t > 0 , the gas is set in motion with a constant speed U 0 . When calculating the Reynolds number, the diameter of the cylinder and the velocity of the undisturbed flow are used as characteristic parameters of the problem. The equations are solved in the region external to the cylinder boundary.
The center of the cylinder with diameter D is located at the origin. In the ( x , y ) plane, the computational domain has the shape of a rectangle ( 20 < x < 80 , 20 < y < 20 ). The extent of the computational domain in the direction of the z axis is 4 D .
The potential flow field is used as the initial conditions, which is modified in the boundary layer region in such a way as to satisfy the no-slip conditions on the cylinder surface. At the inlet boundary, the velocity of the undisturbed flow is set, which is found from the given Reynolds number, and the temperature is equal to 290 K. At the outlet boundary of the computational domain, non-reflecting boundary conditions are used. A static pressure equal to 10 5 Pa is fixed at the upper and lower boundaries of the computational domain. The surface of the cylinder is assumed to be thermally insulated. In the direction of the z axis, periodic boundary conditions are applied.

7.2. Regime 1

The distribution of pressure over time at various points on the surface of the cylinder is shown in Figure 3. The lines correspond to the numerical simulation data, and the symbols correspond to the calculated data of [10]. There is a fairly good agreement between the calculation results and the data of [10]. The amplitude of pressure fluctuations at θ = 90 ° (symbols △, solid line) and θ = + 90 ° (symbols ◻, dotted line) is several times greater than the amplitude pressure fluctuations observed at θ = 0 ° (symbols ♢, dash-dotted line).
The distributions of the lift coefficient, drag coefficient, and pressure coefficient are shown in Figure 4. The lines correspond to the numerical simulation data, the symbols ∘ correspond to the data of [10]. It achieves good agreement with the available data for all characteristics.
The characteristics of the flow in regime 1 are in good agreement with the available numerical and experimental data (Table 1, where C p b , C L , and C D are the mean pressure coefficient, the mean square lift coefficient, and the mean drag coefficient). In particular, in calculations the Strouhal number is 0.186, while the data of a physical experiment give a value of 0.184 [8] and 0.183 [12]. The angular coordinate of the flow separation point is θ s = 108.5 ° . Averaging is carried out over a time interval equal to 150. To find the Strouhal number, the correlation relation St = 0.2684 1.0356 Re 1 / 2 from [15] is used. The Strouhal number equals St = 0.1838 for Re = 150 , which is in good agreement with the obtained value. The results of calculations of the Strouhal number in comparison with experimental data are shown in Figure 5.
The instantaneous pressure distribution near a circular cylinder is shown in Figure 6, which clearly shows oscillations of the wake behind the cylinder.

7.3. Regime 2

The characteristics of the flow in mode 2 are in good agreement with the available numerical and experimental data (Table 2). C p b , C L , C D , C ¯ D , θ s are understood as mean pressure coefficient, mean square lift coefficient, mean square drag coefficient, mean drag coefficient, mean angular coordinate of the lift-off point, average length of the recirculation region. In particular, in the calculations, the Strouhal number is 0.186, while the data of a physical experiment give values of 0.184 [8] and 0.183 [12]. The angular coordinate of the flow separation point is θ s = 108.5 ° . Averaging is carried out over a time interval equal to 150.
The cause of noise generation is the deterministic tone sound with discrete frequencies and turbulent fluctuations in the boundary layer and wake, which are stochastic in nature and have a continuous frequency spectrum. The sound intensity is found from the relation I = p 2 ¯ / ( ρ 0 c 0 ) . The RMS pressure is determined by the relation
p 2 ( x ) = lim T 1 T 0 T p 2 ( x , t ) d t ,
where T = t n t 0 is the averaging period. Under t 0 and t n , we mean the moments of time corresponding to the beginning and end of saving the samples.
RMS pressure distributions obtained using various approaches are shown in Figure 7. The control surface is located at r 0 / R = 1.2 , and the observation point—on a circle with radius r c / D = 250 . Symbols ◻ correspond to simulation results, dash-dotted line—experimental data [12] at Re = 3000 , double dash-dotted line corresponds to experimental data [12] at Re = 5100 . Good agreement between the calculated and experimental data takes place at 30 < L c / D < 150 .

8. Results and Discussion

The structure of the flow behind a transversely streamlined cylinder in a stationary external flow has been well studied. In the case of a transverse flow around a cylinder, a boundary layer is formed on its surface, the thickness of which gradually increases upstream. The main determining parameters of this layer are the Reynolds number and the turbulence of the oncoming flow. The separation of the boundary layer and the formation of vortices in the wake of cylindrical bodies are periodic processes. The frequency of separation is characterized by the Strouhal number, which depends on the Reynolds number, the turbulence of the external flow, the degree of obstruction of the channel, and other factors.
The flow that is formed when the boundary cylinder is separated from the surface loses its stability and interacts with the pipe walls. As a result, self-sustaining pressure oscillations arise, which lead to the presence of intense peaks in the spectrum of the generated noise. The field of the vorticity modulus and the z-component of the vorticity vector in the symmetry plane are shown in Figure 8. The vortices form an extensive wake, grow downstream, and interact with the pipe walls. The level lines of the z-component of the vorticity vector give an idea of the direction of rotation of the vortices.
The distribution of overpressure on the channel walls is shown in Figure 9. The pressure is distributed in a significantly non-uniform way, and therefore, one should expect the occurrence of oscillations of the elastic walls of the pipe.
One of the criteria used to identify eddy flows is the Q-criterion
Q = 1 2 | Ω | 2 | S | 2 ,
where | S | = [ trace ( S S ) ] 1 / 2 , | Ω | = [ trace ( Ω Ω ) ] 1 / 2 . The tensor invariant is
Q = 1 2 λ 1 + λ 2 + λ 3 ,
where λ i are the eigenvalues of the velocity gradient tensor ( i = 1 , 2 , 3 ). A vortex is defined as a flow region in which the inequality Q > 0 is satisfied (the flow region in which the norm of the vorticity tensor exceeds the norm of the strain rate tensor). The connection between the pressure and the Q criterion is established using the Poisson equation for the pressure 2 p = 2 ρ Q . The condition Q > 0 does not necessarily correspond to the pressure maximum (for Q > 0 , the pressure maximum is observed at the boundary).
To identify eddies, the λ 2 -criterion is used, based on the decomposition of the velocity gradient tensor into symmetric (strain rate tensor) and antisymmetric (rotation tensor) tensors [35]. The S 2 + Ω 2 tensor is symmetric and has real eigenvalues ( λ 1 λ 2 λ 3 ), two of which are negative. The vortex flow region is defined as the region in which
λ 2 S 2 + Ω 2 < 0 ,
where λ 2 ( A ) is an eigenvalue of the symmetric tensor A. Under adiabatic conditions, this criterion guarantees an instantaneous pressure minimum in a two-dimensional flow [36].
The level surfaces of the Q and λ 2 criteria, which are different combinations of the velocity gradient tensor invariants [37], colored according to the value of the velocity modulus, make it possible to visualize the vortex structure of the flow (Figure 10). The figure clearly shows both the three-dimensional flow structure and the resulting large-scale vortex component of the Karman street.
The presented fields of physical quantities illustrate the general picture of the separated flow near the cylinder. As quantitative estimates of the accuracy of solving this problem, the most interesting are the spectra of pressure fluctuations on the cylinder surface and acoustic loads on the pipe walls. At each point on the wall, as well as at some selected points, the Fourier transform of the data is performed, expressing the dependence of pressure on time. The time dependence of overpressure on the pipe wall is shown in Figure 11 at x = 0.1 m. Dotted lines show pressure levels corresponding to ± 200 Pa.
The Fourier transform of the pressure signal at the same point is shown in Figure 12. Line 1 shows the results of calculations, line 2—smoothed spectrum, and columns—results of calculations using the 1 / 12 filter. A peak with a frequency corresponding to the fundamental tone of the pulsations stands out well.
The distribution of the spectrum over all points on the channel wall for a certain frequency makes it possible to judge the places of the most intense impact from the flow on the walls at a given frequency (Figure 13).
The power spectral densities of near-wall pressure fluctuations on the streamlined surface and on the channel walls have clear discrete peaks that correspond to the nature of the vortex motion over the surface under study. Discrete peaks of the spectral levels of near-wall pressure fluctuations, which have the most energy-intensive main frequency harmonic, have higher-order harmonics, as well as sub-harmonics, due to the nonlinear interaction of vortex structures (sources of pseudo-sonic fluctuations) with each other and with the surface in a stream. This occurs due to the merging of the vortex structures dominating here with each other (generation of sub-harmonics) and their destruction in accordance with the mechanism of the cascade process of vortex transformation.
As a result of the calculations, the displacement amplitudes were obtained depending on the frequency. The dependence of the displacement amplitude at the point ( 0.5 , 0 , 0 ) is shown in Figure 14. There is a clear peak corresponding to the frequency of the fundamental tone.

9. Conclusions

A mathematical model and computational algorithm to perform coupled simulations of unsteady processes in the circular channel with a round cylinder were developed and validated. At the first stage, the velocity and pressure fields are calculated with delayed detached eddy simulation. At the second stage, the problem of calculating the stress–strain state of the pipe is solved using the finite element method.
Stand-alone CFD simulation of unsteady viscous flowfield around a circular cylinder were performed, and results of CFD calculations were validated against experimental and computational data from the literature. The values of Strouhal numbers and amplitudes of oscillations of the aerodynamic coefficients are determined with the fast Fourier transformation technique.
When fluid flow separates from the cylinder and flow loses its stability, self-sustaining pressure oscillations arise, which lead to the presence of intense peaks in the spectrum of the generated noise. The spectrum of pressure oscillations on the cylinder surface and wall of the pipe have clear discrete peaks. The behavior of spectral dependencies is generally preserved if the flow velocity (Reynolds number) increases.

Funding

The research is partially funded by the Ministry of Science and Higher Education of the Russian Federation as part of World-class Research Center program: Advanced Digital Technologies (contract No. 075-15-2020-903 dated 16 November 2020).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Kudashev, E.B. Experimental studies of flow noise on a pop-up device. Acoust. J. 2005, 51, 488–499. [Google Scholar]
  2. Snarski, S.R.; Lueptow, R.M. Wall pressure and coherent structures in a turbulent boundary layer on a cylinder in a axial flow. J. Fluid Mech. 1995, 286, 137–171. [Google Scholar] [CrossRef]
  3. Khasanov, H.A. Acoustic oscillations of a gas near nested thin-walled cylindrical obstacles in a channel. J. Appl. Mech. Tech. Phys. 2011, 4, 91–99. [Google Scholar] [CrossRef]
  4. Mazo, A.B.; Morenko, I.V. Numerical simulation of viscous separated flow around a rotating circular cylinder. J. Eng. Phys. Thermophys. 2006, 79, 75–81. [Google Scholar] [CrossRef]
  5. Okhotnikov, D.I.; Molochnikov, V.M.; Mazo, A.B.; Malyukov, A.V.; Goltsman, A.E.; Saushin, I.I. Viscous flow in the wake of a circular cylinder near the wall at moderate Reynolds numbers. Thermophys. Aeromech. 2017, 24, 897–907. [Google Scholar] [CrossRef]
  6. Isaev, S.A.; Sudakov, A.G.; Zhukova, Y.V.; Usachov, A.E. Simulation of drag reduction and removal of alternating load on a circular cylinder due to the throttling effect. J. Eng. Phys. Thermophys. 2014, 87, 904–907. [Google Scholar] [CrossRef]
  7. Isaev, S.A.; Zhukova, Y.V.; Popov, I.A.; Sudakov, A.G. Vortex intensification of heat transfer in laminar flow around a circular and elliptical cylinders with air and M-20 oil. J. Eng. Phys. Thermophys. 2018, 91, 664–672. [Google Scholar]
  8. Williamson, C.H.K. Oblique and parallel modes of vortex shedding in the wake of a circular cylinder at low Reynolds numbers. J. Fluid Mech. 1989, 206, 579–627. [Google Scholar] [CrossRef] [Green Version]
  9. Williamson, C.H.K. Vortex dynamics in the cylinder wake. Annu. Rev. Fluid Mech. 1996, 28, 477–539. [Google Scholar] [CrossRef]
  10. Inoue, O.; Hatakeyama, N. Sound generation by a two-dimensional circular cylinder in a uniform flow. J. Fluid Mech. 2002, 471, 285–314. [Google Scholar] [CrossRef]
  11. Etkin, B.; Korbacher, G.K.; Keefe, R.T. Acoustic radiation from a stationary cylinder in a fluid stream (Aeolian tones). J. Acoust. Soc. Am. 1957, 29, 30–36. [Google Scholar] [CrossRef]
  12. Norberg, C. An experimental investigation of the flow around a circular cylinder: Influence of aspect ratio. J. Fluid Mech. 1994, 258, 287–316. [Google Scholar] [CrossRef]
  13. Henderson, R.D. Details of the drag curve near the onset of vortex shedding. Phys. Fluids 1995, 7, 2102–2104. [Google Scholar] [CrossRef] [Green Version]
  14. Zhang, H.-Q.; Fey, U.; Noack, B.R.; Konig, M.; Eckelmann, H. On the transition of the cylinder wake. Phys. Fluids 1995, 7, 779–794. [Google Scholar] [CrossRef] [Green Version]
  15. Fey, U.; Kőnig, M.; Eckelmann, H. A new Strouhal–Reynolds-number relationship for the circular cylinder in the range 47 < Re < 2×105. Phys. Fluids 1998, 10, 1547–1549. [Google Scholar]
  16. Persillon, A.; Braza, M. Physical analysis of the transition to turbulence in the wake of a circular cylinder by three-dimensional Navier–Stokes simulation. J. Fluid Mech. 1998, 365, 23–88. [Google Scholar] [CrossRef]
  17. Ma, X.; Karamanos, G. Dynamics and low-dimensionality of a turbulent wake. J. Fluid Mech. 2000, 410, 29–65. [Google Scholar] [CrossRef] [Green Version]
  18. Posdziech, O.; Grundmann, R. Numerical simulation of the flow around an infinitely long circular cylinder in the transition regime. Theor. Comput. Fluid Dyn. 2001, 15, 121–141. [Google Scholar]
  19. Tremblay, F.; Manhart, M.; Friedrich, R. DNS of flow around a circular cylinder at a subcritical Reynolds number with cartesian grids. In Proceedings of the 8th European Turbulence Conference, Barcelona, Spain, 27–30 June 2002; pp. 659–662. [Google Scholar]
  20. Qu, L.; Norberg, C.; Davidson, L.; Peng, S.-H.; Wang, F. Quantitative numerical analysis of flow past a circular cylinder at Reynolds number between 50 and 200. J. Fluids Struct. 2013, 39, 347–370. [Google Scholar] [CrossRef] [Green Version]
  21. Choi, W.-S.; Choi, Y.; Hong, S.-Y.; Song, J.-S.; Kwon, H.-W.; Jung, C.-M. Turbulence-induced noise of a submerged cylinder using a permeable FW–H method. Int. J. Nav. Archit. Ocean. Eng. 2016, 8, 235–242. [Google Scholar] [CrossRef] [Green Version]
  22. Koval, K.A.; Sukhorukov, A.L.; Chernyshev, I.A. Verification results of the numerical method for calculating the hydrodynamic and hydroacoustic characteristics of a fin propulsor. Fundam. Appl. Hydrophys. 2016, 9, 60–72. [Google Scholar]
  23. Liu, X.; Thompson, D.J.; Hu, Z. Numerical investigation of aerodynamic noise generated by circular cylinders in cross-flow at Reynolds numbers in the upper subcritical and critical regimes. Int. J. Aeroacoust. 2019, 18, 470–495. [Google Scholar] [CrossRef]
  24. Kiselev, M.G.; Nifagin, V.A.; Drozdov, A.V. On one model of interaction of a vibroacoustic system with a rigid base. J. Eng. Phys. Thermophys. 1958, 77, 20–27. [Google Scholar]
  25. Valles, N.G.; Kaplunov, R.M. Mathematical modeling of non-stationary separated flow past a lattice of tubes of heat exchangers. J. Eng. Phys. Thermophys. 1987, 53, 57–64. [Google Scholar]
  26. Khabbouchi, I.; Guellouz, M.S.; Nasrallah, S.B. A study of the effect of the jet-like flow on the near wake behind a circular cylinder close to a plane wall. Exp. Therm. Fluid Sci. 2013, 44, 285–300. [Google Scholar] [CrossRef]
  27. Armellini, A.; Casarsa, L.; Giannattasio, P. Separated flow structures around a cylindrical obstacle in a narrow channel. Exp. Therm. Fluid Sci. 2009, 33, 604–619. [Google Scholar] [CrossRef]
  28. Palkin, E.V.; Mulyadzhanov, R.I. Low-frequency pulsations in the problem of fluid flow around a cylinder in a narrow gap at Re = 3750. Sib. Phys. J. 2017, 12, 43–49. [Google Scholar]
  29. Voskoboynik, V.A.; Grinchenko, V.V.; Makarenkov, A.P. Pseudosound behind an obstacle on a longitudinally streamlined cylinder. Acoust. Bull. 2002, 5, 22–36. [Google Scholar]
  30. Voropaev, G.A.; Voskoboinik, A.V.; Voskoboinik, V.A.; Grinchenko, V.V.; Isaev, S.A.; Rozumnyuk, N.V. Sources of pseudosonic pressure fluctuations in the flow around a spherical dimple. Acoust. Bull. 2008, 11, 27–49. [Google Scholar]
  31. Martin, V.; Cummings, A.; Gronier, C. Discrimination of coupled structural/acoustic duct modes by active control: Principles and experimental results. J. Sound Vib. 2004, 274, 583–603. [Google Scholar] [CrossRef]
  32. Shinde, V.; Marcel, T.; Hoarau, Y.; Deloze, T.; Harran, G.; Baj, F.; Cardolaccia, J.; Magnaud, J.-P.; Longatte, E.; Braza, M. Numerical simulation of the fluid—Structure interaction in a tube array under cross flow at moderate and high Reynolds number. J. Fluids Struct. 2014, 47, 99–113. [Google Scholar] [CrossRef] [Green Version]
  33. Volkov, K.N.; Emelaynov, V.N. Gas Flows Induced by Fluid Injection in Energy Systems; Publishing House of Physical and Mathematical Literature: Moscow, Russia, 2011. [Google Scholar]
  34. Volkov, K.N. Simulation of coupled heat transfer on unstructured meshes. Comput. Math. Math. 2013, 53, 396–416. [Google Scholar]
  35. Nikushchenko, D.; Pavlovsky, V. Fluid motion equations in tensor form. In Advances on Tensor Analysis and their Applications; IntechOpen: London, UK, 2020. [Google Scholar]
  36. Haller, G. An objective definition of a vortex. J. Fluid Mech. 2005, 525, 1–26. [Google Scholar] [CrossRef]
  37. Volkov, K.N.; Zapryagaev, V.I.; Emelyanov, V.N.; Gubanov, D.A.; Kavun, I.N.; Kiselev, N.P.; Teterina, I.V.; Yakovchuk, M.S. Visualization of Data of Physical and Mathematical Simulation in Gas Dynamics; Publishing House of Physical and Mathematical Literature: Moscow, Russia, 2018. [Google Scholar]
Figure 1. Scheme of the computational domain in the plane of symmetry (dimensions are given in millimeters).
Figure 1. Scheme of the computational domain in the plane of symmetry (dimensions are given in millimeters).
Fluids 07 00382 g001
Figure 2. A block-structured mesh.
Figure 2. A block-structured mesh.
Fluids 07 00382 g002
Figure 3. Change in pressure over time for Re = 150 .
Figure 3. Change in pressure over time for Re = 150 .
Fluids 07 00382 g003
Figure 4. Changes in the lift coefficient (dashed line) and drag coefficient (solid line) over time (a), and the distribution of the average pressure coefficient over the cylinder surface (b) for Re = 150 .
Figure 4. Changes in the lift coefficient (dashed line) and drag coefficient (solid line) over time (a), and the distribution of the average pressure coefficient over the cylinder surface (b) for Re = 150 .
Fluids 07 00382 g004
Figure 5. Comparison of calculation results (symbols •) of the Strouhal number with experimental data [15] (solid line).
Figure 5. Comparison of calculation results (symbols •) of the Strouhal number with experimental data [15] (solid line).
Fluids 07 00382 g005
Figure 6. Instantaneous pressure distribution around a cylinder for Re = 150 .
Figure 6. Instantaneous pressure distribution around a cylinder for Re = 150 .
Fluids 07 00382 g006
Figure 7. RMS pressure distributions depending on the length of the cylinder at θ = 90 ° for Re = 3900 .
Figure 7. RMS pressure distributions depending on the length of the cylinder at θ = 90 ° for Re = 3900 .
Fluids 07 00382 g007
Figure 8. Contours of the modulus of the vorticity vector (a) and the z-component of the vorticity vector (b) in the symmetry plane for Re c = 39 , 809 .
Figure 8. Contours of the modulus of the vorticity vector (a) and the z-component of the vorticity vector (b) in the symmetry plane for Re c = 39 , 809 .
Fluids 07 00382 g008
Figure 9. Overpressure field on pipe walls for Re c = 39,809 .
Figure 9. Overpressure field on pipe walls for Re c = 39,809 .
Fluids 07 00382 g009
Figure 10. Contours of Q-criterion (a) and λ 2 -criterion (b) for Re c = 39 , 809 .
Figure 10. Contours of Q-criterion (a) and λ 2 -criterion (b) for Re c = 39 , 809 .
Fluids 07 00382 g010
Figure 11. Time dependence of pressure on the wall at point ( 0.1 , 0 , 0 ) .
Figure 11. Time dependence of pressure on the wall at point ( 0.1 , 0 , 0 ) .
Fluids 07 00382 g011
Figure 12. Spectrum of pressure fluctuations on the wall at point ( 0.1 , 0 , 0 ) .
Figure 12. Spectrum of pressure fluctuations on the wall at point ( 0.1 , 0 , 0 ) .
Fluids 07 00382 g012
Figure 13. Spectrum of pressure fluctuations on the walls at f = 167.5 Hz.
Figure 13. Spectrum of pressure fluctuations on the walls at f = 167.5 Hz.
Fluids 07 00382 g013
Figure 14. The dependence of the displacement amplitude on the frequency at the point ( 0.5 , 0 , 0 ) .
Figure 14. The dependence of the displacement amplitude on the frequency at the point ( 0.5 , 0 , 0 ) .
Fluids 07 00382 g014
Table 1. Comparison of calculated data with available numerical and experimental results (for Re = 150 ).
Table 1. Comparison of calculated data with available numerical and experimental results (for Re = 150 ).
SourceSt C ¯ pb C L C D L z / D
CFD0.1860.8500.3421.22420
[18]0.1840.8481.313140
[13]0.1850.8851.33356
[16]0.1810.3381.26816.46
[14]0.1910.4031.41812
[12]0.1830.356
[8]0.1840.856
[10]0.1841.334
[20]0.1840.8460.3551.305
Table 2. Comparison of calculated data with available numerical and experimental results (for Re = 3900 .
Table 2. Comparison of calculated data with available numerical and experimental results (for Re = 3900 .
SourceSt C ¯ pb C L C D C ¯ D θ ¯ s L r / D L z / D
CFD0.2090.8600.0740.0290.96388.31.5720
[17]0.2190.8401.5920
[19]0.2200.9201.0385.71.3018
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Volkov, K. Coupled Simulation of Fluid Flow and Vibro-Acoustic Processes in the Channel with a Circular Cylinder. Fluids 2022, 7, 382. https://doi.org/10.3390/fluids7120382

AMA Style

Volkov K. Coupled Simulation of Fluid Flow and Vibro-Acoustic Processes in the Channel with a Circular Cylinder. Fluids. 2022; 7(12):382. https://doi.org/10.3390/fluids7120382

Chicago/Turabian Style

Volkov, Konstantin. 2022. "Coupled Simulation of Fluid Flow and Vibro-Acoustic Processes in the Channel with a Circular Cylinder" Fluids 7, no. 12: 382. https://doi.org/10.3390/fluids7120382

APA Style

Volkov, K. (2022). Coupled Simulation of Fluid Flow and Vibro-Acoustic Processes in the Channel with a Circular Cylinder. Fluids, 7(12), 382. https://doi.org/10.3390/fluids7120382

Article Metrics

Back to TopTop