Low-order stabilized finite element for the full Biot formulation in soil mechanics at finite strain     

Lluís Monforte  Pedro Navas  Josep Maria Carbonell  Marcos Arroyo  Antonio Gens 《国际地质力学数值与分析法杂志》2019,43(7):1488-1515

This article presents a novel finite element formulation for the Biot equation using low-order elements. Additionally, an extra degree of freedom is introduced to treat the volumetric locking steaming from the effective response of the medium; its balance equation is also stabilized. The accuracy of the proposed formulation is demonstrated by means of numerical analyses.  相似文献   

Impacts of groundwater flow on the evolution of a thermokarst lake in the permafrost–dominated region on the Qinghai–Tibet plateau     

Jinlong Li  Wei Wang  Xianmin Ke  Jianhua Li  Yating Yu  Zeyong Gao  Fujun Niu 《水文研究》2021,35(12):e14443

As a result of global warming induced permafrost degradation in recent decades, thermokarst lakes in the Qinghai–Tibet plateau (QTP) have been regulating local hydrological and ecological processes. Simulations with coupled moisture–heat numerical models in the Beiluhe basin (located in the hinterland of permafrost regions on the QTP) have provided insights into the interaction between groundwater flow and the freeze–thaw process. A total of 30 modified SUTRA scenarios were established to examine the effects of hydrodynamic forces, permeability, and climate on thermokarst lakes. The results indicate that the hydrodynamic condition variables regulate the permafrost degradation around the lakes. In case groundwater recharges to the lake, a low–temperature groundwater flow stimulates the expansion of the surrounding thawing regions through thermal convection. The thawing rate of the permafrost underlying the lake intensifies when groundwater is discharged from the lake. Under different permeability conditions, spatiotemporal variations in the active layer thickness significantly influence the occurrence of an open talik at the lake bottom. A warmer and wetter climate will inevitably lead to a sharp decrease in the upper limit of the surrounding permafrost, with a continual decrease in the duration of open talik events. Overall, our results underscore that comprehensive consideration of the relevant hydrologic processes is critical for improving the understanding of environmental and ecological changes in cold environments.  相似文献   

An Eocene/Oligocene blueschist‐/greenschist facies P–T loop from the Cycladic Blueschist Unit on Naxos Island,Greece: Deformation‐related re‐equilibration vs. thermal relaxation          下载免费PDF全文

Alexandre Peillod  Uwe Ring  Johannes Glodny  Alasdair Skelton 《Journal of Metamorphic Geology》2017,35(7):805-830

Geothermobarometric and geochronological work indicates a complete Eocene/early Oligocene blueschist/greenschist facies metamorphic cycle of the Cycladic Blueschist Unit on Naxos Island in the Aegean Sea region. Using the average pressure–temperature (P–T) method of thermocalc coupled with detailed textural work, we separate an early blueschist facies event at 576 ± 16 to 619 ± 32°C and 15.5 ± 0.5 to 16.3 ± 0.9 kbar from a subsequent greenschist facies overprint at 384 ± 30°C and 3.8 ± 1.1 kbar. Multi‐mineral Rb–Sr isochron dating yields crystallization ages for near peak‐pressure blueschist facies assemblages between 40.5 ± 1.0 and 38.3 ± 0.5 Ma. The greenschist facies overprint commonly did not result in complete resetting of age signatures. Maximum ages for the end of greenschist facies reworking, obtained from disequilibrium patterns, cluster near c. 32 Ma, with one sample showing rejuvenation at c. 27 Ma. We conclude that the high‐P rocks from south Naxos were exhumed to upper mid‐crustal levels in the late Eocene and early Oligocene at rates of 7.4 ± 4.6 km/Ma, completing a full blueschist‐/greenschist facies metamorphic cycle soon after subduction within c. 8 Ma. The greenschist facies overprint of the blueschist facies rocks from south Naxos resulted from rapid exhumation and associated deformation/fluid‐controlled metamorphic re‐equilibration, and is unrelated to the strong high‐T metamorphism associated with the Miocene formation of the Naxos migmatite dome. It follows that the Miocene thermal overprint had no impact on rock textures or Sr isotopic signatures, and that the rocks of south Naxos underwent three metamorphic events, one more than hitherto envisaged.  相似文献   

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation          下载免费PDF全文

Chongbin Zhao  Thomas Poulet  Klaus Regenauer‐Lieb 《国际地质力学数值与分析法杂志》2015,39(6):571-593

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method          下载免费PDF全文

Mohammad Hossein Bazyar  Abbas Talebi 《国际地质力学数值与分析法杂志》2015,39(1):1-22

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Assessment of flow paths and confluences for saltwater intrusion in a deltaic river network          下载免费PDF全文

Baoshan Cui  Xiaojing Shao  Zhiming Zhang 《水文研究》2015,29(20):4549-4558

Saltwater intrusion is a serious issue in estuarine deltas all over the world due to rapid urban sprawl and water shortage. Therefore, detecting the major flow paths or locations at risk of saltwater intrusion in estuarine ecosystems is important for mitigating saltwater intrusion. In this paper, we introduce a centrality index, the betweenness centrality (BC), to address this problem. Using the BC as the weighted attribute of the river network, we identify the critical confluences for saltwater intrusion and detect the preferential flow paths for saltwater intrusion through the least‐cost‐path algorithm from a graph theory approach. Moreover, we analyse the responses of the BC values of confluences calculated in the river network to salinity. Our results show that the major flow paths and critical confluences for saltwater intrusion in a deltaic river network can be represented by the least cost paths and the BC values of confluences, respectively. In addition, a significant positive correlation between the BC values of confluences and salinity is determined in the Pearl River Delta. Changes in the salinity can produce significant variation in the BC values of confluences. Therefore, freshwater can be diverted into these major flow paths and critical confluences to improve river network management under saltwater intrusion. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

利用跨断层形变资料研究天山中段断层现今活动特征     

韩桂红  滕海涛  高歌  李桂荣  张治广 《大地测量与地球动力学》2019,39(1):72-76

利用2001年以来天山中段多个跨断层形变场地的测量资料，从断层运动方向和断层平均速率2个方面分析天山中段3条断层现今活动特征。结果表明，霍玛吐断裂呼图壁段断层的活动性质为正、逆断和左、右旋交替活动；北轮台-辛格尔断裂库尔勒段断层表现为右旋运动特征，活动性质为逆断性；妖魔山F2断裂乌鲁木齐段的断层活动性质以逆断为主。妖魔山F2断裂乌鲁木齐段断层的活动程度最高，霍玛吐断裂呼图壁段最低。霍玛吐断裂、北轮台-辛格尔断裂、妖魔山F2断裂垂直方向上基本均处于张性变化，前两者水平向处于挤压状态；断层活动速率的时间进程具有阶段性特征，各个时间段转折点有较好的匹配性。  相似文献   

Fingerprinting the 8.2 ka event climate response in a coupled climate model     

Ane P. Wiersma  Didier M. Roche  Hans Renssen 《第四纪科学杂志》2011,26(1):118-127

Using results from coupled climate model simulations of the 8.2 ka climate event that produced a cold period over Greenland in agreement with the reconstructed cooling from ice cores, we investigate the typical pattern of climate anomalies (fingerprint) to provide a framework for the interpretation of global proxy data for the 8.2 ka climate event. For this purpose we developed an analysis method that isolates the forced temperature response and provides information on spatial variations in magnitude, timing and duration that characterise the detectable climate event in proxy archives. Our analysis shows that delays in the temperature response to the freshwater forcing are present, mostly in the order of decades (30 a over central Greenland). The North Atlantic Ocean initially cools in response to the freshwater perturbation, followed in certain parts by a warm response. This delay, occurring more than 200 a after the freshwater pulse, hints at an overshoot in the recovery from the freshwater perturbation. The South Atlantic and the Southern Ocean show a warm response reflecting the bipolar seesaw effect. The duration of the simulated event varies for different areas, and the highest probability of recording the event in proxy archives is in the North Atlantic Ocean area north of 40° N. Our results may facilitate the interpretation of proxy archives recording the 8.2 ka event, as they show that timing and duration cannot be assumed to correspond with the timing and duration of the event as recorded in Greenland ice cores. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

The Scientific Conceptual Framework for Ecological Quality of the Dryland Ecosystem: Concepts,Indicators, Monitoring and Assessment     

WU Rina  CONG Weiwei  LI Yonghua  LI Siyao  WANG Dongfang  JIA Zhiqing  WANG Feng 《资源与生态学报(英文版)》2019,10(2):196-201

The dryland ecosystem is the dominant component of the global terrestrial ecosystem since arid regions occupy 45% of the earth’s land area and feed 38% of the world's population. The stability and sustainable development of the dryland ecosystem are critical for achieving the millennium development goal (MDG) in the arid and semiarid areas. However, there is still no scientific guideline for measuring and conserving the health and productivity of dryland ecosystems. Therefore, the purpose of this study is to develop the scientific conceptual framework of defining, monitoring and evaluating the ecological quality of dryland ecosystems. The ecological quality of dryland ecosystems is represented by a system of comprehensive indicators that are each extracted from the ecological elements, and structural and functional indices of the ecosystem. These indicators can be monitored by integrating satellites and unmanned aerial vehicles with ground-based sensor networks at the scale of either observational sites or regional scales. Finally, the ecological quality can be evaluated by evaluation models based on the normalized index values and their thresholds. This paper presents a conceptual framework for monitoring and evaluating the ecological quality of drylands, which provides a way of advancing the monitoring, diagnosis, and evaluation of the ecological quality of the dryland ecosystems.  相似文献   

青藏高原东南缘构造应力场重构     

侯强  丁小军  赵宏  欧明霖  张博康 《大地测量与地球动力学》2019,39(5):458-463

采用有限元方法，针对青藏高原东南缘建立更细致、更精确的三维有限元弹性模型。选取9种不同的应力边界条件，分别进行优化分析后处理，将对应台站形变模拟值与GPS实测值进行误差分析，最终选取最佳方案作为古构造应力场。结果表明，青藏高原东南缘4 Ma BP的古应力场主要起源于中国大陆周围板块的相互作用，特别是印度板块NNE向强烈碰撞作用，成为中国大陆尤其是西南部青藏高原地区构造应力场最主要的动力来源，控制各个块体相互作用的方式和运动格局。青藏高原东南缘古应力场主要包括几个力源：西北部青藏高原侧向挤压造成的WE向应力约105 MPa；西南部直接来自于印度板块的NE-WS向应力约70 MPa；南部NS向作用力33 MPa；东南部扬子块体侧向NW-SE阻挡力56 MPa；北东部受扬子块体强烈EW向阻挡力90 MPa。这些力源共同作用于青藏高原东南缘，形成现今复杂应力场。
  相似文献