# Search Results

## You are looking at 1 - 10 of 23 items for

- Author or Editor: Hans Burchard x

- All content x

## Abstract

This short paper focuses on the problems caused by the equal contribution of shear and buoyancy production in the length scale equation transporting the product of turbulent kinetic energy *q*
^{2} and macro length scale *l* as suggested by Mellor and Yamada. It is shown that this model has no steady-state solutions for homogeneous shear flows. The concept of the steady-state Richardson number is used for estimating the adequate contribution of the buoyancy production to the equation for *q*
^{2}
*l.* In a simple wind entrainment experiment, the failure of the *q*
^{2}
*l* equation in the case of unlimited *l* is shown. With the new estimate for the buoyancy production contribution to *q*
^{2}
*l,* physically sound results are obtained even for the unlimited case. In applications of the modified *q*
^{2}
*l* equation to three oceanic test cases in the northern North Sea and the northern Pacific, the results of the wind entrainment experiment are confirmed.

## Abstract

This short paper focuses on the problems caused by the equal contribution of shear and buoyancy production in the length scale equation transporting the product of turbulent kinetic energy *q*
^{2} and macro length scale *l* as suggested by Mellor and Yamada. It is shown that this model has no steady-state solutions for homogeneous shear flows. The concept of the steady-state Richardson number is used for estimating the adequate contribution of the buoyancy production to the equation for *q*
^{2}
*l.* In a simple wind entrainment experiment, the failure of the *q*
^{2}
*l* equation in the case of unlimited *l* is shown. With the new estimate for the buoyancy production contribution to *q*
^{2}
*l,* physically sound results are obtained even for the unlimited case. In applications of the modified *q*
^{2}
*l* equation to three oceanic test cases in the northern North Sea and the northern Pacific, the results of the wind entrainment experiment are confirmed.

## Abstract

A universal law of estuarine mixing is derived here, combining the approaches of salinity coordinates, Knudsen relations, total exchange flow, mixing definition as salinity variance loss, and the mixing–exchange flow relation. As a result, the long-term average mixing within an estuarine volume bounded by the isohaline of salinity *S* amounts to *M*(*S*) = *S*
^{2}
*Q*
_{r}, where *Q*
_{r} is the average river runoff into the estuary. Consequently, the mixing per salinity class is *m*(*S*) = ∂_{S}
*M*(*S*) = 2*SQ*
_{r}, which can also be expressed as the product of the isohaline volume and the mixing averaged over the isohaline. The major differences between the new mixing law and the recently developed mixing relation based on the Knudsen relations are threefold: (i) it does not depend on internal dynamics of the estuary determining inflow and outflow salinities (universality), (ii) it is exactly derived from conservation laws (accuracy), and (iii) it calculates mixing per salinity class (locality). The universal mixing law is demonstrated by means of analytical stationary and one-dimensional and two-dimensional numerical test cases. Some possible consequences for the salinity distribution in real estuaries are briefly discussed. Since the mixing per salinity class only depends on the river runoff and the chosen salinity, and not on local processes at the isohaline, low-mixing estuaries must have large isohaline volumes and vice versa.

## Abstract

A universal law of estuarine mixing is derived here, combining the approaches of salinity coordinates, Knudsen relations, total exchange flow, mixing definition as salinity variance loss, and the mixing–exchange flow relation. As a result, the long-term average mixing within an estuarine volume bounded by the isohaline of salinity *S* amounts to *M*(*S*) = *S*
^{2}
*Q*
_{r}, where *Q*
_{r} is the average river runoff into the estuary. Consequently, the mixing per salinity class is *m*(*S*) = ∂_{S}
*M*(*S*) = 2*SQ*
_{r}, which can also be expressed as the product of the isohaline volume and the mixing averaged over the isohaline. The major differences between the new mixing law and the recently developed mixing relation based on the Knudsen relations are threefold: (i) it does not depend on internal dynamics of the estuary determining inflow and outflow salinities (universality), (ii) it is exactly derived from conservation laws (accuracy), and (iii) it calculates mixing per salinity class (locality). The universal mixing law is demonstrated by means of analytical stationary and one-dimensional and two-dimensional numerical test cases. Some possible consequences for the salinity distribution in real estuaries are briefly discussed. Since the mixing per salinity class only depends on the river runoff and the chosen salinity, and not on local processes at the isohaline, low-mixing estuaries must have large isohaline volumes and vice versa.

## Abstract

The aim of this one-dimensional water column study is to combine modified versions of three characteristic parameters for periodic tidal flow under the influence of a longitudinal buoyancy gradient—the horizontal Richardson number, the inverse Strouhal number, and the inverse Ekman number—into a parameter space study, including constant wind forcing from various directions. It is shown how the underlying dynamical equations can be cast into nondimensional form, depending mainly on these three nondimensional parameters plus the relative wind speed and the wind direction. Idealized model simulations are carried out for the whole realistic range of horizontal Richardson and inverse Strouhal numbers for various latitudes, showing the amplitude of the tidally induced stratification for a wide range of scenarios. It is found that classical threshold values for the horizontal Richardson number, indicating the switch from periodic stratification to permanent stratification, are valid only for special cases, and that this switch also strongly depends on the inverse Strouhal number, the inverse Ekman number, and the wind vector. The transverse residual flow is close to the thermal wind balance for a variety of parameters, and nondimensionalized longitudinal residual flow shows the classical estuarine exchange flow pattern with little variation in the near-bed onshore component. Wind straining is confirmed as an important estuarine and coastal process, enhancing estuarine circulation for offshore (down estuary) winds and vice versa. Agreement with field data from Liverpool Bay is good, including the explanation of a flood tide local maximum of the dissipation rate in the upper half of the water column. An equation for the second time derivative of the potential energy anomaly is derived for quantifying the dynamical processes leading to stratification due to the straining of the horizontal density gradient.

## Abstract

The aim of this one-dimensional water column study is to combine modified versions of three characteristic parameters for periodic tidal flow under the influence of a longitudinal buoyancy gradient—the horizontal Richardson number, the inverse Strouhal number, and the inverse Ekman number—into a parameter space study, including constant wind forcing from various directions. It is shown how the underlying dynamical equations can be cast into nondimensional form, depending mainly on these three nondimensional parameters plus the relative wind speed and the wind direction. Idealized model simulations are carried out for the whole realistic range of horizontal Richardson and inverse Strouhal numbers for various latitudes, showing the amplitude of the tidally induced stratification for a wide range of scenarios. It is found that classical threshold values for the horizontal Richardson number, indicating the switch from periodic stratification to permanent stratification, are valid only for special cases, and that this switch also strongly depends on the inverse Strouhal number, the inverse Ekman number, and the wind vector. The transverse residual flow is close to the thermal wind balance for a variety of parameters, and nondimensionalized longitudinal residual flow shows the classical estuarine exchange flow pattern with little variation in the near-bed onshore component. Wind straining is confirmed as an important estuarine and coastal process, enhancing estuarine circulation for offshore (down estuary) winds and vice versa. Agreement with field data from Liverpool Bay is good, including the explanation of a flood tide local maximum of the dissipation rate in the upper half of the water column. An equation for the second time derivative of the potential energy anomaly is derived for quantifying the dynamical processes leading to stratification due to the straining of the horizontal density gradient.

## Abstract

The purpose of this paper is to modify two-equation turbulence models such that they are capable of simulating dynamics in the wave-enhanced layer near the surface. A balance of diffusion of turbulent kinetic energy (TKE) and dissipation is assumed as the surface boundary condition for TKE following the suggestion of Craig and Banner. It is shown that this theory, originally developed under the assumption of a macro length scale linearly increasing down from the surface, fails for two-equation models such as the well-known *k*–ε model. Suggestions are made how to modify such models for overcoming this deficiency. The basic idea is to insert the analytic solution of a model problem suggested by Craig into the dissipation rate equation and solve for the turbulent Schmidt number of the dissipation rate equation, which may be formulated as a function of the production/dissipation ratio. With this modification, the linear behavior of the macro length scale is properly reproduced by the *k*–ε model. It is shown how near-surface dissipation rate measurements under breaking waves can be simulated by an extended *k*–ε model considering a shear-dependent closure for the second moments. Finally, the overall performance of this new model approach is tested with a typical upper mixed layer scenario in the northern North Sea.

## Abstract

The purpose of this paper is to modify two-equation turbulence models such that they are capable of simulating dynamics in the wave-enhanced layer near the surface. A balance of diffusion of turbulent kinetic energy (TKE) and dissipation is assumed as the surface boundary condition for TKE following the suggestion of Craig and Banner. It is shown that this theory, originally developed under the assumption of a macro length scale linearly increasing down from the surface, fails for two-equation models such as the well-known *k*–ε model. Suggestions are made how to modify such models for overcoming this deficiency. The basic idea is to insert the analytic solution of a model problem suggested by Craig into the dissipation rate equation and solve for the turbulent Schmidt number of the dissipation rate equation, which may be formulated as a function of the production/dissipation ratio. With this modification, the linear behavior of the macro length scale is properly reproduced by the *k*–ε model. It is shown how near-surface dissipation rate measurements under breaking waves can be simulated by an extended *k*–ε model considering a shear-dependent closure for the second moments. Finally, the overall performance of this new model approach is tested with a typical upper mixed layer scenario in the northern North Sea.

## Abstract

The interaction of shear, stratification, and turbulence in boundary layers on sloping topography is investigated with the help of an idealized theoretical model, assuming uniform bottom slope and homogeneity in the upslope direction. It is shown theoretically that the irreversible vertical buoyancy flux generated in the boundary layer is directly proportional to the molecular destruction rate of small-scale buoyancy variance, which can be inferred, for example, from microstructure observations. Dimensional analysis of the equations shows that, for harmonic boundary layer forcing and no rotation, the problem is governed by three nondimensional parameters (slope angle, roughness number, and ratio of forcing and buoyancy frequencies). Solution of the equations with a second-moment closure model for the turbulent fluxes reveals the periodic generation of gravitationally unstable boundary layers during upslope flow, consistent with available observations. Investigation of the nondimensional parameter space with the help of this model illustrates a systematic increase of the bulk mixing efficiencies for (i) steep slopes and (ii) low-frequency forcing. Except for very steep slopes, mixing efficiencies are substantially smaller than the classical value of Γ = 0.2.

## Abstract

The interaction of shear, stratification, and turbulence in boundary layers on sloping topography is investigated with the help of an idealized theoretical model, assuming uniform bottom slope and homogeneity in the upslope direction. It is shown theoretically that the irreversible vertical buoyancy flux generated in the boundary layer is directly proportional to the molecular destruction rate of small-scale buoyancy variance, which can be inferred, for example, from microstructure observations. Dimensional analysis of the equations shows that, for harmonic boundary layer forcing and no rotation, the problem is governed by three nondimensional parameters (slope angle, roughness number, and ratio of forcing and buoyancy frequencies). Solution of the equations with a second-moment closure model for the turbulent fluxes reveals the periodic generation of gravitationally unstable boundary layers during upslope flow, consistent with available observations. Investigation of the nondimensional parameter space with the help of this model illustrates a systematic increase of the bulk mixing efficiencies for (i) steep slopes and (ii) low-frequency forcing. Except for very steep slopes, mixing efficiencies are substantially smaller than the classical value of Γ = 0.2.

## Abstract

In this comparative study, four different algebraic second-moment turbulence closure models are investigated in detail. These closure schemes differ in the number of terms considered for the closure of the pressure–strain correlations. These four turbulence closures result in the eddy-diffusivity principle such that the closure assumptions are contained in dimensionless so-called stability functions. Their performance in terms of Prandtl number, Monin–Obukhov similarity theory, and length scale ratios are first tested against data for simple flows. The turbulence closure is then completed by means of a *k*–*ϵ* two-equation model, but other models such as the two-equation model by Mellor and Yamada could also be used. The concept of the steady-state Richardson number for homogeneous shear layers is exploited for calibrating the sensitivity of the four models to shear and stable stratification. Idealized simulations of mixed layer entrainment into stably stratified flow due to surface stress and due to free convection are carried out. For the latter experiment, comparison to recent large eddy simulation data is made. Finally, the well-known temperature profile data at OWS Papa are simulated for an annual cycle. The main result of this paper is that the overall performance of the new second-moment closure model by Canuto et al.—expressed as nondimensional stability functions—is superior compared to the others in terms of physical soundness, predictability, computational economy, and numerical robustness.

## Abstract

In this comparative study, four different algebraic second-moment turbulence closure models are investigated in detail. These closure schemes differ in the number of terms considered for the closure of the pressure–strain correlations. These four turbulence closures result in the eddy-diffusivity principle such that the closure assumptions are contained in dimensionless so-called stability functions. Their performance in terms of Prandtl number, Monin–Obukhov similarity theory, and length scale ratios are first tested against data for simple flows. The turbulence closure is then completed by means of a *k*–*ϵ* two-equation model, but other models such as the two-equation model by Mellor and Yamada could also be used. The concept of the steady-state Richardson number for homogeneous shear layers is exploited for calibrating the sensitivity of the four models to shear and stable stratification. Idealized simulations of mixed layer entrainment into stably stratified flow due to surface stress and due to free convection are carried out. For the latter experiment, comparison to recent large eddy simulation data is made. Finally, the well-known temperature profile data at OWS Papa are simulated for an annual cycle. The main result of this paper is that the overall performance of the new second-moment closure model by Canuto et al.—expressed as nondimensional stability functions—is superior compared to the others in terms of physical soundness, predictability, computational economy, and numerical robustness.

## Abstract

By means of a numerical model of an idealized flat-bottom estuary, the paper studies the hydrodynamic control of the turbidity zone by the combined effect of the salt wedge and tidal movements. The model is of two- dimensional (*x, z*) finite-difference type with high resolution in time and space. It computes momentum, surface elevation, salinity, suspended particulate matter (SPM), turbulent kinetic energy, and dissipation rate as prognostic state variables. At the seaward boundary a tidal forcing is applied. At the landward boundary a weir is situated where a constant freshwater discharge is prescribed. The initial SPM concentration is horizontally homogeneous. After simulating a few tidal periods the model results exhibit the evolution of a stable SPM peak (the estuarine turbidity maximum or ETM) at the tip of the salt wedge. An inspection of the tidal mean velocity profiles around the ETM shows that this trapping of SPM is due to a residual near-bottom upstream current in the region of the salt wedge. Three physical causes for this residual countercurrent are investigated in greater detail by numerical experiments, namely, (i) the *residual gravitational circulation,* (ii) the *tidal velocity asymmetry,* and (iii) *the tidal mixing asymmetry.* The first mechanism is related to the baroclinic part of the longitudinal pressure gradient. The second and third mechanism are based on the differences between the vertical profiles of velocity and SPM, respectively, at flood and ebb tide. For the macrotidal estuary considered here, the consideration of both (i) and (ii) could be shown to be *necessary* for the establishment of an ETM in the considered idealized estuary. It could further be shown that (iii) affects the ETM formation only quantitatively but not qualitatively and appears to be not necessary for the existence of an ETM.

## Abstract

By means of a numerical model of an idealized flat-bottom estuary, the paper studies the hydrodynamic control of the turbidity zone by the combined effect of the salt wedge and tidal movements. The model is of two- dimensional (*x, z*) finite-difference type with high resolution in time and space. It computes momentum, surface elevation, salinity, suspended particulate matter (SPM), turbulent kinetic energy, and dissipation rate as prognostic state variables. At the seaward boundary a tidal forcing is applied. At the landward boundary a weir is situated where a constant freshwater discharge is prescribed. The initial SPM concentration is horizontally homogeneous. After simulating a few tidal periods the model results exhibit the evolution of a stable SPM peak (the estuarine turbidity maximum or ETM) at the tip of the salt wedge. An inspection of the tidal mean velocity profiles around the ETM shows that this trapping of SPM is due to a residual near-bottom upstream current in the region of the salt wedge. Three physical causes for this residual countercurrent are investigated in greater detail by numerical experiments, namely, (i) the *residual gravitational circulation,* (ii) the *tidal velocity asymmetry,* and (iii) *the tidal mixing asymmetry.* The first mechanism is related to the baroclinic part of the longitudinal pressure gradient. The second and third mechanism are based on the differences between the vertical profiles of velocity and SPM, respectively, at flood and ebb tide. For the macrotidal estuary considered here, the consideration of both (i) and (ii) could be shown to be *necessary* for the establishment of an ETM in the considered idealized estuary. It could further be shown that (iii) affects the ETM formation only quantitatively but not qualitatively and appears to be not necessary for the existence of an ETM.

## Abstract

In straight tidal estuaries, residual overturning circulation results mainly from a competition between gravitational forcing, wind forcing, and friction. To systematically investigate this for tidally energetic estuaries, the dynamics of estuarine cross sections is analyzed in terms of the relation between gravitational forcing, wind stress, and the strength of estuarine circulation. A system-dependent basic Wedderburn number

## Abstract

In straight tidal estuaries, residual overturning circulation results mainly from a competition between gravitational forcing, wind forcing, and friction. To systematically investigate this for tidally energetic estuaries, the dynamics of estuarine cross sections is analyzed in terms of the relation between gravitational forcing, wind stress, and the strength of estuarine circulation. A system-dependent basic Wedderburn number

## Abstract

Recent finescale observations of shear and stratification in temperate shelf sea thermoclines show that they are of marginal stability, suggesting that episodes of enhanced shear could potentially lead to shear instability and diapcynal mixing. The bulk shear between the upper and lower boundary layers in seasonally stratified shelf seas shows remarkable variability on tidal, inertial, and synoptic time scales that has yet to be explained. In this paper observations from the seasonally stratified northern North Sea are presented for a time when the water column has a distinct two-layer structure. Bulk shear estimates, based on ADCP measurements, show a bulk shear vector that rotates in a clockwise direction at the local inertial period, with episodes of bulk shear spikes that have an approximately twice daily period, and occur in bursts that last for several days. To explain this observation, a simple two-layer model based on layer averaging of the one-dimensional momentum equation is developed, forced at the surface by wind stress and damped by (tidally dominated) sea bed friction. The two layers are then linked through an interfacial stress term. The model reproduces the observations, showing that the bulk shear spikes are a result of the alignment of the wind stress, tidal bed stress, and (clockwise rotating) bulk shear vectors. Velocity microstructure measurements are then used to confirm enhanced levels of mixing during a period of bulk shear spikes. A numerical study demonstrates the sensitivity of the spike generation mechanism to the local tidal conditions and the phasing and duration of wind events.

## Abstract

Recent finescale observations of shear and stratification in temperate shelf sea thermoclines show that they are of marginal stability, suggesting that episodes of enhanced shear could potentially lead to shear instability and diapcynal mixing. The bulk shear between the upper and lower boundary layers in seasonally stratified shelf seas shows remarkable variability on tidal, inertial, and synoptic time scales that has yet to be explained. In this paper observations from the seasonally stratified northern North Sea are presented for a time when the water column has a distinct two-layer structure. Bulk shear estimates, based on ADCP measurements, show a bulk shear vector that rotates in a clockwise direction at the local inertial period, with episodes of bulk shear spikes that have an approximately twice daily period, and occur in bursts that last for several days. To explain this observation, a simple two-layer model based on layer averaging of the one-dimensional momentum equation is developed, forced at the surface by wind stress and damped by (tidally dominated) sea bed friction. The two layers are then linked through an interfacial stress term. The model reproduces the observations, showing that the bulk shear spikes are a result of the alignment of the wind stress, tidal bed stress, and (clockwise rotating) bulk shear vectors. Velocity microstructure measurements are then used to confirm enhanced levels of mixing during a period of bulk shear spikes. A numerical study demonstrates the sensitivity of the spike generation mechanism to the local tidal conditions and the phasing and duration of wind events.

## Abstract

Tidal straining, which can mathematically be described as the covariance between eddy viscosity and vertical shear of the along-channel velocity component, has been acknowledged as one of the major drivers for estuarine circulation in channelized tidally energetic estuaries. In this paper, the authors investigate the role of lateral circulation for generating this covariance. Five numerical experiments are carried out, starting with a reference scenario including the full physics and four scenarios in which specific key physical processes are neglected. These processes are longitudinal internal pressure gradient forcing, lateral internal pressure gradient forcing, lateral advection, and the neglect of temporal variation of eddy viscosity. The results for the viscosity–shear covariance are correlated across different experiments to quantify the change due to neglect of these key processes. It is found that the lateral advection of vertical shear of the along-channel velocity component and its interaction with the tidally asymmetric eddy viscosity (which is also modified by the lateral circulation) is the major driving force for estuarine circulation in well-mixed tidal estuaries.

## Abstract

Tidal straining, which can mathematically be described as the covariance between eddy viscosity and vertical shear of the along-channel velocity component, has been acknowledged as one of the major drivers for estuarine circulation in channelized tidally energetic estuaries. In this paper, the authors investigate the role of lateral circulation for generating this covariance. Five numerical experiments are carried out, starting with a reference scenario including the full physics and four scenarios in which specific key physical processes are neglected. These processes are longitudinal internal pressure gradient forcing, lateral internal pressure gradient forcing, lateral advection, and the neglect of temporal variation of eddy viscosity. The results for the viscosity–shear covariance are correlated across different experiments to quantify the change due to neglect of these key processes. It is found that the lateral advection of vertical shear of the along-channel velocity component and its interaction with the tidally asymmetric eddy viscosity (which is also modified by the lateral circulation) is the major driving force for estuarine circulation in well-mixed tidal estuaries.