Research Thrusts

Additive Manufacturing

Metal SLM
SLM: Semi-coupled fully resolved CFD with DEM modeling of metallic powder based selective laser melting

Metal Selective Laser Melting

Our lab has gone for metallic poweder based selective laser melting (PB-SLM) for additive manufacturing. We have developed a physically based, high-fidelity approach based on coupling of DEM with fully resolved CFD to tackle the complicated thermal induced phase transition of granular powders into metal liquid and thereafter solidified product as well as the accompanied fluid-particle interactions that may control the finishing quality. Read our recent paper published in CMAME.

Advanced Particulate Modeling

SudoDEM for non-spherical particle modeling

A new generation advanced DEM code, SudoDEM, has been developed Dr Shiwei Zhao with full function for both 2D and 3D non-spherical particles modeling. It offers a rich library of prime particle shapes, including poly-superellipsoids, superellipsoids, cylinders, cones, polyhedrons for 3D and disks and superellipses for 2D, and is further empowered by OpenMP acceleration. It features effecient contact detection algoriths to handle contacts among complex-shaped particles.

We have elected to have it open-sourced for the community to free download and use: Full features of SudoDEM are summarised in a new paper published in Computer Physics Communications.

SudoDEM: rigorous and efficient DEM simulation of non-spherical granular particles
1D consolidation grain crushing
Grain crushing: A coupled Peridynamics with physics engine on continuous grain crushing simulation

Revolutionising Simulation on Grain Crushing

Grain crushing underpins the operation and performance of applications pertaining to multiple engineering and industrial disciplines. Our Lab has developed a whole new approach to simulate grain crushing, based on hybrid Peridynamics and Physics Engine which provides more objective (crushing criteria) and realistical (shape) consideration of major aspects of the physical process of grain crushing in a context of granular packing. We use Peridynamics to treat the crushing of each individual grain, and allow Physics Engine to effectively handle irregularly shaped grains in an assembly during the crushing process (see animation of continuous crushing process in 1D compression), where the continuous evolution and interplay between grain shape and grain size apparently reflect more of the physical nature of this complicated process.

Read our recent paper published in CMAME and an earlier paper appearing in Géotechnique for details.

Rock blasting induced fractures by Peridynamics

Novel Peridynamic models have been further developed to simulate dynamic fractures caused by rock blasting. In the study we developed a new non-ordinary state-based peridynamics approach to consider pressure dependency, strain rate effect, and viscoplasticity in rock fracture under blasting, in conjunction with updated Lagrangian peridynamics with Jones-Wilkins-Lee equation of state to simulate detonation of explosives and gas expansion. For detail please reader our paper published in Journal of the Mechanics and Physics of Solids (see also PDF).

rock blasting
Peridynamic simulation of rock blasting induced fratures in rocks.

Debris Flow Mitigation by Flexible Barriers

Collaborated with government sectors, we have been engaged in developing innovative computational tools to simulate debris flow mitigation by flexible ring barrier. We coupled computational fluid dynamics (CFD) with discrete element method (DEM) to capture the multi-phase, multi-way, multi-physics interactions among the debris solids, the debris fluid and the different deformable components of the flexible ring barrier are captured, including ring-ring and ring-cable frictional contact and sliding. Read our latest paper on this subject.

debris mixture impacts on flexible ring barrier CFD-DEM simulation of debris-fluid mixture impacton a flexible barrie
Rockfall from top Rockfall from side
Calibration of CFD-DEM simulation of flexible barrier: rockfall test

We are able to have the same ring flexible barrier subjected to the impact of different geophysical flows which may occur in a natural setting, ranging from dry rock avalanche, debris avalanche, debris flow, debris flood to mud flow, to investigate the impact mechanisms. We found that the transition from pileup to run up mechanisms is found closely related to flow features (including Froude number, solid volume concentration and fluid rheology) and barrier responses (barrier components including barrier net, cables and brakes as well as their coordinated interactions in response to the impact). The insights on the two impact mechanisms and their transitions gained in the study may provide a useful reference for future design of flexible barriers in mitigating hazardous geophysical flows. For detail please reader our paper published in Engineering Geology (see also PDF).

Hydrodynamic dead zone
Pileup mechanism during the impact of debris flow onto a ring flexible barrier

Multiscale Modeling in Geomechanics

soil pile interaciton
Coupled MPM-DEM simulation of large deformation in seabed soil driven by pipeline drifting

Large Deformation in Geomechanics

Our Lab has developed an innovative computational multiscale modeling framework, based upon FEM/DEM coupling or MPM/DEM coupling, to effectively tackle a wide range of geomechanics boundary value problems to offer cross-scale links between the macroscopic phenomena (such as large deformation, strain localisation, and failures) and their underlying microstructural origins and mechanisms. Large deformation failures in offshore geotechnics are a recent focus of our research (see a MPM/DEM simulation of soil-pipeline interaction in a typical seabed). Read our recent paper on this topic.

soil pile interaction equivalent shear strain A pipeline settles down before drifting laterally by wave and pushes large deformation in seabed soil

Large deformation in offshore anchor pullout

Multiscale modeling of offshore anchor pullout has been performed to predict the pullout bearing capacity of both horizontally and vertically placed anchor buried in loose and dense sand at different depth. Shown in the following animation are contact force network evolutions for three RVEs (from top to bottom) located behind and immediately in front of the anchor plate and at the edge of the front shear banding, respectively. For detail of the study refer to our paper published in JGGE (see also PDF). The work was a collaboration with Prof. Kenichi Soga of UC Berkeley.

anchor pullout Vertically placed plate anchor pulled horizontally to the right: evoluation of RVE force network behind, immediate in front and at the bottow right of the anchor.
anchor pullout
Coupled MPM-DEM simulation of large deformation around an anchor during pullout: variation of deformation pattern with anchor burial depth (H/B).
cyclic thermal mechanical loading
Multiscale modeling of biaxial shearing of granular media subjected to boundary cyclic thermal loads

Multiscale modeling of thermo-mechanical coupling

We have further developed a novel multiscale modeling framework to simulate coupled thermo-mechanical behavior in granular media. We employ dual FEMs, one for mechanical and the other for thermal, which share common Gauss points. At each shared Gauss point, we embed a DEM assembly that considers and solve for heat conduction and force transfer at grain contacts to return required variables for the global FEM solutions. Illustrated are animations of the thermal and displacement fields for a biaxially sheared sample subjected to cyclic thermal loading.Read our paper published in Computer Methods in Applied Mechanics and Engineering.

Thermal multiscael modeling Hierarchical multiscale modeling of coupled thermo-mechanical behavior in granular media

Compaction Front in Granular Rocks

We have further extended the multiscale modeling framework to simulating granular rocks such as sandstone and limestone. We replicated the compaction band in a high-porosity limestone (Tuffeau de Maastricht) observed in X-ray tomography based on coupled 3D FEM/DEM simulations, We confirmed that significant debonding accompanied by collapse of macropores and grain rearrangements may constitute the major microstructural mechanisms causing the formation of compaction bands in high-porosity Tuffeau de Maastricht (Read paper).

We have further performed coupled hydromechanical simulation to show that excess pore pressure can result in transitional evolution of deformation band patterns in saturated sandstone, e.g., from compaction band to shear band as shown above. Read our paper published in Engineering Geology.

Compaction band and flow flux Multiscale modeling of hydro-mechanical coupling during the formation of compaction bands in high-porosity sandsone
Compaction front in high-porosity limestone Compaction front in Tuffeau de Maastricht
Herarchnical multiscale modeling of evolaution of compaction front in Tuffeau de Maastricht (a porous limestone) in comparison with X-ray CT data

Constitutive Modeling of Granular Sand

(To be updated) We focus on anisotropic critical state of granular soils and develop advanced constitutive models centered on the use of fabric tensor and fabric evaluation to describe both quasi-static and cyclic behavior of granular sand.