This study introduces GeoTaichi, an open-source high-performance numerical simulator designed for addressing multiscale geophysical problems. By leveraging the power of the Taichi parallel language, ...GeoTaichi maximizes the utilization of modern computer resources on multicore CPU and GPU architectures. It offers robust and reliable modules for the discrete element method (DEM), material point method (MPM), and coupled material point-discrete element method (MPDEM). These modules enable efficient solving of large-scale problems while being implemented in pure Python. The design philosophy of GeoTaichi focuses on creating a framework that is readable, extensible, and user-friendly. This paper highlights the coupling procedure of MPDEM, the code structures, and the most important features of GeoTaichi. Rigorous benchmark tests have been conducted to verify the validity and robustness of GeoTaichi. Additionally, the performance of GeoTaichi is compared with similar software tools in the field, underscoring a notable improvement in both computational efficiency and memory savings when compared to existing alternatives.
Program title: GeoTaichi
CPC Library link to program files:https://doi.org/10.17632/858bmcf7j6.1
Developer's repository link:https://github.com/Yihao-Shi/GeoTaichi
Licensing provisions: GNU General Public License v3.0
Programming language: Python
Nature of problem: The simulations of large-deformation geophysical flows and their interaction with structures play a crucial role in the field of geophysics. To address the complexities of these nonlinear problems, the discrete element method (DEM), material point method (MPM), and their coupling (MPDEM) have proven to be highly suitable numerical schemes. However, these schemes impose substantial computational demands, necessitating the development of an efficient framework that can harness modern computer resources on multicore CPU and GPU architectures.
Solution method: The open-source code GeoTaichi implements the DEM, MPM, and coupled MPDEM, encompassing a range of constitutive models and contact laws for different geologic materials. The clump particle model is also introduced in DEM to solve granular mechanics involving complex-shaped particles. One significant advantage of GeoTaichi is its utilization of the Taichi parallel language, which is designed to be user-friendly and easily extensible for customized applications.
•A detailed FDEM numerical method to simulate mechanical and fracturing responses of heterogeneous geomaterials with irregular inclusions is systematically developed.•A computational geometry method ...named CWSVM is proposed to control mesh quantity and quality.•A signed-distance-field-based discrete element method (SDF-DEM) is employed to approach the natural allocation and orientation of inclusions.•A combined constitutive model is proposed to consider the shearing hardening behaviour for the cohesive elements.•Effects of the interface strength on the mechanical and fracturing behaviours of inclusion-containing geomaterials are extensively discussed.
In this paper, a detailed FDEM approach to simulate the mechanical and fracturing responses of heterogeneous geomaterials with irregular inclusions is systematically developed. The inclusion surface morphology is first obtained through 3D scanning techniques. A computational geometry method, the curvature-weighted sphere Voronoi method (CWSVM), is adopted to control the mesh quantity and quality and ensure the efficiency and accuracy of the FDEM numerical model. A signed-distance-field-based discrete element method (SDF-DEM) is employed to approximate the natural distribution and orientation of inclusions. Heterogeneous geomaterials with large inclusion contents (such as 60% and 70%) are generated effectively and efficiently through this approach. Next, to model the fracturing process, a finite discrete element method (FDEM) model is developed by integrating cohesive elements into the mesh in a fast and efficient manner. In addition, a combined constitutive model is proposed to consider the shear-hardening behaviour of the cohesive elements. The proposed numerical approach is verified through comparison with experimental results, including the shape of inclusions and mechanical responses of geomaterials. The results demonstrate that both satisfactory precision and low calculation costs can be achieved using the proposed algorithm. The consequent simulation performance is verified through comparisons of observations and numerical results with experimental results for failure patterns and mechanical behaviours. In addition, the effects of the strength of the interfaces between the inclusions and matrix on the mechanical and fracturing characteristics of inclusion-containing geomaterials are analysed quantitatively. The mechanical strength decreases rather than increases with increasing content of inclusions when the interface strength is less than the matrix strength.
•CFD-DEM formulation, including heat and mass transfer and long-range forces is described.•Implementation of CFD-DEM in simulation of different processes is discussed.•Different applications, ...including drying, coating, mixing combustion, gasification and etc. are discussed.
With increasing the computational resources, the number of publications about coupled computational fluid dynamics – discrete element method is in the rise in the recent years. This technique is very useful, especially in simulation of fluid-solid flows in process engineering. This paper provides an introduction to CFD-DEM modeling in process engineering systems, including heat and mass transfer and long range forces, and reviews the major researches in simulation of two-phase processes such as drying, coating, granulation, crystallization, chemical reactions (including combustion, gasification and pyrolysis) and mixing. Details of implementing unresolved CFD-DEM in these applications are explained in details and major assumptions and findings are discussed.
•A modelling framework for multi-track, multi-layer and multi-material SLM.•Simulation of multi-material powder deposition in various patterns.•Effect of process parameters on balling effect, keyhole ...and lack of fusion.•Molten pool evolution of multi-material SLM on the same and across different layers.•Modelling of phase migration at the interface between two different materials.
Selective laser melting (SLM) is a promising powder-based additive manufacturing technology due to its capability to fabricate metallic components with complex geometries. While most previous investigations focus on printing with a single material, recent industry-orientated studies indicate the need for multi-material SLM in several high-value manufacturing sectors including medical devices, aerospace and automotive industries. However, understanding the underlying physics in multi-material SLM remains challenging due to the difficulties of experimental observation. In this paper, an integrated modelling framework for multi-track, multi-layer and multi-material SLM is developed to advance the in-depth understanding of this process. The main novelty is in modelling the molten pool evolvement and track morphology of multiple materials deposited on the same and across different layers. Discrete element method (DEM) is employed to reproduce the powder deposition process of multiple materials in different deposition patterns, with particle size distribution imported from a particle size analyser. Various phenomena including balling effect, keyhole depression, and lack of fusion between layers are investigated with different laser energy inputs. As a result of the different thermal properties, several process parameters including energy density and hatch spacing are optimised for different powder materials to obtain a continuous track profile and improved scanning efficiency. The interface between two layers of different materials is visualised by simulation; it was found that the phase migration at the interface is related to the convection flow inside the molten pool, which contributes to the mixing of the two materials and elemental diffusion. This study significantly contributes to the challenging area of multi-material additive manufacturing by providing a greater in-depth understanding of the SLM process from multi-material powder deposition to laser interaction with powders across multiple scanning tracks and different building layers than can be achieved by experimentation alone.
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.
Summary
The combined finite‐discrete element method (FDEM) was originally developed for fracture and fragmentation of brittle materials, more specifically for cementitious and rock‐like materials. In ...this work, a combination of a discrete crack and plastic deformation has been combined and applied to FDEM simulation of fracture. The deformation is described using a FDEM‐specific mechanistic approach with plastic deformation being formulated in material embedded coordinate systems leading to multiplicative decomposition and plastic flow, that is, resolved in stretch space; this is combined with the FDEM fracture and fragmentation criteria. The result and main novelty of the present work is a robust framework for simulation of large strain solid deformation combined with a multiplicative decomposition‐based model that simultaneously involves elasticity, plasticity, and fracture.
Based on a combined finite-discrete element method (FDEM), this study builds a coupled thermo-mechanical model (termed FDEM-TM) to simulate thermal cracking of rock. The coupled thermo-mechanical ...model consists of two major parts. In the first part the temperature distribution of the system is analyzed based on the heat conduction equation. In the second part the thermal stress caused by temperature change is added to perform mechanical fracture calculation. Based on these two parts, we can model rock fracture driven by thermo-mechanical coupling. Three examples with analytic solutions are used to verify the correctness of the model in dealing with the problems of steady-state heat conduction, unsteady-state conduction, and thermal-mechanic coupling, respectively. In addition, an example of thermal cracking is also given and compared with the experimental results. The simulation results are in excellent agreement with the analytical solutions or experimental results, verifying the correctness of the coupled thermo-mechanical model to simulate thermal cracking. The proposed method provides a new tool for thermal-mechanical coupling problems in geothermal exploitation.
•An FDEM-based model, FDEM-TM, is proposed for simulating thermal cracking of rock.•FDEM-TM is able to describe thermal-mechanic behaviors in rock.•The proposed method provides a new way for thermal-mechanical coupling problems in geothermal exploitation.