Access the full text.
Sign up today, get DeepDyve free for 14 days.
In this work, Dynamic Mode Decomposition (DMD) and Proper Orthogonal Decomposi- tion (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 ow elds is presented with modal decomposition analy- sis. 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 ow elds with all the spatial and temporal frequencies necessary to support accurate noise evaluation. Results show that while DMD is capable to capture ner coherent structures in the wake region for the same amount of employed modes, reconstructed ow elds using POD exhibit smaller magnitudes of global spatiotemporal errors compared with DMD counterparts. Sec- ond, 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 con rm that the predictive character of both re- duced approaches on the ow elds is suciently 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 op- erations in DMD, for the present setup. Finally, a post processing analysis on the evaluation of FWH acoustic signals utilizing reduced uid dynamic elds as input demonstrates that both DMD and PODI data-driven reduced models are ecient and suciently accurate in predicting acoustic noises. Contents 1 Introduction 2 2 Methodology 4 2.1 Full order model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 2.2 Acoustic model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 2.3 Dynamic mode decomposition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7 2.4 Proper orthogonal decomposition with interpolation . . . . . . . . . . . . . . . . . 8 mahmoud.gadalla@sissa.it marta.cianferra@dia.units.it marco.tezzele@sissa.it giovanni.stabile@sissa.it andrea.mola@sissa.it gianluigi.rozza@sissa.it arXiv:2006.14428v2 [math.NA] 28 Dec 2020 3 Numerical results 9 3.1 Full order CFD . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 3.2 Modal decomposition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 3.2.1 Singular values decay . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 3.2.2 Modal representation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 3.2.3 Associated coecients . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 3.3 Flow elds reconstruction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 3.3.1 Global error . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 3.3.2 Fields visualization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 3.3.3 Error statistics in the wake region . . . . . . . . . . . . . . . . . . . . . . . 16 3.4 Fields mid cast using DMD and PODI . . . . . . . . . . . . . . . . . . . . . . . . . 18 3.4.1 Prediction error analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 3.4.2 Coherent structures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 3.4.3 Drag and lift coecients . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 3.5 Spectral analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22 3.5.1 Data probes inside and outside the wake region . . . . . . . . . . . . . . . . 22 3.5.2 Fast Fourier transform . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24 3.6 Acoustic analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26 4 Conclusions and perspectives 28 A ROMs cross-comparison 29 1 Introduction In several engineering elds, there has been recently a growing need to include uid dynamic per- formance evaluation criteria associated with acoustic emissions. This is, for example, the case in aircraft engine design, car manufacturing, and ship design optimization. A particular motiva- tion behind the present research article is to investigate and propose a methodology that can be successively used for the noise level prediction of naval propellers since the early design process [16, 42]. The need for a reduction of the acoustic emissions through ship design optimization usually involves virtual prototyping and parametric high delity simulations. Thanks to the increase of the available computational resources, a deeper insight towards the complex physics associated with hydroacoustic phenomena has become nowadays aordable with unprecedented spatial and tem- poral scales (see, for example, wall-resolving Large Eddy Simulation (LES) [67, 48]). However, the enormous data sizes resulting from such simulations pose several challenges on the input/output operations, post-processing, or the long-term data storage. On the other hand, hybrid techniques such as, among others, Detached-Eddy Simulation or wall-layer model LES (WLES), have allowed obtaining eddy resolving eld data with a reasonable use of computational resources. These tech- niques are, however, still expensive for an early stage design process, when a number of dierent geometric con gurations has to be rapidly analysed to restrict the range of variation of the principal design parameters. Therefore, seeking a suitable data compression strategy that allows extracting the most relevant and revealing information in a reduced order manner, hence providing quick access as well as ecient data storage, becomes a crucial asset. Through multidisciplinary scale, eorts have been made to realize optimal shape design for underwater noise sources, including ship hulls [13, 98, 93] and propellers [62, 99], using ecient geometrical parameterization techniques [95, 32]. On the acoustic side, the development of new gen- eration noise prediction tools was considered a major focus in this work. Particularly, based on the Ffowcs Williams and Hawkings (FWH) analogy [29], several improvements have been developed. For example, Cianferra et al. [15] compared several implementations of the non-linear quadrupole term, highlighting its signi cant contribution to the overall hydrodynamic noise emissions in wide range of frequencies. In a companion paper [14], they showed the eect of shape deformation on the radiated noise for elementary geometries. On a more engineering level, the generated hydro- dynamic noise from a benchmark marine propeller was evaluated in open sea conditions [16, 17] using FWH coupled with LES. The aforementioned hydroacoustic models are typically described by a system of non-linear Partial Dierential Equations (PDEs), the resolution of which results in the uid dynamic elds 2 necessary to reproduce the noise source for an acoustic predictions. In fact, a reliable reconstruction of the noise source is crucially dependent on the ow structures and the resolved spatial and temporal scales of the uid dynamic elds (for a discussion, see [8]). Besides LES works [86, 2, 65, 48, 14, 17, 15], several uid dynamic models have been employed in marine hydroacoustics. To name a few, we mention the potential ow theory [47], boundary element method [85], Reynolds Averaged Navier-Stokes equations (RANS) [36, 42, 1], DES [57], and Direct Numerical Simulation (DNS) [19, 84, 79]. Being regarded as an optimal advance between RANS and DNS, recent literature has reported LES to be the most suitable model which reproduces the noise source with high level of realism [65, 2, 41]. The resolution of the described system of PDEs using standard discretization methods ( nite elements, nite volumes, nite dierences), which we will refer hereafter as the Full Order Model (FOM), allows for high delity acoustic evaluation. Although hybrid uid dynamic techniques constitute a good compromise between accuracy and computational cost for engineering purposes, they are still expensive in the case of parametric analysis and shape optimization. To overcome this issue, the development of a Reduced Order Model (ROM) [77, 74, 39, 73] | which alleviates both the computational complexity and data capacity | becomes essential. Moreover, several ROM developments have been realized to account for various eects including subgrid scales [40], heat transfer [34] and stabilization of the Galerkin projection [89]. One of the necessary assumptions to construct an ecient ROM is that, the solution mani- fold of the underlying problem lies in a low dimensional space and, therefore, can be expressed in terms of a linear combination of a few number of global basis functions (reduced basis functions). Among various techniques to generate such reduced basis set, the Proper Orthogonal Decompo- sition (POD) and the Dynamic Mode Decomposition (DMD) have been widely exploited due to their versatile properties [92]. POD provides a set of orthogonal and optimal basis functions [49], whereas DMD computes a set of modes with an intrinsic temporal behavior, hence it is particu- larly suited for time advancing problems [80]. One of the main goals of this work is to understand whether these two established techniques are able to generate global basis functions which include all the wide range of frequencies associated with hydroacoustic phenomena, and therefore allow for accurate noise predictions. In addition, the interest is in assessing the eciency of the modal de- compositions algorithms, evaluating the amount of modal shapes required to accurately reproduce all high frequencies relevant to acoustic analysis. This investigation is of course a rst, fundamental step towards developing ecient reduced order models for hydroacoustic applications. As will be thoroughly discussed, in the present work we only considered data-driven ROMs, but the POD and DMD eciency assessment reported provides valuable information also for the development of projection based ROMs. In literature, POD has been widely used for the past few decades to identify the coherent structures of turbulent ows (see, among others, [88, 3]), and has been applied towards various ow conditions [59, 71]. Correspondingly, DMD has been also exploited [80, 83]. An intuitive question may arise considering the comparative performance. In that regard, several works have carried out both DMD and POD on various ow con gurations. For instance, Liu et al. [103] conducted comprehensive analysis, concluding that DMD has the ability to clearly separate the ow coherent structure in both spatial and spectral senses, whereas POD was contaminated by other uncorrelated structures. Consistent to the previous, in [4] it was noted that DMD is useful when the main interest is to capture the dominant frequency of the phenomenon, while the optimality of the POD modes prevails for coherent structure identi cation that are energetically ranked. In high- speed train DES, Muld et al. [63] examined the convergence and reported that the most dominant DMD mode requires a longer sample time to converge when compared to the POD counterpart. As discussed earlier, parametric studies such as shape optimization can result in extremely high computational costs. Possible ways to circumvent such issue could be degrading the high delity model, or restricting the design parameter space sampling. Since the only analysed parameter in the unsteady uid dynamic problem is the time, in this work we decided to rely on data- driven ROMs. A signi cant source of error in data-driven ROMs is in fact associated with the interpolation of modal coecients based on the problem parameters. Thus, such methods are an extremely valid alternative in presence of a one dimensional parameter space, in which interpolation errors are rather modest. In particular, POD with Interpolation (noted as PODI here-after) can be an adequate solution in such scenarios [54]. The basic idea is to exploit POD on selected ensemble of high delity solutions in the design space to identify the set of optimal basis functions and associated projection coecients representing the solution dynamics. Such nite set of scalar 3 coecients are then utilized to train a response surface that allows predictions at parameter values that are not in the original high- delity ensemble. It was demonstrated that PODI can be eciently utilized in various events: 1) enhancing the temporal resolution of experimental measurements [7], 2) optimal control [54], 3) reconstruction of incomplete data [27], 4) multi-dimensional parametric analysis [102, 30], 5) inverse design [11], 6) variable delity models [60], and more. Both DMD and PODI are regarded as data-driven reduced models [10] since they operate on the snapshots produced from the FOM, and predict the system dynamics in a non-intrusive manner. In shape optimization context, DMD and PODI have been successively applied as in [20, 93, 21, 22, 96, 31] in naval engineering, [78, 26, 37] in automotive engineering, or [70, 43] in aeronautics. For acoustic analysis, only very few studies have been reported in literature. This can possibly include the DMD application as in [9, 44] or POD as in [35, 55, 87]. As previously demonstrated, there are very limited examples of ROMs speci cally tailored for acoustic analysis. To the best of authors' knowledge, the literature is presently devoid of documents which characterize, in a thorough and systematic fashion, the performance of DMD and POD/I techniques on the hydroacoustic ow elds reconstruction or prediction. The overarching goal of this study is to investigate the use of DMD and PODI on a hydroacoustic dataset corresponding to turbulent incompressible ow past sphere at Re = 5000, and computed using wall-resolving LES for the uid dynamic elds, and FWH analogy with direct integration of the nonlinear quadrupole terms for the acoustic elds. In particular, the objectives of this work are 1) to understand the eect of DMD/POD modal truncation on the local and global reconstruction accuracy of the uid dynamic elds, 2) to compare the eciency of DMD and PODI on recovering the ow and spectral information when half the dataset are utilized, 3) to evaluate the performance of DMD and PODI in terms of data compression and dipole/quadrupole acoustic prediction. The present article is organized as follows: rst, a brief overview on the FOM (LES, FWH) and ROM (DMD, PODI) formulations and respective speci c works are presented in section 2. In subsection 3.1, the FOM results are presented and then followed by modal analysis with DMD and POD in subsection 3.2. The reconstructed and predicted uid dynamic data obtained from both ROMs are discussed in sections 3.3 and 3.4, while their spectral and acoustic performances are addressed in section 3.5. Conclusions are drawn in section 4. 2 Methodology First, the high delity data are generated and uniformly sampled in time using high delity sim- ulations with LES turbulence modelling. The resulting matrix of snapshots is then factorized using Singular Value Decomposition (SVD) which is then used to construct both the DMD and POD spaces. The scalar coecients resulting from the projection of the FOM data onto the POD space are used to train a continuous representation of the system temporal dynamics, i.e. PODI approach. In this context, both DMD and PODI are considered as linear approximation of the dy- namical system of the snapshot matrix and, therefore, are used to predict the data at intermediate timesteps. The predicted snapshots are nally exploited to run a post-processing acoustic analogy and validate the accuracy of the noise generation compared to the FOM data. A summary of the procedure is presented in the ow chart in Figure 1. This section is organized as follows: in 2.1 the FOM is introduced recalling also the utilized LES turbulence model, in 2.2 the acoustic model used to perform the hydro-acoustic analysis is introduced, in 2.3 and 2.4 the non-intrusive pipeline used to perform model order reduction using DMD and PODI, respectively, is described. 2.1 Full order model The full order model, adopted to provide the snapshots dataset as input for the reduced order model, is a Large Eddy Simulation. In LES, the large anisotropic and energy-carrying scales of motion are directly resolved through an unsteady three dimensional (3D) simulation, whereas the more isotropic and dissipative small scales of motion are con ned in the sub-grid space. Scale separation is carried out through a ltering operation of the ow variables. In literature, the contribution of the Sub-Grid Scales (SGS) of motion on noise generation and propagation has been found negligible [66, 86]. This means that the LES model can be considered accurate enough to provide a noise-source ow eld, when compared to DNS. At the same time, the unsteady vortex and coherent structures can be of extreme importance when computing the noise signature. 4 Data acquisition High delity uid dynamic FOM in design space. Basis functions Modal decomposition using SVD. Construct ROM DMD/POD + continuous modal coecients. Mid cast Interrogate intermediate snapshots in time as parameter. Acoustic evaluation Analyse acoustic performance from prediction data. Figure 1: Flow chart of the FOM/ROM operations performed in the present study Indeed, it has been shown (see among others [42]) that, in case of complex con gurations (e.g. marine propellers), the RANS methodology may be unsatisfactory in reproducing the ow elds, adopted as input for the acoustic analysis. The detailed numerical simulation of the ow around sphere was described in a previous work [14], where three dierent blu bodies (a sphere, a cube and a prolate spheroid) were in- vestigated concerning their noise signature. In the mentioned work, validations of the LES solver against experimental and DNS data were presented. In particular, the ow elds resulting from the employed framework reproduce the thin boundary layer developing along the wall-normal di- rection, the leading edge, and in the wake. Moreover, they were able to reproduce distributions of the friction coecients, as well as capturing the underlying ow dynamics with a wide coverage of the spectral content. The ltered Navier-Stokes equations in the incompressible regime are considered. Within this sgs setting, the SGS stress tensor = u u u u , which represents the eect of the unresolved i j i j ij uctuations on the resolved motion, is modeled considering the Smagorinsky eddy-viscosity closure: 1 1 @u @u sgs sgs j i = 2 S ; with S = + ; (1) ij t ij ij ij kk 3 2 @x @x i j and is the isotropic second order tensor, while the overbar denotes the ltering operator. The ij SGS eddy viscosity is expressed as: = (C ) jS j; (2) t s ij where = V is the lter width which is de ned as the cubic root of the cell volume V , and c c the Smagorinsky constant C is computed dynamically using the Lagrangian procedure of [58], averaging over the uid-particle Lagrangian trajectories. 2.2 Acoustic model The acoustic model herein considered is the one proposed by Ffowcs Williams and Hawings [29], which is an extension of the Lighthill theory. The basic idea behind the acoustic analogies is that pressure perturbation originates in the ow eld and propagates in the far- eld where the medium is assumed quiescent and uniform. The integral solution of the acoustic wave equation presented by Ffowcs Williams and Hawkings consists of surface and volume integrals, meaning that the sources of uid-dynamic noise can be 5 found as pressure-velocity uctuations developing in the uid region or as re ected pressure on an immersed solid surface. We consider the original formulation presented in [29], and modi ed according to the works of Naja et al. [64] and Cianferra et al. [14]. The modi cation of the original FWH equation takes into account the advection of acoustic waves. To account for the surrounding uid moving at a constant speed (along the x axis), the advective form of the Green's function must be considered. A derivation of the advective FWH equation is reported in [64], where the authors developed an integral solving formulation for the linear (surface) terms. The advective formulation of the volume term for the particular case of the wind tunnel ow is reported in [14]. In the present work, as done for the uid dynamic part, we describe the formulation without dwelling into details, for which we refer to the previous works [14, 15]. The acoustic pressure p ^, at any point x and time t, is represented by the sum of surface (p ^ ) 2D and volume (p ^ ) integrals, respectively: 3D Z Z 1 @ p ~n ^ r ^ p ~n ^ r ^ i i i 4p ^ (x; t) = dS + dS; (3) 2D c @t r r S S 1 @ r ^ r ^ i j 4p ^ (x; t) = T dV 3D ij c @t r 0 f>0 2r ^ r ^ r ^ r ^ R 1 @ i j i j ij + T + dV (4) ij 2 2 c @t r r f>0 3r ^ r ^ R i j ij + T dV: ij f>0 The pressure perturbation with respect to the reference value p is denoted with p ~ = p p , n ^ 0 0 is the (outward) unit normal vector to the surface element dS, and c is the sound speed. r ^ and r ^ are unit radiation vectors, r and r are the module of the radiation vectors r and r respectively. Their description is given in detail in [14]. Equation (4) contains two second{order tensors: R and the Lighthill stress tensor T , the ij ij latter characterizing the FWH quadrupole term. Under the assumption of negligible viscous eects and iso-entropic transformations for the uid in the acoustic eld, the Lighthill tensor reads as: T = u u + p ~ c ~ ; (5) ij 0 i j ij where ~ is the density perturbation of the ow which, in our case, is equal to zero. The surface integrals in Equation (3) are referred to as linear terms of the FWH equation and represent the loading noise term. The volume integrals in Equation (4) are slightly dierent from the standard FWH (non-advective) equation. For their derivation we consider a uniform ow with velocity U along the streamwise direction. Obviously, the direct integration of the volume terms gives accurate results, however, this method can be used if the calculation of the time delays can be omitted, otherwise the computa- tional burden makes it unfeasible. In fact, the calculation of the time delays requires storing at each time step the pressure and velocity data related to the entire (noise-source) volume, in order to perform an interpolation over all the data. However, for the case herein investigated, the evaluation of the non-dimensional Maximum Frequency Parameter (MFP) [15] (which is greater than unity for every microphone considered) allows to adopt the assumption of compact noise source. This means that, in the investigated case, the time delay is very small and the composition of the signals is not expected to contribute to the radiated noise. Since the evaluation of the time delays may be reasonably omitted, a remarkable saving of the CPU time is achieved, and the direct computation of the quadrupole volume terms becomes feasible. In the two previous works [14, 15], as to perform a validation test for the acoustic model, the solution of the advective FWH equation was compared with the pressure signal provided by LES, considered as reference data. This comparison is useful to verify the ability of the acoustic post processing to accurately reconstruct the pressure eld. Also, it points out the frequency range which is important to consider. 6 2.3 Dynamic mode decomposition Dynamic Mode Decomposition is a data-driven modal decomposition technique for analysing the dynamics of nonlinear systems [81, 82]. DMD was developed in [80], and since then it has attracted a lot of attention and lead to diverse applications. Its popularity is also due to its equation-free nature since it relies only on snapshots of the state of the system at given times. DMD is able to identify spatiotemporal coherent structures, that evolve linearly in time, in order to approximate nonlinear time-dependent systems. For a comprehensive overview we refer to [50]. In the last years, many variants of the classical algorithm arose such as, for instance, the multiresolution DMD [51], DMD with control [68], and compressed DMD [28]. Higher order DMD [52] was developed for the cases that show limited spatial complexity but a very large number of involved frequencies. Sparsity promoting DMD to choose dynamically important DMD modes, can be found in [45]. Kernel DMD was proposed in [101], where a large number of observables were able to be incorporated into DMD, exploiting the kernel trick. Randomized DMD for non-intrusive ROMs can be found in [5, 6]. We also cite a new paradigm for data-driven modelling that uses deep learning and DMD for signal-noise decomposition [75]. For numerical non-intrusive pipelines in applied sciences and industrial application we mention [74, 94, 97]. From a practical point of view, the Python package PyDMD [24] contains all the major and most used aforementioned versions stemmed from the classical DMD in a user friendly way, from which the present work also exploits. Here we present a brief overview of the standard algorithm, suppose a given discrete dataset x 2 R , where k 2 [1; :::; m] denotes the snapshot number, and n denotes the number of de- grees of freedom. The m snapshots can be arranged by columns into the data matrix S = x ;x ; : : : ;x and its temporal evolution S = x ;x ; : : : ;x . Now, the objective is to 1 2 m1 2 3 m seek a Koopman-like operator [72] A, which is a best t linear approximation that minimizes y nn y _ _ kS ASk . Such operator is given as A = SS 2 R , where is the Moore-Penrose pseudoin- verse. For practical reasons, the matrix A is not solved explicitly since it contains n elements, from which n is typically of several order of magnitudes, and is usually much larger than the number of snapshots m. Instead, a reduced operator A is considered which has much lower dimension, yet it preserves the spectral information of A. The procedure is de ned in the following. First, a singular value decomposition (SVD) is performed on the snapshots matrix, S = U V U V ; (6) r r where r (m 1) de nes the SVD truncation rank, U is the truncated POD modes as it nr rr will be further discussed in the following section, and the three terms U 2 C , 2 C r r r(m1) and V 2 C denote the truncated SVD components. Second, the reduced operator A is obtained through an r r projection of A onto the POD modes U , that is 1 rr e _ A = U AU = U SV 2 C : (7) r r r r r The eigenspace of A is denoted as (also called the DMD modes), and it can be revealed through the low-rank projection of the eigenvectors W of A onto the POD modes, = U W; (8) where W is obtained from the eigendecomposition of A, i.e. AW = W ; (9) and is the diagonal matrix of eigenvalues of A which is also considered the same matrix obtained from the eigendecomposition of A. Finally, the DMD modes, , and the diagonal matrix of eigenvalues, , are both used to provide a linear approximation to the solution vector x(t) at any time instance, that is x(t) exp(t) x(0): (10) In this work, we are going to characterize the DMD modes for velocity and pressure elds, and exploit them to make predictions between time instants. The reconstructed elds will be used to approximate the hydroacoustic noise. 7 2.4 Proper orthogonal decomposition with interpolation Proper orthogonal decomposition is a linear dimensionality reduction technique that is widely used to identify the underlying structures within large datasets. In general, a linear dimensionality re- duction assumes that each snapshot in the dataset can be generated as a linear combination of a properly chosen small set of basis functions. In this regards, similar to the previous section, we in- nm troduce the snapshots matrix S 2 R which now contains m snapshots, thus S = [x x x ]. 1 2 m nr Then, we assume it exists the linear combination S U C, where the matrix U 2 R contains r r rm the vectors de ning the basis functions, r m, and C 2 R is the matrix containing the linear coecients. This approximation is equivalent to low-rank matrix approximation [56]. If we choose the basis functions that minimize the sum of the squares of the residual (S U C), in the least square sense, then we recover the Proper Orthogonal Decomposition. In such case, the matrix U contains the so-called POD modes, and they can be recovered with the SVD of the snapshots matrix, as de ned in Equation 6, for a low-rank truncation that retain the rst r modes. In fact, the error introduced by the truncation can measured as [69]: 2 2 kS U V k = ; (11) r r r 2 r+1 kS U V k = ; (12) r r F r i i=r+1 where the subscripts 2 and F refer to the Euclidean norm and to the Frobenius norm respectively, and are the singular values arranged in descending order. To approximate the solution manifold, after the projection of the original snapshots onto the POD space spanned by the POD modes, we interpolate the modal coecients for new input parameters. This non-intrusive data-driven approach is called POD with interpolation [54, 11, 12]. A cou- pling with isogeometric analysis can be found in [33]. For an enhancement of the method using active subspaces property, in the case of small computational budget devoted to the oine phase, see [25]. An open source Python implementation of PODI can be found in the EZyRB package [23]. For an implementation based on OpenFOAM [100] see the freely available ITHACA-FV library [91] from which the present work exploits. The PODI method can be decomposed in two distinct phases: an oine phase, and an online one. In the rst phase, given the solutions snapshots for given parameters , we compute the POD modes (which construct the matrix U) through SVD, and the corresponding modal coecients for all the original snapshot with a projection. This results in m r X X x = ; 8k 2 [1; 2; : : : ; m]; (13) k jk j jk j j=1 j=1 where are the elements of the coecient matrix C. jk After performing the oine phase, the r most energetic modes are now obtained, and the corresponding modal coecients are used to construct the reduced space. Using the pairs ( ; ) we are able to reconstruct the function that maps the input parameters to the modal k k coecients. This function is interpolatory for the snapshots used in the training phase, hence the name POD with interpolation. With this surrogate function, we are able to reconstruct in real-time the solution elds of interest for an arbitrary parameter in the so-called online phase. In the present study, the input parameter is the time while the output elds are the velocity and pressure. Also, a cubic spline interpolation is considered for the modal coecients. The degree of the interpolator was rather a modelling choice, given the exibility provided by the simple, one-dimensional coecient samples. In addition, the fairly easy implementation allowed for a robust performance that is comparative to simpler models, while smoother time response is herein preserved due to the higher order polynomial interpolation. Indeed, a more complicated application involving higher dimensional parameter space would require a careful choice of the hypothesis space, and typically more advanced machine learning (or deep learning) techniques would be more relevant, which is not the case in the present study. In literature, several works have been devoted to study the accuracy of the projection coecients manifold construction and performance. In [61] they explored the full factorial and latin hyper-cube 8 Figure 2: Sketch of the computational domain used for the FOM. The sphere diameter is D = 0:01 m. Successive mesh re nement layers (R2, R3, R4) are performed through cell splitting approach until reaching the nest grid spacing 0:001D in the region R5. sampling techniques of the parameter space. The manifold representation was also studied by [46, 61] mainly comparing between linear regression, polynomial, spline, nite dierence, and radial basis functions (RBF). The high-order singular value decomposition (HOSVD) was demonstrated in [53] as an ecient method for ows with strong discontinuities. Kriging interpolation was utilized in [30] for n-dimensional approximation of potentially complex surfaces. 3 Numerical results 3.1 Full order CFD The uid dynamic elds are solved in the framework of the OpenFOAM library [100] which is based on the Finite Volume Method (FVM). The ltered Navier-Stokes equations are solved using the PISO pressure-velocity coupling algorithm implemented in the pisoFoam solver. The spatial derivatives are discretized through second-order central dierences. Implicit time advancement runs according to the Euler scheme. The numerical algorithm, including the SGS closure, has been customized at the laboratory of Industrial and Environmental Fluid Mechanics (IE-Fluids) of the University of Trieste, and more details can be found in [18]. The uid dynamic full order model simulates a sphere of diameter D = 0:01 m, immersed in a water stream with constant streamwise velocity U = 1 m=s. The kinematic viscosity is 6 2 = 2:0 10 m =s, so that the Reynolds number based on the sphere diameter is Re = 5000. The computational domain, depicted in Figure 2, is a box with dimensions 16D 16D 16D along the x, y and z axes respectively. The sphere is located such that a distance of 12D is attained downstream, along the x-axis, while it is centered with respect to the other axes. A zero-gradient condition is set for the pressure at the domain boundaries, except for the outlet where pressure is set to zero. The velocity is set to U at the inlet, stress-free condition is set at the lateral boundaries, and zero-gradient condition is set for the velocity components at the outlet. The grid, unstructured, and body- tted, consists of about 5 millions of cells. It is created using the OpenFOAM snappyHexMesh utility. The grid spacing normal to the wall for the densest layer + + of cells (indicated as R5) is such to have rst cell center within a wall unit y (y = u y= with u = = and the mean shear stress). An A posteriori analysis showed that about 5 grid w 0 w points are placed within 10 wall units o the wall. The grid spacing is obtained through successive transition re nements (indicated as R3 and R4 in Figure 2). A re nement box around the body (named R2 in Figure 2) is considered so as to obtain, in the wake region, a grid size of less than 0:1D at a distance of 8D. Out of the region of interest, a coarser grid (indicated as R1) allows for possible extension of the domain dimensions, and reducing possible disturbance eects coming from the boundaries. 9 0 velocity pressure −1 −2 −3 0 50 100 150 200 Modal index Figure 3: Normalized singular values (SV) for streamwise velocity and pressure snapshots. SV are arranged in descending order. The absence of any gaps in the plots suggests no speci c truncation rank for the SV-based reduced models. For the time integration, a constant time step is set to t = 10 sec in order to keep the Courant number under the threshold of 0:5. The ow around the sphere is completely developed after about 80 characteristic times D=U . 3.2 Modal decomposition In this section, we report the numerical results concerning the modal decomposition of the full order snapshots. This phase is particularly useful in order to have an insight onto the dominant structures and on the frequencies hidden in the full order dynamical system. We report both an analysis on the eigenvalue decay which is associated with the Kolmogorov width, and the modal representation which permits to visualize the turbulent structures associated with each mode. Moreover, we analyse the time evolution of the temporal coecients in order to identify the time frequencies associated with each mode. All the computations have been carried out using the PyDMD Python package [24] and the ITHACA-FV library [91, 90]. 3.2.1 Singular values decay Figure 3 depicts the rst step of the modal analysis applied to the Navier{Stokes uid dynamic problem considered. The plot shows the magnitude of the normalized singular values (SV) obtained from the SVD factorization of the snapshot matrix, obtained both for the streamwise velocity eld component and for the pressure eld. In the diagram | and in all the following discussion | the modes are arranged in descending order according to the corresponding SV magnitude. Typically, a presence of SV magnitude gaps in such plot provides an indication of a convenient truncation rank for the modal analysis. In the present case, a steady and continuous decay is observed after a steep slope corresponding to the rst 15 to 20 modes. Thus, the absence of SV magnitude gaps in the higher frequencies region suggests that, for both uid dynamic elds considered, there is no speci c truncation rank for the modal analysis. 3.2.2 Modal representation As shown, the SV magnitude observation is not resulting in an obvious indication of the modal truncation rank. To start understanding the eect of modal truncation rank on the uid dynamic solution accuracy, we then resort to considerations based on the spatial and time frequencies which need to be reproduced in the hydro-acoustic simulations. Figure 4 depicts a set of three dimensional modal shapes resulting from the longitudinal velocity eld DMD and POD modal decomposition, respectively. The purple diagrams represent isosurfaces passing through the eld data points of the DMD modes at value 15 10 m=s, while the olive color plots refer to isosurfaces of POD modes at value 6:5 m=s. For each modal decomposition methodology, the images in the Figure are arranged in tabular fashion and refer to modes 1, 2, Normalized singular value Velocity (DMD) Velocity (POD) Figure 4: Isosurfaces of the DMD modes (purple color) at value of 15 10 m=s versus the isosurfaces of POD modes (olive color) at value of 6:5 m=s for the streamwise velocity eld. Modes: 1, 2, 4, 8, 36, and 128. 4 on the rst row, and 8, 36, 128 on the second one. The plots suggest that the DMD modal shapes present a more pronounced tendency to be organized according to spatial frequencies with respect to the POD modes. In fact wider turbulent structures are only found in the very rst DMD modes, while in the case of POD modes they can be identi ed also in higher rank modes, among higher frequency patterns. Along with this tendency, the turbulent structures associated with the low rank DMD modes also appear organized in longitudinal streaks, and gradually become more isotropic for higher rank modes. It should be pointed out though that despite the fact that POD modes are in general not designed to separate the contributes of single harmonic components, the modal shapes obtained for the longitudinal velocity do appear to be at least qualitatively correlated to spatial frequencies. In fact, by a qualitative standpoint, the spatial frequencies appearing in the plots corresponding to higher modes are in general higher with respect to those associated to the rst modes. Similar considerations can been drawn from the observation of similar plots corresponding to the modal decomposition of the pressure elds, presented in Figure 5. Also in this case, the purple 4 2 2 plots refer to isosurfaces passing through all data points of value 5 10 m =s of DMD modal 2 2 shape functions, while the olive diagrams refer to isosurfaces of POD modes at the value 300 m =s . The plots are again arranged, for each decomposition methodology considered, in tabular fashion portraying modes 1, 2, 4 on the rst row, and 8, 36, 128 on the second one. The pressure modes associated with both methods seem again qualitatively arranged according to spatial frequency content. As previously observed for the longitudinal velocity modes, the DMD modes appear more closely correlated to spatial frequencies, and present a less isotropic appearance with respect to 11 Pressure (DMD) Pressure (POD) Figure 5: Isosurfaces of the DMD modes (purple color) at value of 5 10 Pa versus the isosurfaces of POD modes (olive color) at value of 300 Pa for the pressure eld. Modes: 1, 2, 4, 8, 36, and 128. 12 −4 × 10 Velocity Velocity (DMD) 4 (POD) −2 −50 N=2 N=4 −4 −100 N=8 N=36 −6 1.0 1.1 1.2 1.3 1.4 1.0 1.1 1.2 1.3 1.4 −4 × 10 Pressure (POD) −25 −50 −75 Pressure −100 −1 (DMD) 1.0 1.1 1.2 1.3 1.4 1.0 1.1 1.2 1.3 1.4 Time [s] Figure 6: Temporal coecients associated to the DMD and POD (left to right) modes number 2, 4, 8, and 36 of the streamwise velocity and pressure data (top to bottom). their POD counterparts. 3.2.3 Associated coecients After having characterized the spatial frequency content of dierent POD and DMD modes, we now want to analyse the time frequencies associated to each mode. To do this, we decompose each snapshot into its modal components and observe the time evolution of the modal coecients. In fact, as the POD and DMD modes generated from the snapshots are constant in time, the corresponding modal coecients must depend on time to allow for the reconstructed solution to reproduce the correct time variation. The four plots presented in Figure 15 show, the temporal evolution of the modal coecients associated, respectively, to the DMD decomposition of the longitudinal velocity eld (top left), to the POD decomposition of the longitudinal velocity eld (top right), to the DMD decomposition of the pressure eld (bottom left) and to the POD decomposition of the pressure eld (bottom right). Each diagram reports four lines referring to the coecients of modes 2, 4, 8, 36. By a qualitative perspective, the diagrams in Figure 15 suggest that DMD and POD modal coecients are strongly correlated with the time frequencies. In fact, higher frequency harmonics appear in the time evolution of higher order modal coecients, which are not observed in lower ones. As expected, also in this case the frequency-mode association is de nitely stronger for DMD modes, in which a single dominant harmonic can be identi ed in correspondence with each modal coecient. As for POD, the respective frequency content seems to cover wider set of harmonics, associated with higher frequencies as higher modes are utilized. 3.3 Flow elds reconstruction A rst aim of this work is to assess whether the proposed DMD and POD algorithms, are able to accurately reproduce the full order model solutions. A rst step in such assessment will be that of checking the eectiveness of the SVD based modal decomposition strategies of the model reduc- tion algorithm considered. In particular, the reconstructed solution convergence to the snapshots considered will be discussed both through the visualization of single snapshots ow elds and by presenting convergence plots of error averaged among snapshots. Finally, we will present similar plots for the solution predicted by means of both DMD and PODI. Temporal coeﬃcients 6 Velocity Velocity 5 (DMD) (POD) 0 0 50 100 150 200 50 100 150 200 r=20 r=120 r=20 r=120 1.25 3 r=40 r=160 r=40 r=160 r=80 r=200 r=80 r=200 1.00 Pressure Pressure (DMD) (POD) 0.75 0.50 0.25 0 0.00 50 100 150 200 50 100 150 200 Snapshot index Figure 7: Relative error percentage in the Frobenius norm for the velocity and pressure (top and bottom, resp.) reconstructed elds using the DMD and POD (left and right, resp.) modes. As expected, reconstruction accuracy converges with higher SVD truncation rank (r). Figure 7 presents a rst evaluation of the eectiveness of the DMD and POD modal decompo- sition algorithm in reducing the number of degrees of freedom of the uid dynamic problem. The results refer to a reconstruction exercise in which the LES solutions at all the time steps have been used for the modal decomposition. The curves in the plots indicate the relative reconstruction error at each time step for both velocity (top plots) and pressure (bottom plots) when a growing number of modes are considered. Such relative reconstruction error is computed as the Frobenius norm of the dierence between the LES solution vector and the reconstructed one, divided by the Frobenius norm of the LES solution. We here remark that to make the velocity and pressure elds error values comparable, the gauge atmospheric pressure value in the simulations has been set to one. A null value would in fact result in lower LES solution norm, leading in turn to large pressure relative errors compared to the velocity ones, even in presence of comparable absolute errors. In Figure 7, the two plots on the left refer to DMD reconstruction results, while the ones on the right present the POD reconstruction error. As can be appreciated, for both modal de- composition methods the errors presented follow the expected behavior, and reduce as a growing number of modes is used in the reconstruction, until machine precision error is obtained when all the 200 modes available are used to reconstruct the 200 snapshots. More interestingly, the data indicate that both for DMD and POD, a number of modes between 80 and 120 leads to velocity and pressure reconstruction errors which fall under 1% across all the time interval considered. It is worth pointing out that, compared to most typical low Reynolds and RANS ows applications of DMD and POD methodologies, such convergence rate is rather slow, as higher number of modes are needed to obtain comparable accuracy. This should not surprise, as LES resolves more tur- bulent structures than the aforementioned models, resulting in higher spatial frequencies which in turn require a higher number of modal shapes to be accurately reproduced. The peak we see for smaller rank truncations around the 40-th snapshot for the DMD reconstruction is justi ed by the fact that those snapshots present a wider range of frequencies thus we need more DMD modes to properly reconstruct them. Finally, the plots suggest that the reconstruction with POD modes leads to errors that are slightly lower to the corresponding DMD errors. In fact, the results consistently show that for both the pressure and the velocity elds, the POD reconstruction error is approximately half of the DMD error obtained with the same amount of modes. b b b b P − P / P % U − U / U % F F F F U-DMD U-POD P-DMD P-POD 40 80 120 160 200 Modal truncation rank Figure 8: Global error of the reconstructed velocity and pressure elds versus SVD truncation rank. The error denotes spatio-temporal averaging of the ow eld data, rst by evaluating relative error in the Frobenius norm, then by averaging over all the snapshots. 3.3.1 Global error Figure 8 presents further veri cation of the modal reconstruction accuracy. The curves in the plot represent the percentage modal reconstruction error | computed in Frobenius norm and averaged among all snapshots | as a function of the number of modes considered in the reconstruction. The black and red curves con rm that the reconstructed velocity and pressure elds, respectively, converge to the corresponding LES elds as the number of modes is gradually increased. Moreover, these results indicate that an ecient reconstruction, characterized for instance by a 1% relative error, would require more than 80 modes for the velocity elds, and more than 30 for the pressure elds. These values are higher than those typically observed for RANS and low Reynolds simu- lations, probably due to the higher spatial frequencies typically found in eddy-resolving solution elds. 3.3.2 Fields visualization The error indicators considered in the previous sections are extremely useful in con rming that the reconstructed solution is globally converging to the LES one as the number of modes is increased. Yet, they oer little information on the error distribution in the ow eld, and the impact of the reconstruction on local ow characteristics of possible interest. In particular, for the test case considered in the present work, it is quite important to assess whether the reconstruction error does not alter the ow in proximity and in the wake of the sphere, as such regions are crucial both to the evaluation of the uid dynamic forces on the sphere and to the acoustic analysis. To this end, in the present section we present a series of visualization of the reconstructed ow elds, which are compared to their LES counterparts. 2 2 Making use of the Q-criterion, de ned as Q = 0:5(k k kSk ) with and S denoting the vorticity and strain rate tensors respectively, Figure 9 depicts the turbulent structures characteriz- ing the ow in the wake region past the sphere. By de nition, a positive value of Q implies relative 4 2 dominance of the vorticity magnitude over the strain rate [38]. Here, a positive value of 10 1=s is chosen to generate isosurfaces which pass through all data points holding this value. In the gure, the top plot refers to the original LES solution obtained at the last snapshot of the dataset, while the centered and bottom plots refer to the corresponding DMD and POD reconstructions, respec- tively, utilizing 160 modes. The images show that both POD and DMD reconstruction algorithms lead to fairly accurate representation of the turbulent structures shape past the sphere. In fact, the con guration of the wider vortical structures detaching from the sphere appears to be correctly reproduced in the reconstructed solution. As for ner details associated with the smaller turbulent scales, the DMD reconstruction is observed to be in closer agreement with the original LES solution than the POD reconstructed eld. Such observation is consistent with ndings from [103, 4] in which they demonstrated the superiority of DMD to accurately determine spectral and convective Global error [%] Figure 9: Coherent structures represented by the iso-contours of the Q-criterion at non-dimensional 2 2 value QD =U = 4 for the last snapshot (top), compared with corresponding ones obtained from modal reconstruction with DMD (centered) and POD (bottom) using 160 modes. information of the vortical structures in wake regions. For a more signi cant and quantitative assessment, Figure 10 includes a series of of contour plots representing the instantaneous ow eld at time instance corresponding to a maximized relative error, cf. Figure 7. The three plots in the rst row represent contours of the LES streamwise velocity component eld, and of its DMD and POD reconstructed counterparts, respectively, utilizing 160 modes. At a rst glance, the reconstructed elds seem to reproduce the main features of the LES ow. In particular, both the stagnation region ahead of the sphere and the ow detachment past it appear to be correctly reproduced by both modal reconstruction strategies. In addition, the detached vortex, located downstream with respect to the sphere in this particular time instant, is also correctly reproduced. For a better assessment, the two images of the second row represent contours of the local error for the DMD and POD streamwise velocity reconstruction, respectively. Both for the DMD and POD reconstruction, higher local error values are found in the wake region downstream with respect to the sphere. In particular, it observed that the local error peaks in the DMD reconstruction is larger than that for POD. Additionally, the high frequency error pattern and the elevated local error values located in the wake region seem to indicate that, as expected, the modes disregarded in the reconstruction are associated with high spatial frequencies. A similar comparison is presented for the pressure eld in the following rows of the gure. The three plots in the third row represent the instantaneous pressure eld obtained with LES, and its reconstructions computed with DMD and POD, respectively, for the same snapshot. Again, the full order eld appears well reproduced by both DMD and POD reconstruction algorithms, as features like the peak pressure in the stagnation region and the pressure minimum within the vortex detaching past the sphere are correctly reproduced. The plots of pressure reconstruction error for DMD and POD, respectively, presented in the last row con rm that both methods are capable to adequately represent the LES solution. Also here, the disregarded high spatial frequency modes are likely responsible for the high frequency error pattern observed. 3.3.3 Error statistics in the wake region As Figure 10 shows, the highest reconstruction errors in the velocity and pressure elds reconstruc- tion are mostly located in the wake region. Therefore, a follow up analysis is required to quantify the spatial error distribution within such region and to assess the local convergence behavior of both the DMD and POD modal decomposition methodologies. Figure 11 presents a normalized density function plot based on the reconstruction error of both the velocity and pressure elds corresponding to the snapshot illustrated in Figure 10. To generate the plot, the following steps are performed. First, the computational cells in the wake region are identi ed according to the con- dition k! k = kr Uk > 1:0 with ! denoting the streamwise vorticity component. Second, the x x x arrays corresponding to the velocity and pressure error elds in the wake region are computed, i.e. U P ^ ^ e = (U U ) and e = (P P ) wherein the subscript denotes the wake cells. Third, to account w w w w w w for variations in the grid resolution, the mentioned error elds within the wake region are weighted P P n n norm norm norm by the corresponding normalized cell volumes V , i.e. e = ( V e = V ) w w;i c c;i c;i i=1 i=1 wherein n is the cell count in the wake. Finally, the error interval for both the velocity and 16 Figure 10: Left to right: LES, DMD and POD reconstructions with N = 160 modes. Top to bottom: streamwise velocity, corresponding error elds, pressure, corresponding error elds. In- stantaneous snapshot is selected based on errors peak. 17 50 25 Velocity Velocity r=20 r=20 (DMD) (POD) r=40 r=40 40 20 r=80 r=80 r=120 r=120 30 15 r=160 r=160 20 10 10 5 0 0 −0.15 −0.10 −0.05 0.00 0.05 0.10 0.15 −0.15 −0.10 −0.05 0.00 0.05 0.10 0.15 Pressure Pressure 50 50 (DMD) (POD) 40 40 30 30 20 20 10 10 0 0 −0.05 0.00 0.05 −0.05 0.00 0.05 error [-] Figure 11: Error statistics of velocity and pressure eld data of selected snapshot corresponding to maximized relative error. Sampled data points are conditioned by k! k = kr Uk > 1:0 x x to identify the wake region. Resulting density function is weighted by normalized cell volumes. Presented plots correspond to various modal truncation ranks (r). pressure data are uniformly divided into equal-width bins and the density function (normalized histogram) is plotted for each bin. The scattered plots highlight dierences in the behavior of the two modal reconstruction strate- gies considered. The diagrams on the left, which refer to the DMD results for velocity (top) and pressure (bottom), show in fact a clear error reduction as the modal truncation order is increased. The curves corresponding to growing truncation orders tend to get closer to the vertical axis as the error, displayed on the horizontal axis, is progressively reduced. The same convergence rate cannot be observed in the POD plots on the right, as both the velocity (top) and pressure (bot- tom) reconstruction error statistical distribution curves appear less aected by an increase of the modal truncation order. Again, this observation can be explained in the light of the DMD theory and the ability of its modes to be organized according to the eld spatial frequencies, unlike POD ones which can instead become contaminated by uncorrelated structures, according to the claim reported in [103]. It is worth pointing out that for image de nition purposes, the left and right tails of Figure 11 have not been reported. Yet, a cross comparison with Figure 10 readily suggests that POD results showed lower error margin | hence narrower tails | in this regard. Therefore, we could infer that higher DMD modes are capable to capture more frequencies in the wake, resulting in lower mean error but higher peak errors compared to POD modes. 3.4 Fields mid cast using DMD and PODI The previous sections are focused on assessing the accuracy of the POD and DMD modal decom- position strategies. The reconstruction results con rmed that both methods can be considered eective tools for the reduction of the degrees of freedom of the uid dynamic problem. We now want to analyse the ability of the data-driven DMD and PODI reduced order models considered in this work, in predicting the LES solution at time steps that are not included in the original snapshot set. In particular, throughout the remaining analysis of this work, the original 200 LES snapshots are decomposed into two sets, one set comprises the 100 odd snapshots from which they are used to train the ROMs and are called the train dataset hereafter, while the remaining even snapshots within the temporal interval of the train dataset are contained in a test dataset, and are used to compare the full order solution with the DMD and PODI model prediction results (i.e. Normalized density function Table 1: ROM performance at various modal truncation rank. Data compression level is de ned as (GB /GB ) with GB =74. Speedup is de ned as (CPU =CPU ) with CPU = FOM ROM FOM FOM ROM FOM 2:4 10 sec. Compression level is averaged between DMD and PODI. Rank POD cumulative energy Speedup Compress. (r) U P DMD PODI level 10 0.999992 0.994093 29.29 35.67 19.995 20 0.999996 0.997076 28.46 32.46 9.9973 40 0.999998 0.998735 26.41 27.93 4.9986 60 0.999999 0.999295 25.41 24.68 3.3324 80 0.999999 0.999571 23.58 22.48 2.4993 100 1.000000 1.000000 23.18 20.46 1.9995 ROM prediction dataset). Before progressing with analysis, a summary on the POD energetic content, as well as the data compression level and computational speedup for the considered DMD and PODI data-driven models at various modal truncation ranks is listed in Table 1. Here, the kinetic energy content within a ROM is described by cumulative sum of the POD eigenvalues. The compression level, taken as the arithmetic mean between DMD and PODI in the table, considers the size ratio of the FOM data (train and test datasets, totalling 74 GB) to the ROM data (spatial modes of the train dataset, and continuous representation of the temporal dynamics computed via time integration or cubic spline interpolation for DMD or PODI, respectively). Speedup is de ned as (CPU =CPU ) with CPU = 2:4 10 sec corresponding to the time required to solve FOM ROM FOM for the 200 snapshots and to write out the train dataset, while CPU is the time needed to ROM 1) perform modal decomposition on the train dataset, 2) extract the associated dynamics with a continuous representation, and 3) write out the prediction dataset in OpenFOAM format. It is worth noting that, since time is the considered parameter in the present ROM procedure, an oine phase still requires solving for the same temporal window as in the FOM solution, in order to obtain the train dataset. Nevertheless, a computational gain in the generated ROM can be still attained, since the less amount of stored data (i.e. few modes and associated dynamics) compared with FOM are utilized for a swift construction of the uid dynamic elds at arbitrary sample points. Such procedure is considered a lot faster and more economic than recomputing the simulation to write out the uid dynamic eld data at those sample points. Indeed, an extension of the present procedure to consider additional parameters while performing, for instance, multiparameter interpolation using PODI, would result in a further speedup since the oine phase would solve a full order solution only for a subset of the parameter combinations in the parameter space, hence saving up complete FOM computations from being performed. Such multiparameter investigation is considered a follow up study of this work. To this end, now we progress the analysis by considering the temporal evolution of the error due to ROMs prediction. 3.4.1 Prediction error analysis The rst test presented is designed to assess the predictive accuracy of the DMD and PODI methodologies described in this work. Figure 12 depicts the results obtained with the ROM methodologies applied in such alternate snapshots arrangements as previously described. The top plots report the percentage Frobenius norm error between the DMD and PODI predictions with respect to the LES solution for the longitudinal component of the velocity, while the bottom diagrams present similar error norms corresponding to the pressure eld. The left plots refer to the DMD results, while the right ones depict the PODI errors. We point out that, the | signi cantly lower | error values corresponding to the train dataset have been omitted for clarity. Finally, the dierent colors in the plots indicate growing number of modes considered in the DMD and PODI prediction, up to a maximum of 100. The results show a substantial convergence of the data-driven reduced model solutions to the full order one. In particular, selecting a number of modes between 80 and 90 results in errors lower than 2% on the velocity eld for most time steps in both DMD and PODI methods. As for the pressure eld, the bottom left diagram in Figure 12 shows that 1% error goal can be obtained with an even slightly lower amount of DMD modes. It must be pointed out though, that the DMD solution 19 Velocity Velocity (DMD) (PODI) 0 0 0 50 100 150 200 0 50 100 150 200 r=10 r=60 r=10 r=60 1.5 r=20 r=80 r=20 r=80 r=40 r=100 r=40 r=100 Pressure Pressure 1.0 (DMD) (PODI) 0.5 0 0.0 50 100 150 200 50 100 150 200 Snapshot index Figure 12: Relative error percentage in the Frobenius norm for velocity and pressure (top and bottom, resp.) elds, where half sample rate snapshots are used to train the reduced model, DMD (left) or PODI (right), for particular truncation (r) while predicting the intermediate snapshots. Plots show only the snapshot-wise prediction error, while disregarding errors in the training set, which is almost null. in the very rst time steps is not converging to the LES one, and the error grows as the number of modes is increased. This behavior, which could be related to the high frequencies introduced by the higher DMD modes added to the solution, is currently under further investigation. Aside from the rst few time steps, it is generally observed that PODI errors are consistently lower by a factor two with respect to the DMD ones, especially in the cases utilizing fewer modes. It is worth pointing out that a direct comparison between the PODI plots in Figure 8 and those in Figure 12 can indirectly result in a possible estimate of the interpolation error associated with the PODI strategy. The plots suggest that the eect of interpolation on the global spatial error at each time step is rather low, as the errors in Figure 12 present the same behavior and are not signi cantly higher than the ones obtained with pure reconstruction. Now, to summarize the performance of each ROM with respect to the modal truncation level, a global spatio-temporal error is measured against growing number of modes, as depicted in Fig- ure 13. In particular, similar to the previous analysis in Figure 8, the mean absolute error is eval- uated for the percentage Frobenius norm spatial error for all the prediction dataset, cf. Figure 12. Here, it is observed that, for both the velocity and pressure predictions, a signi cant reduction in the global error is achieved in the PODI models with respect to DMD up to a utilization of 60 modes, after which the global error becomes comparable between both ROMs. Additionally, a nite global error is noted in the gure even when a full modal rank is utilized. Indeed, such observation should not be surprising since half the dataset is only employed to train both ROMs while interrogating the remaining sample points. This is not the case in pure eld reconstructions, cf. Figure 8, where the global error vanishes with a full modal rank since dynamic and spectral information become entirely recovered. 3.4.2 Coherent structures The prediction error indicators considered in the previous section indicate whether the reduced models solution is globally converging to the LES one as the number of modes is increased. We now resort to ow visualizations to obtain better information on the error distribution in the ow eld and on the reduced models performance in reproducing local ow characteristics of possible b b b b P − P / P % U − U / U % F F F F 4 U-DMD U-PODI P-DMD P-PODI 20 40 60 80 100 Modal truncation rank (r) Figure 13: Global error of the prediction dataset versus SVD truncation rank. The error denotes spatio-temporal averaging of the ow eld data, rst by evaluating relative error in the Frobenius norm, then by averaging over all the snapshots. Figure 14: Coherent structures of the last snapshot in the prediction dataset, represented by the iso- 2 2 contours of the Q-criterion at non-dimensional value of QD =U = 4, compared with corresponding ones obtained from DMD-100 (centered), and PODI-100 (bottom). interest. 2 2 Again, making use of the Q-criterion isosurfaces at value QD =U = 4, Figure 14 portrays the turbulent structures characterizing the wake ow past the sphere. The top plot refers to the original LES solution obtained at the last snapshot of the prediction dataset, while the centered and bottom plots refer to the respective DMD and PODI predicted solution utilizing 100 modes. Also in this case, the images show that both DMD and PODI reduced order models allow for rather accurate reproduction of the turbulent structures past the sphere. For both ROM solutions, the main vortical structures detaching from the sphere appear in fact very similar to those of the original LES ow eld. Finer details associated with smaller turbulent scales are also in good agreement, hence suggesting that errors in the PODI time interpolation and DMD time integration are not signi cantly higher with respect to the reconstruction error analysed earlier. 3.4.3 Drag and lift coecients Besides coherent structures, it is important to assess the performance of data-driven ROMs in cap- turing hydrodynamic phenomena that could nd particular interest in the engineering community. One of these phenomena is the drag and lift forces on the sphere surface. First, pressure-induced drag and lift forces are evaluated by integrating the pressure in the streamwise and the ow-normal directions, respectively, over the sphere surface. Then, the pressure forces are normalized by the dynamic pressure acting on the sphere, with a sphere projected cross-sectional area being con- sidered for the reference area, to compute coecients of drag (C ) and lift (C ). The temporal D L evolution of both coecients on the sphere surface is depicted in Figure 15 for various ROM pre- dictions compared with the FOM. The presented plots show that both C and C are accurately D L Global error [%] (DMD) (PODI) 0.5 0.5 LES r=40 r=20 r=80 0.4 0.4 1.0 1.1 1.2 1.3 1.4 1.0 1.1 1.2 1.3 1.4 (DMD) (PODI) 0.1 0.1 0.0 0.0 1.0 1.1 1.2 1.3 1.4 1.0 1.1 1.2 1.3 1.4 Time [s] Figure 15: Temporal evolution of the drag and lift coecients corresponding to pressure forces on the immersed sphere surface. The force coecients are computed for various DMD and PODI models in comparison to FOM data. Presented plots show the predictive performance of both data- driven ROMs on estimating the sphere pressure forces, with noted superiority of PODI predictions over DMD. predicted while only half the dataset is employed. It is noted that the accuracy of C predictions is higher in comparison to C predictions, and that employing as low as 40 modes are able to pre- dict the forces coecients with a sucient accuracy. More importantly, the predictive character of PODI models are observed to be superior over DMD models for a particular modal truncation level. The observed discrepancies in the early temporal window of C predictions via DMD at r = 80 are related to the poor DMD predictions of the ow elds at early snapshots, as previously noted in the relative error plots, cf. Figure 12. On a general note, the presented data-driven ROMs are able to predict the pressure dynamics on the surface of the immersed sphere while employing considerably lower levels for modal truncation. We would like to highlight that the framework that we are proposing is solely based on input and output quantities. Therefore, instead of computing engineering quantities from the reconstructed ow elds, one could also construct a reduced order model to directly approximate them. In such a case we would expect that the so generated ROM would produce even more accurate results. 3.5 Spectral analysis 3.5.1 Data probes inside and outside the wake region A rst necessary step in order to evaluate how the frequency content of the FOM solution is reproduced at the ROM level, consists in observing the time dependency of the solution at xed points in the computational domain. To this end, Figure 16 presents the reduced models results for the time evolution of the streamwise velocity component in correspondence with two dierent locations of the ow eld. In the 2 2 tabular arrangement of the gure, the top plots correspond to a point located outside of the sphere turbulent wake, while the bottom plots refer to a point inside the turbulent wake. In addition, the two diagrams on the left refer to the DMD results, and the ones on the right report the PODI result. In the plots, the red lled circles refer to the LES result at prediction sample points, while the dash-dotted lines represent the ROM results, which have been drawn for a growing number of modes. Finally, the black continuous line represents the pure DMD and PODI reconstruction result utilizing 160 modes. C C L D 1.0175 LES r=160 (recon) 1.0150 r=20 r=80 1.0125 1.0100 1.0075 out-wake out-wake (DMD) (PODI) 1.0050 1.0 1.1 1.2 1.3 1.4 1.0 1.1 1.2 1.3 1.4 0.75 in-wake in-wake (DMD) (PODI) 0.50 0.25 0.00 −0.25 −0.50 1.0 1.1 1.2 1.3 1.4 1.0 1.1 1.2 1.3 1.4 Time [s] Figure 16: Time history of streamwise velocity signals at xed points located outside the wake at x = f2D; 0; 2Dg (top) and inside the wake at x = f2D; 0; 0g (bottom) for DMD (left) and A B PODI (right) models. Red symbols represent LES results, while dashed lines represent the ROM results for dierent SVD truncation ranks (r). Black continuous line represents the reconstructed eld utilizing 160 modes. As expected, the plots show that the velocity variation over time becomes more chaotic in the wake region, due to the high frequency uctuations associated with the turbulent structures typically observed in such part of the ow domain. All the plots show that, as the number of DMD and PODI modes used is increased, the reduced solution approaches the full order model one. At a rst glance, the PODI convergence at the point located outside of the wake appears faster than that obtained with DMD. The reduced PODI solution obtained with 20 modes is, in fact, already rather close to the original time signal, while the corresponding DMD solution is still far from the LES one. The behavior of both DMD and PODI reduced models is clearly more accurate when a higher number of modes is selected. The 80 modes curve obtained with both ROMs is practically indistinguishable from the LES one. Given the small estimate of the time interpolation error, as discussed in subsubsection 3.4.1, the faster PODI convergence observed should depend on the fact that single PODI modes are richer in spatial frequency content than the DMD ones. The low order PODI modes might then already include higher frequencies not contained in the corresponding DMD modes. Such spatial frequency content should nally reverberate in the time evolution of the solution including some higher frequencies also when the only lower order PODI modes are used. Despite these favorable characteristics, when the full order model solution presents even higher frequencies, also the PODI solution requires a higher number of modes to obtain satisfactory accuracy. This is clear by the plots of the time history of the velocity at the point within the wake region. Here, both PODI and DMD reduced result curves become suciently close to the original solution one only when 80 modes are considered. This should not come as a surprise, as the in-wake signal presents higher frequency uctuations due to the wake turbulent structures, which can be fully represented only by including high order modes. The reduced pressure local time evolution prediction results presented in Figure 17 exhibit the same behavior observed for the reduced velocity eld. Also this Figure, in which the plots and the curve colors are arranged as described in Figure 16, suggests in fact that both inside and outside the wake region, employing 80 modes allows for both ROMs to obtain pressure predictions that are suciently close to the corresponding LES time signals. Again, the plots also suggest that PODI in the region outside the wake is able to obtain viable pressure predictions with fewer modes. Thus, these results indicate that in presence of fully attached ows, in which wake eects are less dominant, the choice of POD as modal decomposition methodology could result in more economic Velocity [m/s] out-wake out-wake (DMD) (PODI) −0.004 −0.006 LES r=160 (recon) −0.008 r=20 r=80 −0.010 1.0 1.1 1.2 1.3 1.4 1.0 1.1 1.2 1.3 1.4 0.3 in-wake in-wake 0.2 (DMD) (PODI) 0.1 0.0 −0.1 −0.2 −0.3 1.0 1.1 1.2 1.3 1.4 1.0 1.1 1.2 1.3 1.4 Time [s] Figure 17: Time history of pressure signals at xed points located outside the wake at x = f2D; 0; 2Dg (top) and inside the wake at x = f2D; 0; 0g (bottom) for DMD (left) and PODI (right) models. Red symbols represent LES results, while dashed lines represent the ROM results for dierent SVD truncation ranks (r). Black continuous line represents the reconstructed eld utilizing 160 modes. reduced models. 3.5.2 Fast Fourier transform The plots in Figure 16 and Figure 17 seem to indicate that, on a qualitative level, the PODI and DMD solutions are able to recover a growing portion of the full order model frequency content as the number of modes is gradually increased. This aspect should be, of course, investigated in a more accurate way, as the presence of high frequency components in the ROM solution results in their ability to be eective surrogates in acoustic analysis. Thus, we make use of the FFT spectra with the aim of obtaining a quantitative assessment of the impact of the POD and DMD modes considered on the frequency content of the solution. Figure 18 presents the magnitude of the FFT of the local streamwise velocity signals previously presented. Also in this case, the top plot refers to the point located outside of the wake, while the bottom plot refers to the in- wake point spectrum. The black, continuous lines refer to the FFT of the local velocity signal obtained with the LES full order model. The red dash-dotted lines represent the spectrum of the signal composed only by the train dataset which, of course, is characterized by half the sampling frequency with respect to the full LES signal. Finally, the green and yellow dashed lines respectively denote the corresponding plots obtained by means of truncated DMD (left panel) and PODI (right panel) models employing 20 and 80 modes each. Figure 19 | which also employs the line color arrangement just described | displays analogous spectral results obtained when the pressure eld is considered. We point out that, based on experimental results presented in [76], a Re = 5000 ow past a sphere results in the Strouhal number S = 0:2 associated with vortex shedding. This corresponds to f = S U=D = 20 Hz being U = 1:0 m=s and D = 0:01 m. In Figure 18 and Figure 19, it is clearly visible that a 20 Hz peak appears in the LES signal spectra when the velocity and pressure probe is located within the cylinder wake (lower plots). Such peak is accurately reproduced by both PODI and DMD results. The top plots in both Figures 18 and 19 show that all the reduced models considered allow for a suciently good reconstruction of the solution spectra for a point outside of the wake. In fact, all the yellow dashed curves appear very close to the black continuous one representing the LES Pressure [N/m ] −3 LES train r=20 −4 r=80 −5 out-wake out-wake (DMD) (PODI) −6 0 50 100 150 200 250 0 50 100 150 200 250 −1 −2 −3 in-wake in-wake (DMD) (PODI) −4 0 50 100 150 200 250 0 50 100 150 200 250 Frequency [Hz] Figure 18: Normalized amplitude spectrum of the streamwise velocity component, computed as (jt(U U )j=N ) with U denoting mean velocity, in a xed point located outside the wake (top) freq and inside the wake (bottom) for DMD (left) and PODI (right) at dierent truncation ranks (r). Black continuous and red dash-dotted lines refer to LES results of the full dataset and train dataset, respectively. Green and yellow dashed lines mark ROM surrogates employing 20 and 80 modes, respectively. −3 LES train −4 r=20 r=80 −5 −6 out-wake out-wake (DMD) (PODI) 0 50 100 150 200 250 0 50 100 150 200 250 −2 −3 −4 in-wake in-wake (DMD) (PODI) 0 50 100 150 200 250 0 50 100 150 200 250 Frequency [Hz] Figure 19: Normalized amplitude spectrum of the pressure signal, computed as (jt(p p)j=N ) freq with p denoting the mean pressure, in a xed point located outside the wake (top) and inside the wake (bottom) for DMD (left) and PODI (right) at dierent truncation ranks (r). Color arrangements are same as in Figure 18. Normalized amplitude spectrum of pressure Normalized amplitude spectrum of velocity solution spectrum. Thus, this con rms that considering 80 modes or more leads to a spectrum that is indistinguishable from the original, especially for the PODI results in which a relative superiority is again noted in comparison with the spectra from DMD models. Yet, the higher frequency turbulent structures occur in the wake. In such region, as suggested by the bottom plots in Figures 18 and 19, the frequency content of the streamwise velocity and the pressure solutions is de nitely richer, and the accuracy of the reduced PODI and DMD solutions is clearly lower when higher frequencies are considered. Here, for velocity signals, the plot suggests that considering 80 modes both PODI and DMD algorithms lead to good quantitative spectral reconstructions of the velocity signal for frequencies up to approximately 110 Hz. Interestingly, being based on the train signal yellow dashed line, in this case the PODI and DMD algorithms seem able to improve the behavior of the signal interpolated at half the sampling frequency, and somewhat extend its accuracy at frequencies very close to the Nyquist one. It is worth commenting on the ROM sensitivity against grid resolution. When a coarser grid is utilized, larger length scales are modeled. As it is expected to note deviations in the corresponding FOM performance with respect to DNS data, it is also doubtful to expect recovering these missing details via a ROM. Furthermore, ROM predictions closely correlate with the energetic content and the induced frequencies from the employed dataset and respective modes. This implies that, while shifting from LES towards RANS grid resolutions, the underlying eld dynamics are expected to be suciently recovered with less snapshots and less modes for the ROM. 3.6 Acoustic analysis The acoustic analysis carried out in this work is aimed at comparing FWH signal obtained using the full order model data and the corresponding signal obtained with both the PODI and the DMD reduced order models. By a practical standpoint, the FWH post processing described in subsec- tion 2.2 is applied to the pressure and velocity elds obtained from the LES full order simulation, and from both the PODI and DMD reduced models. The FWH integrals are here computed consid- ering two microphones located at x = f0; 2D; 0g (microphone A) and x = f2D; 2D; 0g mic A mic B (microphone B). As mentioned, the sphere is centered at the origin O = f0; 0; 0g and has diameter D. The acoustic pressure time-history is converted to sound Spectrum Level SpL = 20 log(p=p ), ref considering p = 1 Pa. ref We will rst focus on the results of the PODI model. The corresponding SpL from both microphones is reported in Figures 20 and 21, right panels. In the gures, the dierent FWH signal curves denote the dierent number of modes employed in the PODI reconstruction. Also, as seen in previous works where the FWH formulation is used, it is convenient to separate the linear contribution to the acoustic pressure, obtained from the surface integrals in Equation 3, from the nonlinear part obtained from the volume integrals in Equation 4. Such contributions are referred to as dipole and quadrupole terms, respectively. This distinction is particularly relevant, as the former contribution only requires the knowledge of the pressure eld on the body surface, while the latter component takes into consideration the evolution of the pressure and velocity eld in the wake region. More speci cally, we emphasize that to obtain an adequate reconstruction of the non-linear acoustic signal, an accurate reconstruction of the vorticity eld is necessary. In general, the PODI results seem satisfactory especially for what concerns utilizing 80 modes. Figures 20 and 21 (right top plots) show in fact good agreement between the linear acoustic pressure contribution based on the LES pressure eld and the corresponding linear contribution computed with PODI employing 80 modes (yellow lines in the plots). As the plots suggest, this trend is appreciable both at microphones A and B. As for the nonlinear contribution, Figures 20 and 21 (right bottom panels) also show a satisfactory agreement between the acoustic pressure signals obtained from the LES ow elds and those computed based on PODI. In this case, in the plots referring to both microphones, the blue and yellow lines (40 and 80 modes respectively) on the right bottom panels are suciently close to their reference FOM counterpart, represented by the continuous black line. At a closer look, close agreement between PODI and FOM acoustic pressures, both for the linear and non-linear contributions, can be established at least up to frequencies in the range 80{ 100 Hz, which may be considered the meaningful frequency interval related to the considered test case. In fact, we point out that, since the reference pressure is set to 1 Pa, the ambient noise may be considered in the range of 60{100 dB. At higher frequencies though, we note that some spurious 26 LES r=10 r=40 r=80 dipole dipole (DMD) (PODI) 1 2 1 2 10 10 10 10 quadrupole quadrupole (DMD) (PODI) 1 2 1 2 10 10 10 10 Frequency [Hz] Figure 20: Linear dipole (top) and nonlinear quadrupole (bottom) terms of FWH equation evalu- ated from LES data and compared to corresponding DMD (left) and PODI (right) reduced models at dierent truncation ranks (r). Microphone A. LES r=10 r=40 r=80 dipole dipole (DMD) (PODI) 1 2 1 2 10 10 10 10 quadrupole quadrupole (DMD) (PODI) 1 2 1 2 10 10 10 10 Frequency [Hz] Figure 21: Linear dipole (top) and nonlinear quadrupole (bottom) terms of FWH equation evalu- ated from LES data and compared to corresponding DMD (left) and PODI (right) reduced models at dierent truncation ranks (r). Microphone B. Spectrum level for microphone B [dB] Spectrum level for microphone A [dB] oscillations appear in the PODI results, especially in the nonlinear noise terms corresponding to microphone B (Figure 21), at which the value of the integrals in the FWH formulation are strongly conditioned by the presence of the wake. The observed spurious oscillations could be due to an imperfect reconstruction of the vorticity eld, although it must be emphasized that the sampling frequency of the original dataset is not very high, f = 500 Hz, thereby for the train dataset it becomes f = 250 Hz for such an alternative arrangement. Therefore, it is reasonable to obtain less accurate results at frequencies higher than 125 Hz. Further work will be devoted to investigating whether longer temporal records or higher sampling frequencies will eliminate or mitigate the problem. However, considering the 80 modes PODI reconstruction, the maximum error is of the order of 5 dB, observed at microphone B, for the non linear terms (Figure 21 bottom right panel). In addition, the main peaks observed at very low frequencies, up to 20 Hz, are well captured by the PODI signal. In such case, the agreement is also veri ed for the 40 modes signals. As these low frequency peaks describe the main and most energetic features of the acoustic signal, this indicates that a moderate amount of modes might still be a good compromise when the acoustic analysis is only aiming at a general characterization of the principal noise features. Finally, we must point out that the errors observed in the non-linear noise contributions are not in general higher than those introduced in the linear contributions, as might be expected given the interpolation procedure involved in PODI. In some occurrence (Figure 20 top right plot), the linear contribution errors are even quite surprisingly higher than their nonlinear counterparts. This might be related to the fact that the linear noise contributions are only based on pressure evaluations on the body surface, while the non-linear contributions are based on volume integrals. The POD procedure used selects the modal shapes so as to minimize the error in a norm based on the whole volumetric solution, rather than only on a surface restriction. For such reason the modal shape selected might result in non optimal results in the computation of surface integrals. Possible gains might then be obtained, if needed, adding a separate PODI only built on the body surface degrees of freedom, which would result in a modest increase of the computational cost. As regards the acoustic signals provided by the DMD method, a good agreement is observed with respect to the FOM spectrum, up to about 80 Hz, for both microphones, and both linear and non-linear terms, see Figures 20 and 21 (left panels). The most accurate signals are obtained using 40 and 80 modes, while for the signal related to 10 modes we observe, as expected, a considerable discrepancy. In general, the reconstruction of the vorticity eld, obtained by both POD and DMD, is found to provide an adequate input eld for the acoustic model. Particularly, the volume integral of the FWH equation involves both the velocity and the pressure eld in the wake region. Having observed a good match of the ROM signals with respect to the FOM reference signals, we may conclude that the entire spectrum associated with the vortex wake has been adequately reconstructed. 4 Conclusions and perspectives This work discussed the application of data-driven dimensionality reduction algorithms on a hy- droacoustic dataset which was numerically measured using Large Eddy Simulation (LES) uid dynamic turbulence model and the Ffowcs Williams and Hawkings (FWH) acoustic analogy. An extensive set of data was presented to fully characterize the ability of both Dynamic Mode De- composition (DMD) and Proper Orthogonal Decomposition (POD) to reconstruct the ow elds, based on their spectral and energetic contents, with all spatial and temporal frequencies needed to support accurate noise predictions. First, Singular Value Decomposition (SVD) analysis did not indicate the presence of signi cant constraints on the modal truncation rank for such ows. In fact, no signi cant singular value energy gaps were observed. SVD was then used to extract DMD and POD modes with associated coecients, and to employ them for ow elds reconstruction. In general, both DMD and POD algorithms showed ecient reconstruction accuracy. Spatial and temporal error analyses indicated relatively lower error magnitudes in POD-based reconstructed elds. On the other hand, statistical analysis and vortical structures identi cation methods demonstrated the ability of DMD to capture nest wake scales by employing higher modes, compared with POD. Second, two data-driven reduced models based on DMD mid cast and on POD with interpo- lation (PODI) were created utilizing half the LES original dataset. Both DMD and POD based reduced models showed good eciency and accurate ow reconstruction. In addition, the spectral 28 − 3 LES (P signal) − 4 − 2 − 5 − 3 − 6 out-wake in-wake − 4 0 50 100 150 200 250 0 50 100 150 200 250 Frequency [Hz] LES (dipole) LES (quadrupole) mic-A mic-B 1 2 1 2 10 10 10 10 Frequency [Hz] Figure 22: Top: FFT spectrum of the pressure signal at points located out of wake (left) and inside the wake (right). Bottom: dipole and quadrupole acoustic noise spectrum measured at mic-A (left) and mic-B (right) locations. PODI spectra ( ) demonstrate better agreement with LES data ( ) compared with DMD ( ) ones. Both PODI and DMD results in the graphs have been obtained employing r = 40 modes. analysis of the reduced ow solutions at selected points inside and outside vortical wake regions indicated that both models were able to recover most of the model ow frequencies in the ranges of interest for the acoustic analysis. In particular, PODI showed notable capability to capture additional frequencies which were present in the original dataset, but not in the subset employed to train the model. As a consequence, both data-driven reduced models developed proved ecient and suciently accurate in predicting acoustic noises. Given the results obtained, introducing blu body geometrical parameterization as additional dimension for the PODI analysis is an interesting possibility for future works. A further future perspective which is currently being explored, is represented by the development of a reduced model employing the POD modal decomposition for the Galerkin projection of the continuity and momentum equations for the uid dynamic variables. Other studies could be devoted to better understand how the choice of the interpolator for PODI aects the results. Appendix A ROMs cross-comparison For additional clarity, the spectral and acoustic analyses discussed in the previous sections, cf.Figures 18 to 21, are here reported in single graphs directly comparing both DMD and PODI ROMs results against LES data. In particular, we are considering both DMD and PODI ROMs truncated at r = 40. The aim is to provide an easier comparison of the ROMs performance. Figure 22 depicts FFT spectrum of the pressure signal (inside and outside the wake, top plots) and the dipole and quadrupole acoustic noises (microphones A and B, bottom plots). In all plots DMD ( ) and PODI ( ) results are reported alongside LES ( ) results. The plots clearly indicate that, given equal modal truncation level, PODI models appear to accurately recover a wider range of frequencies compared to DMD models. Normalized amplitude spectrum Spectrum level [dB] Acknowledgements This work was partially supported by the project PRELICA, \Advanced methodologies for hydro- acoustic design of naval propulsion", supported by Regione FVG, POR-FESR 2014-2020, Piano Operativo Regionale Fondo Europeo per lo Sviluppo Regionale, and partially funded by European Union Funding for Research and Innovation | Horizon 2020 Program | in the framework of European Research Council Executive Agency: H2020 ERC CoG 2015 AROMA-CFD project 681447 \Advanced Reduced Order Methods with Applications in Computational Fluid Dynamics" P.I. Gianluigi Rozza. We gratefully thank professor Vincenzo Armenio, University of Trieste, for the general supervision, his precious comments and suggestions, during the development of the presented work. References [1] D.-G. Baek, H.-S. Yoon, J.-H. Jung, K.-S. Kim, and B.-G. Paik. Eects of the advance ratio on the evolution of a propeller wake. Computers & Fluids, 118:32{43, Sept. 2015. [2] E. Balaras, S. Schroeder, and A. Posa. Large-Eddy Simulations of Submarine Propellers. Journal of Ship Research, 59(4):227{237, Dec. 2015. [3] G. Berkooz, P. Holmes, and J. L. Lumley. The Proper Orthogonal Decomposition in the Analysis of Turbulent Flows. Annual Review of Fluid Mechanics, 25(1):539{575, Jan. 1993. [4] D. A. Bistrian and I. M. Navon. An improved algorithm for the shallow water equations model reduction: Dynamic Mode Decomposition vs POD. International Journal for Numerical Methods in Fluids, 78(9):552{580, Apr. 2015. [5] D. A. Bistrian and I. M. Navon. Randomized dynamic mode decomposition for nonintru- sive reduced order modelling. International Journal for Numerical Methods in Engineering, 112(1):3{25, 2017. [6] D. A. Bistrian and I. M. Navon. Eciency of randomised dynamic mode decomposition for reduced order modelling. International Journal of Computational Fluid Dynamics, 32(2- 3):88{103, 2018. [7] E. Bouhoubeiny and P. Druault. Note on the POD-based time interpolation from successive PIV images. Comptes Rendus M ecanique, 337(11-12):776{780, Nov. 2009. [8] K. S. Brentner and F. Farassat. Modeling aerodynamically generated sound of helicopter rotors. Progress in Aerospace Sciences, 39(2-3):83{120, Feb. 2003. [9] A. Broatch, J. Garc a-T scar, F. Roig, and S. Sharma. Dynamic mode decomposition of the acoustic eld in radial compressors. Aerospace Science and Technology, 90:388{400, 2019. [10] S. L. Brunton and J. N. Kutz. Data-Driven Science and Engineering: Machine Learning, Dynamical Systems, and Control. Cambridge University Press, 2019. [11] T. Bui-Thanh, M. Damodaran, and K. Willcox. Proper Orthogonal Decomposition Exten- sions for Parametric Applications in Compressible Aerodynamics. In 21st AIAA Applied Aerodynamics Conference, page 4213, 2003. [12] T. Bui-Thanh, M. Damodaran, and K. E. Willcox. Aerodynamic Data Reconstruction and Inverse Design Using Proper Orthogonal Decomposition. AIAA journal, 42(8):1505{1516, [13] E. F. Campana, D. Peri, Y. Tahara, and F. Stern. Shape optimization in ship hydrody- namics using computational uid dynamics. Computer Methods in Applied Mechanics and Engineering, 196(1-3):634{651, Dec. 2006. [14] M. Cianferra, V. Armenio, and S. Ianniello. Hydroacoustic noise from dierent geometries. International Journal Heat Fluid Flow, 70:348{362, 2018. 30 [15] M. Cianferra, S. Ianniello, and V. Armenio. Assessment of methodologies for the solution of the Ffowcs Williams and Hawkings equation using LES of incompressible single-phase ow around a nite-size square cylinder. Journal of Sound and Vibration, 453:1{24, 2019. [16] M. Cianferra, A. Petronio, and V. Armenio. Hydrodynamic Noise from a Propeller in Open Sea Condition. In Technology and Science for the Ships of the Future: Proceedings of NAV 2018: 19th International Conference on Ship & Maritime Research, pages 149{156. IOS Press, 2018. [17] M. Cianferra, A. Petronio, and V. Armenio. Non-linear noise from a ship propeller in open sea condition. Ocean Engineering, 191:106474, Nov. 2019. [18] C. Cintolesi, A. Petronio, and V. Armenio. Large eddy simulation of turbulent buoyant ow in a con ned cavity with conjugate heat transfer. Phys uids, 27, 2015. [19] M. F. de Pando, P. J. Schmid, and D. Sipp. A global analysis of tonal noise in ows around aerofoils. Journal of Fluid Mechanics, 754:5{38, July 2014. [20] N. Demo, M. Tezzele, G. Gustin, G. Lavini, and G. Rozza. Shape Optimization by Means of Proper Orthogonal Decomposition and Dynamic Mode Decomposition. In Technology and Science for the Ships of the Future: Proceedings of NAV 2018: 19th International Conference on Ship & Maritime Research, pages 212{219. IOS Press, 2018. [21] N. Demo, M. Tezzele, A. Mola, and G. Rozza. An Ecient Shape Parametrisation by Free- Form Deformation Enhanced by Active Subspace for Hull Hydrodynamic Ship Design Prob- lems in Open Source Environment. In Proceedings of ISOPE 2018: The 28th International Ocean and Polar Engineering Conference, volume 3, pages 565{572, 2018. [22] N. Demo, M. Tezzele, A. Mola, and G. Rozza. A complete data-driven framework for the ecient solution of parametric shape design and optimisation in naval engineering problems. In R. Bensow and J. Ringsberg, editors, Proceedings of MARINE 2019: VIII International Conference on Computational Methods in Marine Engineering, pages 111{121, 2019. [23] N. Demo, M. Tezzele, and G. Rozza. EZyRB: Easy Reduced Basis method. The Journal of Open Source Software, 3(24):661, 2018. [24] N. Demo, M. Tezzele, and G. Rozza. PyDMD: Python Dynamic Mode Decomposition. The Journal of Open Source Software, 3(22):530, 2018. [25] N. Demo, M. Tezzele, and G. Rozza. A non-intrusive approach for reconstruction of POD modal coecients through active subspaces. Comptes Rendus M ecanique de l'Acad emie des Sciences, DataBEST 2019 Special Issue, 347(11):873{881, November 2019. [26] V. Dolci and R. Arina. Proper Orthogonal Decomposition as Surrogate Model for Aerody- namic Optimization. International Journal of Aerospace Engineering, 2016. [27] P. Druault and C. Chaillou. Use of Proper Orthogonal Decomposition for reconstructing the 3D in-cylinder mean- ow eld from PIV data. Comptes Rendus M ecanique, 335(1):42{47, Jan. 2007. [28] N. B. Erichson, S. L. Brunton, and J. N. Kutz. Compressed dynamic mode decomposition for background modeling. Journal of Real-Time Image Processing, 16(5):1479{1492, Oct 2019. [29] J. E. Ffowcs Williams and D. L. Hawkings. Sound generation by turbulence and surfaces in arbitrary motion. Philosophical Transaction of Royal Society, 264:321{342, 1969. [30] M. Fossati and W. G. Habashi. Multiparameter Analysis of Aero-Icing Problems Using Proper Orthogonal Decomposition and Multidimensional Interpolation. AIAA Journal, 51(4):946{960, Apr. 2013. [31] T. Franz, R. Zimmermann, S. G ortz, and N. Karcher. Interpolation-based reduced-order modelling for steady transonic ows via manifold learning. International Journal of Compu- tational Fluid Dynamics, 28(3-4):106{121, Mar. 2014. 31 [32] M. Gadalla, M. Tezzele, A. Mola, and G. Rozza. BladeX: Python Blade Morphing. The Journal of Open Source Software, 4(34):1203, 2019. [33] F. Garotta, N. Demo, M. Tezzele, M. Carraturo, A. Reali, and G. Rozza. Reduced Order Isogeometric Analysis Approach for PDEs in Parametrized Domains. In M. D'Elia, M. Gun- zburger, and G. Rozza, editors, Quanti cation of Uncertainty: Improving Eciency and Technology: QUIET selected contributions, volume 137 of Lecture Notes in Computational Science and Engineering, pages 153{170. Springer International Publishing, Cham, 2020. [34] S. Georgaka, G. Stabile, K. Star, G. Rozza, and M. J. Bluck. A hybrid reduced order method for modelling turbulent heat transfer problems. Computers & Fluids, 208:104615, 2020. [35] S. A. Glegg and W. J. Devenport. Proper orthogonal decomposition of turbulent ows for aeroacoustic and hydroacoustic applications. Journal of Sound and Vibration, 239(4):767{ 784, 2001. [36] X. Gloerfelt, F. P erot, C. Bailly, and D. Juv e. Flow-induced cylinder noise formulated as a diraction problem for low Mach numbers. Journal of Sound and Vibration, 287(1-2):129{ 151, Oct. 2005. [37] Y. L. Guennec, J.-P. Brunet, F.-Z. Daim, M. Chau, and Y. Tourbier. A parametric and non-intrusive reduced order model of car crash simulation. Computer Methods in Applied Mechanics and Engineering, 338:186{207, Aug. 2018. [38] G. Haller. An objective de nition of a vortex. Journal of Fluid Mechanics, 525:1{26, Feb. [39] J. S. Hesthaven, G. Rozza, B. Stamm, et al. Certi ed Reduced Basis Methods for Parametrized Partial Dierential Equations. Springer, 2016. [40] S. Hijazi, G. Stabile, A. Mola, and G. Rozza. Data-Driven POD{Galerkin reduced order model for turbulent ows. Journal of Computational Physics, 416:109513, 2020. [41] S. Ianniello, R. Muscari, and A. D. Mascio. Ship underwater noise assessment by the acoustic analogy. Part I: nonlinear analysis of a marine propeller in a uniform ow. Journal of Marine Science and Technology, 18(4):547{570, July 2013. [42] S. Ianniello, R. Muscari, and A. D. Mascio. Ship underwater noise assessment by the acoustic analogy. Part II: hydroacoustic analysis of a ship scaled model. J Mar Sci Technol, 18(4):547{ 570, 2013. [43] E. Iuliano and D. Quagliarella. Proper Orthogonal Decomposition, surrogate modelling and evolutionary optimization in aerodynamic design. Computers & Fluids, 84:327{350, Sept. [44] G. Jourdain, L.-E. Eriksson, S. H. Kim, and C. H. Sohn. Application of dynamic mode decomposition to acoustic-modes identi cation and damping in a 3-dimensional chamber with baed injectors. Journal of Sound and Vibration, 332(18):4308{4323, Sept. 2013. [45] M. R. Jovanovi c, P. J. Schmid, and J. W. Nichols. Sparsity-promoting dynamic mode de- composition. Physics of Fluids, 26(2):024103, Feb 2014. [46] S. Karri, J. Charonko, and P. P. Vlachos. Robust wall gradient estimation using radial basis functions and proper orthogonal decomposition (POD) for particle image velocimetry (PIV) measured elds. Measurement Science and Technology, 20(4):045401, Feb. 2009. [47] J. E. Kerwin. Marine Propellers. Annual Review of Fluid Mechanics, 18(1):367{403, Jan. [48] P. Kumar and K. Mahesh. Large eddy simulation of propeller wake instabilities. Journal of Fluid Mechanics, 814:361{396, Feb. 2017. [49] K. Kunisch and S. Volkwein. Proper orthogonal decomposition for optimality systems. ESAIM: Mathematical Modelling and Numerical Analysis, 42(1):1{23, Jan. 2008. 32 [50] J. N. Kutz, S. L. Brunton, B. W. Brunton, and J. L. Proctor. Dynamic Mode Decomposition: Data-Driven Modeling of Complex Systems, volume 149. SIAM, 2016. [51] J. N. Kutz, X. Fu, and S. L. Brunton. Multiresolution Dynamic Mode Decomposition. SIAM Journal on Applied Dynamical Systems, 15(2):713{735, 2016. [52] S. Le Clainche and J. M. Vega. Higher Order Dynamic Mode Decomposition. SIAM Journal on Applied Dynamical Systems, 16(2):882{925, 2017. [53] L. S. Lorente, J. M. Vega, and A. Velazquez. Generation of Aerodynamics Databases Using High-Order Singular Value Decomposition. Journal of Aircraft, 45(5):1779{1788, Sept. 2008. [54] H. V. Ly and H. T. Tran. Modeling and control of physical processes using proper orthogonal decomposition. Mathematical and Computer Modelling, 33(1-3):223{236, Jan. 2001. [55] M. Mancinelli, T. Pagliaroli, R. Camussi, and T. Castelain. On the hydrodynamic and acoustic nature of pressure proper orthogonal decomposition modes in the near eld of a compressible jet. Journal of Fluid Mechanics, 836:998{1008, Dec. 2017. [56] I. Markovsky. Low-Rank Approximation. Communications and Control Engineering. Springer International Publishing, 2 edition, 2019. [57] A. D. Mascio, R. Muscari, and G. Dubbioso. On the wake dynamics of a propeller operating in drift. Journal of Fluid Mechanics, 754:263{307, July 2014. [58] C. Meneveau, T. Lund, and W. Cabot. A Lagrangian dynamic subgrid-scale model of tur- bulence. J Fluid Mech, 319:353{385, 1996. [59] K. E. Meyer, J. M. Pedersen, and O. Ozcan. A turbulent jet in cross ow analysed with proper orthogonal decomposition. Journal of Fluid Mechanics, 583:199{227, July 2007. [60] M. J. Mifsud, D. G. MacManus, and S. Shaw. A variable- delity aerodynamic model using proper orthogonal decomposition. International Journal for Numerical Methods in Fluids, 82(10):646{663, Apr. 2016. [61] M. J. Mifsud, S. T. Shaw, and D. G. MacManus. A high- delity low-cost aerodynamic model using proper orthogonal decomposition. International Journal for Numerical Methods in Fluids, pages n/a{n/a, 2009. [62] A. Mola, M. Tezzele, M. Gadalla, F. Valdenazzi, D. Grassi, R. Padovan, and G. Rozza. Ecient Reduction in Shape Parameter Space Dimension for Ship Propeller Blade Design. In R. Bensow and J. Ringsberg, editors, Proceedings of MARINE 2019: VIII International Conference on Computational Methods in Marine Engineering, pages 201{212, 2019. [63] T. W. Muld, G. Efraimsson, and D. S. Henningson. Flow structures around a high-speed train extracted using Proper Orthogonal Decomposition and Dynamic Mode Decomposition. Computers & Fluids, 57:87{97, 2012. [64] A. Naja -Yazdi, G.A.Bres, and L. Mongeau. An acoustic analogy formulation for moving sources in uniformly moving media. Proc R Soc Lond, A467:144{165, 2011. [65] Z. Nitzkorski and K. Mahesh. A dynamic end cap technique for sound computation using the Ffowcs Williams and Hawkings equations. Physics of Fluids, 26(11):115101, Nov. 2014. [66] U. Piomelli, L. S. Craig, and S. Sarkar. On the computation of sound by large-eddy simula- tions. Journal of Engineering and Mathematics, 32:217{236, 1997. [67] A. Posa, R. Broglia, M. Felli, M. Falchi, and E. Balaras. Characterization of the wake of a submarine propeller via Large-Eddy simulation. Computer and Fluids, 184:138{152, 2019. [68] J. L. Proctor, S. L. Brunton, and J. N. Kutz. Dynamic Mode Decomposition with Control. SIAM Journal on Applied Dynamical Systems, 15(1):142{161, 2016. [69] A. Quarteroni, A. Manzoni, and F. Negri. Reduced Basis Methods for Partial Dierential Equations: An Introduction, volume 92. Springer, 2015. 33 [70] M. Ripepi, M. Verveld, N. Karcher, T. Franz, M. Abu-Zurayk, S. G ortz, and T. Kier. Reduced-order models for aerodynamic applications, loads and MDO. CEAS Aeronautical Journal, 9(1):171{193, 2018. [71] C. W. Rowley, T. Colonius, and R. M. Murray. Model reduction for compressible ows using POD and Galerkin projection. Physica D: Nonlinear Phenomena, 189(1-2):115{129, Feb. [72] C. W. Rowley, I. Mezi c, S. Bagheri, P. Schlatter, and D. S. Henningson. Spectral analysis of nonlinear ows. Journal of uid mechanics, 641:115{127, 2009. [73] G. Rozza, M. W. Hess, G. Stabile, M. Tezzele, and F. Ballarin. Basic Ideas and Tools for Projection-Based Model Reduction of Parametric Partial Dierential Equations . In P. Benner, S. Grivet-Talocia, A. Quarteroni, G. Rozza, W. H. A. Schilders, and L. M. Silveira, editors, Handbook on Model Order Reduction, volume 2, chapter 1. De Gruyter, In Press, 2020. [74] G. Rozza, M. H. Malik, N. Demo, M. Tezzele, M. Girfoglio, G. Stabile, and A. Mola. Ad- vances in Reduced Order Methods for Parametric Industrial Problems in Computational Fluid Dynamics. In R. Owen, R. de Borst, J. Reese, and P. Chris, editors, ECCOMAS ECFD 7 - Proceedings of 6th European Conference on Computational Mechanics (ECCM 6) and 7th European Conference on Computational Fluid Dynamics (ECFD 7), pages 59{76, Glasgow, UK, 2018. [75] S. H. Rudy, J. N. Kutz, and S. L. Brunton. Deep learning of dynamics and signal-noise decomposition with time-stepping constraints. J. Comput. Physics, 396:483{506, 2018. [76] H. Sakamoto and H. Haniu. A Study on Vortex Shedding From Spheres in a Uniform Flow. Journal of Fluids Engineering, 112(4):386{392, 12 1990. [77] F. Salmoiraghi, F. Ballarin, G. Corsi, A. Mola, M. Tezzele, and G. Rozza. Advances in geometrical parametrization and reduced order models and methods for computational uid dynamics problems in applied sciences and engineering: overview and perspectives. ECCO- MAS Congress 2016 - Proceedings of the 7th European Congress on Computational Methods in Applied Sciences and Engineering, 1:1013{1031, 2016. [78] F. Salmoiraghi, A. Scardigli, H. Telib, and G. Rozza. Free-form deformation, mesh mor- phing and reduced-order methods: enablers for ecient aerodynamic shape optimisation. International Journal of Computational Fluid Dynamics, 32(4-5):233{247, 2018. [79] R. D. Sandberg and N. D. Sandham. Direct numerical simulation of turbulent ow past a trailing edge and the associated noise generation. Journal of Fluid Mechanics, 596:353{385, Jan. 2008. [80] P. J. Schmid. Dynamic mode decomposition of numerical and experimental data. Journal of uid mechanics, 656:5{28, 2010. [81] P. J. Schmid. Application of the dynamic mode decomposition to experimental data. Exper- iments in uids, 50(4):1123{1130, 2011. [82] P. J. Schmid, L. Li, M. P. Juniper, and O. Pust. Applications of the dynamic mode decom- position. Theoretical and Computational Fluid Dynamics, 25(1-4):249{259, 2011. [83] A. Seena and H. J. Sung. Dynamic mode decomposition of turbulent cavity ows for self- sustained oscillations. International Journal of Heat and Fluid Flow, 32(6):1098{1110, Dec. [84] J. H. Seo, Y. J. Moon, and B. R. Shin. Prediction of cavitating ow noise by direct numerical simulation. Journal of Computational Physics, 227(13):6511{6531, June 2008. [85] H. Seol, B. Jung, J.-C. Suh, and S. Lee. Prediction of non-cavitating underwater propeller noise. Journal of Sound and Vibration, 257(1):131{156, Oct. 2002. [86] C. Seror, P. Sagaut, C. Bailly, and D. Juv e. On the radiated noise computed by large-eddy simulation. Physics of Fluids, 13(2):476{487, Feb. 2001. 34 [87] W. Shen, T. K. Patel, and S. A. Miller. Extraction of Large-Scale Coherent Structures from Large Eddy Simulation of Supersonic Jets for Shock-Associated Noise Prediction. In AIAA Scitech 2020 Forum. American Institute of Aeronautics and Astronautics, Jan. 2020. [88] L. Sirovich. Turbulence and the dynamics of coherent structures. II. symmetries and trans- formations. Quarterly of Applied Mathematics, 45(3):573{582, Oct. 1987. [89] G. Stabile, F. Ballarin, G. Zuccarino, and G. Rozza. A reduced order variational multiscale approach for turbulent ows. Advances in Computational Mathematics, 45:2349{2368, 2019. [90] G. Stabile, S. Hijazi, A. Mola, S. Lorenzi, and G. Rozza. POD-Galerkin reduced order methods for CFD using Finite Volume Discretisation: vortex shedding around a circular cylinder. Communications in Applied and Industrial Mathematics, 8(1):210{236, 2017. [91] G. Stabile and G. Rozza. Finite volume POD-Galerkin stabilised reduced order methods for the parametrised incompressible Navier{Stokes equations. Computers & Fluids, 173:273{284, sep 2018. [92] P. Stegeman, A. Ooi, and J. Soria. Proper Orthogonal Decomposition and Dynamic Mode Decomposition of Under-Expanded Free-Jets with Varying Nozzle Pressure Ratios. In Insta- bility and Control of Massively Separated Flows, pages 85{90. Springer, 2015. [93] M. Tezzele, N. Demo, M. Gadalla, A. Mola, and G. Rozza. Model Order Reduction by Means of Active Subspaces and Dynamic Mode Decomposition for Parametric Hull Shape Design Hydrodynamics. In Technology and Science for the Ships of the Future: Proceedings of NAV 2018: 19th International Conference on Ship & Maritime Research, pages 569{576. IOS Press, 2018. [94] M. Tezzele, N. Demo, A. Mola, and G. Rozza. An integrated data-driven computational pipeline with model order reduction for industrial and applied mathematics. Special Volume ECMI, Springer, In Press, 2020. [95] M. Tezzele, N. Demo, A. Mola, and G. Rozza. PyGeM: Python geometrical morphing. Software Impacts, page 100047, 2020. [96] M. Tezzele, N. Demo, and G. Rozza. Shape optimization through proper orthogonal decom- position with interpolation and dynamic mode decomposition enhanced by active subspaces. In R. Bensow and J. Ringsberg, editors, Proceedings of MARINE 2019: VIII International Conference on Computational Methods in Marine Engineering, pages 122{133, 2019. [97] M. Tezzele, N. Demo, G. Stabile, A. Mola, and G. Rozza. Enhancing CFD predictions in shape design problems by model and parameter space reduction. Advanced Modeling and Simulation in Engineering Sciences, 7(40), 2020. [98] M. Tezzele, F. Salmoiraghi, A. Mola, and G. Rozza. Dimension reduction in heterogeneous parametric spaces with application to naval engineering shape design problems. Advanced Modeling and Simulation in Engineering Sciences, 5(1), 2018. [99] F. Valdenazzi, F. Conti, S. Gaggero, D. Grassi, C. Vaccaro, and D. Villa. A practical tool for the hydro-acoustic optimization of naval propellers. In R. Bensow and J. Ringsberg, editors, Proceedings of MARINE 2019: VIII International Conference on Computational Methods in Marine Engineering, pages 296{308, 2019. [100] H. G. Weller, G. Tabor, H. Jasak, and C. Fureby. A tensorial approach to computational continuum mechanics using object-oriented techniques. Computers in physics, 12(6):620{631, [101] M. O. Williams, I. G. Kevrekidis, and C. W. Rowley. A Data{Driven Approximation of the Koopman Operator: Extending Dynamic Mode Decomposition. Journal of Nonlinear Science, 25(6):1307{1346, 2015. [102] D. Xiao, F. Fang, C. Pain, and G. Hu. Non-intrusive reduced-order modelling of the Navier- Stokes equations based on RBF interpolation. International Journal for Numerical Methods in Fluids, 79(11):580{595, July 2015. 35 [103] Q. Zhang, Y. Liu, and S. Wang. The identi cation of coherent structures using proper or- thogonal decomposition and dynamic mode decomposition. Journal of Fluids and Structures, 49:53{72, 2014.
Mathematics – arXiv (Cornell University)
Published: Jun 25, 2020
You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!
Read and print from thousands of top scholarly journals.
Already have an account? Log in
Bookmark this article. You can see your Bookmarks on your DeepDyve Library.
To save an article, log in first, or sign up for a DeepDyve account if you don’t already have one.
Copy and paste the desired citation format or use the link below to download a file formatted for EndNote
Access the full text.
Sign up today, get DeepDyve free for 14 days.
All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.