One of the most important aims of this study is to improve the core of the current VVER reactors to achieve more burn-up(or more cycle length)and more intrinsic safety.It is an independent study on the Russian new pro...One of the most important aims of this study is to improve the core of the current VVER reactors to achieve more burn-up(or more cycle length)and more intrinsic safety.It is an independent study on the Russian new proposed FAs,called TVS-2M,which would be applied for the future advanced VVERs.Some important aspects of neutronics as well as thermal hydraulics investigations(and analysis)of the new type of Fas are conducted,and results are compared with the standards PWR CDBL.The TVS-2M FA contains gadolinium-oxide which is mixed with UO_(2)(for different Gd densities and U-235 enrichments which are given herein),but the core does not contain BARs.The new type TVS-2M Fas are modeled by the SARCS software package to find the PMAXS format for three states of CZP and HZP as well as HFP,and then the whole core is simulated by the PARCS code to investigate transient conditions.In addition,the WIMS-D5 code is suggested for steady core modeling including TVS-2M FAs and/or TVS FAs.Many neutronics aspects such as the first cycle length(first cycle burn up in terms of MWthd/kgU),the critical concentration of boric acid at the BOC as well as the cycle length,the axial,and radial power peaking factors,differential and integral worthy of the most reactive CPS-CRs,reactivity coefficients of the fuel,moderator,boric acid,and the under-moderation estimation of the core are conducted and benchmarked with the PWR CDBL.Specifically,the burn-up calculations indicate that the 45.6 d increase of the first cycle length(which corresponds to 1.18 MWthd/kgU increase of burnup)is the best improving aim of the new FA type called TVS-2M.Moreover,thermal-hydraulics core design criteria such as MDNBR(based on W3 correlation)and the maximum of fuel and clad temperatures(radially and axially),are investigated,and discussed based on the CDBL.展开更多
A theoretical “drift-flux based thermal-hydraulic mixture-fluid coolant channel model” is presented. It is the basis to a corresponding digital “Coolant Channel Module (CCM)”. This purpose derived “Separate-Regio...A theoretical “drift-flux based thermal-hydraulic mixture-fluid coolant channel model” is presented. It is the basis to a corresponding digital “Coolant Channel Module (CCM)”. This purpose derived “Separate-Region Mixture Fluid Approach” should yield an alternative platform to the currently dominant “Separate-Phase Models” where each phase is treated separately. Contrary to it, a direct procedure could be established with the objective to simulate in an as general as possible way the steady state and transient behaviour of characteristic parameters of single- and/or (now non-separated) two-phase fluids flowing within any type of heated or non-heated coolant channels. Their validity could be confirmed by a wide range of verification and validation runs, showing very satisfactory results. The resulting universally applicable code package CCM should provide a fundamental element for the simulation of thermal-hydraulic situations over a wide range of complex systems (such as different types of heat exchangers and steam generators as being applied in both conventional but also nuclear power stations, 1D and 3D nuclear reactor cores etc). Thereby the derived set of equations for different coolant channels (distinguished by their key numbers) as appearing in these systems can be combined with other ODE-s and non-linear algebraic relations from additional parts of such an overall model. And these can then to be solved by applying an appropriate integration routine. Within the solution procedure, however, mathematical discontinuities can arise. This due to the fact that along such a coolant channel transitions from single- to two-phase flow regimes and vice versa could take place. To circumvent these difficulties it will in the presented approach be proposed that the basic coolant channel (BC) is subdivided into a number of sub-channels (SC-s), each of them being occupied exclusively by only a single or a two-phase flow regime. After an appropriate nodalization of the BC (and thus its SC-s) and after applying a “modified finite volume method” together with other special activities the fundamental set of non-linear thermal-hydraulic partial differential equations together with corresponding constitutive relations can be solved for each SC separately. As a result of such a spatial discretization for each SC type (and thus the entire BC) the wanted set of non-linear ordinary differential equations of 1st order could be established. Obviously, special attention had to be given to the varying SC entrance or outlet positions, describing the movement of boiling boundaries or mixture levels along the channel. Including even the possibility of SC-s to disappear or be created anew during a transient.展开更多
In the core of a molten salt fast reactor(MSFR),heavy metal fuel and fission products can be dissolved in a molten fluoride salt to form a eutectic mixture that acts as both fuel and coolant.Fission energy is released...In the core of a molten salt fast reactor(MSFR),heavy metal fuel and fission products can be dissolved in a molten fluoride salt to form a eutectic mixture that acts as both fuel and coolant.Fission energy is released from the fuel salt and transferred to the second loop by fuel salt circulation.Therefore,the MSFR is characterized by strong interaction between the neutronics and the thermal hydraulics.Moreover,recirculation flow occurs,and nuclear heat is accumulated near the fertile blanket,which significantly affects both the flow and the temperature fields in the core.In this work,to further optimize the conceptual geometric design of the MSFR,three geometries of the core and fertile blanket are proposed,and the thermal-hydraulic characteristics,including the three-dimensional flow and temperature fields of the fuel and fertile salts,are simulated and analyzed using a coupling scheme between the open source codes OpenMC and OpenFOAM.The numerical results indicate that a flatter core temperature distribution can be obtained and the hot spot and flow stagnation zones that appear in the upper and lower parts of the core center near the reflector can be eliminated by curving both the top and bottom walls of the core.Moreover,eight cooling loops with a total flow rate of0.0555 m3 s-1 ensur an acceptable temperature distribusure an acceptable temperature distribution in the fertile blanket.展开更多
In molten salt reactors(MSRs),the liquid fuel salt circulates through the primary loop and a part of the delayed neutron precursors(DNPs)decays outside the reactor core.To model and analyze the flow field effect of DN...In molten salt reactors(MSRs),the liquid fuel salt circulates through the primary loop and a part of the delayed neutron precursors(DNPs)decays outside the reactor core.To model and analyze the flow field effect of DNPs in channel-type liquid-fueled MSRs,a three-dimensional space-time dynamics code,named ThorCORE3D,that couples neutronics,core thermalhydraulics,and a molten salt loop system was developed and validated with the Molten Salt Reactor Experiment(MSRE)benchmarks.The effects of external loop recirculation time,fuel flow rate,and core flow field distribution on the delayed neutron fraction loss of MSRE at steadystate were modeled and simulated using the ThorCORE3D code.Then,the flow field effect of the DNPs on the system responses of the MSRE in the reactivity insertion transient under different initial conditions was analyzed systematically for the channel-type liquid-fueled MSRs.The results indicate that the flow field condition has a significant effect on the steady-state delayed neutron fractions and will further affect the transient power and temperature responses of the reactor system.The analysis results for the effect of the DNP flow field can provide important references for the design optimization and safety analysis of liquid-fueled MSRs.展开更多
A molten salt reactor(MSR) is one of the six advanced reactor concepts selected by the generation Ⅳ international forum because of its advantages of inherent safety, and the promising capabilities of Th-U breeding an...A molten salt reactor(MSR) is one of the six advanced reactor concepts selected by the generation Ⅳ international forum because of its advantages of inherent safety, and the promising capabilities of Th-U breeding and transuranics transmutation. A dynamics model for the channel-type MSR is developed in this work based on a three-dimensional thermal–hydraulic model(3DTH) and a point reactor model. The 3DTH couples a three-dimensional heat conduction model and a one-dimensional single-phase flow model that can accurately consider the heat conduction between different assemblies. The 3DTH is validated by the RELAP5 code in terms of the temperature and mass flow distribution calculation. A point reactor model considering the drift of delayed neutron precursors is adopted in the dynamics model. To verify the dynamics model, three experiments from the molten salt reactor experiment are simulated. The agreement of the experimental data and simulation results was excellent.With the aid of this model, the unprotected step reactivity addition and unprotected loss of flow of the 2 MWt experimental MSR are modeled, and the reactor power and temperature evolution are analyzed.展开更多
Molten salt-cooled reactor is one of the six GenIV reactors with promising characteristics including safety,reliability,proliferation resistance,physical protection,economics,and sustainability.In this paper,a small i...Molten salt-cooled reactor is one of the six GenIV reactors with promising characteristics including safety,reliability,proliferation resistance,physical protection,economics,and sustainability.In this paper,a small innovative molten chloride-cooled fast reactor(MCCFR)with 30-year core and a target 120-MWt thermal power was presented.For its feasible study,neutronics,thermal-hydraulics,and radiation damage analysis were performed.The key design properties including kinetics parameters,reactivity swing,reactivity feedback coefficients,maximum accumulated displacement per atom(DPA)of reactor pressure vessel(RPV)and fuel cladding,and maximum coolant,cladding,and fuel temperatures were evaluated.The results showed the MCCFR could operate without refueling for 30 years with overall negative reactivity feedback coefficients up the end of its life.During its 30-year life,the excess reactivity was well managed by constantly pulling out the control rods.The maximum accumulated DPA on RPV and fuel cladding were 8.92 dpa and 197.03 dpa,respectively,which are both below the design limits.Similarly,the maximum coolant,cladding and fuel center temperatures were all below the design limits during its entire lifetime.According to these results,the MCCFR core design with long life is feasible.展开更多
To perform nuclear reactor simulations in a more realistic manner,the coupling scheme between neutronics and thermal-hydraulics was implemented in the HNET program for both steady-state and transient conditions.For si...To perform nuclear reactor simulations in a more realistic manner,the coupling scheme between neutronics and thermal-hydraulics was implemented in the HNET program for both steady-state and transient conditions.For simplicity,efficiency,and robustness,the matrixfree Newton/Krylov(MFNK)method was applied to the steady-state coupling calculation.In addition,the optimal perturbation size was adopted to further improve the convergence behavior of the MFNK.For the transient coupling simulation,the operator splitting method with a staggered time mesh was utilized to balance the computational cost and accuracy.Finally,VERA Problem 6 with power and boron perturbation and the NEACRP transient benchmark were simulated for analysis.The numerical results show that the MFNK method can outperform Picard iteration in terms of both efficiency and robustness for a wide range of problems.Furthermore,the reasonable agreement between the simulation results and the reference results for the NEACRP transient benchmark verifies the capability of predicting the behavior of the nuclear reactor.展开更多
Previous analytical results on flow splitting are generalized to consider multiple boiling channels systems. The analysis is consistent with the approximations usually adopted in the use of systems codes (like RELAP5 ...Previous analytical results on flow splitting are generalized to consider multiple boiling channels systems. The analysis is consistent with the approximations usually adopted in the use of systems codes (like RELAP5 and TRACE5, among others) commonly applied to perform safety analyses of nuclear power plants. The problem is related to multiple, identical, parallel boiling channels, connected through common plena. A theoretical model limited in scope explains this flow splitting without reversal. The unified analysis performed and the confirmatory computational results found are summarized in this paper. New maps showing the zones where this behavior is predicted are also shown considering again twin pipes. Multiple pipe systems have been found not easily amenable for analytical analysis when dealing with more than four parallel pipes. However, the particular splitting found (flow along N pipes dividing in one standalone pipe flow plus N -1 identical pipe flows) has been verified up to fourteen pipes, involving calculations in systems with even and odd number of pipes using the RELAP5 systems thermal-hydraulics code.展开更多
In this study,COMSOL v5.2 Multiphysics software was utilized to perform coupled neutronics and thermal–hydraulics simulations of a molten salt fast reactor,and the SCALE v6.1 code package was utilized to generate the...In this study,COMSOL v5.2 Multiphysics software was utilized to perform coupled neutronics and thermal–hydraulics simulations of a molten salt fast reactor,and the SCALE v6.1 code package was utilized to generate the homogenized cross-section data library.The library’s 238 cross-section groups were categorized into nine groups for the simulations in this study.The results of the COMSOL model under no fuel flow conditions were verified using the SCALE v6.1 code results,and the results of the neutronics and thermal–hydraulics simulations were compared to the results of previously published studies.The results indicated that the COMSOL model that includes the cross-section library generated by the SCALE v6.1 code package is suitable for the steady-state analysis and design assessment of molten salt fast reactors.Subsequently,this model was utilized to investigate the neutronics and thermal–hydraulics behaviors of the reactor.Multiple designs were simulated and analyzed in this model,and the results indicated that even if the wall of the core is curved,hot spots occur in the upper and lower portions of the core’s center near the reflectors.A new design was proposed that utilizes a flow rate distribution system,and the simulation results of this design showed that the maximum temperature in the core was approximately 1032 K and no hot spots occurred.展开更多
In this paper,a novel composite heat transfer enhancement technique comprised of louvered fins(LFs)and rectangular wing vortex generators(RWVGs)is proposed to improve the LF side thermal-hydraulic performance of louve...In this paper,a novel composite heat transfer enhancement technique comprised of louvered fins(LFs)and rectangular wing vortex generators(RWVGs)is proposed to improve the LF side thermal-hydraulic performance of louvered fin and flat tube heat exchangers(LFHEs).After validation of the LF side pressure dropΔP and heat transfer coefficient hLFof the baseline by experiments,the numerical method is applied to investigate the influential mechanisms of the RWVG parameters(the number N(7 to 15),attack angleβ(30°to 90°),height H_(VG)(0.8 mm to 2 mm)and width W_(VG)(0.8 mm to 1.2 mm))on the performance of the LFHE in the velocity range of 3 m/s to 10 m/s.Results show that thermal-hydraulic performance of the LFHE is significantly impacted by the RWVGs,and according to the performance evaluation criteria(PEC),the LFHE achieves its optimal thermal-hydraulic performance when N=7,β=45°,H_(VG)=1.8 mm and W_(VG)=1 mm.Compared to the baseline,the maximum,minimum and average increments of PEC for the optimal case are 13.85%,4.67%and 8.39%,respectively.展开更多
Following the previous relevant works [1]-[3], a scaling analysis is performed to derive a set of scaling criteria which were thought to be suitable for reproducing the major thermal-hydraulic phenomena in a scaled-do...Following the previous relevant works [1]-[3], a scaling analysis is performed to derive a set of scaling criteria which were thought to be suitable for reproducing the major thermal-hydraulic phenomena in a scaled-down CANDU moderator tank similar to that in a prototype power plant during a full power steady state condition. The objective of building a scaled-down moderator tank is to generate the experimental data necessary to validate the computer codes which are used to analyze the accident analysis of CANDU-6 plants. The major variables of interests in this paper are moderator flow velocity and temperature of the moderator which is D2O inside the moderator tank during a steady state and transient conditions. The reason is that the local subcooling of the moderator is found to be a critical parameter determining whether the stable film boiling can sustain on the outer surface of the calandria tube if the contact of overheated pressure tube and cold calandria tube should occur due to pressure tube ballooning during LBLOCA with ECC injection failure [4]. The key phenomena involved include the inlet jet development and impingement, buoyancy force driven by the moderator temperature gradient caused by non-uniform direct heating of the moderator, and the pressure drop due to viscous friction of the flow across the calandria tube array. In this paper, the previous researches are reviewed, some concerns or potential problems associated with them implied by comparing CFD analyses results between the CANDU-6 moderator tank and 1/4 scaled-down test facility are described, and as a way to examine the assumption of the scaling analysis is true an order-of-magnitude analyses are performed. Based on the results of these analyses the assumption of neglecting ?and ?terms cannot be justified for the power of 0.5 MW and 1.566 MW for the 1/4 scaled-down facility. Further investigation is thought to be necessary to confirm this result, i.e. if the scaling of the previous work1 is justifiable by some other independent analyses.展开更多
Development and use of nuclear energy is currently growing very rapidly, in order to achieve increasingly advanced technology, both in terms of design, economic factors and safety factors. Thermal-hydraulics aspects o...Development and use of nuclear energy is currently growing very rapidly, in order to achieve increasingly advanced technology, both in terms of design, economic factors and safety factors. Thermal-hydraulics aspects of nuclear reactors should be done with calculation and near-perfect condition. Including today began development of a nuclear reactor with low power below 300 MW, or commonly called the Small Modular Reactor (SMR). One is CAREM-25 developed by Argentina with a power of 25 MW, where in CAREM already using natural circulation system and the use of nanofluid as coolant fluid. In this research, analytic modeling of thermal-hydraulics nuclear reactor SMR CAREM-25, when the nanofluid Al<sub>2</sub>O<sub>3</sub>-Water used as cooling fluid in the cooling system of a nuclear reactor. Further to this analytic modeling will be done on CFD. Analytic modeling with CFD to determine the flow phenomena and distribution as well as the effect of nano-particles of Al<sub>2</sub>O<sub>3</sub>-Water based on the volume fraction (1% and 3%) of the coefficient of heat transfer by natural convection.展开更多
This study investigated the effects of zigzag-flow channel bending angle in printed circuit heat exchangers(PCHEs) using a computational fluid dynamics method with ANSYS-FLUENT simulation.The three-dimensional model o...This study investigated the effects of zigzag-flow channel bending angle in printed circuit heat exchangers(PCHEs) using a computational fluid dynamics method with ANSYS-FLUENT simulation.The three-dimensional model of PCHE with a 15° curved,zigzag channel was conducted for preliminary validation.The comparisons between the CFD simulation results and the experimental data showed good agreement with some discrepancies in the heat transfer and pressure drop results.In addition,different bending angle configurations(0°,3° to 30°) of zigzag channels were analyzed to obtain better thermal-hydraulic performance of the zigzag channel PCHE under different inlet mass flow rates.The criteria of heat transfer and frictional factor were applied to evaluate the thermal-hydraulic performance of the PCHE.The results showed that the 6° and 9°bending channel provided good thermal-hydraulic performance.New correlations were developed using the 6°and 9° bending channel angles in PCHE designs to predict the Nusselt number and friction factor.展开更多
A double reheat ultra-supercritical boiler is an important development direction for high-parameter and large-capacity coal-fired power plants due to its high thermal efficiency and environmental value.China has devel...A double reheat ultra-supercritical boiler is an important development direction for high-parameter and large-capacity coal-fired power plants due to its high thermal efficiency and environmental value.China has developed a 1000 MW double reheat ultra-supercritical boiler with steam parameters of 35 MPa at 605℃/613℃/613℃.Reasonable water wall design is one of the keys to safe and reliable operation of the boiler.In order to examine the thermal-hydraulic characteristics of the double reheat ultra-supercritical boiler,the water wall system was equivalent to a flow network comprising series-parallel circuits,linking circuits and pressure nodes,and a calculation model was built on account of the conservation equations of energy,momentum and mass.Through the iterative solving of nonlinear equations,the prediction parameters of the water wall at boiler maximum continue rate(BMCR),75%turbine heat-acceptance rate(THA)and 30%THA loads,including total pressure drops,flow distribution,outlet steam temperatures,fluid and metal temperatures were gotten.The results of calculation exhibit excellent thermal-hydraulic characteristics and substantiate the feasibility of the water wall design of the double reheat ultra-supercritical boiler.展开更多
In terms of developing supercritical CO_(2)(sCO_(2))coal-fired power plants,enhancing cooling wall performance is one of significant factors to improve system performance.In this paper,a new cooling wall tube structur...In terms of developing supercritical CO_(2)(sCO_(2))coal-fired power plants,enhancing cooling wall performance is one of significant factors to improve system performance.In this paper,a new cooling wall tube structure is proposed to match the non-uniform heat flux(NUH)with the thermal resistance by changing the cooling wall tube eccentricity.A three-dimensional multi-physical coupling model of cooling wall is constructed to compare the novel structure to the conventional structures.The properties of fluid dynamics,thermal stress,coupled heat transfer and cooling wall deformation are analyzed.In contrast to the traditional structure,the maximum temperature and circumferential temperature difference(CTD)of the proposed structure can be reduced by 2%and 27.4%,respectively.The essential working parameters related to the performances of the cooling wall tube are discussed.The maximum temperature of the new structure is reduced by 8-13 K and the maximum thermal stress is reduced by about 10%-15%under all the simulated working conditions when the eccentricity changes from 0 to 0.2.The proposed structure can effectively reduce the maximum temperature and circumferential temperature gradient under NUH.Consequently,a novel insight is put out for the design and optimization of the cooling wall tube in coal-fired power plants.展开更多
A combustion model of a large-scale supercritical circulati ng fluidized bed (CFB) boiler was developed for comprehensive computational-fluid-dynamics analysis. The model incorporates gas-solid hydrodynamics, coal com...A combustion model of a large-scale supercritical circulati ng fluidized bed (CFB) boiler was developed for comprehensive computational-fluid-dynamics analysis. The model incorporates gas-solid hydrodynamics, coal combustion, heat transfer on heat exchange surfaces in the furnace, and heat transfer between fumace and working medium in the heat transfer tubes. In simulating the dense and dilute phases in the fumace, the gas-solid hydrodynamics is based on the Euler-Euler model and energy-minimization multiscale drag model. Coal combustion entails evaporation, devolatilization, char combustion, gas homoge neous reaction, and pollutant emission. The coefficient ofheat transfer between gas-solid and the waterwall is estimated using the cluster renewal model, and for radiation, the discrete ordinate model is used. Moreover, thermohydraulic processes in the membrane wall are also in eluded in the heat transfer process. The model was successfully applied in simulations of a 350-MW supercritical CFB boiler. Detailed distributions of solids concentration, oxygen, heat flux, and working medium temperature in the boiler furnace are presented.展开更多
Printed circuit heat exchangers(PCHEs)are promising candidates for recuperators in supercritical CO2 Brayton cycles.A comparative study is given in this paper on the flow and heat transfer characteristics of PCHEs wit...Printed circuit heat exchangers(PCHEs)are promising candidates for recuperators in supercritical CO2 Brayton cycles.A comparative study is given in this paper on the flow and heat transfer characteristics of PCHEs with sinusoidal and zigzag channels.With mass flow rates of 0.6–1.8 kg/h and the bend angles of 15°–30°,the thermal-hydraulic performance of the PCHEs is discussed.Results show that the sinusoidal channel is superior to the zigzag channel in its comprehensive performance.Larger bend angles result in greater reductions in pressure drop if sinusoidal channels are used instead of zigzag channels and a maximum of 48.4%reduction can be obtained in the considered working conditions.Meanwhile,the inlet sections should be carefully optimized since these sections account for up to 31%and 17%of the total pressure drop in the sinusoidal and zigzag channels,respectively.The corner shape of the zigzag channel can be specially designed to further reduce the pressure drop.The nonuniform density and heat flux distributions in both channels are found to be related to the periodic changes of flow directions and the centrifugal forces should not be ignored when optimizing the sinusoidal and zigzag channels.展开更多
In this study,an experimental system was built to investigate the global performance of an 80-k W zigzag printed circuit heat exchanger(PCHE).It could meet the requirement of the pre-cooler for the supercritical carbo...In this study,an experimental system was built to investigate the global performance of an 80-k W zigzag printed circuit heat exchanger(PCHE).It could meet the requirement of the pre-cooler for the supercritical carbon dioxide(S-CO_(2))Brayton power cycle and the modified effectiveness considering the pinch point is between 61.5%and 79.3%.When the outlet S-CO_(2)temperature is near the pseudo-critical point,the thermo-physical properties have more effects on heat transfer performance compared to flow characteristics.For the local performance,the mass flow rates of both sides have crucial influences on the location where the peak of S-CO_(2)Nusselt number occurs while only the S-CO_(2)flow rate affects the variation of the peak value.In addition,the influence of the radius of curvature on the secondary-flow should not be ignored.In the end,new empirical correlations were proposed considering the drastic variations of the Prandtl number.展开更多
In this article,laminar convective heat transfer of a confined slot impinging jet with nanofluid has been numerically investigated over Reynolds number ranges of 200e1000.Two circular ribs are mounted on the lower-tar...In this article,laminar convective heat transfer of a confined slot impinging jet with nanofluid has been numerically investigated over Reynolds number ranges of 200e1000.Two circular ribs are mounted on the lower-target surface:one rib located right the stagnation point and another one located on the left of the stagnation point.SiO2-water nanofluid with nanoparticles volume fraction ranging from 0 to 4%and nanoparticles diameters of 30 nm has been examined.The two-dimensional governing continuity,momentum and energy equations have been solved using finite volume method based on SIMPLE algorithm.The effect of Reynolds number,nanoparticles volume fraction,rib height and rib location on the flow and thermal characteristics are presented and discussed.Results showed that the average Nusselt number,performance factor,total entropy generation as well as friction factor increase with increasing nanoparticles volume fraction.In addition,it is found that the best thermal-hydraulic performance factor is around 1.89 which is obtained at Reynolds number of 1000,nanoparticles volume fraction of 4%,the rib height of 0.1 and the rib location of 2.展开更多
基金Sincere thanks are due to the Shiraz University Research Council for his financial support.
文摘One of the most important aims of this study is to improve the core of the current VVER reactors to achieve more burn-up(or more cycle length)and more intrinsic safety.It is an independent study on the Russian new proposed FAs,called TVS-2M,which would be applied for the future advanced VVERs.Some important aspects of neutronics as well as thermal hydraulics investigations(and analysis)of the new type of Fas are conducted,and results are compared with the standards PWR CDBL.The TVS-2M FA contains gadolinium-oxide which is mixed with UO_(2)(for different Gd densities and U-235 enrichments which are given herein),but the core does not contain BARs.The new type TVS-2M Fas are modeled by the SARCS software package to find the PMAXS format for three states of CZP and HZP as well as HFP,and then the whole core is simulated by the PARCS code to investigate transient conditions.In addition,the WIMS-D5 code is suggested for steady core modeling including TVS-2M FAs and/or TVS FAs.Many neutronics aspects such as the first cycle length(first cycle burn up in terms of MWthd/kgU),the critical concentration of boric acid at the BOC as well as the cycle length,the axial,and radial power peaking factors,differential and integral worthy of the most reactive CPS-CRs,reactivity coefficients of the fuel,moderator,boric acid,and the under-moderation estimation of the core are conducted and benchmarked with the PWR CDBL.Specifically,the burn-up calculations indicate that the 45.6 d increase of the first cycle length(which corresponds to 1.18 MWthd/kgU increase of burnup)is the best improving aim of the new FA type called TVS-2M.Moreover,thermal-hydraulics core design criteria such as MDNBR(based on W3 correlation)and the maximum of fuel and clad temperatures(radially and axially),are investigated,and discussed based on the CDBL.
文摘A theoretical “drift-flux based thermal-hydraulic mixture-fluid coolant channel model” is presented. It is the basis to a corresponding digital “Coolant Channel Module (CCM)”. This purpose derived “Separate-Region Mixture Fluid Approach” should yield an alternative platform to the currently dominant “Separate-Phase Models” where each phase is treated separately. Contrary to it, a direct procedure could be established with the objective to simulate in an as general as possible way the steady state and transient behaviour of characteristic parameters of single- and/or (now non-separated) two-phase fluids flowing within any type of heated or non-heated coolant channels. Their validity could be confirmed by a wide range of verification and validation runs, showing very satisfactory results. The resulting universally applicable code package CCM should provide a fundamental element for the simulation of thermal-hydraulic situations over a wide range of complex systems (such as different types of heat exchangers and steam generators as being applied in both conventional but also nuclear power stations, 1D and 3D nuclear reactor cores etc). Thereby the derived set of equations for different coolant channels (distinguished by their key numbers) as appearing in these systems can be combined with other ODE-s and non-linear algebraic relations from additional parts of such an overall model. And these can then to be solved by applying an appropriate integration routine. Within the solution procedure, however, mathematical discontinuities can arise. This due to the fact that along such a coolant channel transitions from single- to two-phase flow regimes and vice versa could take place. To circumvent these difficulties it will in the presented approach be proposed that the basic coolant channel (BC) is subdivided into a number of sub-channels (SC-s), each of them being occupied exclusively by only a single or a two-phase flow regime. After an appropriate nodalization of the BC (and thus its SC-s) and after applying a “modified finite volume method” together with other special activities the fundamental set of non-linear thermal-hydraulic partial differential equations together with corresponding constitutive relations can be solved for each SC separately. As a result of such a spatial discretization for each SC type (and thus the entire BC) the wanted set of non-linear ordinary differential equations of 1st order could be established. Obviously, special attention had to be given to the varying SC entrance or outlet positions, describing the movement of boiling boundaries or mixture levels along the channel. Including even the possibility of SC-s to disappear or be created anew during a transient.
基金supported by the Chinese TMSR Strategic Pioneer Science and Technology Project(No.XDA02010000)the Frontier Science Key Program of the Chinese Academy of Sciences(No.QYZDY-SSW-JSC016)。
文摘In the core of a molten salt fast reactor(MSFR),heavy metal fuel and fission products can be dissolved in a molten fluoride salt to form a eutectic mixture that acts as both fuel and coolant.Fission energy is released from the fuel salt and transferred to the second loop by fuel salt circulation.Therefore,the MSFR is characterized by strong interaction between the neutronics and the thermal hydraulics.Moreover,recirculation flow occurs,and nuclear heat is accumulated near the fertile blanket,which significantly affects both the flow and the temperature fields in the core.In this work,to further optimize the conceptual geometric design of the MSFR,three geometries of the core and fertile blanket are proposed,and the thermal-hydraulic characteristics,including the three-dimensional flow and temperature fields of the fuel and fertile salts,are simulated and analyzed using a coupling scheme between the open source codes OpenMC and OpenFOAM.The numerical results indicate that a flatter core temperature distribution can be obtained and the hot spot and flow stagnation zones that appear in the upper and lower parts of the core center near the reflector can be eliminated by curving both the top and bottom walls of the core.Moreover,eight cooling loops with a total flow rate of0.0555 m3 s-1 ensur an acceptable temperature distribusure an acceptable temperature distribution in the fertile blanket.
基金supported by Strategic Pilot Science and Technology Project of Chinese Academy of Sciences (No. XD02001005)
文摘In molten salt reactors(MSRs),the liquid fuel salt circulates through the primary loop and a part of the delayed neutron precursors(DNPs)decays outside the reactor core.To model and analyze the flow field effect of DNPs in channel-type liquid-fueled MSRs,a three-dimensional space-time dynamics code,named ThorCORE3D,that couples neutronics,core thermalhydraulics,and a molten salt loop system was developed and validated with the Molten Salt Reactor Experiment(MSRE)benchmarks.The effects of external loop recirculation time,fuel flow rate,and core flow field distribution on the delayed neutron fraction loss of MSRE at steadystate were modeled and simulated using the ThorCORE3D code.Then,the flow field effect of the DNPs on the system responses of the MSRE in the reactivity insertion transient under different initial conditions was analyzed systematically for the channel-type liquid-fueled MSRs.The results indicate that the flow field condition has a significant effect on the steady-state delayed neutron fractions and will further affect the transient power and temperature responses of the reactor system.The analysis results for the effect of the DNP flow field can provide important references for the design optimization and safety analysis of liquid-fueled MSRs.
基金supported by the Chinese TMSR Strategic Pioneer Science and Technology Project(No.XDA02010000)the Frontier Science Key Program of the Chinese Academy of Sciences(No.QYZDY-SSW-JSC016)the National Natural Science Foundation of China Key Program(No.91326201)
文摘A molten salt reactor(MSR) is one of the six advanced reactor concepts selected by the generation Ⅳ international forum because of its advantages of inherent safety, and the promising capabilities of Th-U breeding and transuranics transmutation. A dynamics model for the channel-type MSR is developed in this work based on a three-dimensional thermal–hydraulic model(3DTH) and a point reactor model. The 3DTH couples a three-dimensional heat conduction model and a one-dimensional single-phase flow model that can accurately consider the heat conduction between different assemblies. The 3DTH is validated by the RELAP5 code in terms of the temperature and mass flow distribution calculation. A point reactor model considering the drift of delayed neutron precursors is adopted in the dynamics model. To verify the dynamics model, three experiments from the molten salt reactor experiment are simulated. The agreement of the experimental data and simulation results was excellent.With the aid of this model, the unprotected step reactivity addition and unprotected loss of flow of the 2 MWt experimental MSR are modeled, and the reactor power and temperature evolution are analyzed.
基金This work was supported by Strategic Pilot Science and Technology Project of Chinese Academy of Sciences(No.XD20191031).
文摘Molten salt-cooled reactor is one of the six GenIV reactors with promising characteristics including safety,reliability,proliferation resistance,physical protection,economics,and sustainability.In this paper,a small innovative molten chloride-cooled fast reactor(MCCFR)with 30-year core and a target 120-MWt thermal power was presented.For its feasible study,neutronics,thermal-hydraulics,and radiation damage analysis were performed.The key design properties including kinetics parameters,reactivity swing,reactivity feedback coefficients,maximum accumulated displacement per atom(DPA)of reactor pressure vessel(RPV)and fuel cladding,and maximum coolant,cladding,and fuel temperatures were evaluated.The results showed the MCCFR could operate without refueling for 30 years with overall negative reactivity feedback coefficients up the end of its life.During its 30-year life,the excess reactivity was well managed by constantly pulling out the control rods.The maximum accumulated DPA on RPV and fuel cladding were 8.92 dpa and 197.03 dpa,respectively,which are both below the design limits.Similarly,the maximum coolant,cladding and fuel center temperatures were all below the design limits during its entire lifetime.According to these results,the MCCFR core design with long life is feasible.
基金supported by the China Postdoctoral Science Foundation(No.2021M703045)the National Natural Science Foundation of China(No.12075067)the National Key R&D Program of China(No.2018YFE0180900).
文摘To perform nuclear reactor simulations in a more realistic manner,the coupling scheme between neutronics and thermal-hydraulics was implemented in the HNET program for both steady-state and transient conditions.For simplicity,efficiency,and robustness,the matrixfree Newton/Krylov(MFNK)method was applied to the steady-state coupling calculation.In addition,the optimal perturbation size was adopted to further improve the convergence behavior of the MFNK.For the transient coupling simulation,the operator splitting method with a staggered time mesh was utilized to balance the computational cost and accuracy.Finally,VERA Problem 6 with power and boron perturbation and the NEACRP transient benchmark were simulated for analysis.The numerical results show that the MFNK method can outperform Picard iteration in terms of both efficiency and robustness for a wide range of problems.Furthermore,the reasonable agreement between the simulation results and the reference results for the NEACRP transient benchmark verifies the capability of predicting the behavior of the nuclear reactor.
文摘Previous analytical results on flow splitting are generalized to consider multiple boiling channels systems. The analysis is consistent with the approximations usually adopted in the use of systems codes (like RELAP5 and TRACE5, among others) commonly applied to perform safety analyses of nuclear power plants. The problem is related to multiple, identical, parallel boiling channels, connected through common plena. A theoretical model limited in scope explains this flow splitting without reversal. The unified analysis performed and the confirmatory computational results found are summarized in this paper. New maps showing the zones where this behavior is predicted are also shown considering again twin pipes. Multiple pipe systems have been found not easily amenable for analytical analysis when dealing with more than four parallel pipes. However, the particular splitting found (flow along N pipes dividing in one standalone pipe flow plus N -1 identical pipe flows) has been verified up to fourteen pipes, involving calculations in systems with even and odd number of pipes using the RELAP5 systems thermal-hydraulics code.
文摘In this study,COMSOL v5.2 Multiphysics software was utilized to perform coupled neutronics and thermal–hydraulics simulations of a molten salt fast reactor,and the SCALE v6.1 code package was utilized to generate the homogenized cross-section data library.The library’s 238 cross-section groups were categorized into nine groups for the simulations in this study.The results of the COMSOL model under no fuel flow conditions were verified using the SCALE v6.1 code results,and the results of the neutronics and thermal–hydraulics simulations were compared to the results of previously published studies.The results indicated that the COMSOL model that includes the cross-section library generated by the SCALE v6.1 code package is suitable for the steady-state analysis and design assessment of molten salt fast reactors.Subsequently,this model was utilized to investigate the neutronics and thermal–hydraulics behaviors of the reactor.Multiple designs were simulated and analyzed in this model,and the results indicated that even if the wall of the core is curved,hot spots occur in the upper and lower portions of the core’s center near the reflectors.A new design was proposed that utilizes a flow rate distribution system,and the simulation results of this design showed that the maximum temperature in the core was approximately 1032 K and no hot spots occurred.
基金supported by the National Natural Science Foundation of China(51875238)。
文摘In this paper,a novel composite heat transfer enhancement technique comprised of louvered fins(LFs)and rectangular wing vortex generators(RWVGs)is proposed to improve the LF side thermal-hydraulic performance of louvered fin and flat tube heat exchangers(LFHEs).After validation of the LF side pressure dropΔP and heat transfer coefficient hLFof the baseline by experiments,the numerical method is applied to investigate the influential mechanisms of the RWVG parameters(the number N(7 to 15),attack angleβ(30°to 90°),height H_(VG)(0.8 mm to 2 mm)and width W_(VG)(0.8 mm to 1.2 mm))on the performance of the LFHE in the velocity range of 3 m/s to 10 m/s.Results show that thermal-hydraulic performance of the LFHE is significantly impacted by the RWVGs,and according to the performance evaluation criteria(PEC),the LFHE achieves its optimal thermal-hydraulic performance when N=7,β=45°,H_(VG)=1.8 mm and W_(VG)=1 mm.Compared to the baseline,the maximum,minimum and average increments of PEC for the optimal case are 13.85%,4.67%and 8.39%,respectively.
文摘Following the previous relevant works [1]-[3], a scaling analysis is performed to derive a set of scaling criteria which were thought to be suitable for reproducing the major thermal-hydraulic phenomena in a scaled-down CANDU moderator tank similar to that in a prototype power plant during a full power steady state condition. The objective of building a scaled-down moderator tank is to generate the experimental data necessary to validate the computer codes which are used to analyze the accident analysis of CANDU-6 plants. The major variables of interests in this paper are moderator flow velocity and temperature of the moderator which is D2O inside the moderator tank during a steady state and transient conditions. The reason is that the local subcooling of the moderator is found to be a critical parameter determining whether the stable film boiling can sustain on the outer surface of the calandria tube if the contact of overheated pressure tube and cold calandria tube should occur due to pressure tube ballooning during LBLOCA with ECC injection failure [4]. The key phenomena involved include the inlet jet development and impingement, buoyancy force driven by the moderator temperature gradient caused by non-uniform direct heating of the moderator, and the pressure drop due to viscous friction of the flow across the calandria tube array. In this paper, the previous researches are reviewed, some concerns or potential problems associated with them implied by comparing CFD analyses results between the CANDU-6 moderator tank and 1/4 scaled-down test facility are described, and as a way to examine the assumption of the scaling analysis is true an order-of-magnitude analyses are performed. Based on the results of these analyses the assumption of neglecting ?and ?terms cannot be justified for the power of 0.5 MW and 1.566 MW for the 1/4 scaled-down facility. Further investigation is thought to be necessary to confirm this result, i.e. if the scaling of the previous work1 is justifiable by some other independent analyses.
文摘Development and use of nuclear energy is currently growing very rapidly, in order to achieve increasingly advanced technology, both in terms of design, economic factors and safety factors. Thermal-hydraulics aspects of nuclear reactors should be done with calculation and near-perfect condition. Including today began development of a nuclear reactor with low power below 300 MW, or commonly called the Small Modular Reactor (SMR). One is CAREM-25 developed by Argentina with a power of 25 MW, where in CAREM already using natural circulation system and the use of nanofluid as coolant fluid. In this research, analytic modeling of thermal-hydraulics nuclear reactor SMR CAREM-25, when the nanofluid Al<sub>2</sub>O<sub>3</sub>-Water used as cooling fluid in the cooling system of a nuclear reactor. Further to this analytic modeling will be done on CFD. Analytic modeling with CFD to determine the flow phenomena and distribution as well as the effect of nano-particles of Al<sub>2</sub>O<sub>3</sub>-Water based on the volume fraction (1% and 3%) of the coefficient of heat transfer by natural convection.
基金supported by the School of Mechanical,Institute of Engineering,Suranaree University of Technology (Thailand),Mechanical and Process System Engineering Program,and Vithedbundit Scholarship,Institute of Engineering,Suranaree University of Technology (Thailand)。
文摘This study investigated the effects of zigzag-flow channel bending angle in printed circuit heat exchangers(PCHEs) using a computational fluid dynamics method with ANSYS-FLUENT simulation.The three-dimensional model of PCHE with a 15° curved,zigzag channel was conducted for preliminary validation.The comparisons between the CFD simulation results and the experimental data showed good agreement with some discrepancies in the heat transfer and pressure drop results.In addition,different bending angle configurations(0°,3° to 30°) of zigzag channels were analyzed to obtain better thermal-hydraulic performance of the zigzag channel PCHE under different inlet mass flow rates.The criteria of heat transfer and frictional factor were applied to evaluate the thermal-hydraulic performance of the PCHE.The results showed that the 6° and 9°bending channel provided good thermal-hydraulic performance.New correlations were developed using the 6°and 9° bending channel angles in PCHE designs to predict the Nusselt number and friction factor.
基金financially supported by the National Key Research&Development Program of China(2018YFB0604400)。
文摘A double reheat ultra-supercritical boiler is an important development direction for high-parameter and large-capacity coal-fired power plants due to its high thermal efficiency and environmental value.China has developed a 1000 MW double reheat ultra-supercritical boiler with steam parameters of 35 MPa at 605℃/613℃/613℃.Reasonable water wall design is one of the keys to safe and reliable operation of the boiler.In order to examine the thermal-hydraulic characteristics of the double reheat ultra-supercritical boiler,the water wall system was equivalent to a flow network comprising series-parallel circuits,linking circuits and pressure nodes,and a calculation model was built on account of the conservation equations of energy,momentum and mass.Through the iterative solving of nonlinear equations,the prediction parameters of the water wall at boiler maximum continue rate(BMCR),75%turbine heat-acceptance rate(THA)and 30%THA loads,including total pressure drops,flow distribution,outlet steam temperatures,fluid and metal temperatures were gotten.The results of calculation exhibit excellent thermal-hydraulic characteristics and substantiate the feasibility of the water wall design of the double reheat ultra-supercritical boiler.
基金financial support provided by the National Natural Science Foundation of China(No.52076075 and No.52130608)。
文摘In terms of developing supercritical CO_(2)(sCO_(2))coal-fired power plants,enhancing cooling wall performance is one of significant factors to improve system performance.In this paper,a new cooling wall tube structure is proposed to match the non-uniform heat flux(NUH)with the thermal resistance by changing the cooling wall tube eccentricity.A three-dimensional multi-physical coupling model of cooling wall is constructed to compare the novel structure to the conventional structures.The properties of fluid dynamics,thermal stress,coupled heat transfer and cooling wall deformation are analyzed.In contrast to the traditional structure,the maximum temperature and circumferential temperature difference(CTD)of the proposed structure can be reduced by 2%and 27.4%,respectively.The essential working parameters related to the performances of the cooling wall tube are discussed.The maximum temperature of the new structure is reduced by 8-13 K and the maximum thermal stress is reduced by about 10%-15%under all the simulated working conditions when the eccentricity changes from 0 to 0.2.The proposed structure can effectively reduce the maximum temperature and circumferential temperature gradient under NUH.Consequently,a novel insight is put out for the design and optimization of the cooling wall tube in coal-fired power plants.
文摘A combustion model of a large-scale supercritical circulati ng fluidized bed (CFB) boiler was developed for comprehensive computational-fluid-dynamics analysis. The model incorporates gas-solid hydrodynamics, coal combustion, heat transfer on heat exchange surfaces in the furnace, and heat transfer between fumace and working medium in the heat transfer tubes. In simulating the dense and dilute phases in the fumace, the gas-solid hydrodynamics is based on the Euler-Euler model and energy-minimization multiscale drag model. Coal combustion entails evaporation, devolatilization, char combustion, gas homoge neous reaction, and pollutant emission. The coefficient ofheat transfer between gas-solid and the waterwall is estimated using the cluster renewal model, and for radiation, the discrete ordinate model is used. Moreover, thermohydraulic processes in the membrane wall are also in eluded in the heat transfer process. The model was successfully applied in simulations of a 350-MW supercritical CFB boiler. Detailed distributions of solids concentration, oxygen, heat flux, and working medium temperature in the boiler furnace are presented.
基金supported by the National Natural Science Foundation of China(Grant Nos.51822606,51806249)Hunan Provincial Natural Science Foundation of China(Grant No.2019JJ50801)。
文摘Printed circuit heat exchangers(PCHEs)are promising candidates for recuperators in supercritical CO2 Brayton cycles.A comparative study is given in this paper on the flow and heat transfer characteristics of PCHEs with sinusoidal and zigzag channels.With mass flow rates of 0.6–1.8 kg/h and the bend angles of 15°–30°,the thermal-hydraulic performance of the PCHEs is discussed.Results show that the sinusoidal channel is superior to the zigzag channel in its comprehensive performance.Larger bend angles result in greater reductions in pressure drop if sinusoidal channels are used instead of zigzag channels and a maximum of 48.4%reduction can be obtained in the considered working conditions.Meanwhile,the inlet sections should be carefully optimized since these sections account for up to 31%and 17%of the total pressure drop in the sinusoidal and zigzag channels,respectively.The corner shape of the zigzag channel can be specially designed to further reduce the pressure drop.The nonuniform density and heat flux distributions in both channels are found to be related to the periodic changes of flow directions and the centrifugal forces should not be ignored when optimizing the sinusoidal and zigzag channels.
基金supported by the National Natural Science Foundation of China(No.51606191)the National Key Research and Development Program-China(2017YFB0601803)Key deployment project of Chinese Academy of Sciences(Y7220112H1)。
文摘In this study,an experimental system was built to investigate the global performance of an 80-k W zigzag printed circuit heat exchanger(PCHE).It could meet the requirement of the pre-cooler for the supercritical carbon dioxide(S-CO_(2))Brayton power cycle and the modified effectiveness considering the pinch point is between 61.5%and 79.3%.When the outlet S-CO_(2)temperature is near the pseudo-critical point,the thermo-physical properties have more effects on heat transfer performance compared to flow characteristics.For the local performance,the mass flow rates of both sides have crucial influences on the location where the peak of S-CO_(2)Nusselt number occurs while only the S-CO_(2)flow rate affects the variation of the peak value.In addition,the influence of the radius of curvature on the secondary-flow should not be ignored.In the end,new empirical correlations were proposed considering the drastic variations of the Prandtl number.
文摘In this article,laminar convective heat transfer of a confined slot impinging jet with nanofluid has been numerically investigated over Reynolds number ranges of 200e1000.Two circular ribs are mounted on the lower-target surface:one rib located right the stagnation point and another one located on the left of the stagnation point.SiO2-water nanofluid with nanoparticles volume fraction ranging from 0 to 4%and nanoparticles diameters of 30 nm has been examined.The two-dimensional governing continuity,momentum and energy equations have been solved using finite volume method based on SIMPLE algorithm.The effect of Reynolds number,nanoparticles volume fraction,rib height and rib location on the flow and thermal characteristics are presented and discussed.Results showed that the average Nusselt number,performance factor,total entropy generation as well as friction factor increase with increasing nanoparticles volume fraction.In addition,it is found that the best thermal-hydraulic performance factor is around 1.89 which is obtained at Reynolds number of 1000,nanoparticles volume fraction of 4%,the rib height of 0.1 and the rib location of 2.