Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.     

水布垭面板堆石坝的三维弹塑性数值分析研究

在建的清江水布垭面板堆石坝高达233 m,是目前同类坝型中最高的.采用MSC.Marc非线性有限元程序,发展了三维子模型法,对该坝进行了三维弹塑性有限元仿真分析,模拟了面板的分缝、坝体材料分区、填筑及蓄水过程,采用双屈服面弹塑性模型模拟堆石体的变形特征.根据数值分析的结果,对坝体和面板的应力变形分布规律进行了探讨.     

Numerical simulation for 1975 Haicheng and 1999 Xiuyan earthquake processes by DDA+FEM

Introduction The Xiuyan earthquake(MS=5.4)on Nov.29,1999,occurred between Hushan town and Pi-anling town at the southeast end of the aftershocks of Haicheng earthquake(MS=7.3).These two earthquakes′sequences are highly similar in the way that the foreshocks are very frequent,focal mechanisms of the main shocks are about the same and fault strikes are approximately in the same direction.It is important to study these earthquakes for further research on seismogenesis to fore-cast earthquake…     

Mercury and methylmercury biogeochemistry in a thawing permafrost wetland complex,Northwest Territories,Canada

In arctic and sub‐arctic environments, mercury (Hg), more specifically toxic methylmercury (MeHg), is of growing concern to local communities because of its accumulation in fish. In these regions, there is particular interest in the potential mobilization of atmospherically deposited Hg sequestered in permafrost that is thawing at unprecedented rates. Permafrost thaw and the resulting ground surface subsidence transforms forested peat plateaus into treeless and permafrost‐free thermokarst wetlands where inorganic Hg released from the thawed permafrost and draining from the surrounding peat plateaus may be transformed to MeHg. This study begins to characterize the spatial distribution of MeHg in a peat plateau–thermokarst wetland complex, a feature that prevails throughout the wetland‐dominated southern margin of thawing discontinuous permafrost in Canada's Northwest Territories. We measured pore water total Hg, MeHg, dissolved organic matter characteristics and general water chemistry parameters to evaluate the role of permafrost thaw on the pattern of water chemistry. A gradient in vegetation composition, water chemistry and dissolved organic matter characteristics followed a toposequence from the ombrotrophic bogs near the crest of the complex to poor fens at its downslope margins. We found that pore waters in poor fens contained elevated levels of MeHg, and the water draining from these features had dissolved MeHg concentrations 4.5 to 14.5 times higher than the water draining from the bogs. It was determined through analysis of historical aerial images that the poor fens in the toposequence had formed relatively recently (early 1970s) as a result of permafrost thaw. Differences between the fens and bogs are likely to be a result of their differences in groundwater function, and this suggests that permafrost thaw in this landscape can result in hotspots for Hg methylation that are hydrologically connected to downstream ecosystems. Copyright © 2016 John Wiley & Sons, Ltd.     

Validating the ability of the discontinuous deformation analysis method to model normal P‐wave propagation across rock fractures

The effects of fractures on wave propagation problems are increasingly abstracting the attention of scholars and engineers in rock engineering field. This study aims to fully validate the ability of discontinuous deformation analysis (DDA) to model normal P‐wave propagation across rock fractures. The effects of a single fracture and multiple parallel fractures are all tested. The results indicate that DDA can accurately reflect the fracture effects, including the fractures stiffness, the fracture spacing and the fracture number, and the effects of incident wave frequency on one‐dimensional P‐wave propagation problems. Thus, DDA is able to deal well with normal incident P‐wave propagation problems. Copyright © 2017 John Wiley & Sons, Ltd.     

A discontinuous finite element suspended sediment transport model for water quality assessments in river networks

Suspended sediment plays an important role in the distribution and transport of many pollutants (such as radionuclides) in rivers. Pollutants may adsorb on fine suspended particles (e.g. clay) and spread according to the suspended sediment movement. Hence, the simulation of the suspended sediment mechanism is indispensable for realistic transport modelling. This paper presents and tests a simple mathematical model for predicting the suspended sediment transport in river networks. The model is based on the van Rijn suspended load formula and the advection–diffusion equation with a source or sink term that represents the erosion or deposition fluxes. The transport equation is solved numerically with the discontinuous finite element method. The model evaluation was performed in two steps, first by comparing model simulations with the measured suspended sediment concentrations in the Grote Nete–Molse Nete River in Belgium, and second by a model intercomparison with the sediment transport model NST MIKE 11. The simulations reflect the measurements with a Nash‐Sutcliffe model efficiency of 0.6, while the efficiency between the proposed model and the NST MIKE 11 simulations is 0.96. Both evaluations indicate that the proposed sediment transport model, that is sufficiently simple to be practical, is providing realistic results.     

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.     

Symplectitic augite: evidence for discontinuous precipitation as an exsolution mechanism in Ca-rich clinopyroxene

Abstract Part of the augite in the Artfjället gabbro consists of symplectitic intergrowths between augite and blebs or lamellae of orthopyroxene. Mineral compositions are consistent with formation of these symplectites by exsolution of orthopyroxene from magmatic augite at a temperature of ca . 900–1000°C. The microstructures indicate that the exsolution mechanism is discontinuous precipitation, whereby the boundary of an augite grain sweeps through a neighbouring augite, leaving the symplectite in its wake. The formation of this symplectitic augite is catalysed by the presence of an intergranular water-rich fluid phase, which promotes grain boundary mobility.     

Three‐parameter discontinuous distributions for hydrological samples with zero values

A consistent approach to the frequency analysis of hydrologic data in arid and semiarid regions, i.e. the data series containing several zero values (e.g. monthly precipitation in dry seasons, annual peak flow discharges, etc.), requires using discontinuous probability distribution functions. Such an approach has received relatively limited attention. Along the lines of physically based models, the extensions of the Muskingum‐based models to three parameter forms are considered. Using 44 peak flow series from the USGS data bank, the fitting ability of four three‐parameter models was investigated: (1) the Dirac delta combined with Gamma distribution; (2) the Dirac delta combined with two‐parameter generalized Pareto distribution; (3) the Dirac delta combined with two‐parameter Weibull (DWe) distribution; (4) the kinematic diffusion with one additional parameter that controls the probability of the zero event (KD3). The goodness of fit of the models was assessed and compared both by evaluation of discrepancies between the results of both estimation methods (i.e. the method of moments (MOM) and the maximum likelihood method (MLM)) and using the log of likelihood function as a criterion. In most cases, the DWe distribution with MLM‐estimated parameters showed the best fit of all the three‐parameter models. Copyright © 2005 John Wiley & Sons, Ltd.     

Nuclear Waste Disposal Simulations: Couplex Test Cases

We give some results obtained for the Couplex test cases proposed by the ANDRA. In this paper our aim is twofold. Firstly, to compute the release of nuclides out of the repository by concentrating on the 3D near field (Couplex 2). The simulation of the transport phenomena takes into account the dissolution of the glass containers and congruent emissions of the radio-nuclides including filiation chains and some simplified chemistry. Secondly, it is to use the near field computations in order to simulate the nuclide migrations in a 2D far field (Couplex 3). Coupling in between the two simulations takes into consideration the periodicity of the disposal modules and the geometry of the repository described in Couplex 1. The mixed finite element and discontinuous Galerkin methods are used to solve the convection–diffusion equations. In order to handle the nonlinear precipitation/dissolution term, we developed a new iterative technique that combines Picard and Newton–Raphson methods.