Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method

Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd.     

Analysis of the contributions of topographic,soil, and vegetation features on the spatial distributions of surface soil moisture in a steep natural forested headwater catchment

Surface soil moisture has been extensively studied for various land uses and landforms. Although many studies have reported potential factors that control surface soil moisture over space or time, the findings have not always been consistent, indicating a need for identification of the main factors. This study focused on the static controls of topographic, soil, and vegetation features on surface soil moisture in a steep natural forested headwater catchment consisting of three hillslope units of a gully area, side slope, and valley‐head slope. Using a simple correlation analysis to investigate the effects of the static factors on surface soil moisture at depths of 0–20 cm at 470 points in 13 surveys, we addressed the characteristics of surface soil moisture and its main controlling factors. The results indicated that the mean of surface soil moisture was in the decreasing order of gully area > valley‐head slope > side slope. The relationship between the mean and standard deviation of surface soil moisture showed a convex‐upward shape in the headwater catchment, a negative curvilinear shape in the gully area, and positive curvilinear shapes at the side and valley‐head slopes. At the headwater catchment and valley‐head slope, positive contributions of soil porosity and negative contributions of slope gradient and saturated hydraulic conductivity were the main controlling factors of surface soil moisture under wetter conditions, whereas positive contributions of topographic wetness index and negative contributions of vegetation density were the main controlling factors of surface soil moisture under drier conditions. At the side slope underlain by fractured bedrocks, only saturated hydraulic conductivity and vegetation density were observed to be the controlling factors. Surface soil moisture in the gully area was mainly affected by runoff rather than were static features. Thus, using hillslope units is effective for approximately estimating the hydrological behaviours of surface moisture on a larger scale, whereas dependency between the main static factors and moisture conditions is helpful for estimating the spatial distributions of surface moisture on a smaller scale.     

陕西洛川不同时代黄土中的碳酸钙的沉积速率及其在地质年代测定中的意义

选择地处西北半干旱区、又具有代表性的黑木沟标准剖面,逐一测试其土层中的碳酸钙与干容重,进而由相关参数推算各土层的碳酸钙的沉积速率(V,下同),这是笔者的尝试。结果表明,该剖面中的碳酸钙的沉积速率,均有自下而上递增的趋势,各层总平均速率为0.20g/cm2·103a     

青藏高原与落基山背风气旋的对比分析(一)——天气结构

本文分别选用了一个东亚(EA)和北美(NA)的背风气旋进行对比分析.分析了两者的天气形势、动力和热力结构的差异.由于青藏高原和落基山地形对气流的不同影响,北美水汽从墨西哥湾向北输送,东亚水汽从南海向东北输送,因而,NA的对流不稳定区成南北带状分布;而EA的则成东西带状分布.EA的初生时,正压不稳定作用是主要的;在发展中,起主导作用的是斜压不稳定.     

青藏高原与落基山背风气旋的对比分析(二)——诊断分析

本文对第一部分所选的东亚和北美的背风气旋个例继续进行诊断分析。利用Q矢量诊断了两者的锋生函数和垂直运动,同时对两者进行了涡度收支、热量和水汽收支的计算,对比分析了两者物理量的分布特征。结果表明:北美的锋生与落基山地形影响有关,在落基山东侧,对流活跃;而东亚的锋生与海洋上感热和水汽蒸发有关,气旋入海后,对流作用加强。     

青藏高原和落基山对气旋的动力影响

本文利用北京气象中心的业务预报模式(B模式),采用全风速调整方案,进行了绕流、爬坡等各种数值试验,研究青藏高原和落基山对气旋的动力影响。结果表明,气流绕青藏高原而行,产生的散度场与东亚的气旋生成区和路径匹配;在青藏高原主体下游,高层辐合,抑制了我国大陆东部的气旋发展。在北美,落基山对大气的动力影响以爬坡为主,由于气流在迎风坡上升,空气绝热冷却,在背风面大规模下沉增温,对流层低层形成暖脊,这对背风锋生和背风气旋的发生发展有着直接的影响;另外爬坡作用也是美国中西部干线形成的主要原因。     

钨、锡流-熔分配实验结果及其矿床成因意义

本文用实验确定了钨、锡在成分不同的花岗质熔体相与共存水热流体相的流-熔分配系数(D_(Me)~(V/L))及其与介质溶液(NaF,KF,HF,NaCl等水溶液)摩尔浓度间的函数关系。实验结果表明,钨、锡的分配行为明显不同。在相同条件下,随体系的不同,D_W~(V/L)比D_(Sn)~(V/L)大几倍至二十倍。钠和钾对钨、锡的流-熔分配行为的影响基本相同,而氟和氯对钨、锡的分配行为的影响相差甚远,花岗质熔体的主成分对D_W~(V/L)和D_(Sn)~(V/L)有复杂的影响。利用这些结果探讨与花岗岩有关的钨、锡矿床的成矿机理,得出了一些与前人不同的新认识。     

冻土的蠕变及蠕变强度

基于试验，本本讨论分析了冻土三轴蠕变规律及温度和围压对冻土蠕变强度的影响，给出了土蠕变及蠕变强度随时间降低的方程式，进而提出了冻土蠕变强度的物型屈服准则。通过改变参数可将冻土的瞬时强度准则与冻土的瞬时强度准则与冻土的蠕变强度准则统一为同一个方程描述。     

山东沂水太古宙麻粒岩区年代学研究成果简介

本文以详细研究本区麻粒岩-斜长角闪岩等表壳岩的岩石学、地球化学和 Sm-Nd 同位素定年为重点。所取得的成果，对华北陆台古老地层的分布和古构造格局的探索提供了实际资料。