SoilCharacteristic.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: Apache-2.0
//
// Copyright 2026 Caleb Buahin
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//     http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.

 
#ifndef OPENSWMM_ENGINE_2D_SOIL_CHARACTERISTIC_HPP
#define OPENSWMM_ENGINE_2D_SOIL_CHARACTERISTIC_HPP
 
#include "SubsurfaceData.hpp"
 
namespace openswmm::twoD::soil {

struct Params {
    SoilChar law     = SoilChar::RUSSO;
    double   Ks      = 1.0e-5;
    double   theta_s = 0.45;
    double   theta_r = 0.10;
    double   alpha   = 2.0;
    double   psi_b   = 0.20;
    double   lambda  = 0.40;
    double   vg_n    = 1.6;
    double   vg_L    = 0.5;
};

double effectiveSaturation(const Params& p, double psi) noexcept;

double waterContent(const Params& p, double psi) noexcept;

double relativeConductivity(const Params& p, double psi) noexcept;

double conductivity(const Params& p, double psi) noexcept;

double suctionAtSaturation(const Params& p, double Se) noexcept;

double diffusivity(const Params& p, double psi) noexcept;

double equilibriumStorage(const Params& p, double L) noexcept;

double rechargeQ0(const Params& p, double L, double hu) noexcept;

double enslavedStorage(const Params& p, double L) noexcept;

double alphaL(const Params& p, double L) noexcept;

double feddesStress(double psi, double psi_w) noexcept;
 
}  // namespace openswmm::twoD::soil
 
#endif  // OPENSWMM_ENGINE_2D_SOIL_CHARACTERISTIC_HPP