Spurious forces are a significant challenge for multi-scale methods,e.g.,the coupled atomistic/discrete dislocation(CADD)method.The assumption of isotropic matter in the continuum domain is a critical factor leading t...Spurious forces are a significant challenge for multi-scale methods,e.g.,the coupled atomistic/discrete dislocation(CADD)method.The assumption of isotropic matter in the continuum domain is a critical factor leading to such forces.This study aims to minimize spurious forces,ensuring that atomic dislocations experience more precise forces from the continuum domain.The authors have already implemented this idea using a simplified and unrealistic slipping system.To create a comprehensive and realistic model,this paper considers all possible slip systems in the face center cubic(FCC)lattice structure,and derives the required relationships for the displacement fields.An anisotropic version of the three-dimensional CADD(CADD3D)method is presented,which generates the anisotropic displacement fields for the partial dislocations in all the twelve slip systems of the FCC lattice structure.These displacement fields are tested for the most probable slip systems of aluminum,nickel,and copper with different anisotropic levels.Implementing these anisotropic displacement fields significantly reduces the spurious forces on the slip systems of FCC materials.This improvement is particularly pronounced at greater distances from the interface and in more anisotropic materials.Furthermore,the anisotropic CADD3D method enhances the spurious stress difference between the slip systems,particularly for materials with higher anisotropy.展开更多
Numerical analyses of earthquake effects on the deformation, stability, and load transfer of a slope covered by deposits are traditionally based on the assumption that the slope is a continuum. It would be problem...Numerical analyses of earthquake effects on the deformation, stability, and load transfer of a slope covered by deposits are traditionally based on the assumption that the slope is a continuum. It would be problematic, however, to extend these approaches to the simulation of the slide, collapse and disintegration of the deposits under seismic loading. Contrary to this, a discrete element method (DEM) provides a means to consider large displacement and rotation of the non-continuum. To take the advantages of both methods of continuum and non- continuum analyses, seismic responses of a slope covered by deposits are studied by coupling a twodimensional (a-D) finite difference method and a 2-D DEM, with the bedrock being modelled by the finite difference grids and the deposits being represented by disks. A smooth transition across the boundaries of the continuous/discontinuous domains is obtained by imposing the compatibility condition and equilibrium condition along their interfaces. In the course of computation, the same time-step value is chosen for both continuous and discontinuous domains. The free-field boundaries are adopted for lateral grids of bedrock domain to eliminate the radiation damping effect. When the static equilibrium under gravity load is obtained, dynamic calculation begins under excitation of the seismic wave input from the continuum model bottom. In this way, responses to the earthquake of a slope covered by deposits are analyzed dynamically. Combined with field monitoring data, deformation and stability of the slope are discussed. The effects of the relevant parameters of spectrum characteristic, duration, andpeak acceleration of seismic waves are further investigated and explained from the simulations.展开更多
In this paper,we develop the residual based a posteriori error estimates and the corresponding adaptive mesh refinement algorithm for atomistic/continuum(a/c)coupling with finite range interactions in two dimensions.W...In this paper,we develop the residual based a posteriori error estimates and the corresponding adaptive mesh refinement algorithm for atomistic/continuum(a/c)coupling with finite range interactions in two dimensions.We have systematically derived a new explicitly computable stress tensor formula for finite range in-teractions.In particular,we use the geometric reconstruction based consistent atomistic/continuum(GRAC)coupling scheme,which is quasi-optimal if the continuum model is discretized by P1 finite elements.The numerical results of the adaptive mesh refinement algorithm is consistent with the quasi-optimal a priori error estimates.展开更多
Weakly bound states often occur in nuclear physics.To precisely understand their properties,the coupling to the continuum should be worked out explicitly.As the first step,we use a simple nuclear model in the continuu...Weakly bound states often occur in nuclear physics.To precisely understand their properties,the coupling to the continuum should be worked out explicitly.As the first step,we use a simple nuclear model in the continuum and on a lattice to investigate the influence of a third particle on a loosely bound state of a particle and a heavy core.Our approach is consistent with the Lüscher formalism.展开更多
文摘Spurious forces are a significant challenge for multi-scale methods,e.g.,the coupled atomistic/discrete dislocation(CADD)method.The assumption of isotropic matter in the continuum domain is a critical factor leading to such forces.This study aims to minimize spurious forces,ensuring that atomic dislocations experience more precise forces from the continuum domain.The authors have already implemented this idea using a simplified and unrealistic slipping system.To create a comprehensive and realistic model,this paper considers all possible slip systems in the face center cubic(FCC)lattice structure,and derives the required relationships for the displacement fields.An anisotropic version of the three-dimensional CADD(CADD3D)method is presented,which generates the anisotropic displacement fields for the partial dislocations in all the twelve slip systems of the FCC lattice structure.These displacement fields are tested for the most probable slip systems of aluminum,nickel,and copper with different anisotropic levels.Implementing these anisotropic displacement fields significantly reduces the spurious forces on the slip systems of FCC materials.This improvement is particularly pronounced at greater distances from the interface and in more anisotropic materials.Furthermore,the anisotropic CADD3D method enhances the spurious stress difference between the slip systems,particularly for materials with higher anisotropy.
基金the National Basic Research Program of China (Grant No. 2008CB425802)
文摘Numerical analyses of earthquake effects on the deformation, stability, and load transfer of a slope covered by deposits are traditionally based on the assumption that the slope is a continuum. It would be problematic, however, to extend these approaches to the simulation of the slide, collapse and disintegration of the deposits under seismic loading. Contrary to this, a discrete element method (DEM) provides a means to consider large displacement and rotation of the non-continuum. To take the advantages of both methods of continuum and non- continuum analyses, seismic responses of a slope covered by deposits are studied by coupling a twodimensional (a-D) finite difference method and a 2-D DEM, with the bedrock being modelled by the finite difference grids and the deposits being represented by disks. A smooth transition across the boundaries of the continuous/discontinuous domains is obtained by imposing the compatibility condition and equilibrium condition along their interfaces. In the course of computation, the same time-step value is chosen for both continuous and discontinuous domains. The free-field boundaries are adopted for lateral grids of bedrock domain to eliminate the radiation damping effect. When the static equilibrium under gravity load is obtained, dynamic calculation begins under excitation of the seismic wave input from the continuum model bottom. In this way, responses to the earthquake of a slope covered by deposits are analyzed dynamically. Combined with field monitoring data, deformation and stability of the slope are discussed. The effects of the relevant parameters of spectrum characteristic, duration, andpeak acceleration of seismic waves are further investigated and explained from the simulations.
基金supported by National Natural Science Foundation of China grant 11861131004,11771040,91430106supported by Natural Science Foundation of China grant 11871339,11861131004,11571314,11471214 and the One Thousand Plan of China for young scientists.
文摘In this paper,we develop the residual based a posteriori error estimates and the corresponding adaptive mesh refinement algorithm for atomistic/continuum(a/c)coupling with finite range interactions in two dimensions.We have systematically derived a new explicitly computable stress tensor formula for finite range in-teractions.In particular,we use the geometric reconstruction based consistent atomistic/continuum(GRAC)coupling scheme,which is quasi-optimal if the continuum model is discretized by P1 finite elements.The numerical results of the adaptive mesh refinement algorithm is consistent with the quasi-optimal a priori error estimates.
基金financial support from the Deutsche Forschungsgemeinschaft (SFB/TRR 110,"Symmetries and the Emergence of Structure in QCD",grant no. TRR 110)by the Chinese Academy of Sciences (CAS) President’s International Fellowship Initiative (PIFI)(2018DM0034)+1 种基金by Volkswagen Stiftung (93562)by the Fundamental Research Funds for the Central Universities
文摘Weakly bound states often occur in nuclear physics.To precisely understand their properties,the coupling to the continuum should be worked out explicitly.As the first step,we use a simple nuclear model in the continuum and on a lattice to investigate the influence of a third particle on a loosely bound state of a particle and a heavy core.Our approach is consistent with the Lüscher formalism.