Numerical network models and entropy principles for isothermal junction flow

We numerically explore network models which are derived for the isothermal Euler equations. Previously we proved the existence and uniqueness of solutions to the generalized Riemann problem at a junction under the conditions of monotone momentum related coupling constant and equal crosssectional areas for all connected pipe sections. In the present paper we extend this proof to the case of pipe sections of different cross-sectional areas. We describe a numerical implementation of the network models, where the flow in each pipe section is calculated using a classical high-resolution Roe scheme. We propose a numerical treatment of the boundary conditions at the pipe-junction interface, consistent with the coupling conditions. In particular, mass is exactly conserved across the junction. Numerical results are provided for two different network configurations and for three different network models. Mechanical energy considerations are applied in order to evaluate the results in terms of physical soundness. Analytical predictions for junctions connecting three pipe sections are verified for both network configurations. Long term behaviour of physical and unphysical solutions are presented and compared, and the impact of having pipes with different cross-sectional area is shown.