Edge flux computation, gradient calculation, and slope limiting. More...

#include "../data/MeshData.hpp"
#include "../data/SurfaceStateData.hpp"
#include "../data/SolverOptions2D.hpp"

Include dependency graph for SurfaceFluxCalculator.hpp:

This graph shows which files directly or indirectly include this file:

Namespaces
namespace	openswmm

namespace	openswmm::twoD

Functions
void	openswmm::twoD::computeUnlimitedGradients (const MeshData &mesh, SurfaceStateData &state, int nthreads=1)
	Compute unlimited gradients for all triangles via Green-Gauss theorem.

void	openswmm::twoD::computeLimitedGradients (const MeshData &mesh, SurfaceStateData &state, double epsilon, int nthreads=1)
	Apply Jawahar-Kamath slope limiter (Eq. [23]–[24]).

void	openswmm::twoD::computeEdgeFluxes (const MeshData &mesh, SurfaceStateData &state, const SolverOptions2D &opts)
	Compute edge fluxes for all triangles.

double	openswmm::twoD::evapSink (double rate, double depth, double dry_depth) noexcept
	Depth-limited evaporation sink rate (m/s) for one cell.

void	openswmm::twoD::assembleRHS (const MeshData &mesh, const SurfaceStateData &state, const SolverOptions2D &opts, double *ydot)
	Assemble the RHS of the ODE system: dψ/dt for each triangle.

void	openswmm::twoD::assembleImplicitRHS (const MeshData &mesh, const SurfaceStateData &state, const SolverOptions2D &opts, double *ydot)
	Implicit half of the IMEX split: the flux divergence only (ARKODE F_I).

void	openswmm::twoD::assembleExplicitRHS (const MeshData &mesh, const SurfaceStateData &state, const SolverOptions2D &opts, const double y, double ydot)
	Explicit half of the IMEX split: source/sink forcing (ARKODE F_E).

void	openswmm::twoD::computeCellContinuity (const MeshData &mesh, SurfaceStateData &state, const SolverOptions2D &opts, double dt)
	Compute the per-cell continuity residual (local mass-balance check).

void	openswmm::twoD::computeFaceVelocity (const MeshData &mesh, SurfaceStateData &state, const SolverOptions2D &opts)
	Reconstruct cell-centred velocity (vx, vy) from edge fluxes (RT0).

Detailed Description

Edge flux computation, gradient calculation, and slope limiting.

Implements the second-order finite volume flux computation:

Unlimited gradients via Green-Gauss (Eq. [25]–[26])
Slope limiter (Jawahar-Kamath, Eq. [23]–[24])
Edge flux with upwind selection (Eq. [15a], [22])

See also: TWO_DIMENSIONAL_SURFACE_ROUTING_IMPLEMENTATION_STRATEGY.md §3, §10

Author: Caleb Buahin caleb.nosp@m..bua.nosp@m.hin@g.nosp@m.mail.nosp@m..com

License\n MIT License

Namespaces

Functions

Detailed Description