Fate of large-scale vortices in idealized tidal lagoons

The generation and evolution of tidally-induced vortices in coastal and estuarine regions can influence water quality and sedimentary processes. These effects must be taken into consideration in the development of coastal reservoirs, barrages and lagoons, among other environmental flow applications. Results are presented here on the fate of large-scale vortices within confined tidally-forced domains. A computational approach is employed using the Thetis depth-averaged coastal ocean modeling framework. Initially, two test cases serve to demonstrate model capability in capturing the formation of dipoles downstream of oscillatory flow channels. Diagnostic quantities of vorticity and localized circulation are used to track the 2-D vortex evolution and dissipation. This approach is then applied to tidal lagoon geometries, where flows through the inlet induce a pair of counter rotating vortices (dipoles). Idealized model geometries and inlet conditions are used to determine the impact of three design parameters on large-scale vortical structures: (a) the lagoon geometry aspect ratio in the horizontal plane, (b) the inlet width and (c) the bathymetry profile as the coastline is approached. The dependence of vortex flushing behavior on the dimensionless ratio Wi/UT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{W_\text{i}}\!/_{UT}$$\end{document} (where Wi\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_\text{i}$$\end{document} is the width of the inlet channel, U is the maximum velocity and T is the tidal period) is reaffirmed, while the side walls and the sloping bathymetry are found to affect the vortex dissipation process.


Introduction
The presence of tidally-induced vortices in coastal flows is a common phenomenon. Vortices tend to form when tidal currents encounter particular geometries, such as when flowing around sharp headlands [1,2] or islands [3][4][5], and within confined estuaries [6], coastal reservoirs or harbors [7]. These vortices can become hazardous for navigation and affect the entrainment and mixing of sand (potentially creating sandbanks), pollutants or nutrients [8]. The formation of large-scale two-dimensional vortices was recently observed [9] downstream of what is currently the world's largest tidal range based power generation scheme in Lake Sihwa, South Korea. In that case, patterns were identified between the occurrence of recirculation zones and the local suspended sediment transport regime. Given how coastal reservoirs and marine energy developments including tidal lagoons are gaining interest [10][11][12], it is timely to consider the fate of such vortices in the context of confined tidal embayments.
The formation of vortices through tidal inlet flows has been investigated analytically, experimentally and computationally, e.g. [13][14][15][16]. As an example, Wells and Heijst [13] applied potential theory by representing the dipoles as a source-sink system to determine whether the vortices will be long-lived (i.e. surviving longer than a tidal cycle and gradually dissipating downstream from the tidal inlet), or whether they would return towards the inlet as the tide recedes and be eventually flushed out from the impounded area with the outgoing flow. This behavior can be examined by considering the ratio of the amount of water that escapes the influence of the sink to the amount of water that enters with the flood tide. This has been found to be dependent on the inlet flow propagation velocity, which needs to be sufficiently high for the vortices to escape from the sink flow during the ebb tide. Experimental and analytical findings show vortices escape from the inlet provided W i ∕ UT is smaller than a critical value of ∼ 0.13 [13] (where W i is the width of the inlet channel, U the maximum velocity and T the tidal period). Subsequently, experimental results on the propagation of tidally-induced vortices were reported using tracer dye injections and particle image velocimetry (PIV) that were applied to a series of physical modeling campaigns of tidal inlet geometries [14,16].
In the assessment of coastal reservoirs and lagoons, computational approaches have been reported that approximate the problem as anything from a 0-D [17], 1-D [18], 2-D [19][20][21] and 3-D [22] system. Due to the wide range of scales present in such engineering applications, hydrodynamic modeling that is capable of resolving the formation of near-field recirculation zones is feasible in only a subset of these approaches. Two-dimensional (2-D) depth-averaged modeling is often selected as it strikes a balance between accuracy and computational complexity, and allows for the consideration of certain regional and far-field effects. The presence and potential implications of large recirculation zones has been discussed in several studies [23]. For example, low-velocity flow areas at the zones' centers can promote sediment accumulation through the tealeaf paradox [24].
In this paper we consider flow within idealized tidal lagoon geometries with a focus on the generation and propagation of two-dimensional vortex hydrodynamic structures. A computational approach is employed utilizing an unstructured mesh coastal ocean model. This is motivated in part by the need to deliver an impact assessment capability that can be transferred to more realistic representations of coastal engineering infrastructure in shallow water flows. Two scenarios are initially considered to validate the model setup and demonstrate its ability to characterize the primary vortex transport.
Sequentially, simplified lagoon designs and inflow conditions are investigated to identify key parameters for the generation of long-enduring vortices.

Methodology
We consider flow simulations through tidal inlets of confined embayments performed using Thetis [25,27], a coastal ocean modeling framework implemented using the Firedrake finite element partial differential equation (PDE) solver engine [26]. Thetis features 3-D modeling capabilities [27] and has already been applied for the simulation and assessment of tidal range structure operation strategies [28]. The configuration of the model presented here is depthaveraged, assuming homogeneity in the vertical direction, and solves the non-conservative form of the 2-D shallow water equations: where is the free surface perturbation, H d = + h m is the total water depth with h m being the mean water depth and is the depth-averaged velocity vector with horizontal components u, v in the x, y directions respectively. represents kinematic viscosity in the model which is typically equal to the aggregate value of molecular and turbulent viscosities. The fluid molecular viscosity is considered negligible and a constant eddy-viscosity turbulence model is applied for simplicity in this work. The bed shear stress can be represented by the quadratic drag coefficient C D as follows [25]: or by imposing Manning's formulation as with n the Manning's coefficient ( s m −1∕3 ). Intertidal areas constitute a characteristic of most tidal estuaries with wetting and drying processes being challenging to represent numerically [29]. Thetis employs the formulation described in Kärnä et al. [30] that uses a negative depth approach where the bathymetry is allowed to move artificially conserving continuity while also being applicable to tracers. The model is implemented using a discontinuous Galerkin finite element discretization (DG-FEM), using the P 1DG -P 1DG velocity-pressure finite element pair. A semi-implicit Crank-Nicolson timestepping approach is applied for temporal discretization with a constant timestep of t . The discretized equations are solved by a Newton nonlinear solver algorithm through the PETSc library [31].
Vorticity, defined in the strong form as = − u y + v x is given in the weak form using a test function as: where is the simulation domain. Integrating by parts and applying the Gauss divergence theorem, this expands as: where represents the boundary of our domain, ds the boundary integral and n x , n y the exterior facets at the boundaries in the x and y directions respectively. A piecewise linear, continuous Galerkin finite element approximation of vorticity can now be found by solving the system of equations formed by (6) and a basis of test functions in the space of piecewise linear functions, for a function (x) in that same function space. In order to track vortices during the simulations we first identify areas of high . For each of these areas, the circulation (i.e. the velocity line integral or effectively the flux of ), denoted here c , is calculated as: where S is a closed surface around the vorticity peaks, with the area being dependent on the scale of the problem. Once a peak in the field is identified, is integrated over a circular area S centered at the peak with a radius r (7). The value of the radius r effectively dictates the size of the vortices that will be tracked and was selected according to the size of the large scale vortices that are examined. For the simulation domains considered in Sects. 3.2 and 4 a value of r = 500 m was applied. Recirculation zone centers were distinguished from shear-induced peaks through a threshold value ( | critical | ) that aims to exclude peaks that feature a high value at the core but a low value around it; such a combination renders these points unlikely to correspond to vortex cores. Determining this value was the matter of a preliminary sensitivity assessment with | critical | ∝ r∕h l where h l is a representative element length over the area S to account for resolution differences. The algorithm evaluates the proximity of each peak to vortex core coordinates from previous timesteps and either appends the new coordinates to known vortex trajectories or establishes a new trajectory.

Numerical model capabilities, validation and sensitivity
Hydrodynamic models need to be calibrated and validated to ensure that the physical processes of interest are simulated with sufficient accuracy [32]. This is typically achieved for regional-scale simulations through comparisons against observed tidal constituents and field measurements such as water elevation records (e.g. tide gauge data) or velocity samples (Acoustic Doppler Current Profiler readings). Depth-averaged numerical modeling methods have proven their accuracy in many coastal modeling applications if the assumptions underlying the shallow water equations remain valid [23,33]. In particular, as long as the coastal flow can be assumed to be independent of depth, i.e. possessing minimal three-dimensionality, depth-averaged models can efficiently and practically simulate water elevations over extended periods of time. Nonetheless, once extending the applications to marine energy and in the context of tidal barrages, certain 2-D modeling limitations have been observed nearfield to the hydraulic structures (turbines and sluice gates) where the flow becomes strongly three-dimensional and a 2-D model is thus expected to fail in capturing the true dynamics [24]. Therefore, computationally demanding modeling approaches become essential if 3-D profiles such as the flow immediately adjacent to the hydraulic structures are to be characterized [18,34,35]. According to physical modeling experiments [22], the predominantly 3-D flow regions are typically contained within a distance of 20D (where D is the turbine diameter). These near-field scales are beyond the scope and capabilities of the regional-scale models considered here. Limitations of 2-D approaches must be acknowledged in the context of vortex modeling. Tide-induced vortices can feature 3-D behavior due to vertical secondary flows. These lead to the upwelling of water from the bottom boundary layers towards the surface and through the center of the horizontal recirculation [36,37]; these effects will also contribute to the aforementioned accumulation of sediments towards the core. Another physical mechanism that might be neglected is transverse advection in the vertical direction discussed in earlier studies [38]. Its effects are more pronounced in regions where stratification is dominant. For sites that feature a tidal range of interest for marine renewable energy developments, the flows are typically well-mixed and therefore weakly, or not, stratified (e.g. in the Bristol Channel in the UK or the Rance river in France). Nonetheless, given the extended slack water times expected within tidal range structures, internal conditions would potentially merit an expansion to 3-D modeling as stratification effects might emerge. In practice, neglecting such processes could lead to an overestimation of vortex endurance [36] and any associated hydro-environmental implications.
Two test cases (one laboratory scale and one full scale) of oscillatory flow through tidal inlets are considered. Case A examines the capability of Thetis as a means to study the formation of tide-induced vortices where we make comparisons against experimental and numerical findings from Kim et al. [14]. Case B follows Wells and Heijst [39] and serves to determine whether vorticity ( ) and circulation ( c ) quantities can be used to diagnose the fate of these vortices. In addition, Case B forms the basis for a grid sensitivity assessment to highlight the dependence of the diagnostic on the localized mesh resolution and its significance for the study findings.

Case A: Oscillatory flow through a tidal inlet physical model
In the study of Kim et al. [14], a series of test cases are presented for the validation of a depth-integrated numerical model based on experimental results. In considering the behavior of oscillatory flow in a shallow water basin containing an inlet structure (denoted as Layout D in [14]), measurements and numerical results were compared. The particular setup is schematically illustrated in Fig. 1a for completeness.
Thetis was configured to replicate the recorded flow conditions beginning from equilibrium with a still water depth of 0. velocities across the mouth of the tidal inlet shown in Fig. 2, which can be seen to agree well with data reported from the experiment. The vorticity field plotted over the left basin (schematically illustrated in Fig. 1a) is in alignment with the results of Kim et al. [14] in Fig. 3. The exact magnitude and timing of the boundary forcing had to be calibrated to match the velocities through the tidal inlet, with the boundary forcing being multiplied by a ramp function in the form of r(t) = tanh((t − t 0 )∕t r ) where t 0 = 0 , and t r = T∕4 to simulate the time required for the experiment to commence from equilibrium conditions. The t r parameter was determined through preliminary calibration. Figure 3 demonstrates good performance in terms of capturing the formation and propagation of the counter-rotating vortices. A better agreement is observed in the three latter instantaneous plots of vorticity, as the simulation becomes independent of the initial equilibrium conditions. Minor deviations in the exact location of the core appear which could be attributed to differences in the meshes used between the two studies (such as a structured vs unstructured mesh based approach). However, these results provide us with confidence in the ability of the model to simulate the generation and propagation of vortices.

Case B: Dipole formation by tidal flow in a channel
The second case is based on the investigation of Wells and Heijst [13,39] into the propagation of vortices formed through idealized tidal inlet flows. Wells and Heijst [13] observed how counter-rotating vortices can form a coherent dipole structure that can self-propagate away from the channel. By treating the dipoles as a source-sink pair, their trajectory was plotted for different W i ∕ UT values and it was shown that the vortices are expected to escape from the central channel and into the confined basin if W i ∕ UT < 0.13. Numerical simulations consider the domain shown in Fig. 1b as per the case study of Wells and Heijst [39]. We initially consider two instances of bottom friction; these are represented through quadratic drag coefficients of (a) C D = 2.5 × 10 −3 and (b) C D = 1.25 × 10 −3 . In addition, we test three unstructured mesh resolutions with element edge lengths of h = 400, 200 and 100 m respectively. A structured mesh of 100 m × 100 m was employed previously [39] based on the justification that this adequately represented The tidal elevation was defined at the left hand boundary via a sinusoidal forcing of 1 m amplitude. All others boundaries, including the inlet channel, were treated with a noslip condition. A tidal period of T = 9 h was set for the tidal forcing as it led to predictions that reflect the values of the channel velocities reported in Wells and Heijst [39]. Simulations over the two coarser meshes were run with t = 50 s. On the finer mesh the timestep was reduced to 25 s to avoid instabilities related to the Courant-Friedrichs-Lewy (CFL) criterion. The simulations were ramped-up over four tidal periods, and the analysis was performed over the subsequent three periods to ensure independence from the initial equilibrium state specified for the model. The field was monitored every 200 s and processed to track vortices using Eq. 7 and the peaks (Eq. 6) across the domain.
The velocity amplitude is ≈ 1 m s −1 , corresponding to a ratio of W i ∕ UT = 0.154 which suggests that vortices will be eventually flushed out of the confined lagoon by the ebbing flow. However, it was observed [13] that bathymetry and bottom friction effects will be influential on the vortex fate even though these factors are not inherently included in the ratio evaluation. The simulation results of Fig. 4 demonstrate how with a drag coefficient of C D = 2.5 × 10 −3 the vortices dissipate rapidly in the flushing process, while with the lower C D = 1.25 × 10 −3 vortices were found to propagate further into the lagoon, enduring longer and dissipating gradually.
Values of the core and coordinates ( c , x c , y c ) were tracked in time within the confined area. Figure 4 illustrates the bottom half of the lagoon where the vortex of positive is transported. The equivalent negative vortex propagating within the top half of the monitoring domain is omitted here on the grounds of symmetry. Maximum and minimum values of are shown in Fig. 5 for the three mesh resolutions, plotted together with c . The greatest magnitude within the monitoring domain is attributed to shear-induced effects at the inlet corners and on the sides of the incoming flow at flood tide. In order to track vortices, the maxima/minima were therefore generally only monitored in the sub-region x > 16.5 km (as seen from the starting location of the trajectories plotted in Fig. 4) because this successfully identifies the values at the primary vortex core. However,  Fig. 5 is that for a lower friction coefficient the vortices dissipate significantly later, and can be traced multiple tidal cycles from their origin (see Fig. 4b). Wells and Heijst [13] argued that vortices should endure longer in deeper water if the same tidal inlet velocities are assumed and for C D = 2.5 × 10 −3 . This is to be expected given that in deeper water bed-shear effects should impact the depth-averaged flow profile to a lesser extent (as estimated by the Manning formulation in Eq. 4) and there will be a higher wave celerity C 0 = √ gH d . Practically, increasing the maximum depth ( H max ) while maintaining the seaward amplitude reduces the flow velocity magnitude by distributing the inlet flow over a greater surface area and this drastically increases W i ∕ UT .
A grid sensitivity analysis based on the the element edge lengths (h) serves to assess whether the mesh resolution is sufficient to capture the formation and evolution of vortices. It is crucial to observe that relies on the curl of the vector, which in turn requires the spatial gradient of the velocity components as expressed in the weak form in Eq. (6). Unlike the velocity field, a finer mesh resolution corresponds to higher values of as Fig. 5 demonstrates. This is not an unexpected outcome and is often seen in diagnostic quantities involving derivatives. Finer resolutions are capable of representing smaller hydrodynamic structures and therefore higher gradients, and in addition, coarser meshes potentially lead to a higher degree of implicit numerical viscosity. Notwithstanding the foregoing, the The figure values on the right are normalized with respect to U, and t using the maximum velocity at the inlet ( U max ), the peak vorticity at the core center in each case ( peak ), and the tidal period (T) of the simulations respectively. In the C D = 1.25 × 10 −3 case evidence of vortices from previous tidal cycles can be seen in the vorticity plots predictive capability of the model in terms of the vortices being dissipated or not remains consistent for all the resolutions tested, even though a coarser resolution could fail to capture detailed features beyond a certain point.
As we are primarily interested in investigating whether vortices are dissipated over the course of the tidal cycles, was normalized according to the peak values ( peak ) as shown in the right hand side of Fig. 5. This is convenient here as uniform h lengths characterize the domain of the test cases. Once normalized, results from the three mesh resolutions appear to yield very similar trends as to how in the domain gets dissipated. Fluctuations in the values monitored at the core (observed in Fig. 5) are attributed to the tracking algorithm that is based on nodal values rather than the entire field (see Fig. 4); as the resolution increases, these fluctuations decrease in magnitude.
Replicating this test case illustrates the model's capability in capturing the formation and propagation of vortices created by tidal flow. In particular, the reproduction of previous findings [39] serves to demonstrate how the new modeling approach of Thetis agrees with the results of established software (in this case the finite difference Delft3D code). The behavior of the vortices can be monitored though the combination of and c , with their fate being dependent on parameters such as the bathymetric depth H, the inlet velocity U, the channel width W i , the tidal period T and bed shear effects (represented for this case study through the quadratic drag coefficient C D ).

Application to idealized tidal lagoons
This section presents simulation results concerning the behavior of vortices created at the inlet of tidal range structures. Techniques applied in the case studies (A and B) are extended to a set of idealized lagoons. Three parameters are varied when designing the test lagoons: the aspect ratio of the lagoon sides (L / W-domain length/width), the width of the inlet ( W i ) representing the turbine housing, and the bathymetry. The aim is to demonstrate how these affect the fate of the vortices created by the inlet water jets.
Lagoons are modeled as rectangular domains with a flux boundary condition at an opening centered at the middle of the lagoon's left hand side. The opening can be perceived as a representation of the open boundary of turbine caissons [24]. The tidal lagoon simulations do not feature an inlet channel and the flow is evenly distributed across the boundary width assuming a uniform distribution of equally sized hydraulic structures that are driven by head differences, such as turbines or sluice gates. This means that curvature and vorticity quantities produced within a conventional inlet channel are not present in this case. The lagoon geometric specifications are illustrated in Fig. 6 and summarized in Table 1. These are defined in a manner that enables assessment of side wall proximity and sloping bathymetry effects, with the right hand wall of the lagoon assumed to coincide with the coastline. Five domain aspect ratios (L / W) were considered, each constrained to yield an overall rectangular area of 25 km 2 . Consequently, all options theoretically deliver the same power output when subject to identical tidal forcing downstream. As the domain width (W) increases so does the amount of the internal surface area affected by intertidal processes; in practice this will influence the power generation performance. The maximum bathymetry was set to 40 m sloping up to 0 m over a distance of 2 km on the landward (right-hand) boundary of the lagoon. This implies that the widest lagoon ( L 4∕25 ) was entirely covered by the sloping bathymetry. Three inlet width cases (of W i = 1.25 , 1.00 and 0.75 km) were successively tested for every instance of domain aspect ratio (L / W). For a tidal energy scheme P ∝ AH 2 , where A is the impounded surface area and H is the head difference between the inner and outer elevations [40]. As a result, for a given area A and under idealized conditions the same number of turbines would be installed and the remaining control factor is the distribution of the hydraulic structures over the lagoon's seaward side wall. In all simulations the same sinusoidal flux was considered to represent a broader or narrower, yet uniform, distribution of incoming flow over the inlet width. Tidal forcing was applied by specifying the flux at the inlet in the x direction and its period was set to T = 6 h starting with an ebb tide. Finally, the maximum flux at the inlet was set to 18,000 m 3 s −1 for the width of the inlet to have a significant impact on the vortex transport.
The lagoon domain meshes correspond to an element edge length of h = 50 m which is consistent with the values used in Sect. 3.2 with respect to the resolution of the inlet geometry. The time step t of the hydrodynamic simulations was set to 10 s, which satisfies the CFL criterion throughout the duration of the runs. Each case was first ramped up This behavior can indeed be observed in Fig. 7. In the case of W i = 1.25 km, changes in the vorticity field over time demonstrate how vortices do in fact dissipate as they return towards the inlet and get flushed within a single T. This is also shown by the tracking algorithm in Fig. 7c, where vortices cease to be traced until the subsequent flood tide that leads to the formation and propagation of a new recirculation zone. The vorticity intensity progressively increases as the width of the inlet is reduced, through an accelerated inlet velocity.
The results of Fig. 7 can be further complemented by the contour plots of the field that also indicate the trajectory of the vortex within the bottom half of the lagoon (Fig. 8). It should be remarked that in Fig. 8 the entire domain is included to demonstrate that despite the unstructured nature of the mesh, symmetry in terms of the vortices formed and the hydrodynamics is effectively maintained. The contour is plotted at t ≈ 3 T which corresponds to the start of an ebb tide, and also the point of maximum propagation of the dipoles that will eventually be flushed away ( W i = 1.25 km, shown in blue). In the other cases, shown in red ( W i = 1.00 km) and black ( W i = 0.75 km), vortices from the preceding tidal cycles are still observed within the lagoon interior.

Sensitivity of vortex trajectory to lagoon aspect ratio (L / W)
Thetis predictions also demonstrate that the overall lagoon shape will influence the fate of internal vortices. For L 25∕4 , the longest domain, vortices tend to dissipate within a single period. This is shown by the drop in vorticity that happens at approximately t = 0.35 T after the identification of the vortex, which appears to be consistent for every case of W i (Fig. 7a). Nonetheless, there are still differences observed for different inlet widths; a wider inflow distribution corresponds to an earlier drop in the vorticity (see Fig. 7a). Visualizing and the c trajectory it becomes apparent that this dissipation is caused by the presence of the (no-slip) side walls close to the inlet (see the blue trajectory in Fig. 9).
Results for L 25∕16 are plotted in Fig. 7b. Unlike L 25∕4 , vortices dissipate only when the inlet width is equal to 1.25 km indicating that their behavior is governed by the dimensionless parameter W i ∕ UT similarly to the square case L 1 , rather than the proximity of the side walls. Indeed for W i = 1.25 km the main vortex structures get flushed within 0.55 T. Moreover, vortices for L 25/16 tend to endure substantially longer than in all other L / W configurations tested. Initially unhindered by side wall effects, the vortices propagate close to the centerline of the lagoon until they are redirected and start to dissipate through the influence of the shallower bathymetry towards the coastal boundary (see the orange trajectory in Fig. 9).
The transient field for L 16∕25 shows that vortices dissipate through a combination of sloping bathymetry and proximity of the inlet to the coastal side of the lagoon. This leads to further mixing, where the dipoles move back towards the inlet and even merge with the next pair of vortices (see the red trajectories in Fig. 9). Figure 7d shows the intermediate L 16∕25 case where for the 0.75 and 1.00 km-wide inlets, trends observed in L 1 are closely followed. For L 4∕25 , vortices tend to remain closer to the inlet, a behavior attributed to the  increased influence of the friction and reduced wave celerity in shallower waters, largely due to the effect of the sloping bathymetry that influences the entire lagoon. The flushing ratio ( W i ∕ UT ), the maximum peak and the average lifetime of the vortices are summarized in Table 2 for completeness.

Sensitivity of vortex trajectory to lagoon bathymetry
Variations in bathymetry can be significant in coastal applications and must be acknowledged for their influence. In Table 2 an interesting pattern emerges with regard to the values of the vortex lifetime t l for W i = 1.00 and W i = 0.75 km. Even though the maximum peak vorticity is consistently encountered for W i = 0.75 km, the associated vortices in L 25∕16 , L 1 and L 16∕25 survive for less time than the W i = 1.00 km case. As a consequence of the greater momentum of the inflow associated with the narrower inlet (for W i = 0.75 km), vortices are advected closer to the shallower section of the lagoon and along the coast. Increased effects of friction become more pronounced and expedite the dissipation of vorticity, eventually contributing to a reduced lifetime t l relative to the W i = 1.00 km case.
We consider an additional case to demonstrate the implications of bathymetry on the fate of the vortices formed by the inlet flow. In this case the bathymetry of lagoon L 1 is represented by a gradient from 40 m at the inlet of the lagoon to 0 m at the coast; hence the sloping bathymetry covers 100 % of the domain instead of 40% from previously. Results obtained with an inlet width of 0.75 km are compared against the original L 1 in Figs. 10 and 11. As suggested by Table 2 in the original L 1 , the vortex lifetime t l is approximately 3.13 T. In the case of the 100% sloping bathymetry, this value reduces to 1.71 T. The trajectory is also significantly altered, with the vortices encountering more resistance as they move closer to the coast.

3
These predictions can be better understood by observing the bed shear stress ( b ) parametrization in Eq. (4). As the depth ( H d ) reduces with the sloping bathymetry, the b term that determines the energy dissipation rate increases. In turn, the depth-averaged velocity magnitude | | also increases, contributing further to friction effects based again on the same Eq. (4). In addition, the sloping bathymetry leads to differences in | | within the vortex. This leads to an uneven rotation rate that diverts the trajectory of the core depending on its clockwise/counter-clockwise rotation. This asymmetry with regard to the depth-averaged rotational velocities corresponds to vortex stretching effects in the propagation to shallower waters and an equivalent contraction as vortices are diverted to deeper waters. In Fig. 11, and consistently with Fig. 8, the entire domain is included illustrating how the trajectories of the counter-rotating vortices are diverted in opposing directions at an earlier stage of their evolution.
Consequently, the sloping bathymetry can have a significant impact on the dissipation and fate of the vortices. Indeed, both cases in Fig. 10 reach comparable peak vorticity

Discussion
The idealized lagoon representation considered in this work has highlighted certain parameters which influence the trajectory and lifetime of large-scale vortices formed by inlet water jets. Using the field and the | c | quantity around the peaks it has been demonstrated how the trajectory and lifetime of vortices could be monitored over time. Beginning from the fundamental study of Wells and Heijst [13,39], it was shown how the estuarine/lagoon geometry surrounding tidal inlets will play a significant role on the fate of vortices in practical cases. In the context of tidal range structure designs, several assumptions were made. For example, the distribution of the hydraulic structures does not have to be in a single location as was the case in all the scenarios considered here.
In fact, the operation of tidal power plants involves the coordination of turbines and sluice gates. These operate at different times of the tidal cycle and therefore may act to exacerbate or dissipate vortices within the lagoons depending on their configuration [23]. There are also a wide range of control sequences for regulating case-specific inflows associated with the power generation strategy. For example, the La Rance barrage in St Malo (France) features significant holding periods at ebb tide to facilitate a sufficient head difference between the inner and outer water levels with a view to generating power efficiently through bulb turbines as the tide recedes. However, these strategies could be altered in time subject to spring-neap tidal conditions [28] and this might lead to significantly variable flow patterns and internal vortex trajectories.
The inlet flow for our idealized cases was assumed to be sinusoidal. In practice, the flow will vary based upon the characteristics of the turbine and sluice gate structures as well as the transient tidal conditions. Improved numerical modeling should therefore feature appropriate parameterizations that reflect the technology installed in potential tidal range structure schemes. Regarding the ensuing effects of these vortices, the large-scale vortex structures predicted here will not necessarily have a noticeable impact on the power generation function of tidal (or more generally hydro-) power plants, as reported in depth-averaged model assessments [24]. The predictions are primarily of interest in terms of their influence on hydro-environmental processes such as sediment transport and scalar mixing. A better understanding of the links that arise between these hydrodynamic structures and water quality or siltation could shed light on the potential long-term environmental effects of such engineering developments. Furthermore, bed geomorphology will lead to spatially varying bed friction, with irregular features such as shallow reefs and headlands (e.g. [9]) dictating local hydrodynamics and by extension the fate of vortices. Given, these assumptions, this paper has focused on presenting a method for determining the fate of vortices formed in close proximity to coastal and offshore developments. This methodology can now be adapted to the assessment of these structures once the practicalities raised in this discussion (which are certainly not exhaustive) are acknowledged using the appropriate numerical methods that have been developed to represent the operation of realistic tidal range structures in coastal and estuarine environments [11,24]. Going forward, it would be interesting to identify correlations between vortex fate, nearfield morphodynamics and ecological implications induced by the introduction of tidal inlets and coastal reservoirs.

Conclusions
This paper has focused on the development, validation and application of numerical tools to identify vortex hydrodynamic structures associated with flows through tidal inlets and coastal reservoirs/lagoons. Such flow features can be encountered in a range of coastal and estuarine conditions; a particular application discussed here is the design of tidal range energy projects. Understanding the implications of long-enduring vortices adjacent to hydraulic structures has been identified as one of the hydrodynamic challenges that needs to be addressed in evaluating the environmental impacts of marine or hydropower infrastructure. Using two test cases, we have demonstrated that the coastal ocean modeling solver, Thetis, is an appropriate tool to simulate the flow through tidal inlets and tidal lagoon inflows in a manner that predicts the formation of large, jetinduced recirculation zones.
We discussed how the non-dimensional ratio W i ∕ UT can be a useful indicator as to whether vortices will propagate further into confined lagoons or be flushed as the oscillatory flow ebbs (or changes direction depending on the side of the tidal inlet we are interested in). However, through a set of idealized scenarios we highlight the significance of bed friction and bathymetry profile, and indicate how they can influence both the trajectory and the dissipation rate of the vortices.
A simulation methodology was presented that traces moving vortex structures by identifying peaks in the vorticity field ( ) and quantifying the circulation ( | c | ) to assess whether these points in time correspond to a distinct vortex trajectory. These methods were used to consider the effect of several design parameters for idealized tidal lagoons or coastal reservoirs. For instance, the inlet width was shown to impact significantly upon the flushing and dissipation behavior of the vortices, with a wider distribution of the inlet flow leading to a greater W i ∕ UT value and an inflow of lower momentum that leads to an expedited dissipation of the primary vortex structures. Design characteristics such as the aspect ratio of the lagoons (L / W) were also shown to impact the dissipation of vortices as well as the propagation trajectory over time. Rapid vortex dissipation was also promoted by the side walls or the proximity of the inlet to the coast, particularly for very large or small lagoon aspect ratios. Finally, it is proposed that the methodology for assessing the fate of vortices should be employed to model more realistic scenarios that consider further technical and operational aspects associated with the design of coastal structures, such as tidal range power generation schemes.