Numerical Method for Modeling the Constitutive Relationship of Sand under Different Stress Paths

A numerical method was used in order to establish the constitutive relationship of sands under different stress paths, Firstly, based on the numerical method modeling the constitutive law of sands, the elastoplastic constitutive relationship of sand was established for three paths： the constant proportion of principle stress path, the conventional triaxial compression （CTC） path, and the p=constant （TC） path. The yield lines of plastic volumetric strain and plastic generalized shear strain were given. Through visualization, the three dimensional surface of the stress-strain relationship in the whole stress field （p, q） obtained under the three paths was plotted. Also, by comparing the stress-strain surfaces and yield locus of the three stress paths, the differences were found to be obvious, which demonstrates that the influence of the stress paths on constitutive law was not neglected. The numerical modeling method overcame the difficulty of finding an analytical expression for plastic potential. The results simulated the experimental data with an accuracy of 90% on average, so the constitutive model established in this paper provides an effective constitutive equation for this kind of engineering, reflecting the effect of practical stress paths that occur in sands.

A Numerical Modelling Method of Fractured Reservoirs with Embedded Meshes and Topological Fracture Projection Configurations

Projection-based embedded discrete fracture model(pEDFM)is an effective numerical model to handle the flow in fractured reservoirs,with high efficiency and strong generalization of flow models.However,this paper points out that pEDFM fails to handle flow barriers in most cases,and identifies the physical projection configuration of fractures is a key step in pEDFM.This paper presents and proves the equivalence theorem,which explains the geometric nature of physical projection configurations of fractures,that is,the projection configuration of a fracture being physical is equivalent to it being topologically homeomorphic to the fracture,by analyzing the essence of pEDFM.Physical projection configurations of fractures may be rigorously established based on this theorem,allowing pEDFM to obtain physical numerical results for many flow models,particularly those with flow barriers.Furthermore,a natural idea emerges of employing flow barriers to flexibly‘cut’the formation to quickly handle the flow problems in the formation with complex geological conditions,and several numerical examples are implemented to test this idea and application of the improved pEDFM.

Decoupling Conditions for Elasto-plastic Consolidation Question Based on Numerical Modeling Method

Elasto-plastic consolidation is one of the classic coupling questions in geomechanics. To solve this problem, an elasto-plastic constitutive model is derived based on the numerical modeling method. The model is applied to Blot＇s consolidation theory. Incremental governing partial differential equations are established using this method. According to the stress path, the decoupling condition of these equations is discussed. Based on these conditions, an incremental diffusion equation and uncoupling governing equations are presented. The method is then applied to numerical analyses of three examples. The results show that （1） the effect of the stress path should be taken into account in the simulation of the soil consolidation question; （2） this decoupling method can predict the evolvement of pore water pressure; （3） the settlement using cam-clay model is less than that using numerical model because of dilatancy.

A High-Order Conservative Numerical Method for Gross-Pitaevskii Equation with Time-Varying Coefficients in Modeling BEC

We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation （BEC）. This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.

Acoustic viscoelastic modeling by frequency-domain boundary element method

Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green＇s function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method.

A review of discrete modeling techniques for fracturing processes in discontinuous rock masses

The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.

3D random Voronoi grain-based models for simulation of brittle rock damage and fabric-guided micro-fracturing

A grain-based distinct element model featuring three-dimensional （3D） Voronoi tessellations （randompoly-crystals） is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rockand allow for numerical replication of crack damage progression through initiation and propagation ofmicro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the pastfor brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi modelshas limited its application to two-dimensional （2D） codes. The proposed approach is implemented inNeper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files thatcan be read directly into 3DEC. A series of Unconfined Compressive Strength （UCS） tests are simulated in3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate therelationship between each micro-parameter and the model＇s macro-response. The possibility of numericalreplication of the classical U-shape strength curve for anisotropic rocks is also investigated innumerical UCS tests by using complex-shaped （elongated） grains that are cemented to one another alongtheir adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models foraccurate replication of the mechanical behaviour of isotropic and anisotropic （containing a fabric） rocks. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Fracture development around deep underground excavations: Insights from FDEM modelling

Over the past twenty years, there has been a growing interest in the development of numerical modelsthat can realistically capture the progressive failure of rock masses. In particular, the investigation ofdamage development around underground excavations represents a key issue in several rock engineeringapplications, including tunnelling, mining, drilling, hydroelectric power generation, and the deepgeological disposal of nuclear waste. The goal of this paper is to show the effectiveness of a hybrid finitediscreteelement method （FDEM） code to simulate the fracturing mechanisms associated with theexcavation of underground openings in brittle rock formations. A brief review of the current state-of-theartmodelling approaches is initially provided, including the description of selecting continuum- anddiscontinuum-based techniques. Then, the influence of a number of factors, including mechanical and insitu stress anisotropy, as well as excavation geometry, on the simulated damage is analysed for threedifferent geomechanical scenarios. Firstly, the fracture nucleation and growth process under isotropicrock mass conditions is simulated for a circular shaft. Secondly, the influence of mechanical anisotropy onthe development of an excavation damaged zone （EDZ） around a tunnel excavated in a layered rockformation is considered. Finally, the interaction mechanisms between two large caverns of an undergroundhydroelectric power station are investigated, with particular emphasis on the rock mass responsesensitivity to the pillar width and excavation sequence. Overall, the numerical results indicate that FDEMsimulations can provide unique geomechanical insights in cases where an explicit consideration offracture and fragmentation processes is of paramount importance. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Workflow to numerically reproduce laboratory ultrasonic datasets

The risks and uncertainties related to the storage of high-level radioactive waste （HLRW） can be reducedthanks to focused studies and investigations. HLRWs are going to be placed in deep geological repositories,enveloped in an engineered bentonite barrier, whose physical conditions are subjected tochange throughout the lifespan of the infrastructure. Seismic tomography can be employed to monitor itsphysical state and integrity. The design of the seismic monitoring system can be optimized via conductingand analyzing numerical simulations of wave propagation in representative repository geometry.However, the quality of the numerical results relies on their initial calibration. The main aim of this paperis to provide a workflow to calibrate numerical tools employing laboratory ultrasonic datasets. The finitedifference code SOFI2D was employed to model ultrasonic waves propagating through a laboratorysample. Specifically, the input velocity model was calibrated to achieve a best match between experimentaland numerical ultrasonic traces. Likely due to the imperfections of the contact surfaces, theresultant velocities of P- and S-wave propagation tend to be noticeably lower than those a prioriassigned. Then, the calibrated model was employed to estimate the attenuation in a montmorillonitesample. The obtained low quality factors （Q） suggest that pronounced inelastic behavior of the clay has tobe taken into account in geophysical modeling and analysis. Consequently, this contribution should beconsidered as a first step towards the creation of a numerical tool to evaluate wave propagation innuclear waste repositories. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Direct numerical simulation of particle-fluid systems by combining time-driven hard-sphere model and lattice Boltzmann method

A coupled numerical method for the direct numerical simulation of particle-fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by the time-driven hard-sphere model, while the hydrodynamic equations governing fluid flow are solved by the lattice Boltzmann method （LBM）, Particle-fluid coupling is realized by an immersed boundary method （IBM）, which considers the effect of boundary on surrounding fluid as a restoring force added to the governing equations of the fluid. The proposed scheme is validated in the classical flow-around-cylinder simulations, and preliminary application of this scheme to fluidization is reported, demonstrating it to be a promising computational strategy for better understanding complex behavior in particle-fluid systems.

Numeric modelling and risk assessment of pollutions in the Chinese Bohai Sea

To simulate the dynamic process of total nitrogen(TN) in seas, numerical modelling combined with the adjoint method is implemented in this study. Because nonpoint source terms(ST) and initial values(IV) of TN are essential but difficult to determine, the adjoint method was applied to a numerical model, and the ST and IV terms of TN were inverted via routine monitoring data in the Bohai Sea. In twin experiments, the adjoint method was capable of inverting the prescribed spatio-temporally distributed ST and the spatial distributed IV. In practical experiments, the results demonstrated that the simulation precision with ST inversion was higher than that with IV inversion and was accurate with joint initial values and source term(IST) inversion. This result indicates that nonpoint source TN is essential for the simulation of TN concentration. Furthermore, the simulated results indicate that the pollution in three bays of the Bohai Sea is rather severe. The model in this study is not specific to the Bohai Sea and can be generalized to other areas, such as the Beibu Gulf. These findings may assist in the development of cost-effective controls for accidental or planned industrial pollutant releases into coastal waters.