• Nuclear Engineering and Technology
• /
• v.25 no.1
• /
• pp.91-101
• /
• 1993

To predict accurately the thermal hydraulic behavior of light water reactors during normal or abnormal operation, the accurate estimation of the void distribution is required. Up to date, many techniques for predicting void fraction of two-phase flow systems have been suggested. Among these techniques, the drift-flux model is widely used because of its exact calculation ability and simplicity. However, to get more accurate prediction of void fraction using drift-flux model, slip and flow regime effects must be considered more properly In the drift-flux method, these two effects are accounted for by two drift-flux parameters ; $C_{o}$ and (equation omitted). At earlier stage, $C_{o}$ is measured in a circular tube. In this study, $C_{o}$ is experimentally determined by measuring local void fraction and vapor velocity distribution in a rectangular subchannel having 4 heating rods which simulates nuclear subchannels. The measurements are peformed with two-electrical conductivity probes which are known to be adequate for measuring local parameters. The experiments are performed at low flow rate and the system pressure less than 3 atmo spheric pressure. In this experiment, (equation omitted), is not measured, but quoted from well-known empirical correlation to formulate $C_{o}$. Finally, $C_{o}$ is expressed as a function of channel averaged void fraction. fraction.

• 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.

• Journal of the Korean Geotechnical Society
• /
• v.39 no.1
• /
• pp.15-25
• /
• 2023

Recent weather changes have led to an increase in heavy rainfall resulting in frequent large-scale slope failures. To minimize damage to life and property, a measurement system is used in slope failure warning systems. However, understanding the slope failure behavior is difficult as the measurement system only measures a specific point. Therefore, numerical analysis must be p erformed with the measurement system. The soil water characteristic curve (SWCC) drying curve and boundary conditions that consider evapotranspiration and precipitation have been applied to numerical analysis, but the hysteresis of SWCC affects the numerical analysis results. To address this, a new evapotranspiration calculation method is proposed and applied to boundary conditions, and the measurement data are compared with the results of the numerical analysis. This method takes into account the different infiltration behaviors on evapotranspiration according to the drying and wetting curves of the SWCC, and allows for a more rational prediction of water movement on unsaturated slopes.

• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.

• Journal of Soil and Groundwater Environment
• /
• v.9 no.1
• /
• pp.79-86
• /
• 2004

Hydraulic conductivity along rock fracture is mainly dependent on fracture geometries such as orientation, aperture, roughness and connectivity. Therefore, it needs to consider fracture geometries sufficiently on a fracture model for a numerical analysis to calculate permeability coefficient in a fracture. This study performed new type of numerical analysis using a homogenization analysis method to calculate permeability coefficient accurately along single fractures with several fracture models that were considered fracture geometries as much as possible. First of all, fracture roughness and aperture variation due to normal stress applied on a fracture were directly measured under a confocal laser scaning microscope (CLSM). The acquired geometric data were used as input data to construct fracture models for the homogenization analysis (HA). Using the constructed fracture models, the homogenization analysis method can compute permeability coefficient with consideration of material properties both in microscale and in macroscale. The HA is a new type of perturbation theory developed to characterize the behavior of a micro inhomogeneous material with a periodic microstructure. It calculates micro scale permeability coefficient at homogeneous microscale, and then, computes a homogenized permeability coefficient (C-permeability coefficient) at macro scale. Therefore, it is possible to analyze accurate characteristics of permeability reflected with local effect of facture geometry. Several computations of the HA were conducted to prove validity of the HA results compared with the empirical equations of permeability in the previous studies using the constructed 2-D fracture models. The model can be classified into a parallel plate model that has fracture roughness and identical aperture along a fracture. According to the computation results, the conventional C-permeability coefficients have values in the range of the same order or difference of one order from the permeability coefficients calculated by an empirical equation. It means that the HA result is valid to calculate permeability coefficient along a fracture. However, it should be noted that C-permeability coefficient is more accurate result than the preexisting equations of permeability calculation, because the HA considers permeability characteristics of locally inhomogeneous fracture geometries and material properties both in microscale and macroscale.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.7 no.1
• /
• pp.861-876
• /
• 1965

During my stay in the Netherlands, I have studied the following, primarily in relation to the Mokpo Yong-san project which had been studied by the NEDECO for a feasibility report. 1. Unit hydrograph at Naju There are many ways to make unit hydrograph, but I want explain here to make unit hydrograph from the- actual run of curve at Naju. A discharge curve made from one rain storm depends on rainfall intensity per houre After finriing hydrograph every two hours, we will get two-hour unit hydrograph to devide each ordinate of the two-hour hydrograph by the rainfall intensity. I have used one storm from June 24 to June 26, 1963, recording a rainfall intensity of average 9. 4 mm per hour for 12 hours. If several rain gage stations had already been established in the catchment area. above Naju prior to this storm, I could have gathered accurate data on rainfall intensity throughout the catchment area. As it was, I used I the automatic rain gage record of the Mokpo I moteorological station to determine the rainfall lntensity. In order. to develop the unit ~Ydrograph at Naju, I subtracted the basic flow from the total runoff flow. I also tried to keed the difference between the calculated discharge amount and the measured discharge less than 1O~ The discharge period. of an unit graph depends on the length of the catchment area. 2. Determination of sluice dimension Acoording to principles of design presently used in our country, a one-day storm with a frequency of 20 years must be discharged in 8 hours. These design criteria are not adequate, and several dams have washed out in the past years. The design of the spillway and sluice dimensions must be based on the maximun peak discharge flowing into the reservoir to avoid crop and structure damages. The total flow into the reservoir is the summation of flow described by the Mokpo hydrograph, the basic flow from all the catchment areas and the rainfall on the reservoir area. To calculate the amount of water discharged through the sluiceCper half hour), the average head during that interval must be known. This can be calculated from the known water level outside the sluiceCdetermined by the tide) and from an estimated water level inside the reservoir at the end of each time interval. The total amount of water discharged through the sluice can be calculated from this average head, the time interval and the cross-sectional area of' the sluice. From the inflow into the .reservoir and the outflow through the sluice gates I calculated the change in the volume of water stored in the reservoir at half-hour intervals. From the stored volume of water and the known storage capacity of the reservoir, I was able to calculate the water level in the reservoir. The Calculated water level in the reservoir must be the same as the estimated water level. Mean stand tide will be adequate to use for determining the sluice dimension because spring tide is worse case and neap tide is best condition for the I result of the calculatio 3. Tidal computation for determination of the closure curve. During the construction of a dam, whether by building up of a succession of horizontael layers or by building in from both sides, the velocity of the water flowinii through the closing gapwill increase, because of the gradual decrease in the cross sectional area of the gap. 1 calculated the . velocities in the closing gap during flood and ebb for the first mentioned method of construction until the cross-sectional area has been reduced to about 25% of the original area, the change in tidal movement within the reservoir being negligible. Up to that point, the increase of the velocity is more or less hyperbolic. During the closing of the last 25 % of the gap, less water can flow out of the reservoir. This causes a rise of the mean water level of the reservoir. The difference in hydraulic head is then no longer negligible and must be taken into account. When, during the course of construction. the submerged weir become a free weir the critical flow occurs. The critical flow is that point, during either ebb or flood, at which the velocity reaches a maximum. When the dam is raised further. the velocity decreases because of the decrease\ulcorner in the height of the water above the weir. The calculation of the currents and velocities for a stage in the closure of the final gap is done in the following manner; Using an average tide with a neglible daily quantity, I estimated the water level on the pustream side of. the dam (inner water level). I determined the current through the gap for each hour by multiplying the storage area by the increment of the rise in water level. The velocity at a given moment can be determined from the calcalated current in m3/sec, and the cross-sectional area at that moment. At the same time from the difference between inner water level and tidal level (outer water level) the velocity can be calculated with the formula $h= \frac{V^2}{2g}$ and must be equal to the velocity detertnined from the current. If there is a difference in velocity, a new estimate of the inner water level must be made and entire procedure should be repeated. When the higher water level is equal to or more than 2/3 times the difference between the lower water level and the crest of the dam, we speak of a "free weir." The flow over the weir is then dependent upon the higher water level and not on the difference between high and low water levels. When the weir is "submerged", that is, the higher water level is less than 2/3 times the difference between the lower water and the crest of the dam, the difference between the high and low levels being decisive. The free weir normally occurs first during ebb, and is due to. the fact that mean level in the estuary is higher than the mean level of . the tide in building dams with barges the maximum velocity in the closing gap may not be more than 3m/sec. As the maximum velocities are higher than this limit we must use other construction methods in closing the gap. This can be done by dump-cars from each side or by using a cable way.e or by using a cable way.

• Journal of Korea Water Resources Association
• /
• v.54 no.6
• /
• pp.419-431
• /
• 2021

In Seoul, it has been confirmed that the duration of rainfall is shortened and the frequency and intensity of heavy rains are increasing with a changing climate. In addition, due to high population density and urbanization in most areas, floods frequently occur in flood-prone areas for the increase in impermeable areas. Furthermore, the Seoul City is pursuing various projects such as structural and non-structural measures to resolve flood-prone areas. A disaster prevention performance target was set in consideration of the climate change impact of future precipitation, and this study conducted to reduce the overall flood damage in Seoul for the long-term. In this study, 29 GCMs with RCP4.5 and RCP8.5 scenarios were used for spatial and temporal disaggregation, and we also considered for 3 research periods, which is short-term (2006-2040, P1), mid-term (2041-2070, P2), and long-term (2071-2100, P3), respectively. For spatial downscaling, daily data of GCM was processed through Quantile Mapping based on the rainfall of the Seoul station managed by the Korea Meteorological Administration and for temporal downscaling, daily data were downscaled to hourly data through k-nearest neighbor resampling and nonparametric temporal detailing techniques using genetic algorithms. Through temporal downscaling, 100 detailed scenarios were calculated for each GCM scenario, and the IDF curve was calculated based on a total of 2,900 detailed scenarios, and by averaging this, the change in the future extreme rainfall was calculated. As a result, it was confirmed that the probability of rainfall for a duration of 100 years and a duration of 1 hour increased by 8 to 16% in the RCP4.5 scenario, and increased by 7 to 26% in the RCP8.5 scenario. Based on the results of this study, the amount of rainfall designed to prepare for future climate change in Seoul was estimated and if can be used to establish purpose-wise water related disaster prevention policies.