One main concern of this work is to develop an efficient particle-tracking-managing algorithm in the framework of a hybrid pressure-based finite-volume/probability-density-function (FV/PDF) Monte-Carlo (MC) solution algorithm to extend the application of FV/PDF MC methods to absolutely incompressible flows and speedup the convergence rate of solving the fluctuating velocity-turbulent frequency joint PDF equation in turbulent flow simulations. Contrary to the density-based algorithms, the pressure-based algorithms have stable convergence rates even in zero-Mach number flows. As another contribution, literature shows that the past developed methods mostly used mesh searching techniques to attribute particles to cells at the beginning of each tracking time-step. Also, they had to calculate the linear basis functions at every time-step to estimate the particle mean fields and interpolate the data. These calculations would be computationally very expensive, time-consuming, and inefficient in computational domains with arbitrary-shaped 3D meshes. As known, the barycentric tracking is a continuous particle tracking method, which provides more efficiency in case of handling 3D domains with general mesh shapes. The barycentric tracking eliminates any mesh searching technique and readily provides the convenient linear basis functions. So, this work benefits from these advantages and tracks the particles based on their barycentric coordinates. It leads to less computational work and a better efficiency for the present method. A bluff-body turbulent flow case is examined to validate the present FV/PDF MC method. From the accuracy perspective, it is shown that the results of the present algorithm are in great agreement with experimental data and available numerical solutions. The present study shows that the number of particle time-steps required to reach the statistically steady-state condition is at least one-sixth less than the previously developed algorithms. This also approves a faster convergence rate for the present hybrid pressure-based algorithm.
Landslides, which are the sources of most catastrophic natural disasters, can be subaerial (dry), submerged (underwater), or semi-submerged (transitional). Semi-submerged or transitional landslides occur when a subaerial landslide enters water and turns to submerged condition. Predicting the behavior of such a highly dynamic multi-phase granular flow system is challenging, mainly due to the water entry effects, such as wave impact and partial saturation (and resulted cohesion). The mesh-free particle methods, such as the moving particle semi-implicit (MPS) method, have proven their capabilities for the simulation of the highly dynamic multiphase systems. This study develops and evaluates a numerical model, based on the MPS particle method in combination with the μ(I) rheological model, to simulate the morphodynamic of the granular mass in semi-submerged landslides in two and three dimensions. An algorithm is developed to consider partial saturation (and resulting cohesion) during the water entry. Comparing the numerical results with the experimental measurements shows the ability of the proposed model to accurately reproduce the morphological evolution of the granular mass, especially at the moment of water entry.
In simulations using the particle finite element method (PFEM) with node-based strain smoothing technique (NS-PFEM) to simulate the incompressible flow, spatial and temporal instabilities have been identified as crucial problems. Accordingly, this study presents a stabilized NS-PFEM-FIC formulation to simulate an incompressible fluid with free-surface flow. In the proposed approach, (1) stabilization is achieved by implementing the gradient strain field in place of the constant strain field over the smoothing domains, handling spatial and temporal instabilities in direct nodal integration; (2) the finite increment calculus (FIC) stabilization terms are added using nodal integration, and a three-step fractional step method is adopted to update pressures and velocities; and (3) a novel slip boundary with the predictor–corrector algorithm is developed to deal with the interaction between the free-surface flow with rigid walls, avoiding the pressure concentration induced by standard no-slip condition. The proposed stabilized NS-PFEM-FIC is validated via several classical numerical cases (hydrostatic test, water jet impinging, water dam break, and water dam break on a rigid obstacle). Comparisons of all simulations to the experimental results and other numerical solutions reveal good agreement, demonstrating the strong ability of the proposed stabilized NS-PFEM-FIC to solve incompressible free-surface flow with high accuracy and promising application prospects.
With the assistance of the moving least-squares (MLS) interpolation functions, a two-dimensional finite element code is developed to consider the effects of a stationary or moving solid body in a flow domain. At the same time, the mesh or grid is independent of the shape of the solid body. We achieve this goal in two steps. In the first step, we use MLS interpolants to enhance the pressure (P) and velocity (V) shape functions. By this means, we capture different discontinuities in a flow domain. In our previous publications, we have named this technique the PVMLS method (pressure and velocity shape functions enhanced by the MLS interpolants) and described it thoroughly. In the second step, we modify the PVMLS method (the M-PVMLS method) to consider the effect of a solid part(s) in a flow domain. To evaluate the new method's performance, we compare the results of the M-PVMLS method with a finite element code that uses boundary-fitted meshes.