Vibro-acoustic analysis of complex systems

A general method is presented for predicting the ensemble average steady-state response of complex vibro-acoustic systems that contain subsystems with uncertain, or random, properties. The method combines deterministic and statistical techniques to produce a non-iterative hybrid method that incorporates equations of dynamic equilibrium and power balance. The method is derived explicitly without reference to statistical energy analysis (SEA); however, it is seen that the wave approach to SEA can be viewed as a special case of the proposed method. The proposed method provides a flexible way to account for necessary deterministic details in a vibro-acoustic analysis without requiring that an entire system be modeled deterministically. The method therefore provides a potential solution to the mid-frequency problem (in which a system is neither entirely deterministic nor entirely statistical). The application of the method is illustrated with a numerical validation example.     

Quantitative estimation of ultrasonic attenuation in a solid in the immersion case with correction of diffraction effects

This paper presents a method of diffraction correction for the log-spectral difference method to estimate quantitatively attenuation of a solid in the immersion case. The correction method is established based on the angular spectrum approach that is used to calculate the echoes from the front and back surfaces of the immersed solid. An example is given of a copper plate submerged in water and inspected by a linear array with a cylindrically curved surface. The correction method is first applied to a theoretical estimation of attenuation which is linearly dependent on frequency. The results have shown that the evaluated attenuation coefficient is in excellent agreement with the exact value. Then the method is applied to a real situation, in which the results have shown that the method yields reasonable evaluated attenuation values. This work has demonstrated that the method is able to correct effectively the diffraction effect so as to achieve a quantitative estimation of attenuation.     

Numerical method satisfying the first two conservation laws for the Korteweg–de Vries equation

In this paper, we develop a finite-volume scheme for the KdV equation which conserves both the momentum and energy. The main ingredient of the method is a numerical device we developed in recent years that enables us to construct numerical method for a PDE that also simulates its related equations. In the method, numerical approximations to both the momentum and energy are conservatively computed. The operator splitting approach is adopted in constructing the method in which the conservation and dispersion parts of the equation are alternatively solved; our numerical device is applied in solving the conservation part of the equation. The feasibility and stability of the method is discussed, which involves an important property of the method, the so-called Jensen condition. The truncation error of the method is analyzed, which shows that the method is second-order accurate. Finally, several numerical examples, including the Zabusky–Kruskal’s example, are presented to show the good stability property of the method for long-time numerical integration.     

FAAS法对三种制备样品方法的研究

在FAAS法分析中，样品预处理始终占有极为重要的位置。本文比较了三种样品消化方式，即传统酸消化法，灰化法及微波消化法。传统的消化方式过程长、速度慢、效率低、而且被测元素易受到损失及易污染等不足，而微波是一种非常快捷、省时、省试剂和无污染的消化方式。通过对校准物质贻贝中的铜、锌、铁的测定结果证明微波消化方式比传统的消化方式优点多。     

Refined multiload method for measuring acoustical source characteristics of an intake or exhaust system

The one-port source characteristics in a duct system, viz., source impedance and strength, can be determined by using the four-load method. In this paper, to avoid the instability problem of the conventional four-load method, a new formulation for the multiload method has been proposed, which employs an error function based on the linear, time-invariant source model. It is shown that the method is less sensitive to input errors compared to the previous methods. For a 10% input error, the proposed method yields a relative error in the source resistance that is about 1/100 times smaller than for the conventional method. The effectiveness of the present method is demonstrated by two test examples, a loudspeaker and a blower, each operating in a duct. It is observed that the conventional and least-squares methods result in large errors, whereas the present method yields far better agreement with the actual source parameters, as measured by the direct method. The present method is then used to obtain the source parameters on the exhaust side of an operating internal combustion engine. The radiated sound spectrum from the exhaust opening is predicted by using the measured source parameters and the calculated result agrees very well with the measured one.     

A finite difference method of solving anisotropic scattering problems

A new method of solving radiative transfer problems is described including a comparison of its speed with that of the doubling method, and a discussion of its accuracy and suitability for computations involving variable optical properties. The method uses a discretization in angle to produce a coupled set of first-order differential equations which are integrated between discrete depth points to produce a set of recursion relations for symmetric and anti-symmetric angular sums of the radiation field at alternate depth points. The formulation given here includes depth-dependent anisotropic scattering, absorption, and internal sources, and allows arbitrary combinations of specular and non-Lambertian diffuse reflection at either or both boundaries. The method is shown to be faster than the doubling method when the number of depth points and angular quadrature points is identical. Numerical tests of the method show that it can return accurate emergent intensities even for large optical depths. The method is also shown to conserve flux to machine accuracy in conservative atmospheres. Finally, several checks are made that demonstrate that the new method can compute accurate radiation fields in atmospheres with variable optical properties.     

Topology optimization using the improved element-free Galerkin method for elasticity

The improved element-free Galerkin(IEFG) method of elasticity is used to solve the topology optimization problems.In this method, the improved moving least-squares approximation is used to form the shape function. In a topology optimization process, the entire structure volume is considered as the constraint. From the solid isotropic microstructures with penalization, we select relative node density as a design variable. Then we choose the minimization of compliance to be an objective function, and compute its sensitivity with the adjoint method. The IEFG method in this paper can overcome the disadvantages of the singular matrices that sometimes appear in conventional element-free Galerkin(EFG) method. The central processing unit(CPU) time of each example is given to show that the IEFG method is more efficient than the EFG method under the same precision, and the advantage that the IEFG method does not form singular matrices is also shown.     

Eigenmode analysis of arbitrarily shaped two-dimensional cavities by the method of point-matching

In the proposed method, as in the boundary element method (BEM), a field problem is solved on its boundary along which nodes are distributed. The distinct feature of the proposed method is related to the fact that no interpolation functions between the nodes are required, so that the basic collocation method is employed to satisfy a given boundary condition. This approach reduces a large amount of numerical calculation induced due to the interpolation functions. Unlike the method of auxiliary sources (MAS), in the proposed method only the regular functions that have finite functional values everywhere including the origins are employed as basis functions. As a result, the proposed method does not need the auxiliary surface, the construction of which is a troublesome work because its form, dimension, shifting, etc. strongly influence the accuracy of the MAS. The eigenvalues calculated by the proposed method rapidly converge to the exact values thanks to the simplicity of the method. Moreover, the method gives mode shapes successfully without using the interpolation functions between the nodes. The efficiency and accuracy of the method are verified through several applications.     

一种晶体光纤基模色散特性的矢量法分析

利用矢量有效折射率方法(矢量法)对光子晶体光纤基模的色散特性进行了数值模拟，并与双正交归一基矢量法以及标量有效折射率方法(标量法)的模拟结果进行了对比.发现所用矢量法的结果与双正交归一基矢量法的结果符合很好，而标量法在低空气填充率f或较高归一化波数A/λ时是一种较好的近似，在空气填充率f较高或归一化波数A/λ较低时，要得到精确的结果必须利用矢量法对光子晶体光纤的特性进行模拟.讨论了光子晶体光纤包层有效折射率与光纤结构的关系. 关键词： 光子晶体光纤 矢量法 有效折射率 色散     

求解辐射传递方程的离散坐标法

阐述了含吸收散射性介质三维空腔内辐射传递方程的离散坐标解法。讨论了入射散射项积分格式的选取,以及假散射和射线效应对解精度的影响。对三维矩形炉膛内辐射传递过程进行了数值模拟,并与区域法和离散传递法进行比较。比较结果表明离散坐标法具有较好的精度,是目前燃烧室内辐射传热过程数值模拟的一种较好的方法。     

Accurate quantification of (1)H spectra: from finite impulse response filter design for solvent suppression to parameter estimation.

A scheme for accurate quantification of (1)H spectra is presented. The method uses maximum-phase finite impulse response (FIR) filters for solvent suppression and an iterative nonlinear least-squares (NLLS) algorithm for parameter estimation. The estimation algorithm takes the filter influence on the metabolites of interest into account and can thereby correctly incorporate a large variety of prior knowledge into the estimation phase. The FIR filter is designed in such a way that no distortion of the important initial samples is introduced. The FIR filter method is compared numerically with the HSVD method for water signal removal in a number of examples. The results show that the FIR method, using an automatic filter design scheme, slightly outperforms the HSVD method in most cases. The good performance and ease of use of the FIR filter method combined with its low computational complexity motivate the use of the proposed method.     

基于孟塞尔系统的阴极射线管特性化新方法

提出了一种新的阴极射线管特性化的方法。该方法的特点是采用“视觉匹配”方法,在反射体表面色和白发光体之间映射一些色貌因素,但没有使用任何复杂的色貌模型。是一种考虑了一些色貌因素的阴极射线管特性化方法。由于该问题的个性因素较多,采用BP神经网络实现色空间的非线性映射。实验结果表明,只要阴极射线管被标定,在办公室环境下,该方法可以改进在不同的阴极射线管上重现的颜色。采用3-7-7-3简单的网络结构;分色相样本训练。训练样本平均色差可以达到3．07L^*u^*v^*色差单位,测试样本平均色差可以达到4．55L^*u^*v^*色差单位,小于阴极射线管的最大可接受色差,结果是令人满意的。这在电子商务和办公自动化方面有广泛的用途。     

Synthetic aperture imaging using sources with finite aperture: deconvolution of the spatial impulse response

A method for ultrasonic synthetic aperture imaging using finite-sized transducers is introduced that is based on a compact, linear, discrete model of the ultrasonic measurement system developed using matrix formalism. Using this model a time-domain algorithm for deconvolution of the transducer's spatial impulse responses (SIRs) is developed that is based on a minimum mean square error (MMSE) criterion. The algorithm takes the form of a spatiotemporal filter that compensates for the SIRs associated with a finite-sized transducer at every point of the processed image. A major advantage of the proposed method is that it can be used for any transducer, provided that its associated SIRs are known. This is in contrast to the synthetic aperture focusing technique (SAFT), which treats the transducer as a point source. The performance of the method is evaluated with simulations and experiments, performed in water using a linear phased array. The results obtained using the proposed method are compared to those obtained with a classical time-domain SAFT algorithm. For a finite aperture source, it is clearly shown that the resolution obtained using the proposed method is superior to that obtained using the SAFT algorithm.     

平面近场声全息中正则化参数的确定

近场声全息的逆向重建过程属于线性病态逆问题,必须进行正则化处理。本文对三种基于Tikhonov正则化的参数选择方法,即离差原理法、广义交叉验证法、L曲线法,在不同全息距离、声源频率和信噪比的条件下进行了比较,结果表明,它们在远距离及低噪声环境下难以获得合适的正则化参数。采用等效噪声方差的方法,对其中较为稳定的离差原理进行了改进,使其在较远全息距离及低噪声环境下仍能获得合适的正则化参数。相应的仿真实验表明,改进后的离差原理法在很宽的信噪比(＞6 dB)和较远的全息距离(~10 cm)均能稳定地找到合适的正则化参数。此外,由于该方法无须对全息声压进行平滑处理,其有效重建孔径和全息孔径相等。      

定向耦合法测量非标准接口元件驻波系数

 利用网络分析仪，根据需要在设定的频率范围进行扫频测量，网络分析仪的端口1输出的微波信号经过微波传输转换器进入与待测非标准元件匹配连接的连接段（如波导）；网络分析仪的端口2连接定向耦合器，用于监测传输段内的反射微波信号。首先，在非标准接口端面连接短路面，通过网络分析仪测量反射信号在设定频点的相对幅度值；然后去掉短路面，在非标准接口端面连接待测元件，再次测量反射信号在特定频点的相对幅度值；最后根据本文推导的公式得出驻波系数。该方法的测量误差与定向耦合器的方向性（方向性系数越大越好）和待测元件驻波系数有关。将该方法运用在非标准接口匹配负载驻波系数的测量中，定向耦合器方向性系数取为40 dB，测得其驻波系数小于1.2，误差小于20%。该方法简便可行，可以用于测量常用的非标准接口元件尤其是非标准低驻波系数元件的驻波系数。     

Method for designing arbitrary two-dimensional continuous phase elements

We propose a method for designing arbitrary two-dimensional continuous phase elements. The phase anomalies that are due to amplitude zeros that occur in conventional design methods can be avoided completely with this method. Therefore the method is also useful for designing holographic elements that are free from amplitude zeros. The method is based on an iterative Fourier-transform algorithm. In the iteration process the method operates on the unwrapped phase. Consequently, a subsequent continuation step is not needed.     

能量色散X射线荧光分析中改进型基本参数法研究

能量色散X射线荧光分析方法是目前常用的一种多元素分析方法,但该方法检出限和分析精度,受到分析基体的影响。基本参数法是目前一种常用的分析方法,但在使用过程必须获取净峰面积和基体所有成分,而在实际使用时,尤其在分析低含量样品时,净峰面积计算、基体中“暗物质”影响了测量精度,制约了基本参数法的应用。针对基本参数法的不足,将谱线解析方法与基本参数法融合,将重叠峰剥离过程嵌入基本参数法迭代过程中。在含量计算过程中,采用分析样品特征X射线分支比的理论系数,对重叠峰进行剥离,解决能量色散X射线荧光测量中净峰面积计算和定量分析问题;在计算过程中,对“暗物质”进行均一化处理。通过对标准物质测量分析,结果表明对于Ni,Cu,Zn三个元素改进型基本参数法(改进型FP)测量结果准确度高于影响系数法。     

An application of a transcorrelated procedure to the calculation of the energy of the hydrogen molecule

The ground-state energy of the hydrogen molecule is calculated using both the standard variational method and a new transcorrelated method. The value obtained with the new method compares favourably with that obtained by the standard method. If the method used in this paper can be extended to more complicated systems, it is expected that it will enable more accurate energies to be calculated for these systems than is possible using the standard variational method.     

聚合物驱试井分析叠加原理适用性的探索

建立考虑聚合物剪切变稀特性的聚驱试井数学模型,基于矩形网格和井筒周围径向加密的复合网格.采用有限体积方法求得数值解,将数值解与叠加原理得到的注聚井停注压降进行对比.结果表明:叠加原理计算得到的注聚井关井井底压力值远低于数值解,说明叠加原理不能用于聚合物驱试井解释.在实测注聚井压降资料解释中,叠加原理得到的渗透率明显偏低...     

Multiquadric Finite Difference (MQ-FD) Method and Its Application

The conventional finite difference (FD) schemes are based on the low order polynomial approximation in a local region. This paper shows that when the polynomial approximation is replaced by the multiquadric (MQ) function approximation in the same region, a new FD method, which is termed as MQ-FD method in this work, can be developed. The paper gives analytical formulas of the MQ-FD method and carries out a performance study for its derivative approximation and solution of Poisson equation and the incompressible Navier-Stokes equations. In addition, the effect of the shape parameter $c$ in MQ on the formulas of the MQ-FD method is analyzed. Derivative approximation in one-dimensional space and Poisson equation in two-dimensional space are taken as model problems to study the accuracy of the MQ-FD method. Furthermore, a lid-driven flow problem in a square cavity is simulated by the MQ-FD method. The obtained results indicate that this method may solve the engineering problem very accurately with a proper choice of the shape parameter $c$.