A Visual Software Implementation of Numerical Simulation for Transient Process of Pipeline Network System of Water Supply Project

: The transient process is the key for the design, operation and maintenance of the pipeline network system of the water supply project. This paper aims to study and develop a new type of visual numerical simulation software for the transient process of pipeline network system. Firstly, the software architecture design is illustrated. Then, three tiers of software architecture, i.e., back-tier, middle-tier and front-tier, are studied and developed. Finally, the practical application procedure of the software is illustrated and the accuracy of calculation results is verified. The results indicate that the visual numerical simulation software is designed based on C/S architecture. The software architecture contains back-tier, middle-tier and front-tier. The back-tier includes the model and the algorithm. The middle-tier realizes the decoding of the topology of the pipeline network system and the interaction between the back-tier and the front-tier. The front-tier integrates the visual interface and pre-post processing. The main interface includes a menu bar, a visual modeling area and a parameter input dialog box. The software has the features of professionalization, visual-ization, generalization, modularization and intellectualization. The software can calculate the transient process for the pipeline network system of both, pressure flow without the pump and pressure flow with the pump. The physical laws reflected from the calculation results are correct. The calculation results are accurate.


Introduction
The water supply project is an infrastructure for the transmission and distribution of water resources for living and production [1,2]. Generally, a water supply project contains reservoirs, pipelines, pumps, valves and water hammer prevention facilities. The types of water supply projects usually include a pressure flow system without a pump, a pressure flow system with a pump and a mixed water supply system [3,4].
The water supply project generally consists of a complex pipeline network system. Based on the theories of transient flow, the hydraulic design and optimization of the pipeline network system are the basis for the construction and operation of water supply projects [5,6]. Moreover, the hydraulic design and optimization are based on the calculation and analysis of transient processes [7][8][9]. The transient process is the key aspect for the design, operation and maintenance of the pipeline network system of the water supply project. According to the calculation and analysis results of the transient process, the design guidelines, optimization methods and operation strategies can be obtained for the pipeline network system. man-machine interaction become higher. The software is expected to be more intelligent and has a powerful pre-post processing function. Moreover, several new types of valves and water hammer prevention facilities are applied into water supply projects. The visual numerical simulation software should follow the development trend of software design and water supply projects. Based on the above considerations, a new type of visual numerical simulation software for transient process of pipeline network system of water supply project needs to be studied and developed.
This paper aims to study and develop a new software for graphical and numerical simulations of the transient processes in the water supply network. It is expected that the new software can realize the efficient calculation for transient processes in the water supply network and contains the following features: professionalization, visualization, generalization, modularization and intellectualization. Moreover, as the novelty and innovation of the software, more friendly man-machine interaction, more powerful pre-post processing function and new types of valves and water hammer prevention facilities are focused and achieved.
The paper structure is as follows. In Section 2, software architecture design is illustrated. In Sections 3-5, three tiers of the software architecture, i.e., back-tier, middle-tier and front-tier, are studied and developed, respectively. In Section 6, two engineering examples, i.e., a pipeline network system of pressure flow without a pump and a pipeline network system of pressure flow with a pump, are selected to illustrate the practical application of software and verify the accuracy of calculation results. In Section 7, the conclusions are given.

Software Architecture Design
The visual numerical simulation software is designed based on the client/server architecture (C/S architecture). The client realizes the interaction tasks with the user, and the server realizes the management of data. The interfaces and operations of C/S architecture are abundant. The interaction and response speed of C/S architecture are rapid. The useroriented system design realizes the interactive information service between man and the machine. The software architecture contains three tiers, i.e., back-tier, middle-tier and front-tier. The connotations and relationships of the three tiers of software architecture are shown in Figure 1. (1) Back-tier The key of back-tier is the model and the algorithm. The models of the components of the pipeline network system are established. The suitable algorithms are adopted for the solving of models. The models and algorithms are made into general calculation modules. For the numerical simulation for the transient process of the pipeline network system of the water supply project, the operating condition should be set first. Then, under an operating condition, the general calculation modules are selected and called. After the instruction of calculation is input, the calculation of transient process starts and lasts in the set time.
(2) Middle-tier The key to middle-tier is the decoding and interaction. The topology of pipeline network system is decoded firstly. After the topology decoding, the properties and connected relations of different graphical calculation modules are determined. Then, for different graphical calculation modules, the parameters are input, read and transmitted. The data reading and transmission exist in the back-tier, middle-tier and front-tier. The data management interface realizes the data reading and transmission rapidly and smoothly.
(3) Front-tier The key to front-tier is the visual interface and pre-post processing. The visual interface of the visual numerical simulation software is integrated in the front-tier. As a general, visual and friendly numerical simulation software, the visual interface is composed of a menu bar, a visual modeling area and a parameter input dialog box. By using the functions of the menu bar, the visual modeling area and the parameter input dialog box, the models and calculation conditions of the pipeline network system can be established and set. After the calculation of the transient process, the post-processing of results is realized in the front-tier. The calculation results of the transient process can be displayed, analyzed and saved.
Under the combination of the three tiers, the visual numerical simulation software can realize the efficient calculation for the transient process of the pipeline network system of the water supply project. The software is developed and integrated in Java. The numerical simulation software contains the following features: professionalization, visualization, generalization, modularization and intellectualization. Specifically, (1) Professionalization: the models of the pipeline network system are extracted from authoritative literature. The classical algorithms for the calculation of the transient process are adopted. Moreover, several new types of valves and water hammer prevention facilities are considered and included. The transient processes in the water supply network can be simulated in detail, and comprehensively.
(2) Visualization: the software is a user-oriented system to realize the interactive information service between man and a machine. The visual interface of the software integrates all the operating menus and functions. In the visual interface, the model of the pipeline network system is established by graphical modules and interface tools.
(3) Generalization: the application of the software is not limited by the types of the water supply project, the layouts of the pipeline network system or the operating conditions of the transient process. The software is general for the numerical simulation for the transient processes in the water supply network.
(4) Modularization: in software, the models and algorithms are made into general modules. The components of the water supply network, such as pipe, pump and valve, are made into general model modules. The operating conditions of the transient process, such as starting up and shutting down, are made into general calculation modules.
(5) Intellectualization: the software contains several intelligent functions that make the use more convenient and efficient. By importing the CAD file of water supply network, the software can automatically identify the topology of the pipeline network system and then establish the graphical model. The calculation results of the transient process can be output as a formal report. In that report, the calculation results are automatically collected and classified. The changing curves of variables are automatically drawn and inserted.

Model of Unsteady Flow in Pipelines and Method of Characteristics
The research object of the pipeline network system of the water supply project is the water flow in the pipelines. In the present study, the pipeline network system is the pressure flow and completely filled with water. There is no free water surface in the pipelines. During the transient process, the water flow is unsteady. It is assumed that the water in the pipeline network system is incompressible. The basic equations that describe the unsteady flow in the pipelines contain the momentum equation and continuity equation [5,6]. Those two basic equations are presented in Equation (1).
where Q is the discharge in the pipeline, H is the pressure head, x is the distance along the pipeline axis, D is the diameter of the pipeline, A is the cross-sectional area of the pipeline, θ is the angle between the pipeline axis and horizontal line, f is the Darcy-Weisbach friction coefficient, a is the speed of water hammer wave, g is the acceleration of gravity, t is time.
The most widely used method for solving Equation (1) is the method of characteristics [5,6]. The method of characteristics is a kind of one-dimensional numerical solution. By using the integration along the characteristic lines, the ordinary differential equations are transformed into finite difference equations. The finite difference equations are algebraic equations and can be easily solved. The schematic diagram of the method of characteristics for solving Equation (1) is shown in Figure 2. The abscissa axis of Figure 2 represents the pipeline axis. The pipeline axis is divided into many micro sections with a length of x ∆ . The vertical axis of Figure 2 represents the duration time of transient process. t ∆ is the time step of difference calculation. There are two characteristic lines, i.e., C + and C − . The characteristic line C + is from point A to point P, the characteristic line C − is from point B to point P.
The equations of characteristic lines are: where V is the flow speed in the pipeline.
After the transformation and integration of Equation (1) on Equation (2), the following finite difference equations are obtained. : : where 1 A C and Based on Equation (3), the discharge and pressure head of point P, i.e., P Q and P H , can be solved by using the discharges and pressure heads of points A and B. By combining the boundary conditions and initial conditions, the discharges and pressure heads of all cross sections of pipeline network system during the transient process can be solved.

Models of Other Network Elements
The other elements of pipeline network system of water supply project include pump, valve and water hammer prevention facility. In the present study, the models of other network elements are established by matching with the method of characteristics.

Pump
Pump is the pressurized equipment of the pipeline network system of the water supply project [38,39]. There are many kinds of pump. The model of constant speed pump is introduced in the following paragraphs. Figure 3 is the schematic diagram of the constant speed pump. The model of pump is composed of the pressure head equilibrium equation and inertia equation [5,9]. The pressure head equilibrium equation describes the balance of pressure head between the inlet section of the pump and the output section of the valve. The valve locates at the near right side of pump. The pressure head equilibrium equation of the pump is 2 1 where H is the output head of the pump,    : : where 1 P C and f H can be determined by the following equation [5,9,38,39].
where 0 f H is the head loss of valve when the discharge of the pump is the rated value, is the relative opening of the valve, θ and R θ are the opening angle and rated opening angle of the valve, respectively, is the relative discharge of the pump.
The inertia equation of the pump is [5,9,38,39] ( ) where R n n α = is the relative rotational speed of the pump, n and R n are the rotational speed and rated rotational speed of the pump, respectively, is the relative torque of the pump, T and R T are the torque and rated torque of the pump, respectively, the subscript t ∆ represents the value of the variable after t ∆ of the present calculation time, 0 C and 1 C are the coefficients, k is an intermediate variable. The models of other kinds of pump, such as variable speed the pump, can be found in [5,9,38,39].

Valve
Valve is the discharge regulation equipment of the pipeline network system of the water supply project [9,40]. There are many kinds of valves. The model of gate valve is introduced in the following paragraphs. Figure 4 is the schematic diagram of the gate valve. The discharge through the valve is [9,40]   : : When the opening of valve is R θ , the discharge through the valve is R Q . Based on where the subscript R represents the rated value of the variable. By using Equations (8) and (10), we can get The models of other kinds of the valve, such as regulating valve, hydraulic control valve, air valve, check valve and safety valve can be found in [9,40].

Water Hammer Prevention Facility
Water hammer prevention facility is the pressure reduction equipment of the pipeline network system of the water supply project [8,9,40]. There are many kinds of water hammer prevention facilities. The model of the air tank is introduced in the following paragraphs. Figure 5 is the schematic diagram of the air tank. The air in the top of the tank satisfies the following air state equation [8,9,40,41].
where p is the pressure of the air in tank, l is the height of the air cushion, m is the polytropic exponent of air and C is the constant.
At the bottom of the air tank, the discharges T Q , S Q and R Q satisfies the following continuity equation.
where T Q is the discharge in the pipe before the air tank, S Q is the discharge that flows into the air tank and R Q is the discharge in the pipe after the air tank.
The relationship among S Q , l and the water level S H in air tank is [8,9,40,41] The piezometric head at the bottom of the air tank, i.e., point S, is denoted as Then we have where γ is the unit weight of water, H ∆ is the standard atmospheric pressure and A α is the discharge coefficient of head loss at the throttling orifice of air tank.
The head loss at the bottom of the air tank is neglected. Therefore, we have  : : The models of other kinds of water hammer prevention facility, such as two-way surge tower and one-way surge tower, can be found in [8,9,40,41].

Boundary and Initial Conditions
The boundary conditions of pipeline network system of water supply project are the reservoirs, which include upstream reservoir and downstream reservoir [5,6]. The equations of upstream reservoir and downstream reservoir are introduced in the following paragraphs. Figure 6 is the schematic diagram of the upstream reservoir and downstream reservoir.   The initial conditions are determined by the initial steady state of pipeline network system of water supply project. For an operating condition, the initial state for the calculation of transient process of pipeline network system is set as the steady state. The values of variables under the initial steady state of pipeline network system are the initial conditions for the calculation of transient process.

Time and Space Discretization
The time and space discretization is about the determination of time step t ∆ and space step x ∆ for the calculation of transient process. Based on the theory of the method of characteristics, the time step t ∆ and space step x ∆ should satisfy the following stability condition [5,6].
In the visual numerical simulation software, we take

Middle-Tier: Decoding and Interaction
The task of middle-tier contains the decoding of topology of the pipeline network system and the interaction between the back-tier and the front-tier. The key technological difficulty is the decoding of topology of the pipeline network system, which is illuminated in this section.
The topology decoding is composed of two steps, i.e., automatic encoding and automatic pipeline identification.

Automatic Encoding
The topology of the pipeline network system is composed of pipelines and nodes. The purpose of automatic encoding is to determine the connection relation between pipelines and nodes. That connection relation should be identified from the graphical model of pipeline network system and expressed by data structure with a specific form. Then the data structure should be transmitted to the computer.
Based on the feature of the topology of pipeline network system, the automatic encoding should realize the following two functions. Firstly, the pipelines and nodes are numbered, respectively. Then, the connection relation among pipelines, upstream nodes and downstream nodes are determined. The automatic encoding is realized by the method of object-oriented programming. The requirements of automatic encoding include: (a) All the data structures of encoding are generated by the computer automatically. In the method of object-oriented programming, the pipelines and nodes are regarded as objects. Three data members are added for each object. That data member is called the encoding data member. Among the three data members for an object, the first encoding data member represents the encoding of the object. The second encoding data member represents the upstream connection relation, i.e., the encoding of the upstream pipeline or node of the object. The third encoding data member represents the downstream connection relation, i.e., the encoding of the downstream pipeline or node of the object. Because every object has those three encoding data members, those three encoding data members can be defined in base class.
The encoding is dynamic. The serial number of the object in data structure is regarded as the encoding of the object. With the change of the graphical model of pipeline network system, the serial number of the object in data structure changes accordingly. Then the dynamic encoding is realized. During the calculation of transient process, the data reading and transmission mainly exist in the pipelines, upstream nodes and downstream nodes.

Automatic Pipeline Identification
The automatic encoding determines the connection relation between pipelines and nodes. Based on the result of automatic encoding, the automatic pipeline identification aims to arrange the pipelines and nodes in the order from upstream to downstream. All the pipelines are divided into two parts, i.e., upstream pipelines and downstream pipelines.
In order to realize the function of automatic pipeline identification, two data members are added for each pipeline. For a pipeline, the first data member is used to store the information of its upstream pipeline, and the second data member is used to store the information of its downstream pipeline. The automatic pipeline identification contains two aspects.
Aspect 1-from upstream to downstream: in the process of automatic pipeline identification, the pipelines connected with upstream reservoir are identified firstly. Those pipelines are denoted as the starting pipelines. According to the starting pipelines, the searching is carried out along pipelines from upstream to downstream. For a pipeline, all the connected pipelines are identified. The searching stops when the upstream pipeline connected with pump is identified.
Aspect 2-from downstream to upstream: in the process of automatic pipeline identification, the pipelines connected with downstream reservoir are identified firstly. Those pipelines are denoted as the starting pipelines. According to the starting pipelines, the searching is carried out along pipelines from downstream to upstream. For a pipeline, all the connected pipelines are identified. The searching stops when the downstream pipeline connected with pump is identified.
The procedure of automatic pipeline identification is shown in Figure 7.

Front-Tier: Visual Interface and Pre-Post Processing
Following the principles of practicability, concision, beauty and friendly man-machine interaction, the main interface of visual numerical simulation software for transient process of pipeline network system of water supply project is designed as Figure 8. The main interface includes three parts, i.e., menu bar, visual modeling area and parameter input dialog box. The main functions and modules are included in the menu bar. The visual modeling area is used to achieve the graphical modeling of pipeline network system. The practical and three-dimensional pipeline network system of the water supply project is converted into two-dimensional topology. The topology is composed of general modules. The parameter input dialog box is used to input the attribute parameters of modules. The operating condition is also set in the parameter input dialog box.
The main functions and modules are emphatically introduce based on the design of the menu bar in the following paragraphs. The two-stage menus are designed for the menu bar. The first stage menu contains six items, i.e., file, modeling, calculation, result view, report output and help. The functions of six items of the first stage menu are shown in Table 1. Introduce the functions of software, provide the tutorial.
The detailed items and functions of the second stage menu are presented in Table 2. After the graphical model of the pipeline network system is established, the attribute parameters of models can be input based on the general modules. By double clicking a module, the parameter input dialog box pops up. For different modules, the different attribute parameters are needed. After the input of attribute parameters, the parameters are saved in modules and can be called for the calculation of the transient process.
The characteristic curves of pump are input in the module of pump. Two methods for the input of characteristic curves of pump are provided in the software. The first method is the manual input of characteristic curves. The user prepares the file of characteristic curves of pump in the format of Excel. Then by inputting the file into the pump module, the characteristic curves of discharge and characteristic curves of moment can be presented in the software. The second method is the automatic generation of characteristic curves. A database of the characteristic curves of pump is embedded in the software. For a pump, the characteristic curves can be obtained from the database based on its specific speed and curve interpolation.

Engineering Examples
After the combination of back-tier, middle-tier and front-tier, the visual numerical simulation software can realize the efficient calculation for transient process of pipeline network system of water supply project. For a pipeline network system, the procedure of the calculation for transient process by using the software includes the following three steps.
Step 1: The actual layout of pipeline network system is divided into several pipelines and nodes. Then the pipeline network system can be abstracted into a two-dimensional topology. Based on the topology, the graphical model of pipeline network system is established in visual modeling area by using the general modules. For all the modules, the attribute parameters are input. As a result, a complete graphical model of pipeline network system is obtained. The graphical model reflects the layout, structure and characteristic of pipeline network system.
Step 2: The operating condition is set in the modules of pumps and valves. The calculation parameters and conditions are set in the first stage menu of calculation. After all the parameters of operating condition and calculation conditions are set, the calculation of transient process can start.
Step 3: After the calculation of transient process, the results can be saved in the format of Excel. The results of steady flow and unsteady flow can be viewed in the software directly.
The changing curves of variables can be drawn based on the calculation results of unsteady flow. According to the results and curves, the characteristics of transient process of pipeline network system can be analyzed and evaluated. The calculation results of transient process can also be output as a formal report.
In the following parts of this section, two engineering examples, i.e., pipeline network system of pressure flow without pump and pipeline network system of pressure flow with pump, are selected to illustrate the practical application of software and verify the accuracy of calculation results. Both of those two pipeline network systems are pressure flow, and there is no free water surface in the pipelines. The two engineering examples are actual water supply projects. The calculation results of the present software are compared with those of one existing software, i.e., HAMMER [33,34].
HAMMER is a widely used commercial software for the numerical simulation for transient process of pipeline network system of water supply project. The modeling theory and method of HAMMER are the same with those of the present software. Moreover, the solution method of HAMMER is also the method of characteristics. However, for the treatment of many technical and programming details, the two kinds of software are different. For example, (1) For the calculation of transient process, the characteristic curves of pump and valve need to be interpolated because the data points of original characteristic curves are very limited. In the present software, the interpolation of characteristic curves is realized by the method of Lagrange interpolation. However, the interpolation method of HAM-MER is not public. Therefore, the interpolation methods of characteristic curves of pump and valve for the two kinds of software cannot keep the same.
(2) For the calculation of transient process, the equations of pump should be solved together with the equations of pipelines based on iteration. In the present software, the iteration solution is realized by the method of Newton-Raphson iteration. However, the iteration method of HAMMER is not public. Therefore, the iteration methods and accuracies of equations of pump and pipelines for the two kinds of software cannot keep the same.
The present software and HAMMER are large and complex. Besides the basic models, solution methods and initial conditions, there are lots of technical and programming details that would lead to the calculation deviations between the two kinds of software. Because many technical and programming details of HAMMER are not public, it is difficult to make all the calculation conditions of the two kinds of software the same. The calculation deviations between the two kinds of software are difficult to avoid. As a new software, the calculation results of the present software remain to be tested. As a widely used commercial software, the calculation results of HAMMER are credible and can be regarded as a reference. Hence, in order to verify the rationality of treatment of details and accuracy of calculation results of the present software, the contrastive analysis is necessary.

Pipeline Network System of Pressure Flow without Pump
The pipeline network system of pressure flow without pump selected in this section is designed for water transfer. From the upstream reservoir to downstream reservoir, the length of the whole pipeline is about 71.07 km. At the end of the pipeline network system, two regulating valves are set in parallel. One regulating valve is used for normal operation, and the other regulating valve is used for emergency reserve. According to the actual layout of pipeline network system of pressure flow without pump, the graphical model of pipeline network system is established by present software and shown in Figure 9. In Figure 9, N1 is the upstream reservoir, N11 and N12 are the downstream reservoirs. P1-P13 are the pipelines. N2-N7 and N14 are the tandem nodes. N9 and N10 are the two regulating valves. N8 is the bifurcated point. N13 is the two-way surge tower. The values of parameters of pipelines are shown in Table 3. Four operating conditions are selected for the calculation and comparison. The information of the four operating conditions is shown in Table 4. Under the four operating conditions in Table 4, the transient process for the selected pipeline network system of the pressure flow without the pump can be calculated by using the present software and the HAMMER. For each operating condition, four variables, i.e., pressure head at the inlet section of valve, pressure head at the output section of valve, discharge through the valve and water level in the two-way surge tower are presented and compared. The comparisons of calculation results under GF1, GF2, GF3 and GF4 are shown in Figures 10-13, respectively. Software 1 represents the present software and software 2 represents the HAMMER. (2) Under the four operating conditions, the maximum relative error of the extremum of pressure head at the inlet section of valve is 3.68% and occurs in GF4. The maximum relative error of the extremum of pressure head at the output section of valve is 4.98% and occurs in GF3. The maximum relative error of the extremum of water level in the two-way surge tower is 4.99% and occurs in GF4. The above results indicate that the errors of calculation results between software 1 and software 2 are small. Therefore, the calculation results of software 1 are accurate.
To sum up, the software 1 can calculate the transient process for the pipeline network system of pressure flow without pump reasonably and accurately.

Pipeline Network System of Pressure Flow with Pump
The pipeline network system of pressure flow with pump selected in this section is designed for water supplement of reservoir. The length of the whole pipeline is about 2.65 km. The water supplement from downstream reservoir to upstream reservoir is realized by one pump. According to the actual layout of pipeline network system of pressure flow with pump, the graphical model of pipeline network system is established by present software and shown in Figure 14. In Figure 14, N1 is the downstream reservoir, and N8 is the upstream reservoir. P1-P7 are the pipelines. N2 is the pump. N3-N6 are the air valves. N7 is the one-way surge tower. The values of parameters of pipelines are shown in Table 5. The basic data of the pump are shown in Table 6. Three operating conditions are selected for calculation and comparison. The information of the three operating conditions is shown in Table 7. Under the three operating conditions in Table 7, the transient process for the selected pipeline network system of pressure flow with pump can be calculated by using the present software and the HAMMER. For each operating condition, one variable, i.e., pressure head at the output section of valve N3, is presented and compared. The comparisons of calculation results under PF1, PF2 and PF3 are shown in Figures 15-17, respectively. Software 1 represents the present software, and software 2 represents the HAMMER.    are consistent with those of software 2. Therefore, the physical laws reflected from the calculation results of software 1 are correct.
(2) Under the three operating conditions of PF1, PF2 and PF3, the maximum relative errors of the extremum of pressure head at the output section of valve N3 are 2.41%, 2.02% and 2.60%, respectively. The above results indicate that the errors of calculation results between software 1 and software 2 are small. Therefore, the calculation results of software 1 are accurate.
To sum up, the software 1 can calculate the transient process for the pipeline network system of pressure flow with pump reasonably and accurately.

Summary and Conclusions
The visual numerical simulation software for transient process of the pipeline network system of the water supply project is studied and developed. The practical application of software is illustrated, and the accuracy of calculation results is verified. The conclusions are: (1) The visual numerical simulation software is designed based on the C/S architecture. The software architecture contains three tiers, i.e., back-tier, middle-tier and fronttier. Under the combination of the three tiers, the visual numerical simulation software can realize the efficient calculation for transient process of pipeline network system of water supply project. The numerical simulation software contains the following features: professionalization, visualization, generalization, modularization and intellectualization.
(2) The method for solving the basic equations of unsteady flow in the pipelines is the method of characteristics. The boundary conditions of pipeline network system of water supply project include pump, valve, water hammer prevention facility and reservoir. The task of middle-tier contains the decoding of topology of pipeline network system and the interaction between back-tier and front-tier. The topology decoding is composed of two steps, i.e., automatic encoding and automatic pipeline identification. The main interface includes three parts, i.e., menu bar, visual modeling area and parameter input dialog box.
(3) The visual numerical simulation software can calculate the transient process for the pipeline network system of both pressure flow without pump and pressure flow with pump reasonably and accurately. The physical laws reflected from the calculation results of software are correct. The calculation results of software are accurate.
The following points should be noted: (1) In the current version of software, the texts and words of the visual interface are in Chinese. An English version will be updated in the future.
(2) At present, the visual numerical simulation software is developed and used by the authors and their teams. The ranges of application of the software include scientific research, technology consulting, and engineering design. The software is not open source. In near future, we aim to promote it to be a commercial software and available for the general public.
(3) As space is limited, only two engineering examples are presented and one existing software, i.e., HAMMER, is compared. The calculation results of software are not compared with the measured data because the measured data for transient process of pipeline network system of water supply project are scarce.
(4) In the visual numerical simulation software, there is not a special module for flowmeter. If the flowmeter has no influence on the flow state of water in the pipeline and has no head loss, such as electromagnetic flowmeter and ultrasonic flowmeter, the flowmeter can be equivalently regarded as a part of pipeline and then simulated in the software. If the flowmeter has head loss, such as differential pressure flowmeter and thermal time of flight flow meter [42], the flowmeter cannot be simulated in the software.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The author declares no conflict of interest.