| Summary: | On the basis of the two-component pressure approach, we developed a numerical model to capture mixed transient flows in close conduit systems. To achieve this goal, an innovative Godunov finite-volume numerical scheme is proposed to suppress the spurious numerical oscillations occurring during rapid pipe pressurization. To dissipate the spurious numerical oscillations, we admit artificial numerical viscosity to the numerical scheme through applying a proposed Harten, Lax, and van Leer (HLL) Riemann solver for calculating the numerical fluxes at the computational cell interfaces. The proposed solver controls the magnitude of the numerical viscosity through adjusting the left and right wave velocities. A wave velocity calculator is proposed to optimally distribute the numerical viscosity over several computational cells around the computational cell in which the pressurization front is located. The proposed solver admits significant artificial numerical viscosity when the pipe pressurization is imminent and automatically reduces it in other places; in this way the numerical diffusion and data smearing is minimized. The validity of the proposed model is justified by the aid of several test cases in which the numerical results are compared with both experimental data and the results obtained from analytical methods. The results reveal that the proposed model succeeds in completely removing the spurious numerical oscillations, even when the pipe acoustic speed is over 1000 m/s. The numerical results also show that the model can successfully capture occurrence of negative pressures during the course of transient flow.
|
|---|
