@article {2021, title = {On the comparison of LES data-driven reduced order approaches for hydroacoustic analysis}, journal = {Computers \& Fluids}, volume = {216}, year = {2021}, pages = {104819}, abstract = {

In this work, Dynamic Mode Decomposition (DMD) and Proper Orthogonal Decomposition (POD) methodologies are applied to hydroacoustic dataset computed using Large Eddy Simulation (LES) coupled with Ffowcs Williams and Hawkings (FWH) analogy. First, a low-dimensional description of the flow fields is presented with modal decomposition analysis. Sensitivity towards the DMD and POD bases truncation rank is discussed, and extensive dataset is provided to demonstrate the ability of both algorithms to reconstruct the flow fields with all the spatial and temporal frequencies necessary to support accurate noise evaluation. Results show that while DMD is capable to capture finer coherent structures in the wake region for the same amount of employed modes, reconstructed flow fields using POD exhibit smaller magnitudes of global spatiotemporal errors compared with DMD counterparts. Second, a separate set of DMD and POD modes generated using half the snapshots is employed into two data-driven reduced models respectively, based on DMD mid cast and POD with Interpolation (PODI). In that regard, results confirm that the predictive character of both reduced approaches on the flow fields is sufficiently accurate, with a relative superiority of PODI results over DMD ones. This infers that, discrepancies induced due to interpolation errors in PODI is relatively low compared with errors induced by integration and linear regression operations in DMD, for the present setup. Finally, a post processing analysis on the evaluation of FWH acoustic signals utilizing reduced fluid dynamic fields as input demonstrates that both DMD and PODI data-driven reduced models are efficient and sufficiently accurate in predicting acoustic noises.

}, keywords = {Dynamic mode decomposition, Ffowcs Williams and Hawkings, Hydroacoustics, Large eddy simulation, Model reduction, Proper orthogonal decomposition}, issn = {0045-7930}, doi = {https://doi.org/10.1016/j.compfluid.2020.104819}, url = {https://www.sciencedirect.com/science/article/pii/S0045793020303893}, author = {Mahmoud Gadalla and Marta Cianferra and Marco Tezzele and Giovanni Stabile and Andrea Mola and Gianluigi Rozza} } @article {2021, title = {Hull Shape Design Optimization with Parameter Space and Model Reductions, and Self-Learning Mesh Morphing}, journal = {Journal of Marine Science and Engineering}, volume = {9}, year = {2021}, pages = {185}, abstract = {

In the field of parametric partial differential equations, shape optimization represents a challenging problem due to the required computational resources. In this contribution, a data-driven framework involving multiple reduction techniques is proposed to reduce such computational burden. Proper orthogonal decomposition (POD) and active subspace genetic algorithm (ASGA) are applied for a dimensional reduction of the original (high fidelity) model and for an efficient genetic optimization based on active subspace property. The parameterization of the shape is applied directly to the computational mesh, propagating the generic deformation map applied to the surface (of the object to optimize) to the mesh nodes using a radial basis function (RBF) interpolation. Thus, topology and quality of the original mesh are preserved, enabling application of POD-based reduced order modeling techniques, and avoiding the necessity of additional meshing steps. Model order reduction is performed coupling POD and Gaussian process regression (GPR) in a data-driven fashion. The framework is validated on a benchmark ship.

}, issn = {2077-1312}, doi = {10.3390/jmse9020185}, url = {https://www.mdpi.com/2077-1312/9/2/185}, author = {Nicola Demo and Marco Tezzele and Andrea Mola and Gianluigi Rozza} } @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} } @article {tezzele2020pygem, title = {PyGeM: Python Geometrical Morphing}, journal = {Software Impacts}, volume = {7}, year = {2021}, pages = {100047}, abstract = {PyGeM is an open source Python package which allows to easily parametrize and deform 3D object described by CAD files or 3D meshes. It implements several morphing techniques such as free form deformation, radial basis function interpolation, and inverse distance weighting. Due to its versatility in dealing with different file formats it is particularly suited for researchers and practitioners both in academia and in industry interested in computational engineering simulations and optimization studies.}, keywords = {Free form deformation, Geometrical morphing, Inverse distance weighting, Python, Radial basis functions interpolation}, issn = {2665-9638}, doi = {10.1016/j.simpa.2020.100047}, author = {Marco Tezzele and Nicola Demo and Andrea Mola and Gianluigi Rozza} } @conference {2020, title = {Advances in reduced order methods for parametric industrial problems in computational fluid dynamics}, booktitle = {Proceedings of the 6th European Conference on Computational Mechanics: Solids, Structures and Coupled Problems, ECCM 2018 and 7th European Conference on Computational Fluid Dynamics, ECFD 2018}, year = {2020}, abstract = {

Reduced order modeling has gained considerable attention in recent decades owing to the advantages offered in reduced computational times and multiple solutions for parametric problems. The focus of this manuscript is the application of model order reduction techniques in various engineering and scientific applications including but not limited to mechanical, naval and aeronautical engineering. The focus here is kept limited to computational fluid mechanics and related applications. The advances in the reduced order modeling with proper orthogonal decomposition and reduced basis method are presented as well as a brief discussion of dynamic mode decomposition and also some present advances in the parameter space reduction. Here, an overview of the challenges faced and possible solutions are presented with examples from various problems.

}, url = {https://www.scopus.com/inward/record.uri?eid=2-s2.0-85075395686\&partnerID=40\&md5=fb0b1a3cfdfd35a104db9921bc9be675}, author = {Gianluigi Rozza and M.H. Malik and Nicola Demo and Marco Tezzele and Michele Girfoglio and Giovanni Stabile and Andrea Mola} } @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 {2020, title = {Enhancing CFD predictions in shape design problems by model and parameter space reduction}, journal = {Advanced Modeling and Simulation in Engineering Sciences}, volume = {7}, year = {2020}, abstract = {

In this work we present an advanced computational pipeline for the approximation and prediction of the lift coefficient of a parametrized airfoil profile. The non-intrusive reduced order method is based on dynamic mode decomposition (DMD) and it is coupled with dynamic active subspaces (DyAS) to enhance the future state prediction of the target function and reduce the parameter space dimensionality. The pipeline is based on high-fidelity simulations carried out by the application of finite volume method for turbulent flows, and automatic mesh morphing through radial basis functions interpolation technique. The proposed pipeline is able to save 1/3 of the overall computational resources thanks to the application of DMD. Moreover exploiting DyAS and performing the regression on a lower dimensional space results in the reduction of the relative error in the approximation of the time-varying lift coefficient by a factor 2 with respect to using only the DMD.

}, doi = {https://doi.org/10.1186/s40323-020-00177-y}, url = {https://arxiv.org/abs/2001.05237}, author = {Marco Tezzele and Nicola Demo and Giovanni Stabile and Andrea Mola and Gianluigi Rozza} } @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 {gadalla19bladex, title = {BladeX: Python Blade Morphing}, journal = {The Journal of Open Source Software}, volume = {4}, number = {34}, year = {2019}, pages = {1203}, doi = {10.21105/joss.01203}, author = {Mahmoud Gadalla and Marco Tezzele and Andrea Mola and Gianluigi Rozza} } @conference {2019, title = {A complete data-driven framework for the efficient solution of parametric shape design and optimisation in naval engineering problems}, booktitle = {8th International Conference on Computational Methods in Marine Engineering, MARINE 2019}, year = {2019}, abstract = {

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 {\textemdash} assuming the topology is inaltered by the deformation {\textemdash}, 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.

}, url = {https://www.scopus.com/inward/record.uri?eid=2-s2.0-85075342565\&partnerID=40\&md5=d76b8a1290053e7a84fb8801c0e6bb3d}, author = {Nicola Demo and Marco Tezzele and Andrea Mola and Gianluigi Rozza} } @conference {2019, title = {A complete data-driven framework for the efficient solution of parametric shape design and optimisation in naval engineering problems}, booktitle = {VIII International Conference on Computational Methods in Marine Engineering}, year = {2019}, abstract = {

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 {\textendash}- assuming the topology is inaltered by the deformation {\textendash}-, 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.

}, url = {https://arxiv.org/abs/1905.05982}, author = {Nicola Demo and Marco Tezzele and Andrea Mola and Gianluigi Rozza} } @conference {2019, title = {Efficient reduction in shape parameter space dimension for ship propeller blade design}, booktitle = {8th International Conference on Computational Methods in Marine Engineering, MARINE 2019}, year = {2019}, abstract = {

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.

}, url = {https://www.scopus.com/inward/record.uri?eid=2-s2.0-85075395143\&partnerID=40\&md5=b6aa0fcedc2f88e78c295d0f437824d0}, author = {Andrea Mola and Marco Tezzele and Mahmoud Gadalla and Valdenazzi, Federica and Grassi, Davide and Padovan, Roberta and Gianluigi Rozza} } @article {TezzeleSalmoiraghiMolaRozza2018, title = {Dimension reduction in heterogeneous parametric spaces with application to naval engineering shape design problems}, journal = {Advanced Modeling and Simulation in Engineering Sciences}, volume = {5}, number = {1}, year = {2018}, month = {Sep}, pages = {25}, abstract = {

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.

}, doi = {10.1186/s40323-018-0118-3}, author = {Marco Tezzele and Filippo Salmoiraghi and Andrea Mola and Gianluigi Rozza} } @proceedings {demo2018efficient, title = {An efficient shape parametrisation by free-form deformation enhanced by active subspace for hull hydrodynamic ship design problems in open source environment}, year = {2018}, publisher = {International Society of Offshore and Polar Engineers}, address = {Sapporo, Japan}, abstract = {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.}, keywords = {Active subspaces, Boundary element method, Dynamic mode decomposition, Fluid structure interaction, Free form deformation, Fully nonlinear potential, Numerical towing tank}, issn = {978-1-880653-87-6}, url = {https://www.onepetro.org/conference-paper/ISOPE-I-18-481}, author = {Nicola Demo and Marco Tezzele and Andrea Mola and Gianluigi Rozza} } @conference {tezzele2018model, title = {Model Order Reduction by means of Active Subspaces and Dynamic Mode Decomposition for Parametric Hull Shape Design Hydrodynamics}, booktitle = {Technology and Science for the Ships of the Future: Proceedings of NAV 2018: 19th International Conference on Ship \& Maritime Research}, year = {2018}, publisher = {IOS Press}, organization = {IOS Press}, address = {Trieste, Italy}, abstract = {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.}, doi = {10.3233/978-1-61499-870-9-569}, url = {http://ebooks.iospress.nl/publication/49270}, author = {Marco Tezzele and Nicola Demo and Mahmoud Gadalla and Andrea Mola and Gianluigi Rozza} } @conference {2018, title = {SRTP 2.0 - The evolution of the safe return to port concept}, booktitle = {Technology and Science for the Ships of the Future - Proceedings of NAV 2018: 19th International Conference on Ship and Maritime Research}, year = {2018}, abstract = {

In 2010 IMO (International Maritime Organisation) introduced new rules in SOLAS with the aim of intrinsically increase the safety of passenger ships. This requirement is achieved by providing safe areas for passengers and essential services for allowing ship to Safely Return to Port (SRtP). The entry into force of these rules has changed the way to design passenger ships. In this respect big effort in the research has been done by industry to address design issues related to the impact on failure analysis of the complex interactions among systems. Today the research activity is working to bring operational matters in the design stage. This change of research focus was necessary because human factor and the way to operate the ship itself after a casualty on board may have a big impact in the design of the ship/systems. Also the management of the passengers after a casualty is becoming a major topic for safety. This paper presents the state of the art of Italian knowledge in the field of system engineering applied to passenger ship address to safety improvement and design reliability. An overview of present tools and methodologies will be offered together with future focuses in the research activity.

}, doi = {10.3233/978-1-61499-870-9-665}, author = {D. Cangelosi and A. Bonvicini and M. Nardo and Andrea Mola and A. Marchese and Marco Tezzele and Gianluigi Rozza} } @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 {Dassi2017, title = {Curvature-adapted remeshing of {CAD} surfaces}, journal = {Engineering with Computers}, volume = {34}, number = {3}, year = {2017}, month = {dec}, pages = {565{\textendash}576}, publisher = {Springer Nature}, doi = {10.1007/s00366-017-0558-2}, url = {https://doi.org/10.1007/s00366-017-0558-2}, author = {Franco Dassi and Andrea Mola and Hang Si} } @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 {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} } @conference {2016, title = {Advances in geometrical parametrization and reduced order models and methods for computational fluid dynamics problems in applied sciences and engineering: overview and perspectives}, booktitle = {Proceedings of the ECCOMAS Congress 2016, VII European Conference on Computational Methods in Applied Sciences and Engineering,}, year = {2016}, month = {06/2016}, publisher = {ECCOMAS}, organization = {ECCOMAS}, address = {Crete, Greece}, abstract = {

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.

}, author = {Filippo Salmoiraghi and F. Ballarin and Giovanni Corsi and Andrea Mola and Marco Tezzele and Gianluigi Rozza}, editor = {Papadrakakis, M. and Papadopoulos, V. and Stefanou, G. and Plevris, V.} } @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 {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 {2014, title = {Curvature-adapted remeshing of CAD surfaces}, journal = {Procedia Engineering}, volume = {82}, number = {Procedia Engineering;volume 82; issue C; pages 253-265;}, year = {2014}, note = {This is an open access article under the CC BY-NC-ND license}, pages = {253{\textendash}265}, publisher = {Elsevier}, abstract = {

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{\'e}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.

}, doi = {10.1016/j.proeng.2014.10.388}, url = {https://doi.org/10.1016/j.proeng.2014.10.388}, author = {Franco Dassi and Andrea Mola and Hang Si} } @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} } @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 {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} } @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 {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 {FormaggiaEtAlRowing2010, title = {A three-dimensional model for the dynamics and hydrodynamics of rowing boats}, journal = {Proceedings of the Institution of Mechanical Engineers, Part P: Journal of Sports Engineering and Technology}, volume = {224}, number = {1}, year = {2010}, pages = {51-61}, abstract = {

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{\textendash}-oars{\textendash}-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.

}, doi = {10.1243/17543371jset46}, url = {https://doi.org/10.1243/17543371jset46}, author = {L. Formaggia and Andrea Mola and N Parolini and M Pischiutta} } @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 {Formaggia2009, title = {A model for the dynamics of rowing boats}, journal = {International Journal for Numerical Methods in Fluids}, volume = {61}, number = {2}, year = {2009}, month = {sep}, pages = {119{\textendash}143}, publisher = {Wiley}, doi = {10.1002/fld.1940}, url = {https://doi.org/10.1002/fld.1940}, author = {L. Formaggia and Edie Miglio and Andrea Mola and Antonio Montano} } @article {Formaggia2008, title = {Fluid{\textendash}structure interaction problems in free surface flows: Application to boat dynamics}, journal = {International Journal for Numerical Methods in Fluids}, volume = {56}, number = {8}, year = {2008}, pages = {965{\textendash}978}, publisher = {Wiley}, doi = {10.1002/fld.1583}, url = {https://doi.org/10.1002/fld.1583}, author = {L. Formaggia and Edie Miglio and Andrea Mola and N Parolini} } @article {marraEtAl2004, title = {Calculation of impulsively started incompressible viscous flows}, journal = {Int. J. Numer. Meth. Fluids}, volume = {46}, year = {2004}, pages = {877{\textendash}902}, author = {Marra, Andrea and Andrea Mola and Quartapelle, Luigi and Riviello, Luca} }