Page Comparison

You can change between the Piecewise-linear Integral (default) and the Backward Euler reservoir routing methods in (default) and the Piecewise-linear Integral in Scenario Options » Storages.

...

And where ε is the mass-balance error. 

Release and spill are the sum of the water flowing downstream in all of the outflow paths. For each path we have two relationships:

...

Table 1. Example of pre-computed solution table

Image Modified
 

These values stay the same for the length of the model run and can then be used at each time step. Using the current value of Orders_j and E_depth we can calculate the corresponding value of ε for each table row (V):

...

At which point storage routing becomes a table lookup and interpolation exercise where we find the value of V that give a ϵ value of zero. Obviously we don't have to calculate the entire table but just use a standard search technique.

...

Figure 2. Multiple outlet priorities, Process overview

...

1. Check to see if the resulting reservoir volume is outside the range allowed (otherwise the solution from step 1 is acceptable and we can stop).
   
   
2. If it is outside of the range allowed then see if it is possible to get back into the acceptable range by either:
   1. If initial solution is below the minimum target then recompute the time step assuming that there are no downstream orders
   2. If initial solution is above the maximum target then recompute the time step assuming that the downstream orders are the greater of the current order and the maximum flow rate allowed on each outlet path for releasing above target water (the safe release).
      
      
3. If the resulting answer from step 3 is still outside the allowed range the result obtained at step 3 becomes the final answer for this time step.
   
   

4. Compute a total release from:
   
    

   
   Where S_target is the appropriate minimum or maximum target volume (are you trying to cut back releases to stop the reservoir dropping below the minimum or increasing releases to stop it going over the maximum).
   
   

5. To calculate the corresponding release rates for the solution obtained at step 5 we have to scale up or down the orders so when we evaluate equation (5) we'll get the required results. The two cases are:
   1. If the initial solution (step 1) fell below the minimum target volume then we are cutting back and S_target at step 5 was S_min. We cut back each path in proportion to how far their order is above fjminSmin.

      1. Calculate the sum of current demands in excess of the path minimum capacity (limiting to the maximum outlet capacity at the minimum target volume):
         
          

      2. Calculate the release cutback required:
         
         

      3. Scale back the orders: 
         

   2. If the initial solution (step 1) fell above the maximum target volume then we are ramping up the releases and S_target at step 5 was S_max. We add onto each path in proportion to how much remaining release capacity each path has.

      1. Calculate maximum we can currently send down each path:
         
         

         
         Where AllowedQ_j is the maximum permitted release down path j for above target water.
         
         

      2. the flow that's already going to go down path j: 
         
         

      3. the remaining release capacity: 
         
         

      4. The required increase in release is then: 
         
         

      5. Modified orders are then calculated from: