Research
Ongoing Projects
Dynamics and deposition morphology of granular column collapses
We investigated the collapse behavior of granular columns to better understand the dynamics and rheology of granular materials. In the preliminary study, we studied the collapse of dry granular columns, which results in various flow phenomena and deposition morphologies, which we link to the inter-granular and particle-boundary frictional coefficients and the initial aspect ratio of the granular columns. With the guidance of experiments, computational studies with the discrete element method (DEM) were performed to investigate the universality of dry granular column collapses. The resulting morphology of a collapsed granular column can be linked to an effective initial aspect ratio, and be described by different collapsing regimes, which are controlled by both the effective aspect ratio and the particle-boundary friction. The finding is important for understanding how granular materials, such as debris flows and landslides, behave and how they deposit when the granular flow halts.link
Granular collapse patterns observed in both DEM simulation and experiments (left). DEM simulation (right).
Size and Cross-section shape effect on granular column collapses
Based on previous works, we further aim to understand the size effect in granular column collapses. Geophysical hazards, such as landslides, debris flows, pyroclastic flows, and rock avalanches, occur frequently worldwide, and it is difficult to predict their initiation and dynamical behavior. Many geophysical flows can be considered as granular materials, where interaction between particles matters. However, it is difficult to upscale a grain-scale particle interaction to a large-scale predicting tool to capture the behavior of many complex granular flows. This paper starts from the collapse of granular columns, a fundamental problem often used to investigate the mobility of granular materials, and carefully navigates to examine the run-out behaviors of columns with different sizes, so that we can obtain a better picture of how different system sizes lead to different collapsing scenarios. We analyzed the characteristic length scale associated with the collapse of granular columns, and the strong force network presenting at the beginning of the collapse, to understand better what drives the size effect of granular systems, and what we can gain from this research to predict the behavior of geophysical flows.link
Response mechanism of sediments adsorption of pollutants to hydrodynamic conditions
Sediment in rivers and lakes play an important role in pollutant transport since pollutant can be absorbed by sediment or be released to water from sediment surface. Increasing fluid velocity can enhance or hinder adsorption which roughly depends on whether sediments suspend or not. The influence of hydrodynamic conditions on adsorption process is still not well understand. One challenging issue on modelling this process is: in conventional methods, the total amount of adsorbed pollutants is at the same magnitude (even lower) as numerical errors due to the moving particles. we introduced the Random Walking Method to solve pollutants transport with exact mass conservation. This improvement allows us to study the adsorption process at grain-scale. link
Solute transport in the fluid-particle system, particles are under periodic body force. Initial setup (left), Solute concentration porfile at different time step (middle), Dimensionless dispersion coefficient as a funtion of volume fraction (right).
Debris flow modelling
Debris flow is a type of mass movement that is generated under certain mechanical conditions of the soil and external factors as weather, seismic and anthropic activity. The water-sediment interaction processes dominate the behaviour of the flow while is driven by gravity. The sediments size have a broad spectrum from fine grains to big boulders (up to ~11 m in diameter). The dynamic of the interaction of these big boulders with the entire flow is not well known. And, reliable predictions of the dynamics of these phenomena is crucial for risk management.Thus, Smooth Particles Hydrodynamics (SPH) coupled to Discrete Element Method (DEM), are employed in this project to model debris flows link.
SPH to represent water and soil components, whereas DEM is used to include large size boulders).
Fluid–solid–particle systems such as rapid mixed flow are very common in steep mountainous regions and have the potential to damage infrastructures. Baffles are typically used as flow-impeding structures to mitigate the destructive power of this natural hazard. However, modelling these complex interactions is still very challenging due to their multi-phase and multi-scale nature. Therefore, we have developed a hybrid model that combines the advantages of the two-phase (solid and liquid) Material Point Method (MPM) on handling continuum materials (fluids and soils) with large deformations and the capability of the Discrete Element Method (DEM) on simulating the mechanical behaviours of rigid particles (rocks). link.
Invariant forms for the permeability and transport properties of discrete fracture networks
Translating fracture statistics to global hydraulic properties in subsurface fractured rocks is a complex but appealing task, as measuring parameters such as permeability can be prohibitively expensive. To tackle this challenge, we developed the finite-size scaling (FSS) hypothesis, drawing inspiration from percolation theory in statistical physics. We created numerous Monte Carlo iterations using our CUDA-based code, cuDFNsys, to simulate flow and transport in discrete fracture networks (DFNs), as shown in the figure below. We generated a wide range of DFN scenarios with vastly different input parameters to ensure generality. We then nondimensionalized the obtained connectivity, permeability, density of fractures, domain sizes, and so on. By analyzing the dimensionless quantities, we found that the FSS hypothesis can accurately predict connectivity and permeability as an invariant function of dimensionless density and domain sizes with several universal critical quantities. Furthermore, the FSS function can identify the transition point of dimensionless density, where permeability remains constant regardless of domain sizes. Around this point, the scale-dependence of permeability changes from negative to positive. Our findings provide a strong theoretical foundation for understanding the relationship between fracture attributes and field-scale hydraulic properties. This research inspires further investigation into applying the FSS function at the field scale, which will improve the information that earth scientists can obtain from fracture statistics link.
Realistic shapes for fluid immersed particle simulations
We developed a Metaball-Imaging (MI) based Discrete Element Lattice Boltzmann Method (DELBM) for simulations of fluid-particle systems with irregular-shaped particles. The MI algorithm was used to capture the real grain shape for simulations with Metaballs, which are coupled in the DELBM framework. The core challenge, contact detection, was tackled robustly by gradient calculation of the closest point with a Newton-Raphson based scheme. The proposed model is rigorously validated with multiple experiments, including analytical solutions and cases of irregular-shaped granular media immersed in fluid systems. Results show that the proposed model is effective and powerful in probing the micromechanics of these systems. Furthermore, the model's potential for studying shape-induced physical processes is investigated, demonstrating its usefulness in fluid-particle system research. The study's main contribution is its ability to simulate complex-shaped particles, which is a difficult task in the field of fluid-particle system simulations.link.
Mass transfer around a rising bubble
In this project, we focused on the mass transfer process around a rising bubble, utilizing a pseudopotential Lattice Boltzmann Model (LBM) simulation. In contrast to previous studies that assumed bubble rising through an effective buoyancy force, we developed an LBM that accounts for the actual buoyancy force stemming from the pressure gradient. By applying thermodynamic analysis to the multi-component pseudopotential LBM, we established an initial (mechanical and chemical potential) equilibrium state and interpretated the Henry law from the viewpoint of chemical potential balance. We implemented our LBM code to model bubble rising in both saturated and unsaturated solutions. Our results unveiled a significant correlation between the diffusion (Figure 1) and advection processes (Figure 2). Specifically, at the lateral side of the rising bubble, the strong advection process continually refreshes the bubble surface, creating a large concentration gradient and, consequently, a strong diffusion process. Conversely, at the bottom side of the bubble, the advection is relatively weak or even negligible, leading to the accumulation of saturated solution in that area and, therefore, weak or negligible diffusion process. Additionally, our LBM code satisfactorily simulates various bubble shapes and their respective distribution of diffusion zones, as demonstrated in Figure 3. We observed a critical point on the bubble surface where the radius of curvature reaches a minimum. This point effectively separates the bubble surface into two parts, with one side experiencing a strong diffusion process and the other side exhibiting negligible diffusion. These findings provide valuable insights into the mass transfer mechanism around a rising bubble and are expected to enhance the penetration theory to account for the effects of bubble shapes and wakes. This work has been published on the journal of Physics of Fluids.link
Submmerged granular column collapses
The occurrence of submarine landslides can have a significant impact on our lives, both directly and indirectly. This large-scale natural phenomenon involves multiple scales and phases and can result in the transportation of large amounts of sediment over long distances, causing damage to submerged structures such as pipelines, and offshore oil wells, and even leading to tsunamis. However, investigating the initiation, propagation, and deposition of submarine landslides directly can be a challenging task.
To address this issue, we developed an LBM-DEM numerical model to replicate submerged granular column collapses at a laboratory scale and examine their underlying mechanisms. The insights gained from this study can be applied to predict macroscopic submarine landslides more accurately.
Molecular dynamic simulations of flow inside carbon nanotubes (CNT)
We investigate nano-mechanical transport mechanisms and intermolecular structure of nanofluids to better understand and develop water desalination based nanomaterials. In a preliminary study, we have investigated the flow enhancement of nanoconfined water in the carbon nanotube (CNT), which is contributed by the various scale and different force fields. We connect the above contributions to various intermolecular structures and different fluid states. According to previous experiments and simulations, we have studied the flow enhancement inconsistencies through molecular dynamics (MD) simulations. MD simulations indicate that the inconsistency of pressure-dependent flow enhancements can be attributed to the liquid structure, and quantitatively analyzed the structural properties through calculate the apparent friction coefficient, which are affected by the external force field. This finding will help us understand the flow enhancement inconsistencies and improve the water permeability efficiency of desalination.
