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.

%B VIII International Conference on Computational Methods in Marine Engineering %G eng %0 Conference Paper %B VIII International Conference on Computational Methods in Marine Engineering %D 2019 %T Efficient Reduction in Shape Parameter Space Dimension for Ship Propeller Blade Design %A Andrea Mola %A Marco Tezzele %A Mahmoud Gadalla %A Valdenazzi, Federica %A Grassi, Davide %A Padovan, Roberta %A Gianluigi Rozza %XIn 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.

%B VIII International Conference on Computational Methods in Marine Engineering %G eng %0 Journal Article %J Advanced Modeling and Simulation in Engineering Sciences %D 2018 %T Dimension reduction in heterogeneous parametric spaces with application to naval engineering shape design problems %A Marco Tezzele %A Filippo Salmoiraghi %A Andrea Mola %A Gianluigi Rozza %XWe 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.

%B Advanced Modeling and Simulation in Engineering Sciences %V 5 %P 25 %8 Sep %G eng %R 10.1186/s40323-018-0118-3 %0 Conference Proceedings %B The 28th International Ocean and Polar Engineering Conference %D 2018 %T An efficient shape parametrisation by free-form deformation enhanced by active subspace for hull hydrodynamic ship design problems in open source environment %A Nicola Demo %A Marco Tezzele %A Andrea Mola %A Gianluigi Rozza %K Active subspaces %K Boundary element method %K Dynamic mode decomposition %K Fluid structure interaction %K Free form deformation %K Fully nonlinear potential %K Numerical towing tank %X 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. %B The 28th International Ocean and Polar Engineering Conference %I International Society of Offshore and Polar Engineers %C Sapporo, Japan %G eng %U https://www.onepetro.org/conference-paper/ISOPE-I-18-481 %0 Conference Paper %B Technology and Science for the Ships of the Future: Proceedings of NAV 2018: 19th International Conference on Ship & Maritime Research %D 2018 %T Model Order Reduction by means of Active Subspaces and Dynamic Mode Decomposition for Parametric Hull Shape Design Hydrodynamics %A Marco Tezzele %A Nicola Demo %A Mahmoud Gadalla %A Andrea Mola %A Gianluigi Rozza %X 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. %B Technology and Science for the Ships of the Future: Proceedings of NAV 2018: 19th International Conference on Ship & Maritime Research %I IOS Press %C Trieste, Italy %G eng %U http://ebooks.iospress.nl/publication/49270 %R 10.3233/978-1-61499-870-9-569 %0 Journal Article %J Advances in Engineering Software %D 2018 %T π-BEM : A flexible parallel implementation for adaptive , geometry aware , and high order boundary element methods %A Nicola Giuliani %A Andrea Mola %A Luca Heltai %B Advances in Engineering Software %V 121 %P 39–58 %G eng %0 Journal Article %J Communication in Applied Industrial Mathematics %D 2017 %T Advances in Reduced order modelling for CFD: vortex shedding around a circular cylinder using a POD-Galerkin method %A Giovanni Stabile %A Saddam Hijazi %A Stefano Lorenzi %A Andrea Mola %A Gianluigi Rozza %K finite volume, CFD %K Reduced order methods %XVortex 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.

%B Communication in Applied Industrial Mathematics %G eng %U https://arxiv.org/abs/1701.03424 %9 reviewed %0 Journal Article %J Engineering with Computers %D 2017 %T Curvature-adapted remeshing of CAD surfaces %A Franco Dassi %A Andrea Mola %A Hang Si %B Engineering with Computers %I Springer Nature %V 34 %P 565–576 %8 dec %G eng %U https://doi.org/10.1007/s00366-017-0558-2 %R 10.1007/s00366-017-0558-2 %0 Journal Article %J Communications in Applied and Industrial Mathematics %D 2017 %T POD-Galerkin reduced order methods for CFD using Finite Volume Discretisation: vortex shedding around a circular cylinder %A Giovanni Stabile %A Saddam Hijazi %A Andrea Mola %A Stefano Lorenzi %A Gianluigi Rozza %B Communications in Applied and Industrial Mathematics %I Walter de Gruyter {GmbH} %V 8 %P 210–236 %8 dec %G eng %U https://doi.org/10.1515/caim-2017-0011 %R 10.1515/caim-2017-0011 %0 Journal Article %J Journal of Ship Research %D 2017 %T Wet and Dry Transom Stern Treatment for Unsteady and Nonlinear Potential Flow Model for Naval Hydrodynamics Simulations %A Andrea Mola %A Luca Heltai %A Antonio DeSimone %XWe 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.

%B Journal of Ship Research %V 61 %P 1–14 %G eng %R https://doi.org/10.5957/JOSR.61.1.160016 %0 Conference Paper %B Proceedings of the ECCOMAS Congress 2016, VII European Conference on Computational Methods in Applied Sciences and Engineering, %D 2016 %T Advances in geometrical parametrization and reduced order models and methods for computational fluid dynamics problems in applied sciences and engineering: overview and perspectives %A Filippo Salmoiraghi %A Francesco Ballarin %A Giovanni Corsi %A Andrea Mola %A Marco Tezzele %A Gianluigi Rozza %E Papadrakakis, M. %E Papadopoulos, V. %E Stefanou, G. %E Plevris, V. %XSeveral 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.

%B Proceedings of the ECCOMAS Congress 2016, VII European Conference on Computational Methods in Applied Sciences and Engineering, %I ECCOMAS %C Crete, Greece %8 06/2016 %G en %1 35466 %2 Mathematics %4 1 %# MAT/08 %$ Submitted by Gianluigi Rozza (grozza@sissa.it) on 2016-05-13T00:28:01Z No. of bitstreams: 1 eccomas2016_AROMA.pdf: 1846196 bytes, checksum: 9636e713df80de178d87fd2feff76f91 (MD5) %0 Conference Paper %B The 26th International Ocean and Polar Engineering Conference %D 2016 %T Ship Sinkage and Trim Predictions Based on a CAD Interfaced Fully Nonlinear Potential Model %A Andrea Mola %A Luca Heltai %A Antonio DeSimone %A Massimiliano Berti %B The 26th International Ocean and Polar Engineering Conference %I International Society of Offshore and Polar Engineers %V 3 %P 511–518 %G eng %0 Journal Article %J Engineering Analysis with Boundary Elements 59 (2015), pp. 8-22 %D 2015 %T FEM SUPG stabilisation of mixed isoparametric BEMs: application to linearised free surface flows %A Nicola Giuliani %A Andrea Mola %A Luca Heltai %A L. Formaggia %XIn 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.

%B Engineering Analysis with Boundary Elements 59 (2015), pp. 8-22 %G en_US %U http://urania.sissa.it/xmlui/handle/1963/34466 %1 34640 %2 Mathematics %4 1 %# MAT/08 %$ Submitted by ngiuliani@sissa.it (ngiuliani@sissa.it) on 2015-05-04T12:32:04Z No. of bitstreams: 1 paper-free-surface.pdf: 888763 bytes, checksum: 68366266786e84b78a356acd6de50840 (MD5) %R 10.1016/j.enganabound.2015.04.006 %0 Journal Article %J Procedia Engineering %D 2014 %T Curvature-adapted remeshing of CAD surfaces %A Franco Dassi %A Andrea Mola %A Hang Si %XA 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.

%B Procedia Engineering %I Elsevier %V 82 %P 253–265 %G en %U https://doi.org/10.1016/j.proeng.2014.10.388 %1 35220 %2 Mathematics %4 1 %$ Submitted by gfeltrin@sissa.it (gfeltrin@sissa.it) on 2015-11-11T16:19:16Z No. of bitstreams: 1 art2014.pdf: 3913220 bytes, checksum: ba78b55b28c0eed195d2abc48bd7ca4b (MD5) %R 10.1016/j.proeng.2014.10.388 %0 Conference Proceedings %B Proceedings of the 24th International Ocean and Polar Engineering Conference, Busan, 2014 %D 2014 %T A fully nonlinear potential model for ship hydrodynamics directly interfaced with CAD data structures %A Andrea Mola %A Luca Heltai %A Antonio DeSimone %K ship hydrodynamics %X 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. %B Proceedings of the 24th International Ocean and Polar Engineering Conference, Busan, 2014 %I SISSA %G en %1 7357 %2 Mathematics %4 1 %# MAT/08 ANALISI NUMERICA %$ Submitted by Luca Heltai (heltai@sissa.it) on 2014-03-22T15:00:37Z No. of bitstreams: 1 MolaHeltaiDeSimone2014.pdf: 1241007 bytes, checksum: 1a56d909a08f1f470c7efc14018e37f8 (MD5) %0 Conference Paper %B The 24th International Ocean and Polar Engineering Conference %D 2014 %T Potential Model for Ship Hydrodynamics Simulations Directly Interfaced with CAD Data Structures %A Andrea Mola %A Luca Heltai %A Antonio DeSimone %A Massimiliano Berti %B The 24th International Ocean and Polar Engineering Conference %I International Society of Offshore and Polar Engineers %V 4 %P 815–822 %G eng %0 Journal Article %J Engineering Analysis with Boundary Elements, 37(1):128 – 143, 2013. %D 2013 %T A stable and adaptive semi-Lagrangian potential model for unsteady and nonlinear ship-wave interactions %A Andrea Mola %A Luca Heltai %A Antonio DeSimone %K Unsteady ship-wave interaction %XWe 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.

%B Engineering Analysis with Boundary Elements, 37(1):128 – 143, 2013. %I SISSA %G en %U http://hdl.handle.net/1963/5669 %1 5457 %2 Mathematics %3 Functional Analysis and Applications %4 -1 %$ Submitted by Maria Pia Calandra (calapia@sissa.it) on 2012-03-16T11:19:51Z\\nNo. of bitstreams: 1\\nMolaHeltaiDeSimone-2012-aa.pdf: 1981527 bytes, checksum: 8a82da009242da0cf884e2c7f04009a6 (MD5) %0 Conference Paper %B 17th Int. Conf. Ships Shipp. Res. %D 2012 %T A stable semi-lagrangian potential method for the simulation of ship interaction with unsteady and nonlinear waves %A Andrea Mola %A Luca Heltai %A Antonio DeSimone %B 17th Int. Conf. Ships Shipp. Res. %G eng %0 Journal Article %J Sports Engineering %D 2011 %T Multi-physics modelling and sensitivity analysis of olympic rowing boat dynamics %A Andrea Mola %A Mehdi Ghommem %A Muhammad R. Hajj %B Sports Engineering %I Springer Nature %V 14 %P 85–94 %8 nov %G eng %U https://doi.org/10.1007/s12283-011-0075-2 %R 10.1007/s12283-011-0075-2 %0 Journal Article %J Proceedings of the Institution of Mechanical Engineers, Part P: Journal of Sports Engineering and Technology %D 2010 %T A three-dimensional model for the dynamics and hydrodynamics of rowing boats %A L. Formaggia %A Andrea Mola %A N Parolini %A M Pischiutta %XThis 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.

%B Proceedings of the Institution of Mechanical Engineers, Part P: Journal of Sports Engineering and Technology %V 224 %P 51-61 %G eng %U https://doi.org/10.1243/17543371jset46 %R 10.1243/17543371jset46 %0 Journal Article %J Journal of Structural Engineering %D 2009 %T Low-Frequency Variations of Force Coefficients on Square Cylinders with Sharp and Rounded Corners %A Andrea Mola %A Giancarlo Bordonaro %A Muhammad R. Hajj %B Journal of Structural Engineering %I American Society of Civil Engineers ({ASCE}) %V 135 %P 828–835 %8 jul %G eng %U https://doi.org/10.1061/(asce)st.1943-541x.0000034 %R 10.1061/(asce)st.1943-541x.0000034 %0 Journal Article %J International Journal for Numerical Methods in Fluids %D 2009 %T A model for the dynamics of rowing boats %A L. Formaggia %A Edie Miglio %A Andrea Mola %A Antonio Montano %B International Journal for Numerical Methods in Fluids %I Wiley %V 61 %P 119–143 %8 sep %G eng %U https://doi.org/10.1002/fld.1940 %R 10.1002/fld.1940 %0 Journal Article %J International Journal for Numerical Methods in Fluids %D 2008 %T Fluid–structure interaction problems in free surface flows: Application to boat dynamics %A L. Formaggia %A Edie Miglio %A Andrea Mola %A N Parolini %B International Journal for Numerical Methods in Fluids %I Wiley %V 56 %P 965–978 %G eng %U https://doi.org/10.1002/fld.1583 %R 10.1002/fld.1583 %0 Journal Article %J Int. J. Numer. Meth. Fluids %D 2004 %T Calculation of impulsively started incompressible viscous flows %A Marra, Andrea %A Andrea Mola %A Quartapelle, Luigi %A Riviello, Luca %B Int. J. Numer. Meth. Fluids %V 46 %P 877–902 %G eng