有限元/边界元耦合方法在海洋生物目标强度特性研究的应用*

目标强度特性是海洋生物声学识别与资源量评估的重要依据,其中,基于近似几何体和声阻抗特性的理论模型法是研究海洋生物目标强度的重要手段。由于对几何形态近似处理以及数值求解方法的限制,传统理论模型对声波频率、入射方位以及目标声阻抗、形态尺寸等均有各自不同的适用范围,单一模型难以满足不同种类或同一种类但不同尺寸海洋生物的目标强度求解。本文尝试将逐渐见诸应用的有限元/边界元耦合方法用于海洋生物目标强度特性研究,分别以球形生物、纺锤形鱼类尾明角灯鱼(Ceratoscopelus warmingii)和细长形浮游动物南极大磷虾(Euphausia superba)为例进行仿真计算,并与相适应的经典理论模型进行对比分析。结果表明,对于球形生物,有限元/边界元耦合方法与解析模型的目标强度频响曲线完全吻合;对于纺锤形鱼类,有限元/边界元耦合方法可有效弥补基于模态级数解的形变圆柱体模型在中低频和两端入射时的准确性问题;对于细长形浮游动物,有限元/边界元耦合方法与畸变波玻恩近似模型高度吻合。综上,有限元/边界元耦合方法对多种海洋生物目标强度求解均具有较好的适用性,未来有待进一步结合实验测定进行验证。     

Comparison of an analytic horn equation approach and a boundary element method for the calculation of sound fields in the human ear canal

The sound field inside a model human ear canal has been computed, to show both longitudinal variations along the canal length and transverse variations through cross-sectional slices. Two methods of computation were used. A modified horn equation approach parametrizes the sound field with a single coordinate, the position along a curved center axis-this approach can accommodate the curvature and varying cross-sectional area of the ear canal but cannot compute transverse variations of the sound field. A boundary element method (BEM) was also implemented to compute the full three-dimensional sound field. Over 2000 triangular mesh elements were used to represent the ear canal geometry. For a plane piston source at the entrance plane, the pressure along the curved center axis predicted by the two methods is in good agreement, for frequencies up to 15 kHz, for four different ear canals. The BEM approach, though, reveals spatial variations of sound pressure within each canal cross section. These variations are small below 4 kHz, but increase with frequency, reaching 1.5 dB at 8 kHz and 4.5 dB at 15 kHz. For source configurations that are more realistic than a simple piston, large transverse variations in sound pressure are anticipated in the vicinity of the source.     

大倾斜角度下基于冷原子重力仪的绝对重力测量

冷原子重力仪的倾斜角会对绝对重力测量产生显著的影响.高精度的绝对重力测量需要对重力仪的倾斜角进行精确的测量、控制及校正.本文从理论上分析了四种不同情况下倾斜对绝对重力测量的影响规律,并在实验上对得到的理论进行了实验验证.基于此,设计了一种基于双倾斜计的绝对重力测量方案,主要是为了解决恶劣测量环境下的冷原子重力仪倾斜漂移问题.此方案利用高精度倾斜计记录放置在被动隔振平台上的拉曼反射镜的倾斜角度,并使用另外一个倾斜计监控真空系统的倾斜,以实现振动噪声的抑制和倾斜的高精度测量.基于自研的小型化冷原子重力仪,对该方案进行了实验验证,并最终实现了车间复杂环境下的高精度绝对重力测量.由于倾斜得到精确测量和补偿,冷原子重力仪的测量精度达到了12.3μGal.本文为复杂环境下的高精度绝对重力测量提供了一种可行的方案,为冷原子重力仪的实用化提供了参考数据.     

Nanotribology of self-assembled monolayer with a probe tip investigated using molecular dynamics simulations

The nanotribology of an alkanethiol self-assembled monolayer (SAM) under tilt contact with a scanning probe tip is studied using molecular dynamics (MD) simulations. The tilt contact is described in terms of the tilt angle and the magnitude of the specimen–tip separation. The effects of tilt angle and magnitude of the specimen–tip separation on the normal force, friction force, friction coefficient, shear strength of the tip–SAM junction, and self-recovery characteristics are evaluated during the scanning probe tip process at a temperature of 300 K. The simulation results clearly show that the magnitudes and periods of the normal force and friction force increase with decreasing magnitude of the specimen–tip separation due to a large change of the tilt angle of the SAM chains during the deformation and recovery stages. For scanning and indentation processes, the effect of the tilt angle of the probe tip on the normal force is more significant than that on the friction force for the SAM. The behaviors of interfacial contact forces, friction coefficient, and shear strength strongly depend on the number of interacting atoms and the contact area, which increases with decreasing magnitude of the specimen–tip separation and increasing tilt angle of the probe tip. The self-recovery of SAM is significantly affected by the magnitude of the specimen–tip separation; the recovery ability of SAM is worse for magnitude of the specimen–tip separation below −0.9 nm with a large tilt angle of the probe tip.     

Deformation study of bicrystalline and nano-polycrystalline structures using phase field crystal method

Deformation behaviors of bicrystalline and nano-polycrystalline structures of various tilt angles and inclination angles in two dimensions are investigated in detail using a two-mode phase field crystal model.The interaction between grain boundary(GB)and dislocation is also examined in bicrystals and nano-polycrystals that both contain asymmetric and symmetric tilt GBs,with energy analysis being carried out to analyze these processes.During deformation simulations,we assume the volume of each simulation cell at every time step is coincident with that of the initial state just before deformation.Our simulation results show that the behaviors of symmetric and asymmetric GBs in bicrystals and nano-polycrystals differ from each other depending on tilt angle and inclination angle.A new dislocation emission mechanism of interest is observed in bicrystals which contain low angle symmetric tilt GBs.Low angle GB has a higher mobility relative to high angle GB in both bicrystalline and nano-polycrystalline structures,as does asymmetric GB to symmetric GB.The generation,motion,pileup and annihilation of dislocations,grain rotation and grain coalescence are observed,which is consistent with the simulation results obtained by molecular dynamics.These simulation results can provide strong guidelines for experimentation.     

Acoustic remote sensing of swimbladder orientation and species mix in the oreo population on the Chatham Rise

A method for combining in situ measurements and theoretical swimbladder-derived estimates of target strength of the deep-water fish, black and smooth oreos, is described. The technique uses Monte Carlo simulation and yields fish length-target strength relationships suitable for use in estimating biomass from echo integration acoustic surveys. The relationships are derived from estimates of the mean and standard deviation of the tilt angle distributions of the wild fish generated by the method. The relationships may also be used to estimate proportions of the two oreo species in the wild. The mean tilt angle of black oreos in the wild was about 10 degrees with a standard deviation of 8 degrees. For smooth oreos it was close to zero with a standard deviation of about 4 degrees. The target strength relationships derived for biomass estimation purposes were TSB= -78.05+25.3 log10(L)+1.62 sin(0.0815L+0.238) and TSS = -82.16+24.6 log10(L)+1.03 sin(0.1165L-1.765), where L is the fish length and TSB and TSS are the target strengths of black and smooth oreos respectively.     

Effects of non-flat contact and interference on self-assembled monolayers under sliding friction

The nanotribology mechanism of alkanethiol self-assembled monolayers (SAM) chemisorbed on a gold surface under a non-flat contact by a tilt plane was studied using molecular dynamics (MD) simulations. The molecular trajectories, tilt angles, normal forces, shear forces, and frictional coefficient of the SAM were evaluated during the friction and relaxation processes for various parameters, including the tilt angle of the slider, interference magnitude, and SAM length. At the nanoscale, the magnitude of interface interactional forces is strongly dependent on the magnitude of the contact area, not on the surface geometry. The contact area and the exerted normal force of the SAM increase with decreasing the tilt angle of the slider at the same contact interference. In contrast, the periods in both normal force and shear force are gradually delayed as the tilt angle of the slider increases. Once the contact interference increases, the normal force and shear force increase together. During the sliding friction process with a smaller tilt slider angle, SAM molecules can maintain a better collective ordered structure. Short SAM molecules are more sensitive to a compressive loading and react to a larger normal force under the same contact interference due to the deformation of a larger tilt angle and decrease in chain length. The friction coefficient of SAM is significantly more dependent on the tilt angle of the slider than the contact interference.     

Demonstration for a two-axis interferometric tilt sensor in KAGRA

Recently, a folded Mach–Zehnder interferometer with homodyne in- and quadrature-phase detection was proposed as a high-precision, wide-dynamic range tilt sensor. By way of a practical application and to validate actual performance, two-axis tilt sensors were developed and installed for one mirror of the input mode cleaner cavity in KAGRA, the large-scale cryogenic gravitational-wave telescope in Kamioka, Gifu, Japan. Building on previous work, we have demonstrated that the two-axis tilt sensor has properly sensed the tilt angle changes of the mirror motion with high precision and without calibration. Compared with our initial angular sensor, an optical lever, which is calibrated by using the interferometer tilt sensor, we found that both sensors showed actual tilt motions of the mirror at low frequencies, and the two-axis interferometer sensor has a better sensitivity at higher frequencies.     

基于倾斜阵的距离与方位联合估计*

针对阵列倾斜引起的阵不变量方法定位精度下降的问题,该文结合倾斜阵声源定位的三维模型提出一种基于倾斜阵的距离与方位联合估计算法。该算法通过联合自校正方位距离联合估计方法对阵列倾角在声源-接收平面的投影量进行估计和补偿,改善目标距离估计误差的同时,利用匹配出的阵列倾角投影量反向估计目标方位,获得对声源二维位置的估计。通过仿真验证了该算法能够在阵列倾角较小时实现目标方位与距离的联合估计。     

Computer calculations of refractive indices in orthoconic SmC<Subscript>A</Subscript><Superscript>*</Superscript>. Is it possible to obtain “isotropic” antiferroelectric liquid crystal (IAFLC)?

Orthoconic antiferroelectric liquid crystals (OAFLC) are promising materials for applications. The saturated value of the tilt angle for such materials is 45°, while for “classical” antiferroelectric liquid crystals (AFLC) it reaches typically 30°. Homogeneously oriented slab of OAFLC, unwound by surface action of cell walls is apparently optically isotropic for the light beam at normal incidence. Such unwound structure of OAFLC one can call isotropic OAFLC (IOAFLC). Optical indicatrix of OAFLC is represented by an oblate sphere (ellipsoid of revolution) with its short axis perpendicular to cell walls. Here the question arises: “Is it possible to create isotropic AFLC (IAFLC)? An effective refractive index of IAFLC should have the same value for the light beam propagating at each direction, hence the optical indicatrix should be a sphere. In the way of theoretical considerations we found that IAFLC can exist for twisted structure of AFLC characterized by the tilt angle as high as 54.3°. Unfortunately, such a big tilt angle is still not achievable in practice. Calculations show that AFLC can be optically isotropic when the helix is not suppressed by cell walls action.     

Backscattering enhancements for tilted solid plastic cylinders in water due to the caustic merging transition: observations and theory

Bulk shear and longitudinal waves give rise to important contributions to the scattering of ultrasound by tilted finite plastic and rubber cylinders in water. This occurs in situations where either the shear or longitudinal speed is less than the speed of sound in the surrounding water. At a certain critical tilt angle, large backscattering enhancements are observed for finite cylinders, where the wave vector can reverse direction upon reflection from the cylinder truncation. The scattering process is analogous to the enhancement produced by the merging of rainbow caustics of primary rainbow rays in the scattering of light by long dielectric cylinders, also known as the caustic merging transition [C. M. Mount, D. B. Thiessen, and P. L. Marston, Appl. Opt. 37, 1534-1539 (1998)]. A ray theory was developed to model the backscattering mechanism at the critical tilt angle. It employs the idea of the Bravais effective refractive index, convenient for constructing ray diagrams for the projections of rays in the base plane of the cylinder. There is general agreement between the theory and the experiment down to relatively low ultrasonic frequencies (ka as small as 10). The enhancement is the most significant backscattering contribution for a wide range of tilt angles.     

Visualization of wave propagation in muffler

This paper presents a prediction technique for wave propagation in muffler using Boundary Element Method (BEM). The results of the numerical calculation are compared with the plane-wave theory and also the experimental results. It is shown that at high frequencies the plane-wave theory does not yield correct results due to the occurrence of two- and three-dimensional wave motion. On the other hand, the BEM is able to predict the acoustic wave behavior in the high frequencies where the planewave theory fails.     

Vibro-acoustic design sensitivity analysis using the wave-based method

Conventional element-based methods, such as the finite element method (FEM) and boundary element method (BEM), require mesh refinements at higher frequencies in order to converge. Therefore, their applications are limited to low frequencies. Compared to element-based methods, the wave-based method (WBM) adopts exact solutions of the governing differential equation instead of simple polynomials to describe the dynamic response variables within the subdomains. As such, the WBM does not require a finer division of subdomains as the frequency increases in order to exhibit high computational efficiency. In this paper, the design sensitivity formulation of a semi-coupled structural-acoustic problem is implemented using the WBM. Here, the results of structural harmonic analyses are imported as the boundary conditions for the acoustic domain, which consists of multiple wave-based subdomains. The cross-sectional area of each beam element is considered as a sizing design variable. Then, the adjoint variable method (AVM) is used to efficiently compute the sensitivity. The adjoint variable is obtained from structural reanalysis using an adjoint load composed of the system matrix and an evaluation of the wave functions of each boundary. The proposed sensitivity formulation is subsequently applied to a two-dimensional (2D) vehicle model. Finally, the sensitivity results are compared to the finite difference sensitivity results, which show good agreement.     

一维反铁磁光子晶体光学双稳态效应研究

基于传递矩阵方法和光局域化原理, 本文研究了一维反铁磁光子晶体共振频率附近光学双稳态效应随电磁波入射角、 外磁场强度及电介质层数的变化.研究发现, 当外磁场为1.0 kG(1 G=10^-4 T)、电磁波小角度入射时, 反铁磁材料高低共振频率附近均可探测到光学双稳态效应; 当电磁波大角度入射时, 仅在高共振频率附近探测到光学双稳态效应. 然而, 当外磁场强度增加到2.0 kG时, 由于反铁磁材料的高低共振频率向两侧移动, 低共振频率附近缺失的光学双稳态频率窗口又被有效激发. 此外, 电磁波小角度入射时, 电介质层数在低共振频率附近对双稳态效应影响较明显.     

低能Cl^-在Al2O3绝缘微孔膜中的输运过程

研究了10 keV Cl^-离子穿越Al2O3绝缘微孔膜的物理过程,发现穿越的Cl^-其分布中心在初束中心即0°附近,Cl^-离子穿透率下降与几何穿透一致,这是典型的直接几何穿越有一定角发散的微孔导致的结果;而出射的Cl0和Cl^+以微孔轴向为中心分布,Cl^+和Cl0穿透率下降慢于几何穿透.模拟计算发现沉积电荷会使出射粒子中Cl^-占主要成分,并使出射Cl^-角分布中心移动到微孔轴向方向而随微孔膜倾角移动;而在不考虑沉积电荷的情况下,计算结果较好地符合了实验结果.通过分析在不同倾角下散射过程对出射粒子的角分布和电荷态分布的影响,发现绝大部分的Cl0是通过一次和两次散射出射的,其中一次散射出射的Cl0占主要成分,从而导致出射的Cl0沿微孔轴向出射而Cl^+主要是经过一次碰撞出射.这导致了随倾角增大,出射的Cl0穿透率下降速度比Cl^+小,Cl0所占比例相对增大较快,从而导致观测到的Cl^+/Cl0的比例下降.本文结果更仔细地描述了低能离子穿越绝缘体微孔的物理机理,印证了之前实验和理论工作的结果,发现在10 keV以上能区的Cl^-离子穿越绝缘微孔膜的过程中,沉积电荷并未起到主要作用,其主要穿透特征是散射过程造成的.     

Deriving correction functions to model the efficiency of noise barriers with complex shapes using boundary element simulations

The boundary element method (BEM) is a commonly used method to compute the insertion loss of noise barriers having arbitrary cross-sections. For large scale three-dimensional problems, however, the BEM is not feasible. On the other hand, standardized calculation methods for noise mapping are efficient, but shapes other than the straight barrier cannot be properly calculated. Attempts to merge these two approaches by using BEM to derive correction functions based on geometrical quantities such as source and target angle as well as the path length elongation between source and receiver caused by the barrier were usually focused on a small set of barrier types, dimensions, absorptive configurations, source or receiver positions. The main objective of this study is to investigate which functions based on the most common geometrical parameters are well suited for approximating the efficiency of different types of barriers, dimensions and absorptive configurations. To achieve this, numerous combinations of 7 different barrier types, different heights and widths as well as 3 different absorptive configurations were simulated using the 2D BEM for 8 different source positions. The octave-band-wise efficiency, i.e. the frequency-dependent gain in insertion loss compared to an equally high, fully reflective straight barrier was used as a basis for the correction functions. Linear as well as polynomial models were compared yielding a polynomial of third degree in the source and fourth degree in the target angle as the best model. Effects on the error using uniform sampling in the target angle instead of a uniform receiver grid as a basis for the correction functions are also investigated. Furthermore, wide-band efficiencies based on standardized traffic emission spectra are calculated showing small errors compared to single-band errors, in particular in the high-frequency range. A linear interpolation scheme is suggested to deal with barriers having dimensions not simulated in this work.     

Spectral integration and wideband analysis in gap detection and overshoot paradigms

Several listening conditions show that energy remote from a target frequency can deleteriously affect sensitivity. One interpretation of such results entails a wideband analysis involving a wide predetection filter. The present study tested the hypothesis that both temporal gap detection and overshoot results are consistent with a wideband analysis, as contrasted with statistical combination of information across independent channels. For gap detection, stimuli were random or comodulated 50-Hz-wide noise bands centered on 1000, 1932, 3569, and 6437 Hz. For overshoot, the masker was an 8-kHz low-pass filtered noise, with 5-ms tone bursts presented at the same center frequencies used for gap detection. Signals were presented with either 0- or 250-ms delay after masker onset. In each paradigm, the target was introduced at only one frequency or at all four frequencies. Results from gap detection conditions did not favor a wideband analysis interpretation: Results in the random condition were consistent with an optimal combination of cues across frequency. An across-channel interference effect was also evident when only one of the four bands contained the gap. Although results from the overshoot conditions were consistent with a wideband analysis interpretation, they were more parsimoniously accounted for in terms of statistical combination of information.     

脉冲激光探测平面目标特性对测距分布的影响

针对脉冲激光近程探测系统平面目标特性影响测距分布的问题,推导了平面目标脉冲响应方程和脉冲激光回波方程,分析了不同平面目标倾斜角、激光发射发散角和激光发射脉宽对激光回波展宽的影响.基于脉冲激光回波方程和恒阈值时刻鉴别方法,推导了平面目标的测距概率密度函数解析式,并理论仿真分析了不同平面目标倾斜角、激光发射功率、激光发射发散角及阈噪比对测距统计特性分布的影响.运用蒙特卡罗算法进行全波形模拟测距实验;搭建脉冲激光测距实验环境,进行20 m测距实验.实验结果表明:理论仿真、蒙特卡罗模拟实验与实际实验的测距概率密度分布基本一致,随着平面目标倾斜角的增大,测距均值和测距方差增大;当倾斜角度为0°,20°,40°,60°时,回波信噪比高于阈噪比,测距分布呈现高斯分布;当倾斜角度为70°时,回波信噪比低于阈噪比,分布不再呈现高斯分布,呈现上升沿缓下降沿陡的分布特性.研究结果为研究目标平面特性对脉冲激光探测测距分布的影响提供了理论依据.     

非对称条形纳磁体的铁磁共振频率和自旋波模式

通过建立微波激励下的非对称条形多铁纳磁体的微磁模型,研究了倾斜角和缺陷角对该形纳磁体的铁磁共振谱和自旋波模式的影响.通过对微磁仿真得到的动态磁化数据进行分析发现,非对称条形纳磁体倾斜角度增加,铁磁共振频率随之增加,而这一现象与纳磁体的缺陷角度无关.倾斜角不变,非对称条形纳磁体的铁磁共振频率与缺陷角度呈单调递增关系,并且不同缺陷角度纳磁体的自旋波模式显示出极大的差异.非对称条形纳磁体与矩形纳磁体相比,它的自旋波模式局部化,具体为非对称条形纳磁体的自旋波模式不对称且高进动区域存在于边缘,表现为非对称边缘模式.倾斜角改变导致纳磁体内部退磁场变化,引起纳磁体边缘模式的移动,而中心模式对倾斜角的变化并不敏感.最后,对建立的模型在高频微波磁场激励下的磁损耗进行了分析,验证了模型的可靠性.这些结论说明缺陷角和倾斜角可用于纳磁体自旋波模式和铁磁共振频率的调谐,所得结果为可调纳磁微波器件的设计提供了重要的理论依据和思路.     

Direct measurement of tilt angle in smectic mesophases by neutron diffraction

We report the results of a direct measurement of the tilt angle of the molecules with respect to the normal layer in the different smectic phases and in the crystalline phase of the TBBA liquid crystal. The values of the tilt angles are obtained by using a new simple method that consists of direct reading of the angular distance between two Bragg peaks in a single neutron rocking curve. The results on the smectic phases are in good agreement with the existing experimental values obtained with other techniques.