COMPLAS 2021 is the 16th conference of the COMPLAS Series.
The COMPLAS conferences started in 1987 and since then have become established events in the field of computational plasticity and related topics. The first fifteen conferences in the COMPLAS series were all held in the city of Barcelona (Spain) and were very successful from the scientific, engineering and social points of view. We intend to make the 16th edition of the conferenceanother successful edition of the COMPLAS meetings.
The objectives of COMPLAS 2021 are to address both the theoretical bases for the solution of nonlinear solid mechanics problems, involving plasticity and other material nonlinearities, and the numerical algorithms necessary for efficient and robust computer implementation. COMPLAS 2021 aims to act as a forum for practitioners in the nonlinear structural mechanics field to discuss recent advances and identify future research directions.
Scope
COMPLAS 2021 is the 16th conference of the COMPLAS Series.
The development of more accurate force prediction models developed through particle-level experiments is required to accurately model non-dilative interfaces from micro to macro. Further, selecting reliable input parameters for DEM remains a challenge. Thus, micromechanical experimental studies are of fundamental importance that can provide insights into microscale aspects for in-depth knowledge of non-dilative interfaces. This study presents custom-designed, reliable, and sensitive equipment that facilitates shear tests for non-dilative interfaces under different configurations that simulate suitable conditions for geotechnical applications. This research offers a logical rationale for the non-dilative interface system's observed shear behavior.
Abstract The development of more accurate force prediction models developed through particle-level experiments is required to accurately model non-dilative interfaces from micro to [...]
Currently existing computational fluid dynamics-discrete element method (CFD-DEM) solvers suffer from computationally expensive coupling between the CFD and DEM as it requires calculating at each fluid time-step the void fraction and the solid-fluid forces such as drag, lift, buoyancy, and undisturbed flow forces. We develop a unified finite element CFD-DEM solver which integrates the CFD and DEM solvers into a single software resulting in faster and cheaper coupling. Our fluid formulation is stabilized using tailored techniques to prevent oscillations in regions of sharp gradients, to enhance the robustness of the formulation and local mass conservation, and to relax the Ladyzhenskaya-Babuska-Brezzi inf-sup condition. The developed solver supports high order finite elements resulting in better accuracy with larger cell sizes. Moreover, our solver supports dynamic load balance parallelization for both the particles and the fluid. This evens the distribution of workloads among processors, resulting in better efficiency and resource exploitation. Additionally, we develop a new spatially and temporally continuous analytical void fraction scheme called Quadrature-Centered Method (QCM). This scheme results in less computational time, better accuracy and convergence, and enhanced mass conservation. It also enables the use of very small CFD time-steps thus achieving better temporal accuracy and the use of mesh sizes smaller than those commonly used in CFD-DEM (< 3 times the particle diameter.) We validate our solver through several cases among which we will discuss a spouted bed test case where particle velocities matched those of the experiments, a particle sedimentation test case where we study the effect of the void fraction scheme choice on the particle’s terminal velocity, and a particle Rayleigh-Taylor Instability where the particles constituted the heavy phase and where we study the evolution of the mixing layer with time.
Abstract Currently existing computational fluid dynamics-discrete element method (CFD-DEM) solvers suffer from computationally expensive coupling between the CFD and DEM as it requires [...]
The hydrodynamic interactions between particles have significant effects in many engineering fields, such as fluidized beds and slurry sedimentation. This is because they can impact the macroscopic process parameters of these systems (e.g., particles' cluster sedimentation speed) [1-2]. To better understand the hydrodynamic interaction between particles and its effects on macroscopic process parameters, it is necessary to understand the hydrodynamic force interactions between individual particles as a function of their relative position and velocity. The first step in this direction is the interaction between a pair of particles, which remains an active area of research. It is established that drag and lift forces applied on a particle change in function of its relative position and velocity to another particle. However, the impact of these changes on their dynamics is limited. We first study the effects of the relative position between the particles and the Reynolds number on the drag and lift forces applied to them to constitute a pairwise fluid force model. This force model is the basis for a new reduced-order particle dynamics model that also includes lubrication, Basset, and added mass forces. We compare the results obtained for a series of sedimentation cases with the reduced-order model with those obtained from a resolved computational fluid dynamic solver coupled with a discrete element method (CFD-DEM) [3]. Comparing the fluid forces obtained between these to the model enables us to assess the impact of the particles' relative motion on the hydrodynamic forces applied to them (e.g., drag, lift, and forces). These sedimentation cases also allow us to evaluate the effects of the density ratio between the particle and the fluid on the virtual mass force and its impact on the dynamic of the particles. This study is a stepping stone toward a complete model for hydrodynamic forces in particle clusters.
Abstract The hydrodynamic interactions between particles have significant effects in many engineering fields, such as fluidized beds and slurry sedimentation. This is because they [...]
C. Kloss, A. Mayrhofer, M. Niemann, A. Aigner, P. Seil, M. Kwakkel, C. Govina
particles2023.
Abstract
Discrete Element Method (DEM) and DEM coupled to Computational Fluid Dynamics (CFD-DEM) are established techniques for optimization and design of particle processes. Its applicability to a wide range of processes has been proven for many different industrial and environmental applications. The extension to new fields and processes has been made possible by continuous improvements of: (i) models, (ii) numerical methods and (iii) computational performance. DEM and CFD-DEM have evolved from pure “particle modelling tools” to numerical tools to model “particulate flow””. Combining the Lagrangian nature of discretization with complex interaction models, the behaviour of viscous pastes, compressible powders, melting polymers just to name a few, has become feasible. Additionally much attention has been given to improvement of numerical aspects, which led to improved stability and therefore applicability of the models. Last but not least, the computational efficiency and possibility to make use of available computational resources has boosted the technology to new levels. The authors give their perspective on some corner-stones and highlights in modelling and development that were made in the past few years, which lead to this break-through and give some concrete examples of current state of the art modelling capability. Based on this solid foundation that has been build, new goals are within reach and the authors will give some insight on future opportunities for this modelling technology.
Abstract Discrete Element Method (DEM) and DEM coupled to Computational Fluid Dynamics (CFD-DEM) are established techniques for optimization and design of particle processes. Its applicability [...]
Interacting particle systems are ubiquitous in nature and engineering. Access to the governing particle interaction law is of fundamental importance for a complete understanding of such systems. However, it is particularly challenging to extract this information from experimental observations due to the intricate configuration complexities involved. Machine learning methods have the potential to learn the behavior of interacting particle systems by combining experiments with data analysis methods. However, most existing algorithms focus on learning the kinetics at the particle level and do not learn the pairwise interactions specifically. Moreover, in reality, interacting particle systems are often heterogeneous, where multiple interaction types coexist simultaneously and relational inference is required. An approach that can simultaneously reveal the hidden pairwise interaction types and infer the unknown heterogeneous governing interactions constitutes a necessary advancement for our understanding of particle systems. However, this task is considerably more challenging than its homogeneous counterpart. Here, we propose the physics-induced graph network for particle interaction (PIG'N'PI) allowing to precisely infer the pairwise interactions that are consistent with underlying physical laws by only being trained to predict the particle acceleration for homogeneous systems. We further propose a novel method for relational inference which combines probabilistic inference and PIG'N'PI to learn different kinds of interactions for heterogeneous systems. We test the proposed methodologies on multiple benchmark datasets and demonstrate that the learnt interactions are consistent with the underlying physics and the proposed relational inference method achieves superior performance in correctly inferring interaction types. In addition, the proposed model is data-efficient and generalizable to large systems when trained on small systems, contrary to previously proposed solutions. The developed methodology constitutes a key element for the discovery of the fundamental laws that determine macroscopic mechanical properties of particle systems.
Abstract Interacting particle systems are ubiquitous in nature and engineering. Access to the governing particle interaction law is of fundamental importance for a complete understanding [...]
This study focused on investigating the deformation behavior of non-spherical particles under high-load compaction, utilizing the multi-contact discrete element method (MC-DEM). To account for the non-spherical shape of the particles, two methods were employed: the bonded multi-sphere method (BMS) and the conventional multi-sphere (CMS). The BMS approach yielded accurate results in predicting the compression behavior of a single rubber sphere, while the CMS method failed to replicate the same behavior. Building on these findings, the BMS method was utilized to study the uniaxial compaction of Avicel® PH 200, a popular choice of excipient due to its ability to enhance the stability, flowability, and compressibility of tablet formulations. The results obtained from this study showed very good agreement with experimental data. To generate realistic 3D models of particles, a novel approach was introduced, which combines 2D projections and deep learning algorithms utilizing a 3D convolutional neural network (3D-CNN) methodology. Surrogate model was used to overcome the computational cost of DEM simulations. The findings of this study offer a valuable tool for researchers and engineers to efficiently and accurately generate 3D models of particles, leading to new insights and innovations in a range of applications such as rock and mineral analysis, battery materials, pharmaceuticals, and space exploration.
Abstract This study focused on investigating the deformation behavior of non-spherical particles under high-load compaction, utilizing the multi-contact discrete element method (MC-DEM). [...]
The Smoothed Particle Hydrodynamics (SPH) is a numerical scheme in which the domain is discretized into Lagrangian particles in the context of continuum mechanics. It has been widely used for fluid dynamics problems, and recently it has also been applied to solid mechanics with reasonable success. In this work, we present a total Lagrangian SPH (TL-SPH) for the application of solid mechanics and contact problems. Total Langrangian stands for the usage of the reference configuration to calculate the spatial derivatives. As a consequence, all particles maintain a perfect distribution, which, in turns, results in highly accurate calculations. In this way, the method is able to eliminate any problem related to tensile instability, which is one of the main shortcomings of SPH. Here, we introduce a simple, yet robust, way to include finite strain elastoplasticity into the TL-SPH method based on the logarithmic strain. Then, the elastic part can be easily defined with the Hencky elastic model, and the plastic part with any yield criteria such as Drucker-Prager. In addition, we develop a contact algorithm capable of simulating solid-solid contact problems. In this way, the simulation of different objects becomes a mix of continuous and discontinuous problems. Finally, we show the applicability of the method with several simple tests to validate the accuracy of the TL-SPH for simulating the elastoplastic solid material, as well as for contact problems including direct impact and friction effects. REFERENCES [1] Morikawa D.S. Toward robust landslide simulations from initiation to post-failure using the Smoothed Particle Hydrodynamics. Kyushu University PhD thesis (2022).
Abstract The Smoothed Particle Hydrodynamics (SPH) is a numerical scheme in which the domain is discretized into Lagrangian particles in the context of continuum mechanics. It has [...]
The Discrete Element Method (DEM) is a numerical approach that deals with the motion and interactions of individual elements. This method is mainly used in particle mechanics because it is overshadowed by other techniques such as the finite element method (FEM) when dealing with continuous problems. For this reason, it is not commonly known among structural engineers. However, its use can be found in cases that combine particle and continuum mechanics problems, such as crack propagation in reinforced concrete members. Calculating and optimizing these types of problems using FEM is challenging due to frequent mesh changes or the requirement for difficult-to-detect parameters in standard practice. This post discusses the possibility of using the DEM method with an extension of the beam-bound model (BBM). In this method, a beam element is inserted between the bounded discrete elements to transfer all types of forces and moments. The problem of this method is to define the correct cross-sectional and material characteristics of the individual beam members. In study, we focus on the determination of these parameters and the verification of this method on model examples of reinforced concrete elements such as simply supported and fixed beams.
Abstract The Discrete Element Method (DEM) is a numerical approach that deals with the motion and interactions of individual elements. This method is mainly used in particle mechanics [...]