Low-speed opening wind tunnel background noise prediction method

Abstract

The invention discloses a low-speed opening wind tunnel background noise prediction method, and relates to the field of aerodynamic noise research, and the method comprises the steps: 1, obtaining the spectral characteristics and frequency of turbulence pulsation in wind tunnel incoming flow, and obtaining the intensity and frequency of vortex waves based on the spectral characteristics and frequency of turbulence pulsation in wind tunnel incoming flow; 2, obtaining the shear layer thickness of a wind tunnel core test section, and obtaining the vortex core radius of the vortex based on the shear layer thickness of the wind tunnel core test section; step 3, calculating vortex wave propagation velocity, vortex convection velocity and vortex intensity based on wind tunnel incoming flow velocity; and 4, performing numerical calculation based on the data obtained in the steps 1 to 3 to obtain a low-speed opening wind tunnel background noise prediction result. Wherein a control equation used in numerical calculation in the step 4 is an improved control equation ignoring a speed gradient item, and the method can accurately predict aerodynamic noise in vortex-dominated complex flow in engineering practice.

Description

Technical Field

[0001] The invention relates to the field of aerodynamic noise research, and in particular to a method for predicting background noise of a low-speed open wind tunnel. Background Art

[0002] The region where aerodynamic noise is generated and propagated can be divided into the near field and the far field. The far field region is the region where sound waves propagate linearly. In this region, the flow field and the sound field are decoupled, and the linear propagation characteristics such as sound wave refraction and scattering can be directly obtained by solving the linearized Euler equations under the condition of obtaining the background flow field in advance. The near field is the region where aerodynamic noise is generated nonlinearly. In this region, the sound field and the flow field are strongly coupled. As part of the flow, the generation process of sound waves, together with the strongly nonlinear and unsteady flow field motion, is dominated by the compressible Navier-Stokes equations. The nonlinear process in the near field not only generates sound waves, but also vortex waves, while in the far field, the only disturbance component in the flow is sound waves. For most aeroacoustic problems, there is a transition region dominated by vortices between the near field and the far field. In this region, the physical process of vortex wave scattering and sound generation will occur, which is widely present in the field of aerospace engineering. For example: the scattering sound of the vortex waves generated in the upstream fan wake in a turbofan engine in the downstream shear layer; the scattering sound of the interaction between the vortex waves in the upstream blade wake and the blade tip vortex of the downstream blade in the counter-rotating propeller flow, etc. The physical process of vortex wave scattering sound is coupled with the nonlinear generation of near-field noise and the far-field sound propagation process, and plays a bridging role in the process from noise generation to outward propagation. In order to further understand the mechanical mechanism and propagation law of noise generation, it is necessary to deeply study the physical process of vortex wave scattering sound in the vortex-dominated area between the near field and the far field.

[0003] However, the existing numerical models of aerodynamic self-noise only explain the mechanism by which the vortex's own motion generates noise in low Mach number flows, but do not clearly give the process by which vortex waves scatter and generate sound in vortex flows, which hinders further quantitative research on the vortex wave scattering and generating sound process using numerical models. As a result, it is impossible to accurately predict the aerodynamic noise in complex flows dominated by vortices in engineering practice. For example, the traditional background noise prediction method for low-speed open wind tunnels only explains the mechanism by which the vortex's own motion generates noise in low Mach number flows, but does not clearly give the process by which vortex waves scatter and generate sound in vortex flows. As a result, the current background noise prediction method for low-speed open wind tunnels has the problem of inaccurate prediction. Summary of the invention

[0004] The purpose of the present invention is to accurately predict the aerodynamic noise in complex flows dominated by vortices in engineering practice, and to achieve accurate prediction of the background noise of a low-speed open wind tunnel.

[0005] To achieve the above object, the present invention provides a method for predicting the background noise of a low-speed open wind tunnel, and the method includes: Step 1: Obtain the spectral characteristics and frequencies of the turbulent pulsations in the wind tunnel inflow, and obtain the intensity and frequencies of the vortex waves based on the spectral characteristics and frequencies of the turbulent pulsations in the wind tunnel inflow; Step 2: Obtain the shear layer thickness of the core test section of the wind tunnel, and obtain the vortex core radius of the vortices based on the shear layer thickness of the core test section of the wind tunnel; Step 3: Calculate and obtain the vortex wave propagation speed, the vortex convection speed, and the vortex intensity based on the wind tunnel inflow velocity; Step 4: Perform numerical calculations based on the data obtained in Steps 1 to 3 to obtain the prediction results of the background noise of the low-speed open wind tunnel; Among them, the governing equation used in the numerical calculation in Step 4 is the improved governing equation with the velocity gradient term ignored.

[0006] Among them, the applicant's research found that in the problem of acoustic scattering in globally unstable flows, the incident vortex waves can further excite exponentially growing unstable waves, generating strong non-physical perturbations in the flow field and contaminating the physical information of the scattering field. That is, it is found that the velocity gradient will generate unstable waves. In view of this situation, in the present invention, the velocity gradient in the linearized equation is ignored, and the governing equation is corrected to ensure the stability of the solution, so as to accurately obtain the scattering field, and further accurately characterize the complex non-linear dynamic system of the spatio-temporal evolution of the vortex waves scattering and generating sound in typical vortex flows, and can accurately predict the aerodynamic noise in complex flows dominated by vortices in engineering practice, realizing the ability to accurately predict the background noise of low-speed open wind tunnels.

[0007] Preferably, Step 4 specifically includes: Step 4.1: Perform non-dimensionalization processing on all the parameters used in the numerical calculation; Step 4.2: Assign values to the parameters used in the numerical calculation based on the data obtained in Steps 1 to 3; Step 4.3: Set a sponge layer around the physical domain and obtain the damping term based on the sponge layer; Step 4.4: Construct and obtain the improved governing equation with the velocity gradient term ignored and the damping term added, and obtain the background flow expression; Step 4.5: Set the computational domain and the physical domain, and perform mesh division; Step 4.6: Set the initial perturbation field; Step 4.7: Set the boundary conditions; Step 4.8: Set the spatial discretization; Step 4.9: Set the time advancement; Step 4.10: Set the numerical filtering; Step 4.11: Set the calculation duration; Step 4.12: Calculate and obtain the prediction result of the background noise of the low-speed open wind tunnel.

[0008] Among them, through the above method, the prediction result of the background noise of the low-speed open wind tunnel can be accurately obtained.

[0009] Preferably, the specific steps of Step 4.12 include: Start the calculation to obtain the calculation result; Subtract the disturbance amount corresponding to the vortex wave at the current moment from the calculation result to obtain the prediction result of the background noise of the low-speed open wind tunnel.

[0010] Preferably, the improved control equation after ignoring the velocity gradient term and adding the damping term is: ; ; ; ; ; Among them, ρ ′, u ′, v ′ and p ′ are respectively the density of the disturbance in the flow, x directional velocity, y directional velocity and pressure; γ represents the adiabatic index, U c is the vortex convection velocity, is the damping term, σ ( x , y ) is the sponge layer strength, is the reference solution, ρ 0, u 0, v 0 and p 0 are respectively the density of the background flow, x directional velocity, y directional velocity and pressure, u ∞ is the vortex wave propagation velocity, x is the spatial coordinate in the mainstream direction, y is the transverse spatial coordinate perpendicular to the mainstream direction.

[0011] Preferably, the spatial discretization adopts a sixth-order enhanced optimization format, the time advancement adopts the fourth-order Runge-Kutta method for time advancement, and the tenth-order central compact filtering method is used for numerical filtering.

[0012] Preferably, the specific setting of the calculation domain in Step 4.5 includes: in the rectangular coordinate system ( x, y) Calculate under the condition that the calculation uses a rectangular calculation domain. At the initial moment, the vortex core is located at the geometric center of the rectangular calculation domain. Introduce a plane vortex wave propagating in the positive direction from infinity through the calculation boundary. x The plane vortex wave propagating in the positive direction.

[0013] Preferably, the plane vortex wave propagating in the positive direction from infinity x has the following expression: ; where are respectively the density of the vortex wave perturbation, x the velocity in the y direction, the velocity in the ε direction and the pressure; c ∞ is the speed of sound at infinity, f 0 is the vortex wave frequency, u ∞ is the vortex wave propagation speed, and the vortex wave wavelength λ = u ∞ / f 0, t is the time, x is the spatial coordinate in the mainstream direction, y is the transverse spatial coordinate perpendicular to the mainstream direction.

[0014] Preferably, the background flow expression in step 4.4 is: ; where ρ 0, u 0, v 0 and p 0 are respectively the density of the background flow, x the velocity in the y direction, the velocity in the u ∞ direction and the pressure, U v is the maximum tangential velocity, x is the spatial coordinate in the mainstream direction, y is the transverse spatial coordinate perpendicular to the mainstream direction, R is the vortex core radius, M v is the vortex strength, γ represents the adiabatic index, ρ ∞ and c ∞ are respectively the density and the speed of sound at infinity, , U c is the vortex convection velocity,t is time, and r represents the distance from the spatial coordinates (x, y) to the vortex center of the vortex.

[0015] Preferably, the calculation method for obtaining the prediction result of the background noise of the low-speed open wind tunnel by subtracting the disturbance amount corresponding to the vortex wave at the current moment from the calculation result is as follows: ; Where is the density of the acoustic wave disturbance, x the velocity in the direction, y the velocity in the direction, and the pressure; ρ ′, u ′, v ′, and p ′ are respectively the density of the disturbance in the flow, x the velocity in the direction, y the velocity in the direction, and the pressure; is the density of the vortex wave disturbance, x the velocity in the direction, y the velocity in the direction, and the pressure, x is the spatial coordinate in the mainstream direction, y is the transverse spatial coordinate perpendicular to the mainstream direction.

[0016] Preferably, the physical domain for calculation is X 1⩽ x ⩽ X 2, Y 1⩽ y ⩽ Y 2. A sponge layer with a thickness of 2.5 λ is set around the physical domain, λ is the wavelength of the vortex wave, X1 and X2 respectively represent the left boundary and the right boundary of the calculation domain, and Y1 and Y2 respectively represent the lower boundary and the upper boundary of the calculation domain.

[0017] One or more technical solutions provided by the present invention have at least the following technical effects or advantages: The present invention adopts a numerical calculation method of spatio-temporal analysis based on a high-order accurate and high-resolution method, which can simulate the physical process of vortex wave scattering and sound generation, and study the influence law of key physical parameters on the scattering sound field.

[0018] Combined with the classical theory of aeroacoustics, this method can accurately characterize the complex non-linear dynamic system of the spatio-temporal evolution of vortex wave scattering and sound generation in typical vortex flows, can accurately predict the aeroacoustic noise in complex flows dominated by vortices in engineering practice, realize the accurate prediction of the background noise of the low-speed open wind tunnel. The wind tunnel background noise result obtained by using the method proposed by the present invention is less than 3 dB compared with the measured result, and is equivalent to the direct aeroacoustic noise calculation result to achieve accurate prediction, and the calculation efficiency is increased by more than 90% compared with the direct calculation result, realizing efficient prediction. Description of the Drawings

[0019] The drawings described herein are used to provide a further understanding of the embodiments of the present invention, and constitute a part of the present invention, but do not limit the embodiments of the present invention; Figure 1 It is a schematic flow chart of a method for predicting the background noise of a low-speed open wind tunnel; Figure 2 Schematic diagram of the numerical model in the present invention; Figure 3 It is a schematic diagram of the sound pressure result. Detailed Embodiments

[0020] In order to more clearly understand the above objects, features and advantages of the present invention, the present invention will be further described in detail below in conjunction with the drawings and specific embodiments. It should be noted that, without conflict, the embodiments of the present invention and the features in the embodiments can be combined with each other.

[0021] Many specific details are set forth in the following description in order to provide a thorough understanding of the present invention. However, the present invention may be implemented in other ways different from those described herein, and therefore, the protection scope of the present invention is not limited by the specific embodiments disclosed below.

[0022] Embodiment 1; Please refer to Figure 1 , Figure 1 It is a schematic flow chart of a method for predicting the background noise of a low-speed open wind tunnel. Embodiment 1 of the present invention provides a method for predicting the background noise of a low-speed open wind tunnel, and the method includes: Step 1: Obtain the spectral characteristics and frequencies of the turbulent fluctuations in the wind tunnel inflow, and obtain the intensity and frequency of the vortex waves based on the spectral characteristics and frequencies of the turbulent fluctuations in the wind tunnel inflow; Among them, the intensity of the turbulent fluctuations in the wind tunnel inflow is the intensity of the vortex waves. Perform one-third octave analysis on the spectral characteristics of the turbulent fluctuations, and the center frequency of the frequency-dominant band obtained is the frequency of the vortex waves.

[0023] Step 2: Obtain the shear layer thickness of the core test section of the wind tunnel, and obtain the vortex core radius of the vortex based on the shear layer thickness of the core test section of the wind tunnel; Among them, the vortex core radius is equal to half of the shear layer thickness of the core test section.

[0024] Step 3: Calculate and obtain the vortex wave propagation speed, vortex convection speed and vortex intensity based on the wind tunnel inflow velocity; Among them, the wind tunnel inflow velocity is the vortex wave propagation speed; the vortex convection speed is 60% of the wind tunnel inflow velocity; the vortex intensity is the ratio of the wind tunnel inflow velocity to the sound speed at infinity.

[0025] Step 4: Perform numerical calculations based on the data obtained in Steps 1 to 3 to obtain the prediction result of the background noise of the low-speed open wind tunnel; Step 5: Perform three-dimensional correction on the prediction result of the background noise of the low-speed open wind tunnel to obtain the actual background noise result of the low-speed open wind tunnel.

[0026] Among them, the control equation used in the numerical calculation in Step 4 is the improved control equation that ignores the velocity gradient term.

[0027] The following is a detailed introduction to this method: The numerical model in the present invention can be applied to the prediction of aerodynamic noise in low Mach number (Mach number less than 0.3) flows dominated by large-scale vortex structures. Taking the prediction of the background noise of a low-speed open wind tunnel as an example, the application of the present invention in the prediction of actual engineering noise is illustrated below.

[0028] The background noise of a low-speed open wind tunnel is mainly generated by the interaction between the turbulent pulsation manifested as vortex waves upstream and the large-scale vortex structures in the downstream shear layer. The specific prediction steps using the present invention are as follows: By experimental measurement or theoretical analysis, obtain the spectral characteristics and intensity of the turbulent pulsation in the wind tunnel inflow, and then the intensity of the vortex waves can be obtained ε and frequency f 0; Then, by experimental measurement or numerical simulation, obtain the shear layer thickness, and then the vortex core radius of the vortex can be obtained R ; Obtain the vortex wave propagation speed u ∞ , vortex convection speed U c and vortex intensity M v ; Directly adopt the numerical calculation method proposed by the present invention to carry out numerical calculations to obtain the far-field noise; Perform three-dimensional correction on the calculation result to obtain the actual far-field noise result.

[0029] Among them, the three-dimensional correction method can be: perform Fourier transform on the calculation result, is the sound pressure corresponding to the frequency f , and the corrected sound pressure can be expressed as: ; The background noise result of the wind tunnel obtained by using the numerical model proposed by the present invention is less than 3 dB compared with the measured result, which is comparable to the direct aerodynamic noise calculation result, and the calculation efficiency is improved by more than 90% compared with the direct calculation result.

[0030] The numerical calculation of the vortex-wave scattering sound generation process proposed by the present invention includes a physical model and a calculation method. The physical model includes: The physical models used in the present invention mainly include the vortex-wave model and the convective isentropic vortex model. The specific derivation processes of the two models are given below.

[0031] Vortex-wave model: The vortex-wave model adopted by the present invention is a form of perturbation propagating with a uniform background flow. Without loss of generality, the two-dimensional uniform background flow adopted in the present invention satisfies the density ρ ∞ , x the velocity in the u ∞ , y the velocity in the v ∞ direction, and the pressure p ∞ are all constants, and v ∞ = 0. In this case, the specific form of the linearized Euler equation is: ; (1) where ρ ′, u ′, v ′, and p ′ are respectively the density of the perturbation in the flow, x the velocity in the y direction, the velocity in the γ direction, and the pressure; t represents the adiabatic index, generally 1.4 in air, x is the time, y is the spatial coordinate in the mainstream direction,

[0032] Assume that the above equation set (1) has a traveling wave solution in the following form: ; (2) where q 1, q 2, q 3, and q 4 are dimensionless constants, ω is the frequency, i is the imaginary unit, k 1 and k 2 are respectively x the wave numbers in the y direction and the c ∞ direction, and q 1, q2、 q 3 and q 4 for the homogeneous linear equations: ; (3) It should be noted that the relationship of ideal gas is used here . The condition for the homogeneous linear equations (3) to have non - zero solutions is that the coefficient determinant is equal to 0, that is: ; (4) Solving the above equation (4) gives: ; (5) That is or . At this time, there are the following three traveling - wave solutions for the equations (1), where corresponds to sound waves, corresponds to entropy waves and vorticity waves, and the specific discussion is as follows.

[0033] The propagation speed of sound waves is , where θ is the angle between the sound - wave propagation direction and the flow direction, is the dispersion relation of sound waves propagating in a uniform background flow with velocity u 0, and the acoustic perturbation quantity at this time is: ; (6) where ε is an infinitesimal quantity.

[0034] Vorticity waves and entropy waves propagate with the background flow, and the propagation speed is u 0, and the expression of vorticity waves is: ; (7) The expression of entropy waves is: ; (8) Sound waves and entropy waves are longitudinal waves, vorticity waves are transverse waves, sound waves propagate at the speed of sound, and entropy waves and vorticity waves propagate at the speed of the background flow.

[0035] The vorticity - wave model specifically used in the present invention is the vorticity wave propagating along the x direction, and the specific expression is as follows: ; (9) where are respectively the density of vorticity - wave perturbation, x the velocity in the y direction, the velocity in the c ∞ is the speed of sound at infinity, f 0 is the vorticity - wave frequency, u∞ is the vorticity wave propagation velocity, and the vorticity wave length λ = u ∞ / f 0, t is time, x is the spatial coordinate in the mainstream direction, y is the transverse spatial coordinate perpendicular to the mainstream direction; ε is the vorticity wave intensity, , indicating that the vorticity wave perturbation amount is much smaller than the hydrodynamic quantity of the background flow.

[0036] Convective isentropic vortex model: The expression of the classical Taylor vortex in polar coordinates ( r , θ ) is: ; (10) where, r represents the distance from the spatial coordinate (x, y) to the vortex center, θ is the angle between the sound wave propagation direction and the flow direction, u rv is the normal velocity of the vortex, u θv is the tangential velocity of the vortex, ρ v is the density of the vortex, p v is the pressure of the vortex; The Taylor vortex is a compact vortex, and its velocity decays exponentially in the far field. Therefore, the far field can be regarded as a uniform field. ρ ∞ and c ∞ are the density and sound speed at infinity, respectively; R is the vortex core radius. At r = R , the tangential velocity reaches the maximum, denoted as U v , M v is the vortex intensity, defined as .

[0037] From the above expression (10), it is not difficult to see that the classical Taylor vortex is a steady vortex. On this basis, by superimposing a uniform flow with a velocity of x in the u ∞ direction, the convective isentropic vortex, that is, the background flow, can be obtained. Its expression in the rectangular coordinate system ( x , y ) is: ; (11) Among them, ρ 0, u 0, v 0 and p 0 are the density, x directional velocity, y directional velocity and pressure of the background flow respectively; ,( x c , y c ) is the vorticity center position coordinate of the Taylor vortex at t time, x c = x c0 + U c t , y c = y c0 , U c is the vortex convection velocity,( x c0 , y c0 ) is the initial vorticity center position coordinate, x is the spatial coordinate in the mainstream direction, y is the transverse spatial coordinate perpendicular to the mainstream direction.

[0038] When U c = u ∞ , the expression (11) of the convective isentropic vortex is the exact solution of the Euler equation, which is the most classical convective isentropic vortex model. However, considering the interaction between vortices in actual vortex flows, the convective velocity of any vortex should include two parts of contributions: the background flow at this vortex (mostly simple flows such as uniform flow and potential flow) and the flow induced by other vortices at this vortex. For each isolated vortex in the actual flow, its convective velocity is not consistent with the velocity of the local background flow (this has been verified in the wind tunnel jet shear layer: the convective velocity of the vortex in the wind tunnel jet shear is about 65% of the incoming flow velocity). Therefore, in the present invention, a vortex convective velocity different from the uniform background flow is set, that is U c ≠ u ∞ . In this case, the equation satisfied by the expression (11) is: ; (12) Compared with the Euler equation, the above system of equations (12) has some additional source terms on the right side, that is, density, xDirectional velocity y The derivatives of the directional velocity, x pressure in the direction, multiplied by the difference between the convective velocity of the uniform background flow and the vortex, is characterized as the net flux of the primitive variables in the convective direction. The above system of equations (12) is written in conservation form: ; (13) where U is the conserved variable, , F and G are respectively x , y the flux vectors in the direction, , , ρ , u , v , p are respectively density, x and y the velocity in the direction, pressure; the total energy density . S is the source term flux vector, representing the additional flux vector brought to the selected isolated vortex structure in the present invention due to the interaction between vortices. The specific expression is: ; (14) respectively represent the additional mass flux, x the momentum flux in the direction, y the momentum flux in the direction and the energy flux. This is due to the action of other vortices, generating a mass flow in the convective direction with a velocity equal to the difference between the uniform background flow and the convective velocity of the vortex u ∞ -[[]]END]] U c into the mass, and the inflowing mass brings in turn x the momentum in the direction, y the momentum in the direction and the energy.

[0039] In addition, for the convenience of solving the above system of equations (14) can also be written in the following form: ; (15) Calculation method: The following introduces the high-precision finite-difference time-domain method for calculating the scattering of convective isentropic vortices by vortex waves used in the invention, including the governing equations, spatial discretization scheme, time marching scheme, numerical filtering scheme, grid stretching function, sponge layer technique, and non-reflecting boundary conditions, etc.

[0040] Governing equations: Considering that both the vortex waves and acoustic waves involved in the present invention are small perturbation waves, the governing equations for the calculation are selected as the linearized equations obtained by linearizing the non-linear equations (12), (13) or (15).

[0041] In the physical process of vortex wave scattering by convective isentropic vortices, the original physical variables ( ρ , u , v , p ) are composed of the background hydrodynamic quantities ( ρ 0, u 0, v 0, p 0) and the perturbation quantities in the flow ( ρ ′, u ′, v ′, p ′), that is: ; (16) where both the flow field and the background flow field satisfy the non-linear equations (12), (13) or (15) obtained previously.

[0042] Substituting (16) into (12), (13) or (15) and neglecting the higher-order small terms, the linearized equations for the small perturbation quantities of the vortex waves and acoustic waves can be obtained, and its conservation form is: ; (17) where: ; ; ; ; (18) Here, ρ 0, u 0, v 0 and p 0 satisfy the expression (11) of the convective isentropic vortex.

[0043] Compared with the standard linearized Euler equations (LEEs for short), the linearized equation (17) obtained for the vortex wave scattering by convective isentropic vortices has an additional term . Here characterizes the influence of the interaction between vortices on the vortex wave scattering.

[0044] The linearized equation (17) is an ideal control equation for the scattering of convective isentropic vortex pairs by vortex waves in the present invention, but it has its applicable conditions. From the above derivation process, it can be seen that the incident vortex wave must satisfy the small perturbation condition, that is, the amplitude of the perturbation is much smaller than the magnitude of the corresponding hydrodynamic quantity. In addition, the linearized equation only supports the propagation of various different waves such as vortex waves, sound waves, and entropy waves at the same time. These different types of waves are coupled together in a strongly non-uniform flow and are extremely difficult to separate. Even worse, in the problem of acoustic scattering in a globally unstable flow, the incident vortex wave can further excite exponentially growing unstable waves, generating strong non-physical perturbations in the flow field and contaminating the physical information of the scattering field. In view of this situation, in the present invention, the method of neglecting the velocity gradient term correction H in the linearized equation (17) is adopted to ensure the stability of the solution, so as to accurately obtain the scattering field. The corrected H term is in the following form: ; (19) Numerical method: In order to better analyze the various perturbation components in the scattering field, when numerically solving the two-dimensional linearized equation (17), the sixth-order enhanced optimized (Enhanced Optimized, abbreviated as EnhOpt) format is used for spatial discretization, and the classical fourth-order Runge-Kutta (Runge-Kutta, abbreviated as R-K) method is used for time advancement. In order to maintain the stability of the calculation, a tenth-order central compact filtering technique is used to suppress the non-physical high-frequency oscillations generated in the calculation.

[0045] The calculation is carried out in a rectangular coordinate system ( x, y ), and a rectangular calculation domain is adopted. At the initial moment, the vortex center of the vortex is located at the geometric center of the calculation domain. A plane vortex wave propagating in the positive direction along x from infinity is introduced through the calculation boundary, and the specific form can be seen in expression (9).

[0046] The physical domain of the calculation is X 1⩽ x ⩽ X 2, Y 1⩽ y ⩽ Y 2. X1 and X2 represent the left and right boundaries of the calculation domain respectively, and Y1 and Y2 represent the lower and upper boundaries of the calculation domain respectively. A sponge layer with a thickness of δ =2.5 λ is set around the physical domain. λ is the wavelength of the incident vortex wave, that is, in this region, a damping term is added to the control equation (17). σ ( x , y) is the strength of the sponge layer. Among them, the strength of the sponge layer can be obtained by converting the thickness of the sponge layer. U ref is the reference solution. Specifically, at the left boundary, σ ( x ) = 0.5[( x − X 1) / δ 2.5 , x ⩽ X 1; at the right boundary, σ ( x ) = 0.5[( x − X 2) / δ 2.5 , x ⩾ X 2; at the upper boundary, σ ( y ) = 0.5[( y − Y 2) / δ 2.5 , y ⩾ Y 2; at the lower boundary, σ ( y ) = 0.5[( y − Y 1) / δ 2.5 , y ⩽ Y 1. U ref It includes the background flow and the planar vortex wave at the current moment, so that the sponge layer only absorbs the scattered wave without affecting the incident wave entering the calculation area. At the boundary, the physical quantities are divided into vortex background flow, incident vortex wave and scattered wave. Among them, the background flow is given by the expression (11) of the convective isentropic vortex, the incident wave is given by the planar vortex wave expression (9), and the scattered wave part is determined by the radiation boundary condition in the x direction and determined by the radiation boundary condition in the y direction.

[0047] A uniform Cartesian grid is adopted in the calculation. The grid spacing meets the requirements of being able to resolve both the vortex dynamic quantities and the acoustic small perturbation quantities simultaneously. Take ∆ x = ∆ y = min{ λ / 20, R / 10}. Y 1 = −max{4 λc ∞ / u ∞ , 8 R ​​​​}, Y2 = max{4 λc ∞ / u ∞ , 8 R}, X 1 = Y 1 - 10 λU c / u ∞ , X 2 = Y 2 + 10 λU c / u ∞ . Courant number c ∞ ∆ t / ∆ x = 0.2. After the calculation starts for a period of time X 2 - X 1 + ( Y 2 - Y 1) + 4 δ / u ∞ ), it enters a quasi - steady - state scattering process with a definite period. At this time, record the data of 10 incident - wave periods λ / u ∞ . At the beginning stage of the calculation, freeze the flow field to keep it in a steady state. After entering the quasi - periodic state, start the normal calculation. That is, when the calculation time t < X 2 - X 1 + ( Y 2 - Y 1) + 4 δ / u ∞ , let U c = 0.

[0048] Subtract the perturbation quantity corresponding to the vortex wave at the current moment from the calculation result to obtain the acoustic quantity. The specific calculation method is: ; (29) Among them, is the density of the acoustic wave perturbation, x the velocity in the y direction, ρ ′, u ′, v ′ and p ′ are respectively the density of the perturbation in the flow, x the velocity in the y direction and the pressure; The density of the vortex wave perturbation, x the direction velocity, y the direction velocity and the pressure.

[0049] The following introduces and illustrates this method with specific examples and data: Non - dimensionalization: First, all the parameters used in the numerical calculation are non - dimensionalized to facilitate numerical calculation programming. Using the vortex core radius R as the reference quantity of the spatial scale, and the sound speed at infinity c ∞ as the reference quantity of the velocity scale, and the density at infinity ρ ∞ as the reference quantity of the density scale, and the time t is non - dimensionalized with , and the pressure p is non - dimensionalized with .

[0050] Parameter assignment: The numerical model is as Figure 2 shown, Figure 2 which can be understood as a closed space formed by wrapping with a sponge layer, and there is a sound source inside. After non - dimensionalization, the parameters of the given numerical model are as follows: The vortex intensity M v = 0.25, the vortex wave intensity is 0.00001, the frequency of the vortex wave is 0.125, the wavelength of the vortex wave λ = 2, the vortex wave propagates along the x direction and the propagation speed u ∞ = 0.25, the convection direction of the vortex is x direction and the velocity U c = 0.15.

[0051] Control equations: First, the in the vortex reference frame can be expressed as: ; (20) The non - dimensional control equation is: ; (21) where, ; (22) ( ρ 0, u 0, v 0, p 0 ) is the background flow, and its specific form is as follows: ; (23) H The specific form of the item is as follows: ; (24) Sponge layer strength σ ( x , y ) The specific form is: ; (25) Reference solution U ref The specific form is: ; (26) Computational domain and mesh generation: The calculation is carried out in a rectangular coordinate system ( x, y ), using a rectangular computational domain. The size of the computational physical domain is X 1, X 2] × Y 1, Y 2], X 1 = -38, X 2 = 38, Y 1 = -32, Y 2 = 32. A sponge layer with a strength of δ = 10 is set around the physical domain. Cartesian uniform grids are used, and the number of grids in the x and y directions are 481 and 421 respectively.

[0052] Initial perturbation field setting: ; (27) Boundary condition setting: At the boundary, the physical quantity is divided into incident vortex waves and scattered waves. The expression of the incident vortex wave is: ; (28) The scattered wave part is determined by the radiation boundary condition in the x direction, and y the scattered wave part is determined by the radiation boundary condition in the

[0053] Space discretization: Sixth-order enhanced optimization format is used for space discretization.

[0054] Time marching: The fourth-order Runge-Kutta method is used for time marching, and the time step is 0.02.

[0055] Numerical filtering: The tenth-order central compact filtering technique is adopted to suppress the non-physical high-frequency oscillations generated in the calculation.

[0056] Calculation duration and data saving: The total calculation duration is 800. When t > 720, data recording starts.

[0057] Result output: Subtract the perturbation corresponding to the vortex wave at the current moment from the calculation result, and the result is the acoustic quantity. The specific calculation method is ; (29) Where is the density of the acoustic wave perturbation, x the velocity in the y direction, ρ ′, u ′, v ′ and p ′ are respectively the density of the perturbation in the flow, x the velocity in the y direction, Figure 3 The sound pressure result is as shown in Figure 3 where E in

[0058] is the symbol of scientific notation.

[0059] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalent technologies, the present invention is also intended to include these changes and modifications.

Claims

1. A prediction method for the background noise of a low-speed open wind tunnel, characterized in that, The method includes: Step 1: Obtain the spectral characteristics and frequency of turbulent pulsations in the wind tunnel inflow, and obtain the intensity and frequency of vortex waves based on the spectral characteristics and frequency of turbulent pulsations in the wind tunnel inflow; Step 2: Obtain the shear layer thickness of the wind tunnel core test section, and obtain the vortex core radius of the vortex based on the shear layer thickness of the wind tunnel core test section; Step 3: Calculate and obtain the vortex wave propagation speed, vortex convection speed, and vortex intensity based on the wind tunnel inflow velocity; Step 4: Perform numerical calculations based on the data obtained in Steps 1 to 3 to obtain the prediction result of the background noise of the low-speed open wind tunnel; Among them, the control equation used in the numerical calculation in Step 4 is the improved control equation that ignores the velocity gradient term.

2. The method for predicting the background noise of a low-speed open wind tunnel according to claim 1, wherein The specific content of Step 4 includes: Step 4.1: Nondimensionalize all the parameters used in the numerical calculation; Step 4.2: Assign values to the parameters used in the numerical calculation based on the data obtained in Steps 1 to 3; Step 4.3: Set a sponge layer around the physical domain and obtain the damping term based on the sponge layer; Step 4.4: Construct and obtain the improved control equation that ignores the velocity gradient term and adds the damping term, and obtain the background flow expression; Step 4.5: Set the computational domain and physical domain and perform mesh division; Step 4.6: Set the initial perturbation field; Step 4.7: Set the boundary conditions; Step 4.8: Set the spatial discretization; Step 4.9: Set the time advancement; Step 4.10: Set the numerical filtering; Step 4.11: Set the calculation duration; Step 4.12: Calculate and obtain the prediction result of the background noise of the low-speed open wind tunnel.

3. A method for predicting the background noise of a low-speed open wind tunnel according to claim 2, characterized in that The specific content of Step 4.12 includes: Start the calculation to obtain the calculation result; Subtract the perturbation amount corresponding to the vortex wave at the current moment from the calculation result to obtain the prediction result of the background noise of the low-speed open wind tunnel.

4. A method for predicting the background noise of a low-speed open wind tunnel according to claim 1, characterized in that, The improved control equation that ignores the velocity gradient term and adds the damping term is: ； ； ； ； ； Among them, ρ ′, u ′, v ′ and p ′ are respectively the density, x directional velocity, y directional velocity and pressure of the perturbation in the flow; γ represents the adiabatic index, U c is the vortex convection velocity, is the damping term, σ ( x , y ) is the sponge layer strength, is the reference solution, ρ 0, u 0, v 0 and p 0 are respectively the density, x directional velocity, y directional velocity and pressure of the background flow, u ∞ is the vorticity wave propagation velocity, x is the spatial coordinate in the mainstream direction, y is the transverse spatial coordinate perpendicular to the mainstream direction.

5. A method for predicting the background noise of a low-speed open wind tunnel according to claim 2, characterized in that, The spatial discretization adopts a sixth-order enhanced optimization format, the time advancement adopts a fourth-order Runge-Kutta method for time advancement, and a tenth-order central compact filtering method is used for numerical filtering.

6. A method for predicting the background noise of a low-speed open wind tunnel according to claim 2, characterized in that, The specific settings of the computational domain in Step 4.5 include: performing calculations in a rectangular coordinate system ( x,y ), using a rectangular computational domain for the calculations. At the initial moment, the vortex core is located at the geometric center of the rectangular computational domain. A plane vortex wave propagating in the x positive direction from infinity is introduced through the computational boundary. x is the spatial coordinate in the mainstream direction, y is the transverse spatial coordinate perpendicular to the mainstream direction.

7. A method for predicting the background noise of a low-speed open wind tunnel according to claim 6, characterized in that The expression of the plane vortex wave propagating in the x positive direction from infinity is as follows: ； Among them, are respectively the density of the vortex wave perturbation, x the direction velocity, y the direction velocity and the pressure; ε is the vortex wave intensity, c ∞ is the sound speed at infinity, f 0 is the vortex wave frequency, u ∞ is the vortex wave propagation velocity, and the vortex wave wavelength λ = u ∞ / f 0, t is the time.

8. A method for predicting the background noise of a low-speed open wind tunnel according to claim 2, characterized in that, The background flow expression in Step 4.4 is: ； Among them, ρ 0, u 0, v 0 and p 0 are respectively the density, x directional velocity, y directional velocity and pressure of the background flow, u ∞ is the vorticity wave propagation velocity, U v is the maximum tangential velocity, x is the spatial coordinate in the mainstream direction, y is the transverse spatial coordinate perpendicular to the mainstream direction, R is the vorticity core radius, M v is the vorticity intensity, γ represents the adiabatic index, ρ ∞ and c ∞ are respectively the density and the speed of sound at infinity, , U c is the vorticity convection velocity, t is the time, r represents the distance from the spatial coordinate ( x , y ) to the vorticity center of the vortex.

9. A method for predicting the background noise of a low-speed open wind tunnel according to claim 3, characterized in that, The calculation method of subtracting the perturbation amount corresponding to the vortex wave at the current moment from the calculation result to obtain the prediction result of the background noise of the low-speed open wind tunnel is: ； Among them, is the density of the acoustic wave perturbation, x the directional velocity, y the directional velocity, and the pressure; ρ ′, u ′, v ′, and p ′ are respectively the density of the perturbation in the flow, x the directional velocity, y the directional velocity, and the pressure; is the density of the vortex wave perturbation, x the directional velocity, y the directional velocity, and the pressure, x is the spatial coordinate in the mainstream direction, y is the transverse spatial coordinate perpendicular to the mainstream direction.

10. A method for predicting the background noise of a low-speed open wind tunnel according to claim 2, characterized in that, The calculated physical domain is X 1⩽ x ⩽ X 2, Y 1⩽ y ⩽ Y 2, and a sponge layer with a thickness of 2.5 is set around the physical domain λ . λ is the wavelength of the vortex wave, X1 and X2 respectively represent the left and right boundaries of the calculation domain, and Y1 and Y2 respectively represent the lower and upper boundaries of the calculation domain.

Images