TY - JOUR T1 - Model hierarchies and higher-order discretisation of time-dependent thin-film free boundary problems with dynamic contact angle JF - Journal of Computational Physics Y1 - 2022 A1 - Dirk Peschka A1 - Luca Heltai VL - 464 ER - TY - JOUR T1 - The \textttdeal.II Library, Version 9.4 JF - Journal of Numerical Mathematics Y1 - 2022 A1 - Daniel Arndt A1 - Wolfgang Bangerth Marco Feder A1 - Marc Fehling A1 - Rene Gassmöller A1 - Timo Heister A1 - Luca Heltai A1 - Martin Kronbichler A1 - Matthias Maier A1 - Peter Munch A1 - Jean-Paul Pelteret A1 - Simon Sticko A1 - Bruno Turcksin A1 - David Wells ER - TY - JOUR T1 - Variational Approach to Fluid–Structure Interaction via GENERIC JF - Journal of Non-Equilibrium Thermodynamics Y1 - 2022 A1 - Dirk Peschka A1 - Andrea Zafferi A1 - Luca Heltai A1 - Marita Thomas ER - TY - JOUR T1 - Vibration Analysis of Piezoelectric Kirchhoff-Love shells based on Catmull-Clark Subdivision Surfaces JF - International Journal for Numerical Methods in Engineering Y1 - 2022 A1 - Zhaowei Liu A1 - Andrew McBride A1 - Prashant Saxena A1 - Luca Heltai A1 - Yilin Qu A1 - Paul Steinmann ER - TY - JOUR T1 - The deal.II Library, Version 9.3 JF - Journal of Numerical Mathematics Y1 - 2021 A1 - Daniel Arndt A1 - Wolfgang Bangerth A1 - Bruno Blais A1 - Marc Fehling A1 - Rene Gassmöller A1 - Timo Heister A1 - Luca Heltai A1 - Uwe Köcher A1 - Martin Kronbichler A1 - Matthias Maier A1 - Peter Munch A1 - Jean-Paul Pelteret A1 - Sebastian Proell A1 - Konrad Simon A1 - Bruno Turcksin A1 - David Wells A1 - Jiaqi Zhang UR - https://doi.org/10.1515/jnma-2021-0081 ER - TY - JOUR T1 - Multiscale coupling of one-dimensional vascular models and elastic tissues JF - Annals of Biomedical Engineering Y1 - 2021 A1 - Luca Heltai A1 - Alfonso Caiazzo A1 - Lucas Müeller ER - TY - JOUR T1 - Propagating geometry information to finite element computations JF - Transactions on Mathematical Software Y1 - 2021 A1 - Luca Heltai A1 - Wolfgang Bangerth A1 - Martin Kronbichler A1 - Andrea Mola VL - 47 IS - 4 ER - TY - JOUR T1 - Quasi-optimal mesh sequence construction through Smoothed Adaptive Finite Element Methods JF - SIAM Journal on Scientific Computing Y1 - 2021 A1 - Ornela Mulita A1 - Stefano Giani A1 - Luca Heltai ER - TY - JOUR T1 - Smoothed-adaptive perturbed inverse iteration for elliptic eigenvalue problems JF - Computational Methods in Applied Mathematics Y1 - 2021 A1 - Stefano Giani A1 - Luka Grubisic A1 - Luca Heltai A1 - Ornela Mulita VL - 21 ER - TY - JOUR T1 - The deal.II finite element library: Design, features, and insights JF - Computers and Mathematics with Applications Y1 - 2020 A1 - Daniel Arndt A1 - Wolfgang Bangerth A1 - Denis Davydov A1 - Timo Heister A1 - Luca Heltai A1 - Martin Kronbichler A1 - Matthias Maier A1 - Jean-Paul Pelteret A1 - Bruno Turcksin A1 - David Wells UR - https://doi.org/10.1016/j.camwa.2020.02.022 ER - TY - JOUR T1 - The deal.II library, Version 9.2 JF - Journal of Numerical Mathematics Y1 - 2020 A1 - Daniel Arndt A1 - Wolfgang Bangerth A1 - Bruno Blais A1 - Thomas C. Clevenger A1 - Marc Fehling A1 - Alexander V. Grayver A1 - Timo Heister A1 - Luca Heltai A1 - Martin Kronbichler A1 - Matthias Maier A1 - Peter Munch A1 - Jean-Paul Pelteret A1 - Reza Rastak A1 - Ignacio Tomas A1 - Bruno Turcksin A1 - Zhuoran Wang A1 - David Wells VL - 28 ER - TY - JOUR T1 - Multiscale modeling of fiber reinforced materials via non-matching immersed methods JF - Computers & Structures Y1 - 2020 A1 - Giovanni Alzetta A1 - Luca Heltai VL - 239 N1 - To appear ER - TY - JOUR T1 - A numerical study of the jerky crack growth in elastoplastic materials with localized plasticity JF - Journal of Convex Analysis Y1 - 2020 A1 - Gianni Dal Maso A1 - Luca Heltai UR - https://arxiv.org/abs/2004.12705 N1 - To appear ER - TY - JOUR T1 - A priori error estimates of regularized elliptic problems JF - Numerische Mathematik Y1 - 2020 A1 - Luca Heltai A1 - Wenyu Lei AB - Approximations of the Dirac delta distribution are commonly used to create sequences of smooth functions approximating nonsmooth (generalized) functions, via convolution. In this work we show a-priori rates of convergence of this approximation process in standard Sobolev norms, with minimal regularity assumptions on the approximation of the Dirac delta distribution. The application of these estimates to the numerical solution of elliptic problems with singularly supported forcing terms allows us to provide sharp \$\$H\^1\$\$and \$\$L\^2\$\$error estimates for the corresponding regularized problem. As an application, we show how finite element approximations of a regularized immersed interface method results in the same rates of convergence of its non-regularized counterpart, provided that the support of the Dirac delta approximation is set to a multiple of the mesh size, at a fraction of the implementation complexity. Numerical experiments are provided to support our theories. VL - 146 SN - 0945-3245 UR - https://doi.org/10.1007/s00211-020-01152-w ER - TY - JOUR T1 - A priori error estimates of regularized elliptic problems JF - Numerische Mathematik Y1 - 2020 A1 - Luca Heltai A1 - Wenyu Lei ER - TY - JOUR T1 - The deal.II Library, Version 9.1 JF - Journal of Numerical Mathematics Y1 - 2019 A1 - Daniel Arndt A1 - Wolfgang Bangerth A1 - Thomas C. Clevenger A1 - Denis Davydov A1 - Marc Fehling A1 - Garcia-Sanchez, Daniel A1 - Harper, Graham A1 - Timo Heister A1 - Luca Heltai A1 - Martin Kronbichler A1 - Maguire Kynch, Ross A1 - Matthias Maier A1 - Jean-Paul Pelteret A1 - Bruno Turcksin A1 - David Wells ER - TY - JOUR T1 - The deal.II Library, Version 9.1 JF - Journal of Numerical Mathematics Y1 - 2019 A1 - Daniel Arndt A1 - Wolfgang Bangerth A1 - Thomas C. Clevenger A1 - Denis Davydov A1 - Marc Fehling A1 - Garcia-Sanchez, Daniel A1 - Harper, Graham A1 - Timo Heister A1 - Luca Heltai A1 - Martin Kronbichler A1 - Maguire Kynch, Ross A1 - Matthias Maier A1 - Jean-Paul Pelteret A1 - Bruno Turcksin A1 - David Wells AB - This paper provides an overview of the new features of the finite element library deal.II, version 9.1. ER - TY - JOUR T1 - Error estimates in weighted Sobolev norms for finite element immersed interface methods JF - Computers & Mathematics with Applications Y1 - 2019 A1 - Luca Heltai A1 - Nella Rotundo PB - Elsevier BV VL - 78 UR - https://doi.org/10.1016/j.camwa.2019.05.029 ER - TY - JOUR T1 - Multiscale modeling of vascularized tissues via non-matching immersed methods JF - International Journal for Numerical Methods in Biomedical Engineering Y1 - 2019 A1 - Luca Heltai A1 - Alfonso Caiazzo VL - 35 UR - https://doi.org/10.1002/cnm.3264 ER - TY - JOUR T1 - Accelerating the iterative solution of convection-diffusion problems using singular value decomposition JF - NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Y1 - 2018 A1 - Giuseppe Pitton A1 - Luca Heltai UR - https://arxiv.org/abs/1807.09467 ER - TY - JOUR T1 - An authenticated theoretical modeling of electrified fluid jet in core–shell nanofibers production JF - JOURNAL OF INDUSTRIAL TEXTILES Y1 - 2018 A1 - Rafiei, S. A1 - Noroozi, B. A1 - Luca Heltai A1 - Haghi, A. K. VL - 47 ER - TY - JOUR T1 - deal2lkit: A toolkit library for high performance programming in deal.II JF - SOFTWAREX Y1 - 2018 A1 - Alberto Sartori A1 - Nicola Giuliani A1 - Mauro Bardelloni A1 - Luca Heltai VL - 7 ER - TY - JOUR T1 - The deal.II Library, Version 9.0 JF - JOURNAL OF NUMERICAL MATHEMATICS Y1 - 2018 A1 - Giovanni Alzetta A1 - Daniel Arndt A1 - Wolfgang Bangerth A1 - Boddu, Vishal A1 - Brands, Benjamin A1 - Denis Davydov A1 - Rene Gassmöller A1 - Timo Heister A1 - Luca Heltai A1 - Kormann, Katharina A1 - Martin Kronbichler A1 - Matthias Maier A1 - Jean-Paul Pelteret A1 - Bruno Turcksin A1 - David Wells UR - https://doi.org/10.1515/jnma-2018-0054 ER - TY - CHAP T1 - A distributed lagrange formulation of the finite element immersed boundary method for fluids interacting with compressible solids T2 - Mathematical and Numerical Modeling of the Cardiovascular System and Applications Y1 - 2018 A1 - Boffi, Daniele A1 - Gastaldi, Lucia A1 - Luca Heltai JF - Mathematical and Numerical Modeling of the Cardiovascular System and Applications PB - Springer International Publishing CY - Cham VL - 16 UR - https://arxiv.org/abs/1712.02545v1 ER - TY - JOUR T1 - Iterative map-making with two-level preconditioning for polarized cosmic microwave background data sets. A worked example for ground-based experiments JF - ASTRONOMY & ASTROPHYSICS Y1 - 2018 A1 - Puglisi, Giuseppe A1 - Poletti, Davide A1 - Fabbian, Giulio A1 - Baccigalupi, Carlo A1 - Luca Heltai A1 - Stompor, Radek VL - 618 UR - https://arxiv.org/abs/1801.08937 ER - TY - JOUR T1 - NURBS-SEM: A hybrid spectral element method on NURBS maps for the solution of elliptic PDEs on surfaces JF - COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING Y1 - 2018 A1 - Giuseppe Pitton A1 - Luca Heltai VL - 338 UR - https://arxiv.org/abs/1804.08271 ER - TY - JOUR T1 - Predicting and Optimizing Microswimmer Performance from the Hydrodynamics of Its Components: The Relevance of Interactions JF - SOFT ROBOTICS Y1 - 2018 A1 - Nicola Giuliani A1 - Luca Heltai A1 - Antonio DeSimone VL - 5 UR - https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6094362/ ER - TY - JOUR T1 - π-BEM : A flexible parallel implementation for adaptive , geometry aware , and high order boundary element methods JF - Advances in Engineering Software Y1 - 2018 A1 - Nicola Giuliani A1 - Andrea Mola A1 - Luca Heltai VL - 121 ER - TY - JOUR T1 - The deal.II Library, Version 8.5 JF - JOURNAL OF NUMERICAL MATHEMATICS Y1 - 2017 A1 - Daniel Arndt A1 - Wolfgang Bangerth A1 - Denis Davydov A1 - Timo Heister A1 - Luca Heltai A1 - Martin Kronbichler A1 - Matthias Maier A1 - Jean-Paul Pelteret A1 - Bruno Turcksin A1 - David Wells VL - 25 UR - https://www.dealii.org/deal85-preprint.pdf ER - TY - JOUR T1 - A natural framework for isogeometric fluid-structure interaction based on BEM-shell coupling JF - COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING Y1 - 2017 A1 - Luca Heltai A1 - Kiendl, J. A1 - Antonio DeSimone A1 - Alessandro Reali VL - 316 UR - http://cdsads.u-strasbg.fr/abs/2017CMAME.316..522H ER - TY - JOUR T1 - Wet and Dry Transom Stern Treatment for Unsteady and Nonlinear Potential Flow Model for Naval Hydrodynamics Simulations JF - Journal of Ship Research Y1 - 2017 A1 - Andrea Mola A1 - Luca Heltai A1 - Antonio DeSimone AB -

We present a model for the fast evaluation of the total drag of ship hulls operating in both wet and dry transom stern conditions, in calm or wavy water, based on the combination of an unsteady semi-Lagrangian potential flow formulation with fully nonlinear free-surface treatment, experimental correlations, and simplified viscous drag modeling. The implementation is entirely based on open source libraries. The spatial discretization is solved using a streamline upwind Petrov‐Galerkin stabilization of an iso-parametric, collocation based, boundary element method, implemented using the open source library deal.II. The resulting nonlinear differential-algebraic system is integrated in time using implicit backward differentiation formulas, implemented in the open source library SUNDIALS. The Open CASCADE library is used to interface the model directly with computer-aided design data structures. The model accounts automatically for hulls with a transom stern, both in wet and dry regimes, by using a specific treatment of the free-surface nodes on the stern edge that automatically detects when the hull advances at low speeds. In this case, the transom stern is partially immersed, and a pressure patch is applied on the water surface detaching from the transom stern, to recover the gravity effect of the recirculating water on the underlying irrotational flow domain. The parameters of the model used to impose the pressure patch are approximated from experimental relations found in the literature. The test cases considered are those of the U.S. Navy Combatant DTMB-5415 and the National Physical Laboratory hull. Comparisons with experimental data on quasi-steady test cases for both water elevation and total hull drag are presented and discussed. The quality of the results obtained on quasi-steady simulations suggests that this model can represent a promising alternative to current unsteady solvers for simulations with Froude numbers below 0.35.

VL - 61 ER - TY - JOUR T1 - The deal.II Library, Version 8.3 JF - ARCHIVE OF NUMERICAL SOFTWARE Y1 - 2016 A1 - Wolfgang Bangerth A1 - Timo Heister A1 - Luca Heltai A1 - G. Kanschat A1 - Martin Kronbichler A1 - Matthias Maier A1 - Bruno Turcksin VL - 4 UR - http://nbn-resolving.de/urn:nbn:de:bsz:16-ans-231226 ER - TY - JOUR T1 - The deal.II library, Version 8.4 JF - JOURNAL OF NUMERICAL MATHEMATICS Y1 - 2016 A1 - Wolfgang Bangerth A1 - Denis Davydov A1 - Timo Heister A1 - Luca Heltai A1 - G. Kanschat A1 - Martin Kronbichler A1 - Matthias Maier A1 - Bruno Turcksin A1 - David Wells VL - 24 UR - https://www.math.clemson.edu/ heister/preprints/deal84-preprint.pdf ER - TY - JOUR T1 - Error Estimates of B-spline based finite-element method for the wind-driven ocean circulation JF - JOURNAL OF SCIENTIFIC COMPUTING Y1 - 2016 A1 - Rotundo, N. A1 - Kim, T. -Y. A1 - Jiang, W. A1 - Luca Heltai A1 - Fried, E. VL - 69 ER - TY - RPRT T1 - Isogeometric analysis-based reduced order modelling for incompressible linear viscous flows in parametrized shapes Y1 - 2016 A1 - Filippo Salmoiraghi A1 - F. Ballarin A1 - Luca Heltai A1 - Gianluigi Rozza AB - In this work we provide a combination of isogeometric analysis with reduced order modelling techniques, based on proper orthogonal decomposition, to guarantee computational reduction for the numerical model, and with free-form deformation, for versatile geometrical parametrization. We apply it to computational fluid dynamics problems considering a Stokes flow model. The proposed reduced order model combines efficient shape deformation and accurate and stable velocity and pressure approximation for incompressible viscous flows, computed with a reduced order method. Efficient offine-online computational decomposition is guaranteed in view of repetitive calculations for parametric design and optimization problems. Numerical test cases show the efficiency and accuracy of the proposed reduced order model. PB - Springer, AMOS Advanced Modelling and Simulation in Engineering Sciences UR - http://urania.sissa.it/xmlui/handle/1963/35199 U1 - 35493 U2 - Mathematics U4 - 1 U5 - MAT/08 ER - TY - JOUR T1 - LinearOperator – a generic, high-level expression syntax for linear algebra JF - COMPUTERS & MATHEMATICS WITH APPLICATIONS Y1 - 2016 A1 - Matthias Maier A1 - Mauro Bardelloni A1 - Luca Heltai VL - 72 ER - TY - CONF T1 - Ship Sinkage and Trim Predictions Based on a CAD Interfaced Fully Nonlinear Potential Model T2 - The 26th International Ocean and Polar Engineering Conference Y1 - 2016 A1 - Andrea Mola A1 - Luca Heltai A1 - Antonio DeSimone A1 - Massimiliano Berti JF - The 26th International Ocean and Polar Engineering Conference PB - International Society of Offshore and Polar Engineers VL - 3 ER - TY - JOUR T1 - Benchmarking the Immersed Finite Element Method for Fluid-Structure Interaction Problems JF - Computers and Mathematics with Applications 69 (2015) 1167–1188 Y1 - 2015 A1 - Roy Saswati A1 - Luca Heltai A1 - Francesco Costanzo AB - We present an implementation of a fully variational formulation of an immersed methods for fluid-structure interaction problems based on the finite element method. While typical implementation of immersed methods are characterized by the use of approximate Dirac delta distributions, fully variational formulations of the method do not require the use of said distributions. In our implementation the immersed solid is general in the sense that it is not required to have the same mass density and the same viscous response as the surrounding fluid. We assume that the immersed solid can be either viscoelastic of differential type or hyperelastic. Here we focus on the validation of the method via various benchmarks for fluid-structure interaction numerical schemes. This is the first time that the interaction of purely elastic compressible solids and an incompressible fluid is approached via an immersed method allowing a direct comparison with established benchmarks. U1 - 34633 U2 - Mathematics U4 - 1 U5 - MAT/08 ER - TY - JOUR T1 - Deal2lkit: a Toolkit Library for High Performance Programming in deal.II Y1 - 2015 A1 - Alberto Sartori A1 - Nicola Giuliani A1 - Mauro Bardelloni A1 - Luca Heltai AB - We present version 1.0.0 of the deal2lkit (deal.II ToolKit) library. deal2lkit is a collection of modules and classes for the general purpose finite element library deal.II. Its principal aim is to provide a high level interface, controlled via parameter files, for those steps that are common in all finite element programs: mesh generation, selection of the finite element type, application of boundary conditions and many others. Each module can be used as a building block independently on the others, and can be integrated in existing finite element codes based on deal.II, drastically reducing the size of programs, rendering their use automatically parametrised, and reducing the overall time-to-market of finite element programming. Moreover, deal2lkit features interfaces with the SUNDIALS (SUite of Nonlinear and DIfferential/ALgebraic equation Solvers) and ASSIMP (Open Asset Import Library) libraries. Some examples are provided which show the aim and scopes of deal2lkit. The deal2lkit library is released under the GNU Lesser General Public License (LGPL) and can be retrieved from the deal2lkit repository https://github.com/mathLab/deal2lkit. PB - SISSA UR - http://urania.sissa.it/xmlui/handle/1963/35006 U1 - 35235 U2 - Mathematics U4 - 1 U5 - MAT/08 ER - TY - JOUR T1 - The deal.II Library, Version 8.2 JF - Archive of Numerical Software, vol. 3, n. 100, (2015), pages : 1-8 Y1 - 2015 A1 - Wolfgang Bangerth A1 - Timo Heister A1 - Luca Heltai A1 - G. Kanschat A1 - Martin Kronbichler A1 - Matthias Maier A1 - Bruno Turcksin A1 - T. D. Young AB - This paper provides an overview of the new features of the finite element library deal.II version 8.2 UR - http://urania.sissa.it/xmlui/handle/1963/34464 U1 - 34637 U2 - Mathematics U4 - 1 U5 - MAT/08 ER - TY - JOUR T1 - FEM SUPG stabilisation of mixed isoparametric BEMs: application to linearised free surface flows JF - Engineering Analysis with Boundary Elements 59 (2015), pp. 8-22 Y1 - 2015 A1 - Nicola Giuliani A1 - Andrea Mola A1 - Luca Heltai A1 - L. Formaggia AB -

In finite element formulations, transport dominated problems are often stabilised through the Streamline-Upwind-Petrov–Galerkin (SUPG) method. Its application is straightforward when the problem at hand is solved using Galerkin methods. Applications of boundary integral formulations often resort to collocation techniques which are computationally more tractable. In this framework, the Galerkin method and the stabilisation may still be used to successfully apply boundary conditions and resolve instabilities that are frequently observed in transport dominated problems. We apply this technique to an adaptive collocation boundary element method for the solution of stationary potential flows, where we solve a mixed Poisson problem in boundary integral form, with the addition of linearised free surface boundary conditions. We use a mixed boundary element formulation to allow for different finite dimensional spaces describing the flow potential and its normal derivative, and we validate our method simulating the flow around both a submerged body and a surface piercing body. The coupling of mixed surface finite elements and strongly consistent stabilisation techniques with boundary elements opens up the possibility to use non conformal unstructured grids with local refinement, without introducing the inconsistencies of other stabilisation techniques based on up-winding and finite difference schemes.

UR - http://urania.sissa.it/xmlui/handle/1963/34466 U1 - 34640 U2 - Mathematics U4 - 1 U5 - MAT/08 ER - TY - JOUR T1 - Reduced Basis Isogeometric Methods (RB-IGA) for the real-time simulation of potential flows about parametrized NACA airfoils JF - Comput Methods Appl Mech Eng. 2015;284:1147–1180 Y1 - 2015 A1 - Andrea Manzoni A1 - Filippo Salmoiraghi A1 - Luca Heltai AB - We present a Reduced Basis (RB) method based on Isogeometric Analysis (IGA) for the rapid and reliable evaluation of PDE systems characterized by complex geometrical features. At the current state of the art, this is the first case of coupling between RB and IGA methods. The construction of the RB method relies on an Isogeometric Boundary Element Method (IGA-BEM) as the high-fidelity technique, allowing a direct interface with Computer Aided Design (CAD) tools. A suitable Empirical Interpolation Method (EIM) ensures an efficient offline/online decomposition between the construction and the evaluation of the RB method. We consider the real-time simulation of potential flows past airfoils, parametrized with respect to the angle of attack and the NACA number identifying their shape, and we provide a validation of our methodology with respect to experimental data and reference numerical codes, showing in both cases a very good agreement.We present a Reduced Basis (RB) method based on Isogeometric Analysis (IGA) for the rapid and reliable evaluation of PDE systems characterized by complex geometrical features. At the current state of the art, this is the first case of coupling between RB and IGA methods. The construction of the RB method relies on an Isogeometric Boundary Element Method (IGA-BEM) as the high-fidelity technique, allowing a direct interface with Computer Aided Design (CAD) tools. A suitable Empirical Interpolation Method (EIM) ensures an efficient offline/online decomposition between the construction and the evaluation of the RB method. We consider the real-time simulation of potential flows past airfoils, parametrized with respect to the angle of attack and the NACA number identifying their shape, and we provide a validation of our methodology with respect to experimental data and reference numerical codes, showing in both cases a very good agreement. U1 - 34587 U2 - Mathematics U4 - 1 U5 - MAT/08 ER - TY - Generic T1 - A fully nonlinear potential model for ship hydrodynamics directly interfaced with CAD data structures T2 - Proceedings of the 24th International Ocean and Polar Engineering Conference, Busan, 2014 Y1 - 2014 A1 - Andrea Mola A1 - Luca Heltai A1 - Antonio DeSimone KW - ship hydrodynamics AB - We present a model for ship hydrodynamics simulations currently under development at SISSA. The model employs potential flow theory and fully nonlinear free surface boundary conditions. The spatial discretization of the equations is performed by means of a collocation BEM. This gives rise to a Differential Algbraic Equations (DAE) system, solved using an implicit BDF scheme to time advance the solution. The model has been implemented into a C++ software able to automatically generate the computational grids from the CAD geometry of the hull. Numerical results on Kriso KCS and KVLCC2 hulls are presented and discussed. JF - Proceedings of the 24th International Ocean and Polar Engineering Conference, Busan, 2014 PB - SISSA U1 - 7357 U2 - Mathematics U4 - 1 U5 - MAT/08 ANALISI NUMERICA ER - TY - JOUR T1 - Nonsingular Isogeometric Boundary Element Method for Stokes Flows in 3D Y1 - 2014 A1 - Luca Heltai A1 - Marino Arroyo A1 - Antonio DeSimone KW - Isogeometric Analysis AB - Isogeometric analysis (IGA) is emerging as a technology bridging Computer Aided Geometric Design (CAGD), most commonly based on Non-Uniform Rational B-Splines (NURBS) surfaces, and engineering analysis. In finite element and boundary element isogeometric methods (FE-IGA and IGA-BEM), the NURBS basis functions that de- scribe the geometry define also the approximation spaces. In the FE-IGA approach, the surfaces generated by the CAGD tools need to be extended to volumetric descriptions, a major open problem in 3D. This additional passage can be avoided in principle when the partial differential equations to be solved admit a formulation in terms of bound- ary integral equations, leading to Boundary Element Isogeometric Analysis (IGA-BEM). The main advantages of such an approach are given by the dimensionality reduction of the problem (from volumetric-based to surface-based), by the fact that the interface with CAGD tools is direct, and by the possibility to treat exterior problems, where the computational domain is infinite. By contrast, these methods produce system matrices which are full, and require the integration of singular kernels. In this paper we address the second point and propose a nonsingular formulation of IGA-BEM for 3D Stokes flows, whose convergence is carefully tested numerically. Standard Gaussian quadrature rules suffice to integrate the boundary integral equations, and carefully chosen known exact solutions of the interior Stokes problem are used to correct the resulting matrices, extending the work by Klaseboer et al. [27] to IGA-BEM. PB - Elsevier UR - http://hdl.handle.net/1963/6326 U1 - 6250 U2 - Mathematics U4 - 1 U5 - MAT/08 ANALISI NUMERICA ER - TY - CONF T1 - Potential Model for Ship Hydrodynamics Simulations Directly Interfaced with CAD Data Structures T2 - The 24th International Ocean and Polar Engineering Conference Y1 - 2014 A1 - Andrea Mola A1 - Luca Heltai A1 - Antonio DeSimone A1 - Massimiliano Berti JF - The 24th International Ocean and Polar Engineering Conference PB - International Society of Offshore and Polar Engineers VL - 4 ER - TY - RPRT T1 - The deal.II Library, Version 8.1 Y1 - 2013 A1 - Wolfgang Bangerth A1 - Timo Heister A1 - Luca Heltai A1 - G. Kanschat A1 - Martin Kronbichler A1 - Matthias Maier A1 - Bruno Turcksin A1 - T. D. Young AB - This paper provides an overview of the new features of the finite element library deal.II version 8.0. PB - SISSA UR - http://hdl.handle.net/1963/7236 N1 - 5 pages U1 - 7272 U2 - Mathematics U4 - 1 U5 - MAT/08 ANALISI NUMERICA ER - TY - JOUR T1 - A stable and adaptive semi-Lagrangian potential model for unsteady and nonlinear ship-wave interactions JF - Engineering Analysis with Boundary Elements, 37(1):128 – 143, 2013. Y1 - 2013 A1 - Andrea Mola A1 - Luca Heltai A1 - Antonio DeSimone KW - Unsteady ship-wave interaction AB -

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-diff erential equations are discretized in space via an adaptive iso-parametric collocation Boundary Element Method, and in time via adaptive implicit Backward Di erentiation 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 fi xed sink and trim.

PB - SISSA UR - http://hdl.handle.net/1963/5669 U1 - 5457 U2 - Mathematics U3 - Functional Analysis and Applications U4 - -1 ER - TY - CHAP T1 - Computing optimal strokes for low reynolds number swimmers T2 - Natural locomotion in fluids and on surfaces : swimming, flying, and sliding / editors Stephen Childress, Anette Hosoi, William W. Schultz, and Z. Jane Wang, editors, Y1 - 2012 A1 - Antonio DeSimone A1 - Luca Heltai A1 - François Alouges A1 - Lefebvre-Lepot Aline KW - Numerical analysis. AB -

We discuss connections between low-Reynolds-number swimming and geometric control theory, and present a general algorithm for the numerical computation of energetically optimal strokes. As an illustration of our approach, we show computed motility maps and optimal strokes for two model swimmers.

JF - Natural locomotion in fluids and on surfaces : swimming, flying, and sliding / editors Stephen Childress, Anette Hosoi, William W. Schultz, and Z. Jane Wang, editors, PB - Springer SN - 9781461439967 UR - http://hdl.handle.net/1963/6445 U1 - 6381 U2 - Mathematics U4 - 1 U5 - MAT/08 ANALISI NUMERICA ER - TY - RPRT T1 - A Fully Coupled Immersed Finite Element Method for Fluid Structure Interaction via the Deal.II Library Y1 - 2012 A1 - Luca Heltai A1 - Saswati Roy A1 - Francesco Costanzo KW - Finite Element Method KW - Immersed Boundary Method KW - Immersed Finite Element Method AB - We present the implementation of a solution scheme for fluid-structure\\r\\ninteraction problems via the finite element software library deal.II. The\\r\\nsolution scheme is an immersed finite element method in which two independent discretizations are used for the fluid and immersed deformable body. In this type of formulation the support of the equations of motion of the fluid is extended to cover the union of the solid and fluid domains. The equations of motion over the extended solution domain govern the flow of a fluid under the action of a body force field. This body force field informs the fluid of the presence of the immersed solid. The velocity field of the immersed solid is the restriction over the immersed domain of the velocity field in the extended equations of motion. The focus of this paper is to show how the determination of the motion of the immersed domain is carried out in practice. We show that our implementation is general, that is, it is not dependent on a specific choice of the finite element spaces over the immersed solid and the extended fluid domains. We present some preliminary results concerning the accuracy of the proposed method. PB - SISSA UR - http://hdl.handle.net/1963/6255 N1 - 28 pages, 9 figures U1 - 6172 U2 - Mathematics U3 - Functional Analysis and Applications U4 - 1 U5 - MAT/08 ANALISI NUMERICA ER - TY - JOUR T1 - Reverse engineering the euglenoid movement JF - Proceedings of the National Academy of Sciences of the United States of America. Volume 109, Issue 44, 30 October 2012, Pages 17874-17879 Y1 - 2012 A1 - Marino Arroyo A1 - Luca Heltai A1 - Daniel Millán A1 - Antonio DeSimone KW - microswimmers AB - Euglenids exhibit an unconventional motility strategy amongst unicellular eukaryotes, consisting of large-amplitude highly concerted deformations of the entire body (euglenoid movement or metaboly). A plastic cell envelope called pellicle mediates these deformations. Unlike ciliary or flagellar motility, the biophysics of this mode is not well understood, including its efficiency and molecular machinery. We quantitatively examine video recordings of four euglenids executing such motions with statistical learning methods. This analysis reveals strokes of high uniformity in shape and pace. We then interpret the observations in the light of a theory for the pellicle kinematics, providing a precise understanding of the link between local actuation by pellicle shear and shape control. We systematically understand common observations, such as the helical conformations of the pellicle, and identify previously unnoticed features of metaboly. While two of our euglenids execute their stroke at constant body volume, the other two exhibit deviations of about 20% from their average volume, challenging current models of low Reynolds number locomotion. We find that the active pellicle shear deformations causing shape changes can reach 340%, and estimate the velocity of the molecular motors. Moreover, we find that metaboly accomplishes locomotion at hydrodynamic efficiencies comparable to those of ciliates and flagellates. Our results suggest new quantitative experiments, provide insight into the evolutionary history of euglenids, and suggest that the pellicle may serve as a model for engineered active surfaces with applications in microfluidics. UR - http://hdl.handle.net/1963/6444 U1 - 6380 U2 - Mathematics U4 - 1 U5 - MAT/08 ANALISI NUMERICA ER - TY - CONF T1 - A stable semi-lagrangian potential method for the simulation of ship interaction with unsteady and nonlinear waves T2 - 17th Int. Conf. Ships Shipp. Res. Y1 - 2012 A1 - Andrea Mola A1 - Luca Heltai A1 - Antonio DeSimone JF - 17th Int. Conf. Ships Shipp. Res. ER - TY - JOUR T1 - Variational implementation of immersed finite element methods JF - Computer Methods in Applied Mechanics and Engineering. Volume 229-232, 1 July 2012, Pages 110-127 Y1 - 2012 A1 - Luca Heltai A1 - Francesco Costanzo KW - Turbulent flow AB -

Dirac-delta distributions are often crucial components of the solid-fluid coupling operators in immersed solution methods for fluid-structure interaction (FSI) problems. This is certainly so for methods like the Immersed Boundary Method (IBM) or the Immersed Finite Element Method (IFEM), where Dirac-delta distributions are approximated via smooth functions. By contrast, a truly variational formulation of immersed methods does not require the use of Dirac-delta distributions, either formally or practically. This has been shown in the Finite Element Immersed Boundary Method (FEIBM), where the variational structure of the problem is exploited to avoid Dirac-delta distributions at both the continuous and the discrete level. In this paper, we generalize the FEIBM to the case where an incompressible Newtonian fluid interacts with a general hyperelastic solid. Specifically, we allow (i) the mass density to be different in the solid and the fluid, (ii) the solid to be either viscoelastic of differential type or purely elastic, and (iii) the solid to be and either compressible or incompressible. At the continuous level, our variational formulation combines the natural stability estimates of the fluid and elasticity problems. In immersed methods, such stability estimates do not transfer to the discrete level automatically due to the non- matching nature of the finite dimensional spaces involved in the discretization. After presenting our general mathematical framework for the solution of FSI problems, we focus in detail on the construction of natural interpolation operators between the fluid and the solid discrete spaces, which guarantee semi-discrete stability estimates and strong consistency of our spatial discretization.

PB - Elsevier UR - http://hdl.handle.net/1963/6462 N1 - 42 pages, 5 figures, Revision 1 U1 - 6389 U2 - Mathematics U4 - 1 U5 - MAT/05 ANALISI MATEMATICA ER - TY - JOUR T1 - Numerical Strategies for Stroke Optimization of Axisymmetric Microswimmers JF - Mathematical Models and Methods in Applied Sciences 21 (2011) 361-387 Y1 - 2011 A1 - François Alouges A1 - Antonio DeSimone A1 - Luca Heltai KW - Optimal swimming AB - 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. PB - World Scientific UR - http://hdl.handle.net/1963/3657 U1 - 648 U2 - Mathematics U3 - Functional Analysis and Applications ER - TY - RPRT T1 - Optimally swimming Stokesian Robots Y1 - 2010 A1 - François Alouges A1 - Antonio DeSimone A1 - Luca Heltai A1 - Aline Lefebvre A1 - Benoit Merlet AB - We study self propelled stokesian robots composed of assemblies of balls, in dimen-\\nsions 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\\nthe 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. UR - http://hdl.handle.net/1963/3929 U1 - 472 U2 - Mathematics U3 - Functional Analysis and Applications ER - TY - RPRT T1 - The role of membrane viscosity in the dynamics of fluid membranes Y1 - 2010 A1 - Marino Arroyo A1 - Antonio DeSimone A1 - Luca Heltai AB - 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. UR - http://hdl.handle.net/1963/3930 U1 - 471 U2 - Mathematics U3 - Functional Analysis and Applications ER - TY - RPRT T1 - Stratos: a code for 3D free surface flows with floating constraints Y1 - 2009 A1 - Antonio DeSimone A1 - B. Bianchi A1 - Luca Heltai AB - This report presents a brief discussion of the theoretical aspects and practical implementation of STRATOS . STRATOS is a 3D code for the simulation\\nof hydrodynamic flows for incompressible fluids, in the presence of a free surface, capable of simulating the interaction between the free surface and a\\nfloating object via Lagrange multipliers...... UR - http://hdl.handle.net/1963/3701 U1 - 604 U2 - Mathematics U3 - Functional Analysis and Applications ER - TY - RPRT T1 - Tools for the Solution of PDEs Defined on Curved Manifolds with deal.II Y1 - 2009 A1 - Antonio DeSimone A1 - Luca Heltai A1 - Cataldo Manigrasso AB - The deal.II finite element library was originally designed to solve partial differential equations defined on one, two or three space dimensions, mostly\\nvia the Finite Element Method. In its versions prior to version 6.2, the user could not solve problems defined on curved manifolds embedded in two or\\nthree spacial dimensions. This infrastructure is needed if one wants to solve, for example, Boundary Integral Equations. UR - http://hdl.handle.net/1963/3700 U1 - 605 U2 - Mathematics U3 - Functional Analysis and Applications ER -