Nonlinear wave dispersion in monoatomic chains with lumped and distributed masses:discrete and continuum models

The main objective of this paper is to investigate the influence of inertia of nonlinear springs on the dispersion behavior of discrete monoatomic chains with lumped and distributed masses.The developed model can represent the wave propagation problem in a non-homogeneous material consisting of heavy inclusions embedded in a matrix.The inclusions are idealized by lumped masses,and the matrix between adjacent inclusions is modeled by a nonlinear spring with distributed masses.Additionally,the model is capable of depicting the wave propagation in bi-material bars,wherein the first material is represented by a rigid particle and the second one is represented by a nonlinear spring with distributed masses.The discrete model of the nonlinear monoatomic chain with lumped and distributed masses is first considered,and a closed-form expression of the dispersion relation is obtained by the second-order Lindstedt-Poincare method(LPM).Next,a continuum model for the nonlinear monoatomic chain is derived directly from its discrete lattice model by a suitable continualization technique.The subsequent use of the second-order method of multiple scales(MMS)facilitates the derivation of the corresponding nonlinear dispersion relation in a closed form.The novelties of the present study consist of(i)considering the inertia of nonlinear springs on the dispersion behavior of the discrete mass-spring chains;(ii)developing the second-order LPM for the wave propagation in the discrete chains;and(iii)deriving a continuum model for the nonlinear monoatomic chains with lumped and distributed masses.Finally,a parametric study is conducted to examine the effects of the design parameters and the distributed spring mass on the nonlinear dispersion relations and phase velocities obtained from both the discrete and continuum models.These parameters include the ratio of the spring mass to the lumped mass,the nonlinear stiffness coefficient of the spring,and the wave amplitude.

Discrete Element Modelling of Damage Evolution of Concrete Considering Meso-Structure of ITZ

The mechanical properties of interfacial transition zones(ITZs)have traditionally been simplified by reducing the stiffness of cement in previous simulation methods.A novel approach based on the discrete element method(DEM)has been developed for modeling concrete.This new approach efficiently simulates the meso-structure of ITZs,accurately capturing their heterogeneous properties.Validation against established uniaxial compression experiments confirms the precision of thismodel.The proposedmodel canmodel the process of damage evolution containing cracks initiation,propagation and penetration.Under increasing loads,cracks within ITZs progressively accumulate,culminating in macroscopic fractures that traverse themortarmatrix,forming the complex,serpentine path of cracks.This study reveals four distinct displacement patterns:tensile compliant,tensile opposite,mixed tensile-shear,and shear opposite patterns,each indicative of different stages in concrete’s damage evolution.The widening angle of these patterns delineates the progression of cracks,with the tensile compliant pattern signaling the initial crack appearance and the shear opposite pattern indicating the concrete model’s ultimate failure.

A discrete Boltzmann model with symmetric velocity discretization for compressible flow

A discrete Boltzmann model(DBM) with symmetric velocity discretization is constructed for compressible systems with an adjustable specific heat ratio in the external force field. The proposed two-dimensional(2D) nine-velocity scheme has better spatial symmetry and numerical accuracy than the discretized velocity model in literature [Acta Aerodyn. Sin.40 98108(2022)] and owns higher computational efficiency than the one in literature [Phys. Rev. E 99 012142(2019)].In addition, the matrix inversion method is adopted to calculate the discrete equilibrium distribution function and force term, both of which satisfy nine independent kinetic moment relations. Moreover, the DBM could be used to study a few thermodynamic nonequilibrium effects beyond the Euler equations that are recovered from the kinetic model in the hydrodynamic limit via the Chapman–Enskog expansion. Finally, the present method is verified through typical numerical simulations, including the free-falling process, Sod’s shock tube, sound wave, compressible Rayleigh–Taylor instability,and translational motion of a 2D fluid system.

Application of shifted lattice model to 3D compressible lattice Boltzmann method

An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D)shifted lattice model is proposed,therefore a shifted lattice model is introduced into D3Q25.Under the finite volume scheme,several typical compressible calculation examples are used to verify whether the numerical stability of the D3Q25 model can be improved by adding the shifted lattice model.The simulation results show that the numerical stability is indeed improved after adding the shifted lattice model.

Influence of heterogeneity on rock strength and stiffness using discrete element method and parallel bond model

The particulate discrete element method(DEM) can be employed to capture the response of rock,provided that appropriate bonding models are used to cement the particles to each other.Simulations of laboratory tests are important to establish the extent to which those models can capture realistic rock behaviors.Hitherto the focus in such comparison studies has either been on homogeneous specimens or use of two-dimensional(2D) models.In situ rock formations are often heterogeneous,thus exploring the ability of this type of models to capture heterogeneous material behavior is important to facilitate their use in design analysis.In situ stress states are basically three-dimensional(3D),and therefore it is important to develop 3D models for this purpose.This paper revisits an earlier experimental study on heterogeneous specimens,of which the relative proportions of weaker material(siltstone) and stronger,harder material(sandstone) were varied in a controlled manner.Using a 3D DEM model with the parallel bond model,virtual heterogeneous specimens were created.The overall responses in terms of variations in strength and stiffness with different percentages of weaker material(siltstone) were shown to agree with the experimental observations.There was also a good qualitative agreement in the failure patterns observed in the experiments and the simulations,suggesting that the DEM data enabled analysis of the initiation of localizations and micro fractures in the specimens.

Microstructure modeling and virtual test of asphalt mixture based on three-dimensional discrete element method

The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.

Numerical Analysis of Ice Rubble with a Freeze-Bond Model in Dilated Polyhedral Discrete Element Method

Freezing in ice rubble is a common phenomenon in cold regions,which can consolidate loose blocks and change their mechanical properties.To model the cohesive effect in frozen ice rubble,and to describe the fragmentation behavior with a large external forces exerted,a freeze-bond model based on the dilated polyhedral discrete element method(DEM)is proposed.Herein,imaginary bonding is initialized at the contact points to transmit forces and moments,and the initiation of the damage is detected using the hybrid fracture model.The model is validated through the qualitative agreement between the simulation results and the analytical solution of two bonding particles.To study the effect of freeze-bond on the floating ice rubble,punch-through tests were simulated on the ice rubble under freezing and nonfreezing conditions.The deformation and resistance of the ice rubble are investigated during indenter penetration.The influence of the internal friction coefficient on the strength of the ice rubble is determined.The results indicate that the proposed model can properly describe the consolidated ice rubble,and the freeze-bond effect is of great significance to the ice rubble properties.

A modified discrete element method for concave granular materials based on energy-conserving contact model

The development of a general discrete element method for irregularly shaped particles is the core issue of the simulation of the dynamic behavior of granular materials.The general energy-conserving contact theory is used to establish a universal discrete element method suitable for particle contact of arbitrary shape.In this study,three dimentional(3D)modeling and scanning techniques are used to obtain a triangular mesh representation of the true particles containing typical concave particles.The contact volumebased energy-conserving model is used to realize the contact detection between irregularly shaped particles,and the contact force model is refined and modified to describe the contact under real conditions.The inelastic collision processes between the particles and boundaries are simulated to verify the robustness of the modified contact force model and its applicability to the multi-point contact mode.In addition,the packing process and the flow process of a large number of irregular particles are simulated with the modified discrete element method(DEM)to illustrate the applicability of the method of complex problems.

Simulation of a Laboratory Scale Ball Mill via Discrete Element Method Modelling

Discrete Element Method (DEM) is a powerful tool for simulating different types of mills. It also used for computing different types of particles such as rocks, grains, and molecules. DEM has been widely used in the field of rock mechanics. In the present work, DEM approach is applied to model the milling media (powder particles and balls) inside a planetary ball mill and to estimate the distribution of particles of a dry powder during milling. In fact, the efficiency of the DEM strongly depends on the input parameters. The DEM simulation results indicated that DEM is a promising tool for the simulation of the dynamic particles motion and interactions within planetary ball mill. These results could be utilized to further develop the synthesis performance, anticipate the reaction, and reduce the wear in the dry milling reactions.

Numerical Methods for Discrete Double Barrier Option Pricing Based on Merton Jump Diffusion Model

As a kind of weak-path dependent options, barrier options are an important kind of exotic options. Because the pricing formula for pricing barrier options with discrete observations cannot avoid computing a high dimensional integral, numerical calculation is time-consuming. In the current studies, some scholars just obtained theoretical derivation, or gave some simulation calculations. Others impose underlying assets on some strong assumptions, for example, a lot of calculations are based on the Black-Scholes model. This thesis considers Merton jump diffusion model as the basic model to derive the pricing formula of discrete double barrier option;numerical calculation method is used to approximate the continuous convolution by calculating discrete convolution. Then we compare the results of theoretical calculation with simulation results by Monte Carlo method, to verify their efficiency and accuracy. By comparing the results of degeneration constant parameter model with the results of previous models we verified the calculation method is correct indirectly. Compared with the Monte Carlo simulation method, the numerical results are stable. Even if we assume the simulation results are accurate, the time consumed by the numerical method to achieve the same accuracy is much less than the Monte Carlo simulation method.

Modeling calving process of glacier with dilated polyhedral discrete element method

Mass loss caused by glacier calving is one of the direct contributors to global sea level rise.Reliable calving laws are required for accurate modelling of ice sheet mass balance.Both continuous and discontinuous methods have been used for glacial calving simulations.In this study,the discrete element method(DEM)based on dilated polyhedral elements is introduced to simulate the calving process of a tidewater glacier.Dilated polyhedrons can be obtained from the Minkowski sum of a sphere and a core polyhedron.These elements can be utilized to generate a continuum ice material,where the interaction force between adjacent elements is modeled by constructing bonds at the joints of the common faces.A hybrid fracture model considering fracture energy is introduced.The viscous creep behavior of glaciers on long-term scales is not considered.By applying buoyancy and gravity to the modelled glacier,DEM results show that the calving process is caused by cracks which are initialized at the top of the glacier and spread to the bottom.The results demonstrate the feasibility of using the dilated polyhedral DEM method in glacier simulations,additionally allowing the fragment size of the breaking fragments to be counted.The relationship between crack propagation and internal stress in the glacier is analyzed during calving process.Through the analysis of the Mises stress and the normal stress between the elements,it is found that geometric changes caused by the glacier calving lead to the redistribution of the stress.The tensile stress between the elements is the main influencing factor of glacier ice failure.In addition,the element shape,glacier base friction and buoyancy are studied,the results show that the glacier model based on the dilated polyhedral DEM is sensitive to the above conditions.

THE APPLICATION OF DISCRETE ELEMENT METHOD IN SOLVING THREE-DIMENTIONAL IMPACT DYNAMICS PROBLEMS

A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)to contact model(for non-continuum),is developed for simulating the mechanical process from continuum to non-continuum.The wave propagation process in a concrete block(as continuum)made of cement grout under impact loading is numer- ically simulated with this code.By comparing its numerical results with those by LS-DYNA,the calculation accuracy of the model and algorithm is proved.Furthermore,the failure process of the concrete block under quasi-static loading is demonstrated,showing the basic dynamic tran- sitional process from continuum to non-continuum.The results of calculation can be displayed by animation.The damage modes are similar to the experimental results.The two numerical examples above prove that our model and its code are powerful and efficient in simulating the dynamic failure problems accompanying the transition from continuum to non-continuum.It also shows that the discrete element method(DEM)will have broad prospects for development and application.

Understanding roof deformation mechanics and parametric sensitivities of coal mine entries using the discrete element method

Although conventional coal mine designs are conservative regarding pillar strength,local failures such as roof-falls and pillar bursts still affect mine safety and operations.Previous studies have identified that discontinuous,layered roof materials have some self-supporting capacity.This research is a preliminary step towards understanding these mechanics in coal-measure rocks.Although others have considered broad conceptual models and simplified analogs for mine roof behavior,this study presents a unique numerical model that more completely represents in-situ roof conditions.The discrete element method(DEM)is utilized to conduct a parametric analysis considering a range of in-situ stress ratios,material properties,and joint networks to determine the parameters controlling the stability of single-entries modeled in two-dimensions.Model results are compared to empirical observations of roof-support effectiveness(ARBS)in the context of the coal mine roof rating(CMRR)system.Results such as immediate roof displacement,overall stability,and statistical relationships between model parameters and outcomes are presented herein.Potential practical applications of this line of research include:(1)roof-support optimization for a range of coal-measure rocks,(2)establishment of a relationship between roof stability and pillar stress,and(3)determination of which parameters are most critical to roof stability and therefore require concentrated evaluation.

Predicting the Dynamic Behavior of Asphalt Concrete Using Three-dimensional Discrete Element Method

A user-defined three-dimensional （3D） discrete element model was presented to predict the dynamic modulus and phase angle of asphalt concrete （AC）. The 3D discrete element method （DEM） model of AC was constructed employing a user-defined computer program developed using the ＂Fish＂ language in PFC3D. Important microstructural features of AC were modeled, including aggregate gradation, air voids and mastic. The irregular shape of aggregate particle was modeled using a clump of spheres. The developed model was validated through comparing with experimental measurements and then used to simulate the cyclic uniaxial compression test, based on which the dynamic modulus and phase angle were calculated from the output stress- strain relationship. The effects of air void content, aggregate stiffness and volumetric fraction on AC modulus were further investigated. The experimental results show that the 3D DEM model is able to accurately predict both dynamic modulus and phase angle of AC across a range of temperature and loading frequencies. The user- defined 3D model also demonstrated significant improvement over the general existing two-dimensional models.

Influence of particle contact models on soil response of poorly graded sand during cavity expansion in discrete element simulation

The discrete element method(DEM) has been extensively adopted to investigate many complex geotechnical related problems due to its capability to incorporate the discontinuous nature of granular materials. In particular, when simulating large deformations or distortion of soil(e.g. cavity expansion),DEM can be very effective as other numerical solutions may experience convergence problems. Cavity expansion theory has widespread applications in geotechnical engineering, particularly to the problems concerning in situ testing, pile installation and so forth. In addition, the behaviour of geomaterials in a macro-level is utterly determined by microscopic properties, highlighting the importance of contact models. Despite the fact that there are numerous contact models proposed to mimic the realistic behaviour of granular materials, there are lack of studies on the effects of these contact models on the soil response.Hence, in this study, a series of three-dimensional numerical simulations with different contact constitutive models was conducted to simulate the response of sandy soils during cylindrical cavity expansion. In this numerical investigation, three contact models, i.e. linear contact model, rolling resistance contact model,and Hertz contact model, are considered. It should be noted that the former two models are linear based models, providing linearly elastic and frictional plasticity behaviours, whereas the latter one consists of nonlinear formulation based on an approximation of the theory of Mindlin and Deresiewicz. To examine the effects of these contact models, several cylindrical cavities were created and expanded gradually from an initial radius of 0.055 m to a final radius of 0.1 m. The numerical predictions confirm that the calibrated contact models produced similar results regarding the variations of cavity pressure, radial stress, deviatoric stress, volumetric strain, as well as the soil radial displacement. However, the linear contact model may result in inaccurate predictions when highly angular soil particles are involved. In addition, considering the excessive soil displacement induced by the pile installation(i.e. cavity expansion), a minimum distance of11 a(a is the cavity radius) is recommend for practicing engineers to avoid the potential damages to the existing piles and adjacent structures.

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.

A review of methods,applications and limitations for incorporating fluid flow in the discrete element method

The past decade has witnessed the substantial growth in research interests and progress on the subject of coupled hydro-mechanical processes in rocks and soils,driven mainly by the surge of research in unconventional hydrocarbon reservoirs and associated hazards.Many coupling techniques have been developed to include the effects of fluid flow in the discrete element method(DEM),and the techniques have been applied to a variety of geomechanical problems.Although these coupling methods have been successfully applied in various engineering fields,no single fluid/DEM coupling method is universal due to the complexity of engineering problems and the limitations of the numerical methods.For researchers and engineers,the key to solve a specific problem is to select the most appropriate fluid/DEM coupling method among these modeling technologies.The purpose of this paper is to give a comprehensive review of fluid flow/DEM coupling methods and relevant research.Given their importance,the availability or unavailability of best practice guidelines is outlined.The theoretical background and current status of DEM are introduced first,and the principles,applications,and advantages and disadvantages of different fluid flow/DEM coupling methods are discussed.Finally,a summary with speculation on future development trends is given.

A particle-resolved heat-particle-fluid coupling model by DEM-IMB-LBM

Multifield coupling is frequently encountered and also an active area of research in geotechnical engineering.In this work,a particle-resolved direct numerical simulation(PR-DNS)technique is extended to simulate particle-fluid interaction problems involving heat transfer at the grain level.In this extended technique,an immersed moving boundary(IMB)scheme is used to couple the discrete element method(DEM)and lattice Boltzmann method(LBM),while a recently proposed Dirichlet-type thermal boundary condition is also adapted to account for heat transfer between fluid phase and solid particles.The resulting DEM-IBM-LBM model is robust to simulate moving curved boundaries with constant temperature in thermal flows.To facilitate the understanding and implementation of this coupled model for non-isothermal problems,a complete list is given for the conversion of relevant physical variables to lattice units.Then,benchmark tests,including a single-particle sedimentation and a two-particle drafting-kissing-tumbling(DKT)simulation with heat transfer,are carried out to validate the accuracy of our coupled technique.To further investigate the role of heat transfer in particle-laden flows,two multiple-particle problems with heat transfer are performed.Numerical examples demonstrate that the proposed coupling model is a promising high-resolution approach for simulating the heat-particle-fluid coupling at the grain level.

Discrete Element Simulation of Asphalt Mastics Based on Burgers Model

In order to investigate the viscoelastic performance of asphalt mastics, a micro-mechanical model for asphalt mastics was built by applying Burgers model to discrete element simulation and constructing Burgers contact model. Then the numerical simulation of creep tests was conducted, and results from the simulation were compared with the analytical solution for Burgers model. The comparision showed that the two results agreed well with each other, suggesting that discrete element model based on Burgers model could be employed in the numerical simulation for asphalt mastics.

A discrete model for prediction of radon flux from fractured rocks

Prediction of radon flux from the fractured zone of a propagating cave mine is basically associated with uncertainty and complexity. For instance, there is restricted access to these zones for field measure- ments, and it is quite difficult to replicate the complex nature of both natural and induced fractures in these zones in laboratory studies. Hence, a technique for predicting radon flux from a fractured rock using a discrete fracture network （DFN） model is developed to address these difficulties. This model quantifies the contribution of fractures to the total radon flux, and estimates the fracture density from a measured radon flux considering the effects of advection, diffusion, as well as radon generation and decay. Radon generation and decay are classified as reaction processes. Therefore, the equation solved is termed as the advection-diffusion-reaction equation （ADRE）. Peclet number （Pe）, a conventional dimensionless parameter that indicates the ratio of mass transport by advection to diffusion, is used to classify the transport regimes. The results show that the proposed model effectively predicts radon flux from a fractured rock. An increase in fracture density for a rock sample with uniformly distributed radon generation rate can elevate radon flux significantly compared with another rock sample with an equivalent increase in radon generation rate. In addition to Pe, two other independent dimensionless parameters （derived for radon transport through fractures） significantly affect radon dimensionless flux. Findings provide insight into radon transport through fractured rocks and can be used to improve radon control measures for proactive mitigation.