Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition   总被引：3，自引：0，他引：3  

Oliver Baur  Gerrit Austen  Jürgen Kusche 《Journal of Geodesy》2008,82(4-5):207-221

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.  相似文献   

A Step Towards Evaluation of the Seismic Response Reduction Factor in Multistorey Reinforced Concrete Frames   总被引：1，自引：0，他引：1  

Barakat  Samar a.  Husein Malkawi  Abdallah I.  Al-Shatnawi  Anis S. 《Natural Hazards》1997,16(1):65-80

A seismic nonlinear time-history analysis was made for four-, six-, and eight-storey reinforced concrete buildings. These buildings were made as three-dimensional space frame structures with shear walls in both orthogonal directions. They have five bays with 4.8 m spacing each in the horizontal direction, and three bays with 4.2 m spacing each in the transversal direction. The frames were designed according to the Jordanian Seismic Code of practice for Seismic Zones 4, 3, 2, and 1 as proposed for Jordan by several authors. Time-history analysis was made using the El Centro (N-S) earthquake record of May 1940 as an actual earthquake excitation. The response reduction factor (R) that primarily consists of two factors that are the ductility reduction (Rµ) and the overstrength (), is obtained. It has been seen that the seismic zoning has a slight effect on the ductility reduction factor for different buildings, since it ranges from Zone 4 to Zone 1 as 2.37 to 2.52, 1.72 to 1.78, and 1.14 to 1.18 for four-, six-, and eight-storey buildings, respectively. Moreover, it is observed that, for different buildings and different seismic zones, the ductility reduction factor (Rµ) is slightly different from the system ductility factor (µ) especially for higher values of µ (i.e., Rµ µ). The response reduction factor, called overstrength (), was evaluated. The overstrength factor was found to vary with seismic zones (Z) , number of stories, and design gravity loads. However, the dependency on seismic zones was the strongest. The average overstrength of these buildings in Zones 4 and 1 was 2.61 and 6.94, respectively. The overstrength increased as the number of storeys decreased: overstrength of a four-storey building was higher than an eight-storey building by 36% in Zone 4, and 39% in Zone 1. Furthermore, buildings of the three heights had an average overstrength 165.9% higher in Zone 1 than in Zone 4. These observations have a significant implications for the seismic design codes which currently do not take into account the variation of the response reduction factor, R (i.e., ductility reduction factor times overstrength).  相似文献   

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

杨彦军  杨宇  张文喜 《矿物岩石》2003,23(2):108-110

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

生态地理建模中的多尺度问题   总被引：29，自引：0，他引：29  

岳天祥  刘纪远 《第四纪研究》2003,23(3):256-261

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

一种在异构系统中实现负载平衡的方法   总被引：2，自引：3，他引：2       下载免费PDF全文

金之雁  王鼎兴 《应用气象学报》2003,14(4):410-418

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

往复荷载作用下矩形钢管混凝土压弯构件承载力计算     

杨有福  韩林海 《地震工程与工程振动》2003,23(5):125-129

本文基于往复荷载作用下矩形钢管混凝土压弯构件的实验结果，比较了国内外比较典型的设计规范(程)，包括英国BS5400(1979)、美国ACI318—99(1999)和AISC-LRFD(1999)、日本AIJ(1997)、欧洲EC4(1994)、中国GJB4142-2000(2001)和福建省地方标准(送审稿)，及本文数值方法在计算往复荷载作用下矩形钢管混凝土压弯构件承载力的差异。结果表明，在进行往复荷载作用下矩形钢管混凝土压弯构件承载力计算时，各种计算方法获得的承载力都偏于安全，其中，数值计算结果与实验结果最接近，GJB4142—2000(2001)的计算结果与实验结果吻合程度稍差，ACI318—99(1999)，EC4(1994)和福建省地方标准(送审稿)的计算结果比实验结果约低20％，而BS5400(1979)、AISC—LRFD(1999)和AIJ(1997)的计算结果比实验结果总体上低30％以上。本文结果可供进行矩形钢管混凝土结构设计时参考。  相似文献   

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

胡振东  达汉桥  夏林元 《武汉大学学报(信息科学版)》1997,(3)

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

自适应空间信息移动服务   总被引：4，自引：0，他引：4  

张海堂罗睿  张永生李波 《测绘学院学报》2004,21(1):34-37,41

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

Clustering and visualization of earthquake data in a grid environment     

Dave A.Yuen  Benjamin J. Kadlec  Evan F. Bollig  Witold Dzwinel  Zachary A. Garbow  Cesar R. S. da Silva 《Visual Geosciences》2005,10(1):1-12

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     

Nathan Shauer  Carlos Armando Duarte 《国际地质力学数值与分析法杂志》2019,43(18):2707-2742

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.  相似文献