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.     

A simulation of condenser microphones in free field by boundary element approach

A boundary element approach to simulate the characteristics of a microphone in free field is presented. To describe the approach, a condenser microphone is considered as a typical model. The usual boundary element technique is applied to the free acoustic field outside the microphone and a new special boundary element is developed to include the effect of a very thin air layer between the diaphragm and the backplate of the microphone. The vibration of the diaphragm is expressed in terms of its normal modes which are combined with the boundary elements to find the frequency response and the directivity. A numerical example is given for the condenser microphone used as the standard microphone in acoustical measurement, in which the diffraction from the preamplifier house and the nose cone is included. The computed results are in good comparison with the measured data in the operating frequency range.     

Focusing performance of the closed-boundary cylindrical microlenses analyzed by the boundary element method

In this paper, we investigate the focusing performance of closed-boundary cylindrical microlenses (CBCMs) based on rigorous electromagnetic theory and the boundary element method. The CBCMs with different incident angles, different quantization-level numbers, different microlens diameters, different f-numbers, and different polarizations of incidence are studied. Several focusing performance measures, such as the focal spot size, the diffraction efficiency, the real focal position, and the normalized transmitted power, are presented. It provides very useful information in designing the CBCMs in micro-optical systems.     

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

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

Terfenol-D超磁致伸缩换能器的有限元模拟

本文研究了利用ＡＮＳＹＳ软件分析磁致伸缩换能器,使这一有限元分析软件能够用于Ｔｅｒｆｅｎｏｌ－Ｄ换能器的分析设计。设计了１４ｋＨｚ Ｔｅｒｆｅｎｏｌ－Ｄ纵向换能器,计算结果与实验结果符合很好。     

A simulation method of aircraft plumes for real-time imaging

Real-time infrared simulation technology can provide a large number of infrared images under different conditions to support the design, test and evaluation of a system having infrared imaging equipment with very low costs. By synthesizing heat transfer, infrared physics, fluid mechanics and computer graphics, a real-time infrared simulation method is proposed based on the method of characteristics to predict the infrared feature of aircraft plumes, which tries to obtain a good balance between simulation precision and computation efficiency. The temperature and pressure distribution in the under-expansion status can be rapidly solved with dynamically changing flight statuses and engine working states. And a modified C–G (Curtis–Godson) spectral band model that combines the plume streamlines with the conventional C–G spectral band model was implemented to calculate the non-uniformly distributed radiation parameters inside a plume field. The simulation result was analyzed and compared with the CFD++, which validates the credibility and efficiency of the proposed simulation method.     

Calculation of the field enhancement on laser-illuminated scanning probe tips by the boundary element method

Received: 9 April 1997/Accepted: 7 October 1997     

Accelerated boundary element method for diffuse optical imaging

The boundary element method (BEM) is a useful tool in diffuse optical imaging (DOI) when modelling large optical regions whose parameters are piecewise constant, but are computationally expensive. We present here an acceleration technique, the single-level fast multipole method, for a highly lossy medium. The enhanced practicability of the BEM in DOI is demonstrated through test examples on single-layer problems, where order of magnitude reduction factors on solution time are achieved and on a realistic three-layer model of the neonatal head. Our experimental results agree very closely with theoretical predictions of computational complexity.     

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.     

Numerical simulation of laser-generated ultrasound in non-metallic material by the finite element method

The use of a pulsed laser for the generation of the elastic waves in non-metallic materials in the thermoelastic regime is investigated by using finite element method (FEM), taking into account not only thermal diffusion and the finite spatial and temporal shape of the laser pulse, but also optical penetration and the temperature dependence of material properties. The optimum finite element model is established based on analysis of two important parameters, meshing size and time step, and the stability of solution. Temperature distributions and temperature gradient fields in non-metallic material for different time steps are obtained, this temperature field is equivalent to a bulk force source to generate ultrasonic wave. The laser-generated ultrasound waveforms at the epicenter and surface acoustic waveforms (SAWs) are obtained and the influence of optical penetration into the material on the temperature field and the ultrasound waveforms are analyzed. The numerical results indicate that the heat penetration into non-metallic material is caused mainly by the optical penetration, and the ultrasound waveforms, especially the shape of the precursor, are strongly dependent on the optical penetration depth into non-metallic material.     

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

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

Three-dimensional simulation of elastic capsules in shear flow by the penalty immersed boundary method

An improved penalty immersed boundary method (pIBM) has been proposed for simulation of flow-induced deformation of three-dimensional (3D) elastic capsules. The motion of the capsule membrane is described in the Lagrangian coordinates. The membrane deformation takes account of the bending and twisting effects as well as the stretching and shearing effects. The method of subdivision surfaces is adopted to generate the mesh of membrane and the corresponding shape functions, which are required to be C¹ continuous. The membrane motion is then solved by the subdivision-surface based finite element method on the triangular unstructured mesh. On the other hand, the fluid motion is defined on the Eulerian domain, and is advanced by the fractional step method on a staggered Cartesian grid. Coupling of the fluid motion and the membrane motion is realized in the framework of the pIBM. Using the proposed method, deformation of 3D elastic capsules in a linear shear flow is studied in detail, and validations are examined by comparing with previous studies. Both the neo-Hookean membrane and the Skalak membrane are tested. For an initially spherical capsule the tank-treading motion is formed under various dimensionless shear rates and reduced bending moduli. It is found that buckling occurs near the equator of the capsule for small shear rates but near the tips for large shear rates, which is suppressed by including the bending rigidity of the membrane. Effects of the Reynolds number and the membrane density are investigated for an initially spherical capsule. For a non-spherical capsule, with the initial shape of the oblate spheroid or the biconcave circular disk as a model of the red blood cell, the swinging motion is observed due to the shape memory effect. By decreasing the dimensionless shear rate or increasing the reduced bending modulus, the swinging motion is transited into the tumbling motion.     

Noise control mechanisms of inside aircraft

World trends in the development of methods and approaches to noise reduction in aircraft cabins are reviewed. The paper discusses the mechanisms of passive and active noise and vibration control, application of “smart” and innovative materials, new approaches to creating all fuselage-design elements, and other promising directions of noise control inside aircraft.     

On the modeling of narrow gaps using the standard boundary element method

Numerical methods based on the Helmholtz integral equation are well suited for solving acoustic scattering and diffraction problems at relatively low frequencies. However, it is well known that the standard method becomes degenerate if the objects that disturb the sound field are very thin. This paper makes use of a standard axisymmetric Helmholtz integral equation formulation and its boundary element method (BEM) implementation to study the behavior of the method on two test cases: a thin rigid disk of variable thickness and two rigid cylinders separated by a gap of variable width. Both problems give rise to the same kind of degeneracy in the method, and modified formulations have been proposed to overcome this difficulty. However, such techniques are better suited for the so-called thin-body problem than for the reciprocal narrow-gap problem, and only the first is usually dealt with in the literature. A simple integration technique that can extend the range of thicknesses/widths tractable by the otherwise unmodified standard formulation is presented and tested. This technique is valid for both cases. The modeling of acoustic transducers like sound intensity probes and condenser microphones has motivated this work, although the proposed technique has a wider range of applications.     

Fan noise reduction of an aircraft engine by inclining the outlet guide vanes

Acoustical Physics - The paper presents the results of an experimental study of the acoustic characteristics of an aircraft engine fan in the design of which the outlet guide vanes are leaned in...     

Effectiveness of combined aircraft engine noise suppressors

We consider the design features of fan noise suppressors in application to air intakes and the bypass duct of a turbofan engine. A combined liner is developed that has increased acoustic efficiency in comparison to conventional honeycomb liner. We demonstrate the important role of the area of the sound-absorbing liner between fan Rotor and Stator ensuring significant noise reduction.     

快速多极边界元法用于扩散光学断层成像研究

在求解扩散光学断层成像中的正向问题时, 目前普遍采用有限元法, 但是随着实际模型规模的增大, 有限元法的计算量问题日益显著, 而边界元法则由于可以降低计算维度使计算量减少而备受关注. 本文以均匀的高散射介质为模型, 研究了将快速多极边界元法用于扩散光学断层成像的正向问题. 快速多极边界元法利用核函数的多极展开, 将常规边界元法中系数矩阵和迭代矢量的乘积项等价为相应四叉树结构的一次递归, 再结合广义最小残量法进行迭代求解. 将计算结果和蒙特卡罗法的模拟结果进行了比较, 表明利用快速多极边界元法的模拟结果和蒙特卡罗法的结果有很好的一致性. 研究结果验证了快速多极边界元法可以用于扩散光学断层成像, 为其大规模和实时成像带来可观的前景. 关键词： 扩散光学断层成像 边界元法 快速多极边界元法     

Stabilization of time domain acoustic boundary element method for the interior problem with impedance boundary conditions

The time domain boundary element method (BEM) is associated with numerical instability that typically stems from the time marching scheme. In this work, a formulation of time domain BEM is derived to deal with all types of boundary conditions adopting a multi-input, multi-output, infinite impulse response structure. The fitted frequency domain impedance data are converted into a time domain expression as a form of an infinite impulse response filter, which can also invoke a modeling error. In the calculation, the response at each time step is projected onto the wave vector space of natural radiation modes, which can be obtained from the eigensolutions of the single iterative matrix. To stabilize the computation, unstable oscillatory modes are nullified, and the same decay rate is used for two nonoscillatory modes. As a test example, a transient sound field within a partially lined, parallelepiped box is used, within which a point source is excited by an octave band impulse. In comparison with the results of the inverse Fourier transform of a frequency domain BEM, the average of relative difference norm in the stabilized time response is found to be 4.4%.     

A direct mixed-body boundary element method for packed silencers

Bulk-reacting sound absorbing materials are often used in packed silencers to reduce broadband noise. A bulk-reacting material is characterized by a complex mean density and a complex speed of sound. These two material properties can be measured by the two-cavity method or calculated by empirical formulas. Modeling the entire silencer domain with a bulk-reacting lining will involve two different acoustic media, air and the bulk-reacting material. Traditionally, the interior silencer domain is divided into different zones and a multi-domain boundary element method (BEM) may be applied to solve the problem. However, defining different zones and matching the elements along each interface is tedious, especially when the zones are intricately connected. In this paper, a direct mixed-body boundary element method is used to model a packed silencer without subdividing it into different zones. This is achieved by summing up all the integral equations in different zones and then adding the hypersingular integral equations at interfaces. Several test cases, including a packed expansion chamber with and without an absorbing center bullet, and a parallel baffle silencer, are studied. Numerical results for the prediction of transmission loss (TL) are compared to experimental data.