AquiferLinkExchange.hpp

Go to the documentation of this file.

 
#ifndef OPENSWMM_ENGINE_HYDRAULICS_AQUIFER_LINK_EXCHANGE_HPP
#define OPENSWMM_ENGINE_HYDRAULICS_AQUIFER_LINK_EXCHANGE_HPP
 
#include <algorithm>
#include <cstddef>
 
#include "../core/SimulationContext.hpp"
 
namespace openswmm::hydraulics {

inline double aquiferSeepFactor(double depth, double gw_head_rel, double d_c) noexcept {
    if (!(d_c > 0.0))
        return 1.0; // no path length: legacy
    const double h_eff = (gw_head_rel > 0.0) ? gw_head_rel : 0.0;
    return std::min((depth - h_eff) / d_c, 1.0);
}

inline double capAquiferGain(double seep_rate, double gain_max) noexcept {
    return (seep_rate < -gain_max) ? -gain_max : seep_rate;
}

inline void applySeepageForcing(const SimulationContext& ctx, std::size_t link,
                                double& seep_rate) noexcept {
    const auto& f = ctx.forcing;
    if (link >= f.link_seepage_mode.size())
        return;
    switch (f.link_seepage_mode[link]) {
    case ForcingMode::OVERRIDE:
        seep_rate = f.link_seepage_value[link];
        break;
    case ForcingMode::ADD:
        seep_rate += f.link_seepage_value[link];
        break;
    default:
        break;
    }
}
 
} // namespace openswmm::hydraulics
 
#endif // OPENSWMM_ENGINE_HYDRAULICS_AQUIFER_LINK_EXCHANGE_HPP