Computer Simulation of Two-dimensional Random Wave Propagation

A straightforward numerical algorithm for simulating the propagationof a wavefield in a two-dimensional randomly varying mediumis described. Both a finite-difference and a fast Fourier transformmethod are described. These methods are well known, but thereis a novel treatment of the random scattering term in the waveequation that allows accurate answers to be economically obtained.The methods are then tested by comparison with known approximatesolutions for the fourth moment of a propagating wavefield.The two approaches show good agreement, thus confirming theusefulness of the analytic results and at the same time indicatingthat the simulation process should be a powerful tool for investigatingthe higher-order statistics of the field. The agreement shouldhold in situations encountered in optical and acoustic scatteringexperiments.     

地形构造中地震波传播的非对称交错网格模拟

提出了一种新的三维空间对称交错网格差分方法,模拟地形构造中弹性波传播过程．通过具有二阶时间精度和四阶空间精度的不规则网格差分算子用来近似一阶弹性波动方程,引入附加差分公式解决非均匀交错网格的不对称问题．该方法无需在精细网格和粗糙网格间进行插值,所有网格点上的计算在同一次空间迭代中完成．使用精细不规则网格处理海底粗糙界面、 断层和空间界面等复杂几何构造, 理论分析和数值算例表明, 该方法不但节省了大量内存和计算时间, 而且具有令人满意的稳定性和精度．在模拟地形构造中地震波传播时,该方法比常规方法效率更高．     

Simulation of Wave Propagation in Complex Structures with the Spectral Element Method

In nondestructive testing, the use of ultrasonic elastic waves has proven as one of the most successful principles to detect structural damage like cracks, delaminations etc. Especially, Structural Health Monitoring (SHM) is characterized by permanently installed embedded or surface–mounted actuators and sensors (e.g. piezoelectric patches). The capability of most approaches strongly depends on adequate choice of parameters like excitation signals and actuator/sensor positions. For this reason there is a growing interest in efficient and accurate simulation tools to shorten time and cost of the necessary pretests. With respect to high frequency excitation a computationally efficient method is required. This contribution presents the theoretical background of the spectral element method including the electro–mechanical coupling of piezo elements. The spectral element method generates a diagonal mass matrix leading to significant savings of memory and to a crucial reduction of complexity of the time integration algorithm. Both in–plane and out–of–plane waves can be handled. Numerical examples for the propagation of waves in stiffened structures are presented. The effect of improper placement of actuators/sensors is shown. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Efficient Simulation of Wave Propagation with Implicit Finite Difference Schemes

Finite difference method is an important methodology in the approximation of waves. In this paper, we will study two implicit finite difference schemes for the simulation of waves. They are the weighted alternating direction implicit (ADI) scheme and the locally one-dimensional (LOD) scheme. The approximation errors, stability conditions, and dispersion relations for both schemes are investigated. Our analysis shows that the LOD implicit scheme has less dispersion error than that of the ADI scheme. Moreover, the unconditional stability for both schemes with arbitrary spatial accuracy is established for the first time. In order to improve computational efficiency, numerical algorithms based on message passing interface (MPI) are implemented. Numerical examples of wave propagation in a three-layer model and a standard complex model are presented. Our analysis and comparisons show that both ADI and LOD schemes are able to efficiently and accurately simulate wave propagation in complex media.     

侧向加热腔体中的多圈型对流斑图

基于流体力学方程组的数值模拟,研究了倾角θ=90°时侧向加热的大高宽比腔体中的对流斑图.对于Prandtl数Pr=6.99的流体,在相对Rayleigh数2≤Ra r≤25的范围内,腔体中发生的是单圈型对流斑图.对于Pr=0.0272的流体,取Ra r=13.9,随着计算时间的发展,腔体中由最初的单圈型对流斑图过渡到多圈型对流斑图,这是出现在侧向加热大高宽比腔体中的新型对流斑图.对不同Ra r情况的计算结果表明,Ra r对对流斑图的形成存在明显的影响.当Ra r≤4.4时是单圈型对流滚动;当Ra r=8.9~11.1时是过渡状态;当Ra r≥13.9时是多圈型对流滚动.对流最大振幅和Nusselt数Nu随着相对Rayleigh数的增加而增加.该对流斑图与Pr=6.99时对流斑图的比较说明,对流斑图的形成依赖于Prandtl数.     

Numerical Simulation of Detonation Excitation at Shock Wave Interaction with Dust Cloud in Channel

The study of detonation ability of reactive particle gas mixtures is necessary to prevent industrial explosions in industries where dispersed powders are used. The present paper focuses on numerical simulation of the shock wave interaction with semiinfinite aluminum dust cloud, which is situated inside a plane channel. The cloud fills entirely or partly the channel cross‐section and has initially a rectangular shape. The scenarios of detonation initiation in the cloud are determined depending on the incident shock wave amplitude values. The processes of transformation and spreading of finite width clouds under weak incident shock wave action (when the particles do not ignite) are investigated. The types of an oblique shock wave reflection from the plane of symmetry in the cloud are analyzed. The processes of particle ignition and detonation structure formation at strong incident shock wave action are investigated. Nonstationary periodic fuctuations take place in the detonation flow due to transversal wave effect. Nevertheless the detonation structure established propagates in quasistationary regime. If the incident shock wave is attenuated with a rarefaction wave then the detonation formation fails at clouds of insufficient width.     

流体饱和多孔隙介质波动方程小波有限差分法

研究流体饱和多孔隙介质中波动方程的数值模拟．针对求解二维弹性波方程问题,提出小波有限差分法．该方法综合了小波多分辨分析计算灵活、计算效率高特性和有限差分易于实现的优点．数值模拟的结果显示,此方法对于求解流体饱和多孔隙介质方程的数值模拟是有效稳定的．     

弹性波方程有限元模拟的曲边地表处理

本文对比研究了关于弹性波模拟中的曲边地表形状处理的两种方法,一种是用给定的实际介质数值划定的地表形状,另一种是用样条插值逼近地表形状.本文采用有限元方法进行弹性波数值模拟,给出了基于这两种方法计算的数值例子,并对结果进行了分析比较.结果表明使用后一种方法对地表进行处理时,地表人工离散产生的干扰明显减少,优于前一种方法.     

Study on Discharge Coefficients of Drain Orifices in Closed Cavities北大核心CSCD

该文利用数值方法模拟了封闭腔体的排液流动,获取了排液孔的流量系数.通过量纲分析法研究了小孔流量系数的主要影响因素,拟合了计算小孔流量系数的经验公式.结果表明：当水头高度小于200 mm时,小孔流量系数随水头高度的增加而减小;当水头高度大于200 mm时,小孔流量系数稳定在0.61附近.不同厚径比的小孔流量系数表现为两种不同的形式：小厚径比的小孔呈现薄孔流动特性,流量系数为0.6左右;大厚径比的小孔呈现厚孔流动特性,流量系数为0.8左右.     

Cavities in a visco-plastic material: A mesoscopic concept

In the present work we investigate the temporal development of arbitrarily distributed voids in a visco-plastic material under different loading regimes. A mesoscopic continuum model is used in order to take the microstructure of the material into account. In particular, we introduce a mesoscopic space representing an extension of the space-time domain of the continuum mechanical fields. This extended domain requires a reformulation of the classical balance equations as well as the consideration of additional constitutive quantities. Furthermore a mesoscopic distribution function can be introduced which follows an own balance. Assuming a special model of porous composites, the spherical shell model, all required steps are elaborated in order to describe load-induced void-growth in a metal-like matrix. We conclude with some exemplary results which show astonishing similarities with co-called LSW-theories. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Absorbing Interface Conditions for the Simulation of Wave Propagation on Non-Uniform Meshes

We proposed absorbing interface conditions for the simulation of linear wave propagation on non-uniform meshes. Based on the superposition principle of second-order linear wave equations, we decompose the interface condition problem into two subproblems around the interface: for the first one the conventional artificial absorbing boundary conditions is applied, while for the second one, the local analytic solutions can be derived. The proposed interface conditions permit a two-way transmission of low-frequency waves across mesh interfaces which can be supported by both coarse and fine meshes, and perform a one-way absorption of high-frequency waves which can only be supported by fine meshes when they travel from fine mesh regions to coarse ones. Numerical examples are presented to illustrate the efficiency of the proposed absorbing interface conditions.     

Some Basic Ideas for the Simulation of Wave Propagation in Microstructures using Proper Orthogonal Decomposition

For the modeling of micro-heterogeneous materials, models are required that take into account the microstructure. More and more, computational homogenization is used to quantify the effective macroscopic response. The large number of resulting microscopic boundary value problems need to be solved efficiently. To reduce the computational effort and simulation time, model reduction methods may be applied. In this contribution the performance of the proper orthogonal decomposition is analyzed for the wave propagation in fiber reinforced materials. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

解洞穴逆散射问题的一种正则化Newton方法

We consider the inverse problem to determine the shape of a open cavity embedded in the infinite ground plane from knowledge of the far-field pattern of the scattering of TM polarization.For its approximate solution we propose a regularized Newton iteration scheme.For a foundation of Newton type methods we establish the Fréchet differentiability of solution to the scattering problem with respect to the boundary of the cavity.Some numerical examples of the feasibility of the method are presented.     

A Regularized Newton Method in Inverse Scattering Problem from Cavities

We consider the inverse problem to determine the shape of a open cavity embedded in the infinite ground plane from knowledge of the far-field pattern of the scattering of TM polarization.For its approximate solution we propose a regularized Newton iteration scheme.For a foundation of Newton type methods we establish the Fréchet differentiability of solution to the scattering problem with respect to the boundary of the cavity.Some numerical examples of the feasibility of the method are presented.     

物体垂直出入水的非定常空泡数学模型

研究了物体从水下向水面高速运动产生的非定常垂直空泡,建立了出水垂直空泡的数学模型,得到了匀速、空泡压力不变条件下的空泡外形、长度、体积随水深变化的解析解,给出了出水通气空泡发展为超空泡的条件．利用类似的方法建立了水平空泡和入水空泡的数学模型,并对3种空泡进行了比较研究．比较研究的结论是,随物体距水面的水深减小,出水空泡体积自身有增大趋势,空泡不容易发生泄气现象．且要保持出水空泡压力不变,空泡内的气体含量就应该增加（可通过人工通气方式）．入水空泡正好相反,随物体入水深度增加,空泡体积自身有收缩的趋势,并挤压空泡内的气体从环境压力较低的空泡尾部喷射而出,导致空泡内的气体含量减少,空泡压力降低．但是当空泡压力低于环境压力后,空泡尾部又将被环境高压所封闭,气体喷射不出来．随着入水深度继续增加,空泡尾部将重复上述过程,形成周期性的喷射 封闭 喷射 封闭的脉动过程,这个脉动喷射过程将在空泡尾部的流体中形成一连串小气泡,并由于空泡内的压力波动而导致空泡形状发生波动现象．