Additional options
Disable SZ solute flux to dummy UZ
The following Extra Parameter is useful, if you are using an alternative UZ model, such as DAISY, in MIKE SHE and you are trying to couple it to the WQ.
In this case, you will be typically using the Negative Precipitation (V1 p. 737) option. If you use this option, then you will not use a MIKE SHE UZ, and the UZ-SZ exchange will pass through a “dummy UZ” layer. When this is coupled to the water quality, solutes will also be passed to this dummy UZ layer and removed from the SZ domain and the model.
To prevent the upflow of solutes from SZ to the dummy UZ, you must specify the following Extra Parameter..
|
Parameter Name |
Type |
Value |
|---|---|---|
|
disable sz transport to dummy uz |
Boolean |
On |
A detailed test of the MIKE SHE WQ engine comparing an SZ model with fixed concentration at an inflow boundary with an analytical solution for a fixed concentration source, showed that MIKE SHE under-estimates the mass flux into the model when the model includes longitudinal dispersion.
The problem is that the SZ transport scheeme (QUICKEST) doesn't include dispersive transport to/from open boundary cells. This is as designed, but apparently not correct. After including the boundary dispersion, the mass input to the model is within 2 % of the analytical solution.
From Release 2011 and onwards, the boundary dispersion has been made optional for backwards compatibility and is activated with the extra-parameter: .
|
Parameter Name |
Type |
Value |
|---|---|---|
|
enable sz boundary dispersion |
Boolean |
On |
However, the SZ boundary dispersion option (above) does not calculate dispersive transport to an inflow boundary correctly. Again, this problem was identified in the tests of MShe_WQ with MIKE ECO Lab vs analytical solution. For example:
· Species 1 enters the model via an inflow (flux) boundary with fixed concentration - including dispersive transport due to the new sz boundary dispersion option.
· Species 1 decays to Species 2 which again decays to Species 3.
· The concentrations of Sp2 & Sp3 are too high, especially close to the inflow boundary.
The analytical solution includes dispersive transport of Sp2 & Sp3 against the flow direction because the concentration of these species are 0 at the boundary. However, this dispersive mass flux to the boundary is not included in the SZ solver due to an old check in the code. When mass flux to/from a boundary point is reversed compared to the flow direction, the mass flux is simply reset to 0.
This made sense before the boundary dispersion was implemented because advective transport against the flow direction would be wrong. But, now, when the boundary dispersion is active, this situation is allowed.