一般曲线坐标系下波浪传播的数值模拟

在曲线坐标系下,建立了缓变水深水域波浪传播的数值模拟模型.模型适宜于复杂变化的边界形状,克服了各种代数坐标变换的局限性.在建立模型时,将原始的椭圆型缓坡方程的近似型式——依赖时间变化的抛物型方程,作为控制方程,既克服了一般抛物近似方法的缺点,又便利了方程的求解;从开边界条件、不同反射特性的固壁边界条件相统一的表达式出发,对边界条件进行处理;用ADI法数值求解控制方程.对模型的验证表明,数值解与物模实验值吻合良好,模型对于具有复杂边界的工程实际有较强的适应性.     

ECMWF细网格模式探空在阿勒泰地区短时强降水预报中的统计检验

采用2013—2018年5—9月ECMWF细网格资料和阿勒泰地区36次短时强降水资料,用Micaps平台的模式探空模块计算T-log P图及其对流参数,运用统计学方法进行了误差检验。结果表明:模式探空T-log P图48 h预报时效内一致性较高(72%),尤其是24 h预报时效内(92%);72 h预报时效内总指数和干暖盖指数及垂直风切变、60 h预报时效内沙氏指数、48 h预报时效内风暴相对螺旋度和36 h预报时效内850与500 hPa温度差、对流温度、最大抬升指数、抬升指数以及24 h预报时效内K指数、700与850 hPa假相当位温差、大风指数等对流参数的3种误差均较小(3.5),相关系数较高(0.60),特别是沙氏指数、K指数、700与850 hPa假相当位温差、垂直风切变和风暴相对螺旋度的3种误差1.5。T-log P图的一致性随时效的延长而减小,对流参数的误差随预报时效的延长变化不一致,在强对流潜势预报业务中注意订正运用。     

Spatial variability of infiltration and related properties of tropical soils

The levels of variance associated with measuring the infiltration process and modelling it by means of a regression model are compared to see which approach yields the best results in terms of effort and accuracy. A nested sampling scheme has been used in the three major physiographic units of central Guyana, South America: ‘White Sands’; (Haplic and Ferralic Arenosols), ‘Brown Sands’ (Haplic Ferrasols) and ‘Laterite’ (Xanthic and Dystric Leptosols). Cluster analysis yields three sample groups that reflect the sharp landscape boundaries between the units. Multiple regression analysis shows that each unit has a different combination of soil properties that explains the variance in final infiltration rate and sorptivity satisfactorily. Nested analysis of variance indicates that clear spatial patterns with distances of variation of several hundred metres exist for final infiltration rate in White Sands and Laterite. Infiltration rate in Brown Sands and sorptivity in all units have large short-distance variabilities and high ‘noise’ levels. The correlated independent variables behave accordingly. For the majority of the soil properties, sampling at distances of 100 to 200 m results in variance levels of more than 80 per cent of the total variance, which indicates that only a detailed investigation can assess spatial variation in soil hydrological behaviour. The use of simple soil properties to predict infiltration is only possible in a very general sense and with the acceptance of high variance levels.     

关于台风非对称结构与台风路径的数值分析

采用一个有限区域原始方程双向相互作用的移动套网格台风路径数值预报模式，对１９９２年的两个台风进行了５个个例预报试验，结果表明预报路径和实况有很好的一致性，通过对预报个例的台风非对称结构进行分析，发现非对称结构随时间变化与台风移动路径有着比较密切的联系，如能了解实际台风的非对称结构，将有助于提高台风路径预报的水平。     

P-T grids for silica-undersaturated granulites in the systems MAS (n+ 4) and FMAS (n+ 3)-tools for the derivation of P-T paths of metamorphism

Abstract Considering the minerals cordierite (Cd), sapphirine (Sa), hypersthene (Hy), garnet (Ga), spinel (Sp), sillimanite (Si) and corundum (Co) in the system FeO-MgO-Al₂O₃-SiO₂ (FMAS), the stable invariant points are [Co], [Ga], [Cd] and [Sa]. Constraints imposed by experimental data for the system MAS indicate that under low P _H2o conditions the invariant points occur at high temperature (> 900° C) and intermediate pressure (7-10 kbar). This temperature is higher than that commonly advocated for granulite facies metamorphism. In granulites Fe-Mg exchange geothermometers may yield temperatures of 100–150° C below peak metamorphic conditions and evidence for peak temperatures is best preserved by relict high-temperature assemblages and by Al-rich cores in orthopyroxene. Application of the FMAS grid to some well-documented granulite occurrences introduces important constraints on their P-T histories. Rocks of different bulk compositions, occurring in close proximity in the field, may record distinct segments of their P-T paths. This applies particularly to rocks with evidence for reaction in the form of coronas, symplectites and zoned minerals. Consideration of curved reaction boundaries and _XMs isopleths may explain apparently contradictory results for the stability of cordierite obtained from low-temperature experiments and thermochemical calculations on the one hand and hightemperature experimental data on the other.     

无结构网格上平面二维水沙模拟的有限体积法

基于无结构网格有限体积法的算法框架,通过引入跨单元界面法向水沙数值通量的逆风分解,将悬沙与床沙交换以及分组挟沙力计算模式自然地嵌入二维水沙运动方程组的数值格式中,形成高精度、守恒性好的二维水沙有限体积算法。最后,利用该算法对谭江樟州河段的水沙输运和河床变形进行了数值模拟。结果表明,该算法能够较好地模拟复杂条件下河道水沙输运的往复特征和河床变形的动态过程,其精度满足河道工程后效分析的要求。     

非静力模式中几种垂直网格计算特性的对比研究

从非静力滞弹性方程组出发，将惯性重力波的水平波长分几百公里、几十公里和几公里三种情况，从频率、群速和各速的垂直尺度范围等方面对现有的几种垂直网格进行比较讨论。结果表明：在非静力模式中，CP网格和LZ网格适合于各种不同水平尺度惯性重力波的模拟；L网格和LY网格适合于水平尺度为几十公里以上的惯性重力波的模拟；LTS网格和CPTS网格比较适合于水平尺度的几十公里以下的惯性重力波的模拟。     

基于正六边形格网的最短路径算法

本文在分析了现有算法的一些不足之处的基础上,结合正六边形的特点及水流扩散思想,提出了基于正六边形格网的最短路径分析算法。该算法在最短路径搜索过程中,对同一正六边形格网而言,它至起点的累计代价值,不需要进行数据比较和修正。与经典的Dijikstra算法相比,该算法大大节约了搜索的时间。     

Adaptive unstructured grid finite element simulation of two-dimensional magnetotelluric fields for arbitrary surface and seafloor topography

We present an adaptive unstructured triangular grid finite element approach for effectively simulating plane-wave diffusive electromagnetic fields in 2-D conductivity structures.
The most striking advantage of irregular grids is their potential to incorporate arbitrary geometries including surface and seafloor topography. Adaptive mesh refinement strategies using an a posteriori error estimator yield most efficient numerical solutions since meshes are only refined where required.
We demonstrate the robustness of this approach by comparison with analytical solutions and previously published numerical simulations. Maximum errors may systematically be reduced to, for example, 0.8 per cent for the apparent resistivity and 0.2° in the phase.
An additional accuracy study of the thickness of the air layer in E-polarization suggests to keep a minimum thickness depending on lateral conductivity contrasts within the earth.
Furthermore, we point out the new quality and flexibility of our simulation technique by addressing two marine magnetotelluric applications. In the first case, we discuss topographic effects associated with a synthetic sinusoidal sea bottom model and in the second case, we show a close-to-reality scenario using real bathymetry data from the East Pacific Rise at 17°S.