Time-domain simulations of sound propagation through screen-induced turbulence

A flow model in combination with a statistical-dynamical turbulence generator and a linearised Euler time-domain model for sound waves were used to simulate the effect of screen-induced turbulence on the noise level in the acoustical shadow of a screen in wind. Instead of simulating a great number of different frozen turbulence realisations, the concept of transient turbulence was successfully tested and applied. This concept is adequate to the time-domain model and reduces the computational demands. Several two-dimensional simulations allowed to isolate the individual effects of wind and screen on the propagation of 500 Hz sound waves over a 4-m high noise barrier. At a distance of 250 m from the source (240 m behind the screen) the sheltering effect of the screen and the refraction effect of the wind are in the order of 6 and 4 dB, respectively. The screen-induced turbulence leads to fluctuations in the noise level with a standard deviation of 1.2 dB and a maximum amplitude of 3 dB. However, the time averaged effect turned out to be in the order of merely 0.2 dB. The effect of the screen-induced turbulence on the average noise level behind the screen is therefore negligible.     

Three-dimensional coupled-mode model and characteristics of low-frequency sound propagation in ocean waveguide with seamount topography

Large-scale topography, such as a seamount, substantially impacts low-frequency sound propagation in an ocean waveguide, limiting the application of low-frequency acoustic detecting techniques. A three-dimensional (3D) coupled-mode model is developed to calculate the acoustic field in an ocean waveguide with seamount topography and analyze the 3D effect. In this model, a correction is introduced in the bottom boundary, theoretically making the acoustic field satisfy the energy conservation. Furthermore, a large azimuth angle calculation range is obtained by using the operator theory and higher-order Padé approximation. Additionally, the model has advantages related to the coupling mode and parabolic equation theory. The couplings corresponding to the effects of range-dependent environment are fully considered, and the numerical implementation is kept feasible. After verifying the accuracy and reliability of the model, low-frequency sound propagation characteristics in the seamount environment are analyzed. The results indicate lateral variability in bathymetry can lead to out-of-plane effects such as the horizontal refraction phenomenon, while the coupling effect tends to restore the abnormal sound field and produces acoustic field diffraction behind the seamount. This model effectively considers the effects of the horizontal refraction and coupling, which are proportional to the scale of the seamount.     

Air-ground interaction in long range propagation of low frequency sound and vibration—field tests and model verification

An extensive program of intermediate and long range impulsive sound propagation field tests have been conducted. The test program and the performed measurements are presented. Particular focus is given on the air-ground interaction and its effect on low frequency sound and vibration propagation. It is found that the pressure wave interaction with the viscoelastic Rayleigh wave in the ground may have a significant effect on the ground impedance and the sound and vibration propagation. This introduces an important mechanism not covered in commonly used ground impedance models. Numerical simulation models have been developed and verified against the test data. The ground impedance does not only effect the sound pressure propagation. If either acoustically induced ground vibration, or ground to building transmitted vibration, is to be considered, the acousto-seismic impedance has a dramatic effect on the level of ground vibration induced by a given sound pressure. For a site where Rayleigh wave interaction appears at the dominant frequencies of the sound pressure, the ground vibration may be greater than a factor 100 (40 dB) than at a site with ground conditions not making the interaction happen.     

深海海底山环境下声传播水平折射效应研究

声波在深海海底山环境中传播时,海底山会对声传播产生重要影响.2016年在南海深海进行了一次海底山环境下的声传播实验,观测到了由海底山引起的三维声传播效应,本文利用BELLHOP射线理论解释了海底山环境下的三维声传播机理.结果表明：声波在传播过程中与海底山作用后破坏了深海会聚区结构,导致传播损失增大,在海底山后形成具有明显边界的声水平折射区,利用二维声传播模型无法解释实验现象,海底山后声水平折射区实验测量的声场结构与N×2D模型计算结果存在明显差异,实验的传播损失比N×2D模型计算结果大10 dB.通过三维射线模型分析N×2D模型计算结果与实验结果存在明显差异产生的原因,发现由于声波水平折射作用,部分声线无法到达接收器,使得三维声传播效应对海底山后一定角度范围内声场影响较为明显.因此,深海海底山会引起明显的三维水平折射效应,应在水下目标探测和定位等应用中给予重视.     

Numerical model to calculate microphysical influences on sound wave propagation in forests

There is a growing need for increasingly accurate and reliable numerical models to predict micrometeorological and turbulence conditions in complex sound wave propagation environments. In this paper, a prototype finite-difference computer model is developed to calculate the microphysical influences on sound speed in forested areas. Several numerical tests are conducted to assess model code capabilities using micrometeorological field data collected in June 2006. For the current analysis, the model domain and total number of grid points are greatly increased from earlier reported versions and the finite-difference numerical schemes for advection are modified to permit different inflow and outflow boundaries. Preliminary results for three cases are encouraging. Micrometeorological profiles and calculated fields of sound speed are presented. Some initial approximations of short-range acoustic transmission loss for the experimental test site are also discussed.     

Prediction-step staggered-in-time FDTD: An efficient numerical scheme to solve the linearised equations of fluid dynamics in outdoor sound propagation

The finite-difference time-domain (FDTD) method to solve the linearised equations of fluid dynamics has shown to be very powerful and useful in outdoor sound propagation. Practical applications are however limited due to the large need for computational resources. The numerical discretisation influences computational efficiency to an important degree. In this paper, some possible ways to discretise temporal derivatives are studied. Two obvious ways of time-discretisation namely staggered-in-time (SIT) and a simple collocated-in-time (CIT) scheme are compared to the prediction-step staggered-in-time (PSIT) scheme. The latter is intended to be used for the calculation of sound propagation in the typical low wind speeds encountered in the outdoor environment at low heights above the earth’s surface. It was shown that the PSIT scheme is more stable than the SIT scheme, so practical calculations are possible. Computational efficiency is increased to an important degree compared to the CIT scheme. The numerical accuracy (more precisely the amplitude error) of the PSIT scheme is an important improvement upon SIT. The CIT scheme on the other hand conserves amplitude better. The amplitude error becomes larger with increasing wind speed because of some simplifications during the numerical discretisation. In low wind speeds, the PSIT algorithm can serve as an interesting compromise between numerical accuracy and the required amount of computing power.     

Resolution effects and scaling in numerical simulations of passive scalar mixing in turbulence

The effects of finite grid resolution on the statistics of small scales in direct numerical simulations of turbulent mixing of passive scalars are addressed in this paper. Simulations at up to 2048³ grid points with grid spacing Δx varied from about 2 to 1/2 Batchelor scales (η_B) show that most conclusions on Schmidt number (Sc) dependence from prior work at less stringent resolution remain qualitatively correct, although simulations at resolution Δx≈η_B are preferred and will give adequate results for many important quantities including the scalar dissipation intermittency exponent and structure functions at moderately high orders. For Sc≥1, since η_B=ηSc^−1/2 (where η is the Kolmogorov scale), the requirement Δx≈η_B is more stringent than the corresponding criterion Δx≈η for the velocity field, which is thus well resolved in simulations aimed at high Schmidt number mixing. A simple argument is given to help interpret the effects of Schmidt and Reynolds numbers on trends towards local isotropy and saturation of intermittency at high Schmidt number. The present results also provide evidence for a trend to isotropy at high Reynolds number with fixed Sc=1.0. This is a new observation apparently not detected in less well resolved simulations in the past, and will require further investigation in the future.     

Parameter study of sound propagation between city canyons with a coupled FDTD-PE model

A parameter study is performed for the case of two-dimensional sound propagation from a (source) city canyon to a nearby, identical (receiver) city canyon. Focus was on sound pressure levels, relative to the free field, in the shielded canyon. An accurate and efficient coupled FDTD-PE model was applied, exploiting symmetry of the source and receiver canyon. With the proposed calculation method, simulations were necessary in only half the sound propagation domain. The shielding in the receiver canyon in case of a coherent line source was compared to the shielding by an incoherent line source, by means of sound propagation calculations in a number of 2D cross-sections through source and receiver. It was found that the shielding is rather insensitive to the width-height ratio of the canyons. The presence of diffusely reflecting façades and balconies lead to an important increase in shielding compared to flat façades. Rigid façades yield significantly lower shielding compared to partly reflecting façades. Effects of a moving atmosphere were modeled in detail. Shielding decreases significantly in case of downwind sound propagation when comparing to sound propagation in a non-moving atmosphere. Refraction is the most important effect in the latter. In case of upwind sound propagation, turbulent scattering plays an important role and the shielding is similar to the one of a non-moving atmosphere for the parameters used in this paper. The combination of effects, as is shown by some examples, is in general not a simple addition of the separate effects.     

Vortex identification in the wake of a model wind turbine array

Vortex identification techniques are used to analyse the wake of a 4 × 3 array of model wind turbines. The Q-criterion, Δ-criterion, and λ₂-criterion are applied to particle image velocimetry data gathered fore and aft of the last row centerline turbine. The Q-criterion and λ₂-criterion provide a clear indication of regions where vortical activity exists, while the Δ-criterion is not successful. Galilean decomposition, Reynolds decomposition, vorticity, and swirling strength are used to further understand the location and behaviour of the vortices. The techniques identify and display the high-magnitude vortices in high-shear zones resulting from the blade tips. Using Galilean and Reynolds decomposition, swirling motions are shown encapsuling vortex regions in agreement with the identification criteria. The Galilean decompositions selected are 20% and 50% of a convective velocity of 7 m/s. As the vortices convect downstream, the strength of the vortices decreases in magnitude, particularly in the far wake of the array, to approximately 25% of those present in the near wake. A high level of vortex activity is visualised as a result of the top tip of the wind turbine blade -- the location where the highest vertical entrainment is present. Analysing the full frame set, the Q-criterion, λ₂-criterion, and swirling strength prove comparable, while the Δ-criterion under-performs in regions of high turbulence activity, namely in the back of the turbine. Entraining flow into the turbine canopy interacting with the turbine generates high-magnitude vortices concentrated at the blade tips. The count of vortices decreases when moving from the top tip down to the wall, as well as their strength for each Galilean technique when a non-zero threshold is applied. Vortex sizes in the near wake are found comparable to turbine blade, hub, and mast dimensions. In the far wake, the resulting size of the vortices is approximately 30% of those in the near wake. These vortices increase in velocity as they convect downstream, following the mean velocity behaviour. The lowest magnitude vortices reside at the hub height in the near-wake region, where they convect at nearly half the speed of those at the blade tips.     

Modification of Spalart-Allmaras model with consideration of turbulence energy backscatter using velocity helicity

The correlation between the velocity helicity and the energy backscatter is proved in a DNS case of 256³-grid homogeneous isotropic decaying turbulence. The helicity is then proposed to be employed to improve turbulence models and SGS models. Then Spalart-Allmaras turbulence model (SA) is modified with the helicity to take account of the energy backscatter, which is significant in the region of corner separation in compressors. By comparing the numerical results with experiments, it can be concluded that the modification for SA model with helicity can appropriately represent the energy backscatter, and greatly improves the predictive accuracy for simulating the corner separation flow in compressors.     

自混合干涉效应的一般模型及其仿真

从自混合干涉效应的等效结构出发,建立了该效应的一般模型。该模型包含了多重光反馈和多模激光的自混合干涉效应等情况,并且可以转化为以往所报道的弱光反馈水平的弱光模型和单模条件下的一般激光自混合干涉模型,明确了模型中各个参数的物理意义,对自混合干涉理论的研究提供了很好的理论指导。对单模的几种情况分别进行了数值仿真,仿真结果很好地验证了模型的正确性。     

Nonclassical effects in a highly nonlinear generalized homogeneous Dicke model

An intensity dependent nonlinear coupling model of N two-level atoms (generalized Dicke model) interacting dispersively with a bimodal cavity field via two-photon transitions is investigated in a scenario where the rotating wave approximation is assumed. The model becomes homogeneous in the sense that the spin transition frequency is the same for all atoms and the coupling constants emerging from the collective interactions of the atomic system with the cavity field depend only on the particular radiation field mode. This allows us to represent the Dicke Hamiltonian entirely in terms of the total angular momentum J. It is assumed that, initially, the atomic system and the field are in a disentangled state where the field modes are in Glauber coherent states and the atomic system is a superposition of states |JM〉 (Dicke states). The model is numerically tested against simulations of normal squeezing variance of the field, squeezing factors based on the Heisenberg uncertainty principle, along with the statistical properties of the light leading to the possible production of nonclassical effects, such as degree of second-order coherence in the modes, degree of intermode correlation, as well as violation of the Cauchy–Schwartz inequality. Analytical expression of the total density operator matrix elements at t>0 shows the present nonlinear model to be strongly entangled, which is reflected in the time evolution of the linear entropy, where the superposition states are reduced to statistical mixtures. Thus, the present generalized Dicke model does not preserve the modulus of the Bloch vector. The computations, performed in the weak coupling and strong field limits, were conducted via second-order Dyson perturbative expansion of the time evolution operator matrix elements for the totality of the angular momentum states of the atomic system.     

Improved model of streamer propagation in liquid dielectrics

Breakdown and pre‐breakdown cause extensive damage in dielectric oils used in power systems subjected to high voltages. Consequently, several research works are devoted to these phenomena. Related investigations are made very difficult by the numerous influencing parameters and the large number of characteristics of these phenomena. The main goal of the present work was to establish a new electrical model that describes the propagation of streamers in liquid dielectrics for a point‐plane configuration subjected to high‐voltage fields. The novelty consists in the presence of an inductance in the equivalent electrical circuit, which represents the dielectric during the partial discharges characterized by fast variations of current. Criticisms were first made concerning previous electrical models that show some disadvantages. The developed model consists of a series RLC circuit whose element values change during propagation. In particular, the value of the resistance decreases in anticipation of breakdown. The simulation of the developed model gave satisfactory results, similar to experimental ones. For instance, the influence of electrical and geometrical parameters on the current pulses amplitude as well as on the streamers propagation velocity and the quantity of charge and the energy injected were studied. Good agreement was also obtained between the obtained results and experimental ones, which reinforces the validity of the model.     

Construction of a Langevin model from time series with a periodical correlation function: Application to wind speed data

A Langevin-type equation for stochastic processes with a periodical correlation function is introduced. A procedure of reconstruction of the equation from time series is proposed and verified on simulated data. The method is applied to geophysical time series–hourly time series of wind speed measured in northern Italy–constructing the macroscopic model of the phenomenon.     

A simple model of ultrasound propagation in a cavitating liquid. Part I: Theory, nonlinear attenuation and traveling wave generation

The bubbles involved in sonochemistry and other applications of cavitation oscillate inertially. A correct estimation of the wave attenuation in such bubbly media requires a realistic estimation of the power dissipated by the oscillation of each bubble, by thermal diffusion in the gas and viscous friction in the liquid. Both quantities and calculated numerically for a single inertial bubble driven at 20 kHz, and are found to be several orders of magnitude larger than the linear prediction. Viscous dissipation is found to be the predominant cause of energy loss for bubbles small enough. Then, the classical nonlinear Caflish equations describing the propagation of acoustic waves in a bubbly liquid are recast and simplified conveniently. The main harmonic part of the sound field is found to fulfill a nonlinear Helmholtz equation, where the imaginary part of the squared wave number is directly correlated with the energy lost by a single bubble. For low acoustic driving, linear theory is recovered, but for larger drivings, namely above the Blake threshold, the attenuation coefficient is found to be more than 3 orders of magnitude larger then the linear prediction. A huge attenuation of the wave is thus expected in regions where inertial bubbles are present, which is confirmed by numerical simulations of the nonlinear Helmholtz equation in a 1D standing wave configuration. The expected strong attenuation is not only observed but furthermore, the examination of the phase between the pressure field and its gradient clearly demonstrates that a traveling wave appears in the medium.     

多参数背景场误差模型在散射计资料台风风场反演中的应用

利用散射计资料反演海面风场时, 台风区域普遍存在降雨使得风场反演误差很大, 引入降雨地球物理模型函数 (GMF+Rain) 及多解方案 (MSS), 结合二维变分(2DVAR) 模糊去除思想风速反演误差很大程度减小, 但风向反演误差仍有待进一步改善, 如何进一步减小风向反演误差有待进一步研究. 文章介绍了2DVAR模糊去 除方法的基本思想, 针对背景场误差较大时, 2DVAR模糊去除风向误差较大, 引入包含若干参数的背景场误差模型. 基于台风个例数值试验结果, 着重从理论分析角度讨论各参数关于2DVAR模糊去除效果的敏感性, 进而提出最优参数设置方案以改善风向模糊去除效果. 2006年“摩羯”台风QuikSCAT数据风场反演数值试验结果结合理论分析表明: 引入多参数误差模型, 通过设置粗糙误差概率等于0, 2DVAR风向模糊去除效果明显改善; 同时, 背景场的影响可通过增大背景场误差方差, 减小背景场误差相关尺度和减小粗糙误差概率而减小, 进而减小在背景场误差较大情况下的风向反演误差. 关键词： 台风风场反演 二维变分 多参数误差模型 散射计资料     

Cascade-like and scaling behavior of wind velocity increments in the atmospheric surface layer

By using a large amount of data collected in the atmospheric surface layer, we analyze the probability density functions (PDFs), the probability of return and the moments of wind velocity increments. Results show that the PDFs change from the non-Gaussian long-tailed distributions to Gaussian with the increase of time scales. This is similar to what has been observed and interpreted as an indication of cascade in the fully developed homogeneous and isotropic turbulence. Besides, both the probability of return and the moments are found to be scaling with time scales. We then compare above results with the truncated Lévy flights and the log-normal PDF model. It is found that although both models show the cascade-like behavior in the PDFs and the scaling behavior in the probability of return and the moments under some conditions, they are not good enough for quantitatively describing the random process of wind velocity increments.     

Exact solutions of the Boltzmann equation in the VHP model with removal interaction

The linear and nonlinear Boltzmann equation for very hard particles (VHP) is considered in the case when the collision between two particles may lead not only to elastic scattering, but also to a removal event with the disappearance of the molecules. The extended transport equation is solved for arbitrary initial distributions. The computations are carried out explicitly for a special class of initial distributions and for various removal rates. The results are demonstrated graphically. Finally, source terms fulfilling physically reasonable conditions are introduced into the VHP model, and the time-dependent particle number is calculated.