SEAWAT 2000: modelling unstable flow and sensitivity to discretization levels and numerical schemes

A systematic analysis shows how results from the finite difference code SEAWAT are sensitive to choice of grid dimension, time step, and numerical scheme for unstable flow problems. Guidelines to assist in selecting appropriate combinations of these factors are suggested. While the SEAWAT code has been tested for a wide range of problems, the sensitivity of results to spatial and temporal discretization levels and numerical schemes has not been studied in detail for unstable flow problems. Here, the Elder-Voss-Souza benchmark problem has been used to systematically explore the sensitivity of SEAWAT output to spatio-temporal resolution and numerical solver choice. A grid size of 0.38 and 0.60% of the total domain length and depth respectively is found to be fine enough to deliver results with acceptable accuracy for most of the numerical schemes when Courant number (Cr) is 0.1. All numerical solvers produced similar results for extremely fine meshes; however, some schemes converged faster than others. For instance, the 3rd-order total variation-diminishing method (TVD3) scheme converged at a much coarser mesh than the standard finite difference methods (SFDM) upstream weighting (UW) scheme. The sensitivity of the results to Cr number depends on the numerical scheme as expected.     

Coarsening of three-dimensional structured and unstructured grids for subsurface flow

We present a generic, semi-automated algorithm for generating non-uniform coarse grids for modeling subsurface flow. The method is applicable to arbitrary grids and does not impose smoothness constraints on the coarse grid. One therefore avoids conventional smoothing procedures that are commonly used to ensure that the grids obtained with standard coarsening procedures are not too rough. The coarsening algorithm is very simple and essentially involves only two parameters that specify the level of coarsening. Consequently the algorithm allows the user to specify the simulation grid dynamically to fit available computer resources, and, e.g., use the original geomodel as input for flow simulations. This is of great importance since coarse grid-generation is normally the most time-consuming part of an upscaling phase, and therefore the main obstacle that has prevented simulation workflows with user-defined resolution. We apply the coarsening algorithm to a series of two-phase flow problems on both structured (Cartesian) and unstructured grids. The numerical results demonstrate that one consistently obtains significantly more accurate results using the proposed non-uniform coarsening strategy than with corresponding uniform coarse grids with roughly the same number of cells.     

水平气压梯度力差分格式误差的定量分析

通过分析发现地形坐标系中气压梯度力差分格式的计算误差应分为二类，其中，第二类误差在水平坐标面倾斜时出现。这种误差可理解为将水平坐标面上的气压(或位势高度)插回到等高面(或等压面)上的插值误差。 构造气压梯度力差分格式的目的正是为了减小这种第二类误差，而不是其他误差。 误差分析表明静力扣除法和Corby 格式的第二类误差都很小，比一般方法的误差小约一个量级。分析还表明第二类误差产生的主要原因是由于气压随高度非线性变化。依据这一分析思路，本文改进了#As(z)#a坐标系中的静力扣除法。 改进后的方法比改进前的误差小约一个量级。     

基于空间信息网格（SIG）的矿产资源评价

随着GIS的普及和深化,空间信息的大数据量、高复杂性和难共享的问题对其研究和应用提出了新的挑战.空间信息网格(Spatial Information Grid,SIG)作为创新性的网络空间信息基础设施和技术体系,为解决这些问题提供了新的契机.本文在给出SIG基本定义和基本特征的基础上,就如何实现基于SIG的矿产资源评价进行了有益的探索和实践,并取得了很好的效果,为矿产资源评价提供了新的方法和思路.     

Interpretation of a two-point flux stencil for skew parallelogram grids

Control-volume formulations for elliptic equations often use two-point flux stencils, even for skew grids. Any two-point flux stencil may be interpreted as a multipoint flux stencil. This yields a definition of the permeability (or conductivity) tensor. Formulas for calculating the permeability tensor, based on the user-specified quantities in the two-point flux stencil, are given. Numerical test examples demonstrate the validity of the derivation.      

虚拟天文台数据访问系统(VO-DAS)任务调度设计与实现

为解决异地、异构、海量天文数据无缝透明地统一访问,中国虚拟天文台(China-VO)研究小组自主设计并开发了一套解决方案——虚拟天文台数据访问系统(VO-DAS)。着重从执行流程、设计模式、Session机制、生命周期、资源销毁、异常处理等模块详细阐述了该系统的任务调度设计。为了验证该系统的可行性和性能,最后描述了多个天文科学范例。     

内蒙古自治区2000国家大地坐标框架的建立及实时服务

启用2000大地坐标系的必要性;内蒙古2000国家大地坐标框架建立。     

基于GIS技术设计实现西安移动公司网格化管理系统

移动公司网格化管理是城市网格化管理的重要内容,设计移动公司网格化管理系统是实现网格管理的必要条件。根据GIS技术的工具型特性,采用GIS集中二次开发方案、Window XP Professional操作系统和Visual C#语言,以Delphi,Visual Studio等通用软件为开发工具,集成开发了西安移动公司网格管理系统。开发的网格管理系统实现了西安移动公司属性数据和空间数据的系统管理,可方便地进行地图显示、数据查询、网络浏览等功能操作,有助于公司管理人员进行信息管理和实施分析决策。     

Capturing the essential spatial variability in distributed hydrological modelling: Infiltration parameters

Selecting the correct resolution in distributed hydrological modelling at the watershed scale is essential in reducing scale-related errors. The work presented herein uses information content (entropy) to identify the resolution which captures the essential variability, at the watershed scale, of the infiltration parameters in the Green and Ampt infiltration equation. A soil map of the Little Washita watershed in south-west Oklahoma, USA was used to investigate the effects of grid cell resolution on the distributed modelling of infiltration. Soil-derived parameters and infiltration exhibit decreased entropy as resolutions become coarser. This is reflected in a decrease in the maximum entropy value for the reclassified/derived parameters vis a vis the original data. Moreover, the entropy curve, when plotted against resolution, shows two distinct segments: a constant section where no entropy was lost with decreasing resolution and another part which is characterized by a sharp decrease in entropy after a critical resolution of 1209 m is reached. This methodology offers a technique for assessing the largest cell size that captures the spatial variability of infiltration parameters for a particular basin. A geographical information system (GIS) based rainfall-runoff model is used to simulate storm hydrographs using infiltration parameter maps at different resolutions as inputs. Model results up to the critical resolution are reproducible and errors are small. However, at resolutions beyond the critical resolution the results are erratic with large errors. A major finding of this study is that a large resolution (1209 m for this basin) yields reproducible model results. When modelling a river basin using a distributed model, the resolution (grid cell size) can drastically affect the model results and calibration. The error structure attributable to grid cell resolution using entropy as a spatial variability measure is shown.     

球面等积网格系统生成算法与软件模型研究

利用全球离散网格系统管理海量空间信息是一种受到普遍关注的新途径,具有等积特性的网格系统在大范围遥感影像数据处理方面具有重要的理论意义和应用前景。文中论述的算法在二十面体展开图上建立坐标系描述不同剖分产生的多分辨率网格,利用施奈德多面体投影映射到球面。针对算法特点设计具有网络模拟功能的软件模型。实验结果表明,该软件模型具有较好的扩展性。