CN114491864B - Nuclear power pipe network model preprocessing method with parameterization and reconfigurable characteristics - Google Patents
Nuclear power pipe network model preprocessing method with parameterization and reconfigurable characteristics Download PDFInfo
- Publication number
- CN114491864B CN114491864B CN202210115993.1A CN202210115993A CN114491864B CN 114491864 B CN114491864 B CN 114491864B CN 202210115993 A CN202210115993 A CN 202210115993A CN 114491864 B CN114491864 B CN 114491864B
- Authority
- CN
- China
- Prior art keywords
- boundary
- flow
- file
- simulation
- pressure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
The invention relates to a simulation preprocessing method, in particular to a preprocessing method of a nuclear power pipe network model with parameterization and reconfigurable characteristics. Adopting a volume control method to disperse a process system, and disassembling a volume attribute and a flow attribute for modeling; determining the boundary type and position according to the coupling requirement of the simulation model and the system flow; drawing a simulation graph according to the discrete process system and the determined boundary and numbering; formulating a format specification of the parameterized data input card; and generating a parameterized data card file according to the design operation parameters and the simulation diagram of the system. The preprocessing process of the invention is provided with a trial calculation process, and the final initial value file after the trial calculation is successful is the steady-state working condition point, thereby greatly reducing the debugging time of the steady-state working condition and improving the debugging efficiency. If the trial calculation fails, no additional burden is added, and the method is consistent with the treatment of the existing mode. Thus, overall, the preprocessing method of the present invention improves debugging efficiency.
Description
Technical Field
The invention relates to a simulation preprocessing method, in particular to a preprocessing method of a nuclear power pipe network model with parameterization and reconfigurable characteristics.
Background
In the real-time fluid network simulation method of the nuclear power plant, a parametric modeling method is often adopted for modeling. In the preprocessing process of parametric modeling, the variables and arrays of each definition are often required to be manually assigned, and although the method can complete the process of parametric modeling, the manual assignment process is complicated, easy to miss filling and misfilling, and difficult to effectively embody the parametric characteristics of the model. And each automatic modeling tool only adopts an automatic assignment technology, so that an assignment process and an assignment object cannot be clearly displayed and cannot be simply modified, and a model is complex to maintain and difficult to update. The connection relation of systems and equipment is required to be changed by each simulation calculation model, modeling needs to be carried out again mostly, and the simulation calculation models hardly have reconfigurable capability.
With the deepening of simulation work in the design of the nuclear power plant, the stages of design, demonstration, check and evaluation are gradually developed. At the above stage, the requirements for the simulation method are further increased, and besides the recurrence, the functions of multi-scheme comparison, design scheme optimization, pipeline setting modification and the like need to be performed. However, in the currently adopted simulation method, a large amount of initial value assignment work of variables is required when a simulation program is generated for multiple times; the automatic modeling software cannot manually and simply change parameters, so that the modeling needs to be carried out again when a plurality of parameters are changed; when the connection relation of the pipeline and the equipment is changed, the topological structure analysis is carried out again and a simulation program is generated. The current simulation program preprocessing method has insufficient parameterization and no reconfigurable characteristic, so that a simulation model cannot finish on-line quick modification, and is difficult to be superior to other simulation applications except training simulation.
Disclosure of Invention
The invention aims to provide a preprocessing method of a nuclear power pipe network model with parameterized and reconfigurable characteristics, which is a preprocessing method in the process of simulating a nuclear power plant process system. Specifically speaking, the method is a parameter preprocessing method of a two-phase pipe network model with parameterized and reconfigurable characteristics in a nuclear power plant pipe network simulation model.
The purpose of the invention is realized as follows:
(1) Adopting a volume control method to disperse a process system, and disassembling the volume attribute and the flow attribute for modeling;
(2) Determining the boundary type and position according to the coupling requirement of the simulation model and the system flow;
(3) Drawing a simulation graph according to the discrete process system and the determined boundary and numbering;
(4) Formulating a format specification of the parameterized data input card;
(5) Generating a parameterized data card file according to design operation parameters and a simulation diagram of the system;
(6) After the parameterized data card is generated, a reading program reads the data card line by line and generates a corresponding calculation file;
(7) And after the data card is read, performing preliminary trial calculation to generate a final initial value file.
The present invention may further comprise:
1. adopting a volume control method to disperse a process system, and disassembling the volume attribute and the flow attribute for modeling:
a. the method comprises the following steps of splitting the volume attribute and the flow attribute of a pipeline, abstracting volume parameters such as pressure, mass, temperature and enthalpy into nodes, abstracting flow parameters such as flow rate and resistance into flow lines, and simulating a section of pipeline in a mode of combining the nodes and the flow lines;
b. and (4) performing process system dispersion according to measuring points, equipment, pipelines flowing by gravity, important concerned positions and simulation requirements. In a practical process system pipeline, the pressure distribution is continuous. The method is adopted to simulate a pipe network system, continuously-changed pressure needs to be dispersed into a plurality of pressure points, each pressure point is provided with a node, and the nodes are connected by utilizing a streamline.
2. Determining the boundary type and position according to the coupling requirement of the simulation model and the system flow:
a. pressure boundary: the fixed pressure participates in the operation process of the pipe network simulation, pressure boundary pressure, enthalpy and liquid level need to be given, and the pressure boundary and the flow boundary are mutually coupled, can be connected and can also exist independently. The pressure boundary is arranged at the constant pressure inlet/outlet position of the pipe network;
b. and (3) flow boundary: taking fixed mass flow to participate in the operation process of pipe network simulation, and giving flow boundary flow and enthalpy value; the flow and pressure boundaries are coupled to one another, may be connected, or may exist separately. The flow boundary is arranged at the inlet/outlet position of the fixed flow of the pipe network, and the flow boundary is additionally arranged at the inlet position of the fixed pressure for testing and debugging;
c. thermal boundary: the fixed heat exchange power participates in the operation process of pipe network simulation, heat boundary heat exchange power needs to be given, and the heat boundary and the temperature boundary are coupled with each other, can be connected and can also exist independently. The thermal boundary is provided at the electric heater/cooler, simplified heater/cooler location;
d. temperature boundary: the fixed temperature, the flowing Reynolds number and the working medium density are used for participating in the operation process of the pipe network simulation, the temperature boundary temperature, the flowing Reynolds number and the working medium density need to be given, and the temperature boundary and the heat boundary are mutually coupled and can be connected or exist independently. The temperature boundary is arranged at the position of the constant-temperature heating or radiating equipment;
e. the boundary described by the method can be calculated by giving fixed parameters; the results obtained by other software or programs through online and offline calculation can also be read as the input of the method; while being able to transfer data to other software or programs.
3. Drawing a simulation graph according to the discrete process system and the determined boundary and numbering:
a. drawing a simulation sketch according to a system flow chart, wherein pipelines, equipment and measuring points are marked in the sketch;
b. abstracting pipelines, equipment and measuring points appointed in the sketch into nodes according to the discrete nodes in the step 2, replacing the equipment in the original sketch, and adding the nodes into the simulated sketch;
c. placing a valve and a centrifugal pump on a streamline to simulate;
d. connecting the nodes with each other in the sketch according to the discrete streamline in the step 2;
e. reasonably determining the boundary according to the boundary determining method in the step 3, and connecting the boundary with the nodes and the streamline in the sketch to form a complete simulation graph;
f. and numbering nodes, streamlines and boundaries in the simulation graph. The numbering rules are as follows: the pipeline nodes, the container nodes, the pressure boundary nodes, the flow boundary nodes and the steam turbine share a set of serial numbers, wherein the serial numbers of the pipeline nodes are from 1 to 220; the container nodes are numbered from 221 to 270; pressure boundaries numbered 271-340; the flow boundaries are numbered 341-370 and the turbines numbered 371-390. Streamlines are numbered individually, from 1 to 450. The thermal boundary and the temperature boundary have no numbering rules.
4. Formulating a format specification of the parameterized data input card:
a. the data card file is a comma separated file in a csv format, can be generated through automatic software, edited by using a text editing tool, and can be intuitively edited through Excel software without the limitation of an operating system and a software environment;
b. all necessary parameters required by the operation of the model are integrated in the data card file, and direct information of the parameters is not saved except the file;
c. the data card is divided into 12 parts according to the sequence, and the 12 parts are an overview information part, a pipeline node parameter part, a container node parameter part, a streamline parameter part, a pressure boundary parameter part, a flow boundary parameter part, a heat boundary parameter part, a temperature boundary parameter part, a valve parameter part, a centrifugal pump parameter part, a steam turbine parameter part and a reserved part in sequence;
d. before each part starts, a line start identification field, such as "//1-begin//", is set to indicate that part 1 starts; after each section is finished, a line termination identification field, such as "//6-end//", is set to indicate that section 6 is finished.
5. Generating a parameterized data card file according to design operating parameters and a simulation diagram of the system:
a. generating a parameterized data card file meeting the specification of the step 5 according to the simulation diagram established in the step 4 and the element number;
b. summary information section the contents of each row are as follows: system name, system profile, gas phase composition, liquid phase composition, ambient temperature;
c. the contents of each row of the pipeline node parameter part are as follows in sequence: node number, node up/down elevation, node volume, hydraulic diameter, flow cross section area, heat dissipation surface area and wall thickness;
d. the contents of each row of the container node parameter part are as follows in sequence: node number, node elevation/volume parameter group number, elevation array, volume array, heat dissipation surface area and wall thickness;
c. the streamline parameter part comprises the following contents in sequence: the number of the streamline, the upstream/downstream nodes, the up/down elevation of the streamline, the setting resistance admittance coefficient of the streamline (the current line fills in the setting working medium density, the setting flow and the setting pressure difference in turn, if a valve is arranged on the streamline, the setting resistance admittance coefficient parameter of the streamline can not be filled), and whether the streamline has the non-return function or not;
d. the pressure boundary parameter part comprises the following contents in turn: pressure boundary number, constant pressure, working medium type (gas, liquid or input liquid level), gas phase component percentage, liquid phase component percentage, two-phase working medium specific enthalpy;
f. the flow boundary parameter part comprises the following contents in sequence: the serial number of the flow boundary, the constant two-phase flow, the percentage of the gas phase component, the percentage of the liquid phase component and the specific enthalpy of the two-phase working medium;
g. the thermal boundary parameter part comprises the following contents in sequence: heat boundary name, action node, gas phase heat exchange power and liquid phase heat exchange power;
h. the temperature boundary parameter part comprises the following contents in sequence: temperature boundary name, action node, heat exchange area, boundary temperature, boundary density and boundary Reynolds number;
i. the contents of each row of the valve parameter part are as follows in sequence: the method comprises the following steps of (1) obtaining a valve name, an action flow line, an initial opening degree, a full opening KV, different opening degree flow characteristic groups, an opening degree array and a relative KV array;
j. the contents of each row of the parameter part of the centrifugal pump are as follows in sequence: the method comprises the following steps of (1) pump name, action streamline, initial rotating speed, rated rotating speed, flow array at rated rotating speed, head array at rated rotating speed and rated characteristic working medium density;
k. the contents of each row of the turbine parameter part are as follows in sequence: the steam turbine serial number, the parameter array number, the pressure array behind the valve, the exhaust pressure array, the steam flow array and the output power array.
The reservation parameter part content is empty.
6. After the parameterized data card is generated, a reading program reads the data card line by line and generates a corresponding calculation file:
a. the card reading program first reads line 1, if it is not part 1 start field "//1-begin//" indicating that the data card is incomplete, stops running and reports an error;
b. except the 1 st part and the 12 th part, the card reading program circularly reads parameters of all parts. When the number of a part is found to be repeated or not accord with the number rule, the program is terminated and an error is reported;
c. reading parameters line by line according to a convention by a card reading program, and terminating the program and reporting an error when the line number does not accord with the convention;
d. if a certain part of the start fields are read, the end fields cannot be read, the data card is not complete, the operation is stopped, and an error is reported;
e. the card reading program has an error correction feature. If the volume is filled with a negative value, when the height of the streamline is not overlapped with the height of the upstream and downstream nodes, when the parameters of the centrifugal pump cannot be matched, when the effective efficiency of the steam turbine exceeds 95 percent and the like, the card reading is stopped, corresponding errors are reported, wrong parameters and reasons of the errors are indicated, and the nodes, the streamline, equipment numbers and row numbers in the data card are indicated;
f. after all data are read and no error exists, a dotting file, a reference module file, a preliminary initial value file and a link file in a specified format are generated at a specified user position, corresponding calculation dynamic and static link libraries are released, and success in generation is prompted.
7. After the data card is read, performing preliminary trial calculation to generate a final initial value file:
a. the initial value file is an assignment file, is not locked and is only used when the trial calculation program assigns for the first time, and each internal parameter value is only applied once and cannot be used as the simulation program assignment;
b. after the card reading program finishes reading the parameters and releases all files, trial calculation is carried out by utilizing the initial value file, the corresponding calculation dynamic state and the static link library. If trial calculation fails, prompting the failure reason, and converting the initial value file into a final initial value file; if the trial calculation is successful, stopping the calculation when the pressure difference of the same node is less than 20kPa in two continuous calculations, and taking the currently output result as a final initial value file; and if the pressure convergence condition can not be reached after 1 ten thousand times of continuous calculation, converting the initial value file into a final initial value file.
The invention provides a pipe network simulation model data preprocessing method with parameterization and reconfigurable characteristics, which is applied to a nuclear power plant process system, wherein a control volume method is used for discretizing an actual process system, determining boundary types according to the coupling and the flow of the system, drawing simulation graphs and numbering the simulation graphs, and completing physical model modeling; establishing the specification of a parameterized data input card, and determining the format and parameter filling information of each section in the input card; generating a parameterized data input card according to the drawn simulation diagram, and completing the release of a calculation file through a card reading program; and finally, generating a final initial value program through trial calculation.
The simulation method provided by the invention comprises a plurality of important steps of dispersing an actual system, determining a boundary, drawing and numbering a simulation graph, establishing a data card specification, generating a data card, reading the data card, releasing a calculation file, performing trial calculation, generating a final initial value file and the like.
The invention can carry out simulation pretreatment different from the prior mode on the real-time two-phase pipe network of the nuclear power device, and has the main advantages that:
(1) The preprocessing method has the characteristics of parameterization and reconfigurability, all parameters required by the model are concentrated in the input card by formulating the format specification of the data input card, and the modification of a small number of parameters can be directly and manually carried out at the formulated position of the input card without re-modeling; the connection relation of system equipment is modified, only streamline parameter fields in the input card need to be changed, the simulation requirements of 'a few times' and 'fine tuning comparison' in the design and evaluation stages are met, the simulation requirements of the existing functions are different, and the calculation efficiency is improved.
(2) The preprocessing process of the invention is provided with a trial calculation process, and a final initial value file after the trial calculation is successful is a steady-state working condition point, thereby greatly reducing the debugging time of the steady-state working condition and improving the debugging efficiency. If the trial calculation fails, no additional burden is added, and the method is consistent with the treatment of the existing mode. Thus, overall, the preprocessing method of the present invention improves debugging efficiency.
Drawings
FIG. 1 is a flow chart of a simulation preprocessing method of the present invention.
FIG. 2 is a flow diagram of an actual process system.
FIG. 3 is a simulated sketch of an actual process system.
FIG. 4 is a simulation diagram of a discrete actual process system.
FIG. 5 is a screenshot of a parameterized simulation data card for an actual process system.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
With reference to fig. 1, the implementation steps of the present invention mainly include:
(1) And drawing a simulation sketch link. The link is that according to the characteristics of the simulation object system, the actual process system is converted into a simulation sketch with labels, and according to the actual structure of the simulation object, the system boundary is selected, and the simulation sketch is perfected. The specific implementation method comprises the following steps:
a. determining a measuring point, equipment, a pipeline flowing by gravity, a valve position with a non-return function paying attention to in a focused manner;
b. drawing a simulation sketch on the basis of a system flow chart, and marking the determined position;
c. determining the nature of the boundary according to the setting of an actual process system and the type of a corresponding interface of a system to be connected with the actual process system, and marking the boundary in a simulated sketch;
d. and erasing other information and determining a simulation sketch.
(2) And (5) establishing a simulation diagram. The link is to split the system into a simulation graph with nodes, streamlines and boundaries interconnected according to the labeling of the simulation sketch, and the specific implementation method is as follows:
a. abstracting a focus point, a pressure measuring point, a temperature measuring point, a gravity flow point and the like into nodes, a steam turbine, a pump and a valve at corresponding positions according to the mode of the step 1 of the invention content, and marking the nodes, the steam turbine, the pump and the valve in a figure;
b. according to the step c of drawing the simulation sketch link, placing a proper boundary at a corresponding position in the drawing;
c. according to the process flow chart, connecting the boundary, the node and the steam turbine by using a streamline, and determining that the streamline exists at the flow measuring point, otherwise, supplementing the streamline;
d. placing equipment such as a valve, a pump and the like on a flow line at a corresponding position;
e. numbering nodes, streamlines, turbines and boundaries in sequence according to the mode of step 3 f of the invention;
f. and thermal parameter information is marked beside the streamline, the node and the boundary, so that debugging is facilitated, and the establishment of the simulation diagram is completed.
(3) And (5) establishing a parameterized data card. The link completes the establishment of the parameterized data card by the design and operation parameters of an actual system based on the topological relation provided by a simulation diagram according to the format specification of the parameterized data card, and the specific implementation method comprises the following steps:
a. establishing a csv-format parameterized data card according to the format specification of the parameterized data card, setting the data input card into 12 sections, wherein each section is provided with a start mark field// n-begin// (n is a field number, the same below) and a termination field// n-end/;
b. determining a section to be filled according to the topological relation provided in the simulation graph;
c. filling in summary information fields according to simulation requirements;
d. filling in pipeline node parameter fields, container parameter node fields, boundary parameter fields and steam turbine parameter fields according to pipelines pointed by discrete nodes in the simulation graph and based on design structure parameters;
e. filling in a streamline parameter field, a valve parameter field and a centrifugal pump parameter field according to the topological connection relation provided by the simulation diagram and the pipeline resistance indicated by the streamline;
f. and (4) the reserved field is not filled in, and the establishment of the parameterized data card is completed.
(4) And (5) a card reading link. The link is to read the parameterized data card established in the step 3 by using an automatic program, and release a calculation file and a preliminary initial value file according to the content of the data card to finish a basic preprocessing link. The specific implementation method comprises the following steps:
a. the automatic card reading program automatically matches the filled data card file positions according to the configuration information;
b. reading in the parameterized data card line by line, and temporarily storing the related parameters in a memory;
c. generating variable names required by all computations meeting the requirements according to the requirements of the real-time simulation platform and the program computation, and writing the variable names into corresponding dotting files;
d. according to the scale of the node, the matrix and the boundary, selecting and releasing proper scale to calculate and solve the dynamic and static link libraries;
e. and c, generating a Fortran code of the interface program according to the variable name generated in the step c, and completing the generation of the subprogram for calling the link library.
f. And outputting each array variable in the memory as a preliminary initial value file. After the card reading program is completed, the memory is not released temporarily.
(5) And (5) trial calculation. The link calls a calculation program in a card reading program by using the read parameters, performs preliminary calculation on the parameters and generates final initial value parameters. The specific implementation method comprises the following steps:
a. the card reading program firstly calls the solving link library released in the step (4) by utilizing data written into the memory in the early stage;
b. calling a link library and then performing operation solving;
c. if trial calculation fails, directly taking the initial value file as a final initial value file for output, finishing preprocessing, terminating the card reading program and releasing the related memory; if the trial calculation is successful, continuously running and recording the calculation steps;
d. when the pressure deviation calculated twice by the same node in all the nodes at certain two adjacent times is less than 20kPa, the calculation is considered to be converged, the current calculation result is output as a final initial value file, the preprocessing is finished, the card reading program is terminated, and the related memory is released;
e. and when the calculation is carried out for 1 ten thousand times and the convergence state is not reached, the calculation is considered to be not converged, the initial value file is directly taken as the final initial value file to be output, the preprocessing is finished, the card reading program is terminated, and the related memory is released.
The application method of the pretreatment method provided by the invention in an actual process system is explained by combining with the figures 2-5. Wherein figure 2 is a flow diagram of an actual process system with piping, instrumentation, valves, and related parameters indicated above. In FIG. 2, F is the flow measurement point, P is the pressure measurement point, and dP is the differential pressure measurement point; the letter next to the valve is the name of the valve.
(1) And drawing a simulation sketch.
a. Determining 4 pressure measuring points, 1 differential pressure measuring point, 3 flow measuring points, 4 pipeline converging/branching points, 1 common valve and 2 check valves in the graph;
b. marking the position information on a system flow chart only keeping the direction of the pipeline, and preliminarily drawing a simulation sketch;
c. according to the actual inlet and outlet states of the process system in FIG. 2, determining that a flow boundary is used as an inlet and a pressure boundary is used as an outlet, and marking the flow boundary and the pressure boundary in a simulated grass diagram;
d. and erasing other information and determining a simulation sketch. The simulation sketch determined according to the system flow diagram 2 is shown in fig. 3.
(2) And establishing a simulation graph.
a. Setting pressure measuring points and converging/branching points as nodes separately; dividing two sides of the valve with the pressure difference measuring points into two nodes, and placing the valve between the two nodes at the response position;
b. 1 flow boundary is arranged at the inlet of the system, and 1 pressure boundary is respectively arranged at 3 outlets of the system;
c. according to the process flow chart, the boundary and the node are connected by using the flow lines, so that the flow lines corresponding to the 4 flow measurement points can be determined, and the flow lines do not need to be supplemented;
d. c, placing three valve devices on the streamline near the position placed in the step a to finish the primary drawing of the simulation diagram;
e. numbering nodes, streamlines and boundaries in sequence, wherein all the nodes are pipeline nodes, and are numbered 1-6; the streamline is numbered from 1 to 10; pressure boundary numbers 271-273; a traffic boundary number 341 filling the above number in the simulation graph;
f. and thermal parameter information is marked beside the streamline, the node and the boundary, so that debugging is facilitated, and the establishment of the simulation diagram is completed. The simulation graph created from the simulation sketch 3 is shown in fig. 4.
(3) And establishing a parameterized data card. The method is realized by an automatic card filling program, and manual errors can be avoided to the maximum extent. However, in order to better illustrate the content of the present invention, the present bar is described in a manual card filling mode.
a. Establishing a TET.csv file under a D: \ Simtools \ Users \ TET \ card path, or copying a parameterized data card template and then changing names, and filling a system name (TET), a system brief description, gas phase components (water vapor, nitrogen and oxygen), liquid phase components (water and dissolved salt) and an environment temperature (18.5) according to prompts according to the running state of an actual system after summarizing an initial identification of an information field;
b. after the initial identification of the pipeline node parameter part, according to the simulated simplified combined pipeline, filling parameters such as the number, the upper and lower elevations, the volume, the hydraulic diameter, the flow section, the heat dissipation surface area, the wall thickness and the like of each node in sequence, and finishing after filling from the node 1 to the node 6;
c. the embodiment does not contain container nodes, and partial parameters of the container nodes are not filled;
d. after the initial mark of the streamline parameter part, sequentially filling a streamline number, an upstream node, a downstream node, an upper elevation, a lower elevation and a resistance admittance coefficient set by the streamline (the current line sequentially fills set working medium density, set flow and set pressure difference, and the current lines 3, 5 and 9 do not fill the current line), and whether the streamline has a non-return effect (filling 0 or 1, in the example, filling 1 on the current lines 3 and 9, and filling 0);
e. after the initial mark of the pressure boundary parameter part, filling a pressure boundary number, constant pressure, a working medium type (in this example, all liquid), gas phase component percentage (in this example, filling 3 numbers, namely steam, nitrogen and oxygen, corresponding to the step a), liquid phase component percentage (in this example, filling 1 number, namely dissolved salt, corresponding to the step a) and two-phase working medium specific enthalpy in sequence;
f. after the initial mark of the flow boundary parameter part, filling the flow boundary number, the constant two-phase flow (gas phase 0 and liquid phase 36), the gas phase component percentage (in this example, 3 numbers are filled, namely, steam, nitrogen and oxygen, corresponding to the step a), the liquid phase component percentage (in this example, 1 number is filled, namely, dissolved salt, corresponding to the step a), and the specific enthalpy of the two-phase working medium;
g. the example does not include a thermal boundary and a temperature boundary, and the parameters of the two parts are not filled;
h. after the initial mark of the valve parameter part, filling the name of the valve, the streamline where the valve is located, the initial opening, the full-open KV, the flow characteristic group number with different opening degrees (for the valves on the streamline 3 and the streamline 9 without adjusting the opening degree, filling '2' representing the two states of full-open and full-close; the adjusting valve data on the streamline 5 is filled as 11 groups), the opening degree array and the relative KV array in sequence;
i. the present example does not include a centrifugal pump and a steam turbine, and the two parts and the reserved part are not filled in.
j. After completion of filling, whether each field is complete or not is checked in detail, and whether the start/end fields are corresponding or not is checked. After confirmation, the filling of the parameterized data card is completed, and the data card style of this example is shown in fig. 5.
(4) And reading the card and releasing the calculation file.
a. The automatic card reading program automatically matches the position where the filled data card file is stored according to the configuration information, wherein the example is D \ Simtools \ Users \ TET;
b. reading the parameterized data card line by line, and temporarily storing the related parameters in a memory;
c. generating variable names required by all calculations which all meet the requirements according to the requirements of a real-time simulation platform and program calculation, writing the variable names into corresponding dotting files, and creating a variable dotting file TET.addv and a constant dotting file TET.addc in a folder D of Simtools \ Users \ TET \ text, wherein the dotting files comprise variable names, variable descriptions, variable types, variable display formats, initial values of the system and single variable points and root nodes;
d. according to the scales of the nodes, the matrix and the boundary of the embodiment, a matrix solver dynamic and static link library of 10 is released under a D: \ Simtools \ Users \ TET path, and a modeling block file incl.TET is created in a D: \ Simtools \ Users \ TET \ text folder;
e. and c, generating an interface program Fortran code file TETmin.for in a folder D: \ Simtools \ Users \ TET \ src according to the variable name generated in the step c, and completing the generation of the subprogram for calling the link library.
f. And outputting the array variables D \ Simtools \ Users \ TET \ temp in the memory as a preliminary initial value file TET. After the card reading program is completed, the memory is not released temporarily.
(5) And (5) trial-calculating and releasing the final initial value file.
a. The card reading program firstly calls the 10 × 10 solving link library by using data written into the memory at the early stage;
b. calling a link library and then carrying out operation solution, wherein trial calculation is successful, and convergence conditions are achieved in 1833 steps;
c. and outputting the current calculation result as a final initial value file, finishing preprocessing, terminating the card reading program and releasing the related memory.
After the preprocessing is finished, ensuring that 1 solving static link library, 1 solving dynamic link library and 1 solving dll file are contained in the corresponding D < Lamtobols \ Users \ TET folder; an initial value folder bd containing a final initial value file bd.TET; a code program folder src which contains an interface program TETmin. The document folder text comprises a module file incl.TET, a variable plus point file TET.addv and a constant plus point file TET.addc; a temporary folder temp containing a preliminary initial value file TET. The data card folder card contains a parameterized data card tet.
The invention provides a preprocessing method of a nuclear power pipe network model with parameterized and reconfigurable characteristics. (1) Adopting a volume control method to disperse a process system, and disassembling the volume attribute and the flow attribute for modeling; (2) Determining the type and position of the boundary according to the coupling requirement of the simulation model and the system flow; (3) Drawing a simulation graph according to the discrete process system and the determined boundary and numbering; (4) formulating a format specification of the parameterized data input card; (5) Generating a parameterized data card file according to design operation parameters and a simulation diagram of the system; (6) After the parameterized data card is generated, a reading program reads the data card line by line and generates a corresponding calculation file; (7) And after the data card is read, performing preliminary trial calculation to generate a final initial value file. The simulation method provided by the invention can complete the preprocessing process of converting the process system flow chart and the design parameters into the pipe network model calculation file, and the generated parameterized data card covers all parameters required by model calculation and has the characteristics of parameterization and modularization; the calculation file has the capability of quickly changing the connection relation among pipelines, equipment and systems, and has reconfigurable characteristics. The parameterized data card provided by the invention is standard and simple, has uniform format, has the advantage of no influence of an operating system and an operating environment, and can be manually filled in or automatically generated by software.
Claims (4)
1. A nuclear power pipe network model preprocessing method with parameterization and reconfigurable characteristics is characterized by comprising the following steps:
(1) Adopting a volume control method to disperse a process system, and disassembling the volume attribute and the flow attribute for modeling;
the method comprises the following steps of splitting the volume attribute and the flow attribute of a pipeline, abstracting the volume parameters of pressure, mass, temperature and enthalpy value into nodes, abstracting the flow parameters of flow and resistance into flow lines, and simulating a section of pipeline in a mode of combining the nodes and the flow lines;
dispersing a process system according to measuring points, equipment, pipelines flowing by gravity, important concerned positions and simulation requirements; dispersing the continuously changed pressure into a plurality of pressure points, arranging a node at each pressure point, and connecting the nodes by using a streamline;
(2) Determining the boundary type and position according to the coupling requirement of the simulation model and the system flow;
the boundary types comprise a pressure boundary, a flow boundary, a heat boundary and a temperature boundary;
pressure boundary: the fixed pressure participates in the operation process of the pipe network simulation, pressure boundary pressure, enthalpy and liquid level need to be given, and the pressure boundary and the flow boundary are mutually coupled, can be connected and can also exist independently; the pressure boundary is arranged at the constant pressure inlet/outlet position of the pipe network;
and (3) flow boundary: taking fixed mass flow to participate in the operation process of pipe network simulation, and giving flow boundary flow and enthalpy value; the flow boundary and the pressure boundary are coupled with each other, can be connected and can also exist independently; the flow boundary is arranged at the inlet/outlet position of the fixed flow of the pipe network, and the flow boundary is additionally arranged at the inlet position of the fixed pressure for testing and debugging;
thermal boundary: the fixed heat exchange power participates in the operation process of the pipe network simulation, heat boundary heat exchange power needs to be given, and the heat boundary and the temperature boundary are coupled with each other, can be connected and can also exist independently; the thermal boundary is provided at the electric heater/cooler, simplified heater/cooler location;
temperature boundary: the method includes the steps that a fixed temperature, a flow Reynolds number and working medium density are used for participating in the operation process of pipe network simulation, temperature boundary temperature, a flow Reynolds number and the working medium density need to be given, and the temperature boundary and a heat boundary are coupled with each other and can be connected or exist independently; the temperature boundary is arranged at the position of the constant-temperature heating or radiating equipment;
(3) Drawing a simulation graph according to the discrete process system and the determined boundary and numbering;
(3.1) drawing a simulation sketch according to the system flow chart, wherein pipelines, equipment and measuring points are marked in the sketch;
(3.2) abstracting pipelines, equipment and measuring points appointed in the sketch into nodes according to the discrete nodes, replacing the equipment in the original sketch, and adding the nodes into the simulated sketch;
(3.3) placing the valve and the centrifugal pump on a streamline to carry out simulation;
(3.4) connecting the nodes to each other in the sketch according to the discrete streamline;
(3.5) reasonably determining the boundary according to the boundary determining method, and connecting the boundary with the nodes and the streamline in the sketch to form a complete simulation graph;
(3.6) numbering nodes, streamlines and boundaries in the simulation graph;
the numbering rules are as follows: the pipeline nodes, the container nodes, the pressure boundary nodes, the flow boundary nodes and the steam turbine share one set of serial numbers, the streamlines are individually numbered, and the heat boundary and the temperature boundary have no serial number rule;
(4) All parameters are integrated in a parameterized data card file, the file is divided into 12 parts, and clear demarcation field lines are arranged among the parts; can be generated by an automatic program or manually filled and modified;
(5) Generating a parameterized data card file according to design operation parameters and a simulation diagram of the system; the parameterized data card has parameterized features;
(6) After the parameterized data card is generated, reading the data card line by a reading program and generating a corresponding calculation file, wherein the calculation file has reconfigurable characteristics;
(7) After the data card is read, performing preliminary trial calculation to generate a final initial value file, and finishing preprocessing;
the initial value file is an assignment file, is not locked and is only used when the trial calculation program assigns for the first time, and each internal parameter value is only applied once and cannot be used as the simulation program assignment;
after the card reading program finishes reading the parameters and releases all files, trial calculation is carried out by utilizing the initial value file, the corresponding calculation dynamic state and the static link library; if trial calculation fails, prompting the failure reason, and converting the initial value file into a final initial value file; if the trial calculation is successful, stopping the calculation when the pressure difference of the same node is less than 20kPa in two continuous calculations, and taking the currently output result as a final initial value file; if the pressure convergence condition can not be achieved after 1 ten thousand times of continuous calculation, converting the initial value file into a final initial value file; after completion, the preprocessing process ends.
2. The method for preprocessing the nuclear power pipe network model with the parameterized and reconfigurable characteristics as claimed in claim 1, wherein the method comprises the following steps: the step (4) specifically comprises:
the data card file is a comma separated file in a csv format, can be generated by automatic software, edited by a text editing tool, and can be intuitively edited by Excel software without the limitation of an operating system and a software environment;
all necessary parameters required by the operation of the model are integrated in the data card file, and direct information of the parameters is not saved except the file;
the data card is divided into 12 parts according to the sequence, and the 12 parts are an overview information part, a pipeline node parameter part, a container node parameter part, a streamline parameter part, a pressure boundary parameter part, a flow boundary parameter part, a heat boundary parameter part, a temperature boundary parameter part, a valve parameter part, a centrifugal pump parameter part, a steam turbine parameter part and a reserved part in sequence;
before each part starts, a line start identification field "//1-begin//" is set to indicate that part 1 starts; after each section is finished, a row of end identification fields "//6-end//" is set to indicate that section 6 is finished.
3. The method for preprocessing the nuclear power pipe network model with the parameterized and reconfigurable characteristics as claimed in claim 1, wherein the method comprises the following steps: the step (5) specifically comprises:
generating an established parameterized data card file according to the established simulation diagram and the element number; the information stored in the parameterized data card file is modular information and has parameterized characteristics;
summary information section the contents of each row are as follows: system name, system profile, gas phase composition, liquid phase composition, ambient temperature;
the contents of each row of the pipeline node parameter part are as follows in sequence: node number, node up/down elevation, node volume, hydraulic diameter, flow cross section area, heat dissipation surface area and wall thickness;
the contents of each row of the container node parameter part are as follows in sequence: node number, node elevation/volume parameter group number, elevation array, volume array, heat dissipation surface area and wall thickness;
the streamline parameter part comprises the following contents in sequence: the streamline number, the upstream/downstream nodes, the streamline upper/lower elevation, the streamline set resistance admittance coefficient and whether the streamline has the non-return function;
the pressure boundary parameter part comprises the following contents in turn: pressure boundary number, constant pressure, working medium type, gas phase component percentage, liquid phase component percentage and two-phase working medium specific enthalpy;
the flow boundary parameter part comprises the following contents in sequence: the serial number of the flow boundary, the constant two-phase flow, the percentage of the gas phase component, the percentage of the liquid phase component and the specific enthalpy of the two-phase working medium;
the contents of each row of the thermal boundary parameter part are as follows in sequence: heat boundary name, action node, gas phase heat exchange power and liquid phase heat exchange power;
the temperature boundary parameter part comprises the following contents in sequence: temperature boundary name, action node, heat exchange area, boundary temperature, boundary density and boundary Reynolds number;
the contents of each row of the valve parameter part are as follows in sequence: the method comprises the following steps of (1) obtaining a valve name, an action flow line, an initial opening degree, a full opening KV, different opening degree flow characteristic groups, an opening degree array and a relative KV array;
the parameters of the centrifugal pump are as follows: the method comprises the following steps of (1) pump name, action flow line, initial rotating speed, rated rotating speed, flow array at rated rotating speed, head array at rated rotating speed and rated characteristic working medium density;
the contents of each row of the turbine parameter part are as follows in sequence: the steam turbine number, the parameter group number, the pressure array behind the valve, the exhaust pressure array, the steam flow array and the output power array;
the content of the reserved parameter part is empty.
4. The method for preprocessing the nuclear power pipe network model with the parameterized and reconfigurable characteristics as claimed in claim 1, wherein the method comprises the following steps: the step (6) specifically comprises:
the card reading program first reads line 1, if it is not part 1 start field "//1-begin//" indicating that the data card is incomplete, stops running and reports an error;
except the 1 st part and the 12 th part, the card reading program circularly reads parameters of all parts; when the number of a part is found to be repeated or not accord with the number rule, the program is terminated and an error is reported;
reading parameters line by line according to a convention by a card reading program, and terminating the program and reporting an error when the line number does not accord with the convention;
after reading a certain part of the start field, failing to read the end field, indicating that the data card is incomplete, stopping running and reporting errors;
the card reading program has an error correction feature; the method comprises the following steps of filling a negative value in the volume, terminating card reading and reporting a corresponding error under the conditions that the height of a streamline is not overlapped with the heights of upstream and downstream nodes, the parameters of a centrifugal pump cannot be matched and the effective efficiency of a steam turbine exceeds 95%, indicating the wrong parameters and the reason of the error, and indicating the nodes, the streamline, the equipment number and the row number in a data card;
after all data are read and no error exists, generating a dotting file, a quoting module file, a preliminary initial value file and a link file in a specified format at a specified user position, releasing corresponding calculation dynamic and static link libraries, and prompting that the generation is successful; the connection relations in the files are symbolic functions and can be changed at any time, so that the reconfigurable character is achieved.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210115993.1A CN114491864B (en) | 2022-01-26 | 2022-01-26 | Nuclear power pipe network model preprocessing method with parameterization and reconfigurable characteristics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210115993.1A CN114491864B (en) | 2022-01-26 | 2022-01-26 | Nuclear power pipe network model preprocessing method with parameterization and reconfigurable characteristics |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114491864A CN114491864A (en) | 2022-05-13 |
| CN114491864B true CN114491864B (en) | 2022-12-13 |
Family
ID=81479242
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210115993.1A Active CN114491864B (en) | 2022-01-26 | 2022-01-26 | Nuclear power pipe network model preprocessing method with parameterization and reconfigurable characteristics |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114491864B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115079592B (en) * | 2022-07-12 | 2023-07-18 | 哈尔滨工程大学 | Pipe network simulation method for thermodynamic system of ship nuclear power device |
| CN115600322B (en) * | 2022-11-11 | 2023-07-14 | 哈尔滨工程大学 | Simulation method and system for turbine pump of ship nuclear power plant system |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2523640B (en) * | 2012-12-20 | 2020-05-27 | Inst Of Modern Physics | Particle flow simulation system and method |
| CN107315858B (en) * | 2017-06-01 | 2020-08-11 | 哈尔滨工程大学 | Centrifugal pump parametric simulation method for flow network simulation application |
| US11687762B2 (en) * | 2018-02-27 | 2023-06-27 | Stmicroelectronics S.R.L. | Acceleration unit for a deep learning engine |
| CN108595752B (en) * | 2018-03-19 | 2021-12-07 | 哈尔滨工程大学 | Single-phase water fluid network simulation method oriented to nuclear power plant simulation application |
| CN108875151A (en) * | 2018-05-07 | 2018-11-23 | 哈尔滨工程大学 | Nuclear power secondary coolant circuit system parameterized simulation model rapid build system and construction method |
| CN109255164B (en) * | 2018-08-23 | 2019-10-29 | 中国水利水电科学研究院 | A kind of earth's surface two dimension underground pipe network hydrodynamics coupling process based on space topological |
| CN109255185A (en) * | 2018-09-11 | 2019-01-22 | 中国水利水电科学研究院 | An a kind of two-dimentional hydrodynamics coupling analytical method based on urban surface underground pipe network |
| CN110705078B (en) * | 2019-09-25 | 2023-04-21 | 哈尔滨工程大学 | Simulation model development system and analysis method of nuclear power plant |
| CN112752347B (en) * | 2019-10-31 | 2023-04-07 | 华为技术有限公司 | Method and apparatus for scheduling data transmission |
-
2022
- 2022-01-26 CN CN202210115993.1A patent/CN114491864B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN114491864A (en) | 2022-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114491864B (en) | Nuclear power pipe network model preprocessing method with parameterization and reconfigurable characteristics | |
| CN108595752B (en) | Single-phase water fluid network simulation method oriented to nuclear power plant simulation application | |
| CN111259535B (en) | Method for developing logistics simulation system of ship pipeline production line | |
| CN106682376A (en) | Whole-process steam turbine modeling and recognizing method of actual characteristics of parameters changing with working conditions | |
| CN107330970A (en) | A kind of threedimensional model processing method, server and computer read/write memory medium | |
| CN112069692B (en) | Optimization solving method for natural gas pipe network transmission difference calculation | |
| CN108763741A (en) | A kind of hydraulic hose fluid structurecoupling Numerical Predicting Method | |
| CN112231876B (en) | Pipeline system fluid dynamic model generation system and method based on geometric model | |
| Pacchin et al. | Analysis of non-iterative methods and proposal of a new one for pressure-driven snapshot simulations with EPANET | |
| JP5020338B2 (en) | Initial value generating apparatus and initial value generating method | |
| CN109359433B (en) | Automation method for reactor core refueling design | |
| CN111680370B (en) | Design method and design device for hydraulic valve | |
| CN105447256A (en) | Excitation enhancement simulation genetic optimization method | |
| CN115630479A (en) | Gas tracking calculation method and device in multi-gas-source pipe network, server and storage medium | |
| CN115374527A (en) | Construction method of regional electricity-heat-cold comprehensive energy dynamic simulation system | |
| CN115079592A (en) | Pipe network simulation method for thermodynamic system of ship nuclear power device | |
| Sorocins et al. | District heating simulation model development to solve optimization problems in the market conditions | |
| JP3236835B2 (en) | One-dimensional plant model unsteady flow evaluation method, evaluation system, and storage medium | |
| Letniowski | Compressor station modeling in networks | |
| JPH05313557A (en) | Automatic plant model generating device | |
| Beez et al. | Automatic generation of bond graph models of process plants | |
| Wook et al. | Development of a power plant simulation tool with GUI based on general purpose design software | |
| PARINI | Object oriented modeling and dynamic optimization of energy systems with application to combined cycle power plant start-up | |
| Larsson | District Heating Network Models for Production Planning | |
| Link et al. | Testing power plant control systems in Modelica |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |