CN106528985B - A kind of partitioning emulation mode of nuclear power unit condenser - Google Patents

A kind of partitioning emulation mode of nuclear power unit condenser Download PDF

Info

Publication number
CN106528985B
CN106528985B CN201610953179.1A CN201610953179A CN106528985B CN 106528985 B CN106528985 B CN 106528985B CN 201610953179 A CN201610953179 A CN 201610953179A CN 106528985 B CN106528985 B CN 106528985B
Authority
CN
China
Prior art keywords
subregion
condenser
pipe side
hot well
pipe
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
Application number
CN201610953179.1A
Other languages
Chinese (zh)
Other versions
CN106528985A (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 CN201610953179.1A priority Critical patent/CN106528985B/en
Publication of CN106528985A publication Critical patent/CN106528985A/en
Application granted granted Critical
Publication of CN106528985B publication Critical patent/CN106528985B/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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Abstract

The present invention is to provide a kind of partitioning emulation modes of nuclear power unit condenser.(1) subregion division is carried out to simulation object, the simulation object includes conventional condenser and common condenser, and it includes uniform segmentation and non-homogeneous subregion that the subregion, which divides,;(2) it inputs the structure size of condenser, physical property, boundary parameter, calculate initial value;(3) Calculation of Heat Transfer is carried out to the one-to-one subregion of shell-side and pipe side;(4) thermal parameter calculating is carried out;(5) simulation result exports.Output data of the simulation model established by the present invention under stable state and transient condition can satisfy the precise requirements of device emulation compared with physical device, can use this emulation mode and carries out operation characteristic research.

Description

A kind of partitioning emulation mode of nuclear power unit condenser
Technical field
The present invention relates to a kind of emulation mode, specifically a kind of emulation mode of nuclear power unit condenser.
Background technique
Nuclear power unit condenser is one of the capital equipment in secondary circuit, is responsible for the wet steaming finished of doing work in steam turbine Vapour is condensed into condensate, and the heat transfer that condensation is released is to recirculated water.Nuclear power unit use condenser be divided into shell-side with Pipe side two parts, are physically separated by heat exchanger tube.The steam and other steam of Steam Turbine discharge are passed since temperature is higher than Heat pipe outer wall face temperature, therefore heat transfer phenomenon occurs, the latent heat of vaporization is discharged by liquefaction, transfers heat to heat exchange pipe outer wall; Since heat exchanger tube outside wall surface temperature is higher than inner wall temperature, Heat Conduction Phenomenon occurs, heat is transmitted to from heat-transfer pipe outside wall surface Heat-transfer pipe inner wall;Since heat-transfer pipe inner wall temperature is higher than circulating water temperature, continue that forced-convection heat transfer phenomenon occurs, Heat is passed into recirculated water from heat transfer inside pipe wall, and nuclear power unit is taken out of by recirculated water, completes the transfer process of energy.Such as Fruit condenser has multiple air intakes, is responsible for the condensation duty of different turbine discharges, then referred to as common condenser.Common condenser Internal flowing and heat exchange mode is more complicated.
With advances in technology, multiple sensors have been arranged in nuclear power unit condenser at present and have carried out pressure survey Amount, the setting position of different sensors, the pressure received is not also identical.Condenser emulation mode based on lumped-parameter method without The subtleties of method consideration equipment internal pressure;Subregion is carried out for axial in current existing subregion emulation mode, is only capable of body Existing sequence of flow, cannot embody the difference of different location thermal parameter.In common condenser, different air intakes, bleeding point Pressure difference is more important parameter, and Traditional calculating methods cannot be distinguished from calculating.Therefore, traditional simulation modeling method without Method meets the needs of high-precision real-time simulation.
Summary of the invention
The purpose of the present invention is to provide a kind of partitioning emulation modes of nuclear power unit condenser that accuracy is high.
The object of the present invention is achieved like this:
(1) subregion division is carried out to simulation object, the simulation object includes conventional condenser and common condenser, described It includes uniform segmentation and non-homogeneous subregion that subregion, which divides,;
(2) it inputs the structure size of condenser, physical property, boundary parameter, calculate initial value;
(3) Calculation of Heat Transfer is carried out to the one-to-one subregion of shell-side and pipe side;
(4) thermal parameter calculating is carried out;
(5) simulation result exports.
The present invention may also include:
1, described that simulation object progress subregion division is specifically included:
A. the non-homogeneous subregion of conventional condenser is divided: subregion is carried out according to the position of practical condenser sensor It divides, the subregion of division includes turbine discharge mouth region domain, air ejector region, hot well region, sensor region, except above-mentioned Other shell-side regions outside region, pipe side circulating water intake water chamber, pipe side circulating water outlet water chamber;Except hot well region, pipe side are followed Outside ring water inlet water chamber, pipe side circulating water outlet water chamber, other each shell-sides are overlapped with pipe side subregion physical region, are corresponded, It should all include portion of heat transfer pipeline;
B. the uniform segmentation of conventional condenser is divided: condenser is divided into addition to hot well 2 × 2,3 × 3 or 4 × 4 subregions, each subregion volume is identical, and includes the heat transfer pipe of identical quantity;And by turbine discharge mouth region domain, pumping Device region, sensor region are divided into wherein some subregion respectively;Hot well part be individually divided into a subregion, pipe side recirculated water into Saliva room is individually divided into a subregion, pipe side circulating water outlet water chamber is individually divided into a subregion;Except hot well subregion, pipe side Outside circulating water intake water chamber, pipe side circulating water outlet water chamber, other each shell-sides are overlapped with pipe side subregion physical region, and one is a pair of It answers;
C. the non-homogeneous subregion of common condenser is divided: according to the position of practical condenser inner sensor and master The position progress subregion division of machine, subsidiary engine air intake, the subregion of division include host steam discharge region, subsidiary engine steam discharge region, pumping Device region, hot well region, sensor region, other shell-side regions in addition to above-mentioned zone, pipe side circulating water intake water chamber, Pipe side circulating water outlet water chamber;It is outer in addition to hot well region, pipe side circulating water intake water chamber, pipe side circulating water outlet water chamber, other Each shell-side is overlapped with pipe side subregion physical region, is corresponded, should all be included portion of heat transfer pipeline;
D. common condenser uniform segmentation is divided: common condenser is divided into 3 × 3 or 4 × 4 points in addition to hot well Area, each subregion volume is identical, and includes the heat transfer pipe of identical quantity;And by host steam drain region, subsidiary engine steam drain Region, air ejector region, sensor region are divided into wherein some subregion respectively;Subregion divide number with host steam drain region with Subsidiary engine steam drain region is not subject to same subregion;Hot well part is individually divided into a subregion, pipe side circulating water intake water chamber Individually it is divided into a subregion, pipe side circulating water outlet water chamber is individually divided into a subregion;Except hot well subregion, pipe side recirculated water Outside import water chamber subregion, pipe side circulating water outlet water chamber subregion, other each shell-sides are overlapped with pipe side subregion physical region, and one is a pair of It answers.
2, the structure size, physical property of the input condenser, boundary parameter, calculating initial value.It specifically includes:
A. the input of full mold structured data includes with the outer dimension of the nuclear power unit condenser of International System of Units description, pipe Road internal diameter, outer diameter tube, duct length, flow resistance data;
B. the input of integer structured data includes nuclear power unit condenser tube number of rows, process number, heat-transfer pipe arrangement mode;Its In, heat-transfer pipe arrangement mode thinks that value is 1 when being triangle, be 0 when square;
C. the input of running boundary data includes pipe side circulating water intake flow, temperature, atmospheric pressure, air ejector pumping pressure Power, pressure before condensate pump;It further include turbine discharge flow, pressure, mass dryness fraction for conventional condenser;For common condenser, It further include host exhaust steam flow, pressure, mass dryness fraction and subsidiary engine exhaust steam flow, pressure, mass dryness fraction;
D. calculating initial value input includes that described calculate uses initial value, specially pipe side and each by stages flow of shell-side, initial pressure Power, initial temperature, rate of setting.
3, the Calculation of Heat Transfer specifically includes:
A. the input of full mold structured data is started in calculating input included data with integer structured data;Repeatedly It is additional to input included data using calculating initial value when first time during generation calculates;
B. Nu Saier film condensation calculation formula is utilized, the amendment of adding tube number of rows, the amendment of process number calculate and remove hot well subregion Outer each shell-side subregion condensation heat transfer and steam condensation rate;
C. forced-convection heat transfer D-B formula is utilized, addition varied property amendment calculates forced-convection heat transfer in the subregion of each pipe side Amount and recirculated water temperature rise;
D. by each shell-side subregion heat-transfer pipe outer wall of metal heat-conducting equation calculation towards pipe side subregion heat-transfer pipe, that is, inner wall Arrive heat.
4, the progress thermal parameter calculating specifically includes:
A. it is directed to shell-side subregion, in addition to hot well, each subregion steam pressure is solved using quality of steam equilibrium equation;For heat Well subregion is calculated liquid level using water quality equilibrium equation, is asked using quality of steam equilibrium equation and water quality equilibrium equation simultaneous Solve hot well zoned pressure;
B. it is directed to each subregion in pipe side, the circulating water pressure of each subregion is calculated using water quality equilibrium equation;
C. it is directed to shell-side subregion, in addition to hot well, using saturation temperature as zone temperature;Hot well vapor (steam) temperature is saturation temperature Degree, coolant-temperature gage are condensate mixing temperature;
D. it is directed to pipe side subregion, circulating water temperature is determined using pressure and enthalpy.
The present invention provides a kind of partitioning emulation modes of nuclear power unit condenser, pass through the specific of analysis condenser Structure arrangement, needs in conjunction with emulation, and the uniform segmentation that can carry out conventional condenser and common condenser divides and non-homogeneous point Division;It is defeated by full mold structured data such as outer dimension input, integer structured data such as process number after the completion of subregion divides Enter, the input of running boundary data such as circulating water flow and the input of calculating initial value such as each zoned pressure, is assigned to variable is calculated Value;The meter of rate of setting and heat exchange amount is carried out to each shell-side heat exchange subregion by the methods of Nu Saier film condensation calculation formula It calculates, heat exchange calculating is carried out to the pipe side subregion that respectively exchanges heat by forced-convection heat transfer D-B formula etc., by metal heat-conducting equation to each Exchanging heat, subregion heat-transfer pipe is thermally conductive to be calculated;Thermal parameter calculating is carried out to all subregions using mass balance method;Calculate knot Fruit is exported by output program.Emulation mode provided by the present invention divides, input parameter, heat exchange calculating, thermal technology comprising subregion Several important steps such as parameter calculates, result exports, output data of the simulation model established under stable state and transient condition It compared with physical device, can satisfy the precise requirements of device emulation, can use this emulation mode progress operation characteristic and grind Study carefully.
The present invention can carry out the emulation different from existing way to condenser, and main advantages of the present invention are embodied in:
(1) present invention uses partitioned mode, and the radial direction of condenser is divided into different subregions, can calculate condenser not With the different parameters at position;
(2) differential method of pressure solves vapor partial pressure using density, accuracy relatively utilizes perfect gas side Method solution is improved;
(3) different external disturbances correspond to different subregions, solve common condenser in calculating process, since difference is into vapour The different computational problem of back pressure caused by mouth throttle flow is different;
(4) method proposed can be very good to realize the real-time simulated animation of condenser inner parameter, may be implemented pair Pressure change in air ejector failure, common condenser when a certain mutation into vapour inside condenser calculates.
Detailed description of the invention
Emulation mode flow diagram Fig. 1 of the invention.
The non-homogeneous subregion emulation mode schematic diagram of Fig. 2 conventional condenser.
Fig. 3 common condenser uniform segmentation emulation mode schematic diagram.
Specific embodiment
Below in conjunction with attached drawing, the present invention is described further.
In conjunction with Fig. 1, realization step of the invention is specifically included that
(1) subregion division is carried out for simulation object, conventional condenser, common condenser two can be divided by simulation object Kind;Two kinds of uniform segmentation, non-homogeneous subregion can be divided by subregion zoned format.Specific implementation method is determined such that
A. it is directed to the non-homogeneous subregion splitting scheme of conventional condenser, is carried out according to the position of practical condenser sensor Subregion divides, and the subregion of division includes 11 turbine discharge mouth region domain of subregion, 12 air ejector region of subregion, 13 hot well area of subregion Domain, sensor region, other shell-side regions in addition to above-mentioned zone, pipe side circulating water intake water chamber, pipe side recirculated water go out Saliva room.In above-mentioned zone, according to the difference of simulation requirements, sensor region can for one, or it is multiple, It can merge with other subregions;In addition to hot well region, pipe side circulating water intake water chamber, pipe side circulating water outlet water chamber, other are each Shell-side is overlapped with pipe side subregion physical region, is corresponded, should all be included portion of heat transfer pipeline;
B. it is directed to the uniform segmentation splitting scheme of conventional condenser, according to the size of practical condenser cross section and emulation It is required that condenser to be divided into 2 × 2,3 × 3 or 4 × 4 subregions in addition to hot well, each subregion volume is identical, it is believed that packet Heat transfer pipe containing identical quantity;And turbine discharge mouth region domain, air ejector region, sensor region are divided into wherein respectively Some subregion;Hot well part is individually divided into a subregion, pipe side circulating water intake water chamber is individually divided into a subregion, pipe side Circulating water outlet water chamber is individually divided into a subregion;Except hot well subregion, pipe side circulating water intake water chamber, pipe side circulating water outlet Outside water chamber, other each shell-sides are overlapped with pipe side subregion physical region, are corresponded;
C. it is directed to the non-homogeneous subregion splitting scheme of common condenser, according to the position of practical condenser inner sensor And the position of host, subsidiary engine air intake carries out subregion division, the subregion of division includes host steam discharge region, subsidiary engine steam discharge area Domain, air ejector region, hot well subregion, sensor region, other shell-side regions in addition to above-mentioned zone, pipe side recirculated water into Saliva room, pipe side circulating water outlet water chamber.In above-mentioned zone, according to the difference of simulation requirements, sensor region can be One, or it is multiple, can also merge with other subregions;Except hot well subregion, pipe side circulating water intake water chamber, pipe side circulation Outside water out water chamber, other each shell-sides are overlapped with pipe side subregion physical region, are corresponded, should all be included portion of heat transfer pipeline;
D. it is directed to common condenser uniform segmentation splitting scheme, according to the size and emulation of practical common condenser cross section Requirement, common condenser is divided into 3 × 3 or 4 × 4 subregions in addition to hot well, each subregion volume is identical, it is believed that includes There is the heat transfer pipe of identical quantity;And by host steam drain region, subsidiary engine steam drain region, air ejector region, sensor region It is divided into wherein some subregion respectively;Subregion divides number with host steam drain region and subsidiary engine steam drain region not in same subregion Subject to;Hot well part is individually divided into a subregion, pipe side circulating water intake water chamber is individually divided into a subregion, pipe side circulation Water out water chamber is individually divided into a subregion;Except hot well subregion, pipe side circulating water intake water chamber subregion, pipe side circulating water outlet Outside water chamber subregion, other each shell-sides are overlapped with pipe side subregion physical region, are corresponded;
(2) data input.Specifically it is implemented:
A. the input of full mold structured data includes with the outer dimension of the nuclear power unit condenser of International System of Units description, pipe Road internal diameter, outer diameter tube, duct length, flow resistance data;
B. the input of integer structured data includes nuclear power unit condenser tube number of rows, process number, heat-transfer pipe arrangement mode;Its In, heat-transfer pipe arrangement mode thinks that value is 1 when being triangle, be 0 when square;
C. the input of running boundary data includes pipe side circulating water intake flow, temperature, atmospheric pressure, air ejector pumping pressure Power, pressure before condensate pump;It further include turbine discharge flow, pressure, mass dryness fraction for conventional condenser;For common condenser, It further include host exhaust steam flow, pressure, mass dryness fraction and subsidiary engine exhaust steam flow, pressure, mass dryness fraction;
D. calculating initial value input includes that described calculate uses initial value, specially pipe side and each by stages flow of shell-side, initial pressure Power, initial temperature, rate of setting.
(3) Calculation of Heat Transfer process is carried out.The process is just for shell-side and the one-to-one subregion in pipe side.Specific implementation method It is such that
A. the input of full mold structured data is started in calculating input included data with integer structured data;Repeatedly It is additional to input included data using calculating initial value when first time during generation calculates;
B. Nu Saier film condensation calculation formula is utilized, the amendment of adding tube number of rows, the amendment of process number calculate and remove hot well subregion Outer each shell-side subregion condensation heat transfer and steam condensation rate;
C. forced-convection heat transfer D-B formula is utilized, addition varied property amendment calculates forced-convection heat transfer in the subregion of each pipe side Amount and recirculated water temperature rise;
D. by each shell-side subregion heat-transfer pipe outer wall of metal heat-conducting equation calculation towards pipe side subregion heat-transfer pipe, that is, inner wall Arrive heat.
(4) thermal parameter calculating is carried out, specific implementation method is as follows:
A. it is directed to shell-side subregion, in addition to hot well, each subregion steam pressure is solved using quality of steam equilibrium equation;For heat Well subregion is calculated liquid level using water quality equilibrium equation, is asked using quality of steam equilibrium equation and water quality equilibrium equation simultaneous Solve hot well zoned pressure;
B. it is directed to each subregion in pipe side, the circulating water pressure of each subregion is calculated using water quality equilibrium equation;
C. it is directed to shell-side subregion, in addition to hot well, using saturation temperature as zone temperature;Hot well vapor (steam) temperature is saturation temperature Degree, coolant-temperature gage are condensate mixing temperature;
D. it is directed to pipe side subregion, circulating water temperature is determined using pressure and enthalpy.
(5) data export.According to form and type that data need, organized data output program exports calculated result.
In above-mentioned steps, the physical parameter in relation to water and vapor is calculated according to IAPWS-IF97 formula.
In step (4), each zoned pressure is calculated using mass balance method;In hot well subregion, water quality is calculated first Measure equilibrium equation:
In formula, Mw--- the water quality in hot well, kg;
τ --- time, s;
ρw--- the water density in hot well, kg/m3
Vw--- the water volume in hot well, m3
Ps--- hot well pressure, MPa;
∑fw,in--- into the sum of the hydrophobic flow of hot well, kg/s;
∑fw,out--- the sum of the condensate flow of outflow hot well, kg/s;
--- the sum of each subregion steam condensation rate, kg/s.
Water volume can be calculated according to above formula, and then water level can be calculated according to hot well size.The steam matter of hot well subregion Measure equilibrium equation:
In formula, Ms--- the quality of steam in hot well, kg;
ρs--- the vapour density in hot well, kg/m3
Vs--- the water volume in hot well, m3
∑fs,in--- into the sum of the steam flow of hot well, kg/s;
∑fs,out--- the sum of the steam flow of outflow hot well, kg/s.
It is as follows that pressure can be calculated according to above-mentioned two equation:
It is illustrated in conjunction with non-homogeneous subregion emulation mode of the Fig. 2 to conventional condenser.It is its structural representation on the left of Fig. 2 Figure, E is air intake, F is air ejector.Right side is its work flow diagram, and wherein A corresponds to turbine discharge, B corresponds to the matter of by stages Amount transmitting, C correspond to mixed gas, the corresponding condensate for falling into hot well of D that air ejector is taken away.Its emulation mode is as follows:
(1) subregion divides: condenser being divided into 3 subregions by shell radial direction in condenser shell-side, on the right side of Fig. 2 Shown in schematic diagram.Subregion 11,12 includes the heat-transfer pipe of different number, and it further includes pipe side-entrance water chamber that subregion 13, which is hot well subregion, With outlet water chamber.Subregion 11 corresponds to the corresponding air ejector of turbine discharge A, subregion 12 and is evacuated C.Except hot well subregion, pipe side-entrance water Outside room and outlet water chamber, shell-side subregion and pipe side subregion are corresponded;
(2) input data: full mold structured data input include with the International System of Units description nuclear power unit include will be cold Outer dimension, internal diameter of the pipeline, outer diameter tube, duct length, the flow resistance data of condenser;The input of integer structured data includes core Power device condenser tube number of rows, process number, heat-transfer pipe arrangement mode;Wherein, think when heat-transfer pipe arrangement mode is triangle Value is 1, is 0 when square;The input of running boundary data includes pipe side circulating water intake flow, temperature, atmospheric pressure, air ejector Suction pressure, pressure before condensate pump, turbine discharge flow, pressure, mass dryness fraction;Calculating initial value input includes the calculating with just Value, specially pipe side and each by stages flow of shell-side, initial pressure, initial temperature, rate of setting.
(3) heat exchange calculates: the heat exchange of subregion 11,12 shell-sides is calculated using modified Nu Saier film condensation calculation relational expression Amount and rate of setting;The heat exchange amount of subregion 11,12 pipe sides is calculated using revised forced-convection heat transfer D-B formula;Utilize gold Belong to thermally conductive formula and calculates subregion 11, metal heat-conducting;
(4) thermal parameter calculates: using boundary parameter and the mass transfer B of by stages, passing through quality of steam EQUILIBRIUM CALCULATION FOR PROCESS The pressure of subregion 11,12 shell-sides;Pass through water quality EQUILIBRIUM CALCULATION FOR PROCESS subregion 11,12 pipe sides and pipe side-entrance water chamber and outlet water chamber Interior circulating water pressure;Pass through hot well water quality balance, the liquid level for the condensate D calculating subregion 13 for falling into hot well;It is steamed by hot well Vapour mass balance and water quality balance simultaneous calculate 13 pressure of subregion;
(5) result exports: the type and format of data as needed export calculated result.
(6) after circulation primary, a time step is pushed ahead, needs to set cycling condition according to the real-time calculating time, Complete simulation work.
It is illustrated in conjunction with uniform segmentation emulation mode of the Fig. 3 to common condenser.It is its structural schematic diagram, H on the left of Fig. 3 For subsidiary engine air intake, G is host air intake, and M is air ejector;Right side is its workflow schematic diagram, wherein emulation mode are as follows: I For subsidiary engine into vapour, J be host into vapour, L be the mixed gas taken away of air ejector, K is hot well.Its emulation mode is as follows:
(1) subregion divides: being evenly dividing common condenser in addition to hot well by radial direction in common condenser shell-side is 3 × 3 subregions, as shown in Fig. 3 right side schematic view.It include the heat-transfer pipe of identical quantity in subregion 1- subregion 9, subregion 10 is hot well Subregion K further includes recirculated water import and export water chamber.Subregion 1 corresponds to subsidiary engine steam discharge I, 3 corresponding host steam discharge J of subregion, and subregion 7 is corresponding Air ejector is evacuated L.In addition to subregion 10, recirculated water import and export water chamber, shell-side subregion and pipe side subregion are corresponded;
(2) input data: full mold structured data input include with the International System of Units description nuclear power unit include will be cold Outer dimension, internal diameter of the pipeline, outer diameter tube, duct length, the flow resistance data of condenser;The input of integer structured data includes core Power device condenser tube number of rows, process number, heat-transfer pipe arrangement mode;Wherein, think when heat-transfer pipe arrangement mode is triangle Value is 1, is 0 when square;The input of running boundary data includes pipe side circulating water intake flow, temperature, atmospheric pressure, air ejector Suction pressure, pressure before condensate pump, host exhaust steam flow, pressure, mass dryness fraction, subsidiary engine exhaust steam flow, pressure, mass dryness fraction;Calculate initial value Input includes that described calculate uses initial value, specially pipe side and each by stages flow of shell-side, initial pressure, initial temperature, condensation speed Rate.
(3) heat exchange calculates: calculating changing for 9 shell-side of subregion 1- subregion using modified Nu Saier film condensation calculation relational expression Heat and rate of setting;The heat exchange amount of 9 pipe side of subregion 1- subregion is calculated using revised forced-convection heat transfer D-B formula;Benefit The metal heat-conducting of subregion 1- subregion 9 is calculated with metal heat-conducting formula;
(4) thermal parameter calculates: passing through the pressure of 9 shell-side of quality of steam EQUILIBRIUM CALCULATION FOR PROCESS subregion 1- subregion;Pass through water quality 9 pipe side of EQUILIBRIUM CALCULATION FOR PROCESS subregion 1- subregion and the indoor circulating water pressure of recirculated water import and export water;It is balanced by hot well water quality Calculate the liquid level of subregion 10;Simultaneous, which is balanced, with water quality by hot well quality of steam balance calculates 10 pressure of subregion;
(5) result exports: the type and format of data as needed export calculated result.
(6) after circulation primary, a time step is pushed ahead, needs to set cycling condition according to the real-time calculating time, Complete simulation work.

Claims (3)

1. a kind of partitioning emulation mode of nuclear power unit condenser, it is characterized in that:
(1) subregion division is carried out to simulation object, the simulation object includes conventional condenser and common condenser, the subregion Divide includes uniform segmentation and non-homogeneous subregion;It specifically includes:
A. the non-homogeneous subregion of conventional condenser is divided: subregion is carried out according to the position of practical condenser sensor and is drawn Point, the subregion of division includes subregion (11) turbine discharge mouth region domain, subregion (12) air ejector region, subregion (13) hot well area Domain, sensor region, other shell-side regions in addition to above-mentioned zone, pipe side circulating water intake water chamber, pipe side recirculated water go out Saliva room;
B. the uniform segmentation of conventional condenser is divided: condenser is divided into 2 × 2,3 × 3 or 4 × 4 in addition to hot well Subregion, each subregion volume is identical, and includes the heat transfer pipe of identical quantity;
C. the non-homogeneous subregion of common condenser is divided: according to the position of practical condenser inner sensor and host, The position of subsidiary engine air intake carries out subregion division;
D. common condenser uniform segmentation is divided: common condenser is divided into 3 × 3 or 4 × 4 subregions in addition to hot well, Each subregion volume is identical, and includes the heat transfer pipe of identical quantity;
(2) it inputs the structure size of condenser, physical property, boundary parameter, calculate initial value;
(3) Calculation of Heat Transfer is carried out to the one-to-one subregion of shell-side and pipe side;It specifically includes:
A. the input of full mold structured data is started in calculating input included data with integer structured data;In iteration mistake It is additional to input included data using calculating initial value when first time in journey calculates;
B. Nu Saier film condensation calculation formula is utilized, the amendment of adding tube number of rows, the amendment of process number calculate in addition to hot well subregion Each shell-side subregion condensation heat transfer and steam condensation rate;
C. utilize forced-convection heat transfer D-B formula, addition varied property amendment, calculate in the subregion of each pipe side forced-convection heat transfer amount with Recirculated water temperature rise;
D. it is arrived by each shell-side subregion heat-transfer pipe outer wall of metal heat-conducting equation calculation towards pipe side subregion heat-transfer pipe, that is, inner wall Heat;
(4) thermal parameter calculating is carried out;It specifically includes:
A. it is directed to shell-side subregion, in addition to hot well, each subregion steam pressure is solved using quality of steam equilibrium equation;For hot well point Area calculates liquid level using water quality equilibrium equation, using quality of steam equilibrium equation and water quality equilibrium equation simultaneous solution heat Well zoned pressure;
B. it is directed to each subregion in pipe side, the circulating water pressure of each subregion is calculated using water quality equilibrium equation;
C. it is directed to shell-side subregion, in addition to hot well, using saturation temperature as zone temperature;Hot well vapor (steam) temperature is saturation temperature, Coolant-temperature gage is condensate mixing temperature;
D. it is directed to pipe side subregion, circulating water temperature is determined using pressure and enthalpy;
(5) simulation result exports.
2. the partitioning emulation mode of nuclear power unit condenser according to claim 1, it is characterized in that described pair of emulation Object carries out in subregion division further include:
A. in the non-homogeneous subregion division of conventional condenser: except hot well subregion (13), pipe side circulating water intake water chamber, pipe side are followed Outside ring water out water chamber, other each shell-sides are overlapped with pipe side subregion physical region, are corresponded, should all be included portion of heat transfer pipe Road;
B. in the uniform segmentation division of conventional condenser: by turbine discharge mouth region domain, air ejector region, sensor region point It is not divided into wherein some subregion;Hot well part is individually divided into a subregion, pipe side circulating water intake water chamber is individually divided into one A subregion, pipe side circulating water outlet water chamber are individually divided into a subregion;Except hot well subregion, pipe side circulating water intake water chamber, pipe Outside the circulating water outlet water chamber of side, other each shell-sides are overlapped with pipe side subregion physical region, are corresponded;
C. in dividing to the non-homogeneous subregion of common condenser: the subregion of division include host steam discharge region, subsidiary engine steam discharge region, Air ejector region, hot well subregion, sensor region, other shell-side regions in addition to above-mentioned zone, pipe side circulating water intake Water chamber, pipe side circulating water outlet water chamber;Except hot well subregion, pipe side circulating water intake water chamber subregion, pipe side circulating water outlet water chamber Outside outer subregion, other each shell-sides are overlapped with pipe side subregion physical region, are corresponded, should all be included portion of heat transfer pipeline;
D. to common condenser uniform segmentation divide in: by host steam drain region, subsidiary engine steam drain region, air ejector region, Sensor region is divided into wherein some subregion respectively;Subregion divides number with host steam drain region and subsidiary engine steam drain region not Subject to same subregion;Hot well part is individually divided into a subregion, pipe side circulating water intake water chamber is individually divided into one point Area, pipe side circulating water outlet water chamber are individually divided into a subregion;Except hot well subregion, pipe side circulating water intake water chamber subregion, pipe Outside the circulating water outlet water chamber subregion of side, other each shell-sides are overlapped with pipe side subregion physical region, are corresponded.
3. the partitioning emulation mode of nuclear power unit condenser according to claim 1, it is characterized in that the input is cold The structure size of condenser, boundary parameter, calculates initial value at physical property, specifically includes:
A. the input of full mold structured data includes in the outer dimension of the nuclear power unit condenser of International System of Units description, pipeline Diameter, outer diameter tube, duct length, flow resistance data;
B. the input of integer structured data includes nuclear power unit condenser tube number of rows, process number, heat-transfer pipe arrangement mode;Wherein, Heat-transfer pipe arrangement mode thinks that value is 1 when being triangle, be 0 when square;
C. the input of running boundary data includes pipe side circulating water intake flow, temperature, atmospheric pressure, and air ejector suction pressure coagulates Pressure before water pump;It further include turbine discharge flow, pressure, mass dryness fraction for conventional condenser;For common condenser, also wrap Include host exhaust steam flow, pressure, mass dryness fraction and subsidiary engine exhaust steam flow, pressure, mass dryness fraction;
D. calculating initial value input includes that described calculate uses initial value, specially pipe side and each by stages flow of shell-side, initial pressure, just Beginning temperature, rate of setting.
CN201610953179.1A 2016-11-03 2016-11-03 A kind of partitioning emulation mode of nuclear power unit condenser Active CN106528985B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610953179.1A CN106528985B (en) 2016-11-03 2016-11-03 A kind of partitioning emulation mode of nuclear power unit condenser

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610953179.1A CN106528985B (en) 2016-11-03 2016-11-03 A kind of partitioning emulation mode of nuclear power unit condenser

Publications (2)

Publication Number Publication Date
CN106528985A CN106528985A (en) 2017-03-22
CN106528985B true CN106528985B (en) 2019-06-25

Family

ID=58326906

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610953179.1A Active CN106528985B (en) 2016-11-03 2016-11-03 A kind of partitioning emulation mode of nuclear power unit condenser

Country Status (1)

Country Link
CN (1) CN106528985B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107784156B (en) * 2017-09-07 2020-11-03 中国船舶重工集团公司第七一九研究所 Method for calculating parameters of steam discharge system of nuclear power plant
CN108763748B (en) * 2018-05-28 2019-04-09 西安交通大学 A method of for heat pipe heap Transient Analysis
CN110059388B (en) * 2019-04-10 2023-01-03 哈尔滨工程大学 Parameter distribution modeling simulation method for condenser of nuclear power plant along working medium flowing direction
CN110705079B (en) * 2019-09-25 2023-01-03 哈尔滨工程大学 Centrifugal compressor structure optimization method based on simulated annealing algorithm
CN111897238B (en) * 2020-08-03 2024-01-26 中国舰船研究设计中心 Semi-physical control loop simulation test device of large power system
CN111897240B (en) * 2020-08-05 2023-03-31 哈尔滨工程大学 Simulation method and system based on nuclear power system operation
CN113821871B (en) * 2021-08-31 2023-02-17 中国船舶重工集团公司第七0三研究所 Dynamic simulation model development method for multi-structure and multi-pressure-node heat exchanger of ship two-loop system

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103017821A (en) * 2012-11-27 2013-04-03 芜湖博耐尔汽车电气系统有限公司 Novel condenser simulating computation experiment correction method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103017821A (en) * 2012-11-27 2013-04-03 芜湖博耐尔汽车电气系统有限公司 Novel condenser simulating computation experiment correction method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Development and analysis of a detailed parametric simulation model of condensers for nuclear power plants;ZHU Hai-shan 等;《Nuclear Safety and Simulation》;20141231;第5卷(第4期);第347-357页
共用凝汽器建模与仿真;薛若军 等;《核动力工程》;20140430;第35卷(第1期);第21-25页
核电站凝汽器仿真模型研究;吴鹏 等;《原子能科学技术》;20140131;第48卷(第1期);第92-98页

Also Published As

Publication number Publication date
CN106528985A (en) 2017-03-22

Similar Documents

Publication Publication Date Title
CN106528985B (en) A kind of partitioning emulation mode of nuclear power unit condenser
CN104049539B (en) A kind of full working scope analogue system of Steam Turbine Regenerative System
CN102779217B (en) Computer simulation performance computation method of refrigeration system under frosting working condition
CN104088771B (en) The accurate determination method of circulating cooling water of power plant system water pump assembly optimum combination operating scheme
CN109670216A (en) Passive residual heat removal condenser locations optimum design method based on CFD technology
CN104834773A (en) Simulation method for heat exchange performance of straight tube type once-through steam generator
CN110096845A (en) A kind of mixed flow type closed cooling tower design and calculation method
CN107632042B (en) Single-phase heat exchange experiment test platform and test method for shell-and-tube heat exchanger
CN105303037A (en) Method for calculating critical air leakage amount causing condenser vacuum deterioration
CN110309591A (en) It exchanges heat under a kind of flat finned heat exchanger air side laminar condition and drag computation method
CN110059388A (en) Parameter distribution modeling and simulating method of nuclear power plant's condenser along Working fluid flow direction
CN109932649A (en) A kind of fired power generating unit generating efficiency monitoring method
CN110186291A (en) A kind of mixed flow type closed cooling tower calculation and check method
CN111159875B (en) Dynamic coupling mathematical model and modeling method for shell side and tube side of power station condenser
CN105160116B (en) The uniform Optimization Design of heat exchanger liquid separation
CN109522644A (en) A kind of enhanced heat exchange surface comprehensive performance evaluation method
CN108595813B (en) A kind of CO2 gas cooler optimum design method based on optimum quality flow velocity
CN110705075A (en) Optimization design method of static pressure box hole plate
CN109344423A (en) A kind of calculation method for closing the practical IP efficiency of cylinder steam turbine
CN206132397U (en) Turbine test device that bloies
CN205909664U (en) Energy -efficient desiccator of developments
CN107664653A (en) Heat exchanger condensation heat transfer experiment test platform and method of testing
CN113486608A (en) Two-loop modular modeling simulation method for nuclear power system
CN207439695U (en) Heat exchanger boiling heat transfer efficiency test platform
CN207439979U (en) Heat exchanger condensation heat transfer experiment test platform

Legal Events

Date Code Title Description
C06 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