复杂介质地震初至波数值模拟

发展了一种地震初至波(first break)数值模拟的旅行时插值波前追踪方法.将复杂变速介质剖分成匀速的矩形或多边形单元,在单元边界上设置一些对波前旅行时采样的节点.首先从震源开始计算波前到达所有节点的旅行时,然后从接收点开始反向确定震源与接收点之间的射线路径.其中,到达任意一点的波前旅行时和射线是通过对该点所在单元其它两两相邻的已算节点旅行时的插值和Fermat原理求极小而获得的.该算法对介质的复杂程度、单元形状和震源与接收点的位置没有限制,能模拟直达波、临界折射波和绕射波以及盲区射线等,具有较强的适应性和较高的精度.用该方法对一些典型近地表模型的初至波进行了数值计算,清晰地显示出这些模型的波前形态和射线路径.     

Three-dimensional shape measurement for an underwater object based on two-dimensional grating pattern projection

In this paper, a practical method using phase tracking and ray tracing algorithms is proposed for measuring the three-dimensional (3D) shape of an underwater object. A 2D projected sinusoidal fringe goes through the water and illuminates the tested object. Firstly, the phase tracking algorithm is employed to identify homologous points in phase distributions of the deformed fringe captured by the camera and these of the fringe pattern projected by the projector. The projector is regarded as a special camera as regards the stereovision principle. In the calibrated system, both ray directions of the homologous points can be easily figured out. Secondly, the ray tracing algorithm is used to trace the propagation path of each ray and to calculate the 3D coordinates of each point on the tested object's surface. Finally, the whole shape of the tested object can be reconstructed.     

Determination of equivalent sound speed profiles for ray tracing in near-ground sound propagation

The determination of appropriate sound speed profiles in the modeling of near-ground propagation using a ray tracing method is investigated using a ray tracing model which is capable of performing axisymmetric calculations of the sound field around an isolated source. Eigenrays are traced using an iterative procedure which integrates the trajectory equations for each ray launched from the source at a specific direction. The calculation of sound energy losses is made by introducing appropriate coefficients to the equations representing the effect of ground and atmospheric absorption and the interaction with the atmospheric turbulence. The model is validated against analytical and numerical predictions of other methodologies for simple cases, as well as against measurements for nonrefractive atmospheric environments. A systematic investigation for near-ground propagation in downward and upward refractive atmosphere is made using experimental data. Guidelines for the suitable simulation of the wind velocity profile are derived by correlating predictions with measurements.     

Analytic ray curve tracing for outdoor sound propagation

Outdoor sound propagation, which propagates sound through inhomogeneous, moving media with complex obstacles, presents challenging scenarios for computational simulation. In this paper, we present a ray-tracing method that uses analytic ray curves as tracing primitives in order to improve the efficiency of outdoor sound propagation in fully general settings. This ray-curve tracer inherits the efficiency and flexibility of rectilinear ray tracers in handling boundary surfaces, and it overcomes the performance limitations imposed by approximating the curved propagation paths in inhomogeneous media with rectilinear rays. Adaptive media traversal, as well as acceleration structures for surfaces intersections, lead to further savings in computation. Our method’s speedup over existing ray models, at least an order of magnitude for simple 2D scenarios and up to two orders of magnitude for 3D complex scenes, is demonstrated on outdoor benchmark scenes.     

An iterative algorithm for computing aeroacoustic integrals with application to the analysis of free shear flow noise

An iterative algorithm is developed for the computation of aeroacoustic integrals in the time domain. It is specially designed for the generation of acoustic images, thus giving access to the wavefront pattern radiated by an unsteady flow when large size source fields are considered. It is based on an iterative selection of source-observer pairs involved in the radiation process at a given time-step. It is written as an advanced-time approach, allowing easy connection with flow simulation tools. Its efficiency is related to the fraction of an observer grid step that a sound-wave covers during one time step. Test computations were performed, showing the CPU-time to be 30 to 50 times smaller than with a classical non-iterative procedure. The algorithm is applied to compute the sound radiated by a spatially evolving mixing-layer flow: it is used to compute and visualize contributions to the acoustic field from the different terms obtained by a decomposition of the Lighthill source term.     

Numerical simulation of transit-time ultrasonic flowmeters: uncertainties due to flow profile and fluid turbulence

Flowmeter measurement using the ultrasonic transit-time method is based on the apparent difference of the sound velocity in the flow direction and in the opposite direction. This method gives a flow velocity averaged along a particular acoustical path. To convert this path velocity to a velocity averaged over the entire cross-section of the flowing medium, the knowledge of the flow velocity profile is essential. However, the acoustical paths joining the two transducers are supposed to be straight and fluid turbulence phenomena are neglected. In this paper, we describe a numerical procedure to estimate the uncertainties due to these approximations in the case of fully developed turbulence. The ultrasonic propagation is modelled in 2-D moving inhomogeneous media via a ray tracing algorithm. Influence of mean profiles of temperature and velocity is studied on simple examples. Fluid temperature fluctuations and fluid velocity turbulence are considered in the stochastic framework to obtain average uncertainties on the measurements of the liquid flow rate.     

基于射线理论的各向异性煤层地震响应

以VTI型煤层为研究对象,从介质中的地震波的射线方向、走时计算等多个方面展开讨论,分析VTI型煤层的试射法射线追踪技术,充分利用VC++的面向对象类框架泛型编程的特点和Fortran数值计算方面的优势,通过混合编程实现VTI型煤层中的地震波试射法射线追踪;通过各向同性煤层和VTI型煤层的射线追踪试算对比,论证对VTI型煤层进行各向异性射线追踪的必要性,同时验证本文VTI型煤层射线追踪方法的正确性.最后,通过计算时间成本比较,表明该方法计算速度较快,能满足地震正演模拟的要求.该方法适用于各向异性煤层的多波射线追踪.     

Analysis of plasma sheath propagation attenuation based on ray tracing

Based on the complex Snell's law of lossy media, a ray tracing method is proposed to study the propagation attenuation characteristics of electromagnetic (EM) waves in plasma sheaths. The plasma sheath is modelled as layered media. This method considers the complex ray characteristics of inhomogeneous plane EM waves, tracks the propagation rays of EM waves in each layer of media, and calculates the propagation attenuation of EM waves in each layer of media according to the propagation direction of the complex rays. The attenuation during numerical cumulative propagation is the total attenuation of EM waves through the plasma sheath. By comparing the results with that obtained from the WKB method, the accuracy of the ray tracing algorithm is proved. The results of the propagation attenuation of a blunt cone model are calculated by the proposed method, and the effects of different parameters on the EM wave propagation attenuation in the plasma sheath are analysed at different heights, velocities, incident angles, and incident positions. Studying the propagation characteristics of EM waves in the plasma sheath is of importance in application for radar target tracking, blackout communication, and other issues.     

Transient coupled radiative and conductive heat transfer in a semitransparent layer of ceramic

This work considers transient radiative and conductive heat transfer in a semitransparent layer of ceramic, submitted to several thermal and radiative boundary conditions. Each side of the layer is exposed to hot or cold radiative surroundings, while each boundary is heated or cooled by convection. The solution procedure must provide accurate temperature distribution in the layer, so a nodal analysis based on Hottel's zonal method extended by ray tracing method is carried out. A finite difference method with non-uniform space and time increments is used to solve the transient energy equation, including a radiative heat source, coupled to a equation of radiative transfer. Variable spacing was used to concentrate grid points in regions with large temperature gradients. The influence of refractive index, optical thicknesses and conduction-radiation parameters is investigated.     

海洋声学环境下水中声源的时延估计法定位精度分析

水中声源的定位精度受到海洋声学环境的重要影响。结合海上试验的实际应用,分析了水下观测平台采用时延估计法对声源的定位精度问题。根据理论分析,计算了时延估计误差、海洋中声速不均匀、平台非稳性、及声传播起伏等因素引起的俯仰角和方位角误差。利用误差传递公式,获得了上述因素引起的不同平台深度下,不同距离声源的定位误差。比较了采用平面阵与立体阵、是否补偿声线弯曲效应等条件下定位误差的变化,并通过海上试验结果进行了部分验证。研究结果表明,海洋声速不均匀对定位误差的贡献最大。采用立体阵代替平面阵、测量海洋声速剖面并补偿声线弯曲引起的定位误差,在1000m距离上可使定位相对误差从最大30%降低到约10%,有效提高了较远距离上的定位精度。研究结果对于采取措施提高水中声源的定位精度有指导意义。      

一种考虑声线弯曲的温度场重建算法

声波在非均匀温度场中传播时会由于声波的折射而产生声线弯曲现象。为提高非均匀温度场声学CT (Computer Tomography)重建精度,提出一种考虑声线弯曲的重建算法。首先用最小二乘法获得一个不考虑声线弯曲的重建温度场,然后用打靶-插值法确定本征声线出射角,用三角形前向展开法追踪声线,获得此温度场中声波发射/接收器间的本征声线,建立本征声线上声波传播时间与温度分布间的关系,进而实现考虑声线弯曲的最小二乘法温度场重建。单峰、双峰温度场仿真重建结果表明:所采用的本征声线追踪法计算简单运行速度快;考虑声线弯曲后,可明显提高温度场重建精度。因此所提重建算法能更好地适应实际温度场重建需求。      

Modelling the interior sound field of a railway vehicle using statistical energy analysis

The sound field in train compartments, treated as a series of connected air cavities, is modelled using statistical energy analysis, SEA. For the case under study, with five cavities in series and the source in the second cavity, a closed-form solution is obtained. An adjusted SEA model is used to predict the rate of spatial decay within a cavity. The SEA model is validated using results from a ray tracing method and from scale model measurements. For the octave bands 500–4000 Hz, good agreement is shown between the results from measurements, the ray tracing and the SEA model, for the two saloons closest to the source cavity (a vestibule). The SEA model was shown to slightly underestimate the rate of spatial decay within a cavity. It is concluded that a reasonable cause is the additional diffusion due to the seating.     

环形子孔径拼接检测中机械误差的分离

为减少环形子孔径拼接干涉检测中机械误差对检测结果造成的影响,分析环形子孔径拼接过程中机械误差作用分量的表现形式,提出了分离机械误差的全局优化的环形子孔径拼接方法。分析根据波像差理论建立的机械误差分离数学模型,然后将其应用于避免误差传递和累积的全局优化的拼接方法中,并提出利用光线追迹的方法在拼接之前除去理想非球面波前与参考球面波前的差别。应用分离机械误差的拼接方法对口径为75mm、顶点曲率半径为100 mm的抛物面面形进行检测,得到的面形峰谷值误差为0.05,均方根值误差为0.003,验证了该拼接方法可有效分离环形子孔径拼接中的机械误差。     

用射线踪迹法解黑表面矩形介质耦合换热

本文用射线踪迹-节点分析法研究了二维黑体表面矩形、各向同性散射半透明介质内辐射与导热瞬态耦合换热。采用全隐格式的有限差分法离散二维瞬态微分能量方程,用辐射传递系数来表示辐射源项,结合谱带模型并采用射线踪迹法求解辐射传递系数。采用Patankar线性化方法将辐射源项及不透明边界条件线性化,并采用附加源项法处理边界条件,运用ADI方法求解名以上的线性化方程组,从而解得二维矩形介质内的瞬态温度分布。     

光栅波像差对脉冲压缩的影响

 采用光线追迹的方法,对脉冲压缩光栅系统波像差进行理论研究,构建了有像差脉冲压缩光栅对分析模型,使用该模型通过模拟计算得出了波像差对脉冲压缩效果的影响。理论分析结果表明：脉冲压缩光栅系统最大波像差不大于λ/4时,脉冲压缩比的最大变化率约为2%。脉冲压缩光栅对系统达到衍射极限时,在选择适当的入射角度时,光栅波像差所引起的脉冲压缩比变化并不大。     

奇对称型相位板波前编码成像系统的点列图分析

基于几何光学理论利用光线追迹方法得到了适用于任意面型相位板的大视场平行光入射情况下的光线像差近似表达式,并在此基础之上分析了奇对称型相位板波前编码系统的点列图的大小、边界以及光线结构等特性,研究了大视场角平行光入射下的光线像差增强及点列图形变等现象,并且给出了相应的近似数学解析式来描述这些特性。通过对奇对称型相位板波前编码成像系统的点列图特性分析,有助于进一步的理解波前编码技术并且指导实际系统的设计。     

Identification of radiative properties of reticulated ceramic porous inert media using ray tracing technique

The radiative properties of reticulated porous inert media are computationally identified using the real three-dimensional structural data of porous media. The computational grids data are reconstructed from three-dimensional computer tomography scans and magnetic resonance image scans of different reticulated porous media. A ray tracing algorithm is used to track the rays inside the grid structure. Statistically large numbers of rays are traced for their path length and incident angle, which are used to find the probability based equivalent extinction coefficient and scattering phase function. The equivalent extinction coefficients are found for porous media with different porosities and pore densities. The dependency of specular and diffuse scattering phase functions on the porous structure and surface reflectance are also studied.     

An efficient method for the incompressible Navier–Stokes equations on irregular domains with no-slip boundary conditions,high order up to the boundary

Common efficient schemes for the incompressible Navier–Stokes equations, such as projection or fractional step methods, have limited temporal accuracy as a result of matrix splitting errors, or introduce errors near the domain boundaries (which destroy uniform convergence to the solution). In this paper we recast the incompressible (constant density) Navier–Stokes equations (with the velocity prescribed at the boundary) as an equivalent system, for the primary variables velocity and pressure. equation for the pressure. The key difference from the usual approaches occurs at the boundaries, where we use boundary conditions that unequivocally allow the pressure to be recovered from knowledge of the velocity at any fixed time. This avoids the common difficulty of an, apparently, over-determined Poisson problem. Since in this alternative formulation the pressure can be accurately and efficiently recovered from the velocity, the recast equations are ideal for numerical marching methods. The new system can be discretized using a variety of methods, including semi-implicit treatments of viscosity, and in principle to any desired order of accuracy. In this work we illustrate the approach with a 2-D second order finite difference scheme on a Cartesian grid, and devise an algorithm to solve the equations on domains with curved (non-conforming) boundaries, including a case with a non-trivial topology (a circular obstruction inside the domain). This algorithm achieves second order accuracy in the L^∞ norm, for both the velocity and the pressure. The scheme has a natural extension to 3-D.     

Model-based adaptive non-null interferometry for freeform surface metrology

A model-based adaptive non-null interferometry(MANI) is proposed for steep optical freeform surfaces in situ testing. The deformable mirror(DM) affording the flexible compensation is monitored with the beam in the interferometer by a wavefront sensor. The residual wavefront aberration in the non-null interferogram is eliminated by the multi-configuration ray tracing algorithm based on the system model, especially the DM surface model. The final figure error can be extracted together with the surface misalignment aberration correction.Experiments proving the feasibility of the MANI are shown.