CN110059388A - Parameter distribution modeling and simulating method of nuclear power plant's condenser along Working fluid flow direction - Google Patents

Parameter distribution modeling and simulating method of nuclear power plant's condenser along Working fluid flow direction Download PDF

Info

Publication number
CN110059388A
CN110059388A CN201910284801.8A CN201910284801A CN110059388A CN 110059388 A CN110059388 A CN 110059388A CN 201910284801 A CN201910284801 A CN 201910284801A CN 110059388 A CN110059388 A CN 110059388A
Authority
CN
China
Prior art keywords
parameter
condenser
simulation
flow direction
power plant
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.)
Granted
Application number
CN201910284801.8A
Other languages
Chinese (zh)
Other versions
CN110059388B (en
Inventor
彭敏俊
龚梅杰
朱海山
夏庚磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201910284801.8A priority Critical patent/CN110059388B/en
Publication of CN110059388A publication Critical patent/CN110059388A/en
Application granted granted Critical
Publication of CN110059388B publication Critical patent/CN110059388B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

The present invention is to provide a kind of nuclear power plant's condensers along the parameter distribution modeling and simulating method in Working fluid flow direction.According to practical condenser structure parameter, initial parameter, the input parameter for condenser along Working fluid flow direction is determined;According to practical condenser structure and internal flow heat transfer process, actual physics boundary, region division is carried out to condenser simulation model;By nuclear power plant's condenser along Working fluid flow direction mean parameter Simulation Calculation and Integral Processing computation model;Condenser is established along working medium directioin parameter distributed simulation computation model and carries out simulation calculation;Calculating parameter is exported, arbitrary parameter of the acquisition condenser along Working fluid flow direction is distributed and its dynamic change at any time.The present invention is not relying on specific type, apparatus-form, the parameter values of any specific nuclear power plant's condenser, can be modeled and be emulated for various nuclear power plant's condenser apparatus under normal circumstances, have the good scope of application and versatility.

Description

Parameter distribution modeling and simulating method of nuclear power plant's condenser along Working fluid flow direction
Technical field
The present invention relates to one kind to carry out modeling and simulating method to nuclear power plant's condenser.
Background technique
Nuclear power plant's condenser is the important component in secondary circuit steam-water circulation system.Condenser comes from vapour for receiving Steam discharge is carried out cooling condensation, steam generation is entered after feed-water heater, oxygen-eliminating device heating for multiple times by the steam discharge of turbine Device.Cold source of nuclear power plant's condenser as secondary coolant circuit system has important role, the operation shape of condenser to nuclear power plant's operation State will significantly impact the generating efficiency of nuclear power plant.Nuclear power plant's condenser internal physical process is complex, generally for core The Simulation Modeling of power plant condenser is all often to carry out calculating research in the form of lumped-parameter method or control volume divide, The division of these positions and control volume is often fixation, leads to position or any interior location required for being difficult to obtain sometimes Parameter.Meanwhile in these simulation models, due to fixed position and control volume zoned format, condensation can only be often calculated The parameter of device specific position, the case where can not reflecting univers parameter distribution character inside condenser.For different model, parameter There may be nuclear power plant's condensers of variation, it is however generally that the simulation model versatility established is poor, can not carry out simply Modification after be used for different condenser apparatus, often require to re-start modeling and emulate, so that it is imitative to will affect modeling Genuine efficiency extends the modeling and simulating period, will also result in influence for nuclear power plant's total system modeling and simulation.
Summary of the invention
The purpose of the present invention is to provide a kind of nuclear power plant's condensers with the good scope of application and versatility along work The parameter distribution modeling and simulating method of matter flow direction.
The object of the present invention is achieved like this:
Step 1: it according to practical condenser structure parameter, initial parameter, determines for condenser along Working fluid flow direction Input parameter;
Step 2: according to practical condenser structure and internal flow heat transfer process, actual physics boundary, condenser is emulated Model carries out region division;Through nuclear power plant's condenser along Working fluid flow direction mean parameter Simulation Calculation and integral Manage computation model;Condenser is established along working medium directioin parameter distributed simulation computation model and carries out simulation calculation;
Step 3: output calculating parameter, obtain condenser along Working fluid flow direction arbitrary parameter distribution and its at any time Dynamic change.
The present invention may also include:
1. condenser structure parameter specifically includes that
Condenser overall dimensions parameter: condenser height, condenser width, condenser length, inside tube bank dimensional parameters: Heat exchanger tube diameter, heat exchanger tube length, heat-exchanging tube bundle height, heat exchange tube spacing, heat exchange area, hot trap dimensional parameters: hot trap height, Hot trap volume;
Condenser initial parameter includes: steam turbine exhaust pressure, turbine discharge mass dryness fraction, turbine discharge enthalpy, steamer Machine exhaust steam flow, cooling water inlet temperature, cooling water inlet enthalpy, cooling water flow;
The determination refers to for input parameter of the condenser along Working fluid flow direction: condenser overall dimensions parameter and cold Condenser initial parameter is used as the boundary parameter of condenser simulation model simultaneously, by changing condenser boundary ginseng in simulation process Number realizes that condenser dynamic parameter calculates, as nuclear power plant along the input of the parameter distribution modeling and simulating in Working fluid flow direction.
Described in 2. according to practical condenser structure and internal flow heat transfer process, actual physics boundary, it is imitative to condenser True mode carries out region division and specifically includes:
It is specifically divided into shell-side steam condensing zone, the cooling pool in pipe side and hot trap area, these three regions are all in simulation model It is corresponding with a part of region in practical nuclear power plant's condenser apparatus, fluid interchange calculating process and other in each region Region is different, selects several parameter calculating positions according to demand respectively in shell-side steam condensing zone and Guan Ce cooling water area, In hot trap area using hot trap liquid level position as parameter calculating position.
By nuclear power plant's condenser along Working fluid flow direction mean parameter Simulation Calculation and Integral Processing described in 3. Computation model, specifically:
In shell-side steam condensing zone using steam inlet position as simulation modeling starting point, it is with condenser tube row's extreme lower position Simulation modeling terminal, wherein selected any position is parameter position to be asked;In known inlet steam boundary parameter, It is handled with mean parameter computation model and Integral Processing computation model, by partial differential equation to the local derviation number form of spatial position Formula is converted into quantic, converts the partial differential equation about two parameters in spatial position and time to about only with respect to the time ODE, wherein contain the quantic about spatial position, to obtain about required state parameter, packet Include power pressure, working medium steam quality, Temperature of Working different spatial expression formula and emulated;Shell-side steam condensing zone Parameter position arbitrarily to be asked is with any indefinite quantity, to emulate shell-side steam condensing zone distribution parameter estimated performance;
Starting point is modeled for the simulation of domain with cooling water inlet position in the cooling pool in pipe side, with cooling water outlet position Terminal is modeled for the simulation of domain, wherein any position is parameter position to be asked, in known import cooling water boundary parameter, with flat Equal parameter computation model and Integral Processing computation model are handled;For the parameter along shell-side flow direction any position, About in the partial differential equation of time and spatial position, algebra will be converted into about the space position parameter along Working fluid flow direction Form carries out dynamic process to the state parameter of pipe side any position, including pipe side Temperature of Working, pressure and other arbitrary parameters Emulation;The cooling water area in pipe side is supercooling water state, is handled at least two state parameters, the cooling pool in pipe side arbitrarily to Asking parameter position is any indefinite quantity, to emulate Guan Ce cooling water area characteristics of distributed parameters;
Starting point is modeled for the simulation of domain with pipe row's extreme lower position in hot trap area, with the outlet of condenser condensate for the region Simulation modeling terminal, hot trap liquid level position are parameter position to be asked;Mainly hot trap area liquid level is emulated in hot trap area, will be closed It is converted into the ODE about time single parameter in the partial differential equation of two parameters of time and liquid level, thus to hot trap Area's liquid level is emulated.
4. described establish condenser along working medium directioin parameter distributed simulation computation model and carry out simulation calculation, specifically:
It is calculated simultaneously comprising shell-side steam condensation in shell-side steam condensing zone, on-condensible gas calculates and liquid film conductometer It calculates, pipe side is calculated comprising the heat exchange of pipe side single-phase flow simultaneously in cooling pool, metal pipe-wall heat exchange calculates, in the two regions Distinguish while carrying out heat transfer flow calculating according to selected and specified spatial position, according to simultaneously bond area zoned format and Each region calculated relationship constructs three simulation of domain computation models to nuclear power plant's condenser simulation model respectively, is being imitated By successively being called when true calculating, wherein there are parameter transitive relation, shells in shell-side steam condensing zone and pipe side cooling pool Side steam condensing zone and hot trap area collectively form Simulation Calculation there are parameter transitive relation.
5. exporting calculating parameter, arbitrary parameter of the acquisition condenser along Working fluid flow direction is distributed and its dynamic at any time Variation specifically includes:
Constructed condenser parameter distribution simulation model, carries out in a time step to each region calculating section Calling, according to set calculating and calling sequence, realize dynamic real-time simulator calculate, and combine spatial position and time structure At output, according to calculating, to obtain shell-side steam condensing zone temperature, pressure, mass dryness fraction, enthalpy, flow dynamic with spatial position and time State Parameters variation situation, the cooling pool temperature in pipe side, pressure, enthalpy, flow etc. become with the dynamic parameter of spatial position and time Change situation, hot trap area liquid level, temperature, pressure, enthalpy, flow are with spatial position and the dynamic parameter situation of change of time.
The present invention relates to the condenser in greenhouse PWR nuclear power plant secondary coolant circuit system, and in particular to can be right to one kind The emulation mode that nuclear power plant's condenser is effectively modeled, it is special which is based on the fluid interchange of condenser internal working medium Property and actual physics process, can it is good, any position parameter inside condenser is accurately calculated, effectively reflect that nuclear power plant is cold Condenser parameter distribution characteristic and dynamic operational behaviour meet for the Design and optimization of nuclear power plant's condenser, assessment, demonstration etc. no Same emulation demand.
It is different from traditional nuclear power plant's condenser modeling and simulating method, the purpose of the present invention is to provide one kind can be improved Nuclear power plant's design of condenser and computational efficiency, the parameter distribution modeling and simulating method based on internal actual physics process.The modeling Emulation mode and constructed simulation model can effectively improve the modeling and simulating efficiency of nuclear power plant's condenser, based on internal practical Physical process, with nuclear power plant's condenser simulation model along Working fluid flow directioin parameter distribution character can be more accurate it is anti- Answer the parameter distribution situation and dynamic state of parameters situation of change of any position inside practical nuclear power plant's condenser.The present invention is not relying on Specific type, apparatus-form, the parameter values of any specific nuclear power plant's condenser, can be for various nuclear powers under normal circumstances Factory's condenser apparatus is modeled and is emulated, and can be convenient use in nuclear power plant system modeling and simulation, is had good suitable With range and versatility.
Beneficial effects of the present invention are as follows:
The present invention is by being divided into shell-side steam condensing zone, pipe side cooling pool, three, hot trap area for nuclear power plant's condenser Region, three regions are divided in three regions all using the practical fluid interchange process in condenser inside, physical boundary as partitioning standards Condenser inner parameter is not calculated according to along Working fluid flow direction, is met to condenser along any position in Working fluid flow direction The calculating demand of the parameters combination spatial position such as pressure, temperature, flow, enthalpy for setting distribution.It can effectively mention through the invention The design research and development and modeling and simulating efficiency of high nuclear power plant's condenser promote being applicable in for established nuclear power plant's condenser simulation model Range.The simulation model established can satisfy the difference such as the Design and optimization for different type nuclear power plant condenser, assessment and need It asks, simulation result effectively the parameter distribution to nuclear power plant's condenser and dynamic characteristic can carry out simulation calculation, obtain condensation Any position parameter and dynamic state of parameters situation of change inside device.The condenser simulation model established through the invention can be used as A basic unit in nuclear power plant's analogue system uses.
Detailed description of the invention
Fig. 1 nuclear power plant condenser is along Working fluid flow direction modeling and simulating method flow diagram.
The division of Fig. 2 nuclear power plant condenser region, calculating position schematic diagram.
Fig. 3 nuclear power plant condenser Simulation Calculation schematic diagram.
Specific embodiment
It illustrates below and the present invention is described in more detail.
Nuclear power plant's condenser of the invention along Working fluid flow direction parameter distribution modeling and simulating method, first according to reality Condenser structure parameter, initial parameter etc. determine the input parameter for condenser along Working fluid flow direction.Then, according to reality Condenser structure and internal flow heat transfer process, actual physics boundary carry out region division to condenser simulation model.Pass through core Power plant condenser along Working fluid flow direction mean parameter Simulation Calculation and Integral Processing computation model, establish condenser edge Working medium directioin parameter distributed simulation computation model.Finally calculating parameter is exported, obtains condenser along Working fluid flow direction Arbitrary parameter distribution and its dynamic change at any time.Particular content is as follows:
1. practical condenser structure parameter, initial parameter input parameter as condenser, comprising:
Condenser structure parameter is mainly contained about condenser overall dimensions parameter: condenser height, condenser width, Condenser length, inside tube bank dimensional parameters: heat exchanger tube diameter, heat-exchanging tube bundle height, heat exchange tube spacing, is changed at heat exchanger tube length Heat area, hot trap dimensional parameters: hot trap height, hot trap volume.Condenser initial parameter includes: steam turbine exhaust pressure, steamer Machine steam discharge mass dryness fraction, turbine discharge enthalpy, turbine discharge flow, cooling water inlet temperature, cooling water inlet enthalpy, cooling water Flow.These initial parameters are used as the boundary parameter of condenser simulation model simultaneously, by changing condenser in simulation process Boundary parameter realizes that condenser dynamic parameter calculates.Parameter distribution modeling above as nuclear power plant along Working fluid flow direction is imitative The input of true method.
2. being drawn according to condenser structure and internal flow heat transfer process, actual physics boundary to condenser simulation model Divide as follows:
According to practical fluid interchange process and actual boundary inside nuclear power plant's condenser, when carrying out modeling and simulating by its institute Corresponding nuclear power plant's condenser simulation model carries out region division.Specifically it is divided into shell-side steam condensing zone, the cooling pool in pipe side With hot trap area.These three regions can be corresponding with a part of region in practical nuclear power plant's condenser apparatus in simulation model, The fluid interchange calculating process in each region and other regions are different, divide also according to actual physics process, structure Boundary, in this, as the foundation of these three region divisions.Wherein according to the demand of calculating, spatial position required for selecting is carried out Parameter calculates.Several parameter calculating positions are selected according to demand respectively in shell-side steam condensing zone and Guan Ce cooling water area, in heat Well region is using hot trap liquid level position as parameter calculating position.
3. being handled according to mean parameter computation model and Integral Processing computation model, specifically along Working fluid flow direction Are as follows:
In shell-side steam condensing zone using steam inlet position as simulation modeling starting point, it is with condenser tube row's extreme lower position Simulation modeling terminal, wherein selected any position is parameter position to be asked.In known inlet steam boundary parameter, It is handled with mean parameter computation model and Integral Processing computation model, in mathematical model, by by partial differential equation pair The partial derivative form of spatial position is converted into quantic, can be by the partial differential equation about two parameters in spatial position and time It is converted into about the ODE only with respect to the time, wherein the quantic about spatial position is contained, to obtain pass In required state parameter, including power pressure, working medium steam quality, Temperature of Working etc. different spatial expression formula simultaneously It is emulated.Shell-side steam condensing zone parameter position arbitrarily to be asked is with any indefinite quantity, to emulate shell-side steam condensing zone Distribution parameter estimated performance.
Starting point is modeled for the simulation of domain with cooling water inlet position in the cooling pool in pipe side, with cooling water outlet position Terminal is modeled for the simulation of domain, wherein any position is parameter position to be asked, in known import cooling water boundary parameter, with flat Equal parameter computation model and Integral Processing computation model are handled.Likewise, for along shell-side flow direction any position Parameter will turn in the partial differential equation about the time and spatial position about the space position parameter along Working fluid flow direction Quantic is turned to, it can be to the state parameter of pipe side any position, including pipe side Temperature of Working, pressure and other arbitrary parameters Carry out dynamic process simulation.With shell-side different from, the cooling water area in pipe side is supercooling water state, it is therefore desirable to at least Two state parameters, such as pressure and temperature carry out treatment thereto, can emulate to other parameters.Pipe side cools down pool Parameter position arbitrarily to be asked is any indefinite quantity, to emulate Guan Ce cooling water area characteristics of distributed parameters.
Starting point is modeled for the simulation of domain with pipe row's extreme lower position in hot trap area, with the outlet of condenser condensate for the region Simulation modeling terminal, hot trap liquid level position are parameter position to be asked.Its treatment process is similar to two above region, in hot trap area Mainly hot trap area liquid level is emulated, converts the partial differential equation about two parameters of time and liquid level to about time list The ODE of one parameter, to be emulated to hot trap area liquid level.
4. condenser is established along working medium directioin parameter distributed simulation computation model and carries out simulation calculation, specifically:
Contain the calculating of shell-side steam condensation, on-condensible gas calculating and liquid film conductometer simultaneously in shell-side steam condensing zone It calculates, contains pipe side single-phase flow heat exchange calculating in the cooling pool in pipe side simultaneously, metal pipe-wall heat exchange calculates.The two regions These heat transfer flow calculating are distinguished while being carried out in the selected and specified spatial position of interior basis, ensure that simulation calculation when Between on consistency.According to and in conjunction with area above zoned format and each region calculated relationship, nuclear power plant's condenser is emulated Model constructs three simulation of domain computation models respectively, and when carrying out simulation calculation by successively being called, wherein shell-side steams There are parameter transitive relation, shell-side steam condensing zone and hot trap area, there are parameter transmitting to close in vapour condensing zone and pipe side cooling pool System, collectively forms Simulation Calculation.
5. a pair calculating parameter exports, condenser is obtained along the parameter distribution in Working fluid flow direction, specifically:
Constructed condenser parameter distribution simulation model, carries out in a time step to each region calculating section Calling, according to set calculating and calling sequence, realize dynamic real-time simulator calculate, and combine spatial position and time structure At output.Shell-side steam condensing zone temperature, pressure, mass dryness fraction, enthalpy, flow etc. are obtained with spatial position and time according to calculating Dynamic parameter situation of change, the cooling pool temperature in pipe side, pressure, enthalpy, flow etc. are with spatial position and the dynamic parameter of time Situation of change, hot trap area liquid level, temperature, pressure, enthalpy, flow etc. are with spatial position and the dynamic parameter situation of change of time.
In conjunction with Fig. 1 modeling and simulating flow diagram to nuclear power plant's condenser along Working fluid flow direction parameter distribution emulate Modeling method is illustrated.
According to actual physics process inside nuclear power plant's condenser, there are some required input ginsengs when carrying out simulation calculation Number, input parameter mainly includes structural parameters, initial parameter, boundary parameter etc., these parameters are set in nuclear power plant's condenser The parameter that the meter stage must provide.The determination for carrying out condenser simulation model input parameter first, inputs parameter list such as table 1 It is shown.In dynamic simulation, boundary parameter, that is, initial parameter can change.
1 nuclear power plant's condenser simulation model of table inputs table
According to practical fluid interchange process and actual physics boundary inside nuclear power plant's condenser, when being emulated, by it It is divided into shell-side vapour condesation zone, pipe side cooling pool, three, hot trap area region, and determines that the modeling along Working fluid flow direction is imitative True starting point and terminal.In the cooling pool of shell-side steam condensing zone and pipe side, if required for being set according to actual emulation demand Dry calculating position, hot trap area is using hot trap liquid level position as calculating position, as shown in Figure 2.In Fig. 2,1 is shell-side steam condensation Area, 2 be condenser shell-side inlet/turbine discharge, 3 be shell-side steam condensing zone starting point, 4 be shell-side any given position, 5 It is the cooling pool in pipe side for shell-side steam condensing zone terminal, 6,7 be cooling water inlet, 8 be the cooling pool starting point in pipe side, 9 is pipe Side any given position, 10 be the cooling pool terminal in pipe side, 11 be cooling water outlet, 12 be hot trap area, 13 be hot trap liquid level, 14 It is condensation-water drain for condenser overall length, 15,16 be shell-side flow direction, 17 is pipe side flow direction.
After carrying out region division and calculating position setting, according to mean parameter computation model and Integral Processing computation model Processing emulation is carried out to required parameter, procedure parameter is handled using mean parameter computation model, at integral Reason computation model handles state parameter, is classified as follows shown in table:
2 condenser procedure parameter of table and state parameter
Mean parameter computation model is as follows:
For the space average parameter along Working fluid flow direction, from working medium entrance location to region in any given position into Row average integral has following processing form:
In formula:To indicate that process averages parameter, the p such as enthalpy, density indicate the actual procedure parameter in position to be asked, LiFor Position along Working fluid flow direction apart from import, for parameter position to be asked, z is Working fluid flow direction coordinate;In subscript, 0 is indicated Zone inlet position, i are indicated apart from imported LiAny position, for indicating condenser interior zone any position.
Integral Processing computation model is as follows:
For the state parameter along Working fluid flow direction any position, from working medium entrance location to region in any position, By original partial differential equation from 0 position (entry position) to the position i (calculating position), carried out at integral for spatial position Reason, is converted into a certain state parameter only with respect to the ODE of time parameter, wherein including the algebra shape about spatial position Formula, processing form are as follows:
In formula: f is about state parameters such as pressure, mass dryness fraction, temperature, liquid levels, and t is simulation time, and Δ t is time step, n For time step calculation times, p1, p2... ... it is calculating process parameter.
In two above computation model, convert the partial differential equation about two parameters of time and spatial position to only About the ODE of time, wherein by quantic is converted into about the space position parameter along Working fluid flow direction, it is right Dynamic process simulation is carried out to any required parameter in any required position, get parms the output changed over time.It is cold The parameter of any position all obtains specific parameter expression form using above-mentioned calculation method in three regions inside condenser.
The calculating of shell-side steam condensation is carried out simultaneously in shell-side steam condensing zone, on-condensible gas calculates, liquid film heat conduction, pipe Single-phase flow calculating, metal pipe-wall heat conduction in pipe are carried out in the cooling pool in side simultaneously.Hot trap area is related to vehicle repair major volume It is calculated with liquid level.Thus three regions divided can the band of position corresponding with practical nuclear power plant's condenser it is corresponding, simultaneously Simulation calculation in each region is in accordance with practical nuclear power plant's condenser internal physical process and carries out, and all calculating are at one Between in step-length respectively simultaneously, synchronous carry out.In conjunction with input parameter, mean parameter computation model, Integral Processing computation model, to shell The cooling pool of side steam condensing zone and pipe side is emulated.For the cooling pool of shell-side steam condensing zone and pipe side, according to emulation Demand, set in condenser simulation model along Working fluid flow direction several positions and set, to wherein procedure parameter and State parameter is emulated.The parameter distribution situation along Working fluid flow direction is obtained by simulation calculation.Change in boundary condition In the case where change, the parameter calculating of different location can all be changed correspondingly, to complete the dynamic of condenser inner parameter at any time Simulation process.Shell-side steam condensing zone, the cooling pool in pipe side and the hot trap divided according to actual physics process inside condenser Three, area region, in conjunction with above-mentioned mean parameter computation model and Integral Processing computation model, wherein to steam condensing zone, hydrophobic Cooling zone calculate along the parameters simulation in Working fluid flow direction, is substantially carried out vehicle repair major liquid level simulating to hot trap area and calculates, In conjunction with the heat transfer flow calculated relationship of each region, collectively as Simulation Calculation, each region calculate between there are parameters Transitive relation, specific calculation process and parameter transitive relation are as shown in Figure 3.
Constructed condenser parameter distribution Simulation Calculation carries out calculating each region in a time step Partial calling, according to set calculating and calling sequence, dynamic real-time simulator calculating can be carried out, in conjunction with spatial position and when Between constitute output.Required parameter etc. is obtained with spatial position and the dynamic parameter situation of change of time, parameter according to calculating It is as shown in table 3 to export table.
3 nuclear power plant's condenser simulation model of table exports table
For with different form, different structure, different parameters nuclear power plant's condenser, the present invention can be to these feelings Condenser with identical fluid interchange process under condition carries out simulation modeling, while having the input parameter, just of the same category Beginning parameter etc., computation sequence and parameter transitive relation are still consistent.Mould is emulated by parameter distribution nuclear power plant condenser The output result that type calculates can in real time, accurately react the parameter distribution situation of corresponding region in practical nuclear power plant's condenser, Especially steam condensing zone and drain cooling zone, calculating output result carry out any along Working fluid flow direction position as needed Parameter calculate, can more suitably react the dynamic operational behaviour and parameter distribution situation of nuclear power plant's condenser.
It should be noted that the above embodiments are only used to illustrate the present invention, and not limitation of the present invention, related skill The those of ordinary skill in art field, without departing from the spirit and scope of the present invention, can also make a variety of changes and Modification, therefore all equivalent technical solutions also belong to scope of the invention, scope of patent protection of the invention should be wanted by right Ask restriction.

Claims (6)

1. a kind of nuclear power plant's condenser is along the parameter distribution modeling and simulating method in Working fluid flow direction, it is characterized in that:
Step 1: according to practical condenser structure parameter, initial parameter, the input for condenser along Working fluid flow direction is determined Parameter;
Step 2: according to practical condenser structure and internal flow heat transfer process, actual physics boundary, to condenser simulation model Carry out region division;By nuclear power plant's condenser along Working fluid flow direction mean parameter Simulation Calculation and Integral Processing meter Calculate model;Condenser is established along working medium directioin parameter distributed simulation computation model and carries out simulation calculation;
Step 3: output calculating parameter obtains condenser along the arbitrary parameter distribution in Working fluid flow direction and its moving at any time State variation.
2. nuclear power plant's condenser according to claim 1 is along the parameter distribution modeling and simulating method in Working fluid flow direction, It is characterized in that condenser structure parameter specifically includes that
Condenser overall dimensions parameter: condenser height, condenser width, condenser length, inside tube bank dimensional parameters: heat exchange Pipe diameter, heat exchanger tube length, heat-exchanging tube bundle height, heat exchange tube spacing, heat exchange area, hot trap dimensional parameters: hot trap height, hot trap Volume;
Condenser initial parameter includes: steam turbine exhaust pressure, turbine discharge mass dryness fraction, turbine discharge enthalpy, steam turbine row Steam flow amount, cooling water inlet temperature, cooling water inlet enthalpy, cooling water flow;
The determination refers to for input parameter of the condenser along Working fluid flow direction: condenser overall dimensions parameter and condenser Initial parameter is used as the boundary parameter of condenser simulation model simultaneously, by changing condenser boundary parameter in simulation process, Realize that condenser dynamic parameter calculates, as nuclear power plant along the input of the parameter distribution modeling and simulating in Working fluid flow direction.
3. nuclear power plant's condenser according to claim 1 is along the parameter distribution modeling and simulating method in Working fluid flow direction, It is characterized in described according to practical condenser structure and internal flow heat transfer process, actual physics boundary, mould is emulated to condenser Type carries out region division and specifically includes:
Specifically be divided into shell-side steam condensing zone, the cooling pool in pipe side and hot trap area, in simulation model these three regions all with reality A part of region in border nuclear power plant condenser apparatus is corresponding, fluid interchange calculating process and other regions in each region It is different, several parameter calculating positions are selected according to demand respectively in shell-side steam condensing zone and Guan Ce cooling water area, in heat Well region is using hot trap liquid level position as parameter calculating position.
4. nuclear power plant's condenser according to claim 1 is along the parameter distribution modeling and simulating method in Working fluid flow direction, It is characterized in the calculating along Working fluid flow direction mean parameter Simulation Calculation and Integral Processing by nuclear power plant's condenser Model, specifically:
It is emulation with condenser tube row's extreme lower position in shell-side steam condensing zone using steam inlet position as simulation modeling starting point Terminal is modeled, wherein selected any position is parameter position to be asked;In known inlet steam boundary parameter, with flat Equal parameter computation model and Integral Processing computation model are handled, and partial derivative form of the partial differential equation to spatial position is turned Quantic is turned to, converts the partial differential equation about two parameters in spatial position and time to about normal only with respect to the time The differential equation, wherein the quantic about spatial position is contained, to obtain about required state parameter, including work Matter pressure, working medium steam quality, Temperature of Working different spatial expression formula and emulated;Shell-side steam condensing zone is any Parameter position to be asked is with any indefinite quantity, to emulate shell-side steam condensing zone distribution parameter estimated performance;
Starting point is modeled for the simulation of domain with cooling water inlet position in the cooling pool in pipe side, being with cooling water outlet position should Simulation of domain models terminal, and wherein any position is parameter position to be asked, in known import cooling water boundary parameter, with average ginseng Number computation model and Integral Processing computation model are handled;For the parameter along shell-side flow direction any position, about In the partial differential equation of time and spatial position, algebra shape will be converted into about the space position parameter along Working fluid flow direction Formula carries out dynamic process and imitates to the state parameter of pipe side any position, including pipe side Temperature of Working, pressure and other arbitrary parameters Very;The cooling water area in pipe side is supercooling water state, is handled at least two state parameters, and the cooling pool in pipe side is arbitrarily wait ask Parameter position is any indefinite quantity, to emulate Guan Ce cooling water area characteristics of distributed parameters;
Starting point is modeled for the simulation of domain with pipe row's extreme lower position in hot trap area, with the outlet of condenser condensate for the simulation of domain Terminal is modeled, hot trap liquid level position is parameter position to be asked;Mainly hot trap area liquid level is emulated in hot trap area, will about when Between and the partial differential equation of two parameters of liquid level be converted into the ODE about time single parameter, thus to hot trap area liquid Position is emulated.
5. nuclear power plant's condenser according to claim 1 is along the parameter distribution modeling and simulating method in Working fluid flow direction, Be characterized in it is described establish condenser along working medium directioin parameter distributed simulation computation model and carry out simulation calculation, specifically:
It is calculated simultaneously comprising shell-side steam condensation in shell-side steam condensing zone, on-condensible gas calculates and liquid film heat conduction, pipe Side is calculated comprising the heat exchange of pipe side single-phase flow simultaneously in cooling pool, metal pipe-wall heat exchange calculates, basis in the two regions Heat transfer flow calculating is distinguished while being carried out in selected and specified spatial position, according to simultaneously bond area zoned format and each Region calculated relationship constructs three simulation of domain computation models to nuclear power plant's condenser simulation model respectively, is carrying out emulation meter By successively being called when calculation, wherein there are parameter transitive relation, shell-side steams in shell-side steam condensing zone and pipe side cooling pool Vapour condensing zone and hot trap area collectively form Simulation Calculation there are parameter transitive relation.
6. nuclear power plant's condenser according to claim 1 is along the parameter distribution modeling and simulating method in Working fluid flow direction, It is characterized in exporting calculating parameter, arbitrary parameter of the acquisition condenser along Working fluid flow direction is distributed and its dynamic change at any time It specifically includes:
Constructed condenser parameter distribution simulation model, carries out the tune to each region calculating section in a time step With according to set calculating and calling sequence, realization dynamic real-time simulator is calculated, and it is defeated to combine spatial position and time to constitute Out, shell-side steam condensing zone temperature, pressure, mass dryness fraction, enthalpy, flow is obtained according to calculating to join with the dynamic of spatial position and time Number situation of change, the cooling pool temperature in pipe side, pressure, enthalpy, flow etc. change feelings with the dynamic parameter of spatial position and time Condition, hot trap area liquid level, temperature, pressure, enthalpy, flow are with spatial position and the dynamic parameter situation of change of time.
CN201910284801.8A 2019-04-10 2019-04-10 Parameter distribution modeling simulation method for condenser of nuclear power plant along working medium flowing direction Active CN110059388B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910284801.8A CN110059388B (en) 2019-04-10 2019-04-10 Parameter distribution modeling simulation method for condenser of nuclear power plant along working medium flowing direction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910284801.8A CN110059388B (en) 2019-04-10 2019-04-10 Parameter distribution modeling simulation method for condenser of nuclear power plant along working medium flowing direction

Publications (2)

Publication Number Publication Date
CN110059388A true CN110059388A (en) 2019-07-26
CN110059388B CN110059388B (en) 2023-01-03

Family

ID=67317519

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910284801.8A Active CN110059388B (en) 2019-04-10 2019-04-10 Parameter distribution modeling simulation method for condenser of nuclear power plant along working medium flowing direction

Country Status (1)

Country Link
CN (1) CN110059388B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113821871A (en) * 2021-08-31 2021-12-21 中国船舶重工集团公司第七0三研究所 Dynamic simulation model development method for multi-structure and multi-pressure-node heat exchanger of ship two-loop system
CN114239432A (en) * 2021-12-02 2022-03-25 中国船舶重工集团公司第七0三研究所 Modelica-based dynamic simulation model development method for heat exchanger of ship two-loop system
WO2023093490A1 (en) * 2021-11-29 2023-06-01 云南电网有限责任公司电力科学研究院 High-temperature steam heat pump dynamic modeling method and system

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4042813A (en) * 1973-02-23 1977-08-16 Westinghouse Electric Corporation Secondary system modeling and method for a nuclear power plant training simulator
US20050165591A1 (en) * 2004-01-27 2005-07-28 Vijay Bahel Cooling system design simulator
CN102279901A (en) * 2011-05-17 2011-12-14 湖北省电力公司电力试验研究院 Modeling method specific to third generation pressurized water reactor nuclear power generating unit
CN104049539A (en) * 2014-05-29 2014-09-17 华中科技大学 Full-work-condition simulation system of turbine regenerative system
CN106528985A (en) * 2016-11-03 2017-03-22 哈尔滨工程大学 Partitioned simulation method for nuclear power unit condenser
CN107704674A (en) * 2017-09-26 2018-02-16 吉林省电力科学研究院有限公司 The method for numerical simulation of air cooling tubes condenser water vapor condensation process
CN107798167A (en) * 2017-09-21 2018-03-13 东南大学 Direct Air-Cooled generating set cold end system modeling and optimization method
CN108875151A (en) * 2018-05-07 2018-11-23 哈尔滨工程大学 Nuclear power secondary coolant circuit system parameterized simulation model rapid build system and construction method

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4042813A (en) * 1973-02-23 1977-08-16 Westinghouse Electric Corporation Secondary system modeling and method for a nuclear power plant training simulator
US20050165591A1 (en) * 2004-01-27 2005-07-28 Vijay Bahel Cooling system design simulator
CN102279901A (en) * 2011-05-17 2011-12-14 湖北省电力公司电力试验研究院 Modeling method specific to third generation pressurized water reactor nuclear power generating unit
CN104049539A (en) * 2014-05-29 2014-09-17 华中科技大学 Full-work-condition simulation system of turbine regenerative system
CN106528985A (en) * 2016-11-03 2017-03-22 哈尔滨工程大学 Partitioned simulation method for nuclear power unit condenser
CN107798167A (en) * 2017-09-21 2018-03-13 东南大学 Direct Air-Cooled generating set cold end system modeling and optimization method
CN107704674A (en) * 2017-09-26 2018-02-16 吉林省电力科学研究院有限公司 The method for numerical simulation of air cooling tubes condenser water vapor condensation process
CN108875151A (en) * 2018-05-07 2018-11-23 哈尔滨工程大学 Nuclear power secondary coolant circuit system parameterized simulation model rapid build system and construction method

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
GE BIN, ET AL.: "Modeling and dynamic simulation of nuclear power plant secondary loop", 《ELECTRIC POWER AUTOMATION EQUIPMENT》 *
GONG MEIJIE, ET AL.: "Research of fast modeling and simulating platform for nuclear power plant secondary loop", 《PROCEEDINGS OF THE 26TH INTERNATIONAL CONFERENCE ON NUCLEAR ENGINEERING》 *
JI-LIN SUN, ET AL.: "Investigation of the thermal characteristics of condensers in nuclear power plant by simulation with zoning model", 《ANNALS OF NUCLEAR ENERGY》 *
吴鹏,等: "核电站凝汽器仿真模型研究", 《原子能科学技术》 *
吴鹏: "凝汽器分布参数动态仿真模型研究", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技II辑》 *
孔夏明,等: "蒸汽排放***蒸汽冷凝器动态特性仿真研究", 《原子能科学技术》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113821871A (en) * 2021-08-31 2021-12-21 中国船舶重工集团公司第七0三研究所 Dynamic simulation model development method for multi-structure and multi-pressure-node heat exchanger of ship two-loop system
WO2023093490A1 (en) * 2021-11-29 2023-06-01 云南电网有限责任公司电力科学研究院 High-temperature steam heat pump dynamic modeling method and system
CN114239432A (en) * 2021-12-02 2022-03-25 中国船舶重工集团公司第七0三研究所 Modelica-based dynamic simulation model development method for heat exchanger of ship two-loop system

Also Published As

Publication number Publication date
CN110059388B (en) 2023-01-03

Similar Documents

Publication Publication Date Title
CN105224735B (en) Generating set energy efficiency analysis method for air
CN110059388A (en) Parameter distribution modeling and simulating method of nuclear power plant's condenser along Working fluid flow direction
CN106682376B (en) Whole-process steam turbine modeling and identification method for actual characteristics of parameters changing along with working conditions
CN105303032B (en) Influence the objective factor analysis method of generating set efficiency
CN104834773A (en) Simulation method for heat exchange performance of straight tube type once-through steam generator
CN104049539B (en) A kind of full working scope analogue system of Steam Turbine Regenerative System
CN104749999B (en) The Turbo-generator Set cold end system optimization operation of assembling wet cooling tower accurately instructs system
CN106528985B (en) A kind of partitioning emulation mode of nuclear power unit condenser
Gao et al. Performance prediction of wet cooling tower using artificial neural network under cross-wind conditions
CN107763610B (en) A kind of structure-based nuclear power plant's high-pressure feed-water heater fining emulation mode
CN104777008B (en) A kind of power-plant flue gas afterheat utilizing system performance simulation experimental apparatus for testing
CN106960113A (en) A kind of divisions of responsibility method on single shaft combined cycle generating unit performances acceptance Zhong Ji Dao Yulu islands
CN105184087B (en) The computational methods and device that variation of ambient temperature influences Thermal generation unit coal consumption
CN109141541A (en) A kind of coal-fired power station boiler reheater working medium flow on-line correction method
Silva et al. Evaluation of the efficiency of a heat recovery steam generator via computational simulations of off-design operation
CN105447256B (en) A kind of enhancing excitation emulation genetic optimization method
CN204514629U (en) A kind of power-plant flue gas bootstrap system performance simulation experimental apparatus for testing
CN106126822A (en) The steam-water separator modeling method of a kind of direct current cooker and system
Laskowski et al. A useful formulas to describe the performance of a steam condenser in off-design conditions
CN109858073A (en) For the real-time quick construction method based on transient test data of plate-fin heat exchanger efficiency Model
CN107247843A (en) The modeling method of solar energy thermal-power-generating station oil-water type model steam generator
CN108595813B (en) A kind of CO2 gas cooler optimum design method based on optimum quality flow velocity
CN116484524A (en) Rapid simulation and design method for performance of multi-row finned tube heat exchanger
CN113609619A (en) Multidimensional coupling simulation method for long blade blast of low-pressure through-flow area of steam turbine
CN114004159A (en) Variable working condition efficiency calculation method for waste heat boiler superheater

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