Tests and applications of an approach to absorbing reflected waves towards incident boundary

If the upstream boundary conditions are prescribed based on the incident wave only, the time-dependent numerical models cannot effectively simulate the wave field when the physical or spurious reflected waves become significant. This paper describes carefully an approach to specifying the incident wave boundary conditions combined with a set sponge layer to absorb the reflected waves towards the incident boundary. Incorporated into a time-dependent numerical model, whose governing equations are the Boussinesq-type ones, the effectiveness of the approach is studied in detail. The general boundary conditions, describing the down-wave boundary conditions are also generalized to the case of random waves. The numerical model is in detail examined. The test cases include both the normal one-dimensional incident regular or random waves and the two-dimensional oblique incident regular waves. The calculated results show that the present approach is effective on damping the reflected waves towards the incident wave boundary.     

浅剖数据解译中淤泥层层界提取方法

港口与航道淤泥层层界确定对于港航建设开发与工程实施有着重大意义。本文针对港口航道浅地层剖面仪原始采集数据模糊的问题,分析了影响原始数据真实性的主要原因,实施了中值滤波法多次波压制与消噪处理,取得了清晰可靠的真实数据图像文件,根据海底声波不同层界反射特征,对比分析了脉冲信号声强图像变化情况,探讨了基于浅剖数据声强图像提取淤泥层厚的方法。试验证明该方法可以达到当前钻孔资料所能达到的精度范围。     

模拟地震波传播的三维逐元并行谱元法

高效地震波场正演模拟对于复杂模型中地震波传播与成像研究至关重要.本文在谱元法原理框架内,对已有逐元谱元法改进,提出一种新的逐元并行谱元法求解三维地震波运动方程,并得到地震波场.逐元并行谱元法的核心思想在于在单元上进行质量矩阵与解向量的乘积运算,并将此运算平均分配至每一个CPU计算核心,此处理有利提升谱元法的并行计算效率.同时,根据Gauss-Lobatto-Legendre (GLL)数值积分点与插值点重合的特点,将稠密单元刚度矩阵的存储转化成单元雅克比矩阵行列式的值及其逆的存储,大幅减少谱元法计算内存开销.此外,在模型边界上利用逐元并行谱元法求解二阶位移形式完美匹配层(PML)吸收边界条件,消除边界截断而引入的虚假反射.通过逐元并行谱元法得到的数值解与解析解对比,以及实际地震波场模拟,数值结果证实了逐元并行谱元法用于地震波场模拟的高效性.

     

Application of Dynamic Subgrid-scale Models for Large-eddy Simulation of the Daytime Convective Boundary Layer over Heterogeneous Surfaces

The sensitivity of large-eddy simulation (LES) to the representation of subgrid-scale (SGS) processes is explored for the case of the convective boundary layer (CBL) developing over surfaces with varying degrees of spatial heterogeneity. Three representations of SGS processes are explored: the traditional constant Smagorinsky–Lilly model and two other dynamic models with Lagrangian averaging approaches to calculate the Smagorinsky coefficient (C _S ) and SGS Prandtl number (Pr). With initial data based roughly on the observed meteorology, simulations of daytime CBL growth are performed over surfaces with characteristics (i.e. fluxes and roughness) ranging from homogeneous, to striped heterogeneity, to a realistic representation of heterogeneity as derived from a recent field study. In both idealized tests and the realistic case, SGS sensitivities are mostly manifest near the surface and entrainment zone. However, unlike simulations over complex domains or under neutral or stable conditions, these differences for the CBL simulation, where large eddies dominate, are not significant enough to distinguish the performance of the different SGS models, irrespective of surface heterogeneity.     

Paleo‐fires and Atmospheric Oxygen Levels in the Latest Permian: Evidence from Maceral Compositions of Coals in Eastern Yunnan,Southern China

Inertinite maceral compositions of the Late Permian coals from three sections in the terrestrial and paralic settings of eastern Yunnan are analyzed in order to reveal the paleo-fire events and the atmospheric oxygen levels in the latest Permian. Although the macerals in the studied sections are generally dominated by vitrinite, the inertinite group makes up a considerable proportion. Its content increases upward from the beginning of the Late Permian to the coal seam near the Permian-Triassic boundary. Based on the microscopic features and the prevailing theory that inertinite is largely a by-product of paleo-fires, we suggest that the increasing upward trend of the inertinite abundance in the latest Permian could imply that the Late Permian peatland had suffered from frequent wildfires. Since ignition and burning depend on sufficient oxygen, a model-based calculation suggests that the O2 levels near the Wuchiapingian/Changhsingian boundary and the Permian-Triassic boundary are 27% and 28% respectively. This output adds supports to other discoveries made in the temporal marine and terrestrial sediments, and challenges the theories advocating hypoxia as a mechanism for the Permian-Triassic boundary crisis.     

On determining the noon polar cap boundary from SuperDARN HF radar backscatter characteristics

Previous work has shown that ionospheric HF radar backscatter in the noon sector can be used to locate the footprint of the magnetospheric cusp particle precipitation. This has enabled the radar data to be used as a proxy for the location of the polar cap boundary, and hence measure the flow of plasma across it to derive the reconnection electric field in the ionosphere. This work used only single radar data sets with a field of view limited to 2 h of local time. In this case study using four of the SuperDARN radars, we examine the boundary determined over 6 h of magnetic local time around the noon sector and its relationship to the convection pattern. The variation with longitude of the latitude of the radar scatter with cusp characteristics shows a bay-like feature. It is shown that this feature is shaped by the variation with longitude of the poleward flow component of the ionospheric plasma and may be understood in terms of cusp ion time-of-flight effects. Using this interpretation, we derive the time-of-flight of the cusp ions and find that it is consistent with approximately 1 keV ions injected from a subsolar reconnection site. A method for deriving a more accurate estimate of the location of the open-closed field line boundary from HF radar data is described.     

TEM study of a special grain boundary in a synthetic K-feldspar bicrystal: Manebach Twin

 A monoclinic KAlSi₃O₈ feldspar Manebach twin boundary was synthesized by diffusion bonding and examined using high-resolution transmission electron microscopy. The sharp (001) twin boundary is straight and free of strain. The boundary width is smaller than d₀₀₁. There is no rigid body shift observed at the twin boundary, and the feldspar structure is arranged symmetrically across (001). The twin boundary structure consists of bridged tetrahedral crankshafts, which are characteristic of the feldspar lattice. The grain boundary structure is in good agreement with the geometrical model of Taylor et al. (1934). The grain boundary composition of K_1/2H_1/2AlSi₃O₈ differs from their model. Received: 13 February 2002 / Accepted: 24 December 2002 Acknowledgements We thank M. Rühle, S. Hutt, J. Mayer, A. Strecker and U. Salzberger at MPI, Stuttgart, for their support and valuable advice in preparing TEM sections of bicrystals.     

Magma underplating and Hannuoba present crust-mantle transitional zone composition: Xenolith petrological and geochemical evidence

Studies on the deep-seated xenoliths from global volcanoes reveal that the present petrological crust-mantle boundary between the lower crust and the upper mantle actually is a transitional layer from mainly mafic granulites to ultramafic spinel lher-zolites[1,2], i.e. a transitional zone distinctive from the seismological Moho[3]. Oceanic lithosphere crust- mantle transitional zone can be established from the study on the exposed ophiolites. However, as for the continental lithosphere, since …     

鄂尔多斯盆地南缘上二叠统——中下三叠统地球化学特征及其古气候、古环境指示意义

为探索陆相湖泊环境记录中二叠纪—三叠纪之交古气候演化的信息,以鄂尔多斯盆地南缘陕西铜川石川河剖面上二叠统石千峰组(P3s)和三叠系的刘家沟组(T1l)、和尚沟组(T1h)、纸坊组(T2z)为研究对象,对界面上下地层开展了系统的矿物学、岩石学和地球化学研究。通过主量和微量元素、碳氧同位素以及TOC测试对古盐度、古氧化还原环境和古气候演化规律进行分析。测试结果显示在二叠系—三叠系(P/T)界线附近主、微量元素和碳、氧同位素发生较明显的变化,气候环境代用化学指标的波动指示了从晚二叠世至早中三叠世鄂尔多斯盆地的古气候、古环境的变化。研究结果表明上二叠统石千峰组形成于河流-三角洲沉积环境,氧化程度相对低,古气候相对温暖湿润;下三叠统刘家沟组与和尚沟组的紫红色砂泥岩代表水体较浅的河流-三角洲相,强氧化环境,气候干旱炎热;而到中三叠统纸坊组沉积期,湖平面上升,元素的迁移作用加强,氧化程度变弱,气温降低,气候转为半干旱半湿润气候。碳同位素分析结果显示,鄂尔多斯盆地陆相P/T界线上δ^13CPDB存在显著负偏,与全球范围内的海相P/T界线具有一致性,说明在华北地台陆相P/T界线上同样存在气候突变和生物灭绝等重大地质事件的沉积记录,与海相地层可对比。     

移动状态边界面法在黏性土不排水蠕变孔压分析中的应用

在黏性土的各种流变现象的数学模拟研究中,不排水蠕变孔压方程几乎是空白。基于临界状态土力学理论,采用与时间相关的移动状态边界面法,给出了适用于正常固结黏性土的不依赖状态边界面具体形态的不排水蠕变孔压公式。正常固结黏性土的不排水蠕变孔压可以看作是体积蠕变势和剪应力水平增大造成的剪缩这两种作用的结果,前者与参数 有关（ 为次压缩指数, 为压缩指数）,后者则与状态边界面的具体形态有关。在此基础上进一步给出了可方便用于试验数据分析的双对数和单对数形式的两个公式。以修正剑桥模型椭圆形状态边界面为例,分析了双对数和单对数坐标系下蠕变孔压曲线的斜率m和m′与土的流变参数 的关系。采用提出的方法,很好地预测了Walkers给出的Leda clay黏土的蠕变孔压系数,并基于以上两种作用解释了Leda clay黏土的蠕变孔压形成的机制。