In the reduced order modeling (ROM) framework, the solution of a parametric partial differential equation is approximated by combining the high-fidelity solutions of the problem at hand for several properly chosen configurations. Examples of the ROM application, in the naval field, can be found in [31, 24]. Mandatory ingredient for the ROM methods is the relation between the high-fidelity solutions and the parameters. Dealing with geometrical parameters, especially in the industrial context, this relation may be unknown and not trivial (simulations over hand morphed geometries) or very complex (high number of parameters or many nested morphing techniques). To overcome these scenarios, we propose in this contribution an efficient and complete data-driven framework involving ROM techniques for shape design and optimization, extending the pipeline presented in [7]. By applying the singular value decomposition (SVD) to the points coordinates defining the hull geometry –- assuming the topology is inaltered by the deformation –-, we are able to compute the optimal space which the deformed geometries belong to, hence using the modal coefficients as the new parameters we can reconstruct the parametric formulation of the domain. Finally the output of interest is approximated using the proper orthogonal decomposition with interpolation technique. To conclude, we apply this framework to a naval shape design problem where the bulbous bow is morphed to reduce the total resistance of the ship advancing in calm water.

JF - VIII International Conference on Computational Methods in Marine Engineering ER - TY - CONF T1 - Efficient Reduction in Shape Parameter Space Dimension for Ship Propeller Blade Design T2 - VIII International Conference on Computational Methods in Marine Engineering Y1 - 2019 A1 - Andrea Mola A1 - Marco Tezzele A1 - Mahmoud Gadalla A1 - Valdenazzi, Federica A1 - Grassi, Davide A1 - Padovan, Roberta A1 - Gianluigi Rozza AB -In this work, we present the results of a ship propeller design optimization campaign carried out in the framework of the research project PRELICA, funded by the Friuli Venezia Giulia regional government. The main idea of this work is to operate on a multidisciplinary level to identify propeller shapes that lead to reduced tip vortex-induced pressure and increased efficiency without altering the thrust. First, a specific tool for the bottom-up construction of parameterized propeller blade geometries has been developed. The algorithm proposed operates with a user defined number of arbitrary shaped or NACA airfoil sections, and employs arbitrary degree NURBS to represent the chord, pitch, skew and rake distribution as a function of the blade radial coordinate. The control points of such curves have been modified to generate, in a fully automated way, a family of blade geometries depending on as many as 20 shape parameters. Such geometries have then been used to carry out potential flow simulations with the Boundary Element Method based software PROCAL. Given the high number of parameters considered, such a preliminary stage allowed for a fast evaluation of the performance of several hundreds of shapes. In addition, the data obtained from the potential flow simulation allowed for the application of a parameter space reduction methodology based on active subspaces (AS) property, which suggested that the main propeller performance indices are, at a first but rather accurate approximation, only depending on a single parameter which is a linear combination of all the original geometric ones. AS analysis has also been used to carry out a constrained optimization exploiting response surface method in the reduced parameter space, and a sensitivity analysis based on such surrogate model. The few selected shapes were finally used to set up high fidelity RANS simulations and select an optimal shape.

JF - VIII International Conference on Computational Methods in Marine Engineering ER - TY - JOUR T1 - Dimension reduction in heterogeneous parametric spaces with application to naval engineering shape design problems JF - Advanced Modeling and Simulation in Engineering Sciences Y1 - 2018 A1 - Marco Tezzele A1 - Filippo Salmoiraghi A1 - Andrea Mola A1 - Gianluigi Rozza AB -We present the results of the first application in the naval architecture field of a methodology based on active subspaces properties for parameters space reduction. The physical problem considered is the one of the simulation of the hydrodynamic flow past the hull of a ship advancing in calm water. Such problem is extremely relevant at the preliminary stages of the ship design, when several flow simulations are typically carried out by the engineers to assess the dependence of the hull total resistance on the geometrical parameters of the hull, and others related with flows and hull properties. Given the high number of geometric and physical parameters which might affect the total ship drag, the main idea of this work is to employ the active subspaces properties to identify possible lower dimensional structures in the parameter space. Thus, a fully automated procedure has been implemented to produce several small shape perturbations of an original hull CAD geometry, in order to exploit the resulting shapes to run high fidelity flow simulations with different structural and physical parameters as well, and then collect data for the active subspaces analysis. The free form deformation procedure used to morph the hull shapes, the high fidelity solver based on potential flow theory with fully nonlinear free surface treatment, and the active subspaces analysis tool employed in this work have all been developed and integrated within SISSA mathLab as open source tools. The contribution will also discuss several details of the implementation of such tools, as well as the results of their application to the selected target engineering problem.

VL - 5 ER - TY - Generic T1 - An efficient shape parametrisation by free-form deformation enhanced by active subspace for hull hydrodynamic ship design problems in open source environment T2 - The 28th International Ocean and Polar Engineering Conference Y1 - 2018 A1 - Nicola Demo A1 - Marco Tezzele A1 - Andrea Mola A1 - Gianluigi Rozza KW - Active subspaces KW - Boundary element method KW - Dynamic mode decomposition KW - Fluid structure interaction KW - Free form deformation KW - Fully nonlinear potential KW - Numerical towing tank AB - In this contribution, we present the results of the application of a parameter space reduction methodology based on active subspaces to the hull hydrodynamic design problem. Several parametric deformations of an initial hull shape are considered to assess the influence of the shape parameters considered on the hull total drag. The hull resistance is typically computed by means of numerical simulations of the hydrodynamic flow past the ship. Given the high number of parameters involved - which might result in a high number of time consuming hydrodynamic simulations - assessing whether the parameters space can be reduced would lead to considerable computational cost reduction. Thus, the main idea of this work is to employ the active subspaces to identify possible lower dimensional structures in the parameter space, or to verify the parameter distribution in the position of the control points. To this end, a fully automated procedure has been implemented to produce several small shape perturbations of an original hull CAD geometry which are then used to carry out high-fidelity flow simulations and collect data for the active subspaces analysis. To achieve full automation of the open source pipeline described, both the free form deformation methodology employed for the hull perturbations and the solver based on unsteady potential flow theory, with fully nonlinear free surface treatment, are directly interfaced with CAD data structures and operate using IGES vendor-neutral file formats as input files. The computational cost of the fluid dynamic simulations is further reduced through the application of dynamic mode decomposition to reconstruct the steady state total drag value given only few initial snapshots of the simulation. The active subspaces analysis is here applied to the geometry of the DTMB-5415 naval combatant hull, which is which is a common benchmark in ship hydrodynamics simulations. JF - The 28th International Ocean and Polar Engineering Conference PB - International Society of Offshore and Polar Engineers CY - Sapporo, Japan UR - https://www.onepetro.org/conference-paper/ISOPE-I-18-481 ER - TY - CONF T1 - Model Order Reduction by means of Active Subspaces and Dynamic Mode Decomposition for Parametric Hull Shape Design Hydrodynamics T2 - Technology and Science for the Ships of the Future: Proceedings of NAV 2018: 19th International Conference on Ship & Maritime Research Y1 - 2018 A1 - Marco Tezzele A1 - Nicola Demo A1 - Mahmoud Gadalla A1 - Andrea Mola A1 - Gianluigi Rozza AB - We present the results of the application of a parameter space reduction methodology based on active subspaces (AS) to the hull hydrodynamic design problem. Several parametric deformations of an initial hull shape are considered to assess the influence of the shape parameters on the hull wave resistance. Such problem is relevant at the preliminary stages of the ship design, when several flow simulations are carried out by the engineers to establish a certain sensibility with respect to the parameters, which might result in a high number of time consuming hydrodynamic simulations. The main idea of this work is to employ the AS to identify possible lower dimensional structures in the parameter space. The complete pipeline involves the use of free form deformation to parametrize and deform the hull shape, the full order solver based on unsteady potential flow theory with fully nonlinear free surface treatment directly interfaced with CAD, the use of dynamic mode decomposition to reconstruct the final steady state given only few snapshots of the simulation, and the reduction of the parameter space by AS, and shared subspace. Response surface method is used to minimize the total drag. JF - Technology and Science for the Ships of the Future: Proceedings of NAV 2018: 19th International Conference on Ship & Maritime Research PB - IOS Press CY - Trieste, Italy UR - http://ebooks.iospress.nl/publication/49270 ER - TY - JOUR T1 - π-BEM : A flexible parallel implementation for adaptive , geometry aware , and high order boundary element methods JF - Advances in Engineering Software Y1 - 2018 A1 - Nicola Giuliani A1 - Andrea Mola A1 - Luca Heltai VL - 121 ER - TY - JOUR T1 - Advances in Reduced order modelling for CFD: vortex shedding around a circular cylinder using a POD-Galerkin method JF - Communication in Applied Industrial Mathematics Y1 - 2017 A1 - Giovanni Stabile A1 - Saddam Hijazi A1 - Stefano Lorenzi A1 - Andrea Mola A1 - Gianluigi Rozza KW - finite volume, CFD KW - Reduced order methods AB -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.

UR - https://arxiv.org/abs/1701.03424 ER - TY - JOUR T1 - Curvature-adapted remeshing of CAD surfaces JF - Engineering with Computers Y1 - 2017 A1 - Franco Dassi A1 - Andrea Mola A1 - Hang Si PB - Springer Nature VL - 34 UR - https://doi.org/10.1007/s00366-017-0558-2 ER - TY - JOUR T1 - POD-Galerkin reduced order methods for CFD using Finite Volume Discretisation: vortex shedding around a circular cylinder JF - Communications in Applied and Industrial Mathematics Y1 - 2017 A1 - Giovanni Stabile A1 - Saddam Hijazi A1 - Andrea Mola A1 - Stefano Lorenzi A1 - Gianluigi Rozza PB - Walter de Gruyter {GmbH} VL - 8 UR - https://doi.org/10.1515/caim-2017-0011 ER - TY - JOUR T1 - Wet and Dry Transom Stern Treatment for Unsteady and Nonlinear Potential Flow Model for Naval Hydrodynamics Simulations JF - Journal of Ship Research Y1 - 2017 A1 - Andrea Mola A1 - Luca Heltai A1 - Antonio DeSimone AB -We present a model for the fast evaluation of the total drag of ship hulls operating in both wet and dry transom stern conditions, in calm or wavy water, based on the combination of an unsteady semi-Lagrangian potential flow formulation with fully nonlinear free-surface treatment, experimental correlations, and simplified viscous drag modeling. The implementation is entirely based on open source libraries. The spatial discretization is solved using a streamline upwind Petrov‐Galerkin stabilization of an iso-parametric, collocation based, boundary element method, implemented using the open source library deal.II. The resulting nonlinear differential-algebraic system is integrated in time using implicit backward differentiation formulas, implemented in the open source library SUNDIALS. The Open CASCADE library is used to interface the model directly with computer-aided design data structures. The model accounts automatically for hulls with a transom stern, both in wet and dry regimes, by using a specific treatment of the free-surface nodes on the stern edge that automatically detects when the hull advances at low speeds. In this case, the transom stern is partially immersed, and a pressure patch is applied on the water surface detaching from the transom stern, to recover the gravity effect of the recirculating water on the underlying irrotational flow domain. The parameters of the model used to impose the pressure patch are approximated from experimental relations found in the literature. The test cases considered are those of the U.S. Navy Combatant DTMB-5415 and the National Physical Laboratory hull. Comparisons with experimental data on quasi-steady test cases for both water elevation and total hull drag are presented and discussed. The quality of the results obtained on quasi-steady simulations suggests that this model can represent a promising alternative to current unsteady solvers for simulations with Froude numbers below 0.35.

VL - 61 ER - TY - CONF T1 - Advances in geometrical parametrization and reduced order models and methods for computational fluid dynamics problems in applied sciences and engineering: overview and perspectives T2 - Proceedings of the ECCOMAS Congress 2016, VII European Conference on Computational Methods in Applied Sciences and Engineering, Y1 - 2016 A1 - Filippo Salmoiraghi A1 - Francesco Ballarin A1 - Giovanni Corsi A1 - Andrea Mola A1 - Marco Tezzele A1 - Gianluigi Rozza ED - Papadrakakis, M. ED - Papadopoulos, V. ED - Stefanou, G. ED - Plevris, V. AB -Several problems in applied sciences and engineering require reduction techniques in order to allow computational tools to be employed in the daily practice, especially in iterative procedures such as optimization or sensitivity analysis. Reduced order methods need to face increasingly complex problems in computational mechanics, especially into a multiphysics setting. Several issues should be faced: stability of the approximation, efficient treatment of nonlinearities, uniqueness or possible bifurcations of the state solutions, proper coupling between fields, as well as offline-online computing, computational savings and certification of errors as measure of accuracy. Moreover, efficient geometrical parametrization techniques should be devised to efficiently face shape optimization problems, as well as shape reconstruction and shape assimilation problems. A related aspect deals with the management of parametrized interfaces in multiphysics problems, such as fluid-structure interaction problems, and also a domain decomposition based approach for complex parametrized networks. We present some illustrative industrial and biomedical problems as examples of recent advances on methodological developments.

JF - Proceedings of the ECCOMAS Congress 2016, VII European Conference on Computational Methods in Applied Sciences and Engineering, PB - ECCOMAS CY - Crete, Greece U1 - 35466 U2 - Mathematics U4 - 1 U5 - MAT/08 ER - TY - CONF T1 - Ship Sinkage and Trim Predictions Based on a CAD Interfaced Fully Nonlinear Potential Model T2 - The 26th International Ocean and Polar Engineering Conference Y1 - 2016 A1 - Andrea Mola A1 - Luca Heltai A1 - Antonio DeSimone A1 - Massimiliano Berti JF - The 26th International Ocean and Polar Engineering Conference PB - International Society of Offshore and Polar Engineers VL - 3 ER - TY - JOUR T1 - FEM SUPG stabilisation of mixed isoparametric BEMs: application to linearised free surface flows JF - Engineering Analysis with Boundary Elements 59 (2015), pp. 8-22 Y1 - 2015 A1 - Nicola Giuliani A1 - Andrea Mola A1 - Luca Heltai A1 - L. Formaggia AB -In finite element formulations, transport dominated problems are often stabilised through the Streamline-Upwind-Petrov–Galerkin (SUPG) method. Its application is straightforward when the problem at hand is solved using Galerkin methods. Applications of boundary integral formulations often resort to collocation techniques which are computationally more tractable. In this framework, the Galerkin method and the stabilisation may still be used to successfully apply boundary conditions and resolve instabilities that are frequently observed in transport dominated problems. We apply this technique to an adaptive collocation boundary element method for the solution of stationary potential flows, where we solve a mixed Poisson problem in boundary integral form, with the addition of linearised free surface boundary conditions. We use a mixed boundary element formulation to allow for different finite dimensional spaces describing the flow potential and its normal derivative, and we validate our method simulating the flow around both a submerged body and a surface piercing body. The coupling of mixed surface finite elements and strongly consistent stabilisation techniques with boundary elements opens up the possibility to use non conformal unstructured grids with local refinement, without introducing the inconsistencies of other stabilisation techniques based on up-winding and finite difference schemes.

UR - http://urania.sissa.it/xmlui/handle/1963/34466 U1 - 34640 U2 - Mathematics U4 - 1 U5 - MAT/08 ER - TY - JOUR T1 - Curvature-adapted remeshing of CAD surfaces JF - Procedia Engineering Y1 - 2014 A1 - Franco Dassi A1 - Andrea Mola A1 - Hang Si AB -A common representation of surfaces with complicated topology and geometry is through composite parametric surfaces as is the case for most CAD modelers. A challenging problem is how to generate a mesh of such a surface that well approximates the geometry of the surface, preserves its topology and important geometric features, and contains nicely shaped elements. In this work, we present an optimization-based surface remeshing method that is able to satisfy many of these requirements simultaneously. This method is inspired by the recent work of Lévy and Bonneel (Proc. 21th International Meshing Roundtable, October 2012), which embeds a smooth surface into a high-dimensional space and remesh it uniformly in that embedding space. Our method works directly in the 3d spaces and uses an embedding space in R6 to evaluate mesh size and mesh quality. It generates a curvatureadapted anisotropic surface mesh that well represents the geometry of the surface with a low number of elements. We illustrate our approach through various examples.

PB - Elsevier VL - 82 UR - https://doi.org/10.1016/j.proeng.2014.10.388 N1 - This is an open access article under the CC BY-NC-ND license U1 - 35220 U2 - Mathematics U4 - 1 ER - TY - Generic T1 - A fully nonlinear potential model for ship hydrodynamics directly interfaced with CAD data structures T2 - Proceedings of the 24th International Ocean and Polar Engineering Conference, Busan, 2014 Y1 - 2014 A1 - Andrea Mola A1 - Luca Heltai A1 - Antonio DeSimone KW - ship hydrodynamics AB - We present a model for ship hydrodynamics simulations currently under development at SISSA. The model employs potential flow theory and fully nonlinear free surface boundary conditions. The spatial discretization of the equations is performed by means of a collocation BEM. This gives rise to a Differential Algbraic Equations (DAE) system, solved using an implicit BDF scheme to time advance the solution. The model has been implemented into a C++ software able to automatically generate the computational grids from the CAD geometry of the hull. Numerical results on Kriso KCS and KVLCC2 hulls are presented and discussed. JF - Proceedings of the 24th International Ocean and Polar Engineering Conference, Busan, 2014 PB - SISSA U1 - 7357 U2 - Mathematics U4 - 1 U5 - MAT/08 ANALISI NUMERICA ER - TY - CONF T1 - Potential Model for Ship Hydrodynamics Simulations Directly Interfaced with CAD Data Structures T2 - The 24th International Ocean and Polar Engineering Conference Y1 - 2014 A1 - Andrea Mola A1 - Luca Heltai A1 - Antonio DeSimone A1 - Massimiliano Berti JF - The 24th International Ocean and Polar Engineering Conference PB - International Society of Offshore and Polar Engineers VL - 4 ER - TY - JOUR T1 - A stable and adaptive semi-Lagrangian potential model for unsteady and nonlinear ship-wave interactions JF - Engineering Analysis with Boundary Elements, 37(1):128 – 143, 2013. Y1 - 2013 A1 - Andrea Mola A1 - Luca Heltai A1 - Antonio DeSimone KW - Unsteady ship-wave interaction AB -We present an innovative numerical discretization of the equations of inviscid potential flow for the simulation of three dimensional unsteady and nonlinear water waves generated by a ship hull advancing in water. The equations of motion are written in a semi-Lagrangian framework, and the resulting integro-diff erential equations are discretized in space via an adaptive iso-parametric collocation Boundary Element Method, and in time via adaptive implicit Backward Di erentiation Formulas (BDF) with variable step and variable order. When the velocity of the advancing ship hull is non-negligible, the semi-Lagrangian formulation (also known as Arbitrary Lagrangian Eulerian formulation, or ALE) of the free surface equations contains dominant transport terms which are stabilized with a Streamwise Upwind Petrov-Galerkin (SUPG) method. The SUPG stabilization allows automatic and robust adaptation of the spatial discretization with unstructured quadrilateral grids. Preliminary results are presented where we compare our numerical model with experimental results on the case of a Wigley hull advancing in calm water with fi xed sink and trim.

PB - SISSA UR - http://hdl.handle.net/1963/5669 U1 - 5457 U2 - Mathematics U3 - Functional Analysis and Applications U4 - -1 ER - TY - CONF T1 - A stable semi-lagrangian potential method for the simulation of ship interaction with unsteady and nonlinear waves T2 - 17th Int. Conf. Ships Shipp. Res. Y1 - 2012 A1 - Andrea Mola A1 - Luca Heltai A1 - Antonio DeSimone JF - 17th Int. Conf. Ships Shipp. Res. ER - TY - JOUR T1 - Multi-physics modelling and sensitivity analysis of olympic rowing boat dynamics JF - Sports Engineering Y1 - 2011 A1 - Andrea Mola A1 - Mehdi Ghommem A1 - Muhammad R. Hajj PB - Springer Nature VL - 14 UR - https://doi.org/10.1007/s12283-011-0075-2 ER - TY - JOUR T1 - A three-dimensional model for the dynamics and hydrodynamics of rowing boats JF - Proceedings of the Institution of Mechanical Engineers, Part P: Journal of Sports Engineering and Technology Y1 - 2010 A1 - L. Formaggia A1 - Andrea Mola A1 - N Parolini A1 - M Pischiutta AB -This paper proposes a new model describing the dynamics of a rowing boat for general three-dimensional motions. The complex interaction between the different components of the rowers–-oars–-boat system is analysed and reduced to a set of ordinary differential equations governing the rigid motion along the six degrees of freedom. To treat the unstable nature of the physical problem, a rather simple (but effective) control model is included, which mimics the main active control techniques adopted by the rowers during their action.

VL - 224 UR - https://doi.org/10.1243/17543371jset46 ER - TY - JOUR T1 - Low-Frequency Variations of Force Coefficients on Square Cylinders with Sharp and Rounded Corners JF - Journal of Structural Engineering Y1 - 2009 A1 - Andrea Mola A1 - Giancarlo Bordonaro A1 - Muhammad R. Hajj PB - American Society of Civil Engineers ({ASCE}) VL - 135 UR - https://doi.org/10.1061/(asce)st.1943-541x.0000034 ER - TY - JOUR T1 - A model for the dynamics of rowing boats JF - International Journal for Numerical Methods in Fluids Y1 - 2009 A1 - L. Formaggia A1 - Edie Miglio A1 - Andrea Mola A1 - Antonio Montano PB - Wiley VL - 61 UR - https://doi.org/10.1002/fld.1940 ER - TY - JOUR T1 - Fluid–structure interaction problems in free surface flows: Application to boat dynamics JF - International Journal for Numerical Methods in Fluids Y1 - 2008 A1 - L. Formaggia A1 - Edie Miglio A1 - Andrea Mola A1 - N Parolini PB - Wiley VL - 56 UR - https://doi.org/10.1002/fld.1583 ER - TY - JOUR T1 - Calculation of impulsively started incompressible viscous flows JF - Int. J. Numer. Meth. Fluids Y1 - 2004 A1 - Marra, Andrea A1 - Andrea Mola A1 - Quartapelle, Luigi A1 - Riviello, Luca VL - 46 ER -