Antonio DeSimone

Luca Heltai, Josef Kiendl, Antonio DeSimone et al.

The interaction between thin structures and incompressible Newtonian fluids is ubiquitous both in nature and in industrial applications. In this paper we present an isogeometric formulation of such problems which exploits a boundary integral formulation of Stokes equations to model the surrounding flow, and a non linear Kirchhoff-Love shell theory to model the elastic behaviour of the structure. We propose three different coupling strategies: a monolithic, fully implicit coupling, a staggered, elasticity driven coupling, and a novel semi-implicit coupling, where the effect of the surrounding flow is incorporated in the non-linear terms of the solid solver through its damping characteristics. The novel semi-implicit approach is then used to demonstrate the power and robustness of our method, which fits ideally in the isogeometric paradigm, by exploiting only the boundary representation (B-Rep) of the thin structure middle surface.

François Alouges, Antonio DeSimone, Luca Heltai et al.

We study self propelled stokesian robots composed of assemblies of balls, in dimensions 2 and 3, and prove that they are able to control their position and orientation. This is a result of controllability, and its proof relies on applying Chow's theorem in an analytic framework, similarly to what has been done in [3] for an axisymmetric system swimming along the axis of symmetry. However, we simplify drastically the analyticity result given in [3] and apply it to a situation where more complex swimmers move either in a plane or in three-dimensional space, hence experiencing also rotations. We then focus our attention on energetically optimal strokes, which we are able to compute numerically. Some examples of computed optimal strokes are discussed in detail.

Andrea Mola, Luca Heltai, Antonio DeSimone

We present an innovative numerical discretization of the equations of inviscid potential flow for the simulation of three dimensional unsteady and nonlinear water waves generated by a ship hull advancing in water. The equations of motion are written in a semi-Lagrangian framework, and the resulting integro-differential equations are discretized in space via an adaptive iso-parametric collocation Boundary Element Method, and in time via adaptive implicit Backward Differentiation Formulas (BDF) with variable step and variable order. When the velocity of the advancing ship hull is non-negligible, the semi-Lagrangian formulation (also known as Arbitrary Lagrangian Eulerian formulation, or ALE) of the free surface equations contains dominant transport terms which are stabilized with a Streamwise Upwind Petrov-Galerkin (SUPG) method. The SUPG stabilization allows automatic and robust adaptation of the spatial discretization with unstructured quadrilateral grids. Preliminary results are presented where we compare our numerical model with experimental results on the case of a Wigley hull advancing in calm water with fixed sink and trim.

Marino Arroyo, Antonio DeSimone, Luca Heltai

Fluid membranes made out of lipid bilayers are the fundamental separation structure in eukaryotic cells. Many physiological processes rely on dramatic shape and topological changes (e.g. fusion, fission) of fluid membrane systems. Fluidity is key to the versatility and constant reorganization of lipid bilayers. Here, we study the role of the membrane intrinsic viscosity, arising from the friction of the lipid molecules as they rearrange to accommodate shape changes, in the dynamics of morphological changes of fluid vesicles. In particular, we analyze the competition between the membrane viscosity and the viscosity of the bulk fluid surrounding the vesicle as the dominant dissipative mechanism. We consider the relaxation dynamics of fluid vesicles put in an out-of-equilibrium state, but conclusions can be drawn regarding the kinetics or power consumption in regulated shape changes in the cell. On the basis of numerical calculations, we find that the dynamics arising from the membrane viscosity are qualitatively different from the dynamics arising from the bulk viscosity. When these two dissipation mechanisms are put in competition, we find that for small vesicles the membrane dissipation dominates, with a relaxation time that scales as the size of the vesicle to the power 2. For large vesicles, the bulk dissipation dominates, and the exponent in the relaxation time vs. size relation is 3.

François Alouges, Antonio DeSimone, Luca Heltai

We propose a computational method to solve optimal swimming problems, based on the boundary integral formulation of the hydrodynamic interaction between swimmer and surrounding fluid and direct constrained minimization of the energy consumed by the swimmer. We apply our method to axisymmetric model examples. We consider a classical model swimmer (the three-sphere swimmer of Golestanian et al.) as well as a novel axisymmetric swimmer inspired by the observation of biological micro-organisms.