• Title/Summary/Keyword: heat equation

Search Result 1,383, Processing Time 0.025 seconds

Verification of neutronics and thermal-hydraulic coupled system with pin-by-pin calculation for PWR core

  • Zhigang Li;Junjie Pan;Bangyang Xia;Shenglong Qiang;Wei Lu;Qing Li
    • Nuclear Engineering and Technology
    • /
    • v.55 no.9
    • /
    • pp.3213-3228
    • /
    • 2023
  • As an important part of the digital reactor, the pin-by-pin wise fine coupling calculation is a research hotspot in the field of nuclear engineering in recent years. It provides more precise and realistic simulation results for reactor design, operation and safety evaluation. CORCA-K a nodal code is redeveloped as a robust pin-by-pin wise neutronics and thermal-hydraulic coupled calculation code for pressurized water reactor (PWR) core. The nodal green's function method (NGFM) is used to solve the three-dimensional space-time neutron dynamics equation, and the single-phase single channel model and one-dimensional heat conduction model are used to solve the fluid field and fuel temperature field. The mesh scale of reactor core simulation is raised from the nodal-wise to the pin-wise. It is verified by two benchmarks: NEACRP 3D PWR and PWR MOX/UO2. The results show that: 1) the pin-by-pin wise coupling calculation system has good accuracy and can accurately simulate the key parameters in steady-state and transient coupling conditions, which is in good agreement with the reference results; 2) Compared with the nodal-wise coupling calculation, the pin-by-pin wise coupling calculation improves the fuel peak temperature, the range of power distribution is expanded, and the lower limit is reduced more.

Numerical study to reproduce a real cable tray fire event in a nuclear power plant

  • Jaiho Lee ;Byeongjun Kim;Yong Hun Jung;Sangkyu Lee;Weon Gyu Shin
    • Nuclear Engineering and Technology
    • /
    • v.55 no.4
    • /
    • pp.1571-1584
    • /
    • 2023
  • In this study, a numerical analysis was performed as part of an international joint research project to reproduce a real cable tray fire that occurred in the heater bay area of the turbine building of a nuclear power plant. A sensitivity analysis was performed on various input parameters to derive results consistent with the sprinkler activation time obtained from the fire event analysis. For all sensitive parameters, the normalized sprinkler activation time correlated well with the power function of the normalized sprinkler height. A correlation equation was developed to identify the sprinkler activation time at any location when determining the slope or fire growth rate under the conditions assuming a linear or t-squared heat release rate (HRR) time curve. Various cable fire growth assumptions were used to determine which assumption was better to provide the prediction coincident with the information given from the fire event analysis in terms of the sprinkler activation time and total energy generated from cables damaged by fire. In the comprehensive analysis of all the sensitive parameters, the standard deviation of the input parameters increased as the sprinkler height decreased. Within the range of the sensitivity parameter values given in this study, when considering all sprinkler heights, the standard deviation of the cable model change was the largest and that of the overhang position change was the smallest.

Experimental and numerical assessment of helium bubble lift during natural circulation for passive molten salt fast reactor

  • Won Jun Choi;Jae Hyung Park;Juhyeong Lee;Jihun Im;Yunsik Cho;Yonghee Kim;Sung Joong Kim
    • Nuclear Engineering and Technology
    • /
    • v.56 no.3
    • /
    • pp.1002-1012
    • /
    • 2024
  • To remove insoluble fission products, which could possibly cause reactor instability and significantly reduce heat transfer efficiency from primary system of molten salt reactor, a helium bubbling method is employed into a passive molten salt fast reactor. In this regard, two-phase flow behavior of molten salt and helium bubbles was investigated experimentally because the helium bubbles highly affect the circulation performance of working fluid owing to an additional drag force. As the helium flow rate is controlled, the change of key thermal-hydraulic parameters was analyzed through a two-phase experiment. Simultaneously, to assess the applicability of numerical model for the analysis of two-phase flow behavior, the numerical calculation was performed using the OpenFOAM 9.0 code. The accuracy of the numerical analysis code was evaluated by comparing it with the experimental data. Generally, numerical results showed a good agreement with the experiment. However, at the high helium injection rates, the prediction capability for void fraction of helium bubbles was relatively low. This study suggests that the multiphaseEulerFoam solver in OpenFOAM code is effective for predicting the helium bubbling but there exists a room for further improvement by incorporating the appropriate drag flux model and the population balance equation.

Stability analysis in BWRs with double subdiffusion effects: Reduced order fractional model (DS-F-ROM)

  • Gilberto Espinosa-Paredes;Ricardo I. Cazares-Ramirez;Vishwesh A. Vyawahare;Erick-G. Espinosa-Martinez
    • Nuclear Engineering and Technology
    • /
    • v.56 no.4
    • /
    • pp.1296-1309
    • /
    • 2024
  • The aim of this work is to explore the effect of the double subdiffusion on the stability in BWRs. A BWR novel reduced order model with double subdiffusion effects: reduced order fractional model (DS-F-ROM) to describe the neutron and heat transfer processes was proposed for this study. The double subdiffusion was developed with a fractional-order two-equation model, and with different fractional-orders and relaxation times. The stability analysis was carried out using the root-locus method and change from the s to the W domain and were confirmed using the time-domain evolution of neutron flux for a unit step change in reactivity. The results obtained using the reduced fractional-order model are presented for different anomalous diffusion coefficient values. Results are compared with normal diffusion and P1 equations, which are obtained straightforwardly with DS-ROM when relaxation time tends to zero, and when the anomalous diffusion coefficient tends to one, respectively.

Theoretical Study on Observed Heat of Ligation for Iron(Ⅱ) and Nickel(Ⅱ) Octahedral Complexes (팔면체형 Fe(Ⅱ)와 Ni(Ⅱ)착물의 실측 리간드화열에 관한 이론적 연구)

  • Kim, Jung Sung;Choi, Jin Tae;Song, Young Dae;Cho, Tae Sub
    • Journal of the Korean Chemical Society
    • /
    • v.43 no.2
    • /
    • pp.141-149
    • /
    • 1999
  • The correlation was investigated between the observed heat of ligation and calculated quantum chemical quantities for octahedral $[M(H_2O)_{6-x}(NH_3)_x]^{2+} (M=Fe(II),\;Ni(II))$ complexes by EHMO(Extended Huckel Molecular Orbital) and ZINDO/1(Zerner's Intermediate Neglected of Differential Overlap)method. The net charge of $Fe^{2+}$ and $Ni^{2+}$ ion of octahedral $[M(H_2O)_{6-x}(NH_3)_x]^{2+}(M=Fe(II),\;Ni(II))$ complexes(x=O, 1, …, 6) decreased with substituting $NH_3$ for $H_2O$ molecules. It has found that a good correlation exists between the observed heat of ligation and the calculated quantum chemical quantities such as net charge of central atom, enthalpy of formation, and total dissociation energy. From this finding, we have obtained the following semiempirical linear equation ${\Delta}H_{obs}=-0.2858_{qFe}+0.8813(r=0.97),\;{\Delta}H_{obs}=-0.8981_{qNi}+1.7929(r=0.95),\;{\Delta}H_{obs}=-0.0031H_{f(Fe)}+0.5725(r=0.97),\;{\Delta}H_{obs}=-0.0095H_{f(Ni)}+0.9193(r=0.97),\;{\Delta}H_{obs}=0.0476E_{diss(Fe)}+0.6434(r=0.94),\;{\Delta}H_{obs}=0.1401E_{diss(Ni)}+1.1393(r=0.93)$.

  • PDF

Adsorption Characteristics of CO2 on Activated Carbons Treated with Alkali-metal Salts (알칼리금속염으로 처리된 활성탄에 대한 CO2의 흡착특성)

  • Ryu, Dong Kwan;Kim, Sung Hyun
    • Applied Chemistry for Engineering
    • /
    • v.9 no.2
    • /
    • pp.286-293
    • /
    • 1998
  • Two methods were used to enhance the adsorption capacity of activated carbons. One is to impregnate activated carbons with chemical compounds which have a good affinity for $CO_2$. The other is to activate by heat-treating after impregnation with KOH on activated carbons(AC). The chemical compounds impregnated on AC were alkali metal, alkaline earth metal, and transition metal chlorides. The adsorption capacity of $CO_2$ on AC impregnated with these metals was less than that of pure AC. These compounds have not the chemical affinity for $CO_2$ and obstruct the micropore of AC. The experiment of breakthrough for $CO_2$ on AC impregnated with KOH showed the increase of the adsorbed amount of $CO_2$ in influent gases containing water vapor. This means that KOH adsorbes $CO_2$ gas. However, the adsorbents impregnated with KOH had not the reproducibility because of the production of $K_2CO_3$ by the reaction of KOH with $CO_2$. The amount of $CO_2$ adsorbed on the heat-treated AC at $800^{\circ}C$ increased with the amount of impregnation. The adsorption capacity of $CO_2$ was the largest when the ratio of weight of KOH to AC equal to 4. The isosteric heat of adsorption was calculated by the equation of Clausius-Clapeyron form adsorption capacity data of $CO_2$ for the temperature change. In addition, the characteristics of $CO_2$ breakthrough curve were surveyed for the change of flow rate and concentration.

  • PDF

Flow and smoke behavior of a longitudinal ventilation tunnel with various velocities using computational fluid dynamics (팬의 운전조건에 따른 종류식환기터널 내의 연기거동에 관한 전산유체역학연구)

  • Lee, J.H.;Kwon, Y.J.;Kim, D.E.
    • Journal of Korean Tunnelling and Underground Space Association
    • /
    • v.16 no.1
    • /
    • pp.105-115
    • /
    • 2014
  • A numerical analysis on the smoke behavior and evacuee safety has been performed with computational fluid dynamics. The purpose of this study is to build computational processes for an evacuation and prevention of a fire disaster of a 3 km-length tunnel in Korea. To save computational cost, 1.5 km of the tunnel that can include a few cross-passing tunnels is considered. We are going to assess the fire safety in a road tunnel according to the smoke level, which consists of the smoke density and the height from the floor. The smoke density is obtained in detail from three-dimensional unsteady CFD analysis. To obtain proper temperature distributions on the tunnel wall, one-dimensional conduction equation is considered instead of an adiabatic wall boundary or a constant heat flux. The tunnel considered in this study equips the cross passing tunnels for evacuees every 250 m. The distance is critical in both safety and economy. The more cross passing tunnels, the more safe but the more expensive. Three different jet fan operations can be considered in this study; under- and over-critical velocities for normal traffic condition and 0-velocoty operation for the traffic congestion. The SE (smoke environment) level maps show a smoke environment and an evacuating behavior every moment.

Characteristics of Isotherm, Kinetic and Thermodynamic Parameters for the Adsorption of Acid Red 66 by Activated Carbon (활성탄에 의한 Acid Red 66의 흡착에 대한 등온선, 동력학 및 열역학적 특성)

  • Lee, Jong-Jib
    • Clean Technology
    • /
    • v.26 no.1
    • /
    • pp.30-38
    • /
    • 2020
  • The kinetic and thermodynamic parameters of Acid Red 66, adsorbed by granular activated carbon, were investigated on areas of initial concentration, contact time, and temperature. The adsorption equilibrium data were applied to Langmuir, Freundlich, Temkin, Redlich-Peterson, and Temkin isotherms. The agreement was found to be the highest in the Freundlich model. From the determined Freundlich separation factor (1/n = 0.125 ~ 0.232), the adsorption of Acid Red 66 by granular activated carbon could be employed as an effective treatment method. Temkin's constant related to adsorption heat (BT = 2.147 ~ 2.562 J mol-1) showed that this process was physical adsorption. From kinetic experiments, the adsorption process followed the pseudo-second order model with good agreement. The results of the intraparticle diffusion equation showed that the inclination of the second straight line representing the intraparticle diffusion was smaller than that of the first straight line representing the boundary layer diffusion. Therefore, it was confirmed that intraparticle diffusion was the rate-controlling step. From thermodynamic experiments, the activation energy was determined as 35.23 kJ mol-1, indicating that the adsorption of Acid Red 66 was physical adsorption. The negative Gibbs free energy change (ΔG = -0.548 ~ -7.802 kJ mol-1) and the positive enthalpy change (ΔH = +109.112 kJ mol-1) indicated the spontaneous and endothermic nature of the adsorption process, respectively. The isosteric heat of adsorption increased with the increase of surface loading, indicating lateral interactions between the adsorbed dye molecules.

A Correlation Study for the Prediction of the Maximum Heat Release Rate in Closed-Compartments of Various Configurations (다양한 형상의 밀폐된 구획에서 최대 열발생률 예측을 위한 상관식 검토)

  • Yun, Hong-Seok;Hwang, Cheol-Hong
    • Fire Science and Engineering
    • /
    • v.32 no.1
    • /
    • pp.16-23
    • /
    • 2018
  • In a closed-compartment with various configurations, the correlation that can predict the maximum heat release rate (HRR) with the changes in internal volume and fire growth rate was investigated numerically. The volume of the compartment was controlled by varying the length ratio based on the bottom surface shape of the ISO 9705 fire room, where the ceiling height was fixed to 2.4 m. As a main result, the effect of a change in ceiling height on the maximum HRR was examined by a comparison with a previous study that considered the change in ceiling height. In addition, a more generalized correlation equation was proposed that could predict the maximum HRR in closed-compartments regardless of the changes in ceiling height. This correlation had an average error of 7% and a maximum error of 19% for various fire growth rates when compared with the numerical results. Finally, the applicability of the proposed correlation to representative fire compartments applied to the domestic performance-based design (PBD) was examined. These results are expected to provide useful information on predicting the maximum HRR caused by flashover in closed-compartments as well as the input information required in a fire simulation.

Optimization Study for the Production of 6-Shogaol-rich Ginger (Zingiber officinale Roscoe) under Conditions of Mild Pressure and High Temperature (가압조건에서 생강 유래 6-shogaol 변환을 위한 가열 조건 최적화)

  • Park, Ho-Young;Ha, Sang Keun;Choi, Jiwon;Choi, Hee-Don;Kim, Yoonsook;Park, Yongkon
    • Korean Journal of Food Science and Technology
    • /
    • v.46 no.5
    • /
    • pp.588-592
    • /
    • 2014
  • Under optimized condition mild pressure in combination with specific temperature for heat treatment transform the 6-gingerol into 6-shogaol. The purpose of this study was to optimize the conditions used for heat treatment under pressure for increasing 6-shogaol content in ginger (Zingiber officinale Roscoe). A central composite experimental design was used to evaluate the effects of application temperature ($70-130^{\circ}C$) and temperature-holding time (95-265 min) on the transformation of 6-shogaol. The experimental values were shown to be in significantly good agreement with the predicted values (adjusted determination coefficient, $R^2{_{Adj}}=0.9857$). 6-Shogaol content increased as the application temperature and temperature-holding time increased. By analyzing the response surface plots, the optimum conditions of heat treatment (temperature and time) for increasing 6-shogaol content were found to be $127^{\circ}C$ and 109 min, respectively. Under these optimal conditions, the predicted 6-shogaol content was 3.98 mg/g dried ginger. The adequacy of the model equation for predicting the optimum response values was effectively verified by the validation data.