G-steps 4 and 5 — see the header for the coordinate, the flux forms and the handover convention. More...

#include "SigmaColumn.hpp"
#include <algorithm>
#include <cmath>
#include <vector>

Include dependency graph for SigmaColumn.cpp:

Functions
double	openswmm::twoD::sigma::columnStorage (const double *theta, int m, int stride, double L) noexcept
	Total water in the column (m of water) — `Σ θ_j·L·Δσ`.

void	openswmm::twoD::sigma::seedHydrostatic (const soil::Params &p, double *theta, int m, int stride, double L) noexcept
	Seed a column to hydrostatic equilibrium above a table at depth `L`.

double	openswmm::twoD::sigma::columnDtLimit (const soil::Params &p, const double *theta, int m, int stride, double L, double Ldot, double c_col, bool capillary) noexcept
	The explicit stability bound of a column, without advancing it (§2.3): `Δt_u ≤ C_col·min_j[ L·Δσ / ( c(θ_j) + \|σ_j·L̇\| + 2·D̄_j/(L·Δσ) ) ]` with celerity `c(θ) = dK/dθ`.

void	openswmm::twoD::sigma::prepareColumn (const soil::Params &p, const double *theta, int m, int stride, const ColumnStep &st, ColumnPreparation &prepared) noexcept
	Prepare reusable trial invariants. Same input preconditions as advanceColumn; m > 0.

void	openswmm::twoD::sigma::advanceColumn (const soil::Params &p, double *theta, int m, int stride, ColumnStep &st) noexcept
	Advance one σ column by `st.dt`.

void	openswmm::twoD::sigma::advanceColumn (const soil::Params &p, double *theta, int m, int stride, ColumnStep &st, const ColumnPreparation &prepared) noexcept

Detailed Description

G-steps 4 and 5 — see the header for the coordinate, the flux forms and the handover convention.

Namespaces