• Nuclear Engineering and Technology
• /
• v.55 no.5
• /
• pp.1802-1813
• /
• 2023

Conjugate heat transfer between liquid metal and solid is a common phenomenon in a liquid-metal-cooled fast reactor's fuel assembly and heat exchanger, dramatically affecting the reactor's safety and economy. Therefore, comprehensively studying the sophisticated conjugate heat transfer in a liquid-metal-cooled fast reactor is profound. However, it has been evidenced that the traditional Simple Gradient Diffusion Hypothesis (SGDH), assuming a constant turbulent Prandtl number (Pr_t,, usually 0.85 - 1.0), is inappropriate in the Computational Fluid Dynamics (CFD) simulations of liquid metal. In recent decades, numerous studies have been performed on the four-equation model, which is expected to improve the precision of liquid metal's CFD simulations but has not been introduced into the conjugate heat transfer calculation between liquid metal and solid. Consequently, a four-equation model, consisting of the Abe k - ε turbulence model and the Manservisi k_𝜃 - ε_𝜃 heat transfer model, is applied to study the conjugate heat transfer concerning liquid metal in the present work. To verify the numerical validity of the four-equation model used in the conjugate heat transfer simulations, we reproduce Johnson's experiments of the liquid lead-bismuth-cooled turbulent pipe flow using the four-equation model and the traditional SGDH model. The simulation results obtained with different models are compared with the available experimental data, revealing that the relative errors of the local Nusselt number and mean heat transfer coefficient obtained with the four-equation model are considerably reduced compared with the SGDH model. Then, the thermal-hydraulic characteristics of liquid metal turbulent pipe flow obtained with the four-equation model are analyzed. Moreover, the impact of the turbulence model used in the four-equation model on overall simulation performance is investigated. At last, the effectiveness of the four-equation model in the CFD simulations of liquid sodium conjugate heat transfer is assessed. This paper mainly proves that it is feasible to use the four-equation model in the study of liquid metal conjugate heat transfer and provides a reference for the research of conjugate heat transfer in a liquid-metal-cooled fast reactor.

• Nuclear Engineering and Technology
• /
• v.55 no.10
• /
• pp.3874-3897
• /
• 2023

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High-Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next-generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds-averaged Navier-Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.13 no.3
• /
• pp.252-259
• /
• 2008

The deep-sea mining system is generally composed of surface vessel, lifting system, buffer, flexible pipe and miner. The mining system can be regarded as a large-scale system in which each subsystem is interconnected to other ones. In order to control a large-scale system, decentralized control approaches have been proposed recently. In this paper, as a basic study on application of decentralized control, firstly, the mining system was modeled in a simplified way. Lifting system and buffer were regarded as a spherical pendulum and the flexible pipe was taken as a two-dimensional linear spring connection. Based on the simplified model dynamics, the mining system can be decentralized two subsystems, the one consisting of surface vessel, lifting system and buffer, and the other, the miner. Next, this paper proposed the design of controller for each decentralized subsystem by regarding the interacting terms as disturbances. The controllers kept the constant distance between two subsystems during the miner was moving on the specified track. Finally, the efficiency of proposed controller was proven through the numerical simulation of the derived model.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.29 no.11B
• /
• pp.980-990
• /
• 2004

In this paper, we propose the dynamic bandwidth allocation algorithm for the IEEE802.3ah Ethernet Passive Optical Network(EPON) system to provide the fairness among terminals, and evaluate the delay-throughput performance by simulation. For the conventional EPON systems, an Optical Line Termination (OLT) schedules the upstream bandwidth for each Optical Network Unit (ONU), based on its buffer state. This scheme can provide a fair bandwidth allocation for each ONU. However, it has a critical problem that it does not guarantee the fair bandwidth among terminals which are connected to ONUs. For an example, we assume that the traffic from a greedy terminal increases at a time. Then, the buffer state of its ONU is instantly reported to the OLT, and finally the OW can get more bandwidth. As a result, the less bandwidth is allocated to the other ONUs, and thus the transfer delay of terminals connected to the ONUs gets inevitably increased. Noting that this unfairness problem exists in the conventional EPON systems, we propose a fair bandwidth allocation scheme by OLT with considering the buffer state of ONU as welt as the number of terminals connected it. For the performance evaluation, we develop the EPON simulation model with SIMULA simulation language. From the result of the throughput-delay performance and the dynamics of buffer state along time for each terminal and ONU, respectively, one can see that the proposed scheme can provide the fairness among not ONUs but terminals. Finally, it is worthwhile to note that the proposed scheme for the public EPON systems might be an attractive solution for providing the fairness among subscriber terminals.

• Journal of Korean Society of Transportation
• /
• v.20 no.5
• /
• pp.113-130
• /
• 2002

In this paper, we developed DNL(Dynamic Network Loading) model based on Moving cell theory to analyze the dynamic characteristics of traffic flow in congested network. In this paper vehicles entered into link at same interval would construct one cell, and the cells moved according to Cell following rule. In the past researches relating to DNL model a continuous single link is separated into two sections such as running section and queuing section to describe physical queue so that various dynamic states generated in real link are only simplified by running and queuing state. However, the approach has some difficulties in simulating various dynamic flow characteristics. To overcome these problems, we present Moving cell theory which is developed by combining Car following theory and Lagrangian method mainly using for the analysis of air pollutants dispersion. In Moving cell theory platoons are represented by cells and each cell is processed by Cell following theory. This type of simulation model is firstly presented by Cremer et al(1999). However they did not develop merging and diverging model because their model was applied to basic freeway section. Moreover they set the number of vehicles which can be included in one cell in one interval so this formulation cant apply to signalized intersection in urban network. To solve these difficulties we develop new approach using Moving cell theory and simulate traffic flow dynamics continuously by movement and state transition of the cells. The developed model are played on simple network including merging and diverging section and it shows improved abilities to describe flow dynamics comparing past DNL models.

• Journal of Advanced Marine Engineering and Technology
• /
• v.40 no.3
• /
• pp.228-234
• /
• 2016

As a college dormitory has the features of high dwelling density and a floating population that becomes crowded during particular times, when a disaster such as a fire occurs, it has the risk of causing much loss of life. In this study, the fire simulation program Fire Dynamics Simulator (FDS), is used to predict the risk when a fire occurs, to analyze the problem, and to suggest an improvement plan for a new cadet dormitory at an university in Korea. The research results are as follows. When a fire occurred in the ironing room inside the cadet dormitory, a smoke detector operated after 65 seconds. Thirteen seconds later, a sprinkler started to operate. The temperature and carbon monoxide density reached the limit value at 241 and 248 seconds, respectively. Because the limit visibility value was reached within 66 seconds after the occurrence of a fire, it is predicted that preparation must be finished and evacuation should begin within 1 minute after the fire occurs, in order to have no casualties. Synthesizing this dormitory fire risk prediction result, the visibility value is considered to be the most dangerous factor for personal safety. Because of this, installing a smoke extraction system is suggested to secure visibility. After the installation of a smoke extraction system, the problem of smoke diffusion in the corridors improved.

• Ocean and Polar Research
• /
• v.35 no.3
• /
• pp.249-258
• /
• 2013

Vertical and horizontal mixing processes in the ocean mixed layer determine sea surface temperature and temperature variability. Accordingly, simulating these processes properly is crucial in order to obtain more accurate climate simulations and more reliable future projections using an ocean general circulation model (OGCM). In this study, by using Modular Ocean Model version 4 (MOM4) developed by Geophysical Fluid Dynamics Laboratory, the upper ocean temperature and mixed layer depth were simulated with two different vertical mixing schemes that are most widely used and then compared. The resultant differences were analyzed to understand the underlying mechanism, especially in the Tropical Pacific Ocean where the differences appeared to be the greatest. One of the schemes was the so-called KPP scheme that uses K-Profile parameterization with nonlocal vertical mixing and the other was the N scheme that was rather recently developed based on a second-order turbulence closure. In the equatorial Pacific, the N scheme simulates the mixed layer at a deeper level than the KPP scheme. One of the reasons is that the total vertical diffusivity coefficient simulated with the N scheme is ten times larger, at maximum, in the surface layer compared to the KPP scheme. Another reason is that the zonal current simulated with the N scheme peaks at a deeper ocean level than the KPP scheme, which indicates that the vertical shear was simulated on a larger scale by the N scheme and it enhanced the mixed layer depth. It is notable that while the N scheme simulates a deeper mixed layer in the equatorial Pacific compared to the KPP scheme, the sea surface temperature (SST) simulated with the N scheme was cooler in the central Pacific and warmer in the eastern Pacific. We postulated that the reason for this is that in the central Pacific atmospheric forcing plays an important role in determining SST and so does a strong upwelling in the eastern Pacific. In conclusion, what determines SST is crucial in interpreting the relationship between SST and mixed layer depth.

• KIEAE Journal
• /
• v.8 no.5
• /
• pp.49-54
• /
• 2008

The experimental of temperature, humidity and velocity was taken from the underground pit which utilized the system of ground heat source quite similar to the cool-pit system. Also, through CFD analysis, one could review the effectiveness of analysis of future alternatives. Furthermore, the temperature range of mock up cool-pit system was analyzed by inputting the weather data of annual average soil temperature provided by KMA(Korea Meteorological Administration) into the fluid simulation of anticipated heat distribution. Firstly, the difference between the temperature of air exhaust of the pit or the temperature of air supply of the compressor room and the experimental data for the month of May from the CFD analysis came out to be $0.6^{\circ}C$ and $0.9^{\circ}C$ respectively with tolerance of 3.1% and 4.7%. Secondly, the difference between the temperature of air exhaust of the Pit or the temperature of air supply of the compressor room and the experimental data for the month of July from the CFD analysis came out to be $0.8^{\circ}C$ and $1.1^{\circ}C$ respectively with tolerance of 3.3% and 4.5%. Thirdly, for the month of May, the difference between the experimental data taken for the air exhaust of the Pit or the air supply of the compressor room and soil temperature provided by KMA for monthly and yearly average temperature of Jeonju region came out be $1.9^{\circ}C$ and $1.8^{\circ}C$ respectively with tolerance of 10.7% and 9.8%. Fourthly, for the month of July, the difference between the experimental data taken for the air exhaust of the Pit or the air supply of the compressor room and soil temperature provided by KMA for monthly and yearly average temperature of Jeonju region came out be $1.1^{\circ}C$ and $1.4^{\circ}C$ respectively with tolerance of 4.5% and 5.8%. The result of above experiments allowed us to establish CFD model set up as a verification tool that is based on experimental data collected within the Pit area. Also, one could confirm the possibility to apply weather data of soil temperature provided by KMA in order to anticipate proper value for CFD analysis.

• Journal of the Korean Geotechnical Society
• /
• v.36 no.11
• /
• pp.149-156
• /
• 2020

Permeation grouting effectively enhances soil strength and decreases permeability of soil; however, the flow of grout is heavily affected by anisotropy of hydraulic conductivity in layers. Therefore, this study investigates the effect of permeability anisotropy on the effective radius of horizontal permeation grout using computational fluid dynamics (CFD). We modeled the horizontal permeation grout flow as a two-phase viscous fluid flow in porous media, and the model incorporated the chemical diffusion and the viscosity variation due to hardening. The numerical simulation reveals that the permeability anisotropy shapes the grout bulb to be elliptic and the dissolution-driven diffusion causes a gradual change in grout pore saturation at the edge of the grout bulb. For the grout pore saturations of 10%, 50% and 90%, the horizontal and vertical radii of grout bulb are estimated when the horizontal-to-vertical permeability ratio varies from 0.01 to 100, and the predictive model equations are suggested. This result contributes to more efficient design of injection strategy in formation layers with permeability anisotropy.

• Composites Research
• /
• v.30 no.5
• /
• pp.323-330
• /
• 2017

In this study, multiscale analysis in which the information obtained from molecular dynamics simulation is applied to the continuum mechanics level is conducted to investigate the effects of clustering of silicon carbide nanoparticles reinforced into polypropylene matrix on mechanical behavior of nanocomposites. The elastic behavior of polymer nanocomposites is observed for various states of nanoparticulate agglomeration according to the model reflecting the degradation of interphase properties. In addition, factors which mainly affect the mechanical behavior of the nanocomposites are identified, and new index 'clustering density' is defined. The correlation between the clustering density and the elastic modulus of nanocomposites is understood. As the clustering density increases, the interfacial effect decreased and finally the improvement of mechanical properties is suppressed. By considering the random distribution of the nanoparticles, the range of elastic modulus of nanocomposites for same value of clustering density can be investigated. The correlation can be expressed in the form of exponential function, and the mechanical behavior of the polymer nanocomposites can be effectively predicted by using the nanoparticulate clustering density.