SurfaceFluxCalculator.hpp

Go to the documentation of this file.

 
#ifndef OPENSWMM_ENGINE_2D_SURFACE_FLUX_CALCULATOR_HPP
#define OPENSWMM_ENGINE_2D_SURFACE_FLUX_CALCULATOR_HPP
 
#include "../data/MeshData.hpp"
#include "../data/SurfaceStateData.hpp"
#include "../data/SolverOptions2D.hpp"
 
namespace openswmm::twoD {

void computeUnlimitedGradients(const MeshData& mesh, SurfaceStateData& state,
                                int nthreads = 1);

void computeLimitedGradients(const MeshData& mesh, SurfaceStateData& state,
                              double epsilon, int nthreads = 1);

void computeEdgeFluxes(const MeshData& mesh, SurfaceStateData& state,
                        const SolverOptions2D& opts);

inline double evapSink(double rate, double depth, double dry_depth) noexcept {
    if (rate <= 0.0 || depth <= 0.0) return 0.0;
    if (depth >= dry_depth) return rate;
    double t = depth / dry_depth;
    return rate * t * t * (3.0 - 2.0 * t);
}

void assembleRHS(const MeshData& mesh, const SurfaceStateData& state,
                  const SolverOptions2D& opts, double* ydot);

void assembleImplicitRHS(const MeshData& mesh, const SurfaceStateData& state,
                          const SolverOptions2D& opts, double* ydot);

void assembleExplicitRHS(const MeshData& mesh, const SurfaceStateData& state,
                          const SolverOptions2D& opts, const double* y,
                          double* ydot);

void computeCellContinuity(const MeshData& mesh, SurfaceStateData& state,
                            const SolverOptions2D& opts, double dt);

void computeFaceVelocity(const MeshData& mesh, SurfaceStateData& state,
                          const SolverOptions2D& opts);
 
} // namespace openswmm::twoD
 
#endif // OPENSWMM_ENGINE_2D_SURFACE_FLUX_CALCULATOR_HPP