南海北部海域声速剖面水平变化环境声学层析

为获取水平变化的海水声速分布,在等效平均海水声速反演基础上,提出了一种声速剖面水平变化环境声学层析方法。首先利用经验正交函数结合匹配场处理方法反演不同距离下的等效平均声速剖面,进而通过逐步迭代重构出随水平距离变化的声速场分布,并通过数值仿真及海上实验两方面验证了声速剖面水平变化环境声学层析方法的可行性及其准确性。结果表明:声速剖面水平变化环境声学层析结果与实验测量声速剖面均方根误差小于0.8 m/s,利用声传播信号基本能较为准确地重构出随距离水平变化的海水声速分布,为进一步层析中尺度三维结构及声呐水声环境保障应用奠定了基础。      

An accurately fast algorithm of calculating reflection/transmission coefficients

For the boundary between transversely isotropic media with a vertical axis of symmetry (VTI media), the interface between a liquid and a VTI medium, and the free-surface of an elastic half-space of a VTI medium, an accurately fast algorithm was presented for calculating reflection/transmission (R/T) coefficients. Specially, the case of post-critical angle incidence was considered. Although we only performed the numerical calculation for the models of the VTI media, the calculated results can be extended to the models of transversely isotropic media with a horizontal axis of rotation symmetry (HTI media). Compared to previous work, this algorithm can be used not only for the calculation of R/T coefficients of the boundary between ellipsoidally anisotropic media, but also for that between generally anisotropic media, and the speed and accuracy of this algorithm are faster and higher. According to the anisotropic parameters of some rocks given by the published literature, we performed the calculation of R/T coefficients by using this algorithm and analyzed the effect of the rock anisotropy on R/T coefficients. We used Snell’s law and the energy balance principle to perform verification for the calculated results. Supported by the Natural Science Foundation of Shaanxi Province, China (Grant No. 2007D15)     

Reflection and transmission of acoustic waves on multilayered porous media

Based on the Boit theory of acoustic wave propagation in fluid-satu-rated porous medium we have studied in this paper the acoustic reflection andtransmission on multilayered porous media,in which the adequate boundaryconditions across the interfaces are taken into account.Numerical calculationsof the reflection and transmission coefficients at different incident angles andfrequencies of the fast compressional wave incident on porous media with threeor four layers are presented.The results indicate that the maximum or mini-mum reflection and transmission coefficients appear at certain ratios of thewavelength to the thickness.The acoustic incident angle and porous mediumproperties are shown to affect significantly these coefficients.As an example,the measured transmission coefficients in a water-saturated fused glass beadsample are in good agreement with theoretical prediction.     

Interaction of ultrasonic wave with void-containing layer in solid

Using a general mode for sound reflection from multilayered media,we present in this paper the expressions for sound reflection and transmissioncoefficients on void-containing interface layer in solid and derive the character-istic equation for symmetric and antisymmetric modes of the interface wavesalong the layer.The method for evaluating the effective elastic moduli ofvoid-containing solid is also introduced.The numerical calculation given in thispaper shows the influences of the void volume concentration and layer thick-ness on the sound reflection coefficients and interface wave velocity,providing atheoretical basis for inverse of the mechanical properties of void-containinglayer based on the ultrasonic measurement.     

Matrix Riccati equation formulation for radiative transfer in a plane-parallel geometry

In this paper, we formulate the radiative transfer problem as an initial value problem via a pair of nonlinear matrix differential equations (matrix Riccati equations or MREs) which describe the reflection (R) and transmission (T) matrices of the specific intensities in a plane-parallel geometry. One first computes R and T matrices of some small but finite layer thickness (equivalent optical thickness τ∼0.01 and then repetitively applies the doubling method to the reflection and transmission matrices R(τ)and T(τ) until reaching the desired layer thickness. The initial matrices R(τ₀)and T(τ₀) can be computed quite accurately by either of the following methods: multiple-order, multiple-scattering approximation, iterative method or fourth-order Runge-Kutta techniques. In addition, the reflection coefficient matrix of a semi-infinite medium satisfies an algebraic matrix equation which can be solved repetitively by a matrix method. MREs offer an alternative way of solving plane-parallel radiative transport problems. This method requires only elementary matrix operations (addition, multiplication and inversion). For vector and/or beam-wave radiative transfer problems, large matrices are required to describe the physics adequately, and the MRE method provides a significant reduction in computer memory and computation time.     

柱面分层流体饱和孔隙地层中的声波测井波场模拟

因钻井作业时泥浆侵入等因素,井外介质沿径向常呈现非均匀性。本文将地层视为柱面分层孔隙介质,采用传递矩阵法模拟声波测井波场。注意到慢纵波波数的实部和虚部都很大,当井外含大厚度层且频率较高时,易因计算大自变量的Hankel函数而溢出,我们在波场表达式中引入了归一化的Hankel函数。论文将前人模拟关于柱面分层固体弹性介质声波测井的广义反射/透射系数计算方法推广到了分层孔隙介质情况,导出了波幅系数的计算公式,并针对含侵入带的柱面分层孔隙地层,计算了声波测井全波波形。      

超声与固体中含气泡层的相互作用

本文利用固体层状媒质声反射模型,给出了固体中含气泡层声反射和透射系数的表达式,并由此导出沿固体中含气泡薄层对称和反对称模式界面波的特征方程式。本文还介绍了含气泡固体有效弹性模量的估算方法,文中给出的数值计算具体说明了气泡体积浓度和气泡层厚度对声反射系数、声透射系数以及反对称模式界面波传播速度的影响,本文的研究为根据声反射系数和界面波的传播速度的测量反演固体间气泡层的力学性能提供了理论依据。 关键词：     

Permutation approach to finite-alphabet stationary stochastic processes based on the duality between values and orderings

The duality between values and orderings is a powerful tool to discuss relationships between various information-theoretic measures and their permutation analogues for discrete-time finite-alphabet stationary stochastic processes (SSPs). Applying it to output processes of hidden Markov models with ergodic internal processes, we have shown in our previous work that the excess entropy and the transfer entropy rate coincide with their permutation analogues. In this paper, we discuss two permutation characterizations of the two measures for general ergodic SSPs not necessarily having the Markov property assumed in our previous work. In the first approach, we show that the excess entropy and the transfer entropy rate of an ergodic SSP can be obtained as the limits of permutation analogues of them for the N-th order approximation by hidden Markov models, respectively. In the second approach, we employ the modified permutation partition of the set of words which considers equalities of symbols in addition to permutations of words. We show that the excess entropy and the transfer entropy rate of an ergodic SSP are equal to their modified permutation analogues, respectively.     

Simulation of polarization-sensitive optical coherence tomography images by a Monte Carlo method

We introduce a new Monte Carlo (MC) method for simulating optical coherence tomography (OCT) images of complex multilayered turbid scattering media. We demonstrate, for the first time of our knowledge, the use of a MC technique to imitate two-dimensional polarization-sensitive OCT images with nonplanar boundaries of layers in the medium like a human skin. The simulation of polarized low-coherent optical radiation is based on the vector approach generalized from the iterative procedure of the solution of Bethe-Saltpeter equation. The performances of the developed method are demonstrated both for conventional and polarization-sensitive OCT modalities.     

浅海声速剖面与移动声源的跟踪定位

在水平非均匀分布的浅海环境中,针对移动声源跟踪时,声速剖面的变化会对声场产生影响,提出了一种利用集合卡尔曼滤波算法的声速剖面跟踪反演和移动声源跟踪定位的方法。首先,将声速剖面进行距离和深度的参数化表示,从而将对声速剖面的跟踪转化为对声速剖面前3阶经验正交函数系数的跟踪;其次,通过将声源状态信息和声速剖面信息表示为状态变量,而将垂直线列阵接收到的声场信息作为测量值建立状态-测量模型,然后利用集合卡尔曼滤波方法对模型状态变量进行跟踪。仿真结果得出:声速剖面跟踪反演的均方根误差和移动声源跟踪定位的绝对误差都非常小,对声源的跟踪定位精度很高。并且通过增加集合样本数、增加接收信号信噪比以及增加接收阵元数目都可以提高跟踪定位结果精度。最后,利用东海实验数据对本方法进行了验证。      

Tracking-positioning of sound speed profiles and moving acoustic source in shallow water

An ensemble Kalman filter(EnKF) approach is proposed to perform sequential tracking of water column sound speed profile(SSP) using a moving acoustic source. First,the SSPs are discretized in depth and range, and are expressed by the empirical orthogonal functions(EOFs). Second, the acoustic source state information and the first three orders of EOF coefficients are expressed as the state variable, and the acoustic field information received by the vertical line array are the measured values. Successively, the state variables and measured values are used to establish the state-measure model. Last, the EnKF is utilized to track the state variables. The simulation results show that the root mean square error of SSP and the absolute error of source are all small, and thus the acoustic source tracking-positioning has high accuracy. Moreover, increasing the number of sample collection, the signal-to-noise ratio and the number of receiving elements can improve the tracking-positioning results. The method is verified using the experimental data of the East China Sea.     

A concise and efficient scattering matrix formalism for stable analysis of elastic wave propagation in multilayered anisotropic solids

This paper presents a concise and efficient scattering matrix formalism for stable analysis of elastic wave propagation in multilayered anisotropic solids. The formalism is capable of resolving completely the numerical instability problems associated with transfer matrix method, thereby obviating the extensive reformulation in its modified versions based on delta operator technique. In contrast to the earlier reflection matrix formalisms, all scattering matrices are obtained in a direct manner without invoking wave-propagator or scatterer operator concepts. Both local and global reflection and transmission matrices corresponding to scatterings in two and more layers are derived. The derivation of global scattering matrices in terms of the local ones is carried out concisely based on physical arguments to provide better insights into scattering mechanism. Another formulation which is even more succinct is also devised for obtaining the global scattering matrices directly from eigensolutions. The resultant expressions and algorithm are terse, efficient and convenient for implementation.     

Making metals transparent for white light by spoof surface plasmons

From first-principles computations we reveal that metallic gratings consisting of narrow slits may become transparent for extremely broad bandwidths under oblique incidence. This phenomenon can be explained by a concrete picture in which the incident wave drives free electrons on the conducting surfaces and part of the slit walls to form spoof surface plasmons (SSPs). The SSPs then propagate on the slit walls but are abruptly discontinued by the bottom edges to form oscillating charges that emit the transmitted wave. This picture explicitly demonstrates the conversion between light and SSPs and indicates clear guidelines for enhancing SSP excitation and propagation. Making structured metals transparent may lead to a variety of applications.     

On an integral equation arising in the transport of radiation through a slab involving internal reflection

In an earlier contribution to this journal [M.M.R. Williams, Eur. Phys. J. B 47, 291 (2005)], we derived an integral equation for the transmission of radiation through a slab of finite thickness which incorporated internal reflection at the surfaces. Here we generalise the problem to the case when there is a source on each face and the reflection coefficients are different at each face. We also discuss numerical and analytic solutions of the equation discussed in [M.M.R. Williams, Eur. Phys. J. B 47, 291 (2005)] when the reflection is governed by the Fresnel conditions. We obtain numerical and graphical results for the reflection and transmission coefficients, the scalar intensity and current and the emergent angular distributions at each face. The incident source is either a mono-directional beam or a smoothly varying distribution which goes from isotropic to a normal beam. Of particular interest is the philosophy of the numerical solution and whether a direct numerical approach is more effective than one involving a more elegant analytical solution using replication and the Hilbert problem. We also develop the solution of this problem using diffusion theory and compare the results with the exact transport solution. An erratum to this article is available at .     

Optimized matrix inversion technique for the T-matrix method

We suggest a new approach to calculate the inverse matrix in scattering calculations using the T-matrix method. Instead of inversion of the full matrix, we suggest the inversion of two matrices, each of which contains half the number of rows. This approach allows significant time savings and a noticeable increase of the precision of scattering calculations due to fewer arithmetical operations. An iterative method can be applied to matrices whose dimension is also divisible by factors of 2, which can further increase the time savings and accuracy.     

Modeling of wave dispersion along cylindrical structures using the spectral method

Algorithm and code are presented that solve dispersion equations for cylindrically layered media consisting of an arbitrary number of elastic and fluid layers. The algorithm is based on the spectral method which discretizes the underlying wave equations with the help of spectral differentiation matrices and solves the corresponding equations as a generalized eigenvalue problem. For a given frequency the eigenvalues correspond to the wave numbers of different modes. The advantage of this technique is that it is easy to implement, especially for cases where traditional root-finding methods are strongly limited or hard to realize, i.e., for attenuative, anisotropic, and poroelastic media. The application of the new approach is illustrated using models of an elastic cylinder and a fluid-filled tube. The dispersion curves so produced are in good agreement with analytical results, which confirms the accuracy of the method. Particle displacement profiles of the fundamental mode in a free solid cylinder are computed for a range of frequencies.     

Nonlinear mean field Fokker-Planck equations. Application to the chemotaxis of biological populations

We study a general class of nonlinear mean field Fokker-Planck equations in relation with an effective generalized thermodynamical (E.G.T.) formalism. We show that these equations describe several physical systems such as: chemotaxis of bacterial populations, Bose-Einstein condensation in the canonical ensemble, porous media, generalized Cahn-Hilliard equations, Kuramoto model, BMF model, Burgers equation, Smoluchowski-Poisson system for self-gravitating Brownian particles, Debye-Hückel theory of electrolytes, two-dimensional turbulence... In particular, we show that nonlinear mean field Fokker-Planck equations can provide generalized Keller-Segel models for the chemotaxis of biological populations. As an example, we introduce a new model of chemotaxis incorporating both effects of anomalous diffusion and exclusion principle (volume filling). Therefore, the notion of generalized thermodynamics can have applications for concrete physical systems. We also consider nonlinear mean field Fokker-Planck equations in phase space and show the passage from the generalized Kramers equation to the generalized Smoluchowski equation in a strong friction limit. Our formalism is simple and illustrated by several explicit examples corresponding to Boltzmann, Tsallis, Fermi-Dirac and Bose-Einstein entropies among others.     

Reflection and transmission of light for a layer of a naturally gyrotropic medium found in an external magnetic field

The propagation of electromagnetic waves through a layer of a naturally gyrotropic medium found in a magnetic field is studied. The problem is solved mathematically using 4×4 Berreman matrices. The Berreman Δ matrix and exact Jones reflection (R?) and transmission (T?) matrices are constructed. The reflectivity, the transmittance, the azimuth, and the ellipticity of polarization of the transmitted wave are calculated. Specific features of the irreversibility of waves are studied. It is shown that multiple reflections in a finite layer lead to a manifold increase in irreversibility, which makes possible the use of such systems as optical diodes working in reflection or as unidirectional reflectors. It is shown that three eigenmodes are excited in such media under certain conditions. The relation between the azimuths of transmitted and incident waves is found and specific features of the relation are analyzed.     

Generalized eigenproblem of hybrid matrix for Floquet wave propagation in one-dimensional phononic crystals with solids and fluids

A method based on the solution to a generalized eigenproblem of hybrid matrix is presented for stable analysis of Floquet wave propagation in one-dimensional phononic crystals with solids and fluids. The method overcomes the numerical instability in the standard eigenproblem of transfer matrix, thus enabling Floquet waves to be determined reliably. The recursion relations of hybrid matrix for periodic multilayered structure of various solid and/or fluid phases are formulated. Dispersion relation and omnidirectional reflection for one-dimensional phononic crystals with solids and fluids are discussed. The frequency-thickness range of phononic bandgap is determined conveniently based on the Floquet wavenumbers.     

Electrodynamic properties of finite-length layered media in a regular waveguide

Electrodynamic properties of finite-length layered media consisting of a sequence of dielectric plates separated from each other by air interlayers are studied. The transmission and reflection coefficients for a wave incident at such a medium are calculated by the technique of the response matrices. It is shown that a set of plates in a waveguide has alternate bands of transparency and opacity. Phenomenon of trapping of a wave excited in a set of plates in a certain frequency range is investigated.