• Nuclear Engineering and Technology
• /
• v.45 no.4
• /
• pp.439-458
• /
• 2013

REX-10 is a fully-passive small modular reactor in which the coolant flow is driven by natural circulation, the RCS is pressurized by a steam-gas pressurizer, and the decay heat is removed by the PRHRS. To confirm design decisions and analyze the transient responses of an integral PWR such as REX-10, a thermal-hydraulic system code named TAPINS (Thermal-hydraulic Analysis Program for INtegral reactor System) is developed in this study. Based on a one-dimensional four-equation drift-flux model, TAPINS incorporates mathematical models for the core, the helical-coil steam generator, and the steam-gas pressurizer. The system of difference equations derived from the semi-implicit finite-difference scheme is numerically solved by the Newton Block Gauss Seidel (NBGS) method. TAPINS is characterized by applicability to transients with non-equilibrium effects, better prediction of the transient behavior of a pressurizer containing non-condensable gas, and code assessment by using the experimental data from the autonomous integral effect tests in the RTF (REX-10 Test Facility). Details on the hydrodynamic models as well as a part of validation results that reveal the features of TAPINS are presented in this paper.

• Journal of the Korean Geotechnical Society
• /
• v.37 no.12
• /
• pp.25-31
• /
• 2021

Boundary conditions for thermal-hydraulic problems of soils play an essential role in the numerical accuracy. This study presents a boundary condition considering the thermo-hydraulic interaction between the ground and the atmosphere. Ground surface energy balance consists of solar radiation, ground radiation, wind convection, latent heat from water evaporation, and heat conduction to the ground. Equations for each heat flux are presented, and numerical analyses are performed in conjunction with the FEM program for the thermal-hydraulic phenomenon of unsaturated soils. Numerical results using the weather data at the Ulsan Meteorological Observatory are similar to the measured surface temperature. Latent heat caused by water evaporation during the daytime lowers the surface temperature of the bare soil, and a thermal equilibrium is reached at nighttime when the effect of the ground condition is significantly reduced. The temperature change of the surface ground is diminished at the deeper ground due to its thermal diffusion. Numerical analysis where the surface ground temperature is the primary concern requires considering the thermo-hydraulic interaction between the ground and the atmosphere.

• Journal of computational fluids engineering
• /
• v.22 no.1
• /
• pp.26-36
• /
• 2017

A numerical study has been performed to evaluate thermal-hydraulic performance of a finger type cooling module with multiple-jet impingement in a divertor of nuclear fusion reactor. To analyze conjugate heat transfer in both solid and fluid domains, numerical analysis of the flow using three-dimensional Reynolds-averaged Navier-Stokes equations has been performed with shear stress transport turbulence model. The computational domain for the cooling module consisted of a single fluid domain and three solid domains; tile, thimble, and cartridge. The numerical results for the temperature variation on the tile were validated in comparison with experimental data under the same conditions. A parametric study was performed with four geometric parameters, i.e., angles between x-axis and centerlines of hole 1, 2, 3 and 4. The results indicate that the heat transfer rate was increased by 2.7% and 0.7% by the angle ${\theta}_1$ and angle ${\theta}_2$, respectively, and that the pressure drop was decreased by up to 1.8% by the angle ${\theta}_3$.

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.271-278
• /
• 2020

The investigation into a full-scale 27 m high, by 6 m wide, thermosyphon loop. The simulation model is based on a one-dimensional axially-symmetrical control volume approach, where the loop is divided into a series of discreet control volumes. The three conservation equations, namely, mass, momentum and energy, were applied to these control volumes and solved with an explicit numerical method. The flow is assumed to be quasi-static, implying that the mass-flow rate changes over time. However, at any instant in time the mass-flow rate is constant around the loop. The boussinesq approximation was invoked, and a reasonable correlation between the experimental and theoretical results was obtained. Experimental results are presented and the flow regimes of the working fluid inside the loop identified. The results indicate that a series of such thermosyphon loops can be used as a cavity cooling system and that the one-dimensional theoretical model can predict the internal temperature and mass-flow rate of the thermosyphon loop.

• Advances in Energy Research
• /
• v.5 no.2
• /
• pp.91-105
• /
• 2017

Thermal hydraulic (TH) analysis of nuclear power reactors is utmost important. In this way, the numerical codes that preparing TH data in reactor core are essential. In this paper, a subchannel analysis of a Russian pressurized water reactor (WWER1000) core with enhanced numerical code is carried out. For this, in fluid domain, the mass, axial and lateral momentum and energy conservation equations for desired control volume are solved, numerically. In the solid domain, the cylindrical heat transfer equation for calculation of radial temperature profile in fuel, gap and clad with finite difference and finite element solvers are considered. The dependence of material properties to fuel burnup with Calza-Bini fuel-gap model is implemented. This model is coupled with Isotope Generation and Depletion Code (ORIGEN2.1). The possibility of central hole consideration in fuel pellet is another advantage of this work. In addition, subchannel to subchannel and subchannel to rod connection data in hexagonal fuel assembly geometry could be prepared, automatically. For a demonstration of code capability, the steady state TH analysis of a the WWER1000 core is compromised with Thermal-hydraulic analysis code (COBRA-EN). By thermal hydraulic parameters averaging Fuel Assembly-to-Fuel Assembly method, the one sixth (symmetry) of the Boushehr Nuclear Power Plant (BNPP) core with regular subchannels are modeled. Comparison between the results of the work and COBRA-EN demonstrates some advantages of the presented code. Using the code the thermal modeling of the fuel rods with considering the fission gas generation would be possible. In addition, this code is compatible with neutronic codes for coupling. This method is faster and more accurate for symmetrical simulation of the core with acceptable results.

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.230-237
• /
• 2020

In this part, an implicit time dependent solution is presented for the Boltzmann transport equation discretized by the analytic coarse mesh finite difference method (ACMFD) over the spatial domain as well as the simplified P₃ (SP₃) for the angular variable. In the first part of this work we proposed a SP₃-ACMFD approach to solve the static eigenvalue equations which provide the initial conditions for temp dependent equations. Having solved the 3D multi-group SP₃-ACMFD static equations, an implicit approach is resorted to ensure stability of time steps. An exponential behavior is assumed in transverse integrated equations to establish a relationship between flux moments and currents. Also, analytic integration is benefited for the time-dependent solution of precursor concentration equations. Finally, a multi-channel one-phase thermal hydraulic model is coupled to the proposed methodology. Transient equations are then solved at each step using the GMRES technique. To show the sufficiency of proposed transient SP₃-ACMFD approximation for a full core analysis, a comparison is made using transport peers as the reference. To further demonstrate superiority, results are compared with a 3D multi-group transient diffusion solver developed as a byproduct of this work. Outcomes confirm that the idea can be considered as an economic interim approach which is superior to the diffusion approximation, and comparable with transport in results.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1997.10a
• /
• pp.599-604
• /
• 1997

One of important inconsistencies in the six-equation model predictions has been found to be the force experienced by a single bubble placed in a convergent stream of liquid. Various sets of governing equations yield different amount of forces to hold the bubble stationary in a convergent nozzle. By using the first order potential flow theory, it is found that the six-equation model can not be used to estimate the force experienced by a deformed bubble. The theoretical value of the particle stress of a bubble in a convergent nozzle flow has been found to be a function of the Weber number when bubble distortion is allowed. This force has been calculated by using different sets of governing equations and compared with the theoretical value. It is suggested in this study that the bubble size distribution function can be used to remove the presented inconsistency by relating the interfacial variables with different moments of the bubble size distribution function. This study also shows that the inconsistencies in the thermal-hydraulic governing equation can be removed by mechanistic modeling of the phasic interface.

• Nuclear Engineering and Technology
• /
• v.52 no.9
• /
• pp.2017-2024
• /
• 2020

In this study, a fractional order proportional integral derivative (FOPID) controller is designed to create the reference power trajectory and to conquer the uncertainties and external disturbances. A fractional nonlinear model was utilized to describe the nuclear reactor dynamic behaviour considering thermal-hydraulic effects. The controller parameters were tuned using optimization method in Matlab/Simulink. The FOPID controller was simulated using Matlab/Simulink and the controller performance was evaluated for Hard variation of the reference power and compared with that of integer order a proportional integral derivative (IOPID) controller by two models of fractional neutron point kinetic (FNPK) and classical neutron point kinetic (CNPK). Also, the FOPID controller robustness was appraised against the external disturbance and uncertainties. Simulation results showed that the FOPID controller has the faster response of the control attempt signal and the smaller tracking error with respect to the IOPID in tracking the reference power trajectory. In addition, the results demonstrated the ability of FOPID controller in disturbance rejection and exhibited the good robustness of controller against uncertainty.

• Nuclear Engineering and Technology
• /
• v.51 no.7
• /
• pp.1721-1728
• /
• 2019

The advanced node core code CORCA-3D is one of the independent developed codes of NPIC for the nuclear reactor core design. CORCA-3D code can calculate the few-group cross section, solve the 3D diffusion equations, consider the thermal-hydraulic feedback, reconstruct the pin-by-pin power. It has lots of functions such as changing core status calculation, critical searching, control rod value calculation, coefficient calculation and so on. The main theory and functions of CORCA-3D code are introduced and validated with a lot of reactor measured data and the SCIENCE system. Now, CORCA-3D code has been applied in ACP type reactor nuclear cores design.

• Nuclear Engineering and Technology
• /
• v.54 no.9
• /
• pp.3526-3539
• /
• 2022

Nuclear thermal-hydraulic system analysis codes have been developed to simulate nuclear reactor systems, which solve simplified governing equations by replacing source terms with constitutive relations for simulating entire reactor systems with low computational resources. For half a century, many efforts have been made for wider versatility and higher accuracy of system codes, but various factors can affect the code analysis results, and it was difficult to isolate these factors and interpret them individually. In this study, two system codes, RELAP5 and TRACE, which have many users and are highly reliable, are selected to analyze only the effects of constitutive relations. The influence of constitutive relations is analyzed using in-house platforms that replicate constitute relations of RELAP5 and TRACE equally to exclude factors that may affect analysis results, such as governing equation solvers and user effects. Among the various constitutive relations, the analysis is performed on the wall variables expected to have the most influence on the analysis results. Part 1 paper presents the methodology and wall friction model comparison, while Part 2 paper shows wall heat transfer comparison of the two selected codes.