• Nuclear Engineering and Technology
• /
• v.55 no.12
• /
• pp.4431-4446
• /
• 2023

A Multiphysics coupling framework, MPCORE, has been developed to analyze safety parameters using the best estimate codes. The framework contains neutron kinetics (NK), thermal hydraulics (TH), and fuel performance (FP) codes to analyze fuel burnup, radial power distribution, and coolant temperature (T_bc). Shuffling and rotation capabilities have been verified on the Watts Bar reactor for three cycles. This study focuses on two coupling approaches for TH and FP modules. The one-way coupling approach involves coupling the FP code with the NK code, providing no data to the TH modules but getting T_bc as boundary condition from TH module. The two-way coupling approach exchanges information from FP to TH modules, so that the simplified heat conduction solver of the TH module is not used. The power profile in both approaches does not differ significantly, but there is an impact on coolant and cladding parameters. The one-way coupling approach tends to over-predict the cladding hydrogen concentration (CHC). This research highlights the difference between one-way and two-way coupling on critical boron concentration, T_bc, CHC, oxide surface temperature, and pellet centerline temperature. Overall, MPCORE framework with two-way coupling provides a more accurate and reliable analysis of safety parameters for nuclear reactors.

• Nuclear Engineering and Technology
• /
• v.56 no.6
• /
• pp.2099-2112
• /
• 2024

This study examined the effects of CRUD on reflood heat transfer behaviors of nuclear fuel rods during a loss-of-coolant-accident (LOCA) in a pressurized water reactor using a best-estimate thermal-hydraulic analysis code. Changes in thermal properties and boiling heat transfer characteristics of the CRUD layer were extensively reviewed, and a set of correction factors to reflect the changes was implemented into the code. A heat structure layer reflecting the effects of CRUDs on the properties was added to the outer surface of the fuel cladding. Numerical simulations were conducted to examine the effects of CRUDs on reflood cooling of overheated fuel rods for representative separate and integral effect tests, FLECHT-SEASET and LOFT. In LOFT analysis, the average cladding temperature was increased due to the low thermal conductivity of CRUD during steady-state operation; however, in both analyses, the peak cladding temperature decreased, and the quenching time was reduced. Obtained results revealed that when the porous CRUD layer is deposited on the fuel cladding, two opposite effects appear. Low thermal conductivity of the CRUD layer always increases fuel temperature during normal operation; however, its hydrophilic porous structures may contribute to accelerated reflood cooling of fuel rods during a LOCA.

• Proceedings of the Korean Institute of Information and Commucation Sciences Conference
• /
• 2007.06a
• /
• pp.227-230
• /
• 2007

If we use the microwave of SAR, we can observe on the ocean in spite of bad weather, day and night time. Sea surface images on the ocean of SAR have a lot of information on the atmospheric phenomena related to surface wind vector. Information of wind speed which is extracted from SAR images is used variously. Wind direction data and sigma nought value are put in the CMOD which can extract wind information in order to estimate sea surface wind from SAR images. Wind spectrum which is extracted from SAR always presents opposed two points of $180^{\circ}$ because of applying to 2D-FFT. These ambiguities should be decided by position of land, wind direction or numerical model. Previously, we converted into sigma nought after extracting Digital Number from RadarSat-1 SAR using ENVI4.0, thus, it took a long time because every process was manual. Therefore, we converted sigma nought by matlab code after making matlab code. After that, we are extracting wind direction from sigma nought. Now, to decide wind direction needs further study because wind direction has $180^{\circ}$ ambiguity.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.24 no.2
• /
• pp.246-252
• /
• 2018

In this paper, the friction drag force of 3D submerged body is investigated by considering the surface roughness, the first grid height, and the Reynolds number using open CFD source code, OpenFOAM 4.0. A procedure for estimating drag components by CFD code is set up and suggested in this study. In the 3D submerged body, because of the form factor in the 3D computations, the friction resistance with the small roughness of $12{\mu}m$ obtains different result with the smooth wall. As the Reynolds number increased, the boundary layer becomes thinner and the fiction resistance tends to decrease. In the computations for the effect of y+, the friction resistance and wall shear stress are excessively predicted when the y+ value deviates from the log layer. This is presumably because the boundary layer becomes thicker and the turbulence energy is excessively predicted in the nose due to the increase in y+ value. As the roughness increases, the boundary layer becomes thicker and the turbulence kinetic energy on the surface increases. From this study, the drag estimation method, considering the roughness by numerical analysis for ships or offshore structures, can be provided by using the suggested the y+ value and surface roughness with wall function.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.5
• /
• pp.545-565
• /
• 2018

A numerical code is developed based on potential flow theory to investigate nonlinear sloshing in rectangular Liquefied Natural Gas (LNG) tanks under forced excitation. Using this code, internal free-surface elevation and sloshing loads on liquid tanks can be obtained both in time domain and frequency domain. In the mathematical model, acceleration potential is solved in the calculation of pressure on tanks and the artificial damping model is adopted to account for energy dissipation during sloshing. The Boundary Element Method (BEM) is used to solve boundary value problems of both velocity potential and acceleration potential. Numerical calculation results are compared with published results to determine the efficiency and accuracy of the numerical code. Sloshing properties in partially filled rectangular and membrane tank under translational and rotational excitations are investigated. It is found that sloshing under horizontal and rotational excitations share similar properties. The first resonant mode and excitation frequency are the dominant response frequencies. Resonant sloshing will be excited when vertical excitation lies in the instability region. For liquid tank under rotational excitation, sloshing responses including amplitude and phase are sensitive to the location of the center of rotation. Moreover, experimental tests were conducted to analyze viscous effects on sloshing and to validate the feasibility of artificial damping models. The results show that the artificial damping model with modifying wall boundary conditions has better applicability in simulating sloshing under different fill levels and excitations.

• Land and Housing Review
• /
• v.6 no.3
• /
• pp.147-152
• /
• 2015

Site coefficient and amplification factor of current domestic Seismic Design Code based on American Seismic Code, have no consideration for the domestic ground condition in which the base rock is normally placed within 30m form the surface. By previous studies, the measured spectral acceleration of the result of dynamic centrifugal test and analysis was smaller than the design spectral acceleration for the period over 1.5 sec. Accordingly, in this study structural analysis and design using dynamic centrifugal test result for pile foundation were achieved, and the quantity of concrete and reinforcement of wall frame was compared with each other. Comparison results of cost using KBC the design spectral acceleration of SC, SD site and SDS, the quantity of reinforcement using SDS for SD site was 17~23% smaller than using the design code SD site.

• Nuclear Engineering and Technology
• /
• v.26 no.3
• /
• pp.355-366
• /
• 1994

The USNRC issued a revised ECCS rule that allows the use of best estimate computer codes for safety analysis. The rule also requires an estimation of uncertainty in calculated system response when applying the best estimate computer codes. A practical realistic evaluation methodology to evaluate the ECCS performance that satisfies the requirements of the ECCS rule has been developed and this paper describes the application of new realistic evaluation methodology to large break LOCA for, the demonstration of the new methodology. The computer code RELAP5/MOD3/KAERI, which was improved from RELAP5/MOD3.1, was used as the best estimate code in the application. The uncertainty of the code was evaluated by assessing several separate and integral effect tests, and for the application to actual plant Kori 3 & 4 was selected as the reference plant. Response surfaces for blowdown and reflood PCTs were generated from the results of the sensitivity analyses and probability distribution functions were established by random sampling or Monte-Carlo method for each response surface. Final uncertainties were quantified at 95% probability level and safety margins for large break LOCA were discussed.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.11
• /
• pp.1094-1104
• /
• 2010

In the previous studies, the validation of numerical codes has been conducted based on the qualitative comparison of predicted ice shapes with experiments, which poses a significant limit on the systematic analysis of ice shapes due to the variation of meteorological conditions. In response to this, the numerical code has been quantitatively validated against available experiment for the ice accretion on cylinders and airfoils in the present study. Ice shapes accumulated on the bodies are systematically investigated with respect to various icing parameters. To this end, maximum thickness, heading direction and ice thickness are quantified and expressed in the polar coordinate system for the comparison with other numerical results. By applying the quantitative analysis, similar shapes are intuitively distinguished. The developed numerical code underestimates the ice accretion area and the ice thickness of lower surface. In order to improve the accuracy, further accurate aerodynamic solver is required for the water droplet trajectories.

• Journal of ILASS-Korea
• /
• v.22 no.4
• /
• pp.210-217
• /
• 2017

A gasoline direct injection engine has an intake air temperature can be lowered by the fuel vaporization in the combustion chamber increase the volume efficiency is high compression ratio. Therefore, study for injection rate and characteristics which influence mixture formation in combustion chamber is important. Movement of the injector needle has a direct effect on the injection of the fuel, such as formation of cavitation, the fuel injection rate, etc. Therefore, recent studies on the dynamic characteristics of the injector considering the movement of the needle have been reported, but it takes a lot of time and cost to experimentally confirm the movement of the needle inside the injector. In this study, AMESim, a commercial 1-D code, and Star-CCM+, a 3-D CFD code, were used to predict the dynamic performance of the injector with needle motion. In order to predict the movement of the needle under the high pressure, the result of the surface pressure distribution according to the movement of the needle was derived by using the morphing technique of flow analysis. In addition, we predicted the injection rate of the injector considering the movement of the needle in conjunction with the 1-D code. The injection rate of the injector was measured by the BOSCH's method and the results were similar to those of the simulation results. This method can predict the injection rate and injection characteristics and this result is expected to be used to predict the performance of gasoline direct injection engines with low cost and time in the future.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.409-418
• /
• 2014

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.