• Nuclear Engineering and Technology
• /
• v.56 no.9
• /
• pp.3835-3850
• /
• 2024

The high-temperature gas-cooled reactor (HTGR) with spherical fuel elements contains complex pebble flow. The flow behavior of pebbles is influenced by various factors, such as pebble density, friction coefficient, wall structure, and discharge port size. Using a GPU-DEM numerical model, the effects of the friction coefficient on the cyclic loading and unloading of pebbles in the full-scale HTR-PM are studied. Numerical simulations with up to 420,000 spherical pebbles are conducted. Four sets of friction coefficient values are determined for comparative analysis based on experimental measurements. Discharging speed, residence time, stress, porosity, and velocity distribution are quantitatively analyzed. In addition, a comparison with the CT-PFD experiment is carried out to validate the numerical model. The results show that near-wall retention phenomena are observed in the reactor core only when using large friction coefficients. However, using friction coefficient values closer to the measured experimental values, the pebble bed in HTR-PM exhibited good flow characteristics. Furthermore, the friction coefficient also influences the porosity and velocity distribution of the pebble bed, with lower friction coefficients resulting in lower overall stress in the bed. The discharge outlet's influence varies with different friction coefficient values. In summary, this study demonstrates that the value of the friction coefficient has a complex influence on the pebble flow in HTR-PM, which provides important insights for future numerical and experimental studies in this field.

• Korean Journal of Environmental Agriculture
• /
• v.27 no.2
• /
• pp.111-119
• /
• 2008

This study was conducted to evaluate water quality and pollutant loads on small agricultural watershed in Jeonbuk province. The EC level of investigated watershed ranged from 0.07 to 0.52 dS/m, BOD level ranged from 0.1 to 5.0 mg/L, and $COD_{Cr}$ level ranged from 0.6 to 17.7 mg/L. As above, contents of water quality indicators covered wide range, but each indicator was alike in mean content every other year. The contents of EC, $Ca^{2+},\;Mg^{2+},\;K^+\;and\;Na^+$ were decreased in rainy season, but the contents of BOD, $COD_{Cr},\;COD_{Mn}$, T-N and T-P were not greatly different as compared to dry season. And high content of SS showed substantial sediments near the surface flow out and influence on water system in rainy season. The pollutant loads measured in terminal of watershed were $9.6{\sim}757.9$ kg/day for BOD, $51.2{\sim}1418.5$ kg/day for T-N and $0.3{\sim}44.7$ kg/day for T-P. The pollutant loads of BOD, T-N and T-P in rainy season increased several times as compared to dry season. In rainy season, watershed with more than 30% in the proportion of paddy field to land showed relatively low discharge and pollutant loads in comparison to watershed with less than 30%. The discharge of watershed in rainy season increased 5.7times compared with the dry season in watershed with less than 30% in the proportion of paddy field to land, whereas was only 2.3times in watershed with more than 30%. The correlation coefficient($R^2$) of regression between discharge and pollutant loads of T-N were higher than those of BOD and T-P.

• Journal of Environmental Impact Assessment
• /
• v.26 no.2
• /
• pp.160-170
• /
• 2017

The purpose of this study is to perform clustering of the habitat types and to identify the characteristics of species in the habitat types using mammal data (70,562) of the 3rd National Ecosystem Survey conducted from 2006 to 2012. The 15 habitat types recorded in the field-paper of the 3rd National ecosystem survey were reclassified, which was followed by the statistical analysis of mammal habitat types. In the habitat types cluster analysis, non-hierarchical cluster analysis (k-means cluster analysis), hierarchical cluster analysis, and non-metric multidimensional scaling method were applied to 14 habitat types recorded more than 30 times. A total of 7 Orders, 16 Families, and 39 Species of mammals were identified in the 3rd National Ecosystem Survey collected nationwide. When 11 clusters were classified by habitat types, the simple structure index was the highest (ssi = 0.07). As a result of the similarities and hierarchies between habitat types suggested by the hierarchical clustering analysis, the residential areas were the most different habitat types for mammals; the next following type was a cluster together with rivers and coasts. The results of the non-metric multidimensional scaling analysis demonstrated that both Mus musculus and Rattus norvegicus restrictively appeared in a residential area, which is the most discriminating habitat type. Lutra lutra restrictively appeared in coastal and river areas. In summary, according to our results, the mammalian habitat can be divided into the following four types: (1) the forest type (using forest as the main habitat and migration route); (2) the river type (using water as the main habitat); (3) the residence habitat (living near residential area); and (4) the lowland type (consuming grain or seeds as the main feeding resource).

• Analytical Science and Technology
• /
• v.28 no.6
• /
• pp.409-416
• /
• 2015

ICP-AES has been used in many laboratories due to the advantages of wide calibration range and multi-element analysis, but it may give erroneous results and suffer from spectral interference due to the large number of emission lines associated with each element. In this study, certified reference materials (CRMs) and field samples were analyzed by ICP-AES and HG-AAS according to the official Korean testing method for soil pollution to investigate analytical problems. The applicability of HG-ICP-AES was also tested as an alternative method. HG-AAS showed good accuracies (90.8~106.3%) in all CRMs, while ICP-AES deviated from the desired range in CRMs with low arsenic and high Fe/Al. The accuracy in CRM030 was estimated as below 39% at the wavelength of 193.696 nm by ICP-AES. Significant partial overlaps and sloping background interferences were observed near to 193.696 nm with the presence of 50 mg/L Fe and Al. Most CRMs were quantified with few or no interferences of Fe and Al at 188.980 nm. ICP-AES properly assessed low and high level arsenic for field samples, at 188.980 nm and 193.696 nm, respectively. The importance of the choice of measurement wavelengths corresponding to relative arsenic level should be noted. Because interferences were affected by the sample matrix, operation conditions and instrument figures, the analysts were required to consider spectral interferences and compare the analytical performance of the recommended wavelengths. HG-ICP-AES was evaluated as a suitable alternative method for ICP-AES due to improvement of the detection limit, wide calibration ranges, and reduced spectral interferences by HG.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.3
• /
• pp.174-183
• /
• 2014

Seabed beneath and near coastal structures may undergo large excess pore water pressure composed of oscillatory and residual components in the case of long durations of high wave loading. This excess pore water pressure may reduce effective stress and, consequently, the seabed may liquefy. If liquefaction occurs in the seabed, the structure may sink, overturn, and eventually increase the failure potential. In this study, to evaluate the liquefaction potential on the seabed, numerical analysis was conducted using the expanded 2-dimensional numerical wave tank to account for an irregular wave field. In the condition of an irregular wave field, the dynamic wave pressure and water flow velocity acting on the seabed and the surface boundary of the composite breakwater structure were estimated. Simulation results were used as input data in a finite element computer program for elastoplastic seabed response. Simulations evaluated the time and spatial variations in excess pore water pressure, effective stress, and liquefaction potential in the seabed. Additionally, the deformation of the seabed and the displacement of the structure as a function of time were quantitatively evaluated. From the results of the analysis, the liquefaction potential at the seabed in front and rear of the composite breakwater was identified. Since the liquefied seabed particles have no resistance to force, scour potential could increase on the seabed. In addition, the strength decrease of the seabed due to the liquefaction can increase the structural motion and significantly influence the stability of the composite breakwater. Due to limitations of allowable paper length, the studied results were divided into two portions; (I) focusing on the dynamic response of structure, acceleration, deformation of seabed, and (II) focusing on the time variation in excess pore water pressure, liquefaction, effective stress path in the seabed. This paper corresponds to (II).

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.3
• /
• pp.160-173
• /
• 2014

Seabed beneath and near coastal structures may undergo large excess pore water pressure composed of oscillatory and residual components in the case of long durations of high wave loading. This excess pore water pressure may reduce effective stress and, consequently, the seabed may liquefy. If liquefaction occurs in the seabed, the structure may sink, overturn, and eventually increase the failure potential. In this study, to evaluate the liquefaction potential on the seabed, numerical analysis was conducted using the expanded 2-dimensional numerical wave tank to account for an irregular wave field. In the condition of an irregular wave field, the dynamic wave pressure and water flow velocity acting on the seabed and the surface boundary of the composite breakwater structure were estimated. Simulation results were used as input data in a finite element computer program for elastoplastic seabed response. Simulations evaluated the time and spatial variations in excess pore water pressure, effective stress, and liquefaction potential in the seabed. Additionally, the deformation of the seabed and the displacement of the structure as a function of time were quantitatively evaluated. From the results of the analysis, the liquefaction potential at the seabed in front and rear of the composite breakwater was identified. Since the liquefied seabed particles have no resistance to force, scour potential could increase on the seabed. In addition, the strength decrease of the seabed due to the liquefaction can increase the structural motion and significantly influence the stability of the composite breakwater. Due to limitations of allowable paper length, the studied results were divided into two portions; (I) focusing on the dynamic response of structure, acceleration, deformation of seabed, and (II) focusing on the time variation in excess pore water pressure, liquefaction, effective stress path in the seabed. This paper corresponds to (I).

• Journal of Mushroom
• /
• v.11 no.1
• /
• pp.36-40
• /
• 2013

Armillaria gallica was ecologically surveyed to investigate its relationship with Gastrodia elata in cultivation areas of Korea in 2012. In the observation made around October 17 in the area of Namyangju (Gyeonggi), Sangju (Gyeongbuk) and Gimcheon (Gyeongbuk) in Korea, the fruit bodies of A. gallica were consistently found near the cultivation areas of Gastrodia elata across these cultivation areas. Since the temperature and rainfall have been considered as important factors of fruiting of A. gallica, we checked the temperature and rainfall around two weeks ago of mushroom fruiting. The average temperature of all cultivation areas was $9.9{\sim}17.5^{\circ}C$ (the lowest temp.: $3{\sim}23^{\circ}C$, the highest temp.: $15{\sim}26^{\circ}C$) and the annual rainfall was 0.5~1.0mm on October 10 which can be compared to 2.5~4.5mm on October 17. Fruiting bodies were generally developed at weedy field in which G. elata were cultivated five years ago. The field was inclined 15 degrees, east-facing and consists of loam. The fruit bodies were formed on both sides of a small ditch between the Gastrodia elata cultivation fields and neighboring fields, and the areas were common weeds and moisture at all times. The fruiting bodies were growing above the soil rather than oak wood. Gregarious, but also occurs to 1-5 fruiting bodies occur sporadically. Results observed in detail the growing areas, some fruit bodies were developed along the black and thick rhizomorph formation.

• The Journal of Industrial Distribution & Business
• /
• v.1 no.1
• /
• pp.13-21
• /
• 2010

Since Korea establish diplomatic ties with China in 1992, korea and China have had rapid progress in most of field as politic, economy, society and culture through basing on cultural commonality and geographical adjacency. Especially, China is the biggest trading partner to korea, and also Korea is third-biggest trading country to China. They become strategic cooperating relation in 2008. Currently, in terms of international trade relation, WTO/DDA negotiation is proceeding in difficulty, but FTA has been growing and extending in the world, and the two country, china and korea, have been competitively trying wide and active FTA negotiation promotion. After Financial crisis in 1997, according to the requirement of local economic cooperation, China has shown the interest to several countries since the conclusion of FTA treaty with ASEAN in 2005. China also makes the active afford to conclude FTA with Korea. Last May 28th, this was mentioned in the meeting between president Lee and Premier Wen Jiabao, so it is anticipated that the negotiation for FTA will be started in the near future. There are many political suggestions and concerns in terms of way of negotiation korea would choose. Some economist said that "'Continuous FTA aimed at long-term protocol should be promoted between korea and China and negotiated includingly'" However, this research claims that commodity exchange, service, and investment areas should be included and it has to be comprehensive package settlement style in negotiation. This research has found out the characteristics of China's negotiation and implications through the China's existed FTA negotiation examples. Currently, China has taken Continuous or a phase-negotiation method to ASEAN, Pakistan, Chile and some other developing country and to advanced countries like New Zealand or Singapore, comprehensive package settlement method is used in FTA negotiation. In consider of the FTA negotiation between Korea and China, Korea has some problems in the commodity change area in agriculture maket's opening. While, for china, the issues would happen in service trade area, especially when encountering finance and communication industries are opened, China's economy could be exposed to some risk. In result, Korea should expand its negotiation range from commodity trade to service trade, in order to exchange both issues, then the negotiation will be concluded more easily. In other word, for FTA, korea should follow comprehensive package settlement way that is similar to New zealand and Singapore case. Through this kind of method, Korea can expect effect of creating trade, conversion of it and preoccupancy of service field in china's market against the advanced countries like Usa, Europe and Japan. Also, to have a successful FTA negotiation, korea should find out china's policy for FTA negotiation. With this information, korea will be able to suggest the way to make a profit. Systematic analysis and comparison about previous negotiation cases of china are needed before the negotiation begin.

• Journal of the Korean GEO-environmental Society
• /
• v.19 no.4
• /
• pp.5-16
• /
• 2018

In the current study, the engineering behaviour of prebored and precast steel pipe piles was examined from a series of full-scale field measurements by conducting static pile load tests, dynamic pile load tests (EOID and restrike tests) and Class-A and Class-C1 type numerical analysis. The study includes the pile load - settlement relations, allowable pile capacity and shear stress transfer mechanism. Compared to the allowable pile capacity obtained from the static pile load tests, the dynamic pile load tests and the numerical simulation showed surprisingly large variations. Overall among these the restrike tests displayed the best results, however the reliability of the predictions from the numerical analysis was lower than those estimated from the dynamic pile load tests. The allowable pile capacity obtained from the EOID tests and the restrike tests indicated 20.0%-181.0% (avg: 69.3%) and 48.2%-181.1% (avg: 92.1%) of the corresponding measured values from the static pile loading tests, respectively. Furthermore, the computed results from the Class-A type analysis showed the largest scatters (37.1%-210.5%, avg: 121.2%). In the EOID tests, a majority of the external load were carried by the end bearing pile capacity, however, similar skin friction and end bearing capacity in magnitude were mobilised in the restrike tests. The measured end bearing pile capacity from the restrike tests were smaller than was measured from the EOID tests. The present study has revealed that if the impact energy is not sufficient in a restrike test, the end bearing pile capacity most likely will be underestimated. The shear stresses computed from the numerical analysis deviated substantially from the measured pile force distributions. It can be concluded that the engineering behaviour of the pile is heavily affected if a slime layer exists near the pile tip, and that the smaller the stiffness of the slime and the thicker the slime, the greater the settlement of the pile.

• Proceedings of the Korea Association of Crystal Growth Conference
• /
• 1996.06a
• /
• pp.179-200
• /
• 1996

The intrinsic instabilities of fluid flow occurred in the melt of the Czochralski crystal growth system Czochralski method, asymmetric flow patterns and temperature profiles in the melt have been studied by many researchers. The idea that the non-symmetric structure of the growing equipment is responsible for the asymmetric profiles is usually accepted at the first time. However further researches revealed that some intrinsic instabilities not related to the non-symmetric equipment structure in the melt could also appear. Ristorcelli had pointed out that there are many possible causes of instabilities in the melt. The instabilities appears because of the coupling effects of fluid flow and temperature profiles in the melt. Among the instabilities, the B nard type instabilities with no or low crucible rotation rates are analyzed by the visualizing experiments using X-ray radiography and the 3-D numerical simulation in this study. The velocity profiles in the Silicon melt at different crucible rotation rates were measured using X-ray radiography method using tungsten tracers in the melt. The results showed that there exits two types of fluid flow mode. One is axisymmetric flow, the other is asymmetric flow. In the axisymmetric flow, the trajectory of the tracers show torus pattern. However, more exact measurement of the axisymmetrc case shows that this flow field has small non-axisymmetric components of the velocity. When fluid flow is asymmetric, the tracers show random motion from the fixed view point. On the other hand, when the observer rotates to the same velocity of the crucible, the trajectory of the tracer show a rotating motion, the center of the motion is not same the center of the melt. The temperature of a point in the melt were measured using thermocouples with different rotating rates. Measured temperatures oscillated. Such kind of oscillations are also measured by the other researchers. The behavior of temperature oscillations were quite different between at low rotations and at high rotations. Above experimental results means that the fluid flow and temperature profiles in the melt is not symmetric, and then the mode of the asymmetric is changed when rotation rates are changed. To compare with these experimental results, the fluid flow and temperature profiles at no rotation and 8 rpm of crucible rotation rates on the same size of crucible is calculated using a 3-dimensional numerical simulation. A finite different method is adopted for this simulation. 50×30×30 grids are used. The numerical simulation also showed that the velocity and flow profiles are changed when rotation rates change. Futhermore, the flow patterns and temperature profiles of both cases are not axisymmetric even though axisymmetric boundary conditions are used. Several cells appear at no rotation. The cells are formed by the unstable vertical temperature profiles (upper region is colder than lower part) beneath the free surface of the melt. When the temperature profile is combined with density difference (Rayleigh-B nard instability) or surface tension difference (Marangoni-B nard instability) on temperature, cell structures are naturally formed. Both sources of instabilities are coupled to the cell structures in the melt of the Czochralski process. With high rotation rates, the shape of the fluid field is changed to another type of asymmetric profile. Because of the velocity profile, isothermal lines on the plane vertical to the centerline change to elliptic. When the velocity profiles are plotted at the rotating view point, two vortices appear at the both sides of centerline. These vortices seem to be the main reason of the tracer behavior shown in the asymmetric velocity experiment. This profile is quite similar to the profiles created by the baroclinic instability on the rotating annulus. The temperature profiles obtained from the numerical calculations and Fourier transforms of it are quite similar to the results of the experiment. bove esults intend that at least two types of intrinsic instabilities can occur in the melt of Czochralski growing systems. Because the instabilities cause temperature fluctuations in the melt and near the crystal-melt interface, some defects may be generated by them. When the crucible size becomes large, the intensity of the instabilities should increase. Therefore, to produce large single crystals with good quality, the behavior of the intrinsic instabilities in the melt as well as the effects of the instabilities on the defects in the ingot should be studied. As one of the cause of the defects in the large diameter Silicon single crystal grown by the