• Magazine of the Korean Society of Agricultural Engineers
• /
• v.22 no.1
• /
• pp.53-66
• /
• 1980

The research work is concerned with the analytical and experimental studies on the heat transfer phenomenon around the underground concrete digester used for biogas production Systems. A mathematical and computational method was developed to estimate heat losses from underground cylindrical concrete digester used for biogas production systems. To test its feasibility and to evaluate thermal parameters of materials related, the method was applied to six physical model digesters. The cylindrical concrete digester was taken as a physical model, to which the model,atical model of heat balance can be applied. The mathematical model was transformed by means of finite element method and used to analyze temperature distribution with respect to several boundary conditions and design parameters. The design parameters of experimental digesters were selected as; three different sizes 40cm by 80cm, 80cm by 160cm and l00cm by 200cm in diameter and height; two different levels of insulation materials-plain concrete and vermiculite mixing in concrete; and two different types of installation-underground and half-exposed. In order to carry out a particular aim of this study, the liquid within the digester was substituted by water, and its temperature was controlled in five levels-35。 C, 30。 C, 25。 C, 20。C and 15。C; and the ambient air temperature and ground temperature were checked out of the system under natural winter climate conditions. The following results were drawn from the study. 1.The analytical method, by which the estimated values of temperature distribution around a cylindrical digester were obtained, was able to be generally accepted from the comparison of the estimated values with the measured. However, the difference between the estimated and measured temperature had a trend to be considerably increased when the ambient temperature was relatively low. This was mainly related variations of input parameters including the thermal conductivity of soil, applied to the numerical analysis. Consequently, the improvement of these input data for the simulated operation of the numerical analysis is expected as an approach to obtain better refined estimation. 2.The difference between estimated and measured heat losses was shown to have the similar trend to that of temperature distribution discussed above. 3.It was found that a map of isothermal lines drawn from the estimated temperature distribution was very useful for a general observation of the direction and rate of heat transfer within the boundary. From this analysis, it was interpreted that most of heat losses is passed through the triangular section bounded within 45 degrees toward the wall at the bottom edge of the digesten Therefore, any effective insulation should be considered within this region. 4.It was verified by experiment that heat loss per unit volume of liquid was reduced as the size of the digester became larger For instance, at the liquid temperature of 35˚ C, the heat loss per unit volume from the 0. 1m$^3$ digester was 1, 050 Kcal/hr m$^3$, while at for 1. 57m$^3$ digester was 150 Kcal/hr m$^3$. 5.In the light of insulation, the vermiculite concrete was consistently shown to be superior to the plain concrete. At the liquid temperature ranging from 15。 C to 350 C, the reduction of heat loss was ranged from 5% to 25% for the half-exposed digester, while from 10% to 28% for the fully underground digester. 6.In the comparison of heat loss between the half-exposed and underground digesters, the heat loss from the former was fr6m 1,6 to 2, 6 times as much as that from the latter. This leads to the evidence that the underground digester takes advantage of heat conservation during winter.

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

• Journal of the Korean Geotechnical Society
• /
• v.35 no.5
• /
• pp.5-19
• /
• 2019

Considering the natural phenomenon in which steep slopes ($65^{\circ}{\sim}85^{\circ}$) consisting of rock mass remain stable for decades, slopes steeper than 1:0.5 (the standard of slope angle for blast rock) may be applied in geotechnical conditions which are similar to those above at the design and initial construction stages. In the process of analysing the stability of a good to fair continuum rock slope that can be designed as a steep slope, a general method of estimating rock mass strength properties from design practice perspective was required. Practical and genealized engineering methods of determining the properties of a rock mass are important for a good continuum rock slope that can be designed as a steep slope. The Genealized Hoek-Brown (H-B) failure criterion and GSI (Geological Strength Index), which were revised and supplemented by Hoek et al. (2002), were assessed as rock mass characterization systems fully taking into account the effects of discontinuities, and were widely utilized as a method for calculating equivalent Mohr-Coulomb shear strength (balancing the areas) according to stress changes. The concept of calculating equivalent M-C shear strength according to the change of confining stress range was proposed, and on a slope, the equivalent shear strength changes sensitively with changes in the maximum confining stress (${{\sigma}^{\prime}}_{3max}$ or normal stress), making it difficult to use it in practical design. In this study, the method of estimating the strength properties (an iso-angle division method) that can be applied universally within the maximum confining stress range for a good to fair continuum rock mass slope is proposed by applying the H-B failure criterion. In order to assess the validity and applicability of the proposed method of estimating the shear strength (A), the rock slope, which is a study object, was selected as the type of rock (igneous, metamorphic, sedimentary) on the steep slope near the existing working design site. It is compared and analyzed with the equivalent M-C shear strength (balancing the areas) proposed by Hoek. The equivalent M-C shear strength of the balancing the areas method and iso-angle division method was estimated using the RocLab program (geotechnical properties calculation software based on the H-B failure criterion (2002)) by using the basic data of the laboratory rock triaxial compression test at the existing working design site and the face mapping of discontinuities on the rock slope of study area. The calculated equivalent M-C shear strength of the balancing the areas method was interlinked to show very large or small cohesion and internal friction angles (generally, greater than $45^{\circ}$). The equivalent M-C shear strength of the iso-angle division is in-between the equivalent M-C shear properties of the balancing the areas, and the internal friction angles show a range of $30^{\circ}$ to $42^{\circ}$. We compared and analyzed the shear strength (A) of the iso-angle division method at the study area with the shear strength (B) of the existing working design site with similar or the same grade RMR each other. The application of the proposed iso-angle division method was indirectly evaluated through the results of the stability analysis (limit equilibrium analysis and finite element analysis) applied with these the strength properties. The difference between A and B of the shear strength is about 10%. LEM results (in wet condition) showed that Fs (A) = 14.08~58.22 (average 32.9) and Fs (B) = 18.39~60.04 (average 32.2), which were similar in accordance with the same rock types. As a result of FEM, displacement (A) = 0.13~0.65 mm (average 0.27 mm) and displacement (B) = 0.14~1.07 mm (average 0.37 mm). Using the GSI and Hoek-Brown failure criterion, the significant result could be identified in the application evaluation. Therefore, the strength properties of rock mass estimated by the iso-angle division method could be applied with practical shear strength.

• Journal of Korean Society of Steel Construction
• /
• v.25 no.1
• /
• pp.61-70
• /
• 2013

A turn-buckle is capable of adjusting the tensile force by left-hand threads and right-hand threads between tension members. There are different types of turn-buckles according to tension member and connection form but the practical and existing turn-buckles are incapable of measuring the tensile force. A turn-buckle for adjusting and measuring tensile force has therefore been developed. This study shows the ultimate strength and reliability for measurement of the new turn-buckles through finite element analysis of the developed ones. From analytic results of the new turn-buckles which have the measurement limit loads of 100kN, 200kN and 300kN, the ultimate strength is approximately five times stronger than the measurement limit capacity. Additionally, a review of the new turn-buckle, which has the measurement limit load of over 300kN, shows that there is a tendency for the size of turn-buckle to become larger. So the connection devices were designed and the loading test was conducted from the concept that the parallel connection of turn-buckle with 300kN capacity can measure the tensile force of 600kN. The results of parallel loading test show the sufficient possibility. Furthermore, the mock-up test was constructed to investigate the release of initial load and corrosion when the new turn-buckle is installed at the outdoor and exposed to rain and atmosphere.

• Composites Research
• /
• v.26 no.4
• /
• pp.207-212
• /
• 2013

In this study, the thickness optimization for uni-directional (UD) glass fiber reinforced polymer (GFRP) laminates of the spar cap of composite tidal blades was performed under the tip deflection constrains. The spar cap was composed of GFRP composites and carbon fiber reinforced polymer (CFRP) composites. The stress distributions in the blade as well as its material costs for the optimized results were additionally investigated. The optimized thickness was obtained by interacting a sequential quadratic programming (SQP) algorithm and an ABAQUS software to calculate an objective function. It was confirmed that the thickness of UD GFRP increased with a decrease of the restrained tip deflection when a thickness of UD CFRP laminates was constrained to 9 mm. The weight of the optimized spar-cap increased up to 96.2% while the maximum longitudinal tensile stress decreased up to 24.6%. The thickness of UD GFRP laminates increased with a decrease of the thickness of UD CFRP laminates when the tip deflection was constrained to 126.83 mm. The weight increased up to 40.1%, but the material cost decreased up to 16.97%. Finally, the relationships among the weight, internal tensile stress, and material costs were presented based on the optimized thicknesses of the spar cap.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.23 no.6
• /
• pp.265-271
• /
• 2013

Czochralski (CZ) growth process is one of the most important techniques for growing high quality sapphire single crystal for LED application. In this study, the inductively-heated CZ growth processes for the sapphire crystal of 300 mm length have been analyzed numerically using finite element method. The hot zone structures were modified with the crucible geometry change and the additional insulation layer installed above the crucible. The results show that the solid-liquid interface height decreased from about 80 mm at initial stage to 40 mm after mid-stage due to achieve the growth speed balance. Also the optimal input power of the modified system was similar with the original one due to the compensation effects of the crucible geometry and additional insulation. The crystal temperature grown by the modified CZ grower was increased about 10 K than the original one. Therefore the sapphire crystal of 300 mm height was grown successfully.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.2
• /
• pp.139-144
• /
• 2015

Recently, to realize sound-absorbing structures, we have to insert sound-absorbing materials into wall. These shapes are taken limitations because sound-absorbing materials should be fixed. Therefore, the sound absorption is changed by environment that used the sound-absorbing materials. On the other hand, we will take same effect without sound-absorbing material, if we change the shape of wall to sound absorbing structure. If we use this sound absorbing structure, we can get benefits by removing limitation of materials. Therefore we suggest perforated plate for effective sound-absorbing structure. We confirmed the function of sound-absorption of this structure using equivalent property. Then, we found the similarity between perforated plate and resonator. Also, we verify these theories through computer simulation by FEM(Finite Element Method). Finally, we validated that perforated plate has function of sound absorption without sound-absorbing material. This perforated plate is used for sound-absorbing material of buildings and transportations such as vehicle, train etc. Also, these results could be further used basic tool for design of sound-absorption structure.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.33 no.3
• /
• pp.171-178
• /
• 2020

In this study, we develop MPI-OpenMP hybrid parallelization for multibody peridynamic simulations. Peridynamics is suitable for analyzing complicated dynamic fractures and various discontinuities. However, compared with a conventional finite element method, nonlocal interactions in peridynamics cost more time and memory. In multibody peridynamic analysis, the costs increase due to the additional interactions that occur when computing the nonlocal contact and ghost interlayer models between adjacent bodies. The costs become excessive when further refinement and smaller time steps are required in cases of high-velocity impact fracturing or similar instances. Thus, high computational efficiency and performance can be achieved by parallelization and optimization of multibody peridynamic simulations. The analytical code is developed using an Intel Fortran MPI compiler and OpenMP in NURION of the KISTI HPC center and parallelized through MPI-OpenMP hybrid parallelization. Further parallelization is conducted by hybridizing with OpenMP threads in each MPI process. We also try to minimize communication operations by model-based decomposition of MPI processes. The numerical results for the impact fracturing of multiple bodies show that the computing performance improves significantly with MPI-OpenMP hybrid parallelization.

• Journal of the Korean Geosynthetics Society
• /
• v.15 no.3
• /
• pp.27-37
• /
• 2016

Although many studies on the Granular Compaction Pile (GCP) have been done by many researchers, the GCP design has not been systematically done due to the absence of the rational design methodology. As the GCP design has been mostly done by engineers' own experiences, some failure cases have been reported to occur. For this reason, it is very difficult to confirm definite causes of the failure and establish the prevention plans for the failure. Therefore, this study aims to investigate the optimal mixing ratio of gravel and sand, the effects of the internal friction angle of the GCP on the stress concentration ratio and the vertical and horizontal settlements. In order to analyze the behavior of the soft ground reinforced with the GCP depending on the different design parameters such as the stress concentration ratio and the internal friction angle, a number of finite element (FE) analyses were performed. From the direct shear test, the optimal mixing ratio of gravel to sand was found to be 70:30. Based on the numerical analyses, as the internal friction angle increased, the stress concentration ratio increased and it converged to a constant value. In addition, the larger the internal friction angle, the smaller the settlements. Consequently, the use of the optimal mixing ratio of gravel and sand can lead to reducing both the lateral flow and the heaving phenomenon.

• Tunnel and Underground Space
• /
• v.20 no.4
• /
• pp.267-276
• /
• 2010

In modern rock engineering practice, fully grouted rock bolting is actively employed as a major supporting system, so that understanding the behavior of fully grouted rock bolts is essential for the precise design of rock bolting. Despite its importance, the supporting mechanism of rock bolts has not been fully understood yet. Since most of existing analytical models for rock bolts were developed by drastically simplifying their boundary conditions, they are not suitable for the bolts of in-situ condition. In this study, 3-D elastic FE analysis of fully grouted rock bolts has been conducted to provide insight into the supporting mechanism of the bolt. The distribution of shear and axial stresses along the bolt are investigated with the consideration of different boundary conditions including three different displacement boundary conditions at the bolt head, the presence of intersecting rock joints, and the variation of elastic modulus of adjacent rock. The numerical result reveals that installation of the faceplate at the bolt head plays an important role in mobilizing the supporting action and enhancing the supporting capabilities of the fully grouted rock bolts.