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

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.14 no.6
• /
• pp.575-593
• /
• 2012

This paper briefly introduces the scope and objectives of SKB Task 8, which is an international cooperative research project. In addition, the hydraulic behaviors of bentonite buffer focusing on the interactions between bentonite and a rock mass with a joint were investigated using TOUGH2 code as part of a sub-mission of Task 8a. The effects of a rock joint and high capillary pressure of bentonite on the re-saturation properties and pressure distribution in a buffer were identified and successfully incorporated in the TOUGH2 code. Based on the numerical results, it was found that the speed of re-saturation in bentonite surrounded by a rock mass with a joint is 2.5 to 12 times faster than that in a condition without a rock joint, while the degree of saturation in the lower part of the buffer material is generally higher than in the upper part in both the cases of with and without a joint. It can be anticipated that the results obtained from this study can be applied to an estimation of the full saturation time and a determination of optimum thickness with regard to the design of the bentonite buffer in a high level waste disposal system.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.6
• /
• pp.614-627
• /
• 2018

Urban railway lines have been constructed adjacent to residential buildings and urban areas. The expansion of transportation networks and reconstruction of residential buildings in highly populated urban areas require deep excavations in areas adjacent to urban railways. Mobilized soil stresses and changes in the groundwater level induced by deep excavations results in track irregularities in urban railways. In this study, a three-dimensional finite difference model using the commercial program FLAC3D was adopted to estimate the horizontal displacements of earth retaining structures, settlements of backfill, the stability of track irregularity and underground box structure based on the criteria of each railway organization and its relationships. In deep excavations, a change in groundwater level induces relatively very small differences for track gauge irregularities, whereas relatively large differences for longitudinal irregularities of 72.5%, twist irregularities of 83.3%, cross level irregularities of 61.9%, and alignments of 43.3% were found to be the maximum differences when the horizontal displacement of earth retaining wall and settlement of backfill were 65.1% and 21.4%, respectively, because the groundwater level (GWL) on the ground surface-mobilized tensile strength of the underground box structure exceeds the allowable value. Therefore, three-dimensional numerical analysis was performed in this study. Overall, real-time monitoring should be carried out to prevent railway accidents in advance when a deep excavation adjacent to urban railway structures is constructed.

• Journal of Korea Water Resources Association
• /
• v.41 no.8
• /
• pp.773-784
• /
• 2008

The present study was carried out to examine flow properties in linear drainage channels such as road surface drainage facilities. The finite difference formulation for the varied flow analysis was solved for flow profiles in the channels. Starting the first step at the control section, the Newton-Raphson method was applied for producing numerical solutions of the equation. We considered two types of linear drainage channels, a channel with one outlet at downstream end and a channel with two outlets at both ends. Moreover, the flow analysis for various channel slopes was performed. However, we considered channels with the two outlets of slopes satisfying the condition that the both ends are the control section. The maximum of those slopes was decided from the relation between the channel slope and the location of control section. The flow of a channel with one outlet was calculated upward and downward from the control section existing in channel or upward from the control section at downstream end. The flow of a channel with two outlets at both ends were calculated for upstream and downstream channel segments divided by the water dividend, respectively and the flow analysis was completed when the water depth at the water dividend calculated from upstream end was equal to that calculated from downstream end. If the slope was larger than the critical slope, the channel with two outlets was likely to behave like the channel with one outlet. The maximum water depth was investigated and compared with that calculated additionally from the uniform flow analysis. The uniform flow analysis was likely to lead a excessive design of a drainage channel with mild slope.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.32 no.10C
• /
• pp.1032-1041
• /
• 2007

The global navigation satellite system (GNSS), which is the core technique for the location based service, adopts the direct sequence/spread spectrum (DS/SS) as its modulation method. The success of a DS/SS system depends on the synchronization between the received and locally generated pseudo noise (PN) signals. As a step in the synchronization process, the tacking scheme performs fine adjustment to bring the phase difference between the two PN signals to zero. The most widely used tracking scheme is the delay locked loop with early minus late discriminator (EL-DLL). In the ideal case, the EL-DLL is the best estimator among various DLL. However, in the band-limited multipath environment, the EL-DLL has tracking bias. In this paper, the timing offset range of correlation function is divided into advanced offset range (AOR) and delayed offset range (DOR) centering around the correct synchronization time point. The tracking bias results from the following two reasons: symmetry distortion between correlation values in AOR and DOR, and mismatch between the time point corresponding to the maximum correlation value and the synchronization time point. The former and latter are named as the type I and type II tracking bias, respectively. In this paper, when the receiver has finite bandwidth in the presence of multipath signals, it is shown that the type II tracking bias becomes a more dominant error factor than the type I tracking bias, and the correlation values in AOR are not almost changed. Exploiting these characteristics, we propose a novel tracking bias mitigation scheme and demonstrate that the tracking accuracy of the proposed scheme is higher than that of the conventional scheme, both in the presence and absence of noise.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.13 no.2
• /
• pp.122-138
• /
• 2001

It is well known that whcn waves propagating on a shallow water suddenly encounter a much deeper water they do not propagate further but are reflected. If we apply this phenomenon to a harbor by making the harbor depth much greater than outside, we could improve the harbor tranquillity by making the waves impinging into the harbor be reflected at the harbor entrance. In the present paper, first we apply the numerical models based' on the mild-slope equation and extended mild-slope equation to calculate the long wave resonances in a rectangular harbor with a very large depth discontinuity at its entrance to find that the difference between the models is almost negligible. By applying the numerical model to a realistic model harbor whose inside is entirely dredged, it is found that the effect of dredging is insignificant when the inside depth is twice the outside one but tripled inside depth significantly reduces the long waves of period of one to five minutes whieh may exert a bad influence on ship motion. Moreover, even when only a portion of the harbor basin is dredged, the cffect of dredging in the dredged area is found to be comparable to that of entire dredging, showing that the dredging of harbor basin can be a countenncasure for harbor resonance.

• Journal of Biomedical Engineering Research
• /
• v.23 no.6
• /
• pp.467-475
• /
• 2002

Hyperthermia using transrectal thermal probes has been used for a noninvasive treatment of prostate diseases. However it is known that heating the rectal wall at excessively high temperature can lead to destruction of the rectal mucous membrane. and it is difficult to maintain an optimum temperature over the entire prostate. Thus, a more accurate understanding of the heat transfer mechanism between prostate and hyperthermia system is needed Numerical analysis was performed to investigate how the cold/warm stimulations on the prostate surface affect the temperature distribution in the prostate model. The general purpose software "FLUENT" was used for obtaining a finite volume solution to the unsteady conduction equation and to calculate the time-varying temperature in the prostate. Effects of the warm/cold stimulations and the stimulation frequency on the temperature distribution were simulated. and we visualized how hyperthermia affected the inside of the prostate. It was found that the effect of hyperthermia by using a typical heating method is limited due to the low thermal conductivity of the prostate. Consecutive repetitions of warm and cold stimulations were considered to provide the thermal irritations inside a prostate. The effects of temperature difference and duration of warm/cold stimulations were investigated, and basic data for the optimum period and effective patterns of stimulations were obtained. A simplified bioheat equation was also solved to describe effects of the blood flow on the blood-tissue heat transfer. The effect of blood flow was not dominant compared to that of warm/cold stimulations. These results might be used as data for design of prostate treating probe, prostatic therapy and thermal stimulation effects on the prostate.

• Journal of Korean Society of Steel Construction
• /
• v.9 no.4 s.33
• /
• pp.659-672
• /
• 1997

Thermal response induced from nonlinear temperature distribution in composite box gilder bridges depends on several variables(environmental conditions, physical and material properties, location and orientation of bridge, and cross-section geometry). In this paper, parametric study are conducted in order to find the effects of variations of seasons, location and orientation of bridge, sectional geometry and some material properties on the axial deformation, curvature and stresses in composite box girder bridge. A two-dimensional transient finite element model to conduct this parametric studies is briefly presented. Firstly, the effects of the parameters on the diurnal variation of curvature are considered, and for the time of maximum curvature, on the distribution of temperature and stresses of composite box girder sectional are considered. Finally, some considerations about the influence of the parameters on the daily maximum values of axial deformation, curvature and stresses are carried out. The influence of thermal effect on structures is important as much as the influence of live or dead load in some cases. In the design of steel composite bridges, the thermal stresses calculated on the supposition that the temperature difference between the concrete slab and steel girder is $10^{\circ}C$ and the temperature distributions are uniform in concrete slab and steel girder can be underestimated.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2013.04a
• /
• pp.835-841
• /
• 2013

Modern solid-state gyroscopes (HRG) with hemispherical resonators from high-purity quartz glass and special surface superfinishing and ultrathin gold coating become the best instruments for precise-grade inertial reference units (IRU) targeting long-term space missions. Designing of these sensors could be a notable contribution into development of Korea as a space nation. In participial, 40mm diameter thin-shell resonator from high-purity fused quartz, fabricated as a single-piece with its supporting stem has been designed, machined, etched, tuned, tested, and delivered by STM Co. (ATS of Ukraine) several years ago; an extremely-high Q-factor (upto 10~20 millions) has been shown. Understanding of the best way how to match such a unique sensor with inner glass assembly of the gyro means how to use the high potential in a maximal extent; and this has become the urgent task. Inner quartz glass assembly has a very thin indium (In) layer soldered the resonator and its silica base (case), but effects of internal resonances between operational modal pair of the shell-cup and its side (parasitic) modes can notable degrade the potential of the sensor as a whole, instead of so low level of resonator's intrinsic losses. Unfortunately, there are special combinations of dimensions of the parts (so-called, "resonant sizes"), when intensive losses of energy occurs. The authors proposed to use the length of stem's fixture as an additional design parameter to avoid such cases. So-called, a cyclic scheme of finite element method (FEM) and ANSYS software were employed to estimate different combinations of gyro assembly parameters. This variant has no mismatches of numerical origin due to FEM's discrete mesh. The optimum length and dangerous "resonant lengths" have been found. The special attention has been paid to analyses of 3D effects in a cup-stem transient zone, including determination of a difference between the positions of geometrical Pole of the resonant hemisphere and of its "dynamical Pole", i.e., its real zone of oscillation node. Boundary effects between the shell (cup) and 3D short "beams" (inner and outer stems) have been ranged. The results of the numerical experiments have been compared with the classic model of a quasi-hemispherical shell band with inextensional midsurface, and the solution using Rayleigh's functions of the $1^{st}$ and $2^{nd}$ kinds. To guarantee the truth of the recommended sizes to a designer of the real device, the analytical and FEM results have been compared with experimental data for a party of real resonators. The consistency of the results obtained by different means has been shown with errors less than 5%. The results notably differ from the data published earlier by different researchers.

• Journal of Environmental Impact Assessment
• /
• v.22 no.3
• /
• pp.219-230
• /
• 2013

Species distribution modeling is one of the most effective habitat analysis methods for wildlife conservation. This study was for evaluating the suitability of species distribution to distance between forest patches in Seoul city using tits. We analyzed the distribution of the four species of tits: varied tit (Parus varius), marsh tit (P. palustris), great tit (P. major) and coal tit (P. ater), using the landscape indexes and connectivity indexes, and compared the resulting suitability indexes from 100m to 1,000m. As factors affecting to the distribution of tits, we calculated landscape indices by separating them into intra-patch indices (i.e. logged patch area (PA), area-weighted mean patch shape index (PSI), tree rate (TR)) and inter-patch indices (i.e. patch degree (PD), patch betweenness (PB), difference probability of connectivity (DPC)), to analyze the internal properties of the patches and their connectivity by tits occurrence data using logistic regression modeling. The models were evaluated by AICc (Akaike Information Criteria with a correction for finite sample sizes) and AUC (Area Under Curve of ROC). The results of AICc and AUC showed DPC, PA, PSI, and TR were important factors of the habitat models for great tit and marsh tit at the level of distance 500~800m. In contrast, habitat models for coal tit and varied tit, which are known as forest interior species, reflected PA, PSI, and TR as intra-patch indices rather than connectivity. These mean that coal tit and varied tit are more likely to find a large circular forest patch than a small and long-shaped forest patch, which are higher rate of forest. Therefore, different strategies are required in order to enhance the habitats of the forest birds, tits, in a region that has fragmented forest patches such as Seoul city. It is important to manage forest interior areas for coal tit and varied tit, which are known as forest interior species and to manage not only forest interior areas but also connectivity of the forest patches in the threshold distance for great tit and marsh tit as adapted species to the urban ecosystem for sustainable ecosystem management.