A mode matching method for modeling dissipative silencers lined with poroelastic materials and containing mean flow

A mode matching method for predicting the transmission loss of a cylindrical shaped dissipative silencer partially filled with a poroelastic foam is developed. The model takes into account the solid phase elasticity of the sound-absorbing material, the mounting conditions of the foam, and the presence of a uniform mean flow in the central airway. The novelty of the proposed approach lies in the fact that guided modes of the silencer have a composite nature containing both compressional and shear waves as opposed to classical mode matching methods in which only acoustic pressure waves are present. Results presented demonstrate good agreement with finite element calculations provided a sufficient number of modes are retained. In practice, it is found that the time for computing the transmission loss over a large frequency range takes a few minutes on a personal computer. This makes the present method a reliable tool for tackling dissipative silencers lined with poroelastic materials.     

Transmission loss predictions for dissipative silencers of arbitrary cross section in the presence of mean flow

A numerical technique is developed for the analysis of dissipative silencers of arbitrary, but axially uniform, cross section. Mean gas flow is included in a central airway that is separated from a bulk reacting porous material by a concentric perforate screen. The analysis begins by employing the finite element method to extract the eigenvalues and associated eigenvectors for a silencer of infinite length. Point collocation is then used to match the expanded acoustic pressure and velocity fields in the silencer chamber to those in the inlet and outlet pipes. Transmission loss predictions are compared with experimental measurements taken for two automotive dissipative silencers with elliptical cross sections. Good agreement between prediction and experiment is observed both without mean flow and for a mean flow Mach number of 0.15. It is demonstrated also that the technique presented offers a considerable reduction in the computational expenditure when compared to a three-dimensional finite element analysis.     

均匀流直通穿孔消声器的声学特性分析

将数值模态匹配法(NMM)扩展应用于计算有均匀流存在时直通穿孔管抗性和阻性消声器的声学特性,编写了相应的计算程序。对于圆形同轴穿孔管抗性和阻性消声器,应用数值模态匹配法计算得到的传递损失结果与实验测量结果吻合良好,从而验证了计算方法和计算程序的正确性。进而应用数值模态匹配法研究了运流效应和穿孔阻抗以及穿孔管偏移对穿孔管抗性和阻性消声器传递损失的影响。研究结果表明,马赫数越高,穿孔管抗性消声器在中高频的消声量越高,阻性消声器在整体频段内的消声性能越差;低马赫数时运流效应对穿孔管抗性消声器的影响可以忽略,马赫数较高时运流效应和穿孔阻抗的影响比较明显;对于穿孔管阻性消声器,穿孔阻抗对消声器声学特性的影响比运流效应的影响小,但是与真实值的差别不可忽略;穿孔管偏移对消声器声学特性的影响与频率和消声器结构均相关。     

均匀流直通穿孔消声器的声学特性分析

将数值模态匹配法(NMM)扩展应用于计算有均匀流存在时直通穿孔管抗性和阻性消声器的声学特性,编写了相应的计算程序。对于圆形同轴穿孔管抗性和阻性消声器,应用数值模态匹配法计算得到的传递损失结果与实验测量结果吻合良好,从而验证了计算方法和计算程序的正确性。进而应用数值模态匹配法研究了运流效应和穿孔阻抗以及穿孔管偏移对穿孔管抗性和阻性消声器传递损失的影响。研究结果表明,马赫数越高,穿孔管抗性消声器在中高频的消声量越高,阻性消声器在整体频段内的消声性能越差;低马赫数时运流效应对穿孔管抗性消声器的影响可以忽略,马赫数较高时运流效应和穿孔阻抗的影响比较明显;对于穿孔管阻性消声器,穿孔阻抗对消声器声学特性的影响比运流效应的影响小,但是与真实值的差别不可忽略;穿孔管偏移对消声器声学特性的影响与频率和消声器结构均相关。     

A three-dimensional finite element approach for predicting the transmission loss in mufflers and silencers with no mean flow

A three-dimensional finite element method has been implemented to predict the transmission loss of a packed muffler and a parallel baffle silencer for a given frequency range. Iso-parametric quadratic tetrahedral elements have been chosen due to their flexibility and accuracy in modeling geometries with curved surfaces. For accurate physical representation, perforated plates are modeled with complex acoustic impedance while absorption linings are modeled as a bulk media with a complex speed of sound and mean density. Domain decomposition and parallel processing techniques are applied to address the high computational and memory requirements. The comparison of the computationally predicted and the experimentally measured transmission loss shows a good agreement.     

Analytic mode matching for a circular dissipative silencer containing mean flow and a perforated pipe

An analytic mode matching scheme that includes higher order modes is developed for a straight-through circular dissipative silencer. Uniform mean flow is added to the central airway and a concentric perforated screen separates the mean flow from a bulk reacting porous material. Transmission loss predictions are compared with experimental measurements and good agreement is demonstrated for three different silencers. Furthermore, it is demonstrated that, when mean flow is present, the axial kinematic matching condition should equate to that chosen for the radial kinematic boundary condition over the interface between the airway and the material. Accordingly, if the radial matching conditions are continuity of pressure and displacement, then the axial matching conditions should also be continuity of pressure and displacement, rather than pressure and velocity as previously thought. When a perforated screen is present the radial pressure condition changes, but the radial kinematic condition should always remain equivalent to that chosen for the axial kinematic matching condition; here, results indicate that continuity of displacement should be retained when a perforated screen is present.     

A finite element method for modelling dissipative mufflers with a locally reactive lining

The transmission loss of a lined expansion chamber muffler with zero mean flow is studied by using a finite element model of the system. In the cases considered the lining is arranged in small cells so that it is reasonable to assume a locally reacting effect. An existing empirical relationship is used to calculate the lining impedance. The results show that the lining always increases the transmission loss over the unlined chamber valuys; with a thin lining its effect is similar to that of a side branch filter. However, as the thickness increases and the lining impedance approaches the characteristic impedance for the absorption material, the transmission loss develops a single maximum similar to the attenuation produced by a thick lining on the primary mode in a lined duct.     

A comparison between Lanczos and moment methods

A simple procedure is suggested for generating the 2n - 1 Lanczos tridiagonal elements from the 2n - 1 moments of any hamiltonian with regard to a given starting vector. The numerical feasibility is studied in the quartic oscillator example.     

A numerical analytic method for electromagnetic radiation accompanying with fracture of rocks

This paper studies Rabinovitch’s compression experiments on granite and chalk and proposes an oscillating dipole model to analyse and simulate the electromagnetic radiation phenomenon caused by fracture of rocks. Our model assumes that the electromagnetic radiation pulses are initiated by vibrations of the charged rock grains on the tips of the crack. The vibrations of the rock grains are stimulated by the pulses of the cracks. Our simulations show comparable results with Rabinovitch’s compression experiments. From the simulation results, it verifies an assumption that the crack width is inversely proportional to the circular frequency electromagnetic radiation, which is presented by Rabinovitch et al. The simulation results also imply that, by using our oscillating dipole model together with Rabinovitch’s two equations about the crack length and crack width, we can quantitatively analyse and simulate the electromagnetic radiation phenomenon, which is induced from the fracture of the rocks.     

A numerical comparison between two unconstrained variational formulations

In an effort to relieve the often cumbersome burden of meeting the requirements of the end conditions and to unify the solution formulation for boundary- and initial-value problems, unconstrained variational statements have been introduced in conjunction with some approximate methods. In the case of a boundary value problem, it is shown in this paper that two different variational statements can be established: one is arrived at by the use of the Lagrange multipliers, the other by energy considerations. The numerical convergence of the solutions associated with finite element schemes involving use of one of these two different variational statements is compared with that of the other. In the case of an initial value problem, both formulations can again be established when the adjoint field variable and the adjoint variational statement are introduced. The numerical data presented here indicate that while both methods generate excellent convergent results for the boundary problem, the method of stiff springs yields results which show much better convergence for the initial value problem than those achieved by Lagrange multipliers.     

High-order shock-fitting methods for direct numerical simulation of hypersonic flow with chemical and thermal nonequilibrium

In recent years, much progress has been made in the direct numerical simulation of laminar-turbulent transition of hypersonic boundary layer flow. However, most of the efforts at the direct numerical simulation of transition previously have been focused on the idealized perfect gas flow or “cold” hypersonic flows. For practical problems in hypersonic flows, high-temperature effects of thermal and chemical nonequilibrium are important and cannot be modeled by a perfect gas model. Therefore, it is necessary to include the real gas models in the numerical simulation of hypersonic boundary layer transition in order to accurately predict flow field parameters. Currently most numerical methods for hypersonic flow with thermo-chemical nonequilibrium are based on shock-capturing approach at relatively low order of accuracy. Shock capturing schemes reduce to first-order accuracy near the shock and have been shown to produce spurious oscillations behind curved strong shocks. There is a need to develop new methods capable of simulating nonequilibrium hypersonic flow fields with uniformly high-order accuracy and avoid spurious oscillations near the shock. This paper presents a fifth-order shock-fitting method for numerical simulation of thermal and chemical nonequilibrium in hypersonic flows. The method is developed based on the state-of-the-art real gas models for thermo-chemical nonequilibrium and transport phenomena. Shock-fitting approach is used because it has the advantage of capturing the entire flow field with high-order accuracy and without any oscillations near the shock. The new method has been tested and validated for a number of test cases over a wide span of free stream conditions. The developed method is applied for the study of receptivity of free stream acoustic waves over a blunt cone for hypervelocity flow. Some preliminary results of the computations of the high order shock fitting method for the above mentioned study have also been presented.     

A Trefftz-based numerical modelling framework for Helmholtz problems with complex multiple-scatterer configurations

The aim of the work presented in this paper is the numerical solution of low- and mid-frequency time-harmonic acoustic multiple-scattering problem. A novel so-called ‘multi-level’ modelling approach is proposed which is applicable to the study of a configuration of well separated obstacles of arbitrary shape on which any type of acoustic boundary condition can be applied. The generic character of the method is obtained by embedding the superposition principle for the multiple-scattering influence in a state-of-the-art acoustic modelling technique, the so-called Wave Based Method. The resulting approach successfully alleviates the geometrical limitations of the underlying Trefftz-based method and preserves the method’s computational efficiency, resulting in a generic multiple-scattering modelling framework with a superior computational efficiency in the low- as well as the mid-frequency range. Several numerical validation examples show that the proposed approach is as accurate as the classical single-scattering Wave Based Method and illustrate the computational efficiency as compared to Boundary Element Methods.     

Qualitative comparison between numerical and experimental results of unsteady flow in a radial diffuser pump

Comparison between numerical simulation and experimental results for unsteady flow field in a radial diffuser pump is presented for the design operating point. The numerical result is obtained by solving three-dimensional, unsteady Reynolds-averaged Navier-Stokes equations by the commercial CFD code CFX-10 withk-ω based shear stress transport turbulence model. Two-dimensional PIV measurements are conducted to acquire the experiment result. The phase-averaged velocity and turbulent kinetic energy fields are compared in detail between the results by the two methods in the impeller, diffuser and return channel regions. The qualitative comparison between CFD and PIV results is quite good in the phase-averaged velocity field. Although the turbulence level by PIV is higher than that by CFD generally, the main turbulence features are nearly the same. Furthermore, the blade orientation effect and other associated unsteady phenomena are also examined, in order to enhance the understanding on impeller-diffuser interaction in a radial diffuser pump.     

A hybrid numerical method for interfacial fluid flow with soluble surfactant

We address a significant difficulty in the numerical computation of fluid interfaces with soluble surfactant that occurs in the physically representative limit of large bulk Peclet number Pe. At the high values of Pe in typical fluid-surfactant systems, there is a transition layer near the interface in which the surfactant concentration varies rapidly, and large gradients at the interface must be resolved accurately to evaluate the exchange of surfactant between the interface and bulk flow. We use the slenderness of the layer to develop a fast and accurate ‘hybrid’ numerical method that incorporates a separate, singular perturbation analysis of the dynamics in the transition layer into a full numerical solution of the interfacial free boundary problem. The accuracy and efficiency of the method is assessed by comparison with a more ‘traditional’ numerical approach that uses finite differences on a curvilinear coordinate system exterior to the bubble, without the separate transition layer reduction. The traditional method implemented here features a novel fast calculation of fluid velocity off the interface.     

Efficient numerical methods for multiple surfactant-coated bubbles in a two-dimensional stokes flow

We present efficient and highly accurate numerical methods to compute the deformation of surfactant-coated, two-dimensional bubbles in a slow viscous flow. Surfactant acts to locally alter the surface tension and thereby change the nature of the interface motion. In this paper, we restrict our attention to the case of a dilute insoluble surfactant. The convection–diffusion equation for the surfactant concentration on the interface is coupled with the Stokes equations in the fluid domain through a boundary condition based on the Laplace-Young condition. The Stokes equations are first recast as an integral equation and then solved using a fast-multipole accelerated iterative procedure. The computational cost per time-step is only O(N log N) operations, with N being the number of discretization points on the interface. The bubble interfaces are described by a spectral mesh and is advected according to the fluid velocity in such a manner so as to preserve equal arc length spacing of marker points. This equal arc length framework has the dual advantage of dynamically maintaining the spatial mesh and allowing efficient, implicit treatment of the stiffest terms in the dynamics. Several phenomenologically different examples are presented.     

低温等离子体数值模拟方法的分析比较

分析比较了低温等离子体模拟中采用的流体模型、粒子模型和混合模拟方法及在放电特性分析中采用的电路模拟方法。给出了每种方法的基本原理,探讨了它们的适应性。利用粒子模型对外磁镜场作用下四阳极装置辉光放电所产生等离子体进行了模拟,分析了磁场对电子密度径向分布的影响。     

Improved mean field studies of SU(N) chiral models and comparison with numerical simulations

SU(N) × SU(N) matrix models are investigated using a systematic expansion around the classical mean field result. One-loop corrections for N = 2?5, and ∞ are computed for the two-, three-, and four-dimensional models. The results are in good agreement with numerical simulations of these models.     

Energy flow model for thin plate considering fluid loading with mean flow

Energy Flow Analysis (EFA) has been developed to predict the vibration energy density of system structures in the high frequency range. This paper develops the energy flow model for the thin plate in contact with mean flow. The pressure generated by mean flow affects energy governing equation and power reflection–transmission coefficients between plates. The fluid pressure is evaluated by using velocity potential and Bernoulli's equation, and energy governing equations are derived by considering the flexural wavenumbers of a plate, which are different along the direction of flexural wave and mean flow. The derived energy governing equation is composed of two kinds of group velocities. To verify the developed energy flow model, various numerical analyses are performed for a simple plate and a coupled plate for several excitation frequencies. The EFA results are compared with the analytical solutions, and correlations between the EFA results and the analytical solutions are verified.     

经验正交函数与广义数值环境模式重构声速剖面的比较

本文利用Argo资料对经验正交函数(EOF)和GDEM模式拟合声速剖面的效果进行了比较。结果表明,两种方法都能较好地重构局地声速剖面。EOF方法的拟合精度与选取的EOF阶数有关,阶数越高,精度越高;GDEM模式的拟合精度介于EOF的前5阶到前10阶之间。EOF易受到样本数量和采样深度的限制,难以提供完整的声速剖面;而GDEM模式能够将剖面扩展到更深的深度,与仅有上层海洋采样数据的声速剖面相比,结合深海水团的信息可使拟合结果达到更高的精度。