• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.1C
• /
• pp.9-17
• /
• 2006

The numerical analysis to predicting the electrokinetic remediation behavior on unsaturated soil is carried out by aiding HERO, Hanyang Unversity Electrokinetic Remediation program, developed from the finite difference method and in the VISUAL FORTRAN environment. The analysis for the pure kaolinite under saturated conditions is performed on the results of the previous study of Acar (1997). Also the predictions to the characteristics of electrokinetic remediation on unsaturated conditions are performed and the conclusions summarized as follows. First, pH of the electrolyte in the reservoirs is not different with the degree of saturation resulted from the changes in electrical efficiency. But the advance of acid front is increased dependent on the degree of saturation in contrary to the transportation of base front. Second, below the degree of saturation of 83%, which is equivalent to the optimum water content, the removal effect increased with the decreasing of degree of saturation. But it have no effect on the efficiency of removal over the degree of saturation of 83%.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.1C
• /
• pp.33-39
• /
• 2009

Vertical soil stress near a pile subjected to negative skin friction (NSF) may be reduced due to shear transfer at the pile-soil interface. A three-dimensional finite difference analysis has been performed to clarify the influence of vertical and horizontal stress reductions on the pile behavour. In addition, a simple equation has been proposed to estimate vertical stress reduction of the soil near the pile. The vertical and horizontal stresses are reduced by substantial amount compared to corresponding stress components at the Greenfield condition. The horizontal extent of vertical stress reduction of the soil near the pile is rather limited to about up to 4-8 D, where D is the pile diameter. The findings from the current research indicate that widely used $\beta$-method may result in over-estimation of dragload (compressive force on piles due to NSF) and hence stress reduction needs to be incorporated in the original equation.

• Korean Journal of Optics and Photonics
• /
• v.34 no.6
• /
• pp.235-240
• /
• 2023

The optical signal injected into an end-fire optical phased array propagates along the waveguides inside the device and is emitted from the edge of the antenna. In general, reflection and scattering occur at the boundary, thereby reducing the emission efficiency of the optical signal. In this article, we propose a silicon nitride (Si₃N₄) tapered waveguide antenna structure whose width is tapered toward the emitting edge, achieving high emission efficiency operating at the 1,550 nm wavelength. The Si₃N₄ tapered waveguide antenna was numerically designed using the 3D finite-difference time-domain method. The optical signal emission efficiency increased from 78% to 96.3%, while reflectance decreased from 22% to 3.7% compared with the untapered waveguide antenna counterpart. This result will not only boost the optical signal intensity but also mitigate optical noise resulting from back reflection along the waveguide in the end-fire optical phased array device.

• Computers and Concrete
• /
• v.32 no.4
• /
• pp.351-363
• /
• 2023

Reducing the self-weight of reinforced concrete structures problem is discussed in this paper by using two types of self-weight reduction, the first is by using lightweight coarse aggregate (crushed brick) and the second is by using styropor block. Experimental and Numerical studies are conducted on (LWAC) lightweight aggregate reinforced concrete slabs, having styropor blocks with various sizes of blocks and the ratio of shear span to the effective depth (a/d). The experimental part included testing eleven lightweight concrete one-way simply supported slabs, comprising three as reference slabs (solid slabs) and eight as styropor block slabs (SBS) with a total reduction in cross-sectional area of (43.3% and 49.7%) were considered. The holes were formed by placing styropor at the ineffective concrete zones in resisting the tensile stresses. The length, width, and thickness of specimen dimensions were 1.1 m, 0.6 m, and 0.12 m respectively, except one specimen had a depth of 85 mm (which has a cross-sectional area equal to styropor block slab with a weight reduction of 49.7%). Two shear spans to effective depth ratios (a/d) of (3.125) for load case (A) and (a/d) of (2) for load case (B), (two-line monotonic loads) are considered. The test results showed under loading cases A and B (using minimum shear reinforcement and the reduction in cross-sectional area of styropor block slab by 29.1%) caused an increase in strength capacity by 60.4% and 54.6 % compared to the lightweight reference slab. Also, the best percentage of reduction in cross-sectional area is found to be 49.7%. Numerically, the computer program named (ANSYS) was used to study the behavior of these reinforced concrete slabs by using the finite element method. The results show acceptable agreement with the experimental test results. The average difference between experimental and numerical results is found to be (11.06%) in ultimate strength and (5.33%) in ultimate deflection.

• Journal of the Institute of Electronics Engineers of Korea CI
• /
• v.48 no.1
• /
• pp.90-100
• /
• 2011

We present a parallel bi-conjugate gradient (Bi-CG) matrix solver for large scale Bio-FET simulations based on recent graphics processing units (GPUs) which can realize a large-scale parallel processing with very low cost. The proposed method is focused on solving the Poisson equation in a parallel way, which requires massive computational resources in not only semiconductor simulation, but also other various fields including computational fluid dynamics and heat transfer simulations. As a result, our solver is around 30 times faster than those with traditional methods based on single core CPU systems in solving the Possion equation in a 3D FDM (Finite Difference Method) scheme. The proposed method is implemented and tested based on NVIDIA's CUDA (Compute Unified Device Architecture) environment which enables general purpose parallel processing in GPUs. Unlike other similar GPU-based approaches which apply usually 32-bit single-precision floating point arithmetics, we use 64-bit double-precision operations for better convergence. Applications on the CUDA platform are rather easy to implement but very hard to get optimized performances. In this regard, we also discuss the optimization strategy of the proposed method.

• Geophysics and Geophysical Exploration
• /
• v.12 no.1
• /
• pp.105-113
• /
• 2009

The interpretation of observed waveform characteristics identified in refraction and wide-angle reflection data increases confidence in the crustal structure model obtained. When calculating traveltimes and raypaths, wavefront methods on a regular grid based on graph theory are robust even with complicated structures, but basically compute only first arrivals. In this paper, we develop new algorithms to compute traveltimes and raypaths not only for first arrivals, but also for fast and later reflection arrivals, later refraction arrivals, and converted waves between P and S, using the modified wavefront method based on slowness network nodes mapped on a multi-layer model. Using the new algorithm, we can interpret reflected arrivals, Pg-later arrivals, strong arrivals appearing behind Pn, triplicated Moho reflected arrivals (PmP) to obtain the shape of the Moho, and phases involving conversion between P and S. Using two models of an ocean-continent transition zone and an oceanic ridge or seamount, we show the usefulness of this algorithm, which is confirmed by synthetic seismograms using the 2D Finite Difference Method (2D-FDM). Characteristics of arrivals and raypaths of the two models differ from each other in that using only first-arrival traveltime data for crustal structure analysis involves risk of erroneous interpretation in the ocean-continent transition zone, or the region around a ridge or seamount.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.137-145
• /
• 2014

We present a design sensitivity analysis(DSA) method for multiscale problems based on bridging scale decomposition. In this paper, we utilize a bridging scale method for the coupled system analysis. Since the analysis of full MD systems requires huge amount of computational costs, a coupled system of MD-level and continuum-level simulation is usually preferred. The information exchange between the MD and continuum levels is taken place at the MD-continuum boundary. In the bridging scale method, a generalized Langevin equation(GLE) is introduced for the reduced MD system and the GLE force using a time history kernel is applied at the boundary atoms in the MD system. Therefore, we can separately analyze the MD and continuum level simulations, which can accelerate the computing process. Once the simulation of coupled problems is successful, the need for the DSA is naturally arising for the optimization of macro-scale design, where the macro scale performance of the system is maximized considering the micro scale effects. The finite difference sensitivity is impractical for the gradient based optimization of large scale problems due to the restriction of computing costs but the analytical sensitivity for the coupled system is always accurate. In this study, we derive the analytical design sensitivity to verify the accuracy and applicability to the design optimization of the coupled system.

• Journal of the Korean Wood Science and Technology
• /
• v.44 no.3
• /
• pp.370-380
• /
• 2016

The black pine logs damaged by pine wilt disease in Jeju-do were heat-treated to extend the utilization of domestic trees damaged by pine wilt disease. The heat-treatment of wood requires wood to be heated to $56^{\circ}C$ for 30 min at the core. The average moisture content and top-diameter of the black pine logs were ranged from 46% to 141% and from 180 mm to 500 mm, respectively. And the basic specific gravity and oven-dry specific gravity of the black pine logs were 0.47 and 0.52, respectively. The time required for heat-treatment at $105^{\circ}C$ temperature was ranged from 7.7 h to 44.2 h, depending on moisture content and top-diameter. The temperature distribution was used to predict the time required for heat-treatment of black pine log with various moisture contents and top-diameters using finite difference method. The thermal properties of wood including the thermal conductivity and specific heat in accordance with moisture content were calculated. Heat transfer coefficient for mixed convection in form of adding natural convection and forced convection was used for heat transfer analysis. The error between the measured and predicted values ranged from 3% to 45%. The predicted times required for heat-treatment of black pine log with 50% moisture content and 200 mm, 300 mm, and 400 mm top-diameter were 10.9 h, 18.3 h, and 27.0 h, respectively. If the initial moisture content of black pine log is 75%, heat treatment times of 13.6 h, 22.5 h, and 32.8 h were predicted in accordance with top-diameter. And if the initial moisture content of black pine log is 100%, heat treatment times of 16.2 h, 26.5 h, and 38.2 h were predicted in accordance with top-diameter. When the physical properties of logs damaged by pine wilt disease are presented, these results can be applicable to the heat-treatment of red pine and Korean pine logs as well.

• Korean Journal of Agricultural Science
• /
• v.12 no.1
• /
• pp.118-127
• /
• 1985

The grain temperature is a crucial factor determining the deterioration rate of stored grain. Therefore, it is used to be predicted in order to evaluate the various stored methods rapidly and inexpensively. In this study, a mathematical model was developed to simulate the temperatures of grain stored in a cylinderical bin. It was formulated for the two dimensional heat transfer by the finite difference method. Then, it was verified statistically using the actual test deta and the predicted. The changes of grain temperature were analyed using the simulated data of one year for a safe stoarge and the following results were obtained: 1. Simulation model developed by the finite difference method was validated with the actual and the predicted grain temperatures and it's result showed that it could predict the grain temperature of storage bin reasonably well. 2. Grain temperature near the wall of storage bin were changed with $6-7^{\circ}C$ higher then average atmospheric temperature from June to September. Therefore, the parts of stored grain near the wall is supposed to be deteriorated fast. 3. When the dimension of bin diameter is about the same as the bed height, the changes of grain temperature of radial direction was higher than the verticals. 4. The predicted temperature showed that the grain temperature of which were from the end of April to mid October were higher than the safe storage limit at Yusung, Korea.

• Journal of the Korean Geotechnical Society
• /
• v.39 no.5
• /
• pp.51-63
• /
• 2023

Numerical analysis was conducted to determine the effect of soil behavior by thawing and freezing of seasonal frozen soil on pile foundations. The analysis was performed using the finite element method (FEM) to simulate soil-pile interaction based on the atmosphere temperature change. Thermomechanical coupled modeling using FEM was applied with the temperature-dependent nonlinear properties of the frozen soil. The analysis model cases were applied to the MCR and HDP models to simulate the elastoplastic behavior of soil. The numerical analysis results were analyzed and compared with various conditions having different length and width sizes of the pile. The results of the numerical analysis showed t hat t he HDP model was relat ively passive, and t he aspect and magnit ude of t he bearing capacit y and displacement of the pile head were similar depending on the length and width of the pile conditions. The vertical displacement of the pile head by thawing and freezing of the ground showed a large variation in displacement for shorter length conditions. In the MCR model, the vertical displacement appeared in the maximum thaw settlement and frost heaving of 0.0387 and 0.0277 m, respectively. In the HDP model, the vertical displacement appeared in the maximum thaw settlement and frost heaving of 0.0367 and 0.0264 m, respectively. The results of the pile bearing capacity for the two elastoplastic models showed a larger difference in the width condition than the length condition of the pile, with a maximum of about 14.7% for the width L condition, a maximum of about 5.4% for M condition, and a maximum of about 5.3% for S condition. The significance of the effect on the displacement of the pile head and the bearing capacity depended on the pile-soil contact area, and the difference depended on the presence or absence of an active layer in the soil and its thickness.