Parallelization of the flow-path network model using a particle-set strategy

High-performance simulation of flow dynamics remains a major challenge in the use of physical-based, fully distributed hydrologic models. Parallel computing has been widely used to overcome efficiency limitation by partitioning a basin into sub-basins and executing calculations among multiple processors. However, existing partition-based parallelization strategies are still hampered by the dependency between inter-connected sub-basins. This study proposed a particle-set strategy to parallelize the flow-path network (FPN) model for achieving higher performance in the simulation of flow dynamics. The FPN model replaced the hydrological calculations on sub-basins with the movements of water packages along the upstream and downstream flow paths. Unlike previous partition-based task decomposition approaches, the proposed particle-set strategy decomposes the computational workload by randomly allocating runoff particles to concurrent computing processors. Simulation experiments of the flow routing process were undertaken to validate the developed particle-set FPN model. The outcomes of hourly outlet discharges were compared with field gauged records, and up to 128 computing processors were tested to explore its speedup capability in parallel computing. The experimental results showed that the proposed framework can achieve similar prediction accuracy and parallel efficiency to that of the Triangulated Irregular Network (TIN)-based Real-Time Integrated Basin Simulator (tRIBS).     

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

一种求气井二项式方程系数的新方法

针对开发早期气田缺乏产能试井资料的实际情况，由于无法用传统的压力降法计算产能方程，故而基于最优化算法，提出了一种利用稳定生产数据和井控储量求取气井产能方程系数的新方法。在只知道稳定生产数据情况下便可进行计算分析。该方法在气田进行了实际应用，并和实际资料以及传统试凑法计算的结果进行了比较，表明方法的计算结果较为可靠。     

生态地理建模中的多尺度问题

本文在分析生态地理建模内涵的基础上,讨论了生态地理建模中的尺度转换问题、跨尺度相互作用问题、空间尺度与时间尺度的关联问题和多尺度数据处理问题.由于生态地理问题的非线性、生态环境的异质性和随机事件,简单的线性尺度转换方法远不能满足生态地理建模的要求.为了从根本上解决生态地理建模中的时空尺度问题,除需要运用微分几何学和等级理论等经典方法外,还需要引入格点生成法和网格计算等现代理论和技术手段.     

一种在异构系统中实现负载平衡的方法

提出了在异构系统实现负载平衡的区域分解算法和实现负载平衡的计算方法，利用它的负反馈性质解决了异构系统处理机计算速度测量误差造成的负载测量不准问题，并对处理机速度变化，速度测量误差、处理机数量、网格点计算量的分布等因素的影响进行了计算，结果表明本方法具有很强的平衡负载能力和较强的适应性；根据计算结果提出了解决模式网格点计算量不易测量问题的解决方案，并用扩散方程和模拟物理过程进行试验，试验表明这种方法是可行的，平衡负载的效果十分显著。     

计算机网络的分布式计算环境及分布式处理

本文对网络计算模式、分布式计算环境以及分布式处理机能作了分析和探讨。     

自适应空间信息移动服务

空间信息移动服务系统是空间信息系统继由传统PC计算环境向有线Web分布式计算环境扩展后，向移动计算环境的新发展。文中在分析移动计算环境概念和特点的基础上，讨论了空间信息移动服务的若干特点及关键技术，研究了自适应空间信息移动服务方案。     

Clustering and visualization of earthquake data in a grid environment

We present a web client-server service WEB-IS, which we have developed for remote analysis and visualization of seismic data consisting of both small magnitude events and large earthquakes. We show that the problem-solving environment (PSE) intended for prediction of large magnitude earthquakes can be based on this WEB-IS idea. The clustering schemes, feature generation, feature extraction techniques and rendering algorithms form a computational framework of this environment. On the other hand, easy and fast access both to the seismic data distributed among distant computing resources and to computational and visualization resources can be realized in a GRID framework. We discuss the usefulness of NaradaBrokering (iNtegrated Asynchronous Real-time Adaptive Distributed Architecture) as a middleware, allowing for flexibility and high throughput for remote visualization of geophysical data. The WEB-IS functionality was tested both on synthetic and the actual earthquake catalogs. We consider the application of similar methodology for tsunami alerts.     

Improved algorithms for generalized finite element simulations of three-dimensional hydraulic fracture propagation

This paper reports improvements to algorithms for the simulation of 3-D hydraulic fracturing with the Generalized Finite Element Method (GFEM). Three optimizations are presented and analyzed. First, an improved initial guess based on solving a 3-D elastic problem with the pressure from the previous step is shown to decrease the number of Newton iterations and increase robustness. Second, an improved methodology to find the time step that leads to fracture propagation is proposed and shown to decrease significantly the number of iterations. Third, reduced computational cost is observed by properly recycling the linear part of the coupled stiffness matrix. Two representative examples are used to analyze these improvements. Additionally, a methodology to include the leak-off term is presented and verified against asymptotic analytical solutions. Conservation of mass is shown to be well satisfied in all examples.     

Boussinesq相位解析的海岸水动力学数学模型研究进展

由于在平衡计算效率和精度上具有优势,Boussinesq相位解析数学模型研究不断取得突破,已成为波浪和水流精细化模拟的较优解析方式,为海岸工程、环境、地质等问题提供了实用和高效的研究手段。本文对已有Boussinesq类模型的研究进行了评述,深入探讨其重要发展、实际应用和理论瓶颈,从高阶非静压修正、GPU准三维高性能算法编译、波浪破碎和泥沙运移沉积等4个方面提出未来可能的科学突破方向。