A 3D discrete FEM iterative algorithm for solving the water pipe cooling problems of massive concrete structures

Water pipe cooling has been widely used for the temperature control and crack prevention of massive concrete structures such as high dams. Because both under‐cooling and over‐cooling may reduce the efficiency of crack prevention, or even lead to great harm to structures, we need an accurate and robust numerical tool for the prediction of cooling effect. Here, a 3D discrete FEM Iterative Algorithm is introduced, which can simulate the concrete temperature gradient near the pipes, as well as the water temperature rising along the pipes. On the basis of the heat balance between water and concrete, the whole temperature field of the problem can be computed exactly within a few iteration steps. Providing the pipe meshing tool for building the FE model, this algorithm can take account of the water pipe distribution, the variation of water flow, water temperature, and other factors, while the traditional equivalent algorithm based on semi‐theoretical solutions can only solve problems with constant water flow and water temperature. The validation and convergence are proved by comparing the simulated results and analytical solutions of two standard second‐stage cooling problems. Then, a practical concrete block with different cooling schemes is analyzed and the influences of cooling factors are investigated. In the end, detailed guidance for pipe system optimization is provided. Copyright © 2015 John Wiley & Sons, Ltd.     

Regional evapotranspiration rate of oasis and surrounding desert

The surface energy balance algorithm for land method was used in this study to calculate the evapotranspiration (ET) rate for the middle reaches of the Heihe River Basin, Gansu Province, China, to analyse ET distribution within the oasis and the surrounding desert and, especially, on the edge zone of the oasis. Five profile graphs were created vertical to the river. Because of the inverse humidity phenomenon, the least amount of evapotranspiration occurred on the desert close to the oasis. The average evapotranspiration rate was roughly proportioned from the edge of the oasis to inside and outside its boundary. Two meteorological ground stations located close to the oasis edge showed a notable difference in net radiation flux that led to the difference found in ET. The primary reason for the significant differences observed in net radiation may be largely the result of differences in flux reflectivity and surface temperature. Meteorological data show that water supply also played an important role. Copyright © 2012 John Wiley & Sons, Ltd.     

Particle manifold method (PMM): A new continuum‐discontinuum numerical model for geomechanics

Particle manifold method (PMM) is a new extension of the numerical manifold method (NMM). PMM uses a mathematical cover system to describe the motion and deformation of a particle‐based physical domain. By introducing the concept of particle into NMM, PMM takes the advantages of easy topological and contact operations with particles. In this article, the methodology, formulations and implementation of the method are presented, together with modelling examples for validation. It is found that good solutions for both continuous and discontinuous problems are obtained by the new developed PMM. Due to the underlying coupled continuum‐discontinuum property of PMM, it has great potential for modelling of geomechanical problems. Copyright © 2012 John Wiley & Sons, Ltd.     

Parallelization of the distinct lattice spring model

The distinct lattice spring model (DLSM) is a newly developed numerical tool for modeling rock dynamics problems, i.e. dynamic failure and wave propagation. In this paper, parallelization of DLSM is presented. With the development of parallel computing technologies in both hardware and software, parallelization of a code is becoming easier than before. There are many available choices now. In this paper, Open Multi‐Processing (OpenMP) with multicore personal computer (PC) and message passing interface (MPI) with cluster are selected as the environments to parallelize DLSM. Performances of these parallel DLSM codes are tested on different computers. It is found that the parallel DLSM code with OpenMP can reach a maximum speed‐up of 4.68× on a quad‐core PC. The parallel DLSM code with MPI can achieve a speed‐up of 40.886× when 256 CPUs are used on a cluster. At the end of this paper, a high‐resolution model with four million particles, which is too big to handle by the serial code, is simulated by using the parallel DLSM code on a cluster. It is concluded that the parallelization of DLSM is successful. Copyright © 2011 John Wiley & Sons, Ltd.     

Mosaic desert pavement influences water infiltration and vegetation distribution on fluvial fan surfaces

Desert pavements (DPs) are critical for maintaining ecological stability and promoting near-surface hydrological cycling in arid regions. However, few studies have focused on eco-hydrological processes of DPs in the ecological systems of fluvial fans. Although DP surfaces appear to be barren and flat, we found that the surfaces are characterized by surface mosaic patterns of desert pavement (mosaic DP) and bare ground (mosaic BG). We investigated the effects of mosaic DP on water infiltration and vegetation distribution at six sites in fluvial fans (one on a hillside and five within the sectors of fans) along a southwest belt transect in northern Linze County, in the central Hexi Corridor (China). We found significant differences in mosaic DP between the hillside and sector sites in terms of pavement thickness and vesicular horizon thickness (Av thickness), particle composition, and bulk density, although significant differences were absent for mass soil water content, gravel coverage, and surface gravel size. The mosaic DP inhibited water infiltration by the pavement layer, where the sorptivity (S), initial infiltration rate (i_int), steady-state infiltration rate (i_sat) and infiltration time (T) averaged 1.19 cm/min^-0.5, 0.64 cm/min, 0.13 cm/min and 12.76 min, respectively. Where the pavement layer was scalped, the S, i_int, and i_sat increased by 0.27 cm/min^-0.5, 0.52 cm/min, and 0.40 cm/min, respectively, and the T reduced by 7.42 min. Water infiltration was mainly controlled by the pavement layer thickness (+), Av thickness (−), surface gravel coverage (−), fine earth (+) and fine gravel (−) in the pavement layer. The DP surfaces only had a sparse covering of shrubs, but an abundance of herbs. Few shrubs were present on the mosaic DP, but a greater number of shrubs and herbs grew on the mosaic BG. It can be concluded that DPs can maintain vegetation stability for different surface mosaic patterns. This study deepens our understanding of the eco-hydrological cycle of DP landscapes in arid regions.     

Mantle dynamics of the Paleoproterozoic North China Craton: A perspective based on seismic tomography

We investigate the mantle dynamics beneath the North China Craton (NCC) and surrounding regions based on a synthesis of recent P-wave mantle tomographic data down to depths of 600–800 km and their correlation with the surface geological features, with particular reference to the Paleoproterozoic tectonic events associated with the incorporation of the NCC within the Columbia supercontinent amalgam. From the tomographic images, we identify a hot corridor in the mantle transition zone beneath the central region of the Western Block of the NCC sandwiched between two cold corridors. This scenario is similar to the donut-shaped high-velocity anomaly surrounding a region of low-velocity anomaly in the lowermost mantle under the Pacific and suggests that the cold regions might represent slab graveyards which provide the fuel for the plumes rising from the center. A tomographic transect along the collisional suture of the NCC with the Columbia supercontinent, covering the Yinshan-Ordos Blocks in the Western Block through the Central Orogenic Belt and into the Eastern Block of the NCC reveals a ca. 250 km thick lithospheric keel below the Ordos Block defined by a prominent high-velocity anomaly. We identify slab break-off and asthenospheric upwelling in this region and suggest that this process probably initiated the thermal and material erosion of the tectosphere beneath the Eastern Block from the Paleoproterozoic, which was further intensified during the Mesozoic when a substantial part of the sub-continental mantle lithosphere was lost. We visualize heat input from asthenosphere and interaction between asthenosphere and overlying carbonated tectosphere releasing CO₂-rich fluids for the preservation of ultra-high temperature (ca. 1000 °C) metamorphic rocks enriched in CO₂ as well as high-pressure mafic granulites as a paired suite in this region. We also identify a hot swell of the asthenosphere rooted to more than 200 km depth and reaching up to the shallow mantle in the tomographic section along 35°N latitude at a depth of 800 km. This zone represents a cross-section through the southern part of the NCC. The surface distribution of Paleoproterozoic Xiong’er lavas and mafic dykes in this region would indicate that this region might have evidenced similar upwellings in the past. Our study has important implications in understanding the evolution of the NCC and suggests that the extensive modification of the mantle architecture and lithospheric structure beneath one of the fundamental Precambrian nuclei of Asia had a prolonged history probably dating from the Paleoproterozoic suturing of the NCC within the Columbia supercontinent amalgam.     

Theoretical analyses of acidization dissolution front instability in fluid‐saturated carbonate rocks

This paper presents an instability theory that can be used to understand the fundamental behavior of an acidization dissolution front when it propagates in fluid‐saturated carbonate rocks. The proposed theory includes two fundamental concepts, namely the intrinsic time and length of an acidization dissolution system, and a theoretical criterion that involves the comparison of the Zhao number and its critical value of the acidization dissolution system. The intrinsic time is used to determine the time scale at which the acidization dissolution front is formed, while the intrinsic length is used to determine the length scale at which the instability of the acidization dissolution front can be initiated. Under the assumption that the acidization dissolution reaction is a fast process, the critical Zhao number, which is used to assess the instability likelihood of an acidization dissolution front propagating in fluid‐saturated carbonate rocks, has been derived in a strictly mathematical manner. Based on the proposed instability theory of a propagating acidization dissolution front, it has been theoretically recognized that: (i) the increase of the mineral dissolution ratio can stabilize the acidization dissolution front in fluid‐saturated carbonate rocks; (ii) the increase of the final porosity of the carbonate rock can destabilize the acidization dissolution front, while the increase of the initial porosity can stabilize the acidization dissolution front in fluid‐saturated carbonate rocks; (iii) the increase of the mineral dissolution ratio can cause an increase in the dimensionless propagation speed of the acidization dissolution front; (iv) the increase of the initial porosity can enable the acidization dissolution front to propagate faster, while the increase of the final porosity can enable the acidization dissolution front to propagate slower in the acidization dissolution system. Copyright © 2012 John Wiley & Sons, Ltd.     

标准系列制备的一种简单方法

一些新的稀释装置已用于原子发射光谱法和原子吸收分光光度法中,改变了以往预先配制浓度依次变化的标准系列的习惯做法。Ko(?)cielniak提出的一种非连续稀释方法,所配制溶液的浓度与稀释次数成指数关系。本文采用了这种技术对稀释装置作了改进,只需要一个小锥形瓶,添加或稀释同一个原始溶液便可制备出一系列标准溶液。设备简单、操作方便,用于原子吸收分光光度法测定钾长石中的钠、钾并与常规法对照结果满意。