1980年代中国地磁正常场图及其数学模式

本文介绍了1980.0年中国地磁正常场图的编绘和地磁正常场数学模式的建立使用不同年代观测的地磁三要素资料2000余个,经过通化改正,统统改正到1980.0这个特定年代。采用泰勒多项式和最小二乘法,分别建立1980.0年中国地磁正常场和地磁场长期变化数学模式,并计算其网格值,用于编绘1980.0年中国地磁正常场图。本文在建立地磁正常场数学模式时,采用三个独立的地磁要素建立模式的方法,解决了地磁场模式在地磁倾角为零的地方、而垂直强度不为零的问题。     

A survey of numerical methods for stochastic differential equations

The development of numerical methods for stochastic differential equations has intensified over the past decade. The earliest methods were usually heuristic adaptations of deterministic methods, but were found to have limited accuracy regardless of the order of the original scheme. A stochastic counterpart of the Taylor formula now provides a framework for the systematic investigation of numerical methods for stochastic differential equations. It suggests numerical schemes, which involve multiple stochastic integrals, of higher order of convergence. We shall survey the literature on these and on the earlier schemes in this paper. Our discussion will focus on diffusion processes, but we shall also indicate the extensions needed to handle processes with jump components. In particular, we shall classify the schemes according to strong or weak convergence criteria, depending on whether the approximation of the sample paths or of the probability distribution is of main interest.     

无单元法在三维断裂力学中的应用

介绍了无单元法在三维断裂力学问题中的应用。无单元法基于滑动最小二乘拟合,只需结点信息就可建立离散模型。采用罚函数法来满足无单元法的本征边界条件。用可视准则来处理裂纹面对高斯点影响域的隔离作用,推导了裂纹两表面的相互作用对整体平衡方程组的贡献方程。利用单点位移公式来计算三维裂纹前缘点的应力强度因子。通过算例验证了应用该法确定三维裂纹前缘应力强度因子的有效性。     

东亚地磁场模型的计算与分析

根据中国、前苏联和蒙古等地区的地磁复测点和地磁台站资料,使用泰勒多项式方法和冠谐分析方法,计算东亚地磁场(X,Y,Z)的泰勒多项式模型和冠谐模型,以及东亚剩余磁场(ΔX,ΔY,ΔZ)的冠谐模型和泰勒多项式模型,并绘制了相应的理论地磁图. 泰勒多项式模型的展开原点位于45°N和100°E,截断阶数为7.地磁场泰勒多项式模型的均方偏差为:X分量是133.0nT,Y分量是107.4nT,Z分量是14.0nT. 球冠极点位于45°N和100°E,球冠半角为42°,冠谐模型的截断阶数为10. 剩余磁场冠谐模型的均方偏差分别为131.2nT(ΔX),112.6nT(ΔY)和13.7nT(ΔZ). 对比分析了上述两种模型. 提出了确定区域模型截断阶数的判据.     

多孔介质中两相可压缩可混溶驱动的GALERKIN方法

多孔介质中两相可压缩可混溶驱动问题由一对耦合的非线性方程来描述.Douglas在忽略弥散效应的情况下,讨论了该方程的半离散Galerkin求解方法,并得到了解的最优的L2模误差估计.文章讨论了在考虑弥散效应的情况下该方程的半离散Galerkin方法,并得到了饱和度方程解的最优L2模误差估计.     

2003年中国及邻区地磁场模型的计算与分析

根据2003年中国地区的135个测点和35个台站的地磁数据, 建立了2003年中国及邻区地磁场泰勒多项式模型和球冠谐模型. 在模型计算过程中, 细致地分析了模型的截断阶数和边界约束对模型计算的影响. 结果表明, 所使用的地磁观测资料是准确可靠的, 模型截断阶数的选取是合理的. 本文所建立的中国区域地磁场球冠谐模型与泰勒多项式模型具有良好的一致性.      

Experimental investigation of stratorotational instability using a thermally stratified system: instability,waves and associated momentum flux

Stratorotational instability (SRI) has been proposed as a mechanism for outward angular momentum transport in Keplerian accretion disks. A particular designed Taylor–Couette laboratory experiment with axial stratification is suitable for studying the instability. Bottom endplate is cooled and top endplate is heated to achieve axial stratification. Due to constructive constraints, endplates are visually unamenable and quantitative measurement techniques in the co-rotating frame can only be done by looking through the outer cylinder. For this purpose, we built a co-rotating mini-PIV (Particle Image Velocimetry) system with a camera having a tilted viewing angle regarding the horizontal laser sheet. The aim of this study is (i) to quantify the uncertainty of the mini-PIV together with the used calibration technique and (ii) to compare experimental findings on SRI with theoretical predictions.

We perform measurements of the azimuthal and radial component of the velocity in axial stably stratified Taylor–Couette flows, consider velocity profiles and do frequency-filtering and flow decomposition. The absolute error of the mini-PIV system is 2% and we realised that stratified Taylor–Couette flows have smaller Ekman endwall effects than homogeneous ones. Still, Ekman pumping has an impact of the flow and might be responsible for differences between the data and theoretical models ignoring the endwalls. Here we focus on the flow structure during transition to SRI, the drift rate of SRI modes and the radial momentum flux as a function of the Reynolds number. Whereas the structure in form of trapped boundary Kelvin modes and the drift rate corresponds well with earlier predictions, the momentum flux shows a nonlinear dependency with respect to the Reynolds number. Away from the region of transition, theoretical models show a linear relationship. Several possible reasons for the mismatch between the experimental and theoretical models are discussed. Most important, we experimentally demonstrated that in the Rayleigh stable flow regime the SRI can provide a significant amount of outward momentum flux which makes this instability interesting in the context of accretion disks and also of atmospheric vortices where rotation and stratification also play a significant role.     

求解双相和黏弹性介质波传播方程的间断有限元方法及其波场模拟

间断有限元（Discontinuous Galerkin：DG）方法具有低数值频散、网格剖分灵活、能模拟地震波在复杂介质中传播等优点.因此,本文将一种新的DG方法推广到双相和黏弹性等复杂介质的地震波场模拟,发展了求解Biot弹性波方程和D'Alembert介质波动方程的DG方法.首先通过引入辅助变量将Biot双相介质弹性波方程和D'Alembert介质波动方程转化为关于时间-空间的一阶偏微分方程组,然后对该方程组进行DG空间离散,得到半离散化的常微分方程组.最后,对此常微分方程组,应用加权的Runge-Kutta格式进行时间推进计算.数值结果表明,DG方法可以有效地求解Biot双相介质弹性波方程和D'Alembert介质波动方程,并能很好地压制因离散求解波动方程而产生的数值频散,获得清晰的各种地震波震相.     

Contaminant transport with non‐equilibrium processes in unsaturated soils and implicit characteristic Galerkin scheme

A non‐equilibrium sorption—advection—diffusion model to simulate miscible pollutant transport in saturated–unsaturated soils is presented. The governing phenomena modelled in the present simulation are: convection, molecular diffusion, mechanical dispersion, sorption, immobile water effect and degradation, including both physical and chemical non‐equilibrium processes. A finite element procedure, based on the characteristic Galerkin method with an implicit algorithm is developed to numerically solve the model equations. The implicit algorithm is formulated by means of a combination of both the precise and the traditional numerical integration procedures. The stability analysis of the algorithm shows that the unconditional stability of the present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results illustrate good performance of the present algorithm in stability and accuracy, and in simulating the effects of all the mentioned phenomena governing the contaminant transport and the concentration distribution. Copyright © 2000 John Wiley & Sons, Ltd.     

二维浅水波方程的间断-连续Galerkin方法

将基于二维简化浅水波模型的间断Galerkin有限元与连续Galerkin有限元耦合方法推广至形式更为复杂的浅水波方程,并给出了误差分析以及模型问题的数值算例。