• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.1
• /
• pp.23-32
• /
• 1999

A new numerical method using solution-adaptive grids (SAG) is developed to solve the Richards' equation (RE) for unsaturated flow in porous media. Using a grid generation technique, the SAG method automatically redistributes a fixed number of grid points during the flow process, so that more grid points are clustered in regions of large solution gradients. The method uses the coordinate transformation technique to employ a new transformed RE, which is solved with the standard finite difference method. The movement of grid points is incorporated into the transformed RE, and therefore all computation is performed on fixed grid points of the transformed domain without using any interpolation techniques. Thus, numerical difficulties arising from the movement of the wetting front during the infiltration process have been substantially overcome by the new method. Numerical experiments for an one-dimensional infiltration problem are presented to compare the SAG method to the modified Picard method using a fixed grid. Results show that accuracy of a SAG solution using 41 nodes is comparable with the solution of the fixed grid method using 201 nodes, while it requires only 50% of the CPU time. The global mass balance and the convergence of SAG solutions are strongly affected by the time step size (Δt) and the weighting parameter (${\gamma}$) used for generating solution-adaptive grids. Thus, the method requires automated readjustment of Δt and ${\gamma}$ to yield mass-conservative and convergent solutions, although it may increase computational costs. The method can be effective especially for simulating unsaturated flow and other transport problems involving the propagation of a sharp-front.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.10 no.4
• /
• pp.550-563
• /
• 1999

A wavelet-Galerkin scheme based on the time-dependent Maxwell's equations is presented. Daubechies wavelet with two vanishing wavelet moments is expanded for basis function in spatial domain and Yee's leap-frog approach is applied. The shifted interpolation property of Daubechies wavelet family leads to the simplified formulations for inhomogeneous media without the additional matrices for the integral or material operator. The stability condition is formulated. The dispersion characteristics are analyzed and compared with those of finite difference time domain and multiresolution time domain methods. The analyses show the excellent trade-off between the regularity and the support width of the basis function. Although the basis function has only two vanishing wavelet moments, it is enough to provide negligible dispersive error in the numerical analysis and its compact support enables only several involved terms per nodes. The storage effectiveness, execution time reduction and accuracy of this scheme are demonstrated by calculating the resonant frequencies of the homogeneous and inhomogeneous cavities.

• Geophysics and Geophysical Exploration
• /
• v.3 no.4
• /
• pp.119-124
• /
• 2000

A new fast algorithm for travel time calculation using mono-chromatic one-way wave equation was developed based on the delta function and the logarithms of the single frequency wavefield in the frequency domain. We found an empirical relation between grid spacing and frequency by trial and error method such that we can minimize travel time error. In comparison with other methods, travel time contours obtained by solving eikonal equation and the wave front edge of the snapshot by the finite difference modeling solution agree with our algorithm. Compared to the other two methods, this algorithm computes travel time of directly transmitted wave. We demonstrated our algorithm on migration so that we obtained good section showing good agreement with original model. our results show that this new algorithm is a faster travel time calculation method of the directly transmitted wave for imaging the subsurface and the transmission tomography.

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 2008.06a
• /
• pp.80-80
• /
• 2008

바이오 센서 응용 연구에 많이 사용되는 금(Au) 나노 입자를 이용한 국소 표면 플라즈몬 공명(Localized Surface Plasmon Resonance, LSPR)에 의한 산란광을 검출하는데 주로 이용되는 암시야(dark field) 현미경 검출 방식에 관한 전산모사를 통하여 입사광의 입사 방식에 따른 산란광 세기를 정량적으로 분석하였다. 전산모사 기법으로는 국소 표면 플라즈몬 공명의 동역학적인 현상을 모사할 수 있는 유한차분시간영역(Finite Difference Time Domain, FDTD) 기법을 이용하였는데, 이러한 기법이 암시야 현미경 전산 모사에 유효함을 우선적으로 검증하였다. 암시야 현미경 검출 방식의 모사에서 입사 광원의 반사 입사 방식과 투과 입사 방식을 비교하였고, 각각의 방식에 서 입사광의 입사각에 따른 산랑광 세기를 계산하였다. 이러한 전산모사를 통하여 프리즘을 통한 내부 전반사(Total Internal Reflection, TIR) 방식에서 입사 광원의 임계각 근처에서 많이 발생하는 에바네슨트 장(evanescent field)을 결합하는 경우 산란광 세기가 증가함을 관찰하였고, 이러한 세기의 변화를 프레넬(Fresnel) 방정식에 의해 계산된 에바네슨트 장의 세기 분포와 비교 분석하였다.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.2
• /
• pp.1-12
• /
• 1996

Hydrodynamic characteristics of shipping water on deck are investigated by using a simplified two-dimensional model. Formulation of the shipping water on deck leads to a nonlinear hyperbolic system of equations based on the shallow-water wave theory. Time-domain solution of these equations are obtained numerically using a finite difference method which adopts predictor-corrector method for time-marching and 2nd upwind differencing method for convection term calculation. To confirm the validity of the present numerical method, we calculated some shallow-water wave problems accompanying a bore and compared the obtained results with the analytic solutions. We found good agreements between them. Though the calculation results of shipping water on deck, we show that the shipping water flows into the deck as a rarefying wave arid grows into a bore after colliding with a deck structure. Also we examined the effects of acceleration and slope of deck and found that they have significant influences on the flow of shipping water.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.22 no.10
• /
• pp.1005-1011
• /
• 2011

In this paper, we present a marching-on-in-degree(MOD) finite difference method(FDM) based on the Helmholtz wave equation for analyzing transient electromagnetic responses in a general dispersive media. The two issues related to the finite difference approximation of the time derivatives and the time consuming convolution operations are handled analytically using the properties of the Laguerre functions. The basic idea here is that we fit the transient nature of the fields, the flux densities, the permittivity with a finite sum of orthogonal Laguerre functions. Through this novel approach, not only the time variable can be decoupled analytically from the temporal variations but also the final computational form of the equations is transformed from finite difference time-domain(FDTD) to a finite difference formulation through a Galerkin testing. Representative numerical examples are presented for transient wave propagation in general Debye, Drude, and Lorentz dispersive medium.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.41 no.3
• /
• pp.45-54
• /
• 2004

Microwave characteristics of ridge type CPW traveling-wave(TW) electroabsorption modulator and photodetector are affected by the thickness of intrinsic layer, width of guiding layer, and the separation of signal and ground electrodes. These factors are determined effective index of microwave and characteristic impedance due to changing of capacitance(C) and inductance(L) of device. However, conventional equivalent circuit of TW-structure is approximated to microstrip and CPW transmission line by distribution of electric and magnetic fields, respectively. In this paper, we analyzed microwave characteristics of TW-structure and found more accurate value of C and L by using finite difference time domain (FDTD) method. These values are adopted circuit element of equivalent circuit. Microwave characteristics obtained by the FDTD and equivalent circuit model show good agreement.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.21 no.10
• /
• pp.1073-1084
• /
• 2010

In this paper, the bio-heat equation including thermoregulatory functions is solved for an anatomically based human head model comprised of 14 tissues to study the thermal implications of high-power exposure to electromagnetic(EM) fields due to half-wave dipole antenna both at 835 and 1,800 MHz. The dipole antenna is located at the side of the ear and the front of the eyes. The FDTD method has been used for the SAR computation. When solving the BHE, the thermoregulation function and sweating effetecs are included in order to predict more exact temperature increase. It is noted that an approximately proportional relationship between the tissues and the maximum temperature increase and the antenna power is not maintained when the thermoregulation and sweating effects are fully accounted for under high power exposure.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.2
• /
• pp.123-136
• /
• 1996

A three-dimensional Galerkin-FEM model which can handle the temporal and spatial variation of density is presented. The hydrostatic approximation is used and density effects are included by means of conservation equation of heat and the equation of state. The finite difference grids are used in the horizontal plane and a set of linear-shape functions is used for the vertical expansion. The similarity transform is introduced to solve resultant matrix equations. The proposed model was first applied to the density-driven circulation in an idealized basin in the presence of the heat exchange between the air and the sea. The advection terms in the momentum equation were ignored, while the convection terms were retained in the heat equation. Coefficients of the vertical eddy viscosity and diffusivity were fixed to be constant. Calculation in a non-rotating idealized basin shows that the difference in heat capacity with depth gives rise to the horizontal gradient of temperature. Consequently, there is a steady new in the upper layer in the direction of increasing depth with compensatory counter flow .in the lower layer. With Coriolis force, geostrophic flow was predominant due to the balance between the pressure gradient and the Coriolis force. As a test in region of irregular topography, the model is applied to the Yellow Sea. Although the resultant flow was very complex, the character of the flow Showed to be geostrophic on the whole.

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