Abstract

Particle methods are those in which the problem is represented by a discrete number of particles. Each particle moves accordingly with its own mass and the external/internal forces applied on it. In this paper the Particle Finite Element Method based on finite element shape functions is used to solve the continuous fluid mechanics equations in the case of heterogeneous density. To evaluate the external applied forces to each particle, the incompressible Navier–Stokes equations are solved at each time step using a Lagrangian formulation. All the information in the fluid is transmitted via the particles. All kinds of density heterogeneous fluids and multiphase flows with internal interfaces including or not free-surfaces, breaking waves and fluid separations may be easily solved with this methodology. 

Abstract

Heterogeneous incompressible fluid flows with jumps in the viscous properties are solved with the particle finite element method using continuous and discontinuous pressure fields. We show the importance of using discontinuous pressure fields to avoid errors in the incompressibility condition near the interface.

Abstract

The purpose of this paper is to put in evidence that the fractional‐step method (FSM) used to solve the incompressible transient Euler and Navier–Stokes equations for free‐surface flows has a problem inherent to the method that may produce unacceptable variations of the domain volume. A simple modification of the free‐surface boundary term is introduced in order to reduce considerably the volume loss and preserve the computational advantages of the FSM.

Abstract

Laplace formulations are weak formulations of the Navier–Stokes equations commonly used in computational fluid dynamics. In these schemes, the viscous terms are given as a function of the Laplace diffusion operator only. Despite their popularity, recently, it has been proven that they violate a fundamental principle of continuum mechanics, the principle of objectivity. It is remarkable that such flaw has not being noticed before, neither detected in numerical experiments. In this work, a series of objectivity tests have been designed with the purpose of revealing such problem in real numerical experiments. Through the tests it is shown how, for slip boundaries or free-surfaces, Laplace formulations generate non-physical solutions which widely depart from the real fluid dynamics. These tests can be easily reproduced, not requiring complex simulation tools. Furthermore, they can be used as benchmarks to check consistency of developed or commercial software. The article is closed with a discussion of the mathematical aspects involved, including the issues of boundary conditions and objectivity.

Abstract

This work proposes a fully Lagrangian formulation for the numerical modeling of free-surface particle-laden flows. The fluid phase is solved using the particle finite element method (PFEM), while the solid particles embedded in the fluid are modeled with the discrete element method (DEM). The coupling between the implicit PFEM and the explicit DEM is performed through a sub-stepping staggered scheme. This work only considers suspended spherical particles that are assumed not to affect the fluid motion. Several tests are presented to validate the formulation. The PFEM–DEM results show very good agreement with analytical solutions, laboratory tests and numerical results from alternative numerical methods.

Abstract

We present a Lagrangian nodal integration method for the simulation of Newtonian and non-Newtonian free-surface fluid flows. The proposed nodal Lagrangian method uses a finite element mesh to discretize the computational domain and to define the (linear) shape functions for the unknown nodal variables, as in the standard Particle Finite Element Method (PFEM). In this approach, however, the integrals are performed over nodal patches and not over elements, and strains/stresses are defined at nodes and not at Gauss points. This allows to limit the drawbacks associated with the remeshing and leads to a more accurate stress computation than in the classical elemental PFEM. Several numerical tests, in 2D and in 3D, are presented to validate the proposed nodal PFEM. In all cases, the method has shown a very good agreement with analytical solutions and with experimental and numerical results from the literature. A thorough comparison between nodal and elemental PFEMs is also presented, focusing on crucial issues, such as solution accuracy, convergence, mass conservation and sensitivity to mesh distortion.