全波场地震波数值模拟中的边界处理实践研究

在用有限差分法或有限元法模拟无界区域中的波动时,需要对计算区域的边界做特殊处理,以消除由于把地震波的传播设定在有限区域而产生的边界反射。为了这一目的,人们研究出了多种人工边界处理方法,完全匹配层（PML）吸收边界条件就是理想的方法之一,现已被广泛应用。本文将PML吸收边界条件应用于全波场地震波的数值模拟,数值计算实验表明,对qP波,匹配层的厚度为5个网格间距即可达到要求,而对qSV波与qSH波,为达到理想的吸收效果,匹配层的厚度应当增大,当厚度为13个网格间距时达到了理想的吸收效果。     

Perfectly matched layers for time-harmonic elastodynamics of unbounded domains: theory and finite-element implementation

One approach to the numerical solution of a wave equation on an unbounded domain uses a bounded domain surrounded by an absorbing boundary or layer that absorbs waves propagating outwards from the bounded domain. A perfectly matched layer (PML) is an unphysical absorbing layer model for linear wave equations that absorbs, almost perfectly, outgoing waves of all non-tangential angles-of-incidence and of all non-zero frequencies. This paper develops the PML concept for time-harmonic elastodynamics in Cartesian coordinates, utilising insights obtained with electromagnetics PMLs, and presents a novel displacement-based, symmetric finite-element implementation of the PML for time-harmonic plane-strain or three-dimensional motion. The PML concept is illustrated through the example of a one-dimensional rod on elastic foundation and through the anti-plane motion of a two-dimensional continuum. The concept is explored in detail through analytical and numerical results from a PML model of the semi-infinite rod on elastic foundation, and through numerical results for the anti-plane motion of a semi-infinite layer on a rigid base. Numerical results are presented for the classical soil–structure interaction problems of a rigid strip-footing on a (i) half-plane, (ii) layer on a half-plane, and (iii) layer on a rigid base. The analytical and numerical results obtained for these canonical problems demonstrate the high accuracy achievable by PML models even with small bounded domains.     

Long-Time Simulations of Nonlinear Schrödinger-Type Equations using Step Size Exceeding Threshold of Numerical Instability

We propose an exponential time-differencing method based on the leapfrog scheme for numerical integration of the generalized nonlinear Schrödinger-type equations. The key advantage of the proposed method over the widely used Fourier split-step method is that in the new method, numerical instability at high wavenumbers is strongly suppressed. This allows one to use time steps that considerably exceed the instability threshold, which leads to a proportional reduction of the computational time. Moreover, we introduce a technique that eliminates numerical instability at low-to-moderate wavenumbers that is common for methods based on the leapfrog scheme. We illustrate the performance of the proposed method with examples from two applications areas: deterministic wave turbulence and solitary waves.     

Development of PML absorbing boundary conditions for computational aeroacoustics: A progress review

Recent advances in the development of perfectly matched layer (PML) as absorbing boundary conditions for computational aeroacoustics are reviewed. The PML methodology is presented as a complex change of variables. In this context, the importance of a proper space-time transformation in the PML technique for Euler equations is emphasized. A unified approach for the derivation of PML equations is offered that involves three essential steps. The three-step approach is illustrated in details for the PML of linear and non-linear Euler equations. Numerical examples are also given that include non-reflecting boundary conditions for a ducted channel flow and mixing layer roll-up vortices.     

Solution of time-domain Maxwell equation with PML by using modified Laguerre polynomials

In this work, a stable numerical algorithm proposed by Chung et al. for the time-domain Maxwell equations is generalized. The time-domain Maxwell equations are solved by expressing the transient behaviors in terms of the modified Laguerre polynomials, and then the original equations of the initial value and boundary value can be transformed into a series of problems independent of the time variable. In this case the method of finite difference （FD）, the finite element method （FEM）, the method of moment （MoM）, etc. or the combination of these methods can be used to solve the problems. Finally, a numerical model is provided for the scattering problem with perfect matched layer （PML） by using FD. The comparison between the results of the proposed method and FDTD is presented to verify the proposed new method.     

On the existence and convergence of the solution of PML equations

In this article we study the mesh termination method in computational scattering theory known as the method of Perfectly Matched Layer (PML). This method is based on the idea of surrounding the scatterer and its immediate vicinity with a fictitious absorbing non-reflecting layer to damp the echoes coming from the mesh termination surface. The method can be formulated equivalently as a complex stretching of the exterior domain. The article is devoted to the existence and convergence questions of the solutions of the resulting equations. We show that with a special choice of the fictitious absorbing coefficient, the PML equations are solvable for all wave numbers, and as the PML layer is made thicker, the PML solution converge exponentially towards the actual scattering solution. The proofs are based on boundary integral methods and a new type of near-field version of the radiation condition, called here the double surface radiation condition. Partly supported by the Finnish Academy, project 37692.     

Application of the image-analysis method to studies of the space-time wave evolution in an annular gas-liquid flow

A secondary instability of the surface of a liquid film sheared by an intense gas flow was found for the first time by our team using the high-speed modified laser-induced fluorescence method. An image-analysis toolkit (in particular, the Canny method) was used for studying space-time evolution of various types of waves. As a result, quantitative characteristics of both the types of waves and the characteristics of the instability phenomena of back slopes of primary waves were obtained.     

Perfectly Matched Layers for Time-Harmonic Second Order Elliptic Problems

The main goal of this work is to give a review of the Perfectly Matched Layer (PML) technique for time-harmonic problems. Precisely, we focus our attention on problems stated in unbounded domains, which involve second order elliptic equations writing in divergence form and, in particular, on the Helmholtz equation at low frequency regime. Firstly, the PML technique is introduced by means of a simple porous model in one dimension. It is emphasized that an adequate choice of the so called complex absorbing function in the PML yields to accurate numerical results. Then, in the two-dimensional case, the PML governing equation is described for second order partial differential equations by using a smooth complex change of variables. Its mathematical analysis and some particular examples are also included. Numerical drawbacks and optimal choice of the PML absorbing function are studied in detail. In fact, theoretical and numerical analysis show the advantages of using non-integrable absorbing functions. Finally, we present some relevant real life numerical simulations where the PML technique is widely and successfully used although they are not covered by the standard theoretical framework.     

基于正交轨迹的共形PML数值仿真方法

从材料的观点看,完全匹配层是各向异性的介质.场在完全匹配层遵从麦克斯韦方程.通过设计磁导率和介电常数张量参数来实现任意形状的共形完全匹配层.为匹配层动态稳定,可应用二维正交轨迹网格和FDTD方法实现了数值共形完全匹配层仿真,数值仿真结果表明基于正交轨迹网格的共形完全匹配层为一个共形FDTD或FEM方法提供了一个有效地吸收边界条件.     

带矢量ABC底面的共形完全匹配层

共形完全匹配层是一种有耗各向异性媒质组成的凸且光滑的壳体,其底面一般是PEC面或PMC面,但是PEC面或PMC面会对原散射场产生反射;为了减少底面反射,将CPML原有的PEC（或PMC）底面改为矢量ABC吸收边界,并给出了带矢量ABC底面的CPML泛函公式。通过数值算例证明,这种带矢量ABC底面的CPML边界不仅减少了底面反射,而且吸收效果好,计算精度高。     

High Order Weighted Essentially Non-oscillation Schemes for Two-Dimensional Detonation Wave Simulations

In this paper, we demonstrate the detailed numerical studies of three classical two dimensional detonation waves by solving the two dimensional reactive Euler equations with species with the fifth order WENO-Z finite difference scheme (Borges et al. in J. Comput. Phys. 227:3101?C3211, 2008) with various grid resolutions. To reduce the computational cost and to avoid wave reflection from the artificial computational boundary of a truncated physical domain, we derive an efficient and easily implemented one dimensional Perfectly Matched Layer (PML) absorbing boundary condition (ABC) for the two dimensional unsteady reactive Euler equation when one of the directions of domain is periodical and inflow/outflow in the other direction. The numerical comparison among characteristic, free stream, extrapolation and PML boundary conditions are conducted for the detonation wave simulations. The accuracy and efficiency of four mentioned boundary conditions are verified against the reference solutions which are obtained from using a large computational domain. Numerical schemes for solving the system of hyperbolic conversation laws with a single-mode sinusoidal perturbed ZND analytical solution as initial conditions are presented. Regular rectangular combustion cell, pockets of unburned gas and bubbles and spikes are generated and resolved in the simulations. It is shown that large amplitude of perturbation wave generates more fine scale structures within the detonation waves and the number of cell structures depends on the wave number of sinusoidal perturbation.     

基于鲁棒自适应Kalman滤波的PET放射性浓度重建

针对正电子发射断层成像重建过程中存在的系统模型误差和投影数据不确定性,提出了基于状态空间体系的鲁棒自适应Kalman滤波法。该方法根据药物动力学先验信息建立状态方程,结合PET测量方程组成状态空间模型。引入虚拟噪声来表示模型的系统矩阵误差之后,通过应用鲁棒自适应Kalman滤波法对未知的系统噪声以及观测噪声进行估计的同时完成PET放射性浓度的重建。实验结果表明,此算法比传统的最大似然法和滤波反投影法更具鲁棒性,适合应用于实际PET系统中。     

Mixed perfectly-matched-layers for direct transient analysis in 2D elastic heterogeneous media

We are concerned with elastic waves arising in plane-strain problems in an elastic semi-infinite arbitrarily heterogeneous medium. Specifically, we discuss the development of a new mixed displacement–stress formulation for forward elastic wave simulations in perfectly-matched-layer (PML)-truncated heterogeneous media.To date, most PML formulations split the displacement and stress fields, resulting in non-physical components for each field. In this work, we favor unsplit schemes, primarily for the relative ease by which the resulting forms can be incorporated into existing codes, the ease by which the resulting semi-discrete forms can be integrated in time, and the ease by which they can be used in adjoint formulations arising in inverse problems, contrary to most past and current developments. We start by following classical lines, and apply complex-coordinate-stretching to the governing equations in the frequency domain, while retaining both displacements and stress quantities as unknowns. With the aid of auxiliary variables the resulting mixed form is rendered second-order in time, thereby allowing the use of standard time integration schemes. We report on numerical simulations demonstrating the stability and efficacy of the approach.     

Convergence of linearized backward Euler–Galerkin finite element methods for the time-dependent Ginzburg–Landau equations with temporal gauge

We present error analysis of fully discrete Galerkin finite element methods for the time-dependent Ginzburg–Landau equations with the temporal gauge, where a linearized backward Euler scheme is used for the time discretization. We prove that the convergence rate is O(τ+h^r) if the finite element space of piecewise polynomials of degree r is used. Due to the degeneracy of the problem, the convergence rate is one order lower than the optimal convergence rate of finite element methods for parabolic equations. Numerical examples are provided to support our theoretical analysis.     

Influence of centrifugal forces on the development of wave packets in boundary layers with a uniformly valid model

The behavior of three-dimensional wave packets in the boundary layer on curved surfaces is analyzed in this study based on a modification of the triple-deck theory referred to as the “criss-cross” interaction model. The equations of the criss-cross interaction describe a particular type of boundary layer instability mode that possesses underlying properties of both the Tollmien-Schlichting waves and Taylor-Görtler vortices. Previous analysis of the criss-cross interaction regime suggests a possibility for upstream propagation of perturbations in the boundary layer and possible absolute instability of the flow. However, these results cannot be considered as conclusive because the initial-value problem for the criss-cross interaction equations is ill-posed. In a recent work [Turkyilmazoglu M, Ruban AI. A uniformly valid well-posed asymptotic approach to the inviscid-viscous interaction theory. Stud Appl Math 2002;108:161-85] a regularized non-asymptotic model to describe criss-cross interaction has been proposed. Whereas in the original version of the theory, perturbations have an unbounded growth rate as the longitudinal wave number ∣k∣ → ∞, in the new model of [Turkyilmazoglu M, Ruban AI. A uniformly valid well-posed asymptotic approach to the inviscid-viscous interaction theory. Stud Appl Math 2002;108:161-85], as physically expected the amplification rate remains bounded for both spatially growing and temporally growing waves. A Fourier transform method is used in the present study to solve the linearized equations for the flow over concave roughness and humps and it is found that disturbances develop and are convected downstream as wave packets. The behavior of the wave packets is consistent with convective instability, and the upstream influence is no longer present at large times.     

Long-time behavior of PML absorbing boundaries for layered periodic structures

In this work we consider a special case of the Perfectly Matched Layer (PML) divergence which is observed by the simulation of the planar periodic structures such as photonic crystal slabs or antenna arrays. This divergence is caused by an excitation of long-living artefact evanescent waves in these structures by an incident external pulse. We study the application of the known remedies to this problem: increasing the distance between the structure and PML, employing the κ parameter, employing non-PML absorbers. We also suggest a new simple and effective solution, where the usual PML is backed by an additional absorbing layer.     

Unconditionally stable higher order perfectly matched layer applied to terminate anisotropic magnetized plasma

An effective and efficient absorbing boundary condition applied to terminate anisotropic magnetized plasma is proposed in this article. As its mathematical derivations are based on the high order perfectly matched layer (HO‐PML) and Crank‐Nicolson finite‐difference time‐domain method, it has the advantages of reducing late‐time reflections, attenuating evanescent waves, absorbing low‐frequency propagation waves and overcoming the Courant‐Friedrichs‐Levy condition. To validate the efficiency and effectiveness of the proposed HO‐PML, a ridge waveguide model composed of vacuum and the anisotropic magnetized plasma, a dielectric waveguide model composed of the dielectric and the anisotropic magnetized plasma are employed. The numerical results show that the proposal can not only maintain the stability of the algorithm with the increment of the time step but also further enhance the absorbing performance.     

Spatial direct numerical simulation of transitional boundary layer over compliant surfaces

An iterative implicit fractional step method is developed and employed for the simulation of transitional boundary layer over compliant surfaces. The three-dimensional perturbation Navier-Stokes equations are discretized by curvilinear finite volumes on a collocated grid system. A multigrid procedure is used for computations associated with the pressure-Poisson equation, while simulation is carried out by MPI-based parallel computation with domain decomposition. Results in the literature for oblique linear waves and the non-linear breakdown of a wave triad over a rigid wall are repeated to check the accuracy of the codes that had been developed. Oblique linear TS (Tollmien-Schlichting) waves over finite-length compliant membranes are generally found to coexist with CIFI (compliance induced flow instability) or FISI (flow induced surface instabilities) waves. The latter waves usually possess longer wavelengths and thus propagate at a larger oblique wave angles than the TS waves. Simulation reveals that compliant surfaces may slow down the development of secondary instabilities during the early stages of laminar-turbulent transition. However, during the later stages of fundamental wave breakdown, interactions with CIFI and edge-generated waves may increase the amplitudes of the original 2D and 3D TS waves, leading to an earlier breakdown on compliant surfaces. Linear interaction between the flow and compliant membranes has been assumed.     

A column generation approach for the maximal covering location problem

This article presents a column generation algorithm to calculate new improved lower bounds to the solution of maximal covering location problems formulated as a p‐median problem. This reformulation leads to instances that are difficult for column generation methods. The traditional column generation method is compared with the new approach, where the reduced cost criterion used at the column selection is modified by a Lagrangean/surrogate multiplier. The efficiency of the new approach is tested with real data, where computational tests were conducted and showed the impact of sparsity and degeneracy on column generation‐based methods.     

Numerical modeling of ground-penetrating radar in 2-D using MATLAB

We present MATLAB codes for finite-difference time-domain (FDTD) modeling of ground-penetrating radar (GPR) in two dimensions. Surface-based reflection GPR is modeled using a transverse magnetic (TM-) mode formulation. Crosshole and vertical radar profiling (VRP) geometries are modeled using a transverse electric (TE-) mode formulation. Matrix notation is used in the codes wherever possible to optimize them for speed in the MATLAB environment. To absorb waves at the edges of the modeling grid, we implement perfectly matched layer (PML) absorbing boundaries. Although our codes are two-dimensional and do not incorporate features such as dispersion in electrical properties, they capture many of the important elements of GPR surveying and run at a fraction of the computational cost of more elaborate algorithms. In addition, the codes are well commented, relatively easy to understand, and can be easily modified for the user's specific purpose.