• Transactions of the Korean hydrogen and new energy society
• /
• v.34 no.5
• /
• pp.439-446
• /
• 2023

In this study, the computational fluid dynamics (CFD) simulations were conducted to verify the cooling capacity of the cryocooler used for pre-cooling of hydrogen gas. Based on the experimental results, the effect of the flow rate on a copper pipe attached to the bottom of the cryocooler was investigated. In this study, the temperature data was calculated through the change of boundary condition for heat flux in the copper pipe. In addition, the cooling capacity of the cryocooler for pre-cooling hydrogen gas was considered by calculating the cooling temperature according to the flow rate in the certified operating range. Consequently the pre-cooing system for hydrogen gas was validated with a reasonable accuracy through CFD simulations.

• Journal of the Korean Geotechnical Society
• /
• v.25 no.4
• /
• pp.91-103
• /
• 2009

Slope failure depends upon the climatic features related to related rainfall, structural geology and geomorphological features as well as the variation of the mechanical behaviors of soil constituting a slope. In this paper, among many variables, effects of soil layer thickness on the slope failure process, and variations of matric suction and volumetric water content were observed. When the soil layer is relatively thick, the descending wetting front decreases matric suction and the observed matric suction reaches to "0" value. When the wetting front reaches to the impermeable boundary, the bottom surface of steel soil box, ascending wetting front was observed. This observation can be postulated to be the effects of various sizes of pores. When macro size pores exist, the capillary effects can be reduced and infilling of pore will be limited. The partially filled pores would be filled with water during the ascending of the wetting front, which bounces from the impermeable boundary. This assumption has been assured from the observation of variation of the volumetric water contents at different depth. When the soil layer is thick (thickness = 20 cm), for granular material, erosion is a cause triggering the slope failure. It has been found that the initiation of erosion occurs when the top soil is fully saturated. Meanwhile, when the soil layer is shallow (thickness = 10 cm), slope slides as en mass. The slope failure for this condition occurs when the wetting front reaches to the interface between the soil layer and steel soil box. As the wetting front approaches to the bottom of soil layer, reduction of shear resistance along the boundary and increase of the unit weight due to the infiltration occur and these produce complex effects on the slope failure processes.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.355-365
• /
• 1994

An analytical model of return flow is presented outside the surf zone. The governing equation is derived from the Navier-Stokes equation and the continuity. Each term of the governing equation is evaluated by the ordering analysis. Then the infinitesimal terms, i.e. the turbulent normal stress, the squared vertical velocity of water particle and the streaming velocity, are neglected. The driving forces of return flow are calculated using the linear wave theory for the shallow water approximation. Especially, the space derivative of local wave heights is described considering a shoaling coefficient. The vertical distribution of eddy viscosity is discussed to the customary types which are the constant, the linear function and the exponential function. Each coefficient of the eddy viscosities which sensitively affect the precision of solutions is uniquely decided from the additional boundary condition which the velocity becomes zero at the wave trough level. Also the boundary conditions at the bottom and the continuity relation are used in the integration of the governing equation. The theoretical solutions of present model are compared with the various experimental results. The solutions show a good agreement with the experimental results in the case of constant or exponential function type eddy viscosity.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.15 no.4
• /
• pp.391-402
• /
• 2017

The Korea Radioactive Waste Agency (KORAD) has developed a concrete cask for the dry storage of spent nuclear fuel that has been generated by domestic light-water reactors. During long-term storage of spent nuclear fuel in concrete casks kept in dry conditions, the integrity of the concrete cask and spent nuclear fuel must be maintained. In addition, the radiation dose rate must not exceed the storage facility's design standards. A suitable shielding design for radiation protection must be in place for the dry storage facilities of spent nuclear fuel under normal and accident conditions. Evaluation results show that the appropriate distance to the annual dose rate of 0.25 mSv for ordinary citizens is approximately 230 m. For a $2{\times}10$ arrangement within storage facilities, rollover accidents are assumed to have occurred while transferring one additional storage cask, with the bottom of the cask facing the controlled area boundary. The dose rates of 12.81 and 1.28 mSv were calculated at 100 m and 230 m from the outermost cask in the $2{\times}10$ arrangement. Therefore, a spent nuclear fuel concrete cask and storage facilities maintain radiological safety if the distance to the appropriately assessed controlled area boundary is ensured. In the future, the results of this study will be useful for the design and operation of nuclear power plant on-site storage or intermediate storage facilities based on the spent fuel management strategy.

• Journal of Korea Soil Environment Society
• /
• v.4 no.1
• /
• pp.109-118
• /
• 1999

Recently, transport parameters of conservative solutes such as KCl in a porous medium have been successfully determined using time domain reflectometry (TDR) . This study was initiated to Investigate the applicability of TDR technique to monitoring the fate of a heavy metal ion in a sandy soil and the distribution of its concentration along travel distance with time. A column test was conducted in a laboratory that consists of monitoring both resident and flux concentrations of $ZnCl_2$in a sandy soil under a breakthrough condition. A tracer of $ZnCl_2$(10 g/L) was injected onto the top surface of the sample as pulse type as soon as a steady-state condition was achieved. Time-series measurements of resistance and electrical conductivity were performed at 10 cm and 20 cm of distances from the inlet boundary by horizontal-positioning of parallel TDR metallic rods and using an EC-meter for the effluent exiting the bottom boundary respectively. In addition. Zn ions of the effluent were analyzed by ICP-AES. Since the mode and position of concentration detected by TDR and effluent were different, comparison between ICP analysis and TDR-detected concentration was made by predicting flux concentration using CDE model accommodating a decay constant with the transport parameters obtained from the resident concentrations. The experimental results showed that the resident concentration resulted in earlier and higher peak than the flux concentration obtained by EC-meter, implying the homogeneity of the packed sandy soil. A close agreement was found between the predicted from the transport parameters obtained by TDR and the measured $ZnCl_2$concentration. This indicates that TDR technique can also be applied to monitoring heavy metal concentrations in the soil once that a decay constant is obtained for a given soil.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.33 no.1
• /
• pp.9-26
• /
• 1997

Stokes drift(SD) and Lagrangian discharge(LD) are important factors for analysis of flushing time, tidal exchange, solute transport and pollutant dispersion. The factors should be calculated using the approached method to flow phenomena. The aim of this paper re-examines the previous procedures for computing the SD and LD, and is to propose the new method approached to stratified flow field in the cross-section of coastal region, e.g. Masan Bay. The intensity of velocity near the bottom boundary layer(BBL) depends on the sea-bed irregularity in the coastal estuaries. So we calculated the depth mean velocity(DMV) considering that of BBL omitted in Kjerfve's calculation method. It revealed that BBL effect resulting in application of the bay acts largely on DMV in half more among 1l stations. The new expression of SD and LD per unit width in the cross-section using the developed DMV and proposed decomposition procedure of current were derived as follow : $$Q=u_0+\frac{1}{2}H_1{U_1cos(\varphi_h-\varphi_u)+U_3cos(\varphi_h-\varphi{ud})} LD ED SD$(Q_{skim}+Q_{sk2}) The third term, $Q_{sk2}$, on the right-hand of the equation is showed newly and arise from vertical oscillatory shear. According to the results applied in 3 cross-sections including 11 stations of the bay, the volume difference between proposed and previous SD was founded to be almost 2 times more at some stations. But their mean transport volumes over all stations are 18% less than the previous SD. Among two terms of SD, the flux of second term, $Q_{skim}$, is larger than third term, $Q_{sk2}$, in the main channel of cross-section, so that $Q_{skim}$ has a strong dependence on the tidal pumping, whereas third term is larger than second in the marginal channel. It means that $Q_{sk2}$ has trapping or shear effect more than tidal pumping phenomena. Maximum range of the fluctuation in LD is 40% as compared with the previous equations, but mean range of it is showed 11% at all stations, namely, small change. It mean that two components of SD interact as compensating flow. Therefore, the computation of SD and LD depend on decomposition procedure of velocity component in obtaining the volume transport of temporal and spacial flow through channels. The calculation of SD and LD proposed here can separate the shear effect from the previous SD component, so can be applied to non-uniform flow condition of cross-section, namely, baroclinic flow field.

• Journal of Korea Water Resources Association
• /
• v.49 no.1
• /
• pp.61-72
• /
• 2016

In Korea, the aquifers at the coastal areas are mostly shallow alluvial unconfined aquifers. To simulate the flow and dispersion in unconfined aquifer, a FDM model has been developed to solve the nonlinear Boussinesq equation. Related analysis and verification have been executed. The iteration method is used to solve the nonlinearity, and the model shows 3-D shape because it is a 2-D y model that consider the undulation of water table and bottom. For the verification of the model, the output of flow module is compared to the 1-D analytic solution of Lee (1989) which have the drawdown or uplift boundary condition, and the two results show almost the same value. and the mass balance of dispersion module shows about 10% error. The developed model can be used for the analysis and design of the flow and dispersion in the unconfined aquifers. The model has been applied to the estuary area of Ssangcheon watershed, and the parameters have been deduced as a result : hydraulic conductivity is 90 m/day, and longitudinal dispersivity is 15 m. And the analysis with these parameters shows that the wells are situated in the influence circle of each others except for No. 7 well. Groundwater discharge to sea is $3700m^3/day$. And the chlorine ion ($cl^-$) concentration at the pumping wells increase at least 1000 mg/L if groundwater dam is not exist, so the groundwater dam plays an important role for the prevention of sea water intrusion.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.9
• /
• pp.69-76
• /
• 2018

The present study analyzed the effect of film hole position of 45 degree ribbed cooling channel on film cooling performance of gas turbine blades. We also investigated the influence of the ribs under the fixed blowing ratio. Three-dimensional numerical model was constructed and extensive simulation was conducted using the commercial code (Fluent ver. 17.0) under steady-state condition. Base on the simulation results, We investigated the cooling effectiveness, flow velocity, streamline, and pressure coefficient. Moreover, We analyzed the effect of cooling hole position on ejection of the secondary flow caused by the rib structure. From the results, It was found that internal flow of the cooling channel forms a vortex pair in the counterclockwise from the top side, and clockwise from the bottom side. For the channels with ribs, the vortex flow generated by the ribs caused a higher pressure difference near the hole outlet, resulting in at least 12% higher cooling effectiveness than the channel without ribs. Additionally, when the hole is located on the left side of the ribbed channel (Rib-Left), it can be found that the secondary flow generated by the ribs hits against wall surface near the hole to form a flow in the direction of the hole inclination angle. Therefore, It is considered that the region where the cooling gas discharged to the blade surface stays in the main flow boundary layer is wider than the other cases. In this case, The largest pressure coefficient difference was observed near the outlet of the hole, and as a result, the discharge of the cooling gas was accelerated and the cooling efficiency was slightly increased.

• The Journal of the Acoustical Society of Korea
• /
• v.35 no.4
• /
• pp.261-271
• /
• 2016

In the prediction of response of a pile in soil, numerical approaches such as a finite element method are generally applied due to complicate nonlinear behaviors of soils. However, the numerical methods based on the finite elements require heavy efforts in pile and soil modelling and also take long computing time. So their usage is limited especially in the early design stage in which principal dimensions and properties are not specified and tend to vary. On the contrary, theoretical approaches adopting linear approximations for soils are relatively simple and easy to model and take short computing time. Therefore, if they are validated to be reliable, they would be applicable in predicting responses of a pile in soil, particularly in early design stage. In case of wind turbines regarded in this study, it is required to assess their natural frequencies in early stages, and in this simulation the supporting pile inserted in soil could be replaced with a simplified elastic boundary condition at the bottom end of the wind turbine tower. To do this, analysis for a pile in soil is performed in this study to extract the spring constants at the top end of the pile. The pile in soil can be modelled as a beam on elastic spring by assuming that the soils deform within an elastic range. In this study, it is attempted to predict pile deformations and influence factors for lateral loads by means of the beam-on-spring model. As two example supporting structures for wind turbines, mono pile and suction pile models with different diameters are examined by evaluating their influence factors and validated by comparing them with those reported in literature. In addition, the deflection profiles along the depth and spring constants at the top end of the piles are compared to assess their supporting features.

• Economic and Environmental Geology
• /
• v.56 no.1
• /
• pp.55-63
• /
• 2023

Variations in geochemical and mineralogical properties of the ferromanganese(Fe-Mn) crust reflect environmental changes. In the present study, geochemical and mineralogical analyses, including micro X-ray fluorescence and X-ray diffraction, were utilized to reconstruct the paleo-ocean environment of western Pacific Magellan seamount cluster. Samples of the Fe-Mn crust were collected using an epibenthic sledge from the open seamount XX (151° 51.12' 7.2" E and 16° 8.16' 9.6" N, 1557 meters below sea level) in the Western Pacific Magellan Seamount. According to the structure and phosphating status, the Fe-Mn crust of the OSM-XX can be divided into the following: phosphatizated (L4-L5), massive non-phosphatizated (L3), and porous non-phosphatizated (L1-L2) portions. All ferromanganese layers contain vernadite, and owing to the presence of carbonate fluorapatite (CFA), the phosphatizated portion (L4-L5) is rich in Ca and P. The massive non-phosphatizated section (L3) contains high Mn, Ni, and Co, whereas the porous non-phosphatizated portion (L1-L2), which comprises detrital quartz and feldspar, is rich in Fe. Variations in properties of the Fe-Mn crust from the OSM-XX reflect changes in the nearby marine environment. The formation of this crust started at approximately 51.87 Ma, and precipitation of the CFA during the global phosphatization event that occurred at approximately 36-32 Ma highlights an elevated sea level and low temperature during the associated period. The high Mn, Ni, and Co concentrations and elevated Mn/Fe ratios of samples from the massive phosphatizated portion indicate that the oxygen minimum zone (OMZ) was enhanced, and reducing conditions prevailed during the crust formation. The high Fe and low Mn/Fe ratios in the porous portion indicate a weak OMZ and dominantly oxidizing conditions. These data reflect environmental changes following the end of the Mi-1 glacial period in the Miocene-Oligocene boundary. Subsequently, Mn/Fe and Co/Mn ratios increased slightly in the outermost part of Fe-Mn crust because of the enhanced bottom current and OMZ associated with the continued cooling from approximately 9 Ma. However, the reduced carbonate dissolution rate in the Pacific Ocean from approximately 6 Ma decreased the growth rate of the Fe-Mn crust.