TMD结构基于双过滤白噪声随机地震激励的系列响应的简明解

针对既有方法在分析TMD结构基于双过滤白噪声激励下结构响应的解表达式复杂而导致计算效率低的问题,提出了一种简明封闭解法。首先,利用双过滤白噪声谱的滤波方程与TMD结构的地震动方程联立,可将TMD结构基于复杂的双过滤白噪声激励准确的表示为易于求解的运动方程;其次,基于复模态法获得TMD耗能结构位移、层间位移的系列响应的复特征值及复模态参与系数;然后基于随机振动理论获得了TMD结构随机地震动系列响应（相对于地面绝对位移和结构层间位移）的功率谱统一形式的二次正交解,进而获得了TMD结构系列随机响应的0-2阶谱矩和方差的简明封闭解。最后研究了基于首超破坏准、Markov过程假设及串联失效模式的TMD结构的体系动力可靠度。通过一算例分析,表明了本文方法的正确性和高效性。因此,本文方法可用于各类线性结构基于复杂的随机地震动响应的分析及其动力可靠度计算。     

Self-similar adiabatic flow headed by a magnetogasdynamic cylindrical shock wave in a rotating non-ideal gas

Similarity solutions are obtained for one-dimensional adiabatic flow behind a magnetogasdynamic cylindrical shock wave propagating in a rotating non-ideal gas in presence of an azimuthal magnetic field. The density of the medium ahead of the shock is assumed to be constant. In order to obtain the similarity solutions the angular velocity of the ambient medium is assumed to be obeying a power law and to be decreasing as the distance from the axis increases. It is found that the similarity solutions exist, in both the cases, when the initial magnetic field is constant or obeying a power law. The effects of an increase in the value of the index for variation of angular velocity of the ambient medium, in the value of the parameter of the non-idealness of the gas and in the strength of the initial magnetic field are obtained.     

Steady flows in solar magnetic solutions structures-a class of exact MHD solutions

Abstract

A class of exact solutions to the steady, two-dimensional magnetohydrodynamic equations ina cylindrical geometry is presented. These may model both closed and open magnetic structures found in the solar atmosphere. For closed structures, it is found that increasing the flow speed causes the summit of the arcade of closed magnetic fieldlines to rise. Parameter ranges also exist where the solution has regions of open and closed field, and so the solutions may be relevant for modelling flows in solar magnetic structures such as coronal streamers, X-ray bright points coronal plumes and coronal holes.     

黏性Cahn-Hilliard方程中时间周期解的存在性

本文研究一维具有变迁移率的黏性Cahn-Hilliard方程的时间周期解问题。采用了能量方法和Leray-Schauder不动点定理的方法,得到了解的一些先验估计值及Holder模估计值,证明了时间周期解的存在性。     

具有任意阶非线性薛定谔方程的新行波解

对具有任意阶非线性薛定谔方程的行波解进行了研究。给出具有任意阶非线性薛定谔方程的精确解。利用行波变换和辅助函数法把具有任意阶非线性薛定谔方程最终转化为一个非线性常微分方程的解,通过对这个微分方程的研究可以得到具有任意阶非线性薛定谔方程的行波解。具有任意阶非线性薛定谔方程存在的孤立波解、三角孤立波解、扭孤立波解、Jacobi椭圆函数解。得到了具有任意阶非线性薛定谔方程新的Jacobi椭圆函数解。     

二维SH波方程的半解析解及其数值模拟

本文以波动理论为基础, 半解析化求解地震勘探中常用的SH波方程. 获得的主要结果包括: 给出了二维均匀介质中SH波方程的解析解; 利用Cagniard-de Hoop方法详细推导了二维双层介质中SH波方程的解析解, 获得了透射波的解析解表达式. 同时, 基于SH波方程的解析表达式, 给出了包含各种波(如直达波、反射波、首波以及透射波)的解析解和波形图. 对于比较复杂的积分型解析解, 利用数值积分方法给出了数值结果, 并与优化的近似解析离散化方法(ONADM)和4阶Lax-Wendroff修正方法(LWC)的数值结果进行了比较, 以验证解析解的正确性. 本文的研究成果有望在检验波动方程数值新方法的有效性、波传播理论分析等方面得到应用.     

Galerkin boundary integral equation method for spontaneous rupture propagation problems: SH-case

We develop a Galerkin finite element boundary integral equation method (GaBIEM) for spontaneous rupture propagation problems for a planar fault embedded in a homogeneous full 2-D space. A 2-D antiplane rupture propagation problem, with a slip-weakening friction law, is simulated by the GaBIEM. This method allows one to eliminate the strong singularities from the integral representation of the traction, and to separate explicitly the expression for the traction into an instantaneous component; static and time-dependent components with weakly (logarithmic) singular kernels; and a dynamic component and a quasi-static component, with continuous, bounded, kernels. Simulated results throw light into the performance of the GaBIEM and highlight differences with respect to that of the traditional, collocation, boundary integral equation method (BIEM). Both methods converge with a power law with respect to grid size, with different exponents. There is no restriction on the CFL stability number for the GaBIEM since an implicit, unconditionally stable method is used for the time integration. The error of the approximation increases with the time step, as expected, and it can remain below that of the BIEM.     

Effect of 3-D viscoelastic structure on post-seismic relaxation from the 2004 M= 9.2 Sumatra earthquake

The 2004 M = 9.2 Sumatra–Andaman earthquake profoundly altered the state of stress in a large volume surrounding the ∼1400 km long rupture. Induced mantle flow fields and coupled surface deformation are sensitive to the 3-D rheology structure. To predict the post-seismic motions from this earthquake, relaxation of a 3-D spherical viscoelastic earth model is simulated using the theory of coupled normal modes. The quasi-static deformation basis set and solution on the 3-D model is constructed using: a spherically stratified viscoelastic earth model with a linear stress–strain relation; an aspherical perturbation in viscoelastic structure; a 'static' mode basis set consisting of Earth's spheroidal and toroidal free oscillations; a "viscoelastic" mode basis set; and interaction kernels that describe the coupling among viscoelastic and static modes. Application to the 2004 Sumatra–Andaman earthquake illustrates the profound modification of the post-seismic flow field at depth by a slab structure and similarly large effects on the near-field post-seismic deformation field at Earth's surface. Comparison with post-seismic GPS observations illustrates the extent to which viscoelastic relaxation contributes to the regional post-seismic deformation.     

New advances in three-dimensional controlled-source electromagnetic inversion

New techniques for improving both the computational and imaging performance of the three-dimensional (3-D) electromagnetic inverse problem are presented. A non-linear conjugate gradient algorithm is the framework of the inversion scheme. Full wave equation modelling for controlled sources is utilized for data simulation along with an efficient gradient computation approach for the model update. Improving the modelling efficiency of the 3-D finite difference (FD) method involves the separation of the potentially large modelling mesh, defining the set of model parameters, from the computational FD meshes used for field simulation. Grid spacings and thus overall grid sizes can be reduced and optimized according to source frequencies and source–receiver offsets of a given input data set. Further computational efficiency is obtained by combining different levels of parallelization. While the parallel scheme allows for an arbitrarily large number of parallel tasks, the relative amount of message passing is kept constant. Image enhancement is achieved by model parameter transformation functions, which enforce bounded conductivity parameters and thus prevent parameter overshoots. Further, a remedy for treating distorted data within the inversion process is presented. Data distortions simulated here include positioning errors and a highly conductive overburden, hiding the desired target signal. The methods are demonstrated using both synthetic and field data.