Numerical Simulations of Wave-Induced Flow Fields around Large-Diameter Surface-Piercing Vertical Circular Cylinder

A computational analysis is performed on the diffraction of water waves induced by large-diameter, surface-piercing, vertical circular cylinder. With reference to linear-wave cases, the phenomenon is preliminarly considered in terms of velocity potential, a simplified theoretical framework in which both hypotheses of inviscid fluid and irrotational flow are incorporated. Then, and as a first-approximation analysis, the Euler equations in primitive variables are considered (a framework in which the fluid is still handled as inviscid, but the field can be rotational). Finally, the real-fluid behavior is analyzed, by numerically integrating the full Navier-Stokes equations (viscous fluid and rotational field) in their velocity-pressure formulation, by following the approach of the Direct Numerical Simulation (DNS, no models are used for the fluctuating portion of the velocity field). For further investigation of the flow fields, the swirling-strength criterion for flow-structure extraction, and the Karhunen-Loève (KL) decomposition technique for the extraction of the most energetic flow modes respectively, are applied to the computed fields. It is found that remarkable differences exist between the wave-induced fields, as derived within the different computing frameworks tested.


Introduction
The phenomenon of diffraction of small-amplitude water waves impinging on large bodies has been studied in the recent past by several authors both numerically and experimentally, where the results related to the case of large-diameter surface-piercing vertical circular cylinder are often compared against the theoretical close-form potential solution devised by MacCamy and Fuchs [1].In the numerical field, Mansour et al. [2] performed a study of the interaction of linear waves with a uniform-bottom-mounted surface-piercing cylinder.Roy and Ghosh [3] presented a numerical solution, in terms of velocity potential, of shallow-water wave force on vertically-submerged circular thin plates in different configurations, where the plates are submerged in water near the shore on a uniformly-sloping bottom.Kim et al. [4] addressed the issue of the diffraction of water waves by a bottom-fixed vertical circular cylinder using the boundary-element method, within the framework of the linear-potential theory.It has to be recognized that the agreement between potential-flow numerical results and potential-theory results, is rather satisfactory.In the experimental field, Hellstrom and Rundgren [5], Laird [6], Bonnefille and Germain [7], and Isaacson [8] studied the wave-elevation distribution around large circular islands and lighthouses within limited ranges of wave interaction parameters, while Galvin and Hallermeier [9], Hallermeier [10], and Haney and Herbich [11] considered piles of various cross sections.Chakrabarti and Tam [12], and Niedzweki and Duggal [13] performed experiments at values of the Keulegan-Carpenter number falling in the diffraction regime.They found increases in the runup, with respect to the potential theory, when either the incident-wave steepness or the scatter parameter were increased.Niedzweki and Duggal [14], and Niedzweki and Huston [15] investigated wave runups and forces on cylindrical bodies in a wide range of wave characteristic parameters.They found that linear diffraction theory satisfactorily predicts the wave forces, but strongly underestimates the wave runups.Subsequently, Martin et al. [16], Mase et al. [17], and Morris-Thomas and Thiagarajan [18] found that linear diffraction theory underpredicts the runup extent by significant margins.De Vos et al. [19], in investigating technical issues related to cylindrical foundations of offshore wind turbines, also found that the wave runups were actually much higher than those predicted by linear theory.
Overall, it can be noticed that most of the numerical works concerning wave phenomena are executed by adopting the velocity-potential approach, which incorporates the two rather restrictive assumptions of inviscid fluid and irrotational flow.A first step toward more reliable investigations can be represented by the numerical integration of the primitive-variable Euler equations (a framework of inviscid fluid, but no necessarily irrotational flow), but it is the numerical integration of the full Navier-Stokes equations (viscous fluid and rotational flow) the physically-consistent approach to be followed in order to obtain flow fields correctly representing the numerical equivalent of the phenomena at hand.In the numerical literature on wave phenomena still rare are the works in which both hypotheses of viscous fluid and rotational flow are assumed.Among the firsts, Choi et al. [20] performed a numerical simulation of the tsunami runup around a conical island by using a RNG RANS-type of closure model for the Reynods-stress term of the Navier-Stokes equations (see Alfonsi [21] for an extensive review of the RANS modeling technique, and refer to the Nomenclature for the meaning of symbols and acronyms used throughout the whole text).
In this work the aforementioned issues are addressed, with reference to the particular case of diffraction of water waves by large-diameter, surface-piercing, vertical circular cylinder.The velocity-potential approach, the numerical integration of the primitive-variable Euler equations, and the numerical integration of the primitive-variable full Navier-Stokes equations are considered for the analysis of given wave cases, and-in the case of the Navier-Stokes equations-the approach is followed of the Direct Numerical Simulation (DNS, no models for the fluctuating portion of the velocity field are used in the calculations, see Alfonsi [22,23] for extensive DNS reviews).For further investigation of the flow fields, the swirling-strength criterion for flow-structure extraction, and the Karhunen-Loève decomposition (KL) technique for the calculation of the most energetic modes of the flow, are applied to the computed fields.Interestingly, a complex configuration of both viscous-flow structures and high-energy flow modes is unveiled in the vicinity of the cylinder wall, both at the free surface and immediately under the free surface.The present work is organized as follows.In Section 2, concise overviews of different wave regimes are given, while in Section 3 the different numerical techniques used are presented, in conjunction with the results of an extensive and systematic work of validation that has been performed on the Navier-Stokes solver and a discussion about the criteria of accuracy that have been adopted in the Navier-Stokes calculations (mainly related to the adoption of the DNS technique).In Section 4, the methods used for the extraction from the computed fields of the flow structures and the most energetic flow modes, are synthetically outlined, and in Section 5 the results are presented.Concluding remarks are given at the end.

Overview of Wave Regimes
It is known (Sarpkaya and Isaacson [24]) that in the case of waves impinging on a vertical circular cylinder, wave diffraction becomes relevant when 0.2 D L > in conjunction with the condition that the value of the Keulegan-Carpenter number: in which: ( ) ( ) and where Equation (1) reflects the fact that the maximum value of the wave-steepness parameter is: By representing the different wave regimes in graphical form in terms of KC and D L (Figure 1), the graphical implementation of Equation (5) translates in the identification of a zone (Zone A) in which wave diffraction is increasingly important at increasing values of D L , and the phenomenon is linear.In Figure 1, a sub-zone (within the former) can be also distinguished (Zone B), in which diffraction is influenced by nonlinear effects, according to the limit introduced by Isaacson [25] for this condition: Zone A of Figure 1 is then characterized by the fact that diffraction is the sole phenomenon occurring in the process of interaction between the wave and the vertical circular cylinder.10), ( 11), ( 14) and ( 15); (∆∆) parameters of the experimental waves [14] falling under the limit mirrored by Equation (5).

Viscous-Fluid Wave Framework
When a circular cylinder is exposed to a viscous-fluid oscillatory flow, the resulting flow fields are influenced by the value of the wave-field Reynolds number: Oscillatory flows have been investigated experimentally in pioneering works due to Honji [26], and Sarpkaya [27].
In Figure 2 the different flow regimes that establish under the influence of both the Keulegan-Carpenter and the Reynolds number, are represented.Different zones can be distinguished: (i) Zone A is a no-separation zone; (ii) Zone B is a zone with no separation, but where the boundary layer becomes turbulent; (iii) Zone C is a separation zone with the establishment of the Honji vortices [26]; (iv) Zone D is a separation zone where a pair of symmetric vortices establishes downstream from the cylinder; (v) Zone E is a separation zone where a pair of symmetric vortices establishes downstream from the cylinder and turbulence exists over the cylinder surface.In Figure 2, the 50 wave cases already represented in Figure 1 10), (11), ( 14) and ( 15); (∆∆) parameters of the experimental waves [14].

Velocity-Potential Solution
A close-form solution in terms of the velocity potential has been devised by MacCamy and Fuchs [1] for the problem of wave diffraction caused by large-diameter, bottom-fixed, surface-piercing, vertical circular cylinder.According to this solution (Sarpkaya and Isaacson [24], Sumer and Fredsøe [28]), the total velocity-potential function ( ϕ ) is given by the following expression: where: ( ) and 1 The velocity field and other relevant quantities can be obtained from Equation (8).In particular, the maximum runup at the cylinder surface ( max R ) and the maximum force on the cylindrical body ( max F ) are given by the following expressions: where: ( ) ( ) ( ) ( ) ( ) ( ) ( ) For a preliminary linear-theory evaluation of the quantities defined in Equations ( 10) and ( 11), 50 wave-cases have been considered, within the ranges of characteristic parameters as reported in Table 1.With reference to Figure 1, it can be noticed that all the 50 selected wave cases are strictly linear (all the points corresponding to the 50 KC − D L couples of values are located under the limit mirrored by Equation ( 5)).The corresponding nondimensional expressions ( max nd R ), ( max nd F ) of the quantities defined in Equations ( 10) and ( 11) are given by: [ ] ( ) As for the computation of the ϕ -derived wave-flow fields, the following procedure has been followed: (i) a wave case has been selected for inviscid-fluid flow-field analysis (the same used in Section 3.2); (ii) the velocity potential ( ϕ ), the related velocity field ( i i u x = ∂ϕ ∂ ), and some other ϕ -derived quantities have been calculated, on the basis of expression Equation (8).The values of these quantities have been evaluated numerically, on the same computing domain used in the subsequent simulations (Figure 3), and the same computational grid used for the Euler wave-flow simulations (Section 3.2).

Numerical Integration of the Euler Equations
The system of the three-dimensional, time-dependent Euler equations in primitive variables with body forces is considered (the fluid is incompressible and inviscid, Einstein summation convention applies to repeated indexes, i, j = 1, 2, 3): For the execution of the calculations, the Flow-3D finite-volume computational code has been used.In this code, the free-surface condition is handled with the Volume Of Fluid method (VOF), as devised by Hirt and Nichols [29], which has extensively proven to be able to accurately tracking a wave interface (see also at Section 3.3).In particular, the time-marching procedure includes three main steps (Flow Science [30]): (i) evaluation of the velocity in each cell using the initial conditions (or previous-time-step values) for the advective pressures (and/or other accelerations) on the basis of appropriate explicit approximations of the governing equations; (ii) adjustment of the pressure in each cell to satisfy the continuity equation; (iii) updating of the fluid free surface to give the new fluid configuration based on the volume-of-fluid value in each cell.For the present Euler wave-flow calculations, a solution scheme based on the Generalized Minimal Residual (GMRES) method has been used.
Boundary conditions of free-slip (and zero wall-normal velocity) have been imposed on the x-y bottom plane of the computing domain (Figure 3), and at the cylinder external surface.On the two x-z lateral boundary planes, and on the y-z end-plane of the computing domain, outflow conditions have been set, while the free-surface condition holds at the wave surface.On the y-z inlet plane, the flow field corresponding to an incident linear wave (with given parameters, see at Table 2) is generated.
Recall that, by applying the curl operator to Equations ( 16) and ( 17), after some manipulations, one obtains the vorticity transport equations for an inviscid-fluid flow: where: ( ) is the vorticity.The terms on the lhs of Equation ( 18) represent the rate-of-change of particle vorticity, while the term on the rhs is the vortex-stretching term (that has no analogous term in the primitive-variable Euler equations).In a condition as depicted by Equations ( 18) and ( 19), the vorticity-if any-only changes by the turning and/or stretching of the vortex lines (as a line stretches the vorticity increases, as it shrinks, the vorticity decreases).Then, if there is a non-zero presence of vorticity in an inviscid-fluid flow, due to the fact that the particles of an inviscid fluid are not subjected to net viscous forces, one faces a situation in which the presence of vorticity is due to the fact that viscous forces have been active previously, i.e., in the past history of the flow (a previous phase in the flow lifetime in which, as an example, the flow field could have interacted viscously with other bodies, and/or could have been simulated within a viscous-fluid framework of hypotheses).When such a case actually verifies, the result is that one faces an inviscid-fluid flow carrying a field of previously-generated vorticity, where stretching and turning are the only mechanisms that influence the evolution of the vorticity field, the latter being governed by Equations ( 18) and ( 19) (no vorticity is diffused in the absence of viscosity).In the case here at hand, a viscous-fluid wave initial condition has been used to start the simulations, and this initial condition actually represents the source of the previously-generated vorticity for the Euler fields (see also at Section 5.3).
For the present-work Euler wave-flow simulations, a linear case has been selected with characteristic parameters included in the range of values reported in Table 1 (so belonging to Zone A of Figure 1).These parameters are extensively reported in Table 2, in conjunction with the related main computational meaningful quantities used in the simulations.As for the grid resolution, a number of tests has been performed before reaching the final configuration of the grid, employing an increasing number of grid points, up to that finally utilized.22 36.0450 360 6 72.90 10 × 0.20 0.10

Numerical Integration of the Navier-Stokes Equations
The system of the three-dimensional, non-dimensional, time-dependent full Navier-Stokes equations in primitive variables (body forces are included) is considered (the fluid is incompressible and viscous, Einstein summation convention applies to repeated indexes, i, j = 1, 2, 3): Variables and operators have been here nondimensionalized by the cylinder diameter ( D ) for lengths, by the group ( Dg ) for velocities, by ( Dg ρ ) for pressures, by ( D g ) for time, where Re NS is the Reynolds number resulting from the above nondimensionalization procedures.The Reynolds number Re NS and the wave-field Reynolds number defined in Equation ( 7) are related as follows: ( ) For the execution of the calculations, the Flow-3D computational code [30] has been further used, that solves the fluid governing equations using the finite-volume technique.In this code the free-surface condition is handled with the VOF method [29], and the Fractional Area Volume Obstacle Representation (FAVOR) technique is used for the modeling of complex geometries.
According to the VOF method, the portion of the volume occupied by the obstacle in each computational cell is defined at the beginning of the work.The fluid fraction in each cell is calculated, and the equations for the fluid-fraction volumes are formulated using the FAVOR function.The finite-volume technique-on staggered cuboid computational cells-is used for the discretization of the governing equations.The meshing technique does not induce any mesh distortion during transients, and a multi-block meshing technique is available, to provide higher resolutions in the calculations, where needed.With respect to time, most terms in the equations are evaluated explicitly, exception made for the pressures.Pressures and velocities are coupled implicitly using time-advanced pressures in the momentum equation and time-advanced velocities in the mass-conservation equation.The coupled sets of equations resulting from this semi-implicit formulation are solved iteratively.For the present Navier-Stokes wave-flow calculations, a modified Alternating Direction Implicit (ADI) scheme has been used [30].
Boundary conditions of no-slip (and zero wall-normal velocity) at the cylinder external walls have been imposed.On the two x-z lateral boundary planes of the computing domain (Figure 3), symmetry conditions have been set, on the y-z end-plane outflow conditions have been enforced, while the free-surface condition holds at the wave surface.At the y-z inlet plane, the flow field corresponding to the incident wave of the case at hand is generated.
For the present Navier-Stokes wave-flow calculations (and also for the other simulations), a specially-assembled computing system has been used, that includes 2 Intel Xeon 5660 exa-core multi-core processors (a total of 12 CPUs available), a maximum of 48 GB of RAM, and up to 1.8 TB of mass memory (see also ).
Recall that, by applying the curl operator to the Navier-Stokes equations, after some manipulations, one obtains the system of the vorticity-transport equations: Now, the first term on the rhs of Equation ( 24) represents the rate-of-deformation of the vortex lines, while the second term on the rhs of Equation ( 24) represents the vorticity rate-of-diffusion, due to fluid viscosity.Equations ( 24) and ( 25) are of utmost importance to fluid mechanics and wave phenomena in particular, in that they mirror, in terms of vorticity field, the effects of the viscous forces.An unbalanced shear stress can exist only when the vorticity is nonzero, so that the existence of vorticity means that a particle is (or in its past history has been) subjected to net viscous forces.As for the representation of the results arising from numerical simulations, the vorticity field has been often used as a tool to gain new insights on flow phenomena, and in particular in the case of the flow around a circular cylinder (among others, Alfonsi [34][35][36], Alfonsi and Giorgini [37,38]).Though, the need of a deep understanding of many complex flow phenomena has soon motivated the consideration of the vorticity field as associated to the flow vortical structures (Alfonsi [39,40]), in a distinction with respect to the background, noncoherent vorticity field (see at Section 4.1).

Validation of the Navier-Stokes Solver
Before actually running the present-work numerical simulations, an extensive work has been performed in order to appropriately validate the Flow-3D Navier-Stokes solver.Flow-3D is a widely-known computational code, and in the recent past several authors have performed scientific works directed to its validation.Among others, Kocaman et al. [41], Cho [42], Jin and Meng [43], Kim et al. [44], and Kayser and Gabr [45] have validated the Flow-3D solver in steady and unsteady problems involving cylinders and free surfaces, always obtaining rather satisfactory results.Though, in the present context, the issue of validation is in some way subtle.In performing the present numerical simulations by using the DNS approach, for the evaluation of the quality of the results, it becomes essential that the fluctuating portions of the velocity fields are correctly caught by the calculations.Thus, for validation, a rather demanding problem under this aspect has been chosen, namely the flow of a viscous fluid in a plane channel (recall that all the wave cases considered in present work fall in the no-separation, shear-flow Zone A of Figure 2).The plane-channel-flow case has a long tradition in the field of the numerical integration of the Navier-Stokes equations with the DNS technique, so becoming a true reference case.In spite of its apparently-simple geometry, it probably represents the most complex shear-flow case one can find, with regard to the possibility for the flow of developing a full cascade of scales, that only depend on the value of the Reynolds number that one is able to reach (see [23] and references therein).Then, the Flow-3D Navier-Stokes code has been validated in the plane-channel flow case against a reference, rather accurate code, namely the mixed spectral-finite difference Navier-Stokes solver originally developed by Alfonsi et al. [46] and further tested by Passoni et al. [47][48][49].This latter solver has extensively proven to be able to accurately calculate several quantities that characterize shear flows (Alfonsi and Primavera [50]).In the framework of the validation procedure, high-order statistical moments of the fluctuating portions of the velocity fields have been computed, so that, according to the fact that the validation process has given rather satisfactory results, one is now confident that the flow scales, the flow structures, and the most energetic flow modes that have been calculated in the present work, are correct.
Results and comparisons are now reported for the computational test-case simulations (see also at Table 3).In Figure 4 the mean-velocity profile (in wall coordinates), as calculated with the Flow-3D solver, is compared with the law of the wall and with the results given by the reference solver.Figure 4 shows that the viscous sublayer is well resolved, yielding the linear velocity distribution for 5 y + < .At larger distances from the wall ( 30 y + > ) the logarithmic velocity distribution with von Kármán constant 0.4 κ = , and additive constant 5.5 C = , is satisfactorily followed.According to the Flow-3D solver, the value of ( ) u S ′ at the position ( peak y + ) of the ( ) rms u′ peak value, is 0.1419 − , rather close to zero (see also Mochizuki and Nieuwstadt [51], Alfonsi [52]), while the flatness factor ( v F ′ ) assumes a value of about 24 approaching the walls, so denoting the highly intermittent character of the normal velocity fluctuation near the solid border (see also Xu et al. [53], Alfonsi [54], Alfonsi et al. [55,56]).Overall, the comparisons between the results given by the Flow-3D Navier-Stokes solver and those given by the reference solver, are rather satisfactory.

Accuracy of Calculations and Computing Procedures
In the Navier-Stokes simulations of present work, no models for the fluctuating portion of the velocity field have been used, so that the viscous-flow governing equations are solved following the DNS approach.Direct Numerical Simulation consists in solving numerically the full Navier-Stokes equations with physically-consistent accuracy in space and time, in such a way as to resolve all the essential flow scales.If the mesh is fine enough, and the numerical scheme is designed to minimize the dispersion and dissipation errors, one obtains an accurate three-dimensional, time-dependent solution of the equations at hand.As concerns the smallest-scale-resolution requirements in DNS calculations, the criterion is usually followed of resolving the Kolmogorov space and time microscales: or, at most, limited multiples of the latter.
As concerns the procedure for the evaluation of the Kolmogorov scales as defined in Equation ( 26), they can be evaluated by estimating the average rate of dissipation of kinetic energy per unit mass ( ε ), as obtained from some mean-flow quantities (Bakewell and Lumley [57]).In the case of present work, one obtains: where ( f F ) is the portion-due-to-friction of the total in-line force exerted on the cylinder.
Recall that, for small values of the Keulegan-Karpenter number and intermediate values of the Reynolds number, it is possible to determine the in-line forces on the basis of the asymptotic theory of Bearman et al. [58] (see also Morison et al. [59]).On this basis, the latter quantity (per unit length) is given by: ( ) so that, for the wave-case of present work (see also at Table 4), one obtains: For the present Navier-Stokes wave-flow simulations, a linear wave case has been selected, with characteristic parameters included in the range of values reported in Table 1 (so belonging to Zone A of Figure 1).These parameters are extensively reported in Table 4, in conjunction with the related main computational parameters.
The computing domain (Figure 3) has been discretized into a 5-block, multiblock computational grid (Figure 9).The partitions of the computing domain occupied by each grid-point block and the grid-point distances ( , , x y z Δ Δ Δ ) within each block or group of blocks, are reported in Table 4.It is to be noted in particular that: (i) along the z-direction, z Δ assumes the same constant value through the different blocks; (ii) as concerns x Δ and y Δ along the xand the y-directions, respectively, the grid point-distances are constant in the peripheral blocks 1, 2, 4, 5, while within the internal block 3, they are distributed non-uniformly according the hyperbolic-tangent grid-stretching law: where 1.7 PP = and 1.9 QQ = ( 3 L′ is the half length of block 3, both along x and y, see at Figure 9).With the aforementioned grid-point configuration, the space resolutions ( x Δ ) and ( y Δ ) at the cylinder wall are able to resolve less than one half of the Kolmogorov spatial microscale, while at the domain external boundaries they are still able to resolve less than two and a half η.These latter values are perfectly in the ranges usually adopted in DNS ( [23]).( ) ) ) 3.9 3.9 blocks 1-5 3.9 3.9 x − < < 3.9 3.9 y − < <

Overview of Swirling-Strength Criterion for Flow-Structure Extraction
Of the existing techniques for the extraction of flow vortical structures (Alfonsi [39,40], Alfonsi and Primavera [60,61]) in the present work the criterion of the imaginary part of complex eigenvalue pair of the velocity-gradient tensor (the swirling-strength or also the ci λ criterion) as devised by Zhou et.al. [62], has been adopted.By considering the system of the flow governing equations, an arbitrary point O can be chosen in the field, and a Taylor-series expansion of each velocity component can be performed in terms of space coordinates with the origin in O, so that the first-order pointwise linear approximation at that point becomes: where ij i j A u x = ∂ ∂ is the velocity-gradient tensor.If O is located at a critical point, the zero-order terms of ( 22) are zero, and from the characteristic equation of ij A , one has: where: ( are the scalar invariants of the velocity-gradient tensor (tr is trace, det is determinant).In the case of incompressible flow, 0 P = , so that Equation (32), the second invariant ( Q ) of Equation (33b), and the discriminant ( Dsc ) of the characteristic equation of ij A , become respectively: when ( 0 Dsc > ), the velocity-gradient tensor has one real eigenvalue ( 1 λ ) and a pair of complex-conjugate eigenvalues ( 2 3 , λ λ ).Zhou et al. [62] adopted the criterion of identifying vortices by visualizing isosurfaces of prescribed values of the imaginary part of the complex-eigenvalue pair of the velocity-gradient tensor.In particular, one has: The swirling strength ( ci λ ) represents a measure of the local swirling rate inside a vortical structure, so that isosurfaces of the imaginary part of the complex eigenvalue pair of the velocity-gradient tensor can be used to visualize vortices.The method is frame independent and, due to the fact that the eigenvalue pair is complex only in regions of local circular or spiralling streamlines, it automatically eliminates regions having no local spiralling motion.The method has proven to give rather satisfactory results in several different flow cases (Alfonsi [63], Alfonsi and Primavera [64,65], Alfonsi et al. [66-70]).

Overview of Karhunen-Loève Decomposition Technique
Within the group of the Proper Orthogonal Decomposition (POD) methods (see Liang et al. [71], Alfonsi [39,40]) a powerful technique for the extraction of the most energetic modes of a flow is the Karhunen-Loèwe (KL) decomposition (see also Lumley [72], Sirovich [73]).By considering an ensemble of temporal realizations of a non-homogeneous, square integrable (with finite kinetic energy), three-dimensional velocity field ( , ) , on a finite domain Λ , one wants to find the most similar function to the elements of the ensemble on average (operator ⋅ ), i.e., to determine the highest mean-square correlated structure with all the elements of the ensemble.This corresponds to find a deterministic vector function ( ) where , , , 1, 2,3 i j l k = .
In Equation ( 36), ( , ) ( , ) is the two-point velocity-correlation tensor, and the maximum value of ( ) corresponds to the largest eigenvalue λ of ij R .When Λ is bounded, based on the Hilbert-Schmidt theory, there exists a denumerable infinity of solutions of Equation (36), where these solutions are called the empirical eigenfunctions [ ( ) ( ) , normalized, i.e., ( ) ( ) The orthonormality of the empirical eigenfunctions implies that structures of different order do not interact which each other in their contribution to second order statistics.To each empirical eigenfunction is associated a real and positive eigenvalue ( ) n λ ( ij R is non-negative by construction), and the eigenfunctions form a complete set.Thus, every member of the originary ensemble can be reconstructed by means of a modal decomposition in the eigenfunctions themselves: that can be seen as a decomposition of the originary random field into deterministic structures ( ) ( ) with random, time-dependent, coefficients.The modal amplitudes (the time series of the coefficients ( ) n a t when related to velocity fields varying with time) are uncorrelated, and their mean-square values are the eigenvalues themselves: where nm δ is the Kroneker's delta.A diagonal decomposition of the two-point velocity-correlation tensor holds: implying that the contribution of each different structure to the kinetic energy of the flow can be separately calculated, as: where E is the total kinetic energy in the domain.Thus, each eigenvalue represents the contribution of each correspondent structure ( ) n ϕ to the total amount of kinetic energy.The KL technique is optimal for modeling or reconstructing a signal in the sense that, for a given number of modes, the projection into the subspace used for modeling contains the most kinetic energy possible on average, or, equivalently, the decay of the tail of the empirical eigenspectrum is always faster (or at most as fast) than the tail of the spectrum based on any other possible basis, Fourier spectrum included.This technique has proven to give rather satisfactory results in several different cases (Alfonsi and Primavera [74][75][76][77], Alfonsi et al. [78-80], Carbone et al. [81], Vecchio et al. [82], and Zhang et al. [83]).
In the context of the problem here at hand, the high-performance computational code for the execution of KL calculations has been used, as originally developed by Alfonsi and Primavera [84,85].The above-mentioned code is general, it can be used in all kind of problems, and operates directly in physical space.More in particular, the two-point velocity-correlation tensor in Equation ( 36) is calculated in its complete form, so that the optimal representation of the velocity field outlined above, is evaluated in all the three directions (x,y,z).Within the computational procedure, Equation ( 36) is evaluated numerically by using the weight functions of the trapezoidal rule, and a scaling transformation is used to recover the symmetry in the resulting eigenvalue problem, that originally is not symmetric due to the fact that a non-uniformly spaced grid of point locations characterizes the three-dimensional computational domain involved in the decomposition operations.
For the execution of the KL decomposition (only Navier-Stokes fields are involved), the following procedure has been followed.Mainly due to limits in the computational resources available, a spatial subdomain has been selected within the whole computing domain in which the Navier-Stokes calculations have been executed.The subdomain involves ( With respect to time, the flow-field instants considered for the analysis are listed in Table 5.As a result of the decomposition procedure, a total of 3 117, 600 three-dimensional eigenfunctions (and corresponding eigenvalues) have been calculated in the subdomain.In Table 6, the individual fractions of energy as in Equation ( 40), and the cumulative energies are listed for the first ten KL modes.From the data reported in Table 6, it can be noted that the first three eigenfunctions exhibit alone a cumulative energy content of more than 99.8% of the total kinetic energy of the original velocity field.

Forces and Runups
Figure 10 reports the values of the maximum nondimensional force on the cylinder ( max nd F ) obtained in the present Euler and Navier-Stokes simulations (the reported values are related to the first wave, once the condition of regular wave train is reached), in a comparison with linear theory and with the experimental data of [14].From Figure 10 it can be noticed that linear theory slightly underestimates the wave force with respect to the results given by both our Navier-Stokes simulations and the experimental results [14].Figure 11 reports the values of the maximum nondimensional runup on the cylinder ( max nd R ) obtained in the present Euler and Navier-Stokes simulations, again in a comparison with linear theory and with the experimental data of [14].From Figure 11, it can be noted that the circumstance that linear theory strongly underestimates the real-fluid wave runup, is confirmed by the results of our Navier-Stokes simulations.This result is also valid for the experimental data obtained by Niedzwecki and Duggal [14] (exception made for a few of their values that fall under the potential-theory line).From Figures 10 and 11 it can be also noted that, in both cases of forces-and runups results, the agreement between the values obtained from the numerical integration of the Euler equations and the potential-theory results is rather good, showing that when inviscid fluids are considered, both potential theory and primitive-variable mass-and momentum-conservation (Euler) equations, give, in practice, the same results.This circumstance demonstrates the fact that all the differences that one finds between inviscid-fluid-and viscous-fluid results, have to be attributed to the effects of the real-fluid, viscosity-generated shear-stress field.

Velocity Potential
As a result of the computing procedures previously outlined at Section 3.1, results related to the ϕ − derived flow fields are presented.In Figures 12 and 13   Moreover, and as expected, the application of the swirling-strength criterion for structure extraction to the ϕ − derived velocity fields, has given no results, due to the pure mathematical nature of the latter fields.

Euler Equations
As a result of the computing procedures previously outlined (Section 3.2), results related to the Euler-equations-derived flow fields are presented.
It is to be preliminarly noted that the following figures incorporate three types of information, namely: (i) they depict the flow fields in terms of flow structures as extracted with the swirling-strength (or also ci λ ) criterion [62], and represented at given threshold values ( ci th λ , see Alfonsi and Primavera [64]), at different instants; (ii) "lines" eventually detectable in the visualizations actually represent borders between structures; (iii) the external surfaces of the flow structures are colored with the spanwise component of the vorticity ( y ω ), so that the sign and the intensity of the local-particle rotation can be inferred from the visualizations (more intense y ω coloring denotes stretching of vortex lines, less intense y ω coloring denotes shrinking of vortex lines, reddish areas denote positive y ω , bluish areas denote negative y ω , greyish areas denote zero y ω ). the Euler-derived wave field exhibits the maximum runup on the upstream external surface of the cylinder), respectively.The fields are rather regular and substantially symmetric.The flow field exhibits a number of tubular vortical structures concentrated at the free surface, where, from the distribution of these structures, one obtains the perception of the free-surface configuration in terms of wave crests, troughs, accumulation of fluid mass upstream from the cylinder, and progressive development of the wave-diffraction phenomenon (see groups of lines a to d in Figure 14).
In particular, at Overall, the nature of the above-mentioned free-surface tubular structures is that of being largely irrotational, exception made for the previously-generated vorticity contained in some of them that propagates (conservatively) across the field, according to the progression of the wave-cylinder interaction process.From Figure 14 it can be also noticed that the areas characterized by positive and/or negative, previously-generated y ω , are well defined and confined, showing that no viscous diffusion of vorticity occurs, as expected (in the absence of viscosity, no vorticity diffusion occurs).Figure 14 further shows that the wave crest-and trough-structures are substantially irrotational, while positive and/or negative vorticity is present in the ascending and/or descending portions of the wave, approximately at half-way between a crest and a trough.Moreover, it can be noted that, at one can also notice that the structures encompassed in the group of lines a exhibit an increase of their positive y ω -field, due to the turning and stretching of vortex lines, as caused by the progression of the wave-to-cylinder approaching process, while the structures encompassed by the group of lines b also increase their y ω field (of either sign), due to the turning and stretching of vortex lines nearby the cylinder.Moreover, the structures encompassed by the group of lines c increase their negative y ω -field, due to the turning of vortex lines right downstream from the cylinder, while those encompassed by the group of lines d are subjected to remarkable turning and stretching between the wave trough and wave crest, downstream from the cylinder.
In Additional structures are present.In Figures 15 and 16, arrow 4 denotes an upstream structure that represents the underwater counterpart of the accumulation of fluid mass upstream from the cylinder, that verifies when the oncoming wave approaches the cylindrical body.In Figures 16 and 17, arrow 5 denotes a downstream structure that represents the underwater remnant of the wave-cylinder interaction process.In Figures 15-17, a right-side structure and a left-side structure are also visible, mirroring the way trough which the structure field propagates from upstream to downstream from the cylinder, in the absence of viscous forces.

Navier Stokes Equations
As a result of the computing procedures previously outlined (Section 3.3), results related to the Navier-Stokes-equations-derived flow fields, are presented (as for the kind of information incorporated in the following figures, one can refer to the previous section).The fields are rather regular and substantially symmetric, and the structure field exhibits a number of tubular and-in this case-mainly flattened vortical structures, concentrated at the free surface.Again, from the distribution of these structures, one obtains the perception of the free-surface configuration in terms of wave crests and troughs, and also of the accumulation of fluid mass upstream from the cylinder, and the progressive development of the wave-diffraction phenomenon related to the wave-cylinder interaction process (lines a to q in Figure 18).The nature of the free-surface flattened (viscous) tubular structures is that of being predominantly irrotational far away from the cylinder, though becoming rotational when they actually interact with the cylindrical body.Moreover, the flattened structures tend to stay irrotational in their free-surface portion, while becoming rotational immediately under the free surface.This is due to the presence of a negligible shear-stress field at the free surface, and to a non-negligible shear-stress-field under the free surface.The structure fields shown in Figure 18 also exhibit a number of strong and not flattened tubular structures immediately under the free surface, mainly located upstream from the cylinder, that represent the underwater consequence of the manifestation of the wave troughs at the free surface in their viscous interaction with the cylindrical body.An inspection to Figure 18 also shows that, through the field, zones characterized by more and/or less intense coloring (either reddish or bluish) are present in the field, and also that, typically at the free surface, a more intense y ω coloring (either reddish or bluish) mainly verifies in the vicinity of the wave troughs.This latter phenomenon is due to the presence of the above-mentioned strongly-rotational underwater tubular structures, whose vorticity field breaks up the predominantly-irrotational character of the structures at the free-surface.Moreover, from Figure 18a, at max NS F t t = , one can notice that the presence of the cylindrical body induces an inversion in the sign of y ω in the trough immediately upstream from the cylinder (from reddish to bluish in front of the cylinder), being the latter bluish zone progressively pushed towards the boundaries of the domain at max NS R t t = (Figure 18b).What happens there is that the presence of the cylindrical body in front of the ascending portion of the wave determines an inversion of the fluid-particle orbital paths, with the consequence that the shear-stress field induces a strong negative    side (exception made for a small central area, Figure 19b), structure 4 has separated in two portions and changed the sign of y ω (Figure 19b), while on the cylinder downstream side, structure 5 still persists, though smaller than previously (Figure 21b).Additional information on the phenomenon at hand can be derived from specific structure-field representations.Figures 22 and 23 show close-ups of the structure field (block 3 of the computational domain, Figure 9), at   The following can be observed.In Figure 22a it is possible first to observe some phenomena already depicted in Figure 18a, namely the existence-upstream from the cylinder-of a spanwise tubular (not flattened) structure, as part of a more complex wave-trough flow structure, that includes the vorticity-inversion zone at center, in front of the cylinder.It also to be noted how the intensity of the vorticity field under the free-surface level is higher upstream from cylinder than downstream).Figure 22b further shows that the peripheral portion of the above-mentioned tubular structure lies more in depth (lighter-blue areas in Figure 22b near the boundaries of the computing domain) than its central part.
Figure 22 also shows the presence of structure 4 (see also at Figures 19 and 20).This structure, at max NS F t t = , is characterized by negative vorticity (Figure 22a), and by a remarkable depth just in front of the cylinder (Figure 22b).Other structures are detectable, especially in Figure 22a, namely those already characterized by lines k, l, m, n, o, p, q in Figure18a.This is the structural character of the flow field in the vicinity of the cylinder, at the time at which the wave exerts the maximum force on the obstacle.
Figure 23 shows the flow-field representation at the time at which the wave exhibits the maximum runup on the external surface of the cylindrical body.Also here (see also at Figure 18b) a complex tubular structure can be recognized upstream form the cylinder in Figure 23a, that exhibits-in its central part-a massive sub-structure just in front of the cylinder.Note how this sub-structure is characterized by less intense vorticity with respect to other portions of the same structure.From Figure 23b can be also observed that the new sub-structure, and its adjacent portions, extend remarkably in deep under the free surface, in a difference with respect to what shown in Figure 22b (at the previous instant), where the central portion of the structure upstream from the cylinder lied, in practice, at the free-surface level.Figure 23 also shows the presence of structure 4 (see also at Figures 19 and 20) as characterized-this time-by positive vorticity (Figure 23a, the structure has broken-up in two parts), and, again, by a remarkable depth under the free-surface, in front of the cylinder (Figure 23b).
Further, as a result of the computing procedures previously outlined (Section 4.2), results in terms of proper orthogonal flow modes are presented.Top views (the free-surface level) of a reconstructed velocity field, based on the three most energetic eignefunctions of the decomposition, are shown in Figures 24-28, in terms of constant wave-flow x-velocity, through a time interval equal to a wave period, in a comparison with the structure field (same computing subdomain).
Figure 24a ( 0.000 T = , T is the wave period) shows that, at the beginning of the wave period (the oncoming wave is only slightly altered by the presence of the cylinder), the free surface is mainly irrotational.On the upstream side of the cylinder, the diffraction tracks are clearly detectable in terms of borders among structures, where track a denotes the ascending portion of the wave, track b denotes the wave crest, and track e denotes the still-water level (note that, due to the particular nature of information incorporated in these figures, the track lettering is now different from that adopted in Figure 18).Only very small zones of positive (reddish) spanwise vorticity are present, immediatly upstream from the cylinder, at the cylinder wall.Overall, the field is regular and symmetric.Figure 24b shows the flow field in terms of KL modes.It can be first noted that there is no direct correspondence between the field visualized in Figure 24a and that shown in Figure 24b, meaning that the energy content of the flow distributes through the field with no direct correspondence with the topological-structure distribution.The only similarity between the two visualizations is mirrored by the position of track a, also indicated in Figure 24b.Further upstream from track a, in Figure 24b, other tracks are present (oblique arrows) that actually represent borders between flow structures with high kinetic-energy content.A further inspection to Figure 24b shows that the field upstream from the cylinder can be seen as divided in two zones, namely a first zone right in front of the obstacle, characterized by remarkably-small flow scales, and a second zone characterized by larger scales.A flow-mode border is also detectable in Figure 24b, between these two zones of high-kinetic-energy scale-organization.In Figures 25a,b the fields at 0.144 T = are represented.Figure 25a shows the evolution of the diffraction tracks with respect to the previous instant.On the upstream side of the cylinder, track d denotes the descending portion of the wave, track e is now interrupted by track f, and track g now denotes the presence in the field of the already mentioned rather complex wave trough.This latter g-structure actually mirrors the previously-outlined vorticity-inversion phenomenon in front of the cylinder (this section).Also in this case, the only detectable similarity between the visualization of Figure 25a and that of Figure 25b is the position of track d, also indicated in Figure 25b.Again, further upstream from track d, in Figure 25b, other tracks are present (oblique arrows), representing borders between flow structures with high kinetic-energy content, and a further inspection to Figure 25b shows that the field upstream from the cylinder is also divided in two zones, i.e., a first zone characterized by small flow scales, and a second zone characterized by larger flow scales.A flow-mode border also appears in Figure 25b, between these two zones of high-kinetic-energy scale-organization.In Figure 26a,b the fields at 0.500 T = are represented.Figure 26a shows again the evolution of the diffraction tracks with respect to the previous instant.On the upstream side of the cylinder, tracks f and g have in practice collapsed into a unique structure (now approaching the cylinder), tracks h and i denote the ascending portion of the wave (toward the crest), while track k indicates the presence of a new crest in the domain.Also in this case, the only detectable similarity between the visualization of Figure 26a and that of Figure 26b, is the position of track f, and, again, further upstream from track f, in Figure 26b, other tracks are present (oblique arrows), representing borders between structures with high kinetic-energy content.A further inspection to Figure 26b shows that the field upstream from the cylinder is again divided in two zones, a first small-scale zone, and a second large-scale zone, and a flow-mode border is again detectable in Figure 26b.It can be noted that the presence of the aforementioned discontinuity between smaller and larger high-kinetic-energy flow scales (and mirrored by the presence of a distinct flow-mode border) lasts from 0.000 T = up to 0.500 T = i.e., in practice, along the whole wave-to-cylinder approaching phase.In Figure 27a,b the fields at 0.750 T = are represented.In Figure 27a, on the cylinder upstream side, track i clearly mirrors the interaction between wave and cylinder, while the wave-crest track k is approaching the cylinder, assuming, at the same time, a complex configuration in front of the cylinder.In this case, the only detectable similarity between the fields represented respectively in Figure 27a,b, is track i, while, in Figure 27b, other tracks are present (oblique arrows), always representing borders between flow structures with high kinetic-energy content.Differently from the previous instants, it clearly appears from Figure 27b, that no more discontinuities between smaller and larger high-kinetic-energy flow scales are present, and the field of high-energy scales continuously evolves from the center to the periphery of the domain, in forming an elongated pattern.In Figure 28a,b the fields at 1.000 T = are represented.In Figure 28a, on the cylinder upstream side, the crest-track k is now approaching the cylinder, while tracks l and m denote the descending portion of the wave.The configuration of the high-energy scales of Figure 28b is similar to that of the previous instant, where the similarity between the two fields can be now recognized in the position of track k.
Overall, the flow fields (actually the iso-u-fields) that have been reconstructed by using the three most energetic eigenfunctions of the KL decomposition, further show that the borders of the flow structures are actually borders among high-energy flow structures of tubular nature, the latter appearing-with respect to their energy content-the most significant type of viscous-flow structures that have been found in the wave-cases considered in the present work.

Concluding Remarks
In the present work, the phenomenon of diffraction of linear waves impinging on large-diameter, surface-piercing, vertical circular cylinder has been investigated numerically, within three different theoretical frameworks of hypotheses, namely the velocity potential, the numerical integration of the primitive-variable Euler equations, and the numerical integration of the primitive-variable Navier-Stokes equations (with the DNS approach).The results obtained in terms of both global wave parameters (forces and runups), and wave-flow fields have shown that remarkable differences exist, depending on the frame of hypotheses that one eventually choses to follow, being though the numerical integration of the full Navier-Stokes equations in primitive variables (the velocity-pressure formulation) the physically-consistent approach to be adopted to obtain flow fields that correctly represent the numerical equivalent of the phenomena at hand.More in particular, as related to the calculation of forces and runups on the cylindrical body, it has been found that the agreement between the values obtained from the numerical integration of the Euler equations and the potential-theory results is rather good, showing that when inviscid fluids are considered, both potential theory and primitive-variable mass-and momentum-conservation equations, give in practice the same results.Thus, all the differences that one finds between inviscid-fluid-and viscous-fluid results (as related in particular to the maximum runup) have to be correctly attributed to the effects of the real-fluid, viscosity-generated shear-stress field.With reference to the wave-flow fields, in the case of the velocity potential approach, no flow structures appear in the field as a result of the wave-cylinder interaction process (as expected), due to the pure-mathematical nature of the latter approach.In the case of the Euler-equations-derived flow fields, flow structures actually appear, in particular contaminated by a previously-generated spanwise vorticity that conservatively evolves with time.Finally, in the case of the Navier-Stokes-equations-derived flow fields, complex flow structures of physical nature appear, as extracted with both the swirling-strength criterion, and the Karhunen-Loève decomposition technique.

Conflicts of Interest
The author declares no conflicts of interest.
in terms of KC D L − are now shown in terms of Re KC − , within the ranges 0.1 0.76 KC = − , and Re 404 2165 = − .It can be noticed that all the 50 previously-selected wave cases fall in the no-separation Zone A (implying that, in a theoretical condition of zero runup, a pure-shear type of flow establishes around the cylindrical body).

Figure 5 Figure 7
Figure 5 reports the rms values of the velocity fluctuations ( ) ( ) ( ) , , rms rms rms u v w ′ ′ ′ in wall coordinates.According to the Flow-3D solver, the position of the maximum value of the streamwise fluctuations results at 14.51 y + =, with a peak value of 2.70 (see also Mochizuki and Nieuwstadt[51], and Alfonsi[52]).Figure6reports the values of the Reynolds shear stress ( )u v ′ ′ −in wall coordinates,

Figure 9 .
Figure 9. Subdivision of computing domain in grid-point blocks (top view).
normalized inner product of the candidate structure with the field.A necessary condition for this problem is that the above-mentioned function ( ) i j x ϕ is an eigenfunction, solution of the eigenvalue problem and first-kind Fredholm integral equation: and is located immediately upstream from the cylindrical body (see at subsequent Figures).
, the velocity-potential-derived flow fields are shown, in terms of free-surface elevation in the domain, at times max F t ϕ , and max R t ϕ , respectively.

Figure 12 .
Figure 12.Velocity-potential-derived free-surface profile.Top view of computing domain, at: (a) max F t t ϕ = ; (b)

Figure 13 .
Figure 13.Velocity-potential-derived free-surface profile.Side view of computing domain, at: (a) max F t t ϕ = ; (b)

Figure
Figure 14a,b show top views of the computing domain at remarkable accumulation of fluid mass-mainly irrotational-upstream from the cylinder, while at max Euler R t t = both positive and negative vorticity contaminates the latter mass of fluid.

Figure 14 .
Figure 14.Flow structures as obtained from the solution of the primitive-variable Euler equations.Top view of free-surface profile: (a) max Euler F t t = ; (b)

Figure 15 .
Figure 15.Close-up at the cylinder surface of inviscid-flow structures as obtained from the solution of the Euler equations.Upstream view: (a) max Euler F t t = ; (b)

Figure 16 .
Figure 16.Close-up at the cylinder surface of inviscid-flow structures as obtained from the solution of the Euler equations.Side view: (a) max Euler F t t = ; (b)

Figure 17 .
Figure 17.Close-up at the cylinder surface of inviscid-flow structures as obtained from the solution of the Euler equations.Downstream view: (a) max Euler F t t = ; (b) Figure 18a,b show top views of the computing domain at max NS F t t = (the instant at which the Navier-Stokes-derived wave field exerts the maximum force on the cylinder)Navier-Stokes-derived wave field exhibits the maximum runup on the external surface of the cylindrical body), respectively.

Figure 18 .
Figure 18.Flow structures as obtained from the solution of the primitive-variable Navier-Stokes equations.Top view of free-surface profile: (a) max NS F t t = ; (b)

yω
-field in the zone right in front of the cylinder.In Figures19-21, upstream-, side-and downstream-close-up views of the structure field are shown in the vicinity of the cylinder external wall, at max From these figures it can be noticed that other types of viscous-flow structures develop under the free surface.The establishment of the latter structures actually represents the more relevant difference between the Navier-Stokes-(and/or eventually Euler-) flow-fields, and the velocity-potential-derived flow fields, so unveiling physical phenomena that are impossible to detect by only analyzing the potential-derived fields.Figures19-21show first-as also previously mentioned-the presence of two main underwater tubular, not flattened structures upstream from the cylinder and under the wave trough, characterized by high positive vorticity.In going from max in dimensions, and the one that is nearest to the cylinder develops a more extended, massive sub-structure, right in front of the cylinder itself.These structures in the whole represent the underwater consequence of the manifestation of the trough at the free surface, where the aforementioned more extended sub-structure represents the underwater counterpart of the wave-to-cylinder approaching process.

Figure 19 .
Figure 19.Close-up at the cylinder surface of viscous-flow structures as obtained from the solution of the Navier-Stokes equations.Upstream view: (a) max NS F t t = ; (b)

Figure 20 .
Figure 20.Close-up at the cylinder surface of viscous-flow structures as obtained from the solution of the Navier-Stokes equations.Side view: (a) max NS F t t = ; (b)

Figure 21 .
Figure 21.Close-up at the cylinder surface of viscous-flow structures as obtained from the solution of the Navier-Stokes equations.Downstream view: (a) max NS F t t = ; (b) of each of these figures represents a bottom close-up view of the structure field, as represented with the same criterion previously used, such that the external surfaces of the flow structures are colored with the values assumed by the spanwise component of the vorticity ( y ω , reddish areas denote positive y ω , bluish areas denote negative y ω , greyish areas denote zero y ω , more intense coloring corresponds to high y ω -values, less intense coloring corresponds to lower y ω -values).The right portions of Figures 22 and 23 also represents the structure field at max now the external surfaces of the structures are colored with the values assumed by the pressure, in a darker-to-lighter color scale (bluish areas denotes low values of the pressure, in the vicinity of the free surface, rosy areas denote higher values of the pressure, more deep under the free surface).

Figure 22 .
Figure 22.Close-up views of structures as obtained from the solution of the Navier-Stokes equations.Bottom view of computing domain at max NS F t t = .The external surfaces of structures are colored with: (a) wave-field spanwise vorticity; (b) wave-field pressure values.

Figure 23 .
Figure 23.Close-up views of structures as obtained from the solution of the Navier-Stokes equations.Bottom view of computing domain at max NS R t t = .The external surfaces of structures are colored with: (a) wave-field spanwise vorticity; (b) wave-field pressure values.

Figure 24 .
Figure 24.Flow-field representation at 0.000 T = (T is wave period), top view (free-surface level): (a) flow-structures (colored with spanwise vorticity); (b) surfaces of constant x-velocity from reconstruction based on the three most energetic KL modes.

Figure 25 .
Figure 25.Flow-field representation at 0.144 T = (T is wave period), top view (free-surface level): (a) flow-structures (colored with spanwise vorticity); (b) surfaces of constant x-velocity from reconstruction based on the three most energetic KL modes.

Figure 26 .
Figure 26.Flow-field representation at 0.500 T = (T is wave period), top view (free-surface level): (a) flow-structures (colored with spanwise vorticity); (b) surfaces of constant x-velocity from reconstruction based on the three most energetic KL modes.

Figure 27 .
Figure 27.Flow-field representation at 0.750 T = (T is wave period), top view (free-surface level): (a) flow-structures (colored with spanwise vorticity); (b) surfaces of constant x-velocity from reconstruction based on the three most energetic KL modes.

Figure 28 .
Figure 28.Flow-field representation at 1.000 T = (T is wave period), top view (free-surface level): (a) flow-structures (colored with spanwise vorticity); (b) surfaces of constant x-velocity from reconstruction based on the three most energetic KL modes.
domain along (x,y,z) in wall units (test case) , , number of grid points in the computing domain along (x,y,z) cartesian coordinates (x is wave direction, z is vertical direction) complex eigenvalue pair (swirling strength) ( )

Table 1 .
Ranges of linear-wave characteristic parameters.

Table 2 .
Characteristic parameters used for the Euler wave-flow calculations.

Table 4 .
Characteristic parameters used for the Navier-Stokes wave-flow.

Table 5 .
Flow-field instants involved in the KL decomposition (T is wave period).

Table 6 .
Energy content of the first ten eigenfuctions of the KL decomposition.
vorticity far away from the cylinder, and by strong positive vorticity nearby the cylinder (caused by the pressure field in the proximity of the cylindrical body), while the structures encompassed by the group of lines d are characterized by remarkably-less intense negative vorticity.
= , the flow structures encompassed in the group of lines a (located in a descending portion of the wave, upstream from a trough) are characterized by strong positive vorticity, the structures encompassed by the group of lines b (right after a trough) are characterized by less intense positive vorticity, those encompassed by the group of lines c (downstream from a crest) are characterized by strong negative Figures 15-17, upstream-, side-, and downstream-close-up views of the structure field are shown in the vicinity of the cylinder external wall, at = (note that, in the subsequent Navier-Stokes-field representations, structures 2 and 3 will appear as distinct).