Everything the residual needs, gathered so the Newton can stay generic. More...

#include <line/solvers/fluid/fluid_dae.h>

Collaboration diagram for line::fluid::FluidDaeSystem:

Public Member Functions
std::size_t	nstaging () const
	u = [x; staging; sigma2; mult]

Public Attributes
const FluidMomentTerms *	terms = nullptr
const FluidDaeConservation *	cons = nullptr
const FluidDaeConstraints *	con = nullptr
const FluidDaeStaging *	stg = nullptr
const FluidDaeGates *	gates = nullptr
std::vector< std::size_t >	cidx
	closable stations
std::vector< std::size_t >	active
	binding capacity constraints
std::size_t	nstate = 0
Matrix< double >	clampT = Matrix<double>(0, 0, 0.0)
	The tangent space of the caps that CLAMP, or an empty matrix.

Detailed Description

Everything the residual needs, gathered so the Newton can stay generic.

Definition at line 947 of file fluid_dae.h.

Member Function Documentation

◆ nstaging()

std::size_t line::fluid::FluidDaeSystem::nstaging ( ) const

inline

u = [x; staging; sigma2; mult]

Definition at line 967 of file fluid_dae.h.

References line::fluid::FluidDaeStaging::n, and stg.

Referenced by line::fluid::fluid_dae_residual().

Member Data Documentation

◆ active

std::vector<std::size_t> line::fluid::FluidDaeSystem::active

binding capacity constraints

Definition at line 954 of file fluid_dae.h.

Referenced by line::fluid::fluid_dae_residual(), and line::fluid::solver_fluid_dae().

◆ cidx

std::vector<std::size_t> line::fluid::FluidDaeSystem::cidx

closable stations

Definition at line 953 of file fluid_dae.h.

Referenced by line::fluid::fluid_dae_residual(), and line::fluid::solver_fluid_dae().

◆ clampT

Matrix<double> line::fluid::FluidDaeSystem::clampT = Matrix<double>(0, 0, 0.0)

The tangent space of the caps that CLAMP, or an empty matrix.

A cap that holds the job upstream or loses it fixes its own combination of the state for as long as it binds, so that combination does not fluctuate: WITHOUT the projection the Lyapunov solve has no solution at all, because with the multiplier held constant the drift along the constrained direction is neutral. A STAGED cap is not clamped – its room is a state and supplies the restoring force – and must not be projected.

Definition at line 965 of file fluid_dae.h.

Referenced by line::fluid::fluid_dae_residual(), and line::fluid::solver_fluid_dae().