Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA     

Quansheng Liu  Weiqin Wang  Hao Ma 《国际地质力学数值与分析法杂志》2020,44(2):208-238

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.  相似文献   

Rainfall–runoff responses on Arctic hillslopes underlain by continuous permafrost,North Slope,Alaska, USA          下载免费PDF全文

Caitlin R. Rushlow  Sarah E. Godsey 《水文研究》2017,31(23):4092-4106

The Arctic hydrologic cycle is intensifying, as evidenced by increased rates of precipitation, evapotranspiration, and riverine discharge. However, the controls on water fluxes from terrestrial to aquatic systems in upland Arctic landscapes are poorly understood. Upland landscapes account for one third of the Arctic land surface and are often drained by zero‐order geomorphic flowpath features called water tracks. Previous work in the region attributed rapid runoff response at larger stream orders to water tracks, but models suggest water tracks are hydrologically disconnected from the surrounding hillslope. To better understand the role of water tracks in upland landscapes, we investigated the surface and subsurface hydrologic responses of 6 water tracks and their hillslope watersheds to natural patterns of rainfall, soil thaw, and drainage. Between storms, both water track discharge and the water table in the hillslope watersheds exhibited diel fluctuations that, when lagged by 5 hr, were temporally correlated with peak evapotranspiration rate. Water track soils remained saturated for more of the summer season than soils in their surrounding hillslope watersheds. When rainfall occurred, the subsurface response was nearly instantaneous, but the water tracks took significantly longer than the hillslopes to respond to rainfall, and longer than the responses previously observed in nearby larger order Arctic streams. There was also evidence for antecedent soil water storage conditions controlling the magnitude of runoff response. Based on these observations, we used a broken stick model to test the hypothesis that runoff production in response to individual storms was primarily controlled by rainfall amount and antecedent water storage conditions near the water track outlet. We found that the relative importance of the two factors varied by site, and that water tracks with similar watershed geometries and at similar landscape positions had similar rainfall–runoff model relationships. Thus, the response of terrestrial water fluxes in the upland Arctic to climate change depends on the non‐linear interactions between rainfall patterns and subsurface water storage capacity on hillslopes. Predicting these interactions across the landscape remains an important challenge.  相似文献   

北斗系统周跳探测与修复方法探究     

李彬 《测绘与空间地理信息》2016,(8):160-165

北斗卫星导航系统是中国自行研制的继美国GPS和俄罗斯GLONASS之后第3个成熟的卫星导航系统。进行GNSS相位数据处理的过程中,由于某种原因,载波相位观测值突然出现整周数的跳变的现象称为周跳。周跳现象的发生将会严重影响后续的数据处理,因此,有必要在解算前完成周跳探测与修复工作。本文就北斗系统周跳探测与修复课题进行了较为全面的探究,先后介绍了高次差法、多普勒法、伪距载波相位组合法、电离层残差法、三频数据组合探测法等一系列周跳探测与修复方法,可分别应用于不同频段接收机、不同观测环境、不同数据处理模式。  相似文献   

L1范数与IQR统计量组合的GNSS坐标序列粗差探测算法     

明锋  曾安敏  景一帆 《测绘科学技术学报》2016,(2):127-132

GNSS坐标时间序列中不可避免地含有粗差,未剔除的粗差将会导致参数估计有偏。因此,粗差探测与剔除是GNSS坐标序列分析中一项重要的数据预处理工作。针对GNSS坐标时间序列特点,提出了一种将L1范数(L1-norm)估计与四分位距统计量IQR(interquartile range)组合的移动开窗粗差探测算法,称之为L1_Mod IQR。该方法的主要思想是,首先利用L1范数估计得到较"真实"的残差,然后再对残差采用IQR统计量进行粗差探测。将L1_Mod IQR法与"3σ"法、基于最小二乘的τ检验法等粗差探测算法进行了模拟计算与对比,验证了该算法的有效性。进一步采用L1_Mod IQR算法对中国区域10个IGS站的高程时间序列进行了分析,结果表明中国区域IGS站高程序列的粗差剔除率最小为0.1%,最大为2.6%。并且以WUHN站为例与SOPAC提供的结果进行了对比,结果表明SOPAC提供的"Clean"数据仍含有大量的粗差,而L1_Mod IQR算法能够有效地剔除粗差。  相似文献   

Low-temperature constraints on the Alpine thermal evolution of the Western Carpathian basement rock complexes     

《Journal of Structural Geology》2016

The Cretaceous to Palaeogene Alpine exhumation of previously buried Variscan basements is recorded in the southern portion of the Western Carpathians in the Gemeric and Veporic units. The Meso-Cenozoic collisional processes resulted in basement exhumation of the Tatric Unit from Palaeogene to Neogene times. According to zircon and apatite fission track data, the Gemeric Unit, an uppermost thick-skinned thrust sheet, cooled from depth levels of ∼10 up to 2.5 km (temperature interval of ∼250–60 °C) about 88–64 Ma ago, after the collapse of overlying Meliata-Turňa-Silica Mesozoic accretionary prism. The middle and lower thick-skinned thrust sheets, Veporic and Tatric units, cooled from the depths of ∼10 up to 2.5 km ∼110–40 Ma ago. The process was controlled by unroofing of footwall from beneath the Gemeric Unit. About 50–20 Ma ago, the internal zone of Tatric Unit gradually exhumed to depth of <2 km and some parts of the unit appeared at the surface level. However, the external zone of Tatric Unit was buried beneath the Eocene to Lower Miocene sedimentary successions and exhumed to the subsurface level at ∼21–8 Ma ago, as a result of oblique collision of the Western Carpathians with the European Platform.  相似文献   

Exhumation history of the Peake and Denison Inliers: insights from low-temperature thermochronology     

J. W. Hall  S. Glorie  A. S. Collins  A. Reid  N. Evans  B. McInnes 《Australian Journal of Earth Sciences》2016,63(7):805-820

Multi-method thermochronology applied to the Peake and Denison Inliers (northern South Australia) reveals multiple low-temperature thermal events. Apatite fission track (AFT) data suggest two main time periods of basement cooling and/or reheating into AFT closure temperatures (～60–120°C); at ca 470–440 Ma and ca 340–300 Ma. We interpret the Ordovician pulse of rapid basement cooling as a result of post-orogenic cooling after the Delamerian Orogeny, followed by deformation related to the start of the Alice Springs Orogeny and orocline formation relating to the Benambran Orogeny. This is supported by a titanite U/Pb age of 479 ± 7 Ma. Our thermal history models indicate that subsequent denudation and sedimentary burial during the Devonian brought the basement rocks back to zircon U–Th–Sm/He (ZHe) closure temperatures (～200–150°C). This period was followed by a renewal of rapid cooling during the Carboniferous, likely as the result of the final pulses of the Alice Springs Orogeny, which exhumed the inlier to ambient surface temperatures. This thermal event is supported by the presence of the Mount Margaret erosion surface, which indicates that the inlier was exposed at the surface during the early Permian. During the Late Triassic–Early Jurassic, the inlier was subjected to minor reheating to AFT closure temperatures; however, the exact timing cannot be deduced from our dataset. Cretaceous apatite U–Th–Sm/He (AHe) ages coupled with the presence of contemporaneous coarse-grained terrigenous rocks suggest a temporally thermal perturbation related with shallow burial during this time, before late Cretaceous exhumation cooled the inliers back to ambient surface temperatures.  相似文献   

Possible sources of forecast errors generated by the global/regional assimilation and prediction system for landfalling tropical cyclones. Part I: Initial uncertainties     

Feifan Zhou  Munehiko Yamaguchi  Xiaohao Qin 《大气科学进展》2016,33(7):841-851

This paper investigates the possible sources of errors associated with tropical cyclone (TC) tracks forecasted using the Global/Regional Assimilation and Prediction System (GRAPES). The GRAPES forecasts were made for 16 landfalling TCs in the western North Pacific basin during the 2008 and 2009 seasons, with a forecast length of 72 hours, and using the default initial conditions (“initials”, hereafter), which are from the NCEP-FNL dataset, as well as ECMWF initials. The forecasts are compared with ECMWF forecasts. The results show that in most TCs, the GRAPES forecasts are improved when using the ECMWF initials compared with the default initials. Compared with the ECMWF initials, the default initials produce lower intensity TCs and a lower intensity subtropical high, but a higher intensity South Asia high and monsoon trough, as well as a higher temperature but lower specific humidity at the TC center. Replacement of the geopotential height and wind fields with the ECMWF initials in and around the TC center at the initial time was found to be the most efficient way to improve the forecasts. In addition, TCs that showed the greatest improvement in forecast accuracy usually had the largest initial uncertainties in TC intensity and were usually in the intensifying phase. The results demonstrate the importance of the initial intensity for TC track forecasts made using GRAPES, and indicate the model is better in describing the intensifying phase than the decaying phase of TCs. Finally, the limit of the improvement indicates that the model error associated with GRAPES forecasts may be the main cause of poor forecasts of landfalling TCs. Thus, further examinations of the model errors are required.  相似文献   

ZQZ TF型风传感器冻结故障自动检测装置设计     

朱世恒  何晓庆  朱平 《气象科技》2016,44(6):902-906

近年来我国大部分地区在冬、春两季多遭受冻雨灾害,自动气象站的风向风速传感器容易因冻结而无法正常工作。针对冻结故障的实时检测问题,本文设计了一种基于ZQZ-TF型风传感器的故障自动检测装置。该装置通过检测风传感器的工作电压和工作电流,实现对传感器实时状态的监控。同时,结合具体的故障诊断算法,能够快速判断风传感器是否被冻结。系统的核心硬件电路、软件算法可分别集成在自动气象站风向风速数据采集电路和嵌入式软件中,具有结构紧凑、操作简单的特点。经验证,该装置能够实时检测风传感器冻结故障,且工作稳定,检测精度高。  相似文献   

Retrieval of effective leaf area index (LAIe) and leaf area density (LAD) profile at individual tree level using high density multi-return airborne LiDAR     

《International Journal of Applied Earth Observation and Geoinformation》2016

As an important canopy structure indicator, leaf area index (LAI) proved to be of considerable implications for forest ecosystem and ecological studies, and efficient techniques for accurate LAI acquisitions have long been highlighted. Airborne light detection and ranging (LiDAR), often termed as airborne laser scanning (ALS), once was extensively investigated for this task but showed limited performance due to its low sampling density. Now, ALS systems exhibit more competing capacities such as high density and multi-return sampling, and hence, people began to ask the questions like—“can ALS now work better on the task of LAI prediction?” As a re-examination, this study investigated the feasibility of LAI retrievals at the individual tree level based on high density and multi-return ALS, by directly considering the vertical distributions of laser points lying within each tree crown instead of by proposing feature variables such as quantiles involving laser point distribution modes at the plot level. The examination was operated in the case of four tree species (i.e. Picea abies, Pinus sylvestris, Populus tremula and Quercus robur) in a mixed forest, with their LAI-related reference data collected by using static terrestrial laser scanning (TLS). In light of the differences between ALS- and TLS-based LAI characterizations, the methods of voxelization of 3D scattered laser points, effective LAI (LAIe) that does not distinguish branches from canopies and unified cumulative LAI (ucLAI) that is often used to characterize the vertical profiles of crown leaf area densities (LADs) was used; then, the relationships between the ALS- and TLS-derived LAIes were determined, and so did ucLAIs. Tests indicated that the tree-level LAIes for the four tree species can be estimated based on the used airborne LiDAR (R² = 0.07, 0.26, 0.43 and 0.21, respectively) and their ucLAIs can also be derived. Overall, this study has validated the usage of the contemporary high density multi-return airborne LiDARs for LAIe and LAD profile retrievals at the individual tree level, and the contribution are of high potential for advancing forest ecosystem modeling and ecological understanding.  相似文献   

利用GNSS干涉信号振荡幅度反演土壤湿度   总被引：2，自引：1，他引：1  

汉牟田  张波  杨东凯  洪学宝  杨磊  宋曙辉 《测绘学报》2016,45(11):1293-1300

根据干涉效应和GNSS接收机信噪比估计方法,推导了利用GNSS干涉信号幅度进行土壤湿度反演的模型,建模过程考虑了天线增益、土壤介电常数和噪声的影响。提出了使用AMPD算法从含有噪声的归一化干涉功率曲线中提取干涉峰值与谷值,进而反演了土壤介电常数与土壤湿度的方法,并对其进行了仿真。结果表明,利用提取出来的干涉谷值进行反演性能比峰值好,相对稳定准确的卫星仰角范围为5°~25°,湿度大于0.06cm~3/cm~3时反演结果更为准确,标准差在0.01cm~3/cm~3左右波动。  相似文献