@unpublished {2023, title = {Applicable Methodologies for the Mass Transfer Phenomenon in Tumble Dryers: A Review}, year = {2023}, abstract = {

Tumble dryers offer a fast and convenient way of drying textiles independent of weather conditions and therefore are frequently used in ordinary households. However, artificial drying of textiles consumes considerable amounts of energy, approximately 8.2 percent of the residential electricity consumption is for drying of textiles in northern European countries (Cranston et al., 2019). Several authors have investigated the aspects of the clothes drying cycle with experimental and numerical methods to understand and improve the process. The first turning point study on understanding the physics of evaporation for tumble dryers was presented by Lambert et al. (1991) in the early 90s. With the aid of Chilton_Colburn analogy, they introduced the concept of area-mass transfer coefficient to address evaporation rate. Afterwards, several experimental or numerical studies were published based on this concept, and furthermore, the model was then developed into 0-dimensional (Deans, 2001) and 1-dimensional (Wei et al., 2017) to gain more accuracy. The evaporation rate is considered to be the main system parameter for dryers with which other performance parameters including drying time, effectiveness, moisture content and efficiency can be estimated. More recent literature focused on utilizing dimensional analysis or image processing techniques to correlate drying indices with system parameters. However, the validity of these regressed models is machine-specific, and hence, cannot be generalized yet. All the previous models for estimating the evaporation rate in tumble dryers are discussed. The review of the related literature showed that all of the previous models for the prediction of the evaporation rate in the clothes dryers have some limitations in terms of accuracy and applicability.

}, author = {Sajad Salavatidezfouli and Sajad Hajisharifi and Michele Girfoglio and Giovanni Stabile and Gianluigi Rozza} } @article {2022, title = {A comparison of reduced-order modeling approaches using artificial neural networks for PDEs with bifurcating solutions}, journal = {ETNA - Electronic Transactions on Numerical Analysis}, volume = {56}, year = {2022}, pages = {52{\textendash}65}, doi = {10.1553/etna_vol56s52}, author = {Martin W. Hess and Annalisa Quaini and Gianluigi Rozza} } @booklet {2022, title = {Data-Driven Enhanced Model Reduction for Bifurcating Models in Computational Fluid Dynamics}, year = {2022}, author = {Martin W. Hess and Annalisa Quaini and Gianluigi Rozza} } @booklet {2022, title = {A Data-Driven Surrogate Modeling Approach for Time-Dependent Incompressible Navier-Stokes Equations with Dynamic Mode Decomposition and Manifold Interpolation}, year = {2022}, author = {Martin W. Hess and Annalisa Quaini and Gianluigi Rozza} } @article {PeschkaHeltai-2022-a, title = {Model hierarchies and higher-order discretisation of time-dependent thin-film free boundary problems with dynamic contact angle}, journal = {Journal of Computational Physics}, volume = {464}, year = {2022}, pages = {111325}, author = {Dirk Peschka and Luca Heltai} } @unpublished {2022, title = {Model Reduction Using Sparse Polynomial Interpolation for the Incompressible Navier-Stokes Equations}, year = {2022}, author = {Martin W. Hess and Gianluigi Rozza} } @article {ArndtBangerthFeder-2022-a, title = {The \textttdeal.II Library, Version 9.4}, journal = {Journal of Numerical Mathematics}, year = {2022}, author = {Daniel Arndt and Wolfgang Bangerth Marco Feder and Marc Fehling and Rene Gassm{\"o}ller and Timo Heister and Luca Heltai and Martin Kronbichler and Matthias Maier and Peter Munch and Jean-Paul Pelteret and Simon Sticko and Bruno Turcksin and David Wells} } @article {PeschkaZafferiHeltai-2022-a, title = {Variational Approach to Fluid{\textendash}Structure Interaction via GENERIC}, journal = {Journal of Non-Equilibrium Thermodynamics}, year = {2022}, author = {Dirk Peschka and Andrea Zafferi and Luca Heltai and Marita Thomas} } @article {LiuMcBrideSaxena-2022-a, title = {Vibration Analysis of Piezoelectric Kirchhoff-Love shells based on Catmull-Clark Subdivision Surfaces}, journal = {International Journal for Numerical Methods in Engineering}, year = {2022}, author = {Zhaowei Liu and Andrew McBride and Prashant Saxena and Luca Heltai and Yilin Qu and Paul Steinmann} } @unpublished {2021, title = {An artificial neural network approach to bifurcating phenomena in computational fluid dynamics}, year = {2021}, author = {Federico Pichi and Francesco Ballarin and Gianluigi Rozza and Jan S Hesthaven} } @article {ArndtBangerthBlais-2021-h, title = {The deal.II Library, Version 9.3}, journal = {Journal of Numerical Mathematics}, year = {2021}, url = {https://doi.org/10.1515/jnma-2021-0081}, author = {Daniel Arndt and Wolfgang Bangerth and Bruno Blais and Marc Fehling and Rene Gassm{\"o}ller and Timo Heister and Luca Heltai and Uwe K{\"o}cher and Martin Kronbichler and Matthias Maier and Peter Munch and Jean-Paul Pelteret and Sebastian Proell and Konrad Simon and Bruno Turcksin and David Wells and Jiaqi Zhang} } @conference {2021, title = {Discontinuous Galerkin Model Order Reduction of Geometrically Parametrized Stokes Equation}, booktitle = {Numerical Mathematics and Advanced Applications ENUMATH 2019}, year = {2021}, month = {2021//}, publisher = {Springer International Publishing}, organization = {Springer International Publishing}, address = {Cham}, abstract = {

The present work focuses on the geometric parametrization and the reduced order modeling of the Stokes equation. We discuss the concept of a parametrized geometry and its application within a reduced order modeling technique. The full order model is based on the discontinuous Galerkin method with an interior penalty formulation. We introduce the broken Sobolev spaces as well as the weak formulation required for an affine parameter dependency. The operators are transformed from a fixed domain to a parameter dependent domain using the affine parameter dependency. The proper orthogonal decomposition is used to obtain the basis of functions of the reduced order model. By using the Galerkin projection the linear system is projected onto the reduced space. During this process, the offline-online decomposition is used to separate parameter dependent operations from parameter independent operations. Finally this technique is applied to an obstacle test problem.The numerical outcomes presented include experimental error analysis, eigenvalue decay and measurement of online simulation time.

}, isbn = {978-3-030-55874-1}, author = {Nirav Shah and Martin W. Hess and Gianluigi Rozza}, editor = {Vermolen, Fred J. and Vuik, Cornelis} } @article {2021, title = {Efficient computation of bifurcation diagrams with a deflated approach to reduced basis spectral element method}, journal = {Advances in Computational Mathematics}, volume = {47}, year = {2021}, abstract = {

The majority of the most common physical phenomena can be described using partial differential equations (PDEs). However, they are very often characterized by strong nonlinearities. Such features lead to the coexistence of multiple solutions studied by the bifurcation theory. Unfortunately, in practical scenarios, one has to exploit numerical methods to compute the solutions of systems of PDEs, even if the classical techniques are usually able to compute only a single solution for any value of a parameter when more branches exist. In this work, we implemented an elaborated deflated continuation method that relies on the spectral element method (SEM) and on the reduced basis (RB) one to efficiently compute bifurcation diagrams with more parameters and more bifurcation points. The deflated continuation method can be obtained combining the classical continuation method and the deflation one: the former is used to entirely track each known branch of the diagram, while the latter is exploited to discover the new ones. Finally, when more than one parameter is considered, the efficiency of the computation is ensured by the fact that the diagrams can be computed during the online phase while, during the offline one, one only has to compute one-dimensional diagrams. In this work, after a more detailed description of the method, we will show the results that can be obtained using it to compute a bifurcation diagram associated with a problem governed by the Navier-Stokes equations.\ 

}, doi = {10.1007/s10444-020-09827-6}, author = {Moreno Pintore and Federico Pichi and Martin W. Hess and Gianluigi Rozza and Claudio Canuto} } @article {HeltaiCaiazzoMueller-2021-a, title = {Multiscale coupling of one-dimensional vascular models and elastic tissues}, journal = {Annals of Biomedical Engineering}, year = {2021}, author = {Luca Heltai and Alfonso Caiazzo and Lucas M{\"u}eller} } @article {HeltaiBangerthKronbichler-2021, title = {Propagating geometry information to finite element computations}, journal = {Transactions on Mathematical Software}, volume = {47}, year = {2021}, pages = {1--30}, chapter = {1}, doi = {https://dx.doi.org/10.1145/3468428}, author = {Luca Heltai and Wolfgang Bangerth and Martin Kronbichler and Andrea Mola} } @booklet {HunJexLan20-Potential, title = {Quantum Systems at The Brink. Existence and Decay Rates of Bound States at Thresholds; Critical Potentials and dimensionality}, howpublished = {arXiv:2107.14128}, year = {2021}, pages = {8}, author = {Dirk Hundertmark and Michal Jex and Markus Lange} } @article {MulitaGianiHeltai-2021-a, title = {Quasi-optimal mesh sequence construction through Smoothed Adaptive Finite Element Methods}, journal = {SIAM Journal on Scientific Computing}, year = {2021}, author = {Ornela Mulita and Stefano Giani and Luca Heltai} } @article {GianiGrubisicHeltai-2021-a, title = {Smoothed-adaptive perturbed inverse iteration for elliptic eigenvalue problems}, journal = {Computational Methods in Applied Mathematics}, volume = {21}, number = {2}, year = {2021}, pages = {385-405}, author = {Stefano Giani and Luka Grubisic and Luca Heltai and Ornela Mulita} } @inbook {RozzaHessStabileTezzeleBallarin2020, title = {Basic ideas and tools for projection-based model reduction of parametric partial differential equations}, booktitle = {Model Order Reduction, Volume 2 Snapshot-Based Methods and Algorithms}, year = {2020}, pages = {1 - 47}, publisher = {De Gruyter}, organization = {De Gruyter}, address = {Berlin, Boston}, isbn = {9783110671490}, doi = {https://doi.org/10.1515/9783110671490-001}, url = {https://www.degruyter.com/view/book/9783110671490/10.1515/9783110671490-001.xml}, author = {Gianluigi Rozza and Martin W. Hess and Giovanni Stabile and Marco Tezzele and F. Ballarin} } @article {han2020gauge, title = {On coherent Hopf 2-algebras}, year = {2020}, month = {05/2020}, abstract = {We construct a coherent Hopf 2-algebra as quantization of a coherent 2-group, which consists of two Hopf coquasigroups and a coassociator. For this constructive method, if we replace Hopf coquasigroups by Hopf algebras, we can construct a strict Hoft 2-algebra, which is a quantisation of 2-group. We also study the crossed comodule of Hopf algebras, which is shown to be a strict Hopf 2-algebra under some conditions. As an example, a quasi coassociative Hopf coquasigroup is employed to build a special coherent Hopf 2-algebra with nontrivial coassociator. Following this we study functions on Cayley algebra basis.}, url = {https://arxiv.org/abs/2005.11207}, author = {Xiao Han} } @article {2020, title = {Data-driven POD-Galerkin reduced order model for turbulent flows}, journal = {Journal of Computational Physics}, volume = {416}, year = {2020}, pages = {109513}, abstract = {

In this work we present a Reduced Order Model which is specifically designed to deal with turbulent flows in a finite volume setting. The method used to build the reduced order model is based on the idea of merging/combining projection-based techniques with data-driven reduction strategies. In particular, the work presents a mixed strategy that exploits a data-driven reduction method to approximate the eddy viscosity solution manifold and a classical POD-Galerkin projection approach for the velocity and the pressure fields, respectively. The newly proposed reduced order model has been validated on benchmark test cases in both steady and unsteady settings with Reynolds up to $Re=O(10^5)$.

}, doi = {10.1016/j.jcp.2020.109513}, url = {https://arxiv.org/abs/1907.09909}, author = {Saddam Hijazi and Giovanni Stabile and Andrea Mola and Gianluigi Rozza} } @article {ArndtBangerthDavydov-2020-a, title = {The deal.II finite element library: Design, features, and insights}, journal = {Computers and Mathematics with Applications}, year = {2020}, url = {https://doi.org/10.1016/j.camwa.2020.02.022}, author = {Daniel Arndt and Wolfgang Bangerth and Denis Davydov and Timo Heister and Luca Heltai and Martin Kronbichler and Matthias Maier and Jean-Paul Pelteret and Bruno Turcksin and David Wells} } @article {ArndtBangerthBlais-2020-a, title = {The deal.II library, Version 9.2}, journal = {Journal of Numerical Mathematics}, volume = {28}, number = {3}, year = {2020}, pages = {131{\textendash}146}, author = {Daniel Arndt and Wolfgang Bangerth and Bruno Blais and Thomas C. Clevenger and Marc Fehling and Alexander V. Grayver and Timo Heister and Luca Heltai and Martin Kronbichler and Matthias Maier and Peter Munch and Jean-Paul Pelteret and Reza Rastak and Ignacio Tomas and Bruno Turcksin and Zhuoran Wang and David Wells} } @article {13850, title = {Efficient computation of bifurcation diagrams with a deflated approach to reduced basis spectral element method}, journal = {Advances in Computational Mathematics}, year = {2020}, abstract = {

The majority of the most common physical phenomena can be described using partial differential equations (PDEs). However, they are very often characterized by strong nonlinearities. Such features lead to the coexistence of multiple solutions studied by the bifurcation theory. Unfortunately, in practical scenarios, one has to exploit numerical methods to compute the solutions of systems of PDEs, even if the classical techniques are usually able to compute only a single solution for any value of a parameter when more branches exist. In this work we implemented an elaborated deflated continuation method, that relies on the spectral element method (SEM) and on the reduced basis (RB) one, to efficiently compute bifurcation diagrams with more parameters and more bifurcation points. The deflated continuation method can be obtained combining the classical continuation method and the deflation one: the former is used to entirely track each known branch of the diagram, while the latter is exploited to discover the new ones. Finally, when more than one parameter is considered, the efficiency of the computation is ensured by the fact that the diagrams can be computed during the online phase while, during the offline one, one only has to compute one-dimensional diagrams. In this work, after a more detailed description of the method, we will show the results that can be obtained using it to compute a bifurcation diagram associated with a problem governed by the Navier-Stokes equations.

}, url = {https://arxiv.org/abs/1912.06089}, author = {Moreno Pintore and Federico Pichi and Martin W. Hess and Gianluigi Rozza and Claudio Canuto} } @conference {HijaziAliStabileBallarinRozza2020, title = {The Effort of Increasing Reynolds Number in Projection-Based Reduced Order Methods: from Laminar to Turbulent Flows}, booktitle = {Lecture Notes in Computational Science and Engineering}, year = {2020}, pages = {245{\textendash}264}, publisher = {Springer International Publishing}, organization = {Springer International Publishing}, address = {Cham}, abstract = {

We present in this double contribution two different reduced order strategies for incompressible parameterized Navier-Stokes equations characterized by varying Reynolds numbers. The first strategy deals with low Reynolds number (laminar flow) and is based on a stabilized finite element method during the offline stage followed by a Galerkin projection on reduced basis spaces generated by a greedy algorithm. The second methodology is based on a full order finite volume discretization. The latter methodology will be used for flows with moderate to high Reynolds number characterized by turbulent patterns. For the treatment of the mentioned turbulent flows at the reduced order level, a new POD-Galerkin approach is proposed. The new approach takes into consideration the contribution of the eddy viscosity also during the online stage and is based on the use of interpolation. The two methodologies are tested on classic benchmark test cases.

}, isbn = {978-3-030-30705-9}, doi = {10.1007/978-3-030-30705-9_22}, author = {Saddam Hijazi and Shafqat Ali and Giovanni Stabile and F. Ballarin and Gianluigi Rozza} } @article {han2020gauge, title = {On the gauge group of Galois objects}, year = {2020}, month = {03/2020}, abstract = {We study the Ehresmann--Schauenburg bialgebroid of a noncommutative principal bundle as a quantization of the classical gauge groupoid of a principal bundle. When the base algebra is in the centre of the total space algebra, the gauge group of the noncommutative principal bundle is isomorphic to the group of bisections of the bialgebroid. In particular we consider Galois objects (non-trivial noncommutative bundles over a point in a sense) for which the bialgebroid is a Hopf algebra. For these we give a crossed module structure for the bisections and the automorphisms of the bialgebroid. Examples include Galois objects of group Hopf algebras and of Taft algebras.}, url = {https://arxiv.org/abs/2002.06097}, author = {Xiao Han and Giovanni Landi} } @article {mine-02-02-011, title = {MicroMotility: State of the art, recent accomplishments and perspectives on the mathematical modeling of bio-motility at microscopic scales}, journal = {Mathematics in Engineering}, volume = {2}, number = {mine-02-02-01}, year = {2020}, pages = {230}, type = {Perspective}, abstract = {

Mathematical modeling and quantitative study of biological motility (in particular, of motility at microscopic scales) is producing new biophysical insight and is offering opportunities for new discoveries at the level of both fundamental science and technology. These range from the explanation of how complex behavior at the level of a single organism emerges from body architecture, to the understanding of collective phenomena in groups of organisms and tissues, and of how these forms of swarm intelligence can be controlled and harnessed in engineering applications, to the elucidation of processes of fundamental biological relevance at the cellular and sub-cellular level. In this paper, some of the most exciting new developments in the fields of locomotion of unicellular organisms, of soft adhesive locomotion across scales, of the study of pore translocation properties of knotted DNA, of the development of synthetic active solid sheets, of the mechanics of the unjamming transition in dense cell collectives, of the mechanics of cell sheet folding in volvocalean algae, and of the self-propulsion of topological defects in active matter are discussed. For each of these topics, we provide a brief state of the art, an example of recent achievements, and some directions for future research.

}, keywords = {active matter, adhesive locomotion, cell motility, cell sheet folding, knotted DNA, topological defects, unicellular swimmers, unjamming transition}, issn = {2640-3501}, doi = {10.3934/mine.2020011}, url = {http://dx.doi.org/10.3934/mine.2020011}, author = {Daniele Agostinelli and Roberto Cerbino and Del Alamo, Juan C and Antonio DeSimone and Stephanie H{\"o}hn and Cristian Micheletti and Giovanni Noselli and Eran Sharon and Julia Yeomans} } @booklet {2020, title = {MicroROM: An Efficient and Accurate Reduced Order Method to Solve Many-Query Problems in Micro-Motility}, year = {2020}, keywords = {FOS: Mathematics, Numerical Analysis (math.NA)}, doi = {10.48550/ARXIV.2006.13836}, url = {https://arxiv.org/abs/2006.13836}, author = {Nicola Giuliani and Martin W. Hess and Antonio DeSimone and Gianluigi Rozza} } @article {AlzettaHeltai-2020-a, title = {Multiscale modeling of fiber reinforced materials via non-matching immersed methods}, journal = {Computers \& Structures}, volume = {239}, year = {2020}, note = {To appear}, pages = {106334}, doi = {10.1016/j.compstruc.2020.106334}, author = {Giovanni Alzetta and Luca Heltai} } @inbook {HijaziStabileMolaRozza2020a, title = {Non-intrusive Polynomial Chaos Method Applied to Full-Order and Reduced Problems in Computational Fluid Dynamics: A Comparison and Perspectives}, booktitle = {Quantification of Uncertainty: Improving Efficiency and Technology: QUIET selected contributions}, year = {2020}, pages = {217{\textendash}240}, publisher = {Springer International Publishing}, organization = {Springer International Publishing}, address = {Cham}, abstract = {

In this work, Uncertainty Quantification (UQ) based on non-intrusive Polynomial Chaos Expansion (PCE) is applied to the CFD problem of the flow past an airfoil with parameterized angle of attack and inflow velocity. To limit the computational cost associated with each of the simulations required by the non-intrusive UQ algorithm used, we resort to a Reduced Order Model (ROM) based on Proper Orthogonal Decomposition (POD)-Galerkin approach. A first set of results is presented to characterize the accuracy of the POD-Galerkin ROM developed approach with respect to the Full Order Model (FOM) solver (OpenFOAM). A further analysis is then presented to assess how the UQ results are affected by substituting the FOM predictions with the surrogate ROM ones.

}, isbn = {978-3-030-48721-8}, doi = {10.1007/978-3-030-48721-8_10}, url = {https://doi.org/10.1007/978-3-030-48721-8_10}, author = {Saddam Hijazi and Giovanni Stabile and Andrea Mola and Gianluigi Rozza} } @article {Dal-MasoHeltai-2020-a, title = {A numerical study of the jerky crack growth in elastoplastic materials with localized plasticity}, journal = {Journal of Convex Analysis}, year = {2020}, note = {To appear}, url = {https://arxiv.org/abs/2004.12705}, author = {Gianni Dal Maso and Luca Heltai} } @article {HL2020, title = {A priori error estimates of regularized elliptic problems}, journal = {Numerische Mathematik}, volume = {146}, number = {3}, year = {2020}, pages = {571{\textendash}596}, abstract = {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.}, isbn = {0945-3245}, doi = {10.1007/s00211-020-01152-w}, url = {https://doi.org/10.1007/s00211-020-01152-w}, author = {Luca Heltai and Wenyu Lei} } @article {2020, title = {A priori error estimates of regularized elliptic problems}, journal = {Numerische Mathematik}, year = {2020}, author = {Luca Heltai and Wenyu Lei} } @booklet {HunJexLan-PhysRev, title = {Quantum Systems at The Brink: Properties of Atomic Bound States at The Ionization Threshold}, year = {2020}, author = {Dirk Hundertmark and Michal Jex and Markus Lange} } @article {2020, title = {Reduced Basis Model Order Reduction for Navier-Stokes equations in domains with walls of varying curvature}, journal = {International Journal of Computational Fluid Dynamics}, volume = {34}, year = {2020}, pages = {119-126}, abstract = {

We consider the Navier-Stokes equations in a channel with a narrowing and walls of varying curvature. By applying the empirical interpolation method to generate an affine parameter dependency, the offline-online procedure can be used to compute reduced order solutions for parameter variations. The reduced order space is computed from the steady-state snapshot solutions by a standard POD procedure. The model is discretised with high-order spectral element ansatz functions, resulting in 4752 degrees of freedom. The proposed reduced order model produces accurate approximations of steady-state solutions for a wide range of geometries and kinematic viscosity values. The application that motivated the present study is the onset of asymmetries (i.e., symmetry breaking bifurcation) in blood flow through a regurgitant mitral valve, depending on the Reynolds number and the valve shape. Through our computational study, we found that the critical Reynolds number for the symmetry breaking increases as the wall curvature increases.

}, doi = {10.1080/10618562.2019.1645328}, url = {https://arxiv.org/abs/1901.03708}, author = {Martin W. Hess and Annalisa Quaini and Gianluigi Rozza} } @article {2020, title = {Reduced basis model order reduction for Navier{\textendash}Stokes equations in domains with walls of varying curvature}, journal = {International Journal of Computational Fluid Dynamics}, volume = {34}, year = {2020}, pages = {119-126}, abstract = {

We consider the Navier{\textendash}Stokes equations in a channel with a narrowing and walls of varying curvature. By applying the empirical interpolation method to generate an affine parameter dependency, the offline-online procedure can be used to compute reduced order solutions for parameter variations. The reduced-order space is computed from the steady-state snapshot solutions by a standard POD procedure. The model is discretised with high-order spectral element ansatz functions, resulting in 4752 degrees of freedom. The proposed reduced-order model produces accurate approximations of steady-state solutions for a wide range of geometries and kinematic viscosity values. The application that motivated the present study is the onset of asymmetries (i.e. symmetry breaking bifurcation) in blood flow through a regurgitant mitral valve, depending on the Reynolds number and the valve shape. Through our computational study, we found that the critical Reynolds number for the symmetry breaking increases as the wall curvature increases.

}, doi = {10.1080/10618562.2019.1645328}, url = {https://www.scopus.com/inward/record.uri?eid=2-s2.0-85085233294\&doi=10.1080\%2f10618562.2019.1645328\&partnerID=40\&md5=e2ed8f24c66376cdc8b5485aa400efb0}, author = {Martin W. Hess and Annalisa Quaini and Gianluigi Rozza} } @conference {2020, title = {A Reduced Order Approach for the Embedded Shifted Boundary FEM and a Heat Exchange System on Parametrized Geometries}, booktitle = {IUTAM Symposium on Model Order Reduction of Coupled Systems, Stuttgart, Germany, May 22{\textendash}25, 2018}, year = {2020}, publisher = {Springer International Publishing}, organization = {Springer International Publishing}, abstract = {

A model order reduction technique is combined with an embedded boundary finite element method with a POD-Galerkin strategy. The proposed methodology is applied to parametrized heat transfer problems and we rely on a sufficiently refined shape-regular background mesh to account for parametrized geometries. In particular, the employed embedded boundary element method is the Shifted Boundary Method (SBM) recently proposed. This approach is based on the idea of shifting the location of true boundary conditions to a surrogate boundary, with the goal of avoiding cut cells near the boundary of the computational domain. This combination of methodologies has multiple advantages. In the first place, since the Shifted Boundary Method always relies on the same background mesh, there is no need to update the discretized parametric domain. Secondly, we avoid the treatment of cut cell elements, which usually need particular attention. Thirdly, since the whole background mesh is considered in the reduced basis construction, the SBM allows for a smooth transition of the reduced modes across the immersed domain boundary. The performances of the method are verified in two dimensional heat transfer numerical examples.

}, doi = {10.1007/978-3-030-21013-7_8}, url = {https://arxiv.org/abs/1807.07753}, author = {Efthymios N Karatzas and Giovanni Stabile and Nabib Atallah and Guglielmo Scovazzi and Gianluigi Rozza}, editor = {Fehr, J{\"o}rg and Bernard Haasdonk} } @article {2020, title = {A spectral element reduced basis method for navier{\textendash}stokes equations with geometric variations}, journal = {Lecture Notes in Computational Science and Engineering}, volume = {134}, year = {2020}, pages = {561-571}, abstract = {

We consider the Navier-Stokes equations in a channel with a narrowing of varying height. The model is discretized with high-order spectral element ansatz functions, resulting in 6372 degrees of freedom. The steady-state snapshot solutions define a reduced order space through a standard POD procedure. The reduced order space allows to accurately and efficiently evaluate the steady-state solutions for different geometries. In particular, we detail different aspects of implementing the reduced order model in combination with a spectral element discretization. It is shown that an expansion in element-wise local degrees of freedom can be combined with a reduced order modelling approach to enhance computational times in parametric many-query scenarios.

}, doi = {10.1007/978-3-030-39647-3_45}, author = {Martin W. Hess and Annalisa Quaini and Gianluigi Rozza} } @article {han2020gauge, title = {Twisted Ehresmann Schauenburg bialgebroids}, year = {2020}, month = {09/2020}, abstract = {We construct an invertible normalised 2 cocycle on the Ehresmann Schauenburg bialgebroid of a cleft Hopf Galois extension under the condition that the corresponding Hopf algebra is cocommutative and the image of the unital cocycle corresponding to this cleft Hopf Galois extension belongs to the centre of the coinvariant subalgebra. Moreover, we show that any Ehresmann Schauenburg bialgebroid of this kind is isomorphic to a 2-cocycle twist of the Ehresmann Schauenburg bialgebroid corresponding to a Hopf Galois extension without cocycle, where comodule algebra is an ordinary smash product of the coinvariant subalgebra and the Hopf algebra (i.e. $\C(B/$\#$_{\sigma}H, H)\simeq \C(B\#H, H)^{\tilde{\sigma}}$). We also study the theory in the case of a Galois object where the base is trivial but without requiring the Hopf algebra to be cocommutative.}, url = {https://arxiv.org/abs/2009.02764}, author = {Xiao Han} } @article {ArndtBangerthClevenger-2019-a, title = {The deal.II Library, Version 9.1}, journal = {Journal of Numerical Mathematics}, year = {2019}, author = {Daniel Arndt and Wolfgang Bangerth and Thomas C. Clevenger and Denis Davydov and Marc Fehling and Garcia-Sanchez, Daniel and Harper, Graham and Timo Heister and Luca Heltai and Martin Kronbichler and Maguire Kynch, Ross and Matthias Maier and Jean-Paul Pelteret and Bruno Turcksin and David Wells} } @article {ArndtBangerthClevenger-2019-a, title = {The deal.II Library, Version 9.1}, journal = {Journal of Numerical Mathematics}, year = {2019}, abstract = {This paper provides an overview of the new features of the finite element library deal.II, version 9.1.}, issn = {15702820}, doi = {10.1515/jnma-2019-0064}, author = {Daniel Arndt and Wolfgang Bangerth and Thomas C. Clevenger and Denis Davydov and Marc Fehling and Garcia-Sanchez, Daniel and Harper, Graham and Timo Heister and Luca Heltai and Martin Kronbichler and Maguire Kynch, Ross and Matthias Maier and Jean-Paul Pelteret and Bruno Turcksin and David Wells} } @article {HeltaiRotundo-2019-a, title = {Error estimates in weighted Sobolev norms for finite element immersed interface methods}, journal = {Computers \& Mathematics with Applications}, volume = {78}, number = {11}, year = {2019}, pages = {3586{\textendash}3604}, publisher = {Elsevier BV}, doi = {10.1016/j.camwa.2019.05.029}, url = {https://doi.org/10.1016/j.camwa.2019.05.029}, author = {Luca Heltai and Nella Rotundo} } @article {2019, title = {A Localized Reduced-Order Modeling Approach for PDEs with Bifurcating Solutions}, journal = {Computer Methods in Applied Mechanics and Engineering}, volume = {351}, year = {2019}, pages = {379-403}, abstract = {

Reduced-order modeling (ROM) commonly refers to the construction, based on a few solutions (referred to as snapshots) of an expensive discretized partial differential equation (PDE), and the subsequent application of low-dimensional discretizations of partial differential equations (PDEs) that can be used to more efficiently treat problems in control and optimization, uncertainty quantification, and other settings that require multiple approximate PDE solutions. In this work, a ROM is developed and tested for the treatment of nonlinear PDEs whose solutions bifurcate as input parameter values change. In such cases, the parameter domain can be subdivided into subregions, each of which corresponds to a different branch of solutions. Popular ROM approaches such as proper orthogonal decomposition (POD), results in a global low-dimensional basis that does no respect not take advantage of the often large differences in the PDE solutions corresponding to different subregions. Instead, in the new method, the k-means algorithm is used to cluster snapshots so that within cluster snapshots are similar to each other and are dissimilar to those in other clusters. This is followed by the construction of local POD bases, one for each cluster. The method also can detect which cluster a new parameter point belongs to, after which the local basis corresponding to that cluster is used to determine a ROM approximation. Numerical experiments show the effectiveness of the method both for problems for which bifurcation cause continuous and discontinuous changes in the solution of the PDE.

}, doi = {10.1016/j.cma.2019.03.050}, url = {https://arxiv.org/abs/1807.08851}, author = {Martin W. Hess and Alla, Alessandro and Annalisa Quaini and Gianluigi Rozza and Max Gunzburger} } @article {2019, title = {A localized reduced-order modeling approach for PDEs with bifurcating solutions}, journal = {Computer Methods in Applied Mechanics and Engineering}, volume = {351}, year = {2019}, pages = {379-403}, abstract = {

Reduced-order modeling (ROM) commonly refers to the construction, based on a few solutions (referred to as snapshots) of an expensive discretized partial differential equation (PDE), and the subsequent application of low-dimensional discretizations of partial differential equations (PDEs) that can be used to more efficiently treat problems in control and optimization, uncertainty quantification, and other settings that require multiple approximate PDE solutions. Although ROMs have been successfully used in many settings, ROMs built specifically for the efficient treatment of PDEs having solutions that bifurcate as the values of input parameters change have not received much attention. In such cases, the parameter domain can be subdivided into subregions, each of which corresponds to a different branch of solutions. Popular ROM approaches such as proper orthogonal decomposition (POD), results in a global low-dimensional basis that does not respect the often large differences in the PDE solutions corresponding to different subregions. In this work, we develop and test a new ROM approach specifically aimed at bifurcation problems. In the new method, the k-means algorithm is used to cluster snapshots so that within cluster snapshots are similar to each other and are dissimilar to those in other clusters. This is followed by the construction of local POD bases, one for each cluster. The method also can detect which cluster a new parameter point belongs to, after which the local basis corresponding to that cluster is used to determine a ROM approximation. Numerical experiments show the effectiveness of the method both for problems for which bifurcation cause continuous and discontinuous changes in the solution of the PDE.

}, doi = {10.1016/j.cma.2019.03.050}, url = {https://www.scopus.com/inward/record.uri?eid=2-s2.0-85064313505\&doi=10.1016\%2fj.cma.2019.03.050\&partnerID=40\&md5=8b095034b9e539995facc7ce7bafa9e9}, author = {Martin W. Hess and Alla, Alessandro and Annalisa Quaini and Gianluigi Rozza and Max Gunzburger} } @article {HeltaiCaiazzo-2018-a, title = {Multiscale modeling of vascularized tissues via non-matching immersed methods}, journal = {International Journal for Numerical Methods in Biomedical Engineering}, volume = {35}, number = {12}, year = {2019}, pages = {e3264}, url = {https://doi.org/10.1002/cnm.3264}, author = {Luca Heltai and Alfonso Caiazzo} } @booklet {HunJexLan19-Helium, title = {Quantum Systems at The Brink. Existence and Decay Rates of Bound States at Thresholds; Helium}, howpublished = {arXiv:1908.04883}, year = {2019}, pages = {25}, author = {Dirk Hundertmark and Michal Jex and Markus Lange} } @booklet {HunJexLan19-NAtom, title = {Quantum Systems at The Brink. Existence and Decay Rates of Bound States at Thresholds; Atoms}, howpublished = {arXiv:1908.05016}, year = {2019}, pages = {14}, author = {Dirk Hundertmark and Michal Jex and Markus Lange} } @booklet {2019, title = {Regularity of minimizers for a model of charged droplets}, year = {2019}, author = {Guido De Philippis and Jonas Hirsch and Giulia Vescovo} } @inbook {2019, title = {A Spectral Element Reduced Basis Method in Parametric CFD}, booktitle = {Numerical Mathematics and Advanced Applications - ENUMATH 2017}, volume = {126}, year = {2019}, publisher = {Springer International Publishing}, organization = {Springer International Publishing}, chapter = {A Spectral Element Reduced Basis Method in Parametric CFD}, abstract = {

We consider the Navier-Stokes equations in a channel with varying Reynolds numbers. The model is discretized with high-order spectral element ansatz functions, resulting in 14 259 degrees of freedom. The steady-state snapshot solu- tions define a reduced order space, which allows to accurately evaluate the steady- state solutions for varying Reynolds number with a reduced order model within a fixed-point iteration. In particular, we compare different aspects of implementing the reduced order model with respect to the use of a spectral element discretization. It is shown, how a multilevel static condensation in the pressure and velocity boundary degrees of freedom can be combined with a reduced order modelling approach to enhance computational times in parametric many-query scenarios.

}, doi = {10.1007/978-3-319-96415-7_64 pages = 693{\textendash}701}, url = {https://arxiv.org/abs/1712.06432}, author = {Martin W. Hess and Gianluigi Rozza}, editor = {Radu, Florin Adrian and Kumar, Kundan and Berre, Inga and Nordbotten, Jan Martin and Pop, Iuliu Sorin} } @article {2019, title = {A spectral element reduced basis method in parametric CFD}, journal = {Lecture Notes in Computational Science and Engineering}, volume = {126}, year = {2019}, pages = {693-701}, abstract = {

We consider the Navier-Stokes equations in a channel with varying Reynolds numbers. The model is discretized with high-order spectral element ansatz functions, resulting in 14,259 degrees of freedom. The steady-state snapshot solutions define a reduced order space, which allows to accurately evaluate the steady-state solutions for varying Reynolds number with a reduced order model within a fixed-point iteration. In particular, we compare different aspects of implementing the reduced order model with respect to the use of a spectral element discretization. It is shown, how a multilevel static condensation (Karniadakis and Sherwin, Spectral/hp element methods for computational fluid dynamics, 2nd edn. Oxford University Press, Oxford, 2005) in the pressure and velocity boundary degrees of freedom can be combined with a reduced order modelling approach to enhance computational times in parametric many-query scenarios.

}, doi = {10.1007/978-3-319-96415-7_64}, url = {https://www.scopus.com/inward/record.uri?eid=2-s2.0-85060005503\&doi=10.1007\%2f978-3-319-96415-7_64\&partnerID=40\&md5=d1a900db8ddb92cd818d797ec212a4c6}, author = {Martin W. Hess and Gianluigi Rozza} } @article {20.500.11767_81643, title = {Accelerating the iterative solution of convection-diffusion problems using singular value decomposition}, journal = {NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS}, year = {2018}, pages = {1{\textendash}21}, doi = {10.1002/nla.2211}, url = {https://arxiv.org/abs/1807.09467}, author = {Giuseppe Pitton and Luca Heltai} } @article {BraHasLan18, title = {On asymptotic expansions in spin-boson models}, journal = {Ann. Henri Poincar{\'e}}, volume = {19}, number = {2}, year = {2018}, pages = {515{\textendash}564}, issn = {1424-0637}, doi = {10.1007/s00023-017-0625-7}, url = {https://doi.org/10.1007/s00023-017-0625-7}, author = {Br{\"a}unlich, Gerhard and Hasler, David and Markus Lange} } @article {20.500.11767_81647, title = {An authenticated theoretical modeling of electrified fluid jet in core{\textendash}shell nanofibers production}, journal = {JOURNAL OF INDUSTRIAL TEXTILES}, volume = {47}, year = {2018}, pages = {1791{\textendash}1811}, doi = {10.1177/1528083717710711}, author = {Rafiei, S. and Noroozi, B. and Luca Heltai and Haghi, A. K.} } @article {2018, title = {Certified Reduced Basis Approximation for the Coupling of Viscous and Inviscid Parametrized Flow Models}, journal = {Journal of Scientific Computing}, volume = {74}, year = {2018}, pages = {197-219}, doi = {10.1007/s10915-017-0430-y}, url = {https://www.scopus.com/inward/record.uri?eid=2-s2.0-85017156114\&doi=10.1007\%2fs10915-017-0430-y\&partnerID=40\&md5=023ef0bb95713f4442d1fa374c92a964}, author = {Immanuel Martini and Bernard Haasdonk and Gianluigi Rozza} } @article {20.500.11767_83906, title = {deal2lkit: A toolkit library for high performance programming in deal.II}, journal = {SOFTWAREX}, volume = {7}, year = {2018}, pages = {318{\textendash}327}, doi = {10.1016/j.softx.2018.09.004}, author = {Alberto Sartori and Nicola Giuliani and Mauro Bardelloni and Luca Heltai} } @article {20.500.11767_81694, title = {The deal.II Library, Version 9.0}, journal = {JOURNAL OF NUMERICAL MATHEMATICS}, year = {2018}, doi = {10.1515/jnma-2018-0054}, url = {https://doi.org/10.1515/jnma-2018-0054}, author = {Giovanni Alzetta and Daniel Arndt and Wolfgang Bangerth and Boddu, Vishal and Brands, Benjamin and Denis Davydov and Rene Gassm{\"o}ller and Timo Heister and Luca Heltai and Kormann, Katharina and Martin Kronbichler and Matthias Maier and Jean-Paul Pelteret and Bruno Turcksin and David Wells} } @inbook {20.500.11767_86398, title = {A distributed lagrange formulation of the finite element immersed boundary method for fluids interacting with compressible solids}, booktitle = {Mathematical and Numerical Modeling of the Cardiovascular System and Applications}, volume = {16}, year = {2018}, pages = {1{\textendash}21}, publisher = {Springer International Publishing}, organization = {Springer International Publishing}, address = {Cham}, doi = {10.1007/978-3-319-96649-6_1}, url = {https://arxiv.org/abs/1712.02545v1}, author = {Boffi, Daniele and Gastaldi, Lucia and Luca Heltai} } @article {20.500.11767_86219, title = {Iterative map-making with two-level preconditioning for polarized cosmic microwave background data sets. A worked example for ground-based experiments}, journal = {ASTRONOMY \& ASTROPHYSICS}, volume = {618}, year = {2018}, pages = {1{\textendash}14}, doi = {10.1051/0004-6361/201832710}, url = {https://arxiv.org/abs/1801.08937}, author = {Puglisi, Giuseppe and Poletti, Davide and Fabbian, Giulio and Baccigalupi, Carlo and Luca Heltai and Stompor, Radek} } @article {20.500.11767_81737, title = {NURBS-SEM: A hybrid spectral element method on NURBS maps for the solution of elliptic PDEs on surfaces}, journal = {COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING}, volume = {338}, year = {2018}, pages = {440{\textendash}462}, doi = {10.1016/j.cma.2018.04.039}, url = {https://arxiv.org/abs/1804.08271}, author = {Giuseppe Pitton and Luca Heltai} } @article {20.500.11767_81735, title = {Predicting and Optimizing Microswimmer Performance from the Hydrodynamics of Its Components: The Relevance of Interactions}, journal = {SOFT ROBOTICS}, volume = {5}, year = {2018}, pages = {410{\textendash}424}, doi = {10.1089/soro.2017.0099}, url = {https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6094362/}, author = {Nicola Giuliani and Luca Heltai and Antonio DeSimone} } @article {doi:10.1142/S0129055X18500204, title = {Principal fibrations over noncommutative spheres}, journal = {Reviews in Mathematical Physics}, volume = {30}, number = {10}, year = {2018}, pages = {1850020}, abstract = {We present examples of noncommutative four-spheres that are base spaces of $SU(2)$-principal bundles with noncommutative seven-spheres as total spaces. The noncommutative coordinate algebras of the four-spheres are generated by the entries of a projection which is invariant under the action of $SU(2)$. We give conditions for the components of the Connes{\textendash}Chern character of the projection to vanish but the second (the top) one. The latter is then a non-zero Hochschild cycle that plays the role of the volume form for the noncommutative four-spheres.}, doi = {10.1142/S0129055X18500204}, url = {https://arxiv.org/abs/1804.07032}, author = {Michel Dubois-Violette and Xiao Han and Giovanni Landi} } @inbook {2018, title = {Reduced Basis Approximation and A Posteriori Error Estimation: Applications to Elasticity Problems in Several Parametric Settings}, booktitle = {Numerical Methods for PDEs}, volume = {15}, number = {SEMA SIMAI}, year = {2018}, doi = {https://doi.org/10.1007/978-3-319-94676-4_8}, url = {https://link.springer.com/chapter/10.1007/978-3-319-94676-4_8}, author = {Huynh, D. B. P. and Federico Pichi and Gianluigi Rozza} } @article {2018, title = {Reduced Basis Approximation and A Posteriori Error Estimation: Applications to Elasticity Problems in Several Parametric Settings}, journal = {SEMA SIMAI Springer Series}, volume = {15}, year = {2018}, pages = {203-247}, abstract = {

In this work we consider (hierarchical, Lagrange) reduced basis approximation and a posteriori error estimation for elasticity problems in affinely parametrized geometries. The essential ingredients of the methodology are: a Galerkin projection onto a low-dimensional space associated with a smooth {\textquotedblleft}parametric manifold{\textquotedblright}{\textemdash}dimension reduction; an efficient and effective greedy sampling methods for identification of optimal and numerically stable approximations{\textemdash}rapid convergence; an a posteriori error estimation procedures{\textemdash}rigorous and sharp bounds for the functional outputs related with the underlying solution or related quantities of interest, like stress intensity factor; and Offline-Online computational decomposition strategies{\textemdash}minimum marginal cost for high performance in the real-time and many-query (e.g., design and optimization) contexts. We present several illustrative results for linear elasticity problem in parametrized geometries representing 2D Cartesian or 3D axisymmetric configurations like an arc-cantilever beam, a center crack problem, a composite unit cell or a woven composite beam, a multi-material plate, and a closed vessel. We consider different parametrization for the systems: either physical quantities{\textemdash}to model the materials and loads{\textemdash}and geometrical parameters{\textemdash}to model different geometrical configurations{\textemdash}with isotropic and orthotropic materials working in plane stress and plane strain approximation. We would like to underline the versatility of the methodology in very different problems. As last example we provide a nonlinear setting with increased complexity.

}, doi = {10.1007/978-3-319-94676-4_8}, url = {https://www.scopus.com/inward/record.uri?eid=2-s2.0-85055036627\&doi=10.1007\%2f978-3-319-94676-4_8\&partnerID=40\&md5=e9c07038e7bcc6668ec702c0653410dc}, author = {D.B.P. Huynh and Federico Pichi and Gianluigi Rozza} } @article {HasLan18, title = {Renormalization analysis for degenerate ground states}, journal = {J. Funct. Anal.}, volume = {275}, number = {1}, year = {2018}, pages = {103{\textendash}148}, issn = {0022-1236}, doi = {10.1016/j.jfa.2018.03.005}, url = {https://doi.org/10.1016/j.jfa.2018.03.005}, author = {Hasler, David and Markus Lange} } @article {giulianiEtAl2018, title = {π-BEM : A flexible parallel implementation for adaptive , geometry aware , and high order boundary element methods}, journal = {Advances in Engineering Software}, volume = {121}, year = {2018}, pages = {39{\textendash}58}, author = {Nicola Giuliani and Andrea Mola and Luca Heltai} } @article {20.500.11767_47950, title = {The deal.II Library, Version 8.5}, journal = {JOURNAL OF NUMERICAL MATHEMATICS}, volume = {25}, year = {2017}, pages = {137{\textendash}145}, doi = {10.1515/jnma-2017-0058}, url = {https://www.dealii.org/deal85-preprint.pdf}, author = {Daniel Arndt and Wolfgang Bangerth and Denis Davydov and Timo Heister and Luca Heltai and Martin Kronbichler and Matthias Maier and Jean-Paul Pelteret and Bruno Turcksin and David Wells} } @book {MR3729265, title = {$hp$-version discontinuous Galerkin methods on polygonal and polyhedral meshes}, series = {SpringerBriefs in Mathematics}, year = {2017}, pages = {viii+131}, publisher = {Springer, Cham}, organization = {Springer, Cham}, isbn = {978-3-319-67671-5; 978-3-319-67673-9}, author = {Andrea Cangiani and Dong, Zhaonan and E.H. Georgoulis and Houston, Paul} } @inbook {ChinestaHuertaRozzaWillcox2017, title = {Model Reduction Methods}, booktitle = {Encyclopedia of Computational Mechanics Second Edition}, year = {2017}, pages = {1-36}, publisher = {John Wiley \& Sons}, organization = {John Wiley \& Sons}, chapter = {Model Reduction Methods}, abstract = {

This chapter presents an overview of model order reduction {\textendash} a new paradigm in the field of simulation-based engineering sciences, and one that can tackle the challenges and leverage the opportunities of modern ICT technologies. Despite the impressive progress attained by simulation capabilities and techniques, a number of challenging problems remain intractable. These problems are of different nature, but are common to many branches of science and engineering. Among them are those related to high-dimensional problems, problems involving very different time scales, models defined in degenerate domains with at least one of the characteristic dimensions much smaller than the others, model requiring real-time simulation, and parametric models. All these problems represent a challenge for standard mesh-based discretization techniques; yet the ability to solve these problems efficiently would open unexplored routes for real-time simulation, inverse analysis, uncertainty quantification and propagation, real-time optimization, and simulation-based control {\textendash} critical needs in many branches of science and engineering. Model order reduction offers new simulation alternatives by circumventing, or at least alleviating, otherwise intractable computational challenges. In the present chapter, we revisit three of these model reduction techniques: proper orthogonal decomposition, proper generalized decomposition, and reduced basis methodologies.} preprint = {http://preprints.sissa.it/xmlui/bitstream/handle/1963/35194/ECM_MOR.pdf?sequence=1\&isAllowed=y

}, doi = {10.1002/9781119176817.ecm2110}, author = {Francisco Chinesta and Antonio Huerta and Gianluigi Rozza and Karen Willcox} } @article {20.500.11767_11953, title = {A natural framework for isogeometric fluid-structure interaction based on BEM-shell coupling}, journal = {COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING}, volume = {316}, year = {2017}, pages = {522{\textendash}546}, doi = {10.1016/j.cma.2016.08.008}, url = {http://cdsads.u-strasbg.fr/abs/2017CMAME.316..522H}, author = {Luca Heltai and Kiendl, J. and Antonio DeSimone and Alessandro Reali} } @article {2017, title = {POD-Galerkin reduced order methods for CFD using Finite Volume Discretisation: vortex shedding around a circular cylinder}, journal = {Communications in Applied and Industrial Mathematics}, volume = {8}, year = {2017}, pages = {210-236}, abstract = {

Vortex shedding around circular cylinders is a well known and studied phenomenon that appears in many engineering fields. In this work a Reduced Order Model (ROM) of the incompressible flow around a circular cylinder, built performing a Galerkin projection of the governing equations onto a lower dimensional space is presented. The reduced basis space is generated using a Proper Orthogonal Decomposition (POD) approach. In particular the focus is into (i) the correct reproduction of the pressure field, that in case of the vortex shedding phenomenon, is of primary importance for the calculation of the drag and lift coefficients; (ii) for this purpose the projection of the Governing equations (momentum equation and Poisson equation for pressure) is performed onto different reduced basis space for velocity and pressure, respectively; (iii) all the relevant modifications necessary to adapt standard finite element POD-Galerkin methods to a finite volume framework are presented. The accuracy of the reduced order model is assessed against full order results.

}, doi = {10.1515/caim-2017-0011}, author = {Giovanni Stabile and Saddam Hijazi and Andrea Mola and Stefano Lorenzi and Gianluigi Rozza} } @article {2017, title = {Time quasi-periodic gravity water waves in finite depth}, number = {arXiv;1708.01517}, year = {2017}, abstract = {We prove the existence and the linear stability of Cantor families of small amplitude time quasi-periodic standing water wave solutions - namely periodic and even in the space variable x - of a bi-dimensional ocean with finite depth under the action of pure gravity. Such a result holds for all the values of the depth parameter in a Borel set of asymptotically full measure. This is a small divisor problem. The main difficulties are the quasi-linear nature of the gravity water waves equations and the fact that the linear frequencies grow just in a sublinear way at infinity. We overcome these problems by first reducing the linearized operators obtained at each approximate quasi-periodic solution along the Nash-Moser iteration to constant coefficients up to smoothing operators, using pseudo-differential changes of variables that are quasi-periodic in time. Then we apply a KAM reducibility scheme which requires very weak Melnikov non-resonance conditions (losing derivatives both in time and space), which we are able to verify for most values of the depth parameter using degenerate KAM theory arguments.}, url = {http://preprints.sissa.it/handle/1963/35296}, author = {P Baldi and Massimiliano Berti and Emanuele Haus and Riccardo Montalto} } @article {MolaHeltaiDeSimone2017, title = {Wet and Dry Transom Stern Treatment for Unsteady and Nonlinear Potential Flow Model for Naval Hydrodynamics Simulations}, journal = {Journal of Ship Research}, volume = {61}, number = {1}, year = {2017}, pages = {1{\textendash}14}, abstract = {

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.

}, doi = {https://doi.org/10.5957/JOSR.61.1.160016}, author = {Andrea Mola and Luca Heltai and Antonio DeSimone} } @article {20.500.11767_14497, title = {The deal.II Library, Version 8.3}, journal = {ARCHIVE OF NUMERICAL SOFTWARE}, volume = {4}, year = {2016}, pages = {1{\textendash}11}, doi = {10.11588/ans.2016.100.23122}, url = {http://nbn-resolving.de/urn:nbn:de:bsz:16-ans-231226}, author = {Wolfgang Bangerth and Timo Heister and Luca Heltai and G. Kanschat and Martin Kronbichler and Matthias Maier and Bruno Turcksin} } @article {20.500.11767_14498, title = {The deal.II library, Version 8.4}, journal = {JOURNAL OF NUMERICAL MATHEMATICS}, volume = {24}, year = {2016}, pages = {135{\textendash}141}, doi = {10.1515/jnma-2016-1045}, url = {https://www.math.clemson.edu/ heister/preprints/deal84-preprint.pdf}, author = {Wolfgang Bangerth and Denis Davydov and Timo Heister and Luca Heltai and G. Kanschat and Martin Kronbichler and Matthias Maier and Bruno Turcksin and David Wells} } @article {20.500.11767_11949, title = {Error Estimates of B-spline based finite-element method for the wind-driven ocean circulation}, journal = {JOURNAL OF SCIENTIFIC COMPUTING}, volume = {69}, year = {2016}, pages = {430{\textendash}459}, doi = {10.1007/s10915-016-0201-1}, author = {Rotundo, N. and Kim, T. -Y. and Jiang, W. and Luca Heltai and Fried, E.} } @article {MR3507270, title = {$hp$-version discontinuous Galerkin methods for advection-diffusion-reaction problems on polytopic meshes}, journal = {ESAIM Math. Model. Numer. Anal.}, volume = {50}, number = {3}, year = {2016}, pages = {699{\textendash}725}, issn = {0764-583X}, doi = {10.1051/m2an/2015059}, url = {https://doi.org/10.1051/m2an/2015059}, author = {Andrea Cangiani and Dong, Zhaonan and E.H. Georgoulis and Houston, Paul} } @article {2016, title = {Isogeometric analysis-based reduced order modelling for incompressible linear viscous flows in parametrized shapes}, number = {AMOS Advanced Modelling and Simulation in Engineering Sciences}, year = {2016}, institution = {Springer, AMOS Advanced Modelling and Simulation in Engineering Sciences}, abstract = {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.}, url = {http://urania.sissa.it/xmlui/handle/1963/35199}, author = {Filippo Salmoiraghi and F. Ballarin and Luca Heltai and Gianluigi Rozza} } @article {20.500.11767_11950, title = {LinearOperator {\textendash} a generic, high-level expression syntax for linear algebra}, journal = {COMPUTERS \& MATHEMATICS WITH APPLICATIONS}, volume = {72}, year = {2016}, pages = {1{\textendash}24}, doi = {10.1016/j.camwa.2016.04.024}, author = {Matthias Maier and Mauro Bardelloni and Luca Heltai} } @inbook {2016, title = {Model Order Reduction: a survey}, booktitle = {Wiley Encyclopedia of Computational Mechanics, 2016}, year = {2016}, publisher = {Wiley}, organization = {Wiley}, url = {http://urania.sissa.it/xmlui/handle/1963/35194}, author = {Francisco Chinesta and Antonio Huerta and Gianluigi Rozza and Karen Willcox} } @inbook {MR3585793, title = {Review of discontinuous Galerkin finite element methods for partial differential equations on complicated domains}, booktitle = {Building bridges: connections and challenges in modern approaches to numerical partial differential equations}, series = {Lect. Notes Comput. Sci. Eng.}, volume = {114}, year = {2016}, pages = {279{\textendash}308}, publisher = {Springer, [Cham]}, organization = {Springer, [Cham]}, author = {Antonietti, Paola F. and Andrea Cangiani and Collis, Joe and Dong, Zhaonan and E.H. Georgoulis and Stefano Giani and Houston, Paul} } @conference {mola2016ship, title = {Ship Sinkage and Trim Predictions Based on a CAD Interfaced Fully Nonlinear Potential Model}, booktitle = {The 26th International Ocean and Polar Engineering Conference}, volume = {3}, year = {2016}, pages = {511{\textendash}518}, publisher = {International Society of Offshore and Polar Engineers}, organization = {International Society of Offshore and Polar Engineers}, author = {Andrea Mola and Luca Heltai and Antonio DeSimone and Massimiliano Berti} } @article {PhysRevD.93.066001, title = {Towards a gauge theory interpretation of the real topological string}, journal = {Phys. Rev. D}, volume = {93}, year = {2016}, month = {Mar}, pages = {066001}, publisher = {American Physical Society}, abstract = {

We consider the real topological string on certain noncompact toric Calabi-Yau three-folds $\mathbb{X}$, in its physical realization describing an orientifold of type IIA on $\mathbb{X}$\ with an O4-plane and a single D4-brane stuck on top. The orientifold can be regarded as a new kind of surface operator on the gauge theory with 8 supercharges arising from the singular geometry. We use the M-theory lift of this system to compute the real Gopakumar-Vafa invariants (describing wrapped M2-brane Bogomol{\textquoteright}nyi-Prasad-Sommerfield (BPS) states) for diverse geometries. We show that the real topological string amplitudes pick up certain signs across flop transitions, in a well-defined pattern consistent with continuity of the real BPS invariants. We further give some preliminary proposals of an intrinsically gauge theoretical description of the effect of the surface operator in the gauge theory partition function.

}, doi = {10.1103/PhysRevD.93.066001}, url = {https://link.aps.org/doi/10.1103/PhysRevD.93.066001}, author = {Hayashi, Hirotaka and Nicol{\`o} Piazzalunga and Uranga, Angel M.} } @article {2015, title = {Benchmarking the Immersed Finite Element Method for Fluid-Structure Interaction Problems}, journal = {Computers and Mathematics with Applications 69 (2015) 1167{\textendash}1188}, number = {;69}, year = {2015}, abstract = {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.}, doi = {10.1016/j.camwa.2015.03.012}, author = {Roy Saswati and Luca Heltai and Francesco Costanzo} } @book {2015, title = {Certified Reduced Basis Methods for Parametrized Partial Differential Equations}, series = {Springer Briefs in Mathematics}, year = {2015}, pages = {135}, publisher = {Springer}, organization = {Springer}, edition = {1}, address = {Switzerland}, abstract = {

This book provides a thorough introduction to the mathematical and algorithmic aspects of certified reduced basis methods for parametrized partial differential equations. Central aspects ranging from model construction, error estimation and computational efficiency to empirical interpolation methods are discussed in detail for coercive problems. More advanced aspects associated with time-dependent problems, non-compliant and non-coercive problems and applications with geometric variation are also discussed as examples.

}, keywords = {a posteriori error bounds, empirical interpolation, parametrized partial differential equations, reduced basis methods, greedy algorithms}, isbn = {978-3-319-22469-5}, issn = {2191-8201}, author = {Jan S Hesthaven and Gianluigi Rozza and Benjamin Stamm} } @article {2015, title = {Deal2lkit: a Toolkit Library for High Performance Programming in deal.II}, year = {2015}, publisher = {SISSA}, abstract = {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.}, url = {http://urania.sissa.it/xmlui/handle/1963/35006}, author = {Alberto Sartori and Nicola Giuliani and Mauro Bardelloni and Luca Heltai} } @article {2015, title = {The deal.II Library, Version 8.2}, journal = {Archive of Numerical Software, vol. 3, n. 100, (2015), pages : 1-8}, year = {2015}, abstract = {This paper provides an overview of the new features of the finite element library deal.II version 8.2}, doi = {10.11588/ans.2015.100.18031}, url = {http://urania.sissa.it/xmlui/handle/1963/34464}, author = {Wolfgang Bangerth and Timo Heister and Luca Heltai and G. Kanschat and Martin Kronbichler and Matthias Maier and Bruno Turcksin and T. D. Young} } @article {2015, title = {FEM SUPG stabilisation of mixed isoparametric BEMs: application to linearised free surface flows}, journal = {Engineering Analysis with Boundary Elements 59 (2015), pp. 8-22}, year = {2015}, abstract = {

In finite element formulations, transport dominated problems are often stabilised through the Streamline-Upwind-Petrov{\textendash}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.

}, doi = {10.1016/j.enganabound.2015.04.006}, url = {http://urania.sissa.it/xmlui/handle/1963/34466}, author = {Nicola Giuliani and Andrea Mola and Luca Heltai and L. Formaggia} } @article {2015, title = {Reduced basis approximation and a-posteriori error estimation for the coupled Stokes-Darcy system}, journal = {Advances in Computational Mathematics}, volume = {special issue for MoRePaS 2012}, year = {2015}, abstract = {

The coupling of a free flow with a flow through porous media has many potential applications in several fields related with computational science and engineering, such as blood flows, environmental problems or food technologies. We present a reduced basis method for such coupled problems. The reduced basis method is a model order reduction method applied in the context of parametrized systems. Our approach is based on a heterogeneous domain decomposition formulation, namely the Stokes-Darcy problem. Thanks to an offline/online-decomposition, computational times can be drastically reduced. At the same time the induced error can be bounded by fast evaluable a-posteriori error bounds. In the offline-phase the proposed algorithms make use of the decomposed problem structure. Rigorous a-posteriori error bounds are developed, indicating the accuracy of certain lifting operators used in the offline-phase as well as the accuracy of the reduced coupled system. Also, a strategy separately bounding pressure and velocity errors is extended. Numerical experiments dealing with groundwater flow scenarios demonstrate the efficiency of the approach as well as the limitations regarding a-posteriori error estimation.

}, keywords = {Domain decomposition, Error estimation, Non-coercive problem, Porous medium equation, Reduced basis method, Stokes flow}, issn = {1019-7168}, doi = { 10.1007/s10444-014-9396-6}, author = {Immanuel Martini and Gianluigi Rozza and Bernard Haasdonk} } @article {2015, title = {Reduced Basis Isogeometric Methods (RB-IGA) for the real-time simulation of potential flows about parametrized NACA airfoils}, journal = {Comput Methods Appl Mech Eng. 2015;284:1147{\textendash}1180}, number = {;284}, year = {2015}, abstract = {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.}, doi = {10.1016/j.cma.2014.11.037}, author = {Andrea Manzoni and Filippo Salmoiraghi and Luca Heltai} } @article {2014, title = {Finite dimensional Kadomtsev-Petviashvili τ-functions. I. Finite Grassmannians}, number = {Journal of Mathematical Physics;volume 55; issue 8; article number 1.4890818;}, year = {2014}, publisher = {American Institute of Physics Inc.}, abstract = {We study τ-functions of the Kadomtsev-Petviashvili hierarchy in terms of abelian group actions on finite dimensional Grassmannians, viewed as subquotients of the Hilbert space Grassmannians of Sato, Segal, and Wilson. A determinantal formula of Gekhtman and Kasman involving exponentials of finite dimensional matrices is shown to follow naturally from such reductions. All reduced flows of exponential type generated by matrices with arbitrary nondegenerate Jordan forms are derived, both in the Grassmannian setting and within the fermionic operator formalism. A slightly more general determinantal formula involving resolvents of the matrices generating the flow, valid on the big cell of the Grassmannian, is also derived. An explicit expression is deduced for the Pl{\"u}cker coordinates appearing as coefficients in the Schur function expansion of the τ-function.}, doi = {10.1063/1.4890818}, url = {http://urania.sissa.it/xmlui/handle/1963/34952}, author = {Ferenc Balogh and Tiago Fonseca and John P. Harnad} } @proceedings {2014, title = {A fully nonlinear potential model for ship hydrodynamics directly interfaced with CAD data structures}, year = {2014}, publisher = {SISSA}, abstract = {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.}, keywords = {ship hydrodynamics}, author = {Andrea Mola and Luca Heltai and Antonio DeSimone} } @article {MR3211116, title = {$hp$-version discontinuous Galerkin methods on polygonal and polyhedral meshes}, journal = {Math. Models Methods Appl. Sci.}, volume = {24}, number = {10}, year = {2014}, pages = {2009{\textendash}2041}, issn = {0218-2025}, doi = {10.1142/S0218202514500146}, url = {https://doi.org/10.1142/S0218202514500146}, author = {Andrea Cangiani and E.H. Georgoulis and Houston, Paul} } @article {2012, title = {Nonsingular Isogeometric Boundary Element Method for Stokes Flows in 3D}, number = {Computer Methods in Applied Mechanics and Engineering}, year = {2014}, publisher = {Elsevier}, abstract = {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.}, keywords = {Isogeometric Analysis}, doi = {10.1016/j.cma.2013.09.017}, url = {http://hdl.handle.net/1963/6326}, author = {Luca Heltai and Marino Arroyo and Antonio DeSimone} } @conference {mola2014ship, title = {Potential Model for Ship Hydrodynamics Simulations Directly Interfaced with CAD Data Structures}, booktitle = {The 24th International Ocean and Polar Engineering Conference}, volume = {4}, year = {2014}, pages = {815{\textendash}822}, publisher = {International Society of Offshore and Polar Engineers}, organization = {International Society of Offshore and Polar Engineers}, author = {Andrea Mola and Luca Heltai and Antonio DeSimone and Massimiliano Berti} } @article {2013, title = {The deal.II Library, Version 8.1}, number = {arXiv:1312.2266;}, year = {2013}, note = {5 pages}, institution = {SISSA}, abstract = {This paper provides an overview of the new features of the finite element library deal.II version 8.0.}, url = {http://hdl.handle.net/1963/7236}, author = {Wolfgang Bangerth and Timo Heister and Luca Heltai and G. Kanschat and Martin Kronbichler and Matthias Maier and Bruno Turcksin and T. D. Young} } @article {2013, title = {Reduced basis approximation and a posteriori error estimation for Stokes flows in parametrized geometries: roles of the inf-sup stability constants}, journal = {Numerische Mathematik, 2013}, number = {SISSA Preprint;43/2012/M}, year = {2013}, publisher = {Springer}, abstract = {In this paper we review and we extend the reduced basis approximation and a posteriori error estimation for steady Stokes flows in a ffinely parametrized geometries, focusing on the role played by the Brezzi\\\'s and Babu ska\\\'s stability constants. The crucial ingredients of the methodology are a Galerkin projection onto a low-dimensional space of basis functions properly selected, an a ne parametric dependence enabling to perform competitive Off ine-Online splitting in the computational\\r\\nprocedure and a rigorous a posteriori error estimation on eld variables.\\r\\nThe combination of these three factors yields substantial computational savings which are at the basis of an e fficient model order reduction, ideally suited for real-time simulation and many-query contexts (e.g. optimization, control or parameter identi cation). In particular, in this work we focus on i) the stability of the reduced basis approximation based on the Brezzi\\\'s saddle point theory and the introduction of a supremizer operator on the pressure terms, ii) a rigorous a posteriori error estimation procedure for velocity and pressure elds based on the Babu ska\\\'s inf-sup constant (including residuals calculations), iii) the computation of a lower bound of the stability constant, and iv) di erent options for the reduced basis spaces construction. We present some illustrative results for both\\r\\ninterior and external steady Stokes flows in parametrized geometries representing two parametrized classical Poiseuille and Couette \\r\\nflows, a channel contraction and a simple flow control problem around a curved obstacle.}, keywords = {parametrized Stokes equations}, url = {http://hdl.handle.net/1963/6339}, author = {Gianluigi Rozza and Phuong Huynh and Andrea Manzoni} } @article {2012, title = {A stable and adaptive semi-Lagrangian potential model for unsteady and nonlinear ship-wave interactions}, journal = {Engineering Analysis with Boundary Elements, 37(1):128 {\textendash} 143, 2013.}, number = {SISSA;06/2012/M}, year = {2013}, publisher = {SISSA}, abstract = {

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.

}, keywords = {Unsteady ship-wave interaction}, url = {http://hdl.handle.net/1963/5669}, author = {Andrea Mola and Luca Heltai and Antonio DeSimone} } @article {2012, title = {A Codazzi-like equation and the singular set for C1 smooth surfaces in the Heisenberg group.}, journal = {Journal fur die Reine und Angewandte Mathematik, Issue 671, October 2012, Pages 131-198}, year = {2012}, publisher = {SISSA}, abstract = {In this paper, we study the structure of the singular set for a C 1 smooth surface in the 3-dimensional Heisenberg group ℍ 1. We discover a Codazzi-like equation for the p-area element along the characteristic curves on the surface. Information obtained from this ordinary differential equation helps us to analyze the local configuration of the singular set and the characteristic curves. In particular, we can estimate the size and obtain the regularity of the singular set. We understand the global structure of the singular set through a Hopf-type index theorem. We also justify the Codazzi-like equation by proving a fundamental theorem for local surfaces in ℍ 1}, doi = {10.1515/CRELLE.2011.159}, url = {http://hdl.handle.net/1963/6556}, author = {Andrea Malchiodi and Paul Yang and Jih-Hsin Cheng and JennFang Hwang} } @inbook {9897, title = {Computing optimal strokes for low reynolds number swimmers}, booktitle = {Natural locomotion in fluids and on surfaces : swimming, flying, and sliding / editors Stephen Childress, Anette Hosoi, William W. Schultz, and Z. Jane Wang, editors,}, number = {The IMA volumes in mathematics and its applications;155}, year = {2012}, publisher = {Springer}, organization = {Springer}, abstract = {

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.

}, keywords = {Numerical analysis.}, isbn = {9781461439967}, doi = {10.1007/978-1-4614-3997-4_13}, url = {http://hdl.handle.net/1963/6445}, author = {Antonio DeSimone and Luca Heltai and Fran{\c c}ois Alouges and Lefebvre-Lepot Aline} } @article {2012, title = {A Fully Coupled Immersed Finite Element Method for Fluid Structure Interaction via the Deal.II Library}, number = {arXiv:1209.2811;}, year = {2012}, note = {28 pages, 9 figures}, institution = {SISSA}, abstract = {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.}, keywords = {Finite Element Method, Immersed Boundary Method, Immersed Finite Element Method}, url = {http://hdl.handle.net/1963/6255}, author = {Luca Heltai and Saswati Roy and Francesco Costanzo} } @article {2012, title = {Reverse engineering the euglenoid movement}, journal = {Proceedings of the National Academy of Sciences of the United States of America. Volume 109, Issue 44, 30 October 2012, Pages 17874-17879}, year = {2012}, abstract = {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.}, keywords = {microswimmers}, doi = {10.1073/pnas.1213977109}, url = {http://hdl.handle.net/1963/6444}, author = {Marino Arroyo and Luca Heltai and Daniel Mill{\'a}n and Antonio DeSimone} } @conference {Mola2012, title = {A stable semi-lagrangian potential method for the simulation of ship interaction with unsteady and nonlinear waves}, booktitle = {17th Int. Conf. Ships Shipp. Res.}, year = {2012}, author = {Andrea Mola and Luca Heltai and Antonio DeSimone} } @article {2012, title = {Variational implementation of immersed finite element methods}, journal = {Computer Methods in Applied Mechanics and Engineering. Volume 229-232, 1 July 2012, Pages 110-127}, number = {arXiv:1110.2063;}, year = {2012}, note = {42 pages, 5 figures, Revision 1}, publisher = {Elsevier}, abstract = {

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.

}, keywords = {Turbulent flow}, doi = {10.1016/j.cma.2012.04.001}, url = {http://hdl.handle.net/1963/6462}, author = {Luca Heltai and Francesco Costanzo} } @article {Mola2011, title = {Multi-physics modelling and sensitivity analysis of olympic rowing boat dynamics}, journal = {Sports Engineering}, volume = {14}, number = {2-4}, year = {2011}, month = {nov}, pages = {85{\textendash}94}, publisher = {Springer Nature}, doi = {10.1007/s12283-011-0075-2}, url = {https://doi.org/10.1007/s12283-011-0075-2}, author = {Andrea Mola and Mehdi Ghommem and Muhammad R. Hajj} } @article {2011, title = {Numerical Strategies for Stroke Optimization of Axisymmetric Microswimmers}, journal = {Mathematical Models and Methods in Applied Sciences 21 (2011) 361-387}, number = {SISSA;33/2009/M}, year = {2011}, publisher = {World Scientific}, abstract = {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.}, keywords = {Optimal swimming}, doi = {10.1142/S0218202511005088}, url = {http://hdl.handle.net/1963/3657}, author = {Fran{\c c}ois Alouges and Antonio DeSimone and Luca Heltai} } @article {2010, title = {Optimally swimming Stokesian Robots}, number = {SISSA;54/2010/M}, year = {2010}, abstract = {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.}, url = {http://hdl.handle.net/1963/3929}, author = {Fran{\c c}ois Alouges and Antonio DeSimone and Luca Heltai and Aline Lefebvre and Benoit Merlet} } @conference {10.1007/978-90-481-9195-6_4, title = {A Phase Field Approach to Wetting and Contact Angle Hysteresis Phenomena}, booktitle = {IUTAM Symposium on Variational Concepts with Applications to the Mechanics of Materials}, year = {2010}, pages = {51{\textendash}63}, publisher = {Springer Netherlands}, organization = {Springer Netherlands}, address = {Dordrecht}, abstract = {

We discuss a phase field model for the numerical simulation of contact angle hysteresis phenomena in wetting. The performance of the model is assessed by comparing its predictions with experimental data on the critical size of drops that can stick on a vertical glass plate.

}, isbn = {978-90-481-9195-6}, author = {Antonio DeSimone and Livio Fedeli and Turco, Alessandro}, editor = {Hackl, Klaus} } @article {2010, title = {The role of membrane viscosity in the dynamics of fluid membranes}, number = {SISSA;55/2010/M}, year = {2010}, abstract = {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.}, url = {http://hdl.handle.net/1963/3930}, author = {Marino Arroyo and Antonio DeSimone and Luca Heltai} } @article {Mola2009, title = {Low-Frequency Variations of Force Coefficients on Square Cylinders with Sharp and Rounded Corners}, journal = {Journal of Structural Engineering}, volume = {135}, number = {7}, year = {2009}, month = {jul}, pages = {828{\textendash}835}, publisher = {American Society of Civil Engineers ({ASCE})}, doi = {10.1061/(asce)st.1943-541x.0000034}, url = {https://doi.org/10.1061/(asce)st.1943-541x.0000034}, author = {Andrea Mola and Giancarlo Bordonaro and Muhammad R. Hajj} } @article {2009, title = {Stratos: a code for 3D free surface flows with floating constraints}, number = {SISSA;41/2009/M}, year = {2009}, abstract = {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......}, url = {http://hdl.handle.net/1963/3701}, author = {Antonio DeSimone and B. Bianchi and Luca Heltai} } @article {2009, title = {Tools for the Solution of PDEs Defined on Curved Manifolds with deal.II}, number = {SISSA;42/2009/M}, year = {2009}, abstract = {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.}, url = {http://hdl.handle.net/1963/3700}, author = {Antonio DeSimone and Luca Heltai and Cataldo Manigrasso} } @article {2007, title = {Asymptotic behaviour of smooth solutions for partially dissipative hyperbolic systems with a convex entropy}, journal = {Comm. Pure Appl. Math. 60 (2007) 1559-1622}, number = {SISSA;83/2005/M}, year = {2007}, abstract = {We study the asymptotic time behavior of global smooth solutions to general entropy dissipative hyperbolic systems of balance law in m space dimensions, under the Shizuta-Kawashima condition.}, doi = {10.1002/cpa.20195}, url = {http://hdl.handle.net/1963/1780}, author = {Stefano Bianchini and Bernard Hanouzet and Roberto Natalini} } @article {2006, title = {Reflection symmetries for multiqubit density operators}, journal = {J. Math. Phys. 47 (2006) 032104}, number = {arXiv.org;quant-ph/0405123v2}, year = {2006}, abstract = {For multiqubit density operators in a suitable tensorial basis, we show that a number of nonunitary operations used in the detection and synthesis of entanglement are classifiable as reflection symmetries, i.e., orientation changing rotations. While one-qubit reflections correspond to antiunitary symmetries, as is known for example from the partial transposition criterion, reflections on the joint density of two or more qubits are not accounted for by the Wigner Theorem and are well-posed only for sufficiently mixed states. One example of such nonlocal reflections is the unconditional NOT operation on a multiparty density, i.e., an operation yelding another density and such that the sum of the two is the identity operator. This nonphysical operation is admissible only for sufficiently mixed states.}, doi = {10.1063/1.2181827}, url = {http://hdl.handle.net/1963/2121}, author = {Claudio Altafini and Timothy F. Havel} } @article {Bertola:Semiiso, title = {Semiclassical orthogonal polynomials, matrix models and isomonodromic tau functions}, journal = {Comm. Math. Phys.}, volume = {263}, number = {2}, year = {2006}, pages = {401{\textendash}437}, issn = {0010-3616}, author = {Marco Bertola and B. Eynard and Harnad, J.} } @article {2006, title = {Semistability vs. nefness for (Higgs) vector bundles}, journal = {Differential Geom. Appl. 24 (2006) 403-416}, number = {arXiv.org;math/0310040v3}, year = {2006}, abstract = {According to Miyaoka, a vector bundle E on a smooth projective curve is semistable if and only if a certain numerical class in the projectivized bundle PE is nef. We establish a similar criterion for the semistability of Higgs bundles: namely, such a bundle is semistable if and only if for every integer s between 0 and the rank of E, a suitable numerical class in the scheme parametrizing the rank s locally-free Higgs quotients of E is nef. We also extend this result to higher-dimensional complex projective varieties by showing that the nefness of the above mentioned classes is equivalent to the semistability of the Higgs bundle E together with the vanishing of the discriminant of E.}, doi = {10.1016/j.difgeo.2005.12.007}, url = {http://hdl.handle.net/1963/2237}, author = {Ugo Bruzzo and Daniel Hernandez Ruiperez} } @article {2005, title = {Minimal surfaces in pseudohermitian geometry}, journal = {Annali della Scuola Normale Superiore di Pisa, Classe di Scienze (5), 4 (2005) 129-177.}, number = {arXiv:math/0401136;}, year = {2005}, publisher = {Scuola Normale Superiore}, abstract = {We consider surfaces immersed in three-dimensional pseudohermitian manifolds. We define the notion of (p-)mean curvature and of the associated (p-)minimal surfaces, extending some concepts previously given for the (flat) Heisenberg group. We interpret the p-mean curvature not only as the tangential sublaplacian of a defining function, but also as the curvature of a characteristic curve, and as a quantity in terms of calibration geometry. As a differential equation, the p-minimal surface equation is degenerate (hyperbolic and elliptic). To analyze the singular set, we formulate some {\em extension} theorems, which describe how the characteristic curves meet the singular set. This allows us to classify the entire solutions to this equation and to solve a Bernstein-type problem (for graphs over the $xy$-plane) in the Heisenberg group $H_1$. In $H_{1}$, identified with the Euclidean space $R^{3}$, the p-minimal surfaces are classical ruled surfaces with the rulings generated by Legendrian lines. We also prove a uniqueness theorem for the Dirichlet problem under a condition on the size of the singular set in two dimensions, and generalize to higher dimensions without any size control condition. We also show that there are no closed, connected, $C^{2}$ smoothly immersed constant p-mean curvature or p-minimal surfaces of genus greater than one in the standard $S^{3}.$ This fact continues to hold when $S^{3}$ is replaced by a general spherical pseudohermitian 3-manifold.}, doi = {10.2422/2036-2145.2005.1.05}, url = {http://hdl.handle.net/1963/4579}, author = {Jih-Hsin Cheng and JennFang Hwang and Andrea Malchiodi and Paul Yang} } @article {2005, title = {Stability of solutions of quasilinear parabolic equations}, journal = {J. Math. Anal. Appl. 308 (2005) 221-239}, number = {SISSA;51/2003/M}, year = {2005}, publisher = {Elsevier}, abstract = {We bound the difference between solutions $u$ and $v$ of $u_t = a\\\\Delta u+\\\\Div_x f+h$ and $v_t = b\\\\Delta v+\\\\Div_x g+k$ with initial data $\\\\phi$ and $ \\\\psi$, respectively, by $\\\\Vert u(t,\\\\cdot)-v(t,\\\\cdot)\\\\Vert_{L^p(E)}\\\\le A_E(t)\\\\Vert \\\\phi-\\\\psi\\\\Vert_{L^\\\\infty(\\\\R^n)}^{2\\\\rho_p}+ B(t)(\\\\Vert a-b\\\\Vert_{\\\\infty}+ \\\\Vert \\\\nabla_x\\\\cdot f-\\\\nabla_x\\\\cdot g\\\\Vert_{\\\\infty}+ \\\\Vert f_u-g_u\\\\Vert_{\\\\infty} + \\\\Vert h-k\\\\Vert_{\\\\infty})^{\\\\rho_p} \\\\abs{E}^{\\\\eta_p}$. Here all functions $a$, $f$, and $h$ are smooth and bounded, and may depend on $u$, $x\\\\in\\\\R^n$, and $t$. The functions $a$ and $h$ may in addition depend on $\\\\nabla u$. Identical assumptions hold for the functions that determine the solutions $v$. Furthermore, $E\\\\subset\\\\R^n$ is assumed to be a bounded set, and $\\\\rho_p$ and $\\\\eta_p$ are fractions that depend on $n$ and $p$. The diffusion coefficients $a$ and $b$ are assumed to be strictly positive and the initial data are smooth.}, doi = {10.1016/j.jmaa.2005.01.026}, url = {http://hdl.handle.net/1963/2892}, author = {Giuseppe Maria Coclite and Helge Holden} } @article {2004, title = {Fredholm modules for quantum euclidean spheres}, journal = {J. Geom. Phys. 49 (2004) 272-293}, number = {SISSA;66/2002/FM}, year = {2004}, publisher = {SISSA Library}, abstract = {The quantum Euclidean spheres, $S_q^{N-1}$, are (noncommutative) homogeneous spaces of quantum orthogonal groups, $\\\\SO_q(N)$. The *-algebra $A(S^{N-1}_q)$ of polynomial functions on each of these is given by generators and relations which can be expressed in terms of a self-adjoint, unipotent matrix. We explicitly construct complete sets of generators for the K-theory (by nontrivial self-adjoint idempotents and unitaries) and the K-homology (by nontrivial Fredholm modules) of the spheres $S_q^{N-1}$. We also construct the corresponding Chern characters in cyclic homology and cohomology and compute the pairing of K-theory with K-homology. On odd spheres (i. e., for N even) we exhibit unbounded Fredholm modules by means of a natural unbounded operator D which, while failing to have compact resolvent, has bounded commutators with all elements in the algebra $A(S^{N-1}_q)$.}, doi = {10.1016/S0393-0440(03)00092-5}, url = {http://hdl.handle.net/1963/1636}, author = {Eli Hawkins and Giovanni Landi} } @article {Bertola:DifferentialCMP, title = {Differential systems for biorthogonal polynomials appearing in 2-matrix models and the associated Riemann-Hilbert problem}, journal = {Comm. Math. Phys.}, volume = {243}, number = {2}, year = {2003}, pages = {193{\textendash}240}, issn = {0010-3616}, author = {Marco Bertola and B. Eynard and Harnad, J.} } @article {Bertola:PartitionJPA, title = {Partition functions for matrix models and isomonodromic tau functions}, journal = {J. Phys. A}, volume = {36}, number = {12}, year = {2003}, note = {Random matrix theory}, pages = {3067{\textendash}3083}, issn = {0305-4470}, author = {Marco Bertola and B. Eynard and Harnad, J.} } @article {Bertola:DualityCMP, title = {Duality, biorthogonal polynomials and multi-matrix models}, journal = {Comm. Math. Phys.}, volume = {229}, number = {1}, year = {2002}, pages = {73{\textendash}120}, issn = {0010-3616}, author = {Marco Bertola and B. Eynard and Harnad, J.} } @article {2002, title = {Relatively stable bundles over elliptic fibrations}, journal = {Math. Nachr. 238 (2002) 23-36}, number = {arXiv.org;math/0109123v2}, year = {2002}, publisher = {Wiley}, abstract = {We consider a relative Fourier-Mukai transform defined on elliptic fibrations over an arbitrary normal base scheme. This is used to construct relative Atiyah sheaves and generalize Atiyah\\\'s and Tu\\\'s results about semistable sheaves over elliptic curves to the case of elliptic fibrations. Moreover we show that this transform preserves relative (semi)stability of sheaves of positive relative degree.}, url = {http://hdl.handle.net/1963/3132}, author = {Claudio Bartocci and Ugo Bruzzo and Daniel Hernandez Ruiperez and Jose M. Munoz Porras} }