Julio Marti

A fast and accurate method to solve the incompressible Navier‐Stokes equations

Purpose – The purpose of this paper is to highlight the possibilities of a novel Lagrangian formulation in dealing with the solution of the incompressible Navier‐Stokes equations with very large time steps.Design/methodology/approach – The design of the paper is based on introducing the origin of this novel numerical method, originally inspired on the Particle Finite Element Method (PFEM), summarizing the previously published theory in its moving mesh version. Afterwards its extension to fixed mesh version is introduced, showing some details about the implementation.Findings – The authors have found that even though this method was originally designed to deal with heterogeneous or free‐surface flows, it can be competitive with Eulerian alternatives, even in their range of optimal application in terms of accuracy, with an interesting robustness allowing to use large time steps in a stable way.Originality/value – With this objective in mind, the authors have chosen a number of benchmark examples and have pr...

Advances in the Particle Finite Element Method (PFEM) for Solving Coupled Problems in Engineering

We present some developments in the formulation of the Particle Finite Element Method (PFEM) for analysis of complex coupled problems on fluid and solid mechanics in engineering accounting for fluid-structure interaction and coupled thermal effects, material degradation and surface wear. The PFEM uses an updated Lagrangian description to model the motion of nodes (particles) in both the fluid and the structure domains. Nodes are viewed as material points which can freely move and even separate from the main analysis domain representing, for instance, the effect of water drops. A mesh connects the nodes defining the discretized domain where the governing equations are solved, as in the standard FEM. The necessary stabilization for dealing with the incompressibility of the fluid is introduced via the finite calculus (FIC) method. An incremental iterative scheme for the solution of the non linear transient coupled fluid-structure problem is described. The procedure for modelling frictional contact conditions at fluid-solid and solidsolid interfaces via mesh generation are described. A simple algorithm to treat soil erosion in fluid beds is presented. An straight forward extension of the PFEM to model excavation processes and wear of rock cutting tools is described. Examples of application of the PFEM to solve a wide number of coupled problems in engineering such as the effect of large waves on breakwaters and bridges, the large motions of floating and submerged bodies, bed erosion in open channel flows, the wear of rock cutting tools during excavation and tunneling and the melting, dripping and burning of polymers in fire situations are presented.

An improved enrichment method for weak discontinuities for thermal problems

Purpose The purpose of this paper is to propose a new elemental enrichment technique to improve the accuracy of the simulations of thermal problems containing weak discontinuities. Design/methodology/approach The enrichment is introduced in the elements cut by the materials interface by means of adding additional shape functions. The weak form of the problem is obtained using Galerkin approach and subsequently integrating the diffusion term by parts. To enforce the continuity of the fluxes in the “cut” elements, a contour integral must be added. These contour integrals named here the “inter-elemental heat fluxes” are usually neglected in the existing enrichment approaches. The proposed approach takes these fluxes into account. Findings It has been shown that the inter-elemental heat fluxes cannot be generally neglected and must be included. The corresponding method can be easily implemented in any existing finite element method (FEM) code, as the new degrees of freedom corresponding to the enrichment are local to the elements. This allows for their static condensation, thus not affecting the size and structure of the global system of governing equations. The resulting elements have exactly the same number of unknowns as the non-enriched finite element (FE). Originality/value It is the first work where the necessity of including inter-elemental heat fluxes has been demonstrated. Moreover, numerical tests solved have proven the importance of these findings. It has been shown that the proposed enrichment leads to an improved accuracy in comparison with the former approaches where inter-elemental heat fluxes were neglected.