The Brinkman’s equation simplifies the numerical modelling of karst reservoirs by allowing the use of a single transport equation to model the flow of fluids in both the free flow and porous regions, in effect reducing the error arising from improper modelling of the interface between the two regions. However, most of the equations available to model flow within karst reservoirs deal with steady flow conditions. This approach however may not be accurate in reservoirs where unsteady conditions exist. We considered the effects of unsteady flow conditions in karst reservoirs by adding an unsteady flow term to the Brinkman’s equation. We solved the coupled conservation-transport equations that models unsteady fluid transport in karst reservoirs and then studied the effects of unsteady flow conditions on tracer transport in two different sample reservoirs. The solution method adopted is sequential and involves solving the unsteady Brinkman’s model first, followed by advection-diffusion-adsorption equation using the cell-centred finite volume approach. The same problems were also solved using a steady flow Brinkman’s model, and the results obtained were compared. were compared. The results show that, inside the caves, the unsteady Brinkman’s model yielded lower tracer concentrations at early times when compared to the steady flow model.


Article metrics loading...

Loading full text...

Full text loading...


  1. Borghi, A., Renard, P., & Cornaton, F.
    (2016). Can one identify karst conduit networks geometry and properties from hydraulic and tracer test data?Advances in Water Resources, 90, 99–115. http://doi.org/10.1016/j.advwatres.2016.02.009
    [Google Scholar]
  2. Cabras, S., De Waele, J., & Sanna, L.
    (2008). Caves and Karst Aquifer Drainage of Supramonte (Sardinia, Italy): A Review. Acta Carsologica, 37(2–3). http://doi.org/10.3986/ac.v37i2.148
    [Google Scholar]
  3. Göppert, N., & Goldscheider, N.
    (2007). Solute and Colloid Transport in Karst Conduits under Low-and High-Flow Conditions. Ground Water, http://doi.org/10.1111/j.1745-6584.2007.00373.x
    [Google Scholar]
  4. Kincaid, T. R., Hazlett, T. J., & Davies, G. J.
    (2002). Quantitative groundwater tracing and effective numerical modeling in karst: an example from the Woodville Karst Plain of North Florida. Ground Water, (850), 1–8. http://doi.org/10.1061/40796(177)13
    [Google Scholar]
  5. Oehlmann, S., Geyer, T., Licha, T., & Sauter, M.
    (2015). Reducing the ambiguity of karst aquifer models by pattern matching of flow and transport on catchment scale. Hydrology and Earth System Sciences, 19(2), 893–912. http://doi.org/10.5194/hess-19-893-2015
    [Google Scholar]
  6. Rivard, C., & Delay, F.
    (2004). Simulations of solute transport in fractured porous media using 2D percolation networks with uncorrelated hydraulic conductivity fields. Hydrogeology Journal, 12(6), 613–627. http://doi.org/10.1007/s10040-004-0363-z
    [Google Scholar]
  7. Smart, C.
    (1988). Artificial Tracer Techniques for the Determination of the Structure of Conduit Aquifers. Ground Water, 26(4), 445–453. http://doi.org/10.1111/j.1745-6584.1988.tb00411.x
    [Google Scholar]
  8. Staut, M., & Auersperger, P.
    (2006). Tracing of the Stream Flowing Through the Cave Ferranova Buža, Central Slovenia. Acta Carsologica, 35(2–3). http://doi.org/10.3986/ac.v35i2-3.231
    [Google Scholar]
  9. Weeks, S. W., & Sposito, G.
    (1998). Mixing and stretching efficiency in steady and unsteady groundwater flows. Water Resources Research, 34(12), 3315–3322. http://doi.org/10.1029/98WR02535
    [Google Scholar]
  10. Zhu, T., Waluga, C., Wohlmuth, B., & Manhart, M.
    (2014). A Study of the Time Constant in Unsteady Porous Media Flow Using Direct Numerical Simulation. Transport in Porous Media, 104(1), 161–179. http://doi.org/10.1007/s11242-014-0326-3
    [Google Scholar]

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error