Partitioned solution to fluid-structure interaction problems in application to free-surface flows

International audienceIn this work we discuss a way to compute the impact of free-surface flow on nonlinear structures. The approach chosen rely on a partitioned strategy that allows to solve strongly coupled fluid-structure interaction problem. It is then possible to re-use existing and validated strategy for each sub-problem. The structure is formulated in a Lagrangian way and solved by the finite element method. The free-surface flow approach considers a Volume-Of-Fluid (VOF) strategy formulated in an Arbitrary Lagrangian-Eulerian (ALE) framework, and the finite volume are used to discrete and solve this problem. The software coupling is ensured in an efficient way using the Communication Template Library (CTL). Numerical examples presented herein concern 2D validations case but also 3D problems with a large number of equations to be solved

Nonlinear fluid-structure interaction problem. Part II: space discretization, implementation aspects, nested parallelization and application examples

International audienceThe main focus of the present article is the development of a general solution framework for coupled and/or interaction multi-physics problems based upon re-using existing codes into software products. In particular, we discuss how to build this software tool for the case of fluid-structure interaction problem, from finite element code Feap for structural and finite volume code OpenFOAM for fluid mechanics. This is achieved by using the Component Template Library (CTL) to provide the coupling between the existing codes into a single software product. The present CTL code-coupling procedure accepts not only different discretization schemes, but different languages, with the solid component written in Fortran and fluid component written in \Cpp. Moreover, the resulting CTL-based code also accepts the nested parallelization. The proposed coupling strategy is detailed for explicit and implicit fixed-point iteration solver presented in the Part I of this paper, referred to Direct Force-Motion Transfer/Block-Gauss-Seidel. However, the proposed code-coupling framework can easily accommodate other solution schemes. The selected application examples are chosen to confirm the capability of the code-coupling strategy to provide a quick development of advanced computational tools for demanding practical problems, such as 3D fluid models with free-surface flows interacting with structures

Unified semi-analytical wall boundary conditions applied to 2-D incompressible SPH

International audienceThis work aims at improving the 2-D incompressible SPH model (ISPH) by adapting it to the unified semi-analytical wall boundary conditions proposed by Ferrand et al. [10]. The ISPH algorithm considered is as proposed by Lind et al. [25], based on the projection method with a divergence-free velocity field and using a stabilising procedure based on particle shifting. However, we consider an extension of this model to Reynolds-Averaged Navier-Stokes equations based on the k- turbulent closure model, as done in [10]. The discrete SPH operators are modified by the new description of the wall boundary conditions. In particular, a boundary term appears in the Laplacian operator, which makes it possible to accurately impose a von Neumann pressure wall boundary condition that corresponds to impermeability. The shifting and free-surface detection algorithms have also been adapted to the new boundary conditions. Moreover, a new way to compute the wall renormalisation factor in the frame of the unified semi-analytical boundary conditions is proposed in order to decrease the computational time. We present several verifications to the present approach, including a lid-driven cavity, a water column collapsing on a wedge and a periodic schematic fish-pass. Our results are compared to Finite Volumes methods, using Volume of Fluids in the case of free-surface flows. We briefly investigate the convergence of the method and prove its ability to model complex free-surface and turbulent flows. The results are generally improved when compared to a weakly compressible SPH model with the same boundary conditions, especially in terms of pressure prediction

Nonlinear fluid-structure interaction problem. Part I: implicit partitioned algorithm, nonlinear stability proof and validation examples

International audienceIn this work we consider the fluid-structure interaction in fully nonlinear setting, where different space discretization can be used. The model problem considers finite elements for structure and finite volume for fluid. The computations for such interaction problem are performed by implicit schemes, and the partitioned algorithm separating fluid from structural iterations. The formal proof is given to find the condition for convergence of this iterative procedure in the fully nonlinear setting. Several validation examples are shown to confirm the proposed convergence criteria of partitioned algorithm. The proposed strategy provides a very suitable basics for code-coupling implementation as discussed in Part II

Un couplage fluide-fluide de type Boussinesq-SPH pour des calculs à proximité de structures hydroliennes

National audienceAu cours de notre exposé, nous présenterons une stratégie de couplage pour des problèmes d'écoulements à surface libre. L'objectif est ici, par exemple, de prendre en compte la propagation de vagues sur des distances moyennes à grandes, et leur déferlement sur les côtes à une échelle plus réduite

Coupling SPH with a 1-D Boussinesq-type wave model

International audienceThe high computational cost of SPH remains problematic in dealing with wave propagation, especially when the domains considered are large. In order to overcome this difficulty, we propose to couple 2-D SPH with a 1-D Finite Difference Boussinesq-type model. The latter deals with wave propagations for most of the spatial domain, whereas SPH computations focus on the shoreline or close to off-shore structures, where a complex description of the free-surface is required. The re-use of existing codes is achieved using a generic implementation based on Component Technology. The communication between software is ensured by the middleware Component Template Library (CTL). In order to deal with open domains, open-boundaries have to be implemented for SPH, with water height and velocity varying in space and time. These velocity and water height values are then driven by the Boussinesq-type model. As an illustration of the one way coupling, we present herein two simple examples of water waves, the first one with a flat bottom, the other one representing a schematic coastal protection

Transport of isotropic particles in a partially obstructed channel flow

The transport of particles in turbulent flows is a common problem in hydraulic engineering. In this paper, the focus is set on a numerical model used to simulate the transport of small bodies (debris, algae, etc.) along a coastline. In this problem, the particles are larger than the small turbulent eddies, but smaller than the large turbulent eddies, and sufficiently diluted within the flow so that each particle does not affect the flow or the motion of other particles. A mixed Eulerian-Lagrangian approach was chosen in order to model a large flow area with sufficient information for the turbulent diffusion. This model is validated through an experiment on particles released into a partially obstructed channel flow. The measurements are then compared with simulations using two Eulerian industrial codes, Telemac-2D and OpenFoam. Finally, an application to algae bloom transport in a harbour is presented

Functional immune responses against SARS-CoV-2 variants of concern after fourth COVID-19 vaccine dose or infection in patients with blood cancer

Summary Patients with blood cancer continue to have a greater risk of inadequate immune responses following three COVID-19 vaccine doses and risk of severe COVID-19 disease. In the context of the CAPTURE study (NCT03226886) we report immune responses in 80 patients with blood cancer who received a fourth dose of BNT162b2. We measured neutralising antibody titres (NAbT) using a live virus microneutralization assay against wild-type (WT), Delta, Omicron BA.1 and BA.2 and T cell responses against WT and Omicron BA.1 using an activation-induced marker (AIM) assay. The proportion of patients with detectable NAb titres and T cell responses after the fourth vaccine dose increases compared to those after the third vaccine dose. Patients who received B cell-depleting therapies within 12 months before vaccination have the greatest risk of not having detectable NAbT. In addition, we report immune responses in 57 patients with breakthrough infections after vaccination