Abstract

One of the main drawbacks of all the time integration algorithms using an Eulerian formulations in Coupled Problems is the restricted time-step to be used to have acceptable results.

For the case of fluid-structure interactions (FSI) with or without free-surfaces or for the case of fluid with moving internal interfaces (multi-fluids), it is well known that in the explicit integrations, the time-step to be used in the solution is stable only for time-step smaller than two critical values: the Courant-Friedrichs-Lewy (CFL) number and the Fourier number. The first one is concerning with the convective terms and the second one with the diffusive ones. Both numbers must be less than one to have stable algorithms. For convection dominant problems the condition CFL<1 becomes crucial and limit the use of explicit methods or outdistance its to be efficient. On the other hand, implicit integrations using Eulerian formulations are restricted in the time-step size due to the lack of convergence of the non-linear terms. Both time integrations, explicit or implicit are, in most cases, limited to CFL no much larger than one.

In this lecture we will present a Particle Method to solve coupled problems like FSI or multi-fluid problems that use in all the domain (solid and fluid) a Lagrangian formulation with explicit or implicit time integration without the CFL<1 restriction. This allows large time-steps, independent of the spatial discretization, having equal or better precision that an Eulerian integration.

The proposal will be tested numerically for FSI and multi-fluid flows problems using the Particle Finite Element Method second generation (PFEM2). The results show than this Particle Method is largely more efficient compared as well in accuracy as in computing time with other more standard Eulerian formulations.

Abstract

Lagrangian finite element methods emerged in fluid dynamics when the deficiencies of the Eulerian
methods in treating free surface flows (or generally domains undergoing large shape deformations)
were faced. Their advantage relies upon natural tracking of boundaries and interfaces, a feature
particularly important for interaction problems. Another attractive feature is the absence of the
convective term in the fluid momentum equations written in the Lagrangian framework resulting
in a symmetric discrete system matrix, an important feature in case iterative solvers are utilized.
Unfortunately, the lack of the control over the mesh distortions is a major drawback of Lagrangian
methods. In order to overcome this, a Lagrangian method must be equipped with an efficient
re-meshing tool.

This work aims at developing formulations and algorithms where maximum advantage of using
Lagrangian finite element fluid formulations can be taken. In particular we concentrate our attention
at fluid-structure interaction and thermally coupled applications, most of which originate from
practical “real-life” problems. Two fundamental options are investigated - coupling two Lagrangian
formulations (e.g. Lagrangian fluid and Lagrangian structure) and coupling the Lagrangian and
Eulerian fluid formulations.
In the first part of this work the basic concepts of the Lagrangian fluids, the so-called Particle
Finite Element Method (PFEM) [1], [2] are presented. These include nodal variable storage, mesh
re-construction using Delaunay triangulation/tetrahedralization and alpha shape-based method for
identification of the computational domain boundaries. This shall serve as a general basis for all the
further developments of this work.

Next we show how an incompressible Lagrangian fluid can be used in a partitioned fluid-structure
interaction context. We present an improved Dirichlet-Neumann strategy for coupling the incompressible
Lagrangian fluid with a rigid body. This is finally applied to an industrial problem dealing
with the sea-landing of a satellite capsule.

In the following, an extension of the method is proposed to allow dealing with fluid-structure
problems involving general flexible structures. The method developed takes advantage of the symmetry
of the discrete system matrix and by introducing a slight fluid compressibility allows to treat
the fluid-structure interaction problem efficiently in a monolithic way. Thus, maximum benefit from
using a similar description for both the fluid (updated Lagrangian) and the solid (total Lagrangian)
is taken. We show next that the developed monolithic approach is particularly useful for modeling
the interaction with light-weight structures. The validation of the method is done by means of comparison with experimental results and with a number of different methods found in literature.

The second part of this work aims at coupling Lagrangian and Eulerian fluid formulations. The
application area is the modeling of polymers under fire conditions. This kind of problem consists
of modeling the two subsystems (namely the polymer and the surrounding air) and their thermomechanical
interaction. A compressible fluid formulation based on the Eulerian description is used for
modeling the air, whereas a Lagrangian description is used for the polymer. For the surrounding air
we develop a model based upon the compressible Navier-Stokes equations. Such choice is dictated by
the presence of high temperature gradients in the problem of interest, which precludes the utilization
of the Boussinesq approximation. The formulation is restricted to the sub-sonic flow regime, meeting
the requirement of the problem of interest. The mechanical interaction of the subsystems is modeled
by means of a one-way coupling, where the polymer velocities are imposed on the interface elements
of the Eulerian mesh in a weak way. Thermal interaction is treated by means of the energy equation
solved on the Eulerian mesh, containing thermal properties of both the subsystems, namely air and
polymer. The developments of the second part of this work do not pretend to be by any means
exhaustive; for instance, radiation and chemical reaction phenomena are not considered. Rather we
make the first step in the direction of modeling the complicated thermo-mechanical problem and
provide a general framework that in the future can be enriched with a more detailed and sophisticated
models. However this would affect only the individual modules, preserving the overall architecture
of the solution procedure unchanged.

Each chapter concludes with the example section that includes both the validation tests and/or
applications to the real-life problems. The final chapter highlights the achievements of the work and
defines the future lines of research that naturally evolve from the results of this work.

Abstract

En algunos problemas relevantes de la ingeniería el campo mecánico están fuertemente acoplados con el de temperatura y el de transporte de masa. La solución de estos problemas es compleja, sobre todo, cuando las deformaciones son grandes. Por ello, a menudo se acude a aproximaciones numéricas que pueden encontrar todos los campos involucrados en el problema acoplado. En este artículo se propone una discretización por el método de los elementos finitos del problema de tres campos en el rango no lineal de deformación. La formulación tiene en cuenta todos los acoplamientos posibles y se proporcionan todos los detalles necesarios para su completa implementación, llenando un hueco de la literatura de estos métodos. Además, se describe la fundamentación termodinámica de la formulación, con lo que puede emplearse como base para simular procesos de interés como el comportamiento de baterías, la fragilización de metales por hidrógeno, etc.

Abstract