用有限电流元法识别HL-2A 装置的等离子体边界

在电流丝模型对HL-2A装置等离子体边界重建取得满意结果[3]、[5]的基础上，本文用有限电流元模型对HL-2A等离子体边界进行了重建研究。计算和实验结果表明，在通常情况下，有限电流元模型比固定电流丝模型的重建误差稍小，前者误差小于3mm，后者误差小于6mm。当部分有限电流元分布在等离子体边界之外时，用有限电流元模型仍然可以成功重建边界。有限电流元的位置在一定范围内变化时，重建的误差都很小。用普通奔腾4 2.4GHz PC机计算一组等离子体放电边界的时间不超过1ms。有限电流元法能准确而快速地识别等离子体的偏滤器位形，这对于HL-2A装置的实时位形控制是基本和重要的。     

Development of a simplified field method of measuring sound insulation

A simplified field method of measuring the sound insulation of partitions within the framework of the British building regulations has been proposed and its development is described in this paper. The selection of the source type and the weighting function for the measurement of the overall sound pressure levels and the effects of loudspeaker and room response are discussed.A simple method of normalising for receiving room absorption is proposed and evaluated. Some results of experiments, carried out both in the laboratory and in field conditions, are reported. The results of this study show that the simplified field method has considerable potential for bringing about wide-scale enforcement of the building codes, which have so far been ineffective.     

An analysis of a two cylinder-fluid interaction at critical gap spacing by a cell boundary element method

Incompressible viscous, uniform flow past two parallel cylinders of equal diameter at a critical gap spacing g^*=0.90 to 1.10 is investigated using the viscous cell boundary element method. In this critical flow zone an irregular interchange of in-phase and anti-phase wake flows is observed and the switching mechanism between the states is investigated. Numerical studies undertaken at space increments of 0.05, including the central case g^*=1.0, detail the flow characteristics in time and space domains. Probability and power spectral analyses are presented to illustrate the statistical characteristics of the drag and lift flow parameters and vortex shedding patterns assessed through examination of their frequencies.     

运动声源的边界元声全息识别方法研究

提出了一种可实现任意形状的运动结构噪声源识别的声全息方法。通过结合移动框架技术与边界元声全息技术两种算法的特点,提出利用移动框架技术将存在多普勒效应的时域数据转换成边界元声全息所需的双平面全息数据,然后由边界元法声全息公式重构任意结构表面的声学信息,实现运动结构噪声源定位。该方法既具有移动框架技术处理运动问题的快速简便,又具有边界元方法可处理任意形状问题的特点。最后在半消声水池中,对运动速度为9.96 cm/s的带帽圆柱壳体进行了试验验证,结果表明:在低速条件下,该方法能够准确反演得到该结构的表面有功声强以及声压等声场信息,从而实现噪声源定位。由于条件有限,高速验证需进一步验证。     

The use of the boundary element method in solving the problems of sound scattering by an elastic noncircular cylinder

A semianalytic solution of the problem of sound scattering by an elastic cylinder with a noncircular cross section is proposed. The problem is solved by the analytic-numerical boundary element method in terms of the ideal fluid theory and the linear theory of elasticity.     

Assessing the effect of a barrier between two rooms subjected to low frequency sound using the boundary element method

The boundary element method (BEM) has been used to compute the acoustic wave propagation through a single vertical panel, which separates two rooms, made of concrete, when one of the rooms is excited by a steady-state, spatially sinusoidal, harmonic line load pressure at low frequencies. This work focuses on how the connection of the panel to the ceiling affects the acoustic insulation provided by the wall. Perfect double-fixed partitions and acoustic barrier-type structures with differently-sized gaps between the ceiling and the barrier are studied. The BEM model is formulated in the frequency domain and takes the air-solid interaction fully into account. Insulation dips are localised in the frequency domain and identified with dips associated with both the wall's natural dynamic vibration modes and with those associated with the air in the rooms. The influence of the wall's thickness on acoustic insulation is also analysed. The computed results obtained with the acoustic barrier type structure are compared with those obtained by a rigid model. The importance of the rooms' surface conditions is assessed, modelling the rooms with cork.     

小尺度阻尼边界封闭空间声传递函数的有限元素法计算

如何求解阻尼边界封闭空间中声源点到接收点的低频声传递函数已成为目前小尺度封闭空间可听化技术研究的关键技术，能处理任意形状及复杂边界条件的有限元素法可作为求解该问题的适合方法，以室内声声有源Helmholtz方程及其相应边界方程为基础，本文推导出了用于小尺度阻尼边界封闭空间声传递函数的有限元素求解方法，并编制了相应的计算机程序，在算例中，首先通过与模态叠加法计算结果进行比较，验证了该方法的正确性。最后计算了某型车体内腔中任意两点间声传递函数。     

Light intensity fluctuations on a semiconductor microsphere calculated by boundary element method

For a semiconductor microsphere irradiated by monochromatic unpolarized plane light wave, the governing Maxwell equations, which are transformed into Schrödinger equations through Debye potentials, are solved by the boundary element method, one of the integral formulations. The resultant light intensities on the particle surface show noise-like fluctuations depending on various parameters such as the light wavelength, the particle size, the numerical surface element size, etc. The more the numerical surface elements used, the greater the noise extent of the intensities. We think that this noise is related to the fluctuations happening in the real world and that they are somehow made into shape numerically by Green’s function and surface integration. One can consider a numerical surface element as a crystalline or amorphous unit cell. Actually a few experiments with photon energy conversion devices give the consistent results with ours that the energy conversion efficiency is on the increase with the decreasing size of unit cells. It is thus proposed here that this noise from numerical computation may be modeled to be the real thermal fluctuations of photon density on particle surface for photovoltaic cells, photocatalysts, photoluminescence devices, etc.     

求解随机振动结构声辐射的统计边界元方法

机器的机械振动在音频范围内大多是具有明显随机性的随机过程,普通边界元方法不能计算这类具有随机边界条件的结构声辐射问题。本文提出将边界元方法和统计方法相结合的统计边界元方法,导出了求解随机振动物体所辐射噪声声场的计算公式。通过求解随机声场算例的计算,表明本文所提出的方法是正确的。     

一种简化维里型状态方程预测高温甲烷PVT关系

为满足描述爆轰环境下高温气体高温、中高压状态的需求,本文提出了一种基于Lennard-Jones(LJ)势能函数的简化维里型状态方程Han-Long(HL).应用HL状态方程计算了甲烷1000 K以上112组理论和实验数据,计算所得体积平均绝对偏差约为1%,最大误差为3.28%,远低于DMW状态方程和BS状态方程的计算偏差.采用HL状态方程计算了甲烷冲击试验的热力学数据,计算所得体积偏差均小于3%.结果表明,HL状态方程能够很好的描述高温甲烷的热力学状态.     

Analysis of sound absorption by periodic structures using a hybrid-boundary element/mode expansion method

An improved hybrid method is proposed for analyzing sound scattering by a periodic structure. Part of the scattered field formulated with the mode expansion method is combined with other components of the field formulated with the boundary integration method in one period of the structure. Structures treated by this method can have arbitrary periodic forms made of locally reactive boundaries and porous materials. The oblique incident absorption coefficient of the structure is obtained simply from the reflection factor calculated for each elemental wave of the scattered field. The accuracy of the method is demonstrated by the agreement between calculated and measured values of the normal absorption coefficient of some test structures.     

Refraction and reflection of sound at the boundary of a bubbly liquid

Reflection and refraction of acoustic waves at the interface between pure water and bubbly water are investigated for the case of oblique incidence. From an analysis of analytic solutions, it is concluded that, for a wave incident on the interface from the side of a bubbly liquid, a critical angle of incidence, which depends on the frequency and the parameters of the disperse system, exists, so that, at angles of incidence exceeding the critical one, the wave is totally reflected from the interface.

     

用边界无法与四负载法联合预测排气消声器的插入损失

本文提出了将边界元法与四负载法相结合用于预测排气消声器插入损失的新方法,建立了具有低马赫数流时消声器四极参数求值的边界元计算方法,并将四负载法扩展到排气消声系统声源阻抗的求值。文中对两个消声器的插入损失进行了数值计算和实验测量,其预测结果与实测结果吻合良好。     

Noise simulation of aircraft engine fans by the boundary element method

Numerical simulation results of the civil aircraft engine fan stage noise in the far field are presented. Non-steady-state rotor–stator interaction is calculated the commercial software that solves the Navier–Stokes equations using differentturbulence models. Noise propagation to the far acoustic field is calculated by the boundary element method using acoustic Lighthill analogies without taking into account the mean current in the air inlet duct. The calculated sound pressure levels at points 50 m from the engine are presented, and the directional patterns of the acoustic radiation are shown. The use of the eddy resolving turbulence model to calculate rotor–stator interaction increases the accuracy in predicting fan stage noise.     

Nonlinear nonuniform torsional vibrations of bars by the boundary element method

In this paper a boundary element method is developed for the nonuniform torsional vibration problem of bars of arbitrary doubly symmetric constant cross-section taking into account the effect of geometrical nonlinearity. The bar is subjected to arbitrarily distributed or concentrated conservative dynamic twisting and warping moments along its length, while its edges are supported by the most general torsional boundary conditions. The transverse displacement components are expressed so as to be valid for large twisting rotations (finite displacement-small strain theory), thus the arising governing differential equations and boundary conditions are in general nonlinear. The resulting coupling effect between twisting and axial displacement components is considered and torsional vibration analysis is performed in both the torsional pre- or post-buckled state. A distributed mass model system is employed, taking into account the warping, rotatory and axial inertia, leading to the formulation of a coupled nonlinear initial boundary value problem with respect to the variable along the bar angle of twist and to an “average” axial displacement of the cross-section of the bar. The numerical solution of the aforementioned initial boundary value problem is performed using the analog equation method, a BEM based method, leading to a system of nonlinear differential-algebraic equations (DAE), which is solved using an efficient time discretization scheme. Additionally, for the free vibrations case, a nonlinear generalized eigenvalue problem is formulated with respect to the fundamental mode shape at the points of reversal of motion after ignoring the axial inertia to verify the accuracy of the proposed method. The problem is solved using the direct iteration technique (DIT), with a geometrically linear fundamental mode shape as a starting vector. The validity of negligible axial inertia assumption is examined for the problem at hand.     

Relationships between arithmetic averages of sound pressure level calculated in octave bands and Zwicker’s loudness level

Previously it has been found through a series of psychoacoustical experiments that the arithmetic average of sound pressure level calculated in octave bands is a good estimator of loudness for various kinds of environmental noise. Remarkably, the arithmetic average of sound pressure level in octave bands from 63 Hz to 4 kHz, L_m,1/1(63-4k), strongly correlates with the loudness level specified in ISO 532B, LL(Z), as well as with loudness assessment. To investigate this relationship further, a numerical study has been carried out based on Zwicker’s loudness model. As a result, practical expressions to estimate the loudness levels of general environmental noises from the sound pressure levels in octave bands from 63 Hz or 125 Hz to 4 kHz are proposed.     

Measurement of sound insulation with a sound level meter

The present technique for measuring the airborne sound insulation of walls and floors, involving measurements in 16 one-third octave bands, is tedious and expensive. The method provides more information than is needed for most purposes, and is more suited to research.Several investigators have proposed the measurement of the overall A-weighted sound level difference using a sound level meter, with a broad band source of white or pink noise. Consistent results have been obtained but their relation to accepted rating methods such as STC is rather empirical.The reference curves used for airborne sound insulation, i.e. STC and HPGW are very similar to the A weighting curve, and if the latter were adopted as the reference curve, there would be a firm theoretical basis for measurement with a sound level meter. Measurement of the difference between the linear sound level of a source of white noise, and the A weighted received level, would in practice be a test of the conformity to the A weighting curve of the transmission loss curve of the partition. Adverse deviations would show as a higher received level. Favourable deviations would have little effect.A study of practical walls and floors, taken from National Building Studies Research Paper 33, showed that there was good correlation between the sound level difference calculated as if it had been measured directly with a sound level meter, and a proposed rating method similar to ISO R717, but using the A weighting curve as the reference curve. Ninety-six per cent of results were within ±1 dB.The practical difficulties of achieving a reasonably flat transmitted spectrum, and of correcting for room absorption will reduce this precision, but bearing in mind the practical success of other short tests, the proposed test should provide a rapid test which is adequate for approval purposes.     

有限元+边界元——修正的模态分解法预报水下旋转薄壳的振动和声辐射

本文提出了一种有限元+边界元——修正的模态分解法用于预报水下结构的振动和声辐射问题.通过引入结构的剩余模态计入被忽略的高阶模态的准静态响应,使模态分解法的收敛性大大改善.文中以轴对称力作用下水下弹性薄球壳振动和声辐射问题为例进行了数值计算,给出了ka=1,2,3和4时球壳表面位移、声压和远场辐射声的计算结果.与原有的方法比较可见,该方法提高了模态分解法的收敛性,特别对远场辐射声指向性尤为明显。