Simulation of transcritical flow in pipe/channel networks

Using finite difference methods in conjunction with the reduced momentum equation and applying boundary condition structure inherent to subcritical flow to all regimes, is an approach that enables efficient numerical simulation of supercritical and transcritical flows in pipe/channel systems. However, as well as certain errors within a single channel due to incomplete equations, this technique also may introduce unwanted effects propagating across a network in both upstream and downstream directions. These may include: unrealistic backwater effects due to improper boundary conditions, nonamplifying oscillations due to jerky jump movement, and other computational instabilities. Practical implications of these are analyzed in detail and are illustrated using a set of examples. Sensitivity analyzes and comparisons with analytical solutions and laboratory experiments are made. The measures to reduce the inaccuracies inevitable in simulation of transcritical flows are discussed.