Functionally graded materials for impedance matching in elastic media

When functionally graded material layers are inserted between two impedance mismatching media, passbands with extremely large bandwidths can appear in these layered systems. An accurate and effective iterative method is developed to deal with these layered systems with extremely large layer number.     

The propagation of ultrasound in porous media

The equations of Ying and Truell, and Waterman and Truell, describing the propagation of ultrasound in two-phase materials are solved numerically for porous solids, and are found to give unphysical results for high porosity. A new self-consistent theory, which can be solved analytically, is presented and is shown to have reasonable behaviour at high porosity.     

Guided wave propagation in multilayered piezoelectric structures

A general formulation of the method of the reverberation-ray matrix (MRRM) based on the state space formalism and plane wave expansion technique is presented for the analysis of guided waves in multilayered piezoelectric structures. Each layer of the structure is made of an arbitrarily anisotropic piezoelectric material. Since the state equation of each layer is derived from the three-dimensional theory of linear piezoelectricity, all wave modes are included in the formulation. Within the framework of the MRRM, the phase relation is properly established by excluding exponentially growing functions, while the scattering relation is also appropriately set up by avoiding matrix inversion operation. Consequently, the present MRRM is unconditionally numerically stable and free from computational limitations to the total number of layers, the thickness of individual layers, and the frequency range. Numerical examples are given to illustrate the good performance of the proposed formulation for the analysis of the dispersion characteristic of waves in layered piezoelectric structures. Supported by the National Natural Science Foundation of China (Grant Nos. 10725210 and 10832009), the Specialized Research Fund for the Doctoral Program of Higher Education (Grant No. 20060335107), the National Basic Research Program of China (Grant No. 2009CB623204), and the Scientific Research Foundation for Tsuiying Talents of Lanzhou University     

一种处理色散介质问题的通用时域有限差分方法

色散介质的介电系数是频率的函数,使本构关系在时域成为卷积关系.这就给用时域有限差分方法计算色散介质中波的散射和传播带来了困难.现有算法往往要针对不同色散介质模型推导相应的递推公式,算法的通用性较差.本文完善和发展了移位算子-时域有限差分方法,使之成为一种处理色散介质电磁问题的通用方法.首先,证明了常见的三种色散介质模型(德拜模型、洛伦兹模型和德鲁模型)的介电系数均可以写成适于移位算子法计算的有理分式函数形式.然后,用/t代替jω,过渡到时域,再引入时域移位算子z_t代替时间微分算子来处理有理分式函数形式的介电系数,给出离散时域本构关系的表示式,进而导出时域有限差分方法当中电位移矢量和电场强度之间的关系.最后,计算了几种色散介质的电磁散射,数值结果表明了本文方法和程序的通用性和正确有效性. 关键词： 时域有限差分方法 色散介质 移位算子     

Transition in the propagation mechanism during flame acceleration in porous media

The combustion of stoichiometric hydrogen-air at various initial pressures was investigated in a 7.62 cm square cross-section channel filled with 1.27 cm diameter beads. The flame time-of-arrival and pressure time history along the channel were obtained by ionization probes and piezoelectric pressure transducers. Flame acceleration was found to be very rapid, e.g. at an initial pressure of 45 kPa the flame achieves a velocity of over 600 m/s in roughly 0.3 m. It was determined that at this high speed a well defined planar shock wave precedes a thick reaction zone. It was also shown that there is a transition in the flame propagation mechanism, similar to that observed in an obstacle laden channel [G. Ciccarelli and C. Johansen, The role of shock-flame interactions on flame acceleration in an obstacle laden channel, Proc. 22nd International Colloquium on the Dynamics of Explosions and Reactive Systems, Minsk, 2009]. By varying the initial pressure of the mixture, changes in the axial location of the transition between combustion propagation regimes was also observed. A soot foil technique was used to identify the transition in the propagation mechanism, as well as to provide information concerning the local flow field around the beads and the overall average flow direction.     

Multiscale finite-volume method for parabolic problems arising from compressible multiphase flow in porous media

The multiscale finite-volume (MSFV) method was originally developed for the solution of heterogeneous elliptic problems with reduced computational cost. Recently, some extensions of this method for parabolic problems have been proposed. These extensions proved effective for many cases, however, they are neither general nor completely satisfactory. For instance, they are not suitable for correctly capturing the transient behavior described by the parabolic pressure equation. In this paper, we present a general multiscale finite-volume method for parabolic problems arising, for example, from compressible multiphase flow in porous media. Opposed to previous methods, here, the basis and correction functions are solutions of full parabolic governing equations in localized domains. At the same time, to enhance the computational efficiency of the scheme, the basis functions are kept pressure independent and do not have to be recalculated as pressure evolves. This general approach requires no additional assumptions and its good efficiency and high accuracy is demonstrated for various challenging test cases. Finally, to improve the quality of the results and also to extend the scheme for highly anisotropic heterogeneous problems, it is combined with the iterative MSFV (i-MSFV) method for parabolic problems. As one iterates, the i-MSFV solutions of compressible multiphase problems (parabolic problems) converge to the corresponding fine-scale reference solutions in the same way as demonstrated recently for incompressible cases (elliptic problems). Therefore, the proposed MSFV method can also be regarded as an efficient linear solver for parabolic problems and studies of its efficiency are presented for many test cases.     

Stochastic model in microwave propagation

Further experimental results of delay time in microwave propagation are reported in the presence of a lossy medium (wood). The measurements show that the presence of a lossy medium makes the propagation slightly superluminal. The results are interpreted on the basis of a stochastic (or path integral) model, showing how this model is able to describe each kind of physical system in which multi-path trajectories are present.     

非局域克尔介质中厄米高斯光束传输的变分研究

在非局域非线性克尔介质中,通过对介质实对称响应函数的泰勒展开,简化了非局域非线性薛定谔方程所对应的Lagrange密度,进而利用变分法对光束的传输问题进行了分析.求出试探解各个参量的演化方程并得到了自聚焦介质中的厄米高斯型光束的精确解析解,当输入功率达到临界功率时,即形成高阶空间光孤子(厄米高斯孤子),其最低阶(基模光孤子)就是高斯孤子.通过数值模拟发现解析解与数值解符合得很好. 关键词： 非局域克尔介质 变分法 厄米高斯光束 空间光孤子     

Nonlinear pulse propagation in periodic media with and without defects

The influence of single defects within a periodic structure on the nonlinear transmission of optical pulses through the structure is investigated numerically. A stack of alternating linear and cubically nonlinear layers of submicron thicknesses is considered. The simulations are based on a modified equation for the pulse envelope. Diffraction of the pulse has been taken into account, too. The results of simulations have been demonstrated as an essential influence of defects within the periodic structure on the nonlinear propagation of optical pulses, when the carrier frequency is chosen within the stop band of the structure with the defect.     

Radiation propagation in random media: From positive to negative correlations in high-frequency fluctuations

We survey research on radiation propagation or ballistic particle motion through media with randomly variable material density, and we investigate the topic with an emphasis on very high spatial frequencies. Our new results are based on a specific variability model consisting of a zero-mean Gaussian scaling noise riding on a constant value that is large enough with respect to the amplitude of the noise to yield overwhelmingly non-negative density. We first generalize known results about sub-exponential transmission from regular functions, which are almost everywhere continuous, to merely “measurable” ones, which are almost everywhere discontinuous (akin to statistically stationary noises), with positively correlated fluctuations. We then use the generalized measure-theoretic formulation to address negatively correlated stochastic media without leaving the framework of conventional (continuum-limit) transport theory. We thus resolve a controversy about recent claims that only discrete-point process approaches can accommodate negative correlations, i.e., anti-clustering of the material particles. We obtain in this case the predicted super-exponential behavior, but it is rather weak. Physically, and much like the alternative discrete-point process approach, the new model applies most naturally to scales commensurate with the inter-particle distance in the material, i.e., when the notion of particle density breaks down due to Poissonian—or maybe not-so-Poissonian—number-count fluctuations occur in the sample volume. At the same time, the noisy structure must prevail up to scales commensurate with the mean-free-path to be of practical significance. Possible applications are discussed.     

Time domain numerical modeling of wave propagation in 2D heterogeneous porous media

This paper deals with the numerical modeling of wave propagation in porous media described by Biot’s theory. The viscous efforts between the fluid and the elastic skeleton are assumed to be a linear function of the relative velocity, which is valid in the low-frequency range. The coexistence of propagating fast compressional wave and shear wave, and of a diffusive slow compressional wave, makes numerical modeling tricky. To avoid restrictions on the time step, the Biot’s system is splitted into two parts: the propagative part is discretized by a fourth-order ADER scheme, while the diffusive part is solved analytically. Near the material interfaces, a space–time mesh refinement is implemented to capture the small spatial scales related to the slow compressional wave. The jump conditions along the interfaces are discretized by an immersed interface method. Numerical experiments and comparisons with exact solutions confirm the accuracy of the numerical modeling. The efficiency of the approach is illustrated by simulations of multiple scattering.     

A method to identify acoustic reverberation in multilayered homogeneous media

The presence of reverberation is a source of artifacts that can hinder the analysis of ultrasound signals and images. Besides compromising image generation, these artifacts can introduce errors in the quantitative parameter estimation in fields such as material and biological tissue characterization. A method that allows the separation between the first reflection on an interface and all the other reflections from the same interface (reverberation) could improve the quality of these images as well as the precision and accuracy in quantitative results. This work presents an algorithm for the identification of reverberating echoes in multilayered media, based on the comparison of their power spectra (estimated via FFT), through a least mean square approach, and on the temporal relationship among them. It considers that the global effect from attenuation, reflection and transmission coefficients for each layer causes spectral differences that could differentiate echoes that pass through one layer or another. The results of 10 simulations and of 60 experiments, carried out with 6 different phantoms (10 experiments with each one) are presented and discussed. It was found that the algorithm provided a correct identification for 85% of the simulated and 86.6% of the experimental echoes collected from the 60 experiments.     

一种简化非Gauss场模型在随机介质中光传播问题中的应用

为了研究大气湍流间歇性的光传播效应，构造出一种比较简单的非Gauss场模型(Poission场 )用于描述大气介电常数(或折射率)随机起伏.模型特征泛函含有四个待定函数，根据大气湍 流的统计均匀性，介电起伏的单点概率分布函数，以及介电起伏能谱可以选择或确定它们. 对在这种简化湍流中传播的光波平均场及二阶统计矩性质进行了理论分析，并给出数值模拟 的一个简单例子. 关键词： 光波传播 大气湍流 间歇性     

Development of a finite element model for coupled radiative and conductive heat transfer in participating media

Considering the geometrical applicability, a finite element model (FEM) for coupled radiative-conductive heat transfer has been developed which is applicable to enclosures of arbitrary geometry in present research. The present work provides a solution of coupled heat transfer in a rectangular, cylindrical or annulus enclosure with black or gray walls containing an absorbing-emitting-scattering medium. It is also applied to study the influence of conductive/radiation coefficient, albedo and wall emissivity on the temperature distribution in the medium. Compared with the results available in other references, the present FEM has no limitation with respect to geometry and can predict the coupled radiative-conductive heat transfer in participating media accurately.     

运动物体在颗粒物质中的动力学过程及最大穿透深度仿真研究

利用离散元法仿真了运动物体在颗粒物质中的三维动力学过程, 仿真采用周期边界条件, 并考虑了重力、接触力、阻尼力、摩擦力的影响. 将仿真结果和相关的三维实验结果进行了对比, 两者符合较好. 仿真结果表明穿透深度与运动物体的冲击速度、运动物体质量、颗粒介质床的密度均有关系. 运动物体质量越大, 速度越快, 则穿透越深, 而且穿透深度和质量呈线性关系. 仿真过程较为真实地再现了小颗粒的飞溅现象. 关键词： 颗粒物质 动力学过程 仿真 离散元法     

Dynamics of a stochastic rumor propagation model incorporating media coverage and driven by Lévy noise

With the development of information technology, rumors propagate faster and more widely than in the past. In this paper, a stochastic rumor propagation model incorporating media coverage and driven by Lévy noise is proposed. The global positivity of the solution process is proved, and further the basic reproductive number R₀ is obtained. When R₀ < 1, the dynamical process of system with Lévy jump tends to the rumor-free equilibrium point of the deterministic system, and the rumor tends to extinction; when R₀ > 1, the rumor will keep spreading and the system will oscillate randomly near the rumor equilibrium point of the deterministic system. The results show that the oscillation amplitude is related to the disturbance of the system. In addition, increasing media coverage can effectively reduce the final spread of rumors. Finally, the above results are verified by numerical simulation.     

Dynamics and near-optimal control in a stochastic rumor propagation model incorporating media coverage and Lévy noise

The appearance of rumors intensifies people's panic and affects social stability. How to control the spread of rumors has become an important issue which is worth studying. In order to more accurately reflect the actual situation in the real world, a stochastic model incorporating media coverage and Lévy noise is proposed to describe the dynamic process of rumor propagation. By introducing two control strategies of popular science education and media coverage in an emergency event, an near-optimal control problem that minimizes the influence and control cost of rumor propagation is proposed. Sufficient conditions for near-optimal control of the model are established by using a Hamiltonian function. Then the necessary conditions for near-optimal control are obtained by using the Pontryagin maximum principle. Finally, the effect of popular science education, media coverage and Lévy noise on rumor propagation process control is verified by numerical simulation.     

用元胞自动机模型研究二维激发介质中的非线性波

采用Greenberg-Hasting元胞自动机模型研究了激发介质中非线性波的有关性质：在零流边界条件下根据计算机数值模拟研究了行波波速随邻域半径、激发阈值的变化关系,分析得到激发介质的激发条件;同时模拟研究了螺旋波波头轨迹受不应态数和激发阈值的影响并分析产生这些影响的内在微观机理. 关键词： 元胞自动机 激发介质 非线性波     

A hierarchical model for random walks in random media

We show that the random walk generated by a hierarchical Laplacian in ^d has standard diffusive behavior. Moreover, we show that this behavior is stable under a class of random perturbations that resemble an off-diagonal disordered lattice Laplacian. The density of states and its asymptotic behavior around zero energy are computed: singularities appear in one and two dimensions.