2. Materials and Methods
In two-dimensional river mixing, mixing data can be classified into two types; time-variant data and steady data. The former can be acquired from a transient concentration situation by injecting a tracer instantaneously. The latter is acquired under steady concentration conditions from a continuous injection of a tracer. Most field studies based on continuous tracer injection had been conducted during 1970–80s. When the channel discharge and the solute inflow rate are both steady, the downstream plume also eventually becomes steady; time derivative terms vanish. Furthermore, under this condition the longitudinal mixing term may be dropped in the two-dimensional mass transport model. It is well known that dispersive transport is small relative to convective transport in a unidirectional flow [
4]. Moment-based methods were derived in such conditions and were commonly used to calculate the transverse dispersion coefficient by researchers [
4,
5,
6,
7]. Of course, the moment-based methods also may be applied to the conditions of time-variant concentration using Beltaos’ conversion procedure [
8]. In such cases, time variations within the concentration data were vanished by using a defined tracer dosage.
where
is the dosage;
C is the depth-averaged concentration of a tracer;
t is time; and
x and
y are the longitudinal transverse coordinates respectively.
The generalized moment method (generalized MM) by Holley et al. [
6], which can reflect the effects of the transverse velocities and tracer impinging on the banks, for the
-equation yields.
where
is the transverse dispersion coefficients;
is the channel width;
is the second moment of the transverse distribution of the concentration data;
and
are the terms which reflect the tracer impinging on the banks and the transverse velocities, respectively;
and
are the depth-averaged longitudinal and transverse velocities, respectively. We can calculate the transverse dispersion coefficient using the slope of the straight line which is fitted to the plot of variance (
) against longitudinal distance (
). The simple moment method (simple MM) by Sayre and Chang [
4] can be derived from the generalized MM, Equation (2). Neglecting the transverse velocity and the tracer impinging on banks, and assuming a constant longitudinal velocity, Equation (2) becomes:
where
is the reach-averaged longitudinal velocity. The stream–tube moment method (stream–tube MM) by Beltaos [
7] can reflect irregularities of depth and width at rivers by using the stream–tube concept. The method is:
where
and
are the second and first moments of the
distribution, respectively;
;
q and
Q are the cumulative flow and the total flow discharge, respectively;
;
;
and
are the normalized dosages at the left and right banks, respectively;
is the mean depth; and
is a normalized shape–velocity factor that has a range of 1.0–3.6 [
7] and is defined by:
These moment-based methods, however, have some restrictions in which the skewed concentration profile induced by river irregularities makes it difficult to compute a meaningful value for the second moments. This may lead to an inaccurate dispersion coefficient. Seo et al. [
9] developed the stream–tube routing procedure (stream–tube RP), a routing procedure combined with the stream–tube concept to overcome the weak point of the moment-based methods. The equation is:
where
is the observed dosage distribution at an upstream
;
is the predicted dosage distribution at a downstream
;
; and
is a transverse distance variable for the integration. The merit of the routing procedure is that it is relatively simple and easy in comparison to constructing numerical simulation models. The routing procedure is, of course, a type of inverse method where the dispersion coefficient is determined by trial and error until the model results agree with tracer tests or other observations.
When it is needed to acquire both the transverse dispersion and longitudinal dispersion coefficients from tracer data under time-variant concentration situation, other routing procedures should be employed to the observed dispersion coefficients. Baek et al. [
10] expanded the one-dimensional equation by Fischer [
11] into two-dimensional equation. The equation (2D RP) is:
where
is the predicted concentration at a downstream section,
;
is the observed concentration at a upstream section,
;
and
are the mean times of passage in sections
and
, respectively;
is a time variable for the integration;
is a distance variable for the integration;
is the longitudinal dispersion coefficient. The calculated concentration by Equation (7) is matched repeatedly with the measured concentration by changing the dispersion coefficient until the differences between two concentrations become minimized. The value for the dispersion coefficient with the least error in concentration is regarded as the observed dispersion coefficient [
12].
2D RP (Equation (7)) also has limitations. It cannot reflect the irregularities of a river—such as non-uniformities of the bed and the channel width, a skewed velocity profile, meandering—and so on. To consider such irregularities, Baek and Seo [
13] derived another routing procedure from Equation (7) combining with the stream–tube concept. The equation (2D stream-tube RP) is:
where
is the predicted concentration at
; and
is the measured concentration at
. The above methods, including the moment-based method and the routing procedure, to observe the dispersion coefficients are summarized in
Table 1.
3. Results
The strengths and limitations of the methods described in the previous section are obvious, as the methods were derived under specific original assumptions. Application examples may be found in previous works, such as [
2,
3,
10,
13]. Baek et al. [
10] showed comparisons of the observed dispersion coefficients calculated by both the moment-based methods such as the simple MM and generalized MM and the routing procedure, 2D RP, based on tracer mixing data sets which were acquired from the artificial meandering channel with rectangular cross section. The results are depicted in
Figure 1. We can see that observed values of the dispersion coefficient from the moment-based methods and the routing procedure are in the same range. More than half values resulting from the generalized MM are higher than others because this method accounts for the effect of the transverse velocity and the tracer impinging on the banks. We also deduced from this figure that the time-variant concentration condition included in the routing procedure does not affect significantly the mechanism of transverse dispersion. This is because unsteadiness was incorporated into the longitudinal dispersion term, not the transverse dispersion term as shown in Equation (7).
The observed dispersion coefficients in this laboratory channel were further validated by Seo et al. [
2], in which they calculated the dispersion coefficients at each section based on a transport numerical modeling. According to Seo et al. [
2], a finite element scheme (specifically the Petrov–Galerkin type) was used to construct a numerical model, and the simulation results from the model were verified and compared with the tracer experimental data in the laboratory curved channel. The comparison is depicted in
Figure 2. Seo et al. [
2] showed that the dispersion coefficient acquired directly from velocity profiles (DC1 as shown
Figure 2b) supplied a more precise solution than the coefficient based on the dispersion tensor (DC2 as shown
Figure 2c). Additionally, Seo et al. [
2] inversely calculated the dispersion coefficient at each section based on the numerical model to reveal the spatial variations along the curved channel. The results are shown in
Figure 3a. Baek and Seo [
3] calculated inversely the dispersion coefficient along the identical channel by using 2D RP, as illustrated in
Figure 3b. As shown in this figure, the inversely calculated values of dispersion coefficient by the numerical model and 2D RP are in good agreement with each other along the channel.
Baek and Seo [
13] presented comparisons of the observed dispersion coefficients by 2D stream–tube RP and the stream–tube RP based on tracer data sets which were acquired from experiments conducted at bends of natural streams. In this study, several experimental results from tracer field tests were also added to verify and compare the two methods. Nine cases of tracer experiments were conducted at five streams in Korea. Rhodamine WT was selected as a tracer—a fluorescent dye that has been widely used for mixing study in streams—as it is conservative material and easily detectable with low background concentration [
14]. Using combined dispersion data sets (Baek and Seo [
13], and this study), observed dispersion coefficients by two routing procedures are compared in
Figure 4. Theoretically, the values for the transverse dispersion coefficient by the two routing procedures should be identical [
13]. From
Figure 4, it can be seen that there is generally good agreement between values calculated by 2D stream–tube RP and those calculated by stream–tube RP, although there are somewhat scattered trends.
Recently, Baek and Seo [
3] showed comparison of the observed values by the simple MM, the stream–tube MM, 2D RP, and 2D stream–tube MM based on nine experimental cases which were conducted at natural rivers, including Baek and Seo [
13]. They concluded that both the stream–tube RP and 2D stream–tube RP generated reasonable dispersion coefficients values that reflected the irregularities of the river’s geometry and hydraulics.