• Journal of Korean Society of Environmental Engineers
• /
• v.32 no.7
• /
• pp.690-698
• /
• 2010

Contamination of groundwater resources by organic chemicals has become an issue of increasing environmental concern. Surfactant-enhanced aquifer remediation (SEAR) is widely recognized as one of the most promising techniques to remediate organic contaminations in-situ. Solutions of surfactant or surfactant with polymer are used to dramatically expedite the process, which in turn, may reduce the treatment time of a site compared to use of water alone. In the design of surfactant-based technologies for remediation of organic contaminated aquifers, it is very important to have a considerable analysis using extensive numerical simulations prior to full-scale implementation. This study investigated the formation and flow of microemulsions during SEAR of organic-contaminated aquifer using the finite difference model UTCHEM, a three-dimensional, multicomponent, multiphase, compositional model. The remediation process variables considered in this study were the sequence of injection fluids, the injection and extraction rate, the concentrations of polymer in surfactant slug and chase water, and the duration of surfactant injection. For each variable, temporal changes in injection and production wells and spatial distributions of relative saturations in the organic phase were compared. Cleanup time and cumulative organic recovery were also quantified. The study would provide useful information to design strategies for the remediation of nonaqueous phase liquid-contaminated aquifers.

• Korean Journal of Optics and Photonics
• /
• v.14 no.5
• /
• pp.545-554
• /
• 2003

Design of the optical modulator composed of a three-waveguide coupler and CPW traveling-wave electrodes was carried out. Switching phenomena of three-waveguide couplers were analyzed by using the coupled mode theory, and the coupling-lengths of the devices were calculated by means of the FDM. CPW traveling-wave electrodes were analysed by the CMM and SOR simulation technique in order to find the conditions of phase-velocity and impedance matching. Traveling-wave modulators were fabricated on z-cut LiNbO$_3$ substrate. Ti was in-diffused in LiNbO$_3$ to make waveguides and Au electrodes were built on the waveguides by the electrolyte technique. The fabricated modulator chip was end-polished, pig-tailed and packaged in a brass mount with K-connector. The insertion loss and the switching voltage of the optical modulator were about 4㏈ and 19V, respectively. Network analyzer was used to obtain the S parameter and the corresponding RF response. From the measurement, parameters of the traveling-wave electrodes were extracted to be Z$_{c}$= 45 Ω, N$_{eff}$=2.20, and $\alpha$$_{0}$=0.055/cm√GHZ. The measured optical response R($\omega$) was compared with the theoretically estimated one, showing both responses agree well. The measurement results revealed that 3㏈ bandwidth turned out to be about 13 GHz.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.5
• /
• pp.337-343
• /
• 2014

This paper presents a continuum-based shape design sensitivity analysis(DSA) method for crack propagation problems using a reproducing kernel method(RKM), which facilitates the remeshing problem required for finite element analysis(FEA) and provides the higher order shape functions by increasing the continuity of the kernel functions. A linear elasticity is considered to obtain the required stress field around the crack tip for the evaluation of J-integral. The sensitivity of displacement field and stress intensity factor(SIF) with respect to shape design variables are derived using a material derivative approach. For efficient computation of design sensitivity, an adjoint variable method is employed tather than the direct differentiation method. Through numerical examples, The mesh-free and the DSA methods show excellent agreement with finite difference results. The DSA results are further extended to a shape optimization of crack propagation problems to control the propagation path.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.3B
• /
• pp.293-303
• /
• 2011

An efficient diagonalized approximate factorization algorithm (DAF) is developed for the solution of three-dimensional incompressible viscous flows. The pressure-based, artificial compressibility (AC) method is used for calculating steady incompressible Navier-Stokes equations. The AC form of the governing equations is discretized in space using a second-order-accurate finite volume method. The present DAF method is applied to derive a second-order accurate splitting of the discrete system of equations. The primary objective of this study is to investigate the computational efficiency of the present DAF method. The solutions of the DAF method are evaluated relative to those of well-known four-stage Runge-Kutta (RK4) method for fully developed and developing laminar flows in curved square ducts and a laminar flow in a cavity. While converged solutions obtained by DAF and RK4 methods on the same computational meshes are essentially identical because of employing the same discrete schemes in space, both algorithms shows significant discrepancy in the computing efficiency. The results reveal that the DAF method requires substantially at least two times less computational time than RK4 to solve all applied flow fields. The increase in computational efficiency of the DAF methods is achieved with no increase in computational resources and coding complexity.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.35 no.5
• /
• pp.291-298
• /
• 2015

An array type probe for monitoring metal tubes is proposed in this paper which utilizes peak value and peak time of a pulsed eddy current(PEC) signal. The probe consists of an array of encircling coils along a tube and the outside of coils is shielded by ferrite to prevent source magnetic fields from directly affecting sensor signals since it is the magnetic fields produced by eddy currents that reflect the condition of metal tubes. The positions of both exciter and sensor coils are consecutively moved automatically so that manual scanning is not necessary. At one position of send-receive coils, peak value and peak time are extracted from a sensor PEC signal and these data are accumulated for all positions to form an array type peak value signal and an array type peak time signal. Numerical simulation was performed using the backward difference method in time and the finite element method for spatial analysis. Simulation results showed that peak value increases and the peak appears earlier as the defect depth or length increases. The proposed array signals are shown to be excellent in reflecting the defect location as well as variations of defect depth and length within the array probe.

• Tunnel and Underground Space
• /
• v.30 no.5
• /
• pp.473-483
• /
• 2020

This paper examines the vertical stress change concentrated on mine pillar which occurs due to the stress disturbance from opening excavation at room and pillar mine by FLAC3D, a finite difference method (FDM) software. The mesh size combination is decided with a careful consideration of relative error and run-time, then its performance is verified. A series of numerical analyses is conducted and the vertical stress at central pillar was observed for the test cases of 1×1 to 11×11 mine pillars, 40 m to 320 m depth with 40 m difference. The results show that the vertical stress of pillar approaches to the similar value with the value estimated by tributary area theory(TAT) when the development area (N_P) is increased or the height of overburden (H_OB) is decreased, while it is overestimated in the opposite case. Furthermore, it also represents that the vertical stress factor (VSF) converges to a specific value when the depth is increased whille keeping the development area identical.

• Journal of Korea Soil Environment Society
• /
• v.1 no.1
• /
• pp.89-102
• /
• 1996

An integrated model is presented to describe underground flow and mass transport, using a multicomponent multiphase approach. The comprehensive governing equation is derived considering mass and force balances of chemical species over four phases(water, oil, air, and soil) in a schematic elementary volume. Compact and systemati notations of relevant variables and equations are introduced to facilitate the inclusion of complex migration and transformation processes, and variable spatial dimensions. The resulting nonlinear system is solved by a multidimensional finite element code. The developed code with dynamic array allocation, is sufficiently flexible to work across a wide spectrum of computers, including an IBM ES 9000/900 vector facility, SP2 cluster machine, Unix workstations and PCs, for one-, two and three-dimensional problems. To reduce the computation time and storage requirements, the system equations are decoupled and solved using a banded global matrix solver, with the vector and parallel processing on the IBM 9000. To avoide the numerical oscillations of the nonlinear problems in the case of convective dominant transport, the techniques of upstream weighting, mass lumping, and elementary-wise parameter evaluation are applied. The instability and convergence criteria of the nonlinear problems are studied for the one-dimensional analogue of FEM and FDM. Modeling capacity is presented in the simulation of three dimensional composite multiphase TCE migration. Comprehesive simulation feature of the code is presented in a companion paper of this issue for the specific groundwater or flow and contamination problems.

• Geotechnical Engineering
• /
• v.10 no.1
• /
• pp.47-62
• /
• 1994

It has been well known that the rainfall-triggered rise of groundwater levels is one of the most important factors resulting the instability of the hillside slopes. Thus, the prediction of porewater pressure is an essential step in the evaluation of landslide hazard. This study involves the development and verification of numerical groundwater flow model for the prediction of groundwater flow fluctuations accounting for both of unsatu나toed flow and saturated flow on steep hillside slopes. The first part of this study is to develop a nomerical groundwater flow model. The numerical technique chosen for this study is the finitro element method in combination with the finite difference method. The finite element method is used to transform the space derivatives and the finite difference method is used to discretize the time domain. The second part of this study is to estimate the unknown model parameters used in the proposed numerical model. There were three parameters to be estimated from input -output record $K_e$, $\psi_e$, b. The Maximum -A-Posteriori(MAP) optimization method is utilized for this purpose, . The developed model is applied to a site in Korea where two debris avalanches of large scale and many landslides of small scale were occurred. The results of example analysis show that the numerical groundwater flow model has a capacity of predicting the fluctuation of groundwater levels due to rainfall reasonably well.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2006.04a
• /
• pp.166-170
• /
• 2006

본 연구에서는 폐광산 주변에 산재되어 있는 광미/광폐석을 처리하기 위하여 고형화 실증 실험을 수행하였다. 고형화 공정에서 흔히 사용하는 포틀랜드 시멘트와 MSG-E, MSG-N을 고화제로 사용하였으며 현장 광미 및 광폐석을 대상으로 고화체를 양생하고 고화체의 압축강도 및 중금속 용출 정도를 측정하였다. 고화체의 물리/화학적 특성을 비교하기 위해 광미/고화제 비율, 배합수/고화제 비율 그리고 고화체 양생기간을 실험인자로 설정하였다. 실험 결과 광미/고화제의 비율 1:1 만을 고려하더라도 중금속 용출의 급격한 감소가 이루어지는 것을 확인할 수 있었으며 광미/고화제의 비율을 3:1 이하로 유지시키는 경우, 고화체의 압축강도가 현행 폐기물 관리법(20조 관련)에서 규정하고 있는 차단형 매립시설 내부막의 압축강도 기준인 $0.21kgf/mm^2$ 보다 높은 것으로 나타났다. 다양한 pH를 갖는 수용액에 대하여 시간에 따른 고화체의 중금속 용출률을 측정한 결과, 수용액의 pH가 1과 13인 강산/강염기 용액에서 일부 중금속의 용출 농도가 지하수 생활용수 기준치를 초과하였으나, pH와 3 - 11인 경우에는 중금속 용출률이 급격히 감소하여 모두 기준치 이하를 나타내었다. 또한, pH가 1과 13인 수용액의 경우에도 고화체와 반응하는 시간이 증가할수록 고화체의 buffering 효과에 의해 수용액의 pH가 감소하였다. 이러한 결과는 현장에서 접촉수의 pH가 강산이나 강염기라 하여도, 고화체의 buffering 효과에 의해 시간이 지남에 따라 수용액의 pH가 낮아져 고화체로부터의 중금속 용출은 매우 감소할 것임을 의미한다.ss of an active application defined using the model. The technique is developed in a platform- and language-independent way, and it is algorithmic and can be automated by computer program. We give an example dealing with network auction to illustrate the use of the model and the verification technique.품으로 내부 온도분포를 측정하였으며, 유한차분법 프로그램으로 대류열전달계수를 결정하였다. 대류열전달계수는 792에서 2,107 W/m$^2$로 분석되었다. 대류열전달 계수는 액상식품과의 상대속도가 증가함에 따라서 증가하였고, 점도가 증가함에 따라서는 감소하였다.ce of precision/recall of 90.99%/92.52%, and 93.39%/93.41% respectively. 의한 변성에 부분적으로 보호 작용을 나타 낼 것으로 추정된다.경(製麴72時間頃)의 활성(活性)은 보리쌀국(麴), 밀가루국(麴), 찹쌀국(麴), 고구마국(麴)의 순이었다.험 결과 오전용 사료는 관행적인 산란계 배합사료에서 Ca공급제를 제외한 것을 급여하고, 오후용 사료는 Ca공급제를 3배 첨가한 T2처리로 15:00~16:00시에 교체급여를 하면 사료섭취량 감소와 사료비 절감면에서 바람직할 것으로 사료되며, 고에너지-고단백질-저Ca의 분말사료와 저에너지-저단백질-고Ca의 펠렛사료를 혼합급여하면 산란계의 사료

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

Using an isogeometric approach, a continuum-based shape design optimization method is developed for steady state power flow problems at high frequencies. In case the isogeometric method is employed to the shape design optimization, the NURBS basis functions used in CAD geometric modeling are directly utilized to embed the exact geometry into the computational framework so that the design parameterization for shape optimization is much easier than that in the finite element method and consequently provides the enhanced smoothness of design perturbations. Thus, exact geometric models can be used in both the response and the shape sensitivity analyses, where normal vector and curvature are continuous over the whole design space so that enhanced shape sensitivity can be expected. Through numerical examples, the developed isogeometric sensitivity is compared with finite difference one to provide excellent agreement. Also, it turns out that the proposed method works very well in the shape optimization problems.