• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.153-160
• /
• 2008

Seismic modeling is used to simulate wave propagation in the earth. Although the earth's subsurface is usually semi-infinite, we cannot handle the semi-infinite model in seismic modeling because of limited computational resources. For this reason, we usually assume a finite-sized model in seismic modeling. In that case, we need to eliminate the edge reflections arising from the artificial boundaries introducing a proper boundary condition. In this study, we changed three kinds of boundary conditions (sponge boundary condition, Clayton and Engquist's absorbing boundary condition, and Higdon's transparent boundary condition) so that they can be applied in elastic wave modeling for anisotropic media. We then apply them to several models whose Poisson's ratios are different. Clayton and Engquist's absorbing boundary condition is unstable in both isotropic and anisotropic media, when Poisson's ratio is large. This indicates that the absorbing boundary condition can be applied in anisotropic media restrictively. Although the sponge boundary condition yields good results for both isotropic and anisotropic media, it requires too much computational memory and time. On the other hand, Higdon's transparent boundary condition is not only inexpensive, but also reduce reflections over a wide range of incident angles. We think that Higdon's transparent boundary condition can be a method of choice for anisotropic media, where Poisson's ratio is large.

• Proceedings of the Acoustical Society of Korea Conference
• /
• 1998.06c
• /
• pp.445-448
• /
• 1998

고유음향 임피던스가 다른 물질이 층상구조를 이루고 있는 매질내에서 임의의 한 층에 형성된 이차원 초음파 음장이 다른 매질로 전파될 때에의 특성을 각스펙트럼법에 의해 해석하였다. 음원으로서는 baffle된 원형 트랜스듀서를, 층상구조를 매질로서는 물과 생체조직(지방)의 경계를 각각 고려하였다. 해석에 있어서는 먼저 각스펙트럼 분해요소, 즉, 음원의 크기 및 거리에 의존하는 기준음장의 최적 크기와 샘플링 간격을 균질매질을 이용하여 결정한 다음, 그 요소를 층상매질에 적용하였다. 그 결과, 경계면에서의 입사각(공간주파수)의 의존하는 투과계수에 입사음장의 곱으로 나타내어지는 투과음장은 각스펙트럼법에 의해 짧은 시간내에 정도 높게 산출될 수 있음을 보인다.

• The Optical Journal
• /
• s.97
• /
• pp.58-61
• /
• 2005

복굴절이란 하나의 매질로부터 다른 매질로 진입하는 파동이 그 경계면에서 나가는 방향을 바꾸는 현상으로 등방성 매질에서 이방성 매질로 나아갈 때는 경계면에서 굴절파가 두 개로 나뉘어 굴절하게 되는 것을 말한다. 휴대폰용 카메라모듈 렌즈의 사출성형공정에서 발생하는 렌즈 내의 잔류응력에 의한 복굴절률은 사출공정에서의 렌즈불량 발생의 원인이 되며, 잔류응력은 또한 사출성형 이후에 완화되어 렌즈의 외관 치수의 변형을 가져오게 된다. 그러나 플라스틱류를 사용하는 정밀가공 생산 공정 등 실생활에서는 제품의 제조에 있어서 복굴절은 광학적 특성을 결정하는 가장 중요한 요소 중의 하나로서 성형시 복굴절을 최소화할 필요가 있다. 이를 어느 정도 달성했는가의 수준에 따라 플라스틱 정밀 성형 제품의 등급이 결정되는 중요한 측정요소이다.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.10 no.2
• /
• pp.180-186
• /
• 1999

Conformal methods such as staircase approximation and effective medium theory have been used to model arbitrary dielectric surfaces in Cartesian FDTD(finite-difference time-domain) grids. However, the staircase approximation doesn't guarantee accurate results and the effective medium theory cannot be applied for frequency dispersive medium. In this paper, a simple conformal method is introduced for analyzing arbitrary dielectric surfaces. These arbitrarily shaped dielectric surfaces make nonuniformly filled FDTD cells. E(H)-field in a nonuniformly filled FDTD cell is calculated by considering the cell as the combination of two kinds of uniformly filled cells whose material boundaries are normal and parallel to the E(H)-field in the cell. TM scattering from a rotated 2-D dielectric and a ferrite square cylinder is analyzed to show the accuracy of the proposed method.

• Korean Journal of Optics and Photonics
• /
• v.15 no.2
• /
• pp.89-99
• /
• 2004

At the interface of two materials with frequency-dependent material-parameters of permittivity and permeability, there may exist two kinds of surface polaritons: surface electric-polaritons(SEPs) and surface magnetic-polaritons(SMPs). Possible combinations of the material-parameters to support propagation of the two surface polaritons are suggested at the interface between metals and metamaterials such as a left-handed material. Dispersion relations are also derived in order to characterize frequency dependence of propagation of the SEP and SMP. It is found that only one propagation mode of SEP or SMP is allowed at a given set of four material parameters, and that counter-propagation of the phase and group velocities of the propagation mode can be observed even in the case when there are no double negative(or, negative-index) materials. Physical origin of the counter-propagation of the group velocity is proposed by evaluating the ratio of two electromagnetic-energy densities of a surface polariton propagating along within the two interface media, and it is confirmed by the dispersion relations.

• Economic and Environmental Geology
• /
• v.44 no.2
• /
• pp.161-170
• /
• 2011

The multiphase flow simulator, CHEMPS, was developed based on the fractional flow approach reported in the petroleum engineering literature considering fully three phase flow in physically and chemically heterogeneous media. It is a extension of MPS developed by Suk and Yeh (2008) to include the effect of wettability on the migration of NAPL. The fractional flow approach employs water, total liquid saturation and total pressure as the primary variables. Most existing models are limited to two-phase flow and specific boundary conditions when considering physically heterogeneous media. In addition, these models focused mainly on the water-wet media. However, in a real system, variations in wettability between water-wet and oil-wet media often occur. Furthermore, the wetting of porous media by oil can be heterogeneous, or fractional, rather than uniform due to the heterogeneous nature of the subsurface media and the factors that affect the wettability. Therefore, in this study, the chemically heterogeneous media considering fractional wettability as well as physically heterogeneous media were simulated using CHEMPS. In addition, the general boundary conditions were considered to be a combination of two types of boundaries of individual phases, flux-type and Dirichlet type boundaries.

• Geophysics and Geophysical Exploration
• /
• v.6 no.2
• /
• pp.77-86
• /
• 2003

We designed a new time-domain, finite-difference, elastic wave modeling technique, based on a displacement formulation. which yields nearly correct solutions to Lamb's problem. Unlike the conventional, displacement-based, finite-difference method using a node-based grid set (where both displacements and material properties such as density and Lame constants are assigned to nodal points), in our new finite-difference method, we use a cell-based grid set (where displacements are still defined at nodal points but material properties within cells). In the case of using the cell-based grid set, stress-free conditions at the free surface are naturally described by the changes in the material properties without any additional free-surface boundary condition. Through numerical tests, we confirmed that the new second-order finite differences formulated in the cell-based grid let generate numerical solutions compatible with analytic solutions unlike the old second-order finite-differences formulated in the node-based grid set.

• Nuclear Engineering and Technology
• /
• v.21 no.4
• /
• pp.267-276
• /
• 1989

A numerical model for nuclide migration through fractured rock media has been developed. Nuclide transport with groundwater in rock fissures and the diffusion of nuclides into rock matrix are considered one-dimensionally . In the safety assessment of the repository for radioactive waste, this one-dimensional model by the finite-difference scheme, which enables us not only to use more realistic boundary conditions but also to model the nonhomogeneous rock medium as the multilayered media, can be used effectively with the analytical mode. The solution by the numerical model will be verified analytically, and then extended to the double-layered rock medium transport model.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.318-325
• /
• 2008

Assigning an exact boundary condition is of great importance in three-dimensional (3D) magnetotelluric (MT) modeling, in which no source is considered in a computing domain. This paper presents a 3D MT modeling algorithm utilizing a Dirichlet condition for a 2D host. To compute boundary values for a model with a 2D host, we need to conduct additional 2D MT modeling. The 2D modeling consists of transverse magnetic and electric modes, which are determined from the relationship between the polarization of plane wave and the strike direction of the 2D structure. Since the 3D MT modeling algorithm solves Maxwell's equations for electric fields using the finite difference method with a staggered grid that defines electric fields along cell edges, electric fields are calculated at the same place in the 2D modeling. The algorithm developed in this study can produce reliable MT responses for a 3D model with a 2D host.

• Proceedings of the Acoustical Society of Korea Conference
• /
• 1996.10a
• /
• pp.8-13
• /
• 1996

균질 매질내에 형성되는 근거리 음장의 해석에는 종래 회절 이론에서의 Rayleigh 적분식에 근거한 Lommel 근사식이 주로 이용되어져 왔다. 그러나, 불균질한 매질에서는 그 방법의 적용이 어렵고, 유한요소법등 새로운 방법이 적용되어야 한다. 본 연구에서는 불균질 매질에서의 음장을 해석하는 한 방법으로 Angular Spectrum법을 제안하고, 경계를 가상한 등수온의 물에 대하여 Lommel 근사식에 의해 직접 산출한 음장과 Angular Spectrum법을 적용하여 산출한 음장을 서로 비교함으로써 그 유효성을 검토하였다.