Rapid flow net modeling and simulation method for supercritical carbon dioxide Brayton cycle

Abstract

The invention discloses a supercritical carbon dioxide Brayton cycle-oriented flow net rapid modeling simulation method, which comprises the following steps of: constructing an input card, and reading parameters of a flow net system; based on the equipment nodal part of the input card, carrying out node division on a flow network system, and generating a connection relation matrix; constructing a pressure matrix and an enthalpy value matrix of the flow net based on the three conservation equations and the calculation result of the previous time step length; solving a pressure matrix and an enthalpy value matrix through a parallel algorithm; and a calculation result is output in real time, and simulation of different speeds is realized by setting a super real-time multiple. According to the method, various arrangements and parameter settings can be carried out, and the modeling efficiency of technicians can be improved through a carbon dioxide Brayton cycle flow net program which can be modeled through different equipment combinations; and real-time and super-real-time simulation can be carried out, so that the purposes of design and teaching are achieved.

Description

Technical Field

[0001] The present invention belongs to the technical field of real-time simulation of nuclear power plants, and particularly relates to a method for quickly modeling and simulating a flow network for a supercritical carbon dioxide Brayton cycle. Background Technique

[0002] The supercritical carbon dioxide Brayton cycle system is a nuclear power plant loop system based on supercritical carbon dioxide. This system utilizes the unique properties of carbon dioxide under supercritical conditions, such as high density and low viscosity, to achieve higher energy conversion efficiency and smaller equipment size. Currently, there are already multiple schemes for the layout of the supercritical carbon dioxide Brayton cycle, and its simulation research mainly focuses on the exploration and comparison of cycle characteristics, efficiency, and control strategies. In the field of simulation, traditional water reactors usually adopt a method combining experiments and calculations, analyzing operating characteristics, discovering the mechanism of phenomena, and drawing conclusions to optimize system design and operating strategies.

[0003] However, the operating pressure of the supercritical carbon dioxide Brayton cycle system is relatively high, and the development of related equipment is not yet mature, resulting in relatively few experimental studies. Therefore, current research mostly relies on simulation means. When facing different usage environments or proposing new cycle parameters, technicians need to repeatedly develop simulation software specifically for a particular cycle layout, which is inefficient and time-consuming. In addition, existing simulation methods lack flexibility and generality, and it is difficult to quickly adapt to diverse cycle design requirements. Therefore, developing a program that can efficiently and flexibly perform modeling and simulation is an urgent problem for technicians in this field. Summary of the Invention

[0004] To solve the above technical problems, the present invention proposes a method for quickly modeling and simulating a flow network for a supercritical carbon dioxide Brayton cycle to solve the problems existing in the above-mentioned prior art.

[0005] To achieve the above object, in the first aspect, the present invention provides a method for quickly modeling and simulating a flow network for a supercritical carbon dioxide Brayton cycle, including:

[0006] Construct an input card and read the parameters of the flow network system;

[0007] Based on the equipment node part of the input card, divide the nodes of the flow network system and generate a connection relationship matrix;

[0008] Based on the three major conservation equations and the calculation results of the previous time step, construct the pressure matrix and enthalpy matrix of the flow network;

[0009] Solve the pressure matrix and enthalpy matrix through a parallel algorithm;

[0010] Output the calculation results in real time, and realize simulations at different speeds by setting the super real-time multiple.

[0011] Preferably, the steps of constructing the input card include:

[0012] Design a specific input card format to allow users to input various parameters of the system, including node size, pipeline connection relationship, fluid physical property parameters, flow rate data, external energy and momentum input;

[0013] The input card is provided in the form of a text file, and the program automatically reads and parses the input data;

[0014] The types of devices that can be input through the input card include reactors, turbines, compressors, heat exchangers, valves, pipelines, temperature boundaries, pressure boundaries, flow rate boundaries, power boundaries, and individual nodes.

[0015] Preferably, the steps of dividing nodes for the flow network system include:

[0016] Divide the models of turbines and compressors into two nodes plus one streamline, and longitudinally node the model of the heat exchanger in a multi-segmented manner, with each longitudinal node using the method of two heat component nodes;

[0017] Assign an identification number to each node, store it in the streamline, and generate a connection relationship matrix with the identification number as the coordinate and the position of the node in the node matrix as the value for identifying the connection relationship during subsequent flow network calculations.

[0018] Preferably, the pressure drop of the turbine is calculated using the Stodola ellipse method.

[0019] Preferably, the relationship between the flow rate, rotational speed and pressure ratio of the compressor is obtained through calculation of the characteristic curve.

[0020] Preferably, the heat components of the heat exchanger are divided into nodes and calculated in the manner of two nodes, and a two-dimensional matrix of the energy equations for the two heat component nodes is constructed and solved.

[0021] Preferably, the steps of constructing the pressure matrix and enthalpy matrix of the flow network include:

[0022] Calculate the equipment parameters, linear terms and constant terms of the energy equation, and linear terms and constant terms of the mass equation;

[0023] Through the semi-implicit difference format of the momentum equation, gradually calculate the pressure change and enthalpy change of each node in the flow network at the next time step.

[0024] Preferably, the parallel algorithm uses the Gaussian elimination method with partial pivoting to solve the matrix, and the specific steps include:

[0025] Perform partial pivoting on the target matrix;

[0026] Perform LU decomposition on the target matrix. First, calculate the first row of the U matrix and the first column of the L matrix in parallel;

[0027] Calculate the rows of the U matrix and the columns of the L matrix corresponding to the diagonal elements in parallel;

[0028] Calculate the intermediate quantity of the result variable and finally calculate the result.

[0029] Preferably, the steps of real-time outputting the calculation result include:

[0030] Real-time output the calculation result in the form of a txt file, recording the pressure, temperature, flow rate and equipment parameters of the system;

[0031] Output the total duration calculation result in the form of a csv file, recording the changes in the system state during the entire simulation process.

[0032] Preferably, the steps of implementing simulations at different speeds by setting the super-real-time multiple include:

[0033] The user sets the super-real-time multiple in the input card, and the program accelerates the simulation process according to the set multiple;

[0034] In the real-time simulation mode, the calculation is synchronized with the time of the physical reality; in the super-real-time simulation mode, the system accelerates the calculation process by increasing the multiple.

[0035] Compared with the prior art, the present invention has the following advantages and technical effects:

[0036] The present invention provides a fast flow network modeling and simulation method for a supercritical carbon dioxide Brayton cycle. First, construct an input card and read the parameters of the flow network system; second, based on the device node part of the input card, divide the nodes of the flow network system and generate a connection relationship matrix; then, based on the three major conservation equations and the calculation results of the previous time step, construct the pressure matrix and enthalpy matrix of the flow network; again, solve the pressure matrix and enthalpy matrix through a parallel algorithm; finally, real-time output the calculation result and implement simulations at different speeds by setting the super-real-time multiple.

[0037] The present invention can perform various arrangements and parameter settings, and can model the carbon dioxide Brayton cycle flow network program through different device combinations, which can improve the modeling efficiency of technicians; and can perform real-time and super-real-time simulations to achieve the purposes of design and teaching. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] The drawings forming a part of this application are used to provide a further understanding of this application. The schematic embodiments of this application and their descriptions are used to explain this application and do not constitute an improper limitation to this application. In the drawings:

[0039] Figure 1 It is the program structure division diagram of the embodiment of the present invention;

[0040] Figure 2 It is the calculation flow chart of the semi-implicit difference method adopted in the embodiment of the present invention;

[0041] Figure 3 It is the schematic diagram of the turbine node model of the embodiment of the present invention;

[0042] Figure 4 It is the schematic diagram of the compressor node model of the embodiment of the present invention;

[0043] Figure 5 It is the schematic diagram of the calculation input card of the embodiment of the present invention;

[0044] Figure 6 It is the diagram of the txt output file of the embodiment of the present invention;

[0045] Figure 7 It is the diagram of the csv output file of the embodiment of the present invention. Specific implementation manners

[0046] It should be noted that, without conflict, the embodiments in this application and the features in the embodiments can be combined with each other. The following will describe this application in detail with reference to the drawings and in combination with the embodiments.

[0047] It should be noted that the steps shown in the flowchart of the drawings can be executed in a computer system such as a set of computer executable instructions, and although the logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in a different order from here.

[0048] Embodiment 1

[0049] As Figure 1 shown, in this embodiment, a fast flow network modeling and simulation method for a supercritical carbon dioxide Brayton cycle is provided, including:

[0050] S1. Construct an input card and read the parameters of the flow network system;

[0051] Further, the steps of constructing the input card include:

[0052] Design a specific input card format to allow users to input various parameters of the system, including node size, pipeline connection relationship, fluid physical property parameters, flow rate data, external energy and momentum input;

[0053] The input card is provided in the form of a text file, and the program automatically reads and parses the input data.

[0054] Specifically, analyze the supercritical carbon dioxide Brayton cycle system and write input cards. The types of devices that can be input in the input cards include reactors, turbines, compressors, heat exchangers, valves, pipelines, temperature boundaries, pressure boundaries, flow rate boundaries, power boundaries, and individual nodes. Among them, the identification number of the reactor is hprt, the turbine is turb, the compressor is comp, the heat exchanger is htch, the valve is valv, the flow rate boundary is flbd, the pressure boundary is prbd, the power boundary is pwbd, the temperature boundary is tmbd, the individual node is sgnd, and the pipeline is pipe. The input card ends with end.

[0055] As an innovative implementation method, provide characteristic curve parameter files for turbines and compressors. The writing format of the turbine characteristic curve file is flow rate (kg / s), pressure ratio, rotational speed (r / min), and efficiency (%). The writing format of the compressor characteristic curve file is flow rate (kg / s), pressure ratio, and efficiency (%).

[0056] The program startup program, the basic information reading part and the device reading part in the input card module read the parameters of the flow network system. These parameters include the dimensions of each node in the flow network system, the pipeline connection relationship, the fluid physical property parameters, the flow rate data, and the external energy and momentum input. The device parameter processing part of the program will parse and process these input data to provide a basis for subsequent calculation processes.

[0057] S2. Based on the device node division part of the input card, divide the nodes of the flow network system and generate a connection relationship matrix;

[0058] Furthermore, the steps for dividing the nodes of the flow network system include:

[0059] Divide the models of the turbine and the compressor into two nodes plus a streamline, and longitudinally node the model of the heat exchanger in a multi-segmented manner. Each longitudinal node adopts two heat component nodes (heating end and cold end flow channel nodes);

[0060] Assign identification numbers to each node, store them in the streamline, and generate a connection relationship matrix with the identification numbers as coordinates and the positions of the nodes in the node matrix as values for identifying the connection relationship during subsequent flow network calculations.

[0061] Specifically, the device node division part in the input card module automatically divides the nodes of the flow network system and generates a connection relationship matrix. It includes automatically dividing the nodes of the flow network system according to the input parameter information. The program will judge the positions of each node in the flow network and their mutual connection relationships according to the input card. According to the connection method of the nodes, the program will automatically generate a connection relationship matrix for solving the subsequent flow network system.

[0062] As an innovative implementation method, the connection relationship matrix is read by the program calculation module to identify the connection relationships of the flow network. Its specific format is as follows: The position of each node in the flow network is identified by a unique identification number, which is used to characterize the information of this node (specifically composed of the device identification number + the device quantity number + the node number in the device). The identification number of each node is used as the coordinate of the matrix and stored in the streamline. The value of the matrix under the corresponding coordinate is the position of this node in the pressure matrix and the enthalpy matrix.

[0063] If there are places with incorrect input card formats, such as duplicate device names, pipes not connected to devices, etc., the model checking part of the program will generate error messages.

[0064] The input card passes node, streamline, device, and connection relationship information to the program calculation module. This module first reads the connection relationships and generates retrieval variables for nodes and streamlines to facilitate subsequent calculations.

[0065] As an innovative implementation method, the node retrieval variable directly stores the position of the streamline in the pressure matrix, which is used to extract the inlet streamline flow rate of the node in the case of single-streamline connection inside the device; the streamline retrieval variable stores the identification number of the node, which is used to extract the node parameters at the inlet and outlet of the streamline.

[0066] The device calculation part of the program reads device parameters such as the pressure drop of the turbine and the heat exchange amount of the heat exchanger, and processes them into the heat exchange amount of the node, the head and admittance of the streamline. For turbine, compressor, heat exchanger, and valve models, the input and output of the module, as well as the processing flow, are as follows:

[0067] (1) The turbine draws on the idea of the lumped parameter method. As Figure 3 shown, it is modeled by means of two nodes and one streamline. The first node only provides input parameters. The volume of the second node is the volume of the turbine, which is used to reflect the volumetric inertia phenomenon and the output power during the operation of the turbine. The middle streamline is used to reflect the pressure drop of the turbine. In the turbine module, the input is the structural parameters of the turbine, the inlet flow rate, and the fluid physical property parameters at the inlet and outlet. The output is the pressure drop and energy loss of the turbine.

[0068] The pressure drop of the turbine is calculated using the Stodola ellipse method. The calculation formula is:

[0069]

[0070] In the formula, the subscript in is the turbine inlet, out is the turbine outlet, d is the design value, and od is the non-design value. The calculation formula of

[0071]

[0072] In the formula, F is the flow rate, kg / s; t is the node temperature, K; p is the node pressure, Pa.

[0073] The calculation formula for the energy loss of the turbine model is as follows:

[0074]

[0075] In the formula, υ in and υ out are the velocities at the inlet and outlet of the node; h out，j，s is the isentropic stagnation outlet enthalpy, obtained by querying the physical property function:

[0076] h out，j，s = f(P out，j , s in，j ) (4)

[0077] The stagnation entropy s j is obtained by querying the physical property function, and the stagnation pressure P j and are obtained from the following formula:

[0078]

[0079] In the formula, γ is the specific heat ratio of the working fluid; υ is the velocity.

[0080] (2) As an innovative implementation method, the compressor model abandons the ideal gas assumption and draws on the idea of the lumped parameter method. As Figure 4 shown, it is modeled by means of two nodes and a streamline. The first node only provides input parameters, the volume of the second node is the volume of the compressor, used to reflect the volumetric inertia phenomenon during the operation of the compressor, and the middle streamline is used to reflect the head of the compressor. In the compressor module, the inputs are the structural parameters of the compressor, the inlet flow rate, and the fluid physical property parameters at the inlet and outlet, and the output is the head of the compressor.

[0081] The relationship between the flow rate, rotational speed and pressure ratio of the compressor model is as follows:

[0082]

[0083] In the formula, a, b, c are experimental coefficients, which need to be calculated through the provided characteristic curves; n d is the design rotational speed.

[0084] The head of the compressor is:

[0085]

[0086] The relationship between the efficiency and the flow rate is:

[0087] η = a + bF + cF 2 (8)

[0088] (3) As an innovative implementation method, the heat components of the heat exchanger are divided and calculated in a two-node manner. The input of the heat exchanger module is the structural parameters of the heat exchanger and the fluid physical property parameters of the heat exchanger nodes, and the output is the heat transfer amount of the heat exchanger. The heat balance equations of the two nodes are respectively:

[0089]

[0090] In the formula, A is the heat transfer area of the node, m 2 ; λ s is the thermal conductivity of the node, W / (m·K); x is the length of the node, m; t w1 , t w2 are the temperatures of the left and right nodes respectively, K; t hfi , t cfi are the temperatures of the hot and cold end fluids respectively, K; h h , h c are the heat transfer coefficients of the hot and cold end fluids respectively, W / (m 2 ·K);

[0091] Therefore, there is a matrix:

[0092]

[0093] Solving this matrix can obtain the boundary temperature, and then calculate the heat transfer amount.

[0094] The convective heat transfer model adopts the Gnielinski heat transfer model, and the calculation formula is as follows:

[0095]

[0096] In the formula, Nu is the Nusselt number; Pr is the Prandtl number; d is the hydraulic diameter of the flow channel, m; l is the length of the pipe section, m; t f is the fluid temperature, K; t w is the wall temperature, K; f c is the Darcy-Weisbach friction coefficient.

[0097] (4) The valve model is calculated using the critical flow model. First, it is judged whether the valve is in critical flow. The critical admittance is used for calculation. The critical admittance calculation formula is:

[0098]

[0099] In the formula, V cr is the speed of sound; A is the cross-sectional area at the narrowest part of the valve; P in is the inlet pressure of the valve; P out is the outlet pressure of the valve.

[0100] When the critical admittance is less than the streamline admittance, it is proved that critical flow has occurred at the valve. The admittance of this streamline adopts the critical admittance; otherwise, the normal admittance is adopted.

[0101] S3. Based on the three major conservation equations and the calculation results of the previous time step, construct the pressure matrix and enthalpy matrix of the flow network;

[0102] The specific forms of the three major conservation equations for solving the flow network are as follows:

[0103] (1) The expression of the momentum equation is:

[0104]

[0105] In the formula, F i is the fluid flow rate of the streamline, kg / s; A is the cross-sectional area of the flow channel, m 2 ; P i is the pressure drop, Pa; is the identification factor, equal to 0 or 1; P s，i is the head term, Pa; P f，i is the frictional resistance, Pa; P s，i is the local resistance, Pa; is the gravity term.

[0106] After integration, we get:

[0107]

[0108] Among them, the frictional resistance P f，i can be expressed by the Darcy formula:

[0109]

[0110] In the formula, λ i is the frictional resistance coefficient, L i is the node length, m, ρ i is the fluid density, kg / m 3 , d i is the node diameter, m, The frictional resistance coefficient can be calculated according to the Nikuradse experiment diagram. The Zigrang-Sylvester resistance model is used in the program.

[0111] The local resistance is defined by the following formula:

[0112]

[0113] In the formula, ξ is the local resistance coefficient, which is obtained through experiments. Generally, it takes 0.5 at the inlet of the heat exchanger and 1 at the outlet.

[0114] In summary, the total resistance and the resistance coefficient are expressed by the following formula:

[0115]

[0116] Assuming that the time variation term is 0, a semi-implicit difference process is performed, and we can get:

[0117]

[0118] Where,

[0119] (2) The energy conservation equation is expressed as:

[0120]

[0121] Where M i is the node mass, kg; h i is the node average specific enthalpy, kJ / kg; h i，in and h i.out are the node inflow and outflow specific enthalpies, kJ / kg; F i，in and F i，out are the node inflow and outflow mass flow rates, kg / s; Q f is the node heat transfer, W; W shaft The work done by the node is W.

[0122] (3) The mass conservation equation is expressed as:

[0123]

[0124] Where V i is the node volume, m 3 ρ i is the node fluid density, kg / m 3 ; F i，in and F i.out are the node inflow and outflow flow rates, kg / s, respectively.

[0125] Based on the three conservation equations and the calculation results of the previous time step, the pressure matrix and enthalpy matrix are constructed.

[0126] S4. Solve the pressure matrix and enthalpy matrix using a parallel algorithm;

[0127] As an innovative implementation, a parallel algorithm is used to solve the pressure and enthalpy matrices. This includes the use of parallel computing algorithms. Through multi-threaded distributed computing, the program can handle multiple computational tasks simultaneously, significantly reducing matrix solution time to achieve the computational speed requirements of real-time and beyond real-time simulations.

[0128] Assuming there is a linear equation system Ax=b, the parallel algorithm method is as follows:

[0129] (1) Perform pivoting on the target matrix A;

[0130] (2) Perform LU decomposition on the target matrix A. First, calculate the first row of the U matrix and the first column of the L matrix in parallel. The calculation formulas are: u 1i = a 1i , i = 2, 3, …, n and

[0131] (3) Calculate the diagonal elements of the U matrix

[0132] (4) Calculate the rows of the U matrix and the columns of the L matrix corresponding to the diagonal elements in parallel. The calculation formulas are:

[0133] and

[0134] (5) Calculate the intermediate quantity of the result variable. Assume Ly = b, then the calculation formula is:

[0135]

[0136] (6) Calculate the result x. The calculation formula is:

[0137]

[0138] Transfer the calculated node, streamline, and device data to the output card module. The output card module sequentially reads the node and streamline parameters and processes them into output parameters.

[0139] S5. Output the calculation results in real time and implement simulations at different speeds by setting the super-real-time multiple.

[0140] Further, the steps of outputting the calculation results in real time include:

[0141] Output the calculation results in real time in the txt file format, recording the pressure, temperature, flow rate, and device parameters of the system;

[0142] Output the total duration calculation results in the csv file format, recording the changes in the system state during the entire simulation process.

[0143] Specifically, outputting the calculation results in real time in the txt file format includes outputting the calculation results in the txt file format after one output time step. This file will record the parameter changes of the system, including the pressure, temperature, flow rate, and device parameters of the system. Users can read and view this text file in read-only form to understand the operating state of the flow network system in real time.

[0144] The file status check section checks whether the csv file exists. If it does not exist, a new file is created and the header is written; if it exists, the output parameters are printed on a new line. The total duration calculation result is output in the csv file format, including all system change processes during the calculation, in the csv file format. This file details the start and end times of the entire simulation process and the system states corresponding to the times of each step, helping the user to plot the system parameter change graph.

[0145] Real-time and super-real-time simulations are achieved by setting the super-real-time multiple. The real-time simulation control module allows the user to set the super-real-time multiple (such as 1 times, 10 times, 100 times, etc.) in the input card to control the simulation speed. With this setting, the system can accelerate the simulation process as needed. In the real-time simulation mode, the calculation is synchronized with the time of the physical reality; while in the super-real-time simulation mode, the system accelerates the calculation process by increasing the multiple.

[0146] Combined with Figure 5 and Figure 7 Describe the specific implementation manner of the program.

[0147] Figure 5 This is the input card for the supercritical carbon dioxide simple Brayton cycle. The numbers in the first row are the simulation time, time step, output time step, and the program is default set to real-time calculation. The subsequent devices are the reactor, turbine, heat exchanger, compressor, precooler, the boundary of the cold end of the precooler, and the pipeline connecting the devices. The input card ends with "end" to indicate the end of the input.

[0148] Click to run the program, and the program starts to calculate. By reading the Figure 6 txt file, the system parameter changes are displayed in real-time. All the parameters calculated by the program are stored in the Figure 7 csv file shown, which is convenient for the user to view the system parameter change situation later.

[0149] The above is only the preferred specific implementation manner of the present application, but the protection scope of the present application is not limited thereto. Any changes or substitutions that can be easily thought of by those skilled in the art within the technical scope disclosed by the present application should be covered by the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A rapid flow network modeling and simulation method for supercritical carbon dioxide Brayton cycle, characterized in that It includes the following steps: Construct an input card and read the parameters of the flow network system; Based on the device node part of the input card, divide the nodes of the flow network system and generate a connection relationship matrix; Based on the connection relationship matrix and device parameters, calculate the linear term and constant term of the conservation equation; Based on the three major conservation equations and the calculation results of the previous time step, construct the pressure matrix and enthalpy matrix of the flow network; Solve the pressure matrix and enthalpy matrix through a parallel algorithm; Output the calculation results in real time and achieve simulations at different speeds by setting the super real-time multiple.

2. The method according to claim 1, wherein The steps of constructing the input card include: Design a specific input card format to allow users to input various parameters of the system, including node size, pipeline connection relationship, fluid physical property parameters, flow rate data, external energy and momentum input; The input card is provided in the form of a text file, and the program automatically reads and parses the input data; The types of devices that can be input on the input card include reactors, turbines, compressors, heat exchangers, valves, pipelines, temperature boundaries, pressure boundaries, flow rate boundaries, power boundaries, and single nodes.

3. The method according to claim 2, wherein The steps of dividing the nodes of the flow network system include: Divide the models of turbines and compressors into two nodes plus a streamline, and longitudinally node the model of the heat exchanger in a multi-segmented manner, with each longitudinal node adopting the method of two heat component nodes; Assign an identification number to each node, store it in the streamline, and generate a connection relationship matrix with the identification number as the coordinate and the position of the node in the node matrix as the value for identifying the connection relationship during subsequent flow network calculations.

4. The method according to claim 2, wherein The pressure drop of the turbine is calculated using the Stodola ellipse method.

5. The method according to claim 2, wherein The relationship between the flow rate, rotational speed and pressure ratio of the compressor is obtained by calculating the characteristic curve.

6. The method according to claim 2, wherein The heat components of the heat exchanger are divided and calculated in the way of two nodes, and a two-dimensional matrix of the energy equations of the two heat component nodes is constructed and solved.

7. The method according to claim 1, wherein The steps of constructing the pressure matrix and enthalpy matrix of the flow network include: Calculate the device parameters, linear term and constant term of the energy equation, linear term and constant term of the mass equation; Through the semi-implicit difference format of the momentum equation, gradually calculate the pressure change and enthalpy change of each node in the flow network at the next time step.

8. The method according to claim 1, wherein The parallel algorithm uses the Gaussian elimination method with pivoting to solve the matrix. The specific steps include: Perform pivoting on the target matrix; Perform LU decomposition on the target matrix. First, calculate the first row of the U matrix and the first column of the L matrix in parallel; Calculate the rows of the U matrix and the columns of the L matrix corresponding to the diagonal elements in parallel; Calculate the intermediate quantity of the result variable and finally calculate the result.

9. The method according to claim 1, wherein The steps of outputting the calculation results in real time include: Output the calculation results in real time in txt file format, recording the pressure, temperature, flow rate and equipment parameters of the system; Output the total duration calculation results in csv file format, recording the system state changes during the entire simulation process.

10. The method according to claim 1, wherein The simulation steps with different speeds are realized by setting the super real-time multiple, including: The user sets the super real-time multiple in the input card, and the program accelerates the simulation process according to the set multiple; In the real-time simulation mode, the calculation is synchronized with the time of the physical reality; in the super real-time simulation mode, the system accelerates the calculation process by increasing the multiple.

Images