Implementation of a hydraulic check valve. Note that the small reverse flows can still occur with this model.
The basic flow function
ṁ = sign(Δp) k √ Δp ,
with regularization near the origin, is used to compute the pressure drop. The flow coefficient
k = ṁ ⁄ √ Δp
is increased from l*KV_Si to KV_Si,
where KV_Si is equal to Kv but in SI
units. Therefore, the flow coefficient k is set to a
value close to zero for negative pressure differences, thereby
restricting reverse flow to a small value. The flow coefficient
k saturates to its maximum value at the pressure
dpValve_closing. For larger pressure drops, the
pressure drop is a quadratic function of the flow rate.
The parameters m_flow_nominal and
dpValve_nominal determine the flow coefficient of the
check valve when it is fully opened. A typical value for a nominal
flow rate of 1 m/s is dpValve_nominal = 3400
Pa. The leakage ratio l determines the minimum
flow coefficient, for negative pressure differences. The parameter
dpFixed_nominal allows to include a series pressure
drop with a fixed flow coefficient into the model. The parameter
dpValve_closing determines when the flow coefficient
starts to increase, which is typically in the order of
dpValve_nominal.
The check valve implementation approximates the physics where a
forward pressure difference opens the valve such that the valve
opening increases, causing a growing orifice area and thus
increasing the flow coefficient. Near
dp=dpValve_closing, the valve is fully open and the
flow coefficient saturates to the flow coefficient value determined
by dpValve_nominal and m_flow_nominal.
For typical valve diameters, the check valve is only fully open
near nominal mass flow rate. Therefore, the model sets
dpValve_closing=dpValve_nominal/2 by default.