Numerical analysis of dipole sound source around high speed trains

As the maximum speed of high speed trains increases, the effect of aeroacoustic noise on the sound level on the ground becomes increasingly important. In this paper, the distribution of dipole sound sources at the bogie section of high speed trains is predicted numerically. The three-dimensional unsteady flow around a train is solved by the large eddy simulation technique. The time history of vortices shows that unstable shear layer separation at the leading edge of the bogie section sheds vortices periodically. These vortices travel downstream while growing to finally impinge upon the trailing edge of the section. The wavelength of sound produced by these vortices is large compared to the representative length of the bogie section, so that the source region can be regarded as acoustically compact. Thus a compact Green's function adapted to the shape can be used to determine the sound. By coupling the instantaneous flow properties with the compact Green's function, the distribution of dipole sources is obtained. The results reveal a strong dipole source at the trailing edge of the bogie section where the shape changes greatly and the variation of flow with time is also great. On the other hand, the bottom of the bogie section where the shape does not change, or the leading edge and boundary layer where the variation of flow with time is small, cannot generate a strong dipole source.     

Gyroscopic dynamics of antiferromagnetic vortices in the orthoferrite domain wall

The method of generation of antiferromagnetic vortices on the supersound domain wall in the orthoferrites was proposed. Moving antiferromagnetic vortices were accompanied by the solitary deflection waves. These waves were used for investigation of generation and nonlinear dynamics of the antiferromagnetic vortices on a moving domain wall with the help of two- and three-fold digital high-speed photography and Faraday rotation in the orthoferrites plates cut perpendicular to the optical axis. The full velocity of antiferromagnetic vortex nonlinearly increases and saturates on the spin velocity level c. The vortices with smallest topological charges saturate earlier than with big one. The vortices velocity along the domain wall u increases up to the maximum and goes to the dependence u²+v²=c². Vortex dynamics is quasirelativistic on quasirelativistic domain wall. The theory of gyroscopic force in the domain wall of orthoferrites was elaborated by Zvezdin et al. and was confirmed our earlier experimental results.     

Basic sound generation mechanisms in inviscid vortex interactions at low Mach number

The sound generation mechanisms during finite core vortex interactions at low Mach number are investigated in the present study. The theoretical deductions show clearly that the basic sound generation mechanisms are associated with the vortex core deformation and the vorticity centroid dynamics, independent of the vortex system. Such deductions are substantiated by numerical experiments with the interactions of two-dimensional vortices, vortex pairs and vortex rings. Detailed discussions on the similarities and differences between the sound generation processes of the two-dimensional and axisymmetric vortex systems are given. The relative importance of the two sound generation mechanisms in these vortex systems, their characteristics and interactions, which are hardly found in existing literature, are also examined. The present findings have also generalized and substantiated the previous results of the authors on the topic.     

Quantized vortex stability and interaction in the nonlinear wave equation

The stability and interaction of quantized vortices in the nonlinear wave equation (NLWE) are investigated both numerically and analytically. A review of the reduced dynamic law governing the motion of vortex centers in the NLWE is provided. The second order nonlinear ordinary differential equations for the reduced dynamic law are solved analytically for some special initial data. Using 2D polar coordinates, the transversely highly oscillating far field conditions can be efficiently resolved in the phase space, thus giving rise to an efficient and accurate numerical method for the NLWE with non-zero far field conditions. By applying this numerical method to the NLWE, we study the stability of quantized vortices and find numerically that the quantized vortices with winding number m=±1 are dynamically stable, and resp. |m|>1 are dynamically unstable, in the dynamics of NLWE. We then compare numerically quantized vortex interaction patterns of the NLWE with those from the reduced dynamic law qualitatively and quantitatively. Some conclusive findings are obtained, and discussions on numerical and theoretical results are made to provide further understanding of vortex stability and interactions in the NLWE. Finally, the vortex motion under an inhomogeneous potential in the NLWE is also studied.     

Spectrum of the non-abelian phase in Kitaev’s honeycomb lattice model

The spectral properties of Kitaev’s honeycomb lattice model are investigated both analytically and numerically with the focus on the non-abelian phase of the model. After summarizing the fermionization technique which maps spins into free Majorana fermions, we evaluate the spectrum of sparse vortex configurations and derive the interaction between two vortices as a function of their separation. We consider the effect vortices can have on the fermionic spectrum as well as on the phase transition between the abelian and non-abelian phases. We explicitly demonstrate the 2ⁿ-fold ground state degeneracy in the presence of 2n well separated vortices and the lifting of the degeneracy due to their short-range interactions. The calculations are performed on an infinite lattice. In addition to the analytic treatment, a numerical study of finite size systems is performed which is in exact agreement with the theoretical considerations. The general spectral properties of the non-abelian phase are considered for various finite toroidal systems.     

Analytical and numerical study of acoustic intensity field in irregularly shaped room

The modal expansion method has been used to formulate expressions for real and imaginary parts of the complex sound intensity inside enclosures. Based on theoretical results, the computer program has been developed to simulate the acoustic intensity vector field inside the irregular room whose shape resembles the capital letter L. Calculation results have shown that a low-frequency distribution of the acoustic intensity is strongly influenced by the modal localization and the characteristic objects in the active intensity field are energy vortices and saddle points positioned irregularly inside the room. It was found that for small sound damping the vortex centers lie exactly on the lines corresponding to zeros of the eigenfunction for a dominant mode. An increase in a sound attenuation results in the change of vortex positions and can cause the formation of new vortices. Finally, an influence of the wall impedance on the quantitative relation between the acoustic and reactive intensities was studied and it was concluded that for very small sound damping the behavior of the sound intensity is basically only oscillatory.     

Numerical solutions of the acoustic eigenvalue equation in the rectangular room with arbitrary (uniform) wall impedances

Two numerical procedures for finding the acoustic eigenvalues in the rectangular room with arbitrary (uniform) wall impedances are developed. One numerical procedure applies Newton's method. Here, starting with soft walls, the eigenvalues are found by increasing the impedances of each wall pair in small increments up to the terminal impedances. Another procedure poses the eigenvalue problem as one of homotopic continuation from a non-physical reference configuration in which all eigenvalues are known and obvious. The continuation is performed by the numerical integration of two differential equations. The latter procedure was found to be faster and finds all possible solutions. The set of eigenvalues allowed the room modal natural frequencies and damping constants to be obtained. From sound decays measured in a hard-walled rectangular room, and from the collective-modal-decay curve, the impedances of the hard walls are estimated. These are then used to find the reverberation times of the modes in the room with the floor lined with sound absorbing material of known acoustic impedance. It was found that a single reverberation time, for all modes, is only supported in the rectangular room with hard walls and at the higher frequency bands, consistent with Sabine's theory, which assumes a diffuse sound field. In the rectangular room with hard walls and at the lower frequency bands, and in the rectangular room with the floor lined with sound absorbing material and for all frequency bands, modes with rather distinctive reverberation times may produce sound decays not always consistent with Sabine's prediction.     

Sound generated by a vortex convected past an elastic sheet

We study the motion and sound generated when a line vortex is convected in a uniform low-Mach flow parallel to a thin elastic sheet. The linearized sheet motion is analyzed under conditions where the unforced sheet (in the absence of the line vortex) is stationary. The vortex passage above the sheet excites a resonance mode of motion, where the sheet oscillates at its least stable eigenmode. The sources of sound in the acoustic problem include the sheet velocity and fluid vorticity. It is shown that the release of trailing-edge vortices, resulting from the satisfaction of the Kutta condition, has two opposite effects on sound radiation: while trailing-edge vortices act to reduce the pressure fluctuations occurring owing to the direct interaction of the line vortex with the unperturbed sheet, they extend and amplify the acoustic signal produced by the motion of the sheet. The sheet motion radiates higher sound levels as the system approaches its critical conditions for instability, where the effect of resonance becomes more pronounced. It is argued that the present theory describes the essential mechanism by which sound is generated as a turbulent eddy is convected in a mean flow past a thin elastic airfoil.     

三维边界层内定常横流涡的感受性研究

层流向湍流转捩的预测与控制一直是研究的前沿热点问题之一,其中感受性阶段是转捩过程中的初始阶段,它决定着湍流产生或形成的物理过程.但是有关三维边界层内感受性问题的数值和理论研究都比较少;实际工程问题中大部分转捩过程都是发生在三维边界层流中,所以研究三维边界层中的感受性问题显得尤为重要.本文以典型的后掠角45?无限长平板为例,数值研究了在三维壁面局部粗糙作用下的三维边界层感受性问题,探讨了三维边界层感受性问题与三维壁面局部粗糙长、宽和高之间的关系;然后,考虑在后掠平板上设计不同的三维壁面局部粗糙的分布状态、几何形状、距离后掠平板前缘的位置以及流向和展向设计多个三维壁面局部粗糙对三维边界层感受性问题有何影响;最后,讨论两两三维壁面局部粗糙中心点之间的距离以及后掠角的改变对三维边界层感受性的物理过程将会发生何种影响等.这一问题的深入研究将为三维边界层流中层流向湍流转捩过程的认识和理解提供理论依据.     

A numerical analysis on flow field and heat transfer by interaction between a pair of vortices in rectangular channel flow

This paper is a numerical study of the effect of flow field and heat transfer created by interactions between a pair of vortices generated by a vortex generator in a rectangular channel flow. In order to analyze the vortices produced by the vortex generator, the pseudo-compressibility method is introduced into the Navier–Strokes (NS) equation of a three-dimensional unsteady, incompressible viscous flow. A two-layer k–ε turbulence model is used on the flat plate three-dimensional turbulence boundary to predict the turbulence characteristics of the vortices. The computational results accurately predict the vortex characteristics, which are related to Reynolds stress, turbulent kinetic energy, and flow field. Also, in the prediction of thermal boundary layers, skin friction characteristics, and heat transfers, the present results are reasonably close to the experimental results obtained by other researchers.     

A computational and experimental study of resonators in three dimensions

In a previous work by the present authors, a computational and experimental investigation of the acoustic properties of two-dimensional slit resonators was carried out. The present paper reports the results of a study extending the previous work to three dimensions. This investigation has two basic objectives. The first is to validate the computed results from direct numerical simulations of the flow and acoustic fields of slit resonators in three dimensions by comparing with experimental measurements in a normal incidence impedance tube. The second objective is to study the flow physics of resonant liners responsible for sound wave dissipation. Extensive comparisons are provided between computed and measured acoustic liner properties with both discrete frequency and broadband sound sources. Good agreements are found over a wide range of frequencies and sound pressure levels. Direct numerical simulation confirms the previous finding in two dimensions that vortex shedding is the dominant dissipation mechanism at high sound pressure intensity. However, it is observed that the behavior of the shed vortices in three dimensions is quite different from those of two dimensions. In three dimensions, the shed vortices tend to evolve into ring (circular in plan form) vortices, even though the slit resonator opening from which the vortices are shed has an aspect ratio of 2.5. Under the excitation of discrete frequency sound, the shed vortices align themselves into two regularly spaced vortex trains moving away from the resonator opening in opposite directions. This is different from the chaotic shedding of vortices found in two-dimensional simulations. The effect of slit aspect ratio at a fixed porosity is briefly studied. For the range of liners considered in this investigation, it is found that the absorption coefficient of a liner increases when the open area of the single slit is subdivided into multiple, smaller slits.     

Non-nearest-neighbor interactions in nonlinear dynamical lattices

We revisit the theme of non-nearest-neighbor interactions in nonlinear dynamical lattices, in the prototypical setting of the discrete nonlinear Schrödinger equation. Our approach offers a systematic way of analyzing the existence and stability of solutions of the system near the so-called anti-continuum limit of zero coupling. This affords us a number of analytical insights such as the fact that, for instance, next-nearest-neighbor interactions allow for solutions with nontrivial phase structure in infinite one-dimensional lattices; in the case of purely nearest-neighbor interactions, such phase structure is disallowed. On the other hand, such non-nearest-neighbor interactions can critically affect the stability of unstable structures, such as topological charge S=2 discrete vortices. These analytical predictions are corroborated by numerical bifurcation and stability computations.     

E-CUSP scheme for the equations of ideal magnetohydrodynamics with high order WENO Scheme

An E-CUSP (energy-convective upwind and split pressure) scheme is developed to solve the equations of magnetohydrodynamics. Fifth order WENO reconstructions are employed to calculate the fluxes in order to achieve high order spacial accuracy. A characteristic speed of sound by averaging the fast wave speed and the acoustic speed of sound is suggested to evaluate the Mach number, which will yield robust and accurate solutions. The numerical experiments have demonstrated the accuracy and the capability of the new scheme to capture complex interactions of multiple shocks and vortices.     

Charged vortices in an abelian Higgs model with Chern-Simons term

The existence is shown of charged vortices of finite energy in a (2+1)-dimensional abelian Higgs model with Chern-Simons (C-S) term. Further, lower bounds on the energy and angular momentum of these vortices are obtained. Finally it is shown that the “glueballs” are absent in (2+1)-dimensional QED or QCD with C-S term for any compact gauge group.     

On vortices and rings in extended Abelian models

A numerical search for straight superconducting vortices in a U(1) model with a Ginzburg-Landau potential containing a cubic term is presented. Such vortices exist in a small numerically determined region. The reasons of their existence in that narrow region of the parameter space, as well as of their instability in the rest of the parameter space, are explained. Then, the results of a numerical search for axially symmetric solitons in a U(1)×U(1) model with higher derivative terms, which is based on [C.G. Doudoulakis, Physica D 228 (2007) 159], are presented and discussed.     

激波冲击V形界面重气体导致的壁面与旋涡作用及其对湍流混合的影响

基于多组分混合物质量分数模型,采用色散最小耗散可控的高分辨率有限体积方法,数值模拟了弱激波冲击V形空气/SF_6界面后,界面不稳定性生成的旋涡与固体壁面作用问题.激波冲击V形界面之后,因斜压效应诱导涡量沉积在界面附近,形成沿界面规则排列的多个涡对结构.旋涡的诱导作用使界面不断变形和卷起,同时旋涡之间不断发生相互并对,诱导更多更小尺度的旋涡产生.旋涡诱导作用的叠加效应,使界面尖端处的初始涡对向上下壁面发展.随后,涡结构开始与壁面发生复杂的相互作用.旋涡与壁面作用后沿壁面加速,使得物质界面沿壁面伸展,随后,旋涡从壁面回弹,并诱导二次旋涡产生.旋涡与壁面相互作用的过程,能够明显加剧物质混合.本文从物质混合的角度研究了该过程的机理,分析了旋涡与壁面作用对物质混合的影响.     

Roles of initial condition and vortex pairing in jet noise

This is a study of the effect of initial condition on sound generated by vortex pairing in a low Mach number, cold air jet (0·15 ⩽ M ⩽ 0·35). Data has been taken, both flow velocity fields and sound pressure far fields, in a quality anechoic facility, with careful documentation of the effect of initial condition on the sound field of jets of two different geometries (i.e., circular and elliptic). Explanations are presented for most of the observed effects by applying Möhring's theory of vortex sound to vortex filament models of coherent structures in the jets. The explanations also draw upon experience with coherent structure dynamics. The sound source of interest here is that associated with the pairing of shear layer vortices. The evolution of these vortices is greatly affected by the initial condition as is their resultant sound field. The elliptic jets with laminar boundary layers show azimuthal directivity, namely, sound pressure levels in the minor axis plane were greater than in the major axis plane. This difference decreases as the nozzle boundary layer undergoes natural transition with increasing jet speed. When the nozzle boundary layer is tripped, making it fully turbulent and removing the shear layer mode of pairing, the elliptic jet sound fields become nearly axisymmetric. What appears to be the most acoustically active phase of vortex pairing has been modeled, and the resulting sound field calculated for the circular jet. Supporting evidence is found in the experimental data for the validity of this model. The model explains the connection between the initial condition and the far field sound of jets. Interestingly, a general result of Möhring's theory is that motions of vortex rings (of any arbitrary shape) can produce only axisymmetric sound fields if the rings remain in a plane. This implies that the observed asymmetric directivity of the laminar elliptic jet sound field must be due to non-planar ring motions of the vortical structures. The primary contribution of this paper is to examine quantitatively the role of vortex pairing in the production of jet noise; the results are used to reemphasize that “pairing noise” cannot be dominant in most practical jet sound fields, contrary to claims by other researchers.     

Direct numerical simulation of a fully developed compressible wall turbulence over a wavy wall

Compressible turbulent channel flow over a wavy surface is investigated by direct numerical simulations using high-resolution finite difference schemes. The Reynolds number considered in the present paper is 3380 based on the bulk velocity, the channel half-width and the kinetic viscosity at the wall. Four test cases are simulated and analysed at Ma_m = 0.33, 0.8, 1.2, 1.5 based on the bulk velocity and the speed of sound at the wall. We mainly focus on the curvature and the Mach number effects on the compressible turbulent flows. Numerical results show that although the wavy wall has effects on the mean and fluctuation quantities, log law still exists in the distribution of the wave-averaged streamwise velocity if the roughness effects are taken into consideration in the scaling of it. Near-wall streaks are broken by the wavy surface and near-wall quasi-streamwise vortices mostly begin at the upslope of the wave and pass over the crest of it. The wavy wall makes the turbulence more active and the flow easier to be blended. From the viewpoint of turbulent kinetic budgets, curvature effects strengthen both the diffusion terms and the dissipation terms. At the same time, they change the properties of the compressibility-related terms and promote more inner energy transferring into turbulent kinetic energy. As the Mach number increases, the reattachment of the mean flow is delayed, which indicates the mean separation bubble becomes larger. Concerning the near-wall coherent structures, the vortices are more sparsely distributed with the increasing of the Mach number. For the supersonic cases, shock waves appear. Though they have little effects on the mean turbulent quantities, they change the structures of the flow fields and induce local separations at the upper wall of the channel.