A fast volume integral equation method for the direct/inverse problem in elastic wave scattering phenomena

A fast method for solving the volume integral equation is introduced for the solution of forward and inverse multiple scattering problems in an elastic 3-D full space. For both forward and inverse scattering analysis, the volume integral equation in the wavenumber domain is used. By means of the discrete Fourier transform, the volume integral equation in the wavenumber domain can be dealt with as a Fredholm equation of the 2nd kind with respect to a non-Hermitian operator on a finite dimensional vector space. The Bi-CGSTAB method is employed to construct the Krylov subspace in the wavenumber domain. The current procedure establishes a fast and simplified method without requiring the derivation of a coefficient matrix. Several numerical results validate the accuracy and effectiveness of the current method for both forward and inverse scattering analysis. According to the numerical results, the reconstruction of inhomogeneities of the wave field is successful, even for multiple scattering of several cubes.     

任意区域上的粒子均匀分布方法

给出了任意区域内粒子的Voronoi面积均匀划分的方法,对给定的区域边界点,先由边界点得到边界粒子,再利用分布函数结合Delaunay加点算法得到内部粒子,这种方法可以得到一般的区域上粒子的均匀分布。这里给出了一个流体力学中利用SPH方法模拟内爆问题的例子,算例表明,均匀分布可以得到更高的无网格方法数值模拟精度。     

SH波散射与界面圆孔附近的动应力集中

建立了求解含有界面圆孔的二种不同弹性组合介质中ＳＨ波的散射和界面圆孔附近的动应力集中问题的Ｇｒｅｅｎ函数法给出了一个具有半圆形缺口的弹性半空间水平表面上任意一点承受时间谐和的出平面线源荷载作用时位移函数的基本解取基本解作为Ｇｒｅｅｎ函数，建立起问题的定解积分方程最后给出了界面圆孔的动应力集中的算例和结果，并讨论了不同介质参数的组合对动应力集中的影响     

Three-dimensional formulation of time-dependent wave propagation in a gyroanisotropic medium by Bergeron''s method

This paper discusses a numerical vector analysis method for wave propagation and scattering in three-dimensional space and the time domain. The approach is based on both the equivalent-circuit representation of Maxwell's equations and the formulation by Bergeron's method in the time domain. As an example of the application of this technique to a complicated medium, we consider time-dependent electromagnetic wave propagation in a gyroanisotropic medium, specifically the magnetized ferrite and the magnetized plasma. We demonstrate that the numerical approach accurately models relevant physical phenomena, including the Faraday rotation effect.     

Rate dependent properties of perovskite type tetragonal piezoelectric materials using micromechanical model

In this paper, a three-dimensional micromechanical model is presented for simulation of the rate dependent properties of certain perovskite type tetragonal piezoelectric materials. The model is based on linear constitutive, nonlinear domain switching, and linear kinetics theories. The simulation starts with a virgin bulk material of randomly oriented grains. Then the material is electrically loaded with an alternating voltage of various frequencies, which are in the order of 0.01 Hz to 1 Hz. An energy equation in combination with a probability function is used to determine the onset of the domain switching inside the grains. Such a probability function leads to a better phenomenological model for the domain switching even for electrical loadings, which are in a range far below the coercive fields. The propagation of the domain wall during the domain switching process in grains is modeled by means of linear kinetics relations after domain nucleation. The response of the bulk ceramic is predicted by averaging the response of individual grains using Euler angles for the transformation from local coordinates of the grains to global coordinate. Electric displacement hysteresis loops for different frequencies and amplitudes of the alternating electric fields are simulated. A simple micromechanical model without the probabilistic approach is compared with the one that takes it into account. Both models give important insights into the rate dependency of piezoelectric materials, which was observed in some experiments reported in the literature.     

Overall constitutive description of symmetric layered media based on scattering of oblique SH waves

This papers investigates the scattering of oblique shear horizontal (SH) waves off finite periodic media made of elastic and viscoelastic layers. It further considers whether a Willis-type constitutive matrix (in temporal and spatial Fourier domain) may reproduce the scattering matrix (SM) of such a system. In answering this question the procedure to determine the relevant overall constitutive parameters for such a medium is presented. To do this, first the general form of the dispersion relation and impedances for oblique SH propagation in such coupled Willis-type media are developed. The band structure and scattering of layered media are calculated using the transfer matrix (TM) method. The dispersion relation may be derived based on the eigen-solutions of an infinite periodic domain. The wave impedances associated with the exterior surfaces of a finite thickness slab are extracted from the scattering of such a system. Based on reciprocity and available symmetries of the structure and each constituent layer, the general form of the dispersion and impedances may be simplified. The overall quantities may be extracted by equating the scattering data from TM with those expected from a Willis-type medium. It becomes evident that a Willis-type coupled constitutive tensor with components that are assumed independent of wave vector is unable to reproduce all oblique scattering data. Therefore, non-unique wave vector dependent formulations are introduced, whose SM matches that of the layered media exactly. It is further shown that the dependence of the overall constitutive tensors of such systems on the wave vector is not removable even at very small frequencies and incidence angles and that analytical considerations significantly limit the potential forms of the spatially dispersive constitutive tensors.     

Flow simulation on moving boundary‐fitted grids and application to fluid–structure interaction problems

We present a method for the parallel numerical simulation of transient three‐dimensional fluid–structure interaction problems. Here, we consider the interaction of incompressible flow in the fluid domain and linear elastic deformation in the solid domain. The coupled problem is tackled by an approach based on the classical alternating Schwarz method with non‐overlapping subdomains, the subproblems are solved alternatingly and the coupling conditions are realized via the exchange of boundary conditions. The elasticity problem is solved by a standard linear finite element method. A main issue is that the flow solver has to be able to handle time‐dependent domains. To this end, we present a technique to solve the incompressible Navier–Stokes equation in three‐dimensional domains with moving boundaries. This numerical method is a generalization of a finite volume discretization using curvilinear coordinates to time‐dependent coordinate transformations. It corresponds to a discretization of the arbitrary Lagrangian–Eulerian formulation of the Navier–Stokes equations. Here the grid velocity is treated in such a way that the so‐called Geometric Conservation Law is implicitly satisfied. Altogether, our approach results in a scheme which is an extension of the well‐known MAC‐method to a staggered mesh in moving boundary‐fitted coordinates which uses grid‐dependent velocity components as the primary variables. To validate our method, we present some numerical results which show that second‐order convergence in space is obtained on moving grids. Finally, we give the results of a fully coupled fluid–structure interaction problem. It turns out that already a simple explicit coupling with one iteration of the Schwarz method, i.e. one solution of the fluid problem and one solution of the elasticity problem per time step, yields a convergent, simple, yet efficient overall method for fluid–structure interaction problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Parallel computing of radiative transfer in relativistic jets using Monte Carlo method

We present numerical attempts of radiative transfer in a relativistic scattering flow that can produce gamma rays using a three-dimensional Monte Carlo code. We prepared an initial background flowfield obtained from hydrodynamical simulation of a relativistic jet in which Thomson scattering dominates compared to absorption, and solved the radiative transfer equation for the background evolved by a simple expansion model. Since a large number of sample particles is required for an accurate computation, we have parallelized the Monte Carlo code in order to obtain solutions in a practical computational time even for a long-term simulation coupled with a time-dependent flowfield. Using this code, higher parallel efficiency is achieved with larger number of particles. The obtained light curve from the simple model shows a signal of the transition from the opaque post-shock flow to the transparent regime as the flow expands, and the high-energy photons are generated by not only the Doppler boosting but also the inverse Compton scattering.     

Numerical study of flows of complex fluids between eccentric cylinders using transformation functions

In this paper, we investigate fluid flows between eccentric cylinders by means of two stream‐tube analyses. The first method considers a one‐to‐one global transformation function that allows the physical domain to be transformed into a mapped domain, used as computational domain, that involves concentric streamlines. The second approach uses local transformations and domain decomposition techniques to deal with mixed flow regimes. Both formulations are particularly adapted for handling time‐dependent constitutive equations, since particle‐tracking problems are avoided. Mass conservation is verified in both formulations and the relevant numerical procedure can be carried out using simple meshes built on the mapped streamlines. Fluids obeying anelastic and viscoelastic constitutive equations are considered in the calculations. The numerical results are consistent with those in the literature for the flow rates tested. Application of the method to the K‐BKZ memory‐integral constitutive equation highlights significant differences between the model predictions and those provided by more simple rheological models. Copyright © 2002 John Wiley & Sons, Ltd.     

弹性静力问题的无网格弱-强形式结合法

基于改进的移动最小二乘(MLS)二阶导数近似,建立了一种求解弹性静力问题的无网格弱-强形式结合法(MLS-MWS)。该方法采用节点离散求解域,通过MLS构造形函数,将求解域划分为边界域和内部域,并分别使用控制方程的局部弱形式和强形式来建立离散系统方程。对强形式中涉及的近似函数二阶导数计算,提出了一种将其转化为求两次一阶导数的方法,与传统方法相比,该方法计算简单、精度高。MLS-MWS法结合了弱、强形式无网格法的优点,Neumann边界条件容易满足,并且只需在边界区域进行积分。文中应用该方法分析了两个弹性力学平面问题,分析结果表明本文方法具有良好的精度和收敛性。     

A MESH-FREE ALGORITHM FOR DYNAMIC IMPACT ANALYSIS OF HYPERELASTICITY

A mesh-free method based on local Petrov-Galerkin formulation is presented to solve dynamic impact problems of hyperelastic material.In the present method,a simple Heaviside test function is chosen for simplifying domain integrals.Trial function is constructed by using a radial basis function(RBF)coupled with a polynomial basis function,in which the shape function possesses the kronecker delta function property.So,additional treatment is not required for imposing essential boundary conditions.Governing equations of impact problems are established and solved node by node by using an explicit time integration algorithm in a local domain,which is very similar to that of the collocation method except that numerical integration can be implemented over local domain in the present method.Numerical results for several examples show that the present method performs well in dealing with the dynamic impact problem of hyperelastic material.     

Propagation of Transient Acoustic Waves in Layered Porous Media: Fractional Equations for the Scattering Operators

Acoustic waves scattering from a rigid air-saturated porous medium is studied in the time domain. The medium is one dimensional and its physical parameters are depth dependent, i.e., the medium is layered. The loss and dispersion properties of the medium are due to the fluid-structure interaction induced by wave propagation. They are modeled by generalized susceptibility functions which express the memory effects in the propagation process. The wave equation is then a fractional telegraphist’s equation. The two relevant quantities are the scattering operators—transmission and reflection operators—which give the scattered fields from the incident wave. They are obtained from Volterra equations which are fractional equations for the scattering operators.     

A Note on Perturbation Formulae for the Surface-Wave Speed Due To Perturbations in Material Properties

In a recent paper by Tanuma and Man, a two-term asymptotic formula was derived for the speed of surface waves propagating in an anisotropic elastic half-space whose elastic moduli differ only slightly from those for a (base) isotropic elastic material. This formula disagrees with that derived by Delsanto and Clark in an earlier paper using a different method. In this short note, we use a simple procedure to derive another two-term asymptotic formula for the surface-wave speed. Our formula takes the same compact form even if the base material is generally anisotropic. We show that when an error in the work of Delsanto and Clark is corrected, the three different methods do give equivalent results.      

根据散斑统计性质测量运动速度的一种新方法

本文提出了一种根据散斑统计性质测量物体运动速度的新方法.在夫朗和费衍射区一点记录运动物体的漫射散斑场随时间的变化,并将记录时城信号作快速付氏变换,转换到频率域上,从功率谱宽度或其所覆面积可推算出物体运动速度.本文从理论上作了分析,推导出物体的运动速度和记录信号功率谱函数的简单关系.实验验证与理论分析吻合一致.该法可作实时测量.     

Solution of acoustic scattering problems by means of second kind integral equations

In the present paper, the problem of scattering of acoustic waves from a fluid inclusion in two dimensions is solved numerically. The boundary conditions are those of continuous pressure and normal displacement. First, the problem in the frequency domain is reduced to a pair of second kind Fredholm integral equations on the boundary of the scatterer. Then, the integral equations are discretized by means of the Nyström algorithm, and the resulting linear systems are solved iteratively. Finally, the time domain solution is obtained from a sequence of frequency domain values. The integral equations of the present paper possess a very simple physical interpretation which guarantees the stability of their numerical solution and rapid convergence of the iterative solver. The resulting algorithm is an efficient tool for solving relatively large scale two-dimensional scattering problems.     

无网格局部强弱法求解不规则域问题

无网格局部彼得洛夫-伽辽金(meshless local Petrov-Galerkin,MLPG)法是一种具有代表性的无网格方法,在计算力学领域得到广泛应用.然而,这种方法在边界上需执行积分运算,通常很难处理不规则求解域问题.为了克服MLPG法的这种局限性,提出了无网格局部强弱(meshless local strong-weak,MLSW)法.MLSW法采用MLPG法离散内部求解域,采用无网格介点(meshless intervention-point,MIP)法施加自然边界条件,并采用配点法施加本质边界条件,避免执行边界积分运算,可适用于求解各类复杂的不规则域问题.从理论上讲,这种结合式方法,既保持了MLPG法稳定而精确计算的优势,同时兼备配点型方法在处理复杂结构问题时简洁而灵活的优势,实现了弱式法和强式法的优势互补.此外,MLSW法采用移动最小二乘核(moving least squares core,MLSc)近似法来构造形函数,是对传统移动最小二乘(moving least squares,MLS)近似法的一种改进.MLSc使用核基函数代替通常的基函数,有利于数值求解的精确性和稳定性,而且其导数近似计算变得更为简单.数值算例结果初步表明:这种新方法实施简单,求解稳定、精确,表现出适合工程运用的潜力.     

Dispersion and attenuation in thermoelastic multisize particulate composites

This paper deals with the problem of multiple scattering by a random distribution of spherical solid particles in a solid. The material properties of both media are taken as thermoelastic. The radii of the inclusions may be different. The self-consistent method in its variant of the effective medium is used to find the dispersion and attenuation of quasi-elastic, quasi-thermal and shear waves. The single scattering problem required by this technique is solved approximately by means of the Galerkin method applied to an integral equation using the Green function. Numerical results display a characteristic resonance phenomena which appears in the interval where the results are approximately valid, that is, for very long waves down to wavelengths about twice the largest diameter of the spheres. Examples are shown, for composites with two sets of inclusions, which have either a very similar or dissimilar size. Comparisons are made with the elastic counterpart. Among the material properties, the mass density ratio, inclusion to matrix, seems to play an important and simple role. Frequency intervals are distinguished and shown to depend on that ratio, where the attenuation and dispersion of quasi-elastic and P-waves are either very close to each other or not at all. The same applies to shear waves in either composite. The mass density ratio also displays a simple monotonic decreasing behaviour as a function of the frequency at the first attenuation maximum and velocity minimum. These results may be of interest for the nondestructive testing characterization of particulate composites.     

界面脱胶圆夹杂对SH波散射的远场解

采用Green函数方法和复变函数法研究了SH波对界面脱胶圆夹杂的散射问题,并给出 了远场解答. 首先,沿双质材料界面将整个空间分成上下两部分, 在下半空间,给出了在水 平表面上任意一点承受时间谐和的出平面线源载荷作用时的位移函数,取该位移函数作为 Green函数. 其次,在下半空间,利用相关文献给出的Green函数,在上下空间连接时在双质材料界面处满足 连续性条件,构造出半圆形脱胶裂纹,进而求出应力和位移的表达式,建立积分方 程组, 给出了散射波远场位移模式和散射截面的解答, 分析了在不同参数变化时SH 波散射的远场特性. 结果表明,脱胶结构的存在对位移和散射截面有较大的放大作用.     

有限散射信号下二维缺陷形状识别的罚函数方法

研究在有限照射角度和频带宽度下二维缺陷的形状识别问题。首先,通过引进介质参数扰动函数,建立介质参数扰动函数和弹性波散射场之间的非线性关系,并将所关心的缺陷的形状识别问题转化为关于扰动函数的反演;然后,利用变分技术和优化方法求解,为了弥补散射数据的不足,在总的目标函数中,采用附加度量函数作为罚函数;最后,对后场散射远场测量时有限照射角度和频带宽度下几种典型缺陷进行了模拟识别,表明了;表明了罚函数法的有效性。     

DYNAMIC INTERACTION OF PLANE WAVES WITH A UNILATERALLY FRICTIONALLY CONSTRAINED INCLUSION-TIME DOMAIN BOUNDARY ELEMENT ANALYSIS

A 2D time domain boundary element method (BEM) is developed to solve the transient scattering of plane waves by a unilaterally frictionally constrained inclusion. Coulomb friction is assumed along the contact interface. The incident wave is assumed strong enough so that localized slip and separation take place along the interface. The present problem is in effect a nonlinear boundary value problem since the mixed boundary conditions involve unknown intervals (slip, separation and stick regions). In order to determine the unknown intervals, an iterative technique is developed. As an example, we consider the scattering of a circular cylinder embeddedin an infinite solid.