Modeling wave-mud interaction on the central chenier-plain coast,western Louisiana Shelf,USA

The strong coupling between hydrodynamics and seafloors on shallow muddy shelves, and resulting bed reworking, have been extensively documented. On these shelves, spectral wave transformation is driven by a complex combination of forcing mechanisms that include nonlinear wave interactions and wave energy dissipation induced by fluid-mud at a range of frequencies. Wave-mud interaction is investigated herein by using a previously validated nonlinear spectral wave model and observations of waves and near-bed conditions on a mildly-sloping seafloor off the muddy central chenier-plain coast, western Louisiana Shelf, United States. Measurements were made along a cross-shelf transect spanning 1 km between 4 and 3 m water depths. The high-resolution observations of waves and near-bed conditions suggest presence of a fluid mud layer with thickness sometimes exceeding 10 cm under strong long wave action (1 meter wave height with 7 s peak period at 4 meter depth). Spectral wave transformation is modeled using the stochastic formulation of the nonlinear Mild Slope Equation, modified to account for wave-breaking and mud-induced dissipation. The model is used in an inverse manner in order to estimate the viscosity of the fluid mud layer, which is a key parameter controlling mud-induced wave dissipation but complicated to measure in the field during major wave events. Estimated kinematic viscosities vary between 10⁻⁴-10⁻³ m²/s. Combining these results of the wave model simulations with in-depth analysis of near-bed conditions and boundary layer modeling allows for a detailed investigation of the interaction of nonlinear wave propagation and mud characteristics. The results indicate that mud-induced dissipation is most efficient when the wave-induced resuspensions of concentrations  > 10 g/L settle due to relatively small bottom stresses to form a fluid mud layer that is not as thin and viscous as a consolidated seafloor in absence of wave action but also not as thick and soft as a near-bed high concentration layer that forms during strong wave action.     

The influence of sea state on formation speed of alongshore variability in surf zone sand bars

The formation time of alongshore morphological variability in surf zone sand bars has long been known to differ from one beach to the other and from one post-storm period to another. Here we investigate whether the type of sea state, i.e. distant swell waves or locally generated short period wind sea, affects the formation time of the emerging alongshore topographic variability.A numerical modeling approach is used to examine the emergence of alongshore variability under different shore-normal wave forcing. A research version of Delft3D, operating on the time-scale of wave groups, is applied to a schematised bathymetry with a single bar. The model is then used to investigate several wave scenarios, examining the impact of peak period, frequency spread and directional spread on the formation time of alongshore variability.Results show that an increase in wave period has a large effect, changing the formation time up to O (250%) in case the wave period is changed from a representative value for the Dutch coast (T_p ~ 5–6 s) to an Australian South East coast value (T_p ~ 10–12 s). In contrast, modifications in the directional and frequency spread of the wave field result only in a minor change in the formation time.Examination of hydrodynamics and potential sediment transport shows that the variations in formation time are primarily related to changes in the magnitude of the time-averaged flow conditions. Variations in the magnitude of very low frequency (f < 0.004 Hz) or infragravity (0.004 < f < 0.04 Hz) surf zone flow velocities do not affect the mean sediment transport capacity. Consequently the formation speed of patterns is primarily governed by positive feedback between mean flow and morphology, and low frequency flow fluctuations are of minor importance.These findings indicate that the development of alongshore topographic variability may be faster at swell dominated open coasts, primarily due to the occurrence of longer period swell. Also, at a given site, the arrival of a long wave period swell after a storm can accelerate the emergence of variability.     

A numerical study of the effect of hurricane wind asymmetry on storm surge and inundation

The influence of the asymmetric structure of hurricane wind field on storm surge is studied with five types of numerical experiments using a three-dimensional storm surge model. The results from the case of Hurricane Floyd (1999) show that Floyd-induced peak surge would have been much higher had the region of maximum wind rotated 30–90° counterclockwise. The idealized cases (the hypothetical hurricanes) with a wind speed asymmetry of 20 m s^?1 show that the peak (negative) surge varied from 4.7 to 6.0 m (?5 to ?5.7 m) or equivalent to ?8.8% and 16.3% (2.8% and ?10.4%) differences as compared to the control experiment. The area of flooding varied from 3552 to 3660 km². The results from two other idealized cases of varying degree of wind speed asymmetry further show that with decreasing (increasing) asymmetry of wind fields, the variations of peak surge and peak negative surge caused by the rotation of wind fields decrease (increase) accordingly. The results suggest that in storm surge simulations forced by winds derived from balanced models, considerable uncertainty in storm surge and inundation can result from wind asymmetries. This is true even if all other storm parameters, including maximum wind speed, the radius of maximum winds (storm size), minimum central pressure, storm translation speed, drag coefficient, and model settings (domain size and resolution) are identical. Thus, when constructing ensemble and probabilistic storm surge forecasts, uncertainty in wind asymmetry should be considered in conjunction with variations in storm track, storm intensity and size.     

Probability distribution of wave heights attenuated by salt marsh vegetation during tropical cyclone

Zero-crossing wave heights, obtained from the field measurement of random waves propagating through salt marsh vegetation (Spartina alterniflora) during a tropical storm, were analyzed to examine their probability distribution. Wave data (significant wave heights up to 0.4 m in 0.8 m depth) were collected over a two-day period along a 28 m transect using three pressure transducers sampling at 10 Hz. Wave height distribution was observed to deviate from the Rayleigh distribution. The observed probability densities of the larger wave heights were reduced significantly by vegetation, producing wave heights lower than those predicted by the Rayleigh distribution. Assuming Rayleigh distributed wave heights for the incident waves to the vegetation patch, existing vegetation-induced wave attenuation formulations are used to derive a special form of two-parameter Weibull distribution for wave heights in the inundated wetland. The scale parameter of the distribution is theoretically shown to be a function of the shape parameter, which agrees with the measurements, effectively reducing the proposed distribution to a one-parameter type. The derived distribution depends on the local parameters only and fits well to the observed distribution of wave heights attenuated by vegetation. Empirical relationships are developed to estimate the shape parameter from the local wave parameters.     

Offshore wind farm impacts on surface waves and circulation in Eastern Lake Ontario

A coupled wave and hydrodynamic model was applied to the Kingston Basin of eastern Lake Ontario, a region with bathymetric variability due to channels and shoals, to assess the potential impacts on surface waves and wind-driven circulation of an offshore wind farm. The model was used to simulate a series of storm events with time-varying wind forcing and validated against wave, current and water level observations. The wind farm was simulated by adding semi-permeable structures in the surface wave model to represent the turbine monopiles, and by adding an energy loss term to the fluid momentum equations in the hydrodynamic model to represent the added drag of the monopiles on the flow. The results suggest that the wind farm would have a small influence on waves and circulation throughout the wind farm area, with spatial variability due to focussing of wave energy and re-direction of the flow. Overall, the results indicate that the wave height in coastal areas will be minimally affected with changes in significant wave height predicted to be < 3%. Larger changes to the strength of circulation occur inside the wind farm region with localized changes in current magnitude of up to 8 cm s^− 1. The results of this study may help to understand the impacts of future offshore wind farms and other offshore structures in the Great Lakes.     

Coastal defence through wave farms

The possibility of using wave farms for coastal defence warrants investigation because wave energy is poised to become a major renewable in many countries over the next decades. The fundamental question in this regard is whether a wave farm can be used to reduce beach erosion under storm conditions. If the answer to this question is positive, then a wave farm can have coastal defence as a subsidiary function, in addition to its primary role of producing carbon-free energy. The objective of this work is to address this question by comparing the response of a beach in the face of a storm in two scenarios: with and without the wave farm. For this comparison a set of ad hoc impact indicators is developed: the bed level impact (BLI), beach face eroded area (FEA), non-dimensional erosion reduction (NER), and mean cumulative eroded area (CEA); and their values are determined by means of two coupled models: a high-resolution wave propagation model (SWAN) and a coastal processes model (XBeach). The study is conducted through a case study: Perranporth Beach (UK). Backed by a well-developed dune system, Perranporth has a bar between − 5 m and − 10 m. The results show that the wave farm reduces the eroded volume by as much as 50% and thus contributes effectively to coastal protection. This synergy between marine renewable energy and coastal defence may well contribute to improving the viability of wave farms through savings in conventional coastal protection.     

Wave climate scatter performance of a cycloidal wave energy converter

A lift based cycloidal wave energy converter (WEC) was investigated using potential flow numerical simulations in combination with viscous loss estimates based on published hydrofoil data. This type of wave energy converter consists of a shaft with one or more hydrofoils attached eccentrically at a radius. The main shaft is aligned parallel to the wave crests and submerged at a fixed depth. The operation of the WEC as a wave-to-shaft energy converter interacting with straight crested waves was estimated for an actual ocean wave climate. The climate chosen was the climate recorded by a buoy off the north-east shore of Oahu/Hawaii, which was a typical moderate wave climate featuring an average annual wave power P_W = 17 kWh/m of wave crest. The impact of the design variables radius, chord, span and maximum generator power on the average annual shaft energy yield, capacity factor and power production time fraction were explored. In the selected wave climate, a radius R = 5 m, chord C = 5 m and span of S = 60 m along with a maximum generator power of P_G = 1.25 MW were found to be optimal in terms of annual shaft energy yield. At the design point, the CycWEC achieved a wave-to-shaft power efficiency of 70%. In the annual average, 40% of the incoming wave energy was converted to shaft energy, and a capacity factor of 42% was achieved. These numbers exceeded the typical performance of competing renewables like wind power, and demonstrated that the WEC was able to convert wave energy to shaft energy efficiently for a range of wave periods and wave heights as encountered in a typical wave climate.     

Force and flow characteristics of a circular cylinder with uniform surface roughness at subcritical Reynolds numbers

The main purpose of this study is to establish a better understanding of the relationship between drag reduction and surface roughness. Experiments were conducted to measure the force and flow characteristics of a circular cylinder with different types of artificial surface roughness over the range 6 × 10³ < Re < 8 × 10⁴ (Re is based on the cylinder diameter D). The roughness cylinder was formed by covering the exterior surface of the cylinder with uniformly distributed (1) sandpaper, (2) netting, and (3) dimples. The roughness coefficient ranged from k/D = 0.0028 to 0.025 (k is the roughness height). A detailed quantitative measurement of the flow field around the cylinder using Particle Imaging Velocimetry (PIV) was carried out. The hydrodynamic force coefficients (drag and lift) of the rough cylinders are compared against those of a smooth cylinder measured under the same flow conditions. It is found that certain configuration of surface roughness significantly reduces the mean drag coefficient of the cylinder, particularly at large Reynolds numbers. In addition, the root-mean-square (r.m.s.) lift coefficient of the rough cylinders is considerably lower than that of a smooth cylinder.     

Numerical study of turbulence and wave damping induced by vegetation canopies

Vegetation canopies control mean and turbulent flow structure as well as surface wave processes in coastal regions. A non-hydrostatic RANS model based on NHWAVE (Ma et al., 2012) is developed to study turbulent mixing, surface wave attenuation and nearshore circulation induced by vegetation. A nonlinear k − ϵ model accounting for vegetation-induced turbulence production is implemented to study turbulent flow within the vegetation field. The model is calibrated and validated using experimental data from vegetated open channel flow, as well as nonbreaking and breaking random wave propagation in vegetation fields. It is found that the drag-related coefficients in the k − ϵ model C_fk and C_fϵ can greatly affect turbulent flow structure, but seldom change the wave attenuation rate. The bulk drag coefficient C_D is the major parameter controlling surface wave damping by vegetation canopies. Using the empirical formula of Mendez and Losada (2004), the present model provides accurate predictions of vegetation-induced wave energy dissipation. Wave propagation through a finite patch of vegetation in the surf zone is investigated as well. It is found that the presence of a finite patch of vegetation may generate strong pressure-driven nearshore currents, with an onshore mean flow in the unvegetated zone and an offshore return flow in the vegetated zone.     

Total kinetic energy in four global eddying ocean circulation models and over 5000 current meter records

We compare the total kinetic energy (TKE) in four global eddying ocean circulation simulations with a global dataset of over 5000, quality controlled, moored current meter records. At individual mooring sites, there was considerable scatter between models and observations that was greater than estimated statistical uncertainty. Averaging over all current meter records in various depth ranges, all four models had mean TKE within a factor of two of observations above 3500 m, and within a factor of three below 3500 m. With the exception of observations between 20 and 100 m, the models tended to straddle the observations. However, individual models had clear biases. The free running (no data assimilation) model biases were largest below 2000 m. Idealized simulations revealed that the parameterized bottom boundary layer tidal currents were not likely the source of the problem, but that reducing quadratic bottom drag coefficient may improve the fit with deep observations. Data assimilation clearly improved the model-observation comparison, especially below 2000 m, despite assimilated data existing mostly above this depth and only south of 47 °N. Different diagnostics revealed different aspects of the comparison, though in general the models appeared to be in an eddying-regime with TKE that compared reasonably well with observations.     

Solitary wave generation dynamics at Luzon Strait

A high resolution modeling study is undertaken, with a 2.5-dimensional nonhydrostatic model, of the generation of internal waves induced by tidal motion over the ridges in Luzon Strait. The model is forced by the barotropic tidal components K1, M2, and O1. These tidal components, along with the initial density field, were extracted from data and models. As the barotropic tide moves over the Luzon Strait sills, there is a conversion of barotropic tidal energy into baroclinic tidal energy. Depressions are generated that propagate towards the Asian Seas International Acoustics Experiment (ASIAEX) test site on the Chinese continental shelf. Nonlinear effects steepen the depressions, frequency and amplitude dispersion set in, and disintegration into large amplitude solitary waves occurs. The effects of varying the initial density field, tidal component magnitudes, as well as adding a steady background current to represent the occasional excursions of the Kuroshio Current into the strait, are considered.Depressions are generated at each of the two sills in Luzon Strait which radiate away, steepening and evolving into internal solitary wave trains. Baroclinic fluxes of available potential energy, kinetic energy and linear are calculated for various parameter combinations. The solitary wave trains produced in the simulations generally consist of large amplitude wave trains alternating with small amplitude wave trains. During strong tidal flow, Kelvin–Helmholtz type instabilities can develop over the taller double-humped sill. The solitary waves propagating towards the ASIAEX test site have been observed to reach amplitudes of 120–250 m, depending on the tidal strength. ASIAEX observations indicate amplitudes up to 150 m and the Windy Island Experiment (WISE) measurements contain magnitudes over 200 m. The model results yield solitary wave amplitudes of 70–300 m and half widths of 0.60–3.25 km, depending on parameter values. These are in the range of observations. Measurements by Klymak et al. (2006), in the South China Sea, exhibit amplitudes of 170 m, half widths of 3 km and phase speeds of 2.9 m s^?1. Model predictions indicate that the solitary waves making up the wave packet each experience different background currents with strong near surface shear.The energy in the leading soliton of the large amplitude wave trains ranges between 1.8 and 9.0 GJ m^?1. The smaller value, produced using barotropic tidal currents based on the Oregon State University data base, is the same as the energy estimated to be in a solitary wave observed by Klymak et al. (2006). Estimates of the conversion of barotropic tidal energy into radiating internal wave energy yield conversion rates ranging between 3.6% and 8.3%.     

Wave attenuation in mangroves: A quantitative approach to field observations

Coastal mangroves, dwelling at the interface between land and sea, provide an important contribution to reducing risk from coastal hazards by attenuating incident waves and by trapping and stabilizing sediments. This paper focusses on relations between vegetation densities, wave attenuation rates, sediment characteristics and sedimentation rates in mangroves. These processes were studied along two cross-shore transects through mangroves fringing estuaries in the southern Andaman region of Thailand. Volumetric vegetation densities in these mangroves were ranging up to 32‰, depending on the water depth. Generalized total wave attenuation rates increased from 0.002 m^− 1 in the sparsely vegetated forest fringes with Avicennia and Sonneratia species, up to 0.012 m^− 1 in the dense Rhizophora vegetation in the back of the forests. The total wave attenuation rates integrate effects of shoaling and energy losses due to various bio-physical interactions within the mangrove ecosystem. Wave attenuation in the mangroves is presumably dominated by energy losses due to vegetation drag, since wave attenuation due to bottom friction and viscous dissipation on the bare mudflats is significantly lower than those inside the mangrove vegetation.Additionally, wave attenuation in the mangroves was found to facilitate enhanced net sediment deposition and a gradual fining of the bed material. These findings corroborate the coastal defence function of mangroves by quantifying their contribution to wave attenuation and sediment trapping. The explicit linking of these properties to vegetation composition and structure facilitates modelling studies investigating the mechanisms determining the coastal defence capacities of mangroves.     

Quantifying the potential impact of land cover changes due to sea-level rise on storm surge on lower Texas coast bays

In this study we investigated the impacts of potential changes of land cover due to sea-level rise (SLR) on storm surge (i.e., the rise of water above normal sea level, namely mean-sea level and the astronomical tide, caused by hurricane winds and pressure) response inside bays on the lower Texas coast. We applied a hydrodynamic and wave model (ADCIRC + SWAN) forced by hurricane wind and pressure fields to quantify the importance of SLR-induced land cover changes, considering its impacts by changing bottom friction and the transfer of wind momentum to the water column, on the peak surge inside coastal bays. The SLR increments considered, 0.5 m to 2.0 m, significantly impacted the surge response inside the bays. The contribution of land cover changes due to SLR to the surge response, on average, ranged from a mean surge increase of 2% (SLR of 0.5 m) to 15% (SLR of 2.0 m), in addition to the SLR increments. The increase in surge response strongly depended on storm condition, with larger increases for more intense storms, and geographical location. Although land cover changes had little impact on the surge increase for SLR increments lower than 1.0 m, intense storms resulted in surge increase of up to 10% even for SLR below 1.0 m, but in most cases, the geometry changes were the major factor impacting the surge response due to SLR. We also found a strong relationship between changes in bottom friction and the surge response intensification; demonstrating the importance of considering land cover changes in coastal regions that are highly susceptible to SLR when planning for climate change.     

Influence of surface forcing on near-surface and mixing layer turbulence in the tropical Indian Ocean

An autonomous upwardly-moving microstructure profiler was used to collect measurements of the rate of dissipation of turbulent kinetic energy (ε) in the tropical Indian Ocean during a single diurnal cycle, from about 50 m depth to the sea surface. This dataset is one of only a few to resolve upper ocean ε over a diurnal cycle from below the active mixing layer up to the air–sea interface. Wind speed was weak with an average value of ~5 m s⁻¹ and the wave field was swell-dominated. Within the wind and wave affected surface layer (WWSL), ε values were on the order of 10⁻⁷–10⁻⁶ W kg⁻¹ at a depth of 0.75 m and when averaged, were almost a factor of two above classical law of the wall theory, possibly indicative of an additional source of energy from the wave field. Below this depth, ε values were closer to wall layer scaling, suggesting that the work of the Reynolds stress on the wind-induced vertical shear was the major source of turbulence within this layer. No evidence of persistent elevated near-surface ε characteristic of wave-breaking conditions was found. Profiles collected during night-time displayed relatively constant ε values at depths between the WWSL and the base of the mixing layer, characteristic of mixing by convective overturning. Within the remnant layer, depth-averaged values of ε started decaying exponentially with an e-folding time of 47 min, about 30 min after the reversal of the total surface net heat flux from oceanic loss to gain.     

Enhancement of particulate organic carbon export flux induced by atmospheric forcing in the subtropical oligotrophic northwest Pacific Ocean

The effects of extreme atmospheric forcing on the export flux of particulate organic carbon (POC) in the warm oligotrophic nitrogen-limited northwest Pacific Ocean were examined in 2007 during the spring Asian dust storm period. Several strong northeast monsoon events (maximum sustained wind speeds approaching 16.7 m s^? 1, and gusts up to 19.0 m s^? 1) accompanied by dust storms occurred during a 1-month period. The cold stormy events decreased surface water temperature and induced strong wind-driven vertical mixing of the water column, resulting in nutrient entrainment into the mixed layer from subsurface waters. As a result, the export flux of POC ranged from 49 to 98 (average value = 71 ± 16) mg m^? 2 day^? 1, approximately 2–3 times greater than average values in other seasons. As dry and wet deposition of nitrogen attributable to Asian dust storm events does not account for the associated increases in POC stocks in this N-limited oligotrophic oceanic region, the enhancement of POC flux must have been caused by nutrient entrainment from subsurface waters because of the high winds accompanying the dust storm events.     

Biogeochemical responses to late-winter storms in the Sargasso Sea,II: Increased rates of biogenic silica production and export

Previous studies measuring biogenic silica production in the Sargasso Sea, all conducted when no phytoplankton bloom was in progress, have reported a mean rate of 0.4 mmol Si m^?2 d^?1 and maximum rate of 0.9 mmol Si m^?2 d^?1, the lowest rates yet recorded in any ocean habitat. During February/March of 2004 and 2005 we studied the effects of late-winter storms prior to seasonal stratification on the production rate, standing stock and vertical export of biogenic silica in the Sargasso Sea. In 2004, alternating storm and stratification events provided pulsed input of nutrients to the euphotic zone. In contrast, nearly constant storm conditions in 2005 caused the mixed layer to deepen to ～350 m toward the end of the cruise. Biogenic silica production rates in the upper 140 m were statistically indistinguishable between years, averaging ～1.0 mmol Si m^?2 d^?1. In early March 2004, a storm event entrained nutrients into the euphotic zone and, upon stabilization, vertically integrated biogenic silica in the upper 140 m nearly doubled in 2 days. Within 4 days, 75–100% of the accumulated biogenic silica was exported, sustaining a flux to 200 m of ～0.5 mmol Si m^?2 d^?1 (4× greater than export measured during February and March in the mid-1990s). In 2005, destabilization without stratification increased biogenic silica flux at 200 m up to two-fold above previously measured export in late winter, with little or no increase in water-column biogenic silica. Despite comprising <5% of total chlorophyll, diatoms accounted for an estimated 25–50% of the nitrate uptake in the upper 140 m and 35–97% of the particulate organic nitrogen export from the upper 200 m during both cruise periods. These previously unobserved brief episodes of diatom production and export in response to late-winter storms increase the estimated production and export of diatom-derived material in the Sargasso Sea in late winter by >150%, and increase estimated annual biogenic silica production in this region by ～8%.     

Simulation of storm surge,wave, and coastal inundation in the Northeastern Gulf of Mexico region during Hurricane Ivan in 2004

Hurricane-induced storm surge, waves, and coastal inundation in the northeastern Gulf of Mexico region during Hurricane Ivan in 2004 are simulated using a fine grid coastal surge model CH3D (Curvilinear-grid Hydrodynamics in 3D) coupled to a coastal wave model SWAN, with open boundary conditions provided by a basin-scale surge model ADCIRC (Advanced CIRCulation) and a basin-scale wave model WW3 (WaveWatch-III). The H1wind, a reanalysis 10-m wind produced by the NOAA/AOML Hurricane Research Division (HRD), and a relatively simple analytical wind model are used, incorporating the effect of land dissipation on hurricane wind. Detailed comparison shows good agreement between the simulated and measured wind, waves, surge, and high water marks. Coastal storm surge along the coast is around 2–3 m, while peak surge on the order of 3.5 m is found near Pensacola, which is slightly to the east of the landfall location on Dauphin Island. Wind waves reach 20 m at the Mobile South station (National Data Buoy Center buoy 42040) on the shelf and 2 m inside the Pensacola/Escambia Bay. Model results show that wave-induced surge (total surge subtracted by the meteorologically-induced surge due to wind and pressure) accounts for 20–30% of the peak surge, while errors of the simulated surge and waves are generally within 10% of measured data. The extent of the simulated inundation region is increased when the effects of waves are included. Surge elevations simulated by the 3D model are generally up to 15% higher than that by the 2D model, and the effects of waves are more pronounced in the 3D results. The 3D model results inside the Pensacola/Escambia Bay show significant vertical variation in the horizontal currents. While the estuary has little impact on the surge elevation along the open coastal water, surge at the head of Escambia Bay is more than 50% higher than that at the open coast with 1.5 h delay.     

Hydrodynamic performance of double rows of piles suspending horizontal c shaped bars

The wave transmission, reflection, and energy dissipation of the double rows of vertical piles suspending horizontal steel C shaped bars are experimentally and theoretically studied under normal regular waves. Different wave and structural parameters are investigated e.g. the wave length, the C shaped bars draft and spacing, the supporting piles diameter and spacing, and the space between the double rows. Also, the theoretical model based on an eigenfunction expansion method is developed to study the hydrodynamic breakwater performance. In order to examine the validity of the theoretical model, the theoretical results are compared with the experimental and theoretical results obtained by different authors. Comparison between experiments and predictions showed that theoretical model provides a good estimate to the different hydrodynamic coefficients when the friction factors of the upper and the lower parts are f_U = 1.5 and f_L = 0.75. The present breakwater physical model gives efficiency near other similar systems of different shapes.     

Computations of the energy flux to mixing processes via baroclinic wave drag on barotropic tides

The geographical distribution of barotropic to baroclinic transfer of tidal energy by baroclinic wave drag in the abyssal ocean is estimated. Using tidal velocities from a state-of-the-art numerical tidal model, the total loss of barotropic tidal energy in the deep ocean (between 70°S and 70°N and at depths greater than 1000 m) is estimated to be about 0.7 TW (M₂) corresponding to a mean value of the energy flux (e) of 2.4×10⁻³ W/m². The distribution of e is however highly skewed with a median of about 10⁻⁶ W/m². Only 10% of the area is responsible for more than 97% of the total energy transfer.To assess the possible influence of the relatively coarse bathymetry representation upon the present estimate, complementary calculations using better resolved sea floor topography are carried out over a control area around the Hawaiian Ridge. There are no major differences between the results achieved using the two different bathymetry databases. Fluxes of about 16 GW or 6×10⁻³ W/m² are computed in both cases, and the main contributions to the total fluxes originate in the same range of e-values and cover equally large parts of the total area.It is not clear whether the present model is valid at flat or subcritical bottom slopes. However, for the Hawaiian region, only 2% of the total energy flux as calculated in the present study originates in areas of critical and subcritical slopes.     

Large eddy simulations of a circular cylinder at Reynolds numbers surrounding the drag crisis

Large eddy simulations of the flow around a circular cylinder at high Reynolds numbers are reported. Five Reynolds numbers were chosen, such that the drag crisis was captured. A total of 18 cases were computed to investigate the effect of gridding strategy, turbulence modelling, numerical schemes and domain width on the results. It was found that unstructured grids provide better resolution of key flow features, when a ‘reasonable’ grid size is to be maintained.When using coarse grids for large eddy simulation, the effect of turbulence models and numerical schemes becomes more pronounced. The dynamic mixed Smagorinsky model was found to be superior to the Smagorinsky model, since the model coefficient is allowed to dynamically adjust based on the local flow and grid size. A blended upwind-central convection scheme was also found to provide the best accuracy, since a fully central scheme exhibits artificial wiggles, due to dispersion errors, which pollute the solution.Mean drag, fluctuating lift Strouhal number and base pressure are compared to experiments and empirical estimates for Reynolds numbers ranging from 6.31 × 10⁴ to 5.06 × 10⁵. In terms of the drag coefficient, the drag crisis is well captured by the present simulations, although the other integral quantities (rms lift and Strouhal number) show larger discrepancies. For the lowest Reynolds number, the drag is seen to be more sensitive to the domain width than the spanwise grid spacing, while at the higher Reynolds numbers the grid resolution plays a more important role, due to the larger extent of the turbulent boundary layer.