• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.19 no.2
• /
• pp.121-129
• /
• 2008

In this paper, we present a fast analysis of multilayer microstrip fractal structure by using the fast multipole method (FMM). In the analysis, accurate spatial green's functions from the real-axis integration method(RAIM) are employed to solve the mixed potential integral equation(MPIE) with FMM algorithm. MoM's iteration and memory requirement is $O(N^2)$ in case of calculation using the green function. the problem is the unknown number N can be extremely large for calculation of large scale objects and high accuracy. To improve these problem is fast algorithm FMM. FMM use the addition theorem of green function. So, it reduce the complexity of a matrix-vector multiplication and reduce the cost of calculation to the order of $O(N^{1.5})$, The efficiency is proved from comparing calculation results of the moment method and Fast algorithm.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.97-106
• /
• 1989

This paper is on the modeling of two-dimensional groundwater flow, which is the first step of the development of Dynamic System Model for groundwater flow and pollutant transport in subsurface porous media. The particular features of the model are its versatility and flexibility to deal with as many real-world problems as possible. Points as well as distributed sources/sinks are included to represent recharges/pumping and rainfall infiltrations. All sources/sinks can be transient or steady state. Prescribed hydraulic head on the Dirichlet boundaries and fluxes on Neumann or Cauchy boundaries can be time-dependent or constant. Sources/sinks strength over each element and node, hydraulic head at each Dirichlet boundary node and flux at each boundary segment can vary independently of each other. Either completely confined or completely unconfined aquifers, or partially confined and partially unconfined aquifers can be dealt with effectively. Discretization of a compound region with very irregular curved boundaries is made easy by including both quadrilateral and triangular elements in the formulation. Large-field problems can be solved efficiently by including a pointwise iterative solution strategy as an optional alternative to the direct elimination solution methed for the matrix equation approximating the partial differential equation of groundwater flow. The model also includes transient flow through confining leaky aquifers lying above and/or below the aquifer of interest. The model is verified against three simple cases to which analytical solutions are available. The groundwater flow model shall be combined with the model of pollutant transport in subsurface porous media. Then the combined model, with the applications of the Eigenvalue technique and the Dynamic system theory, shall be improved to the Dynamic System Model which can simulate the real groundwater flow and the pollutant transport accurately and effectively for the analyses and predictions.

• Korean Journal of Soil Science and Fertilizer
• /
• v.14 no.4
• /
• pp.179-184
• /
• 1982

To describe a mathematical heat transfer model in saturated paddy soils, an analytical solution of the heat flow equation incorporated with the heat transfer by mass flow of water was obtained under the assumptions: 1) the diurnal (or annual) changes in temperature at a depth follow harmonic curves, 2) the temperature at the infinite depth be constant and 3) the temperatures of soil and water at the one depth be identical. The calculation of thermal diffusivities of the soil is possible with the known values of the physical parameters of each component in the soil matrix (heat capacity, density and porosity), percolation rate and the minimum and maximum temperatures at two different depths. The calculated thermal diffusivities using the solution were $9.5cm^2/hr$ for the loam soil with the percolation rate of 0.88cm/day and $13.9cm^2/hr$ for the sandy loam soil with the percolation rate of 2.64 cm/day.

• Geophysics and Geophysical Exploration
• /
• v.19 no.3
• /
• pp.111-120
• /
• 2016

In order to enhance the connectivity of fracture network as fluid path in enhanced/engineered geothermal system (EGS), the exact locating of hydraulic fractured zone is very important. Hydraulic fractures can be tracked by locating of microseismic events which are occurred during hydraulic fracture stimulation at each stage. However, since the subsurface velocity is changed due to hydraulic fracturing at each stage, in order to find out the exact location of microseismic events, we have to consider the velocity change due to hydraulic fracturing at previous stage when we perform the mapping of microseimic events at the next stage. In this study, we have modified 3D locating algorithm of microseismic data which was developed by Kim et al. (2015) and have developed 3D velocity update algorithm using occurred microseismic data. Eikonal equation which can efficiently calculate traveltime for complex velocity model at anywhere without shadow zone is used as forward engine in our inversion. Computational cost is dramatically reduced by using Fresnel volume approach to construct Jacobian matrix in velocity inversion. Through the numerical test which simulates the geothermal survey geometry, we demonstrated that the initial velocity model was updated by using microseismic data. In addition, we confirmed that relocation results of microseismic events by using updated velocity model became closer to true locations.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.27 no.11
• /
• pp.119-125
• /
• 1990

This paper describes fiber optic voltage sensor using EOM-BSO (Electro-Optic Modulator-Bismuth Silicon Oxcide). Transceiver has an electical/optical converter and an optical/electrical converter which consist of light emitting diode, PIN-PD, and electronic circuits. Multimode fiber cable of $100/140{\mu}m$ core/clad diameter is used for connecting the transceiver to fiber cable and fiber optic voltage sensor. Before our experiments, by applying the Maxwell equations and wave equations, We derive matrix equation on wave propagation in the BSO single crystal. And also we derive optimal equation on intensity modulation arising through an analyzer. According to experi-mental results, fiber optic voltage sensor has maximum $2.5{\%}$ error within the applied AC voltage of 800V. As the applied voltage increases, saturation values of voltage sensor also increase. This phenomenon is caused by optical rotatory power of BSO single crystal. And temperature dependence of sensitivity for fiber optical rotatory power of BSO single crystal. And temperature dependence of sensitivity for fiber optic voltage sensor in the temperature range from$-20^{\circ}C\to\60^{\circ}C$ are measured within ${\pm}0.6{\%}$. And frequency characteristics of the voltage sensor has good frequency characteristics from DC to 100kHz.

• Korean Journal of Optics and Photonics
• /
• v.9 no.5
• /
• pp.333-341
• /
• 1998

The novel integrated device, 1.55 ${\mu}{\textrm}{m}$ DR-LD(distrbuted reflector laser diode) integrated EA-MOD (electro-absorption modulator) as light source, is proposed to improve the device yield and its operational performances. This device can be easily fabricated by the selective MOVPE technique and its fabrication processes are almost the same as the reported 1.55 ${\mu}{\textrm}{m}$ DFB-LD(distributed feedback laser diode) integrated EA-MOD except the asymmetric gratings. The static and dynamic properties are investigated simultaneously by solving the transfer matrix method for light propagation, the time-dependent rate equation for carrier change and schr$\"{o}$dinger equation for QCSE (Quantum-Confined Stark Effect). The performances of the proposed device such as output power, chirp, and extinction ratio are compared with those of DFB-LD integrated EA-MOD. Under 10Gb/s NRZ modulation, we obtain that DR-LD integrated EA-MOD. is 30% higher in output power on the on-state, about 50% lower in chirp, and slightly larger in extinction ratio than DFB-LD integrated EA-MOD.-MOD.

• Geophysics and Geophysical Exploration
• /
• v.1 no.2
• /
• pp.127-135
• /
• 1998

The Born approximation is widely used for solving the complex scattering problems in electromagnetics. Approximating total internal electric field by the background field is reasonable for small material contrasts as long as scatterer is not too large and the frequency is not too high. However in many geophysical applications, moderate and high conductivity contrasts cause both real and imaginary part of internal electric field to differ greatly from background. In the extended Born approximation, which can improve the accuracy of Born approximation dramatically, the total electric field in the integral over the scattering volume is approximated by the background electric field projected to a depolarization tensor. The finite difference and elements methods are usually used in EM scattering problems with a 2D model and a 3D source, due to their capability for simulating complex subsurface conductivity distributions. The price paid for a 3D source is that many wavenumber domain solutions and their inverse Fourier transform must be computed. In these differential equation methods, all the area including homogeneous region should be discretized, which increases the number of nodes and matrix size. Therefore, the differential equation methods need a lot of computing time and large memory. In this study, EM modeling program for a 2D model and a 3D source is developed, which is based on the extended Born approximation. The solution is very fast and stable. Using the program, crosshole EM responses with a vertical magnetic dipole source are obtained and the results are compared with those of 3D integral equation solutions. The agreement between the integral equation solution and extended Born approximation is remarkable within the entire frequency range, but degrades with the increase of conductivity contrast between anomalous body and background medium. The extended Born approximation is accurate in the case conductivity contrast is lower than 1:10. Therefore, the location and conductivity of the anomalous body can be estimated effectively by the extended Born approximation although the quantitative estimate of conductivity is difficult for the case conductivity contrast is too high.

• Economic and Environmental Geology
• /
• v.27 no.2
• /
• pp.191-199
• /
• 1994

We have developed an algorithm based on the method of integral equations to simulate the magnetotelluric (MT) responses of three-dimensional (3-D) bodies in a layered half-space. The inhomogeneities are divided into a number of cells and are replaced by an equivalent current distribution which is approximated by pulse basis functions. A matrix equation is constructed using the electric Green's tensor function appropriate to a layered earth, and is solved for the vector current in each cell. Subsequently, scattered fields are found by integrating electric and magnetic Green's tensor functions over the scattering current About a 3-D conductive body near the earth's surface, interpretation using 2-D transverse electric modeling schemes can imply highly erratic low resistivities at depth. This is why these routines do not account for the effect of boundary charges. However, centrally located profiles across elongate 3-D prisms may be modeled accurately with a 2-D transverse magnetic algorithm, which implicitly includes boundary charges in its formulation. Multifrequency calculations show that apparent resistivity and impedance phase are really two complementary parameters. Hence, they should be treated simultaneously in broadband MT interpretation.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.25 no.5
• /
• pp.601-606
• /
• 2014

The IPO(Iterative Physical Optics) method repeatedly applies the well-known PO(Physical Optics) approximation to calculate the scattered field by a large object. Thus, the IPO method can consider the multiple scattering in the object, which is ignored for the PO approximation. This kind of iteration can improve the final accuracy of the induced current on the scatterer, which can result in the enhancement of the accuracy of the RCS(Radar Cross Section) of the scatterer. Since the IPO method can not exactly but approximately solve the required integral equation, however, the convergence of the IPO solution can not be guaranteed. Hence, we apply the famous techniques used in the inversion of a matrix to the IPO method, which include Jacobi, Gauss-Seidel, SOR(Successive Over Relaxation) and Richardson methods. The proposed IPO methods can efficiently calculate the RCS of a large scatterer, and are numerically verified.

• The Journal of the Institute of Internet, Broadcasting and Communication
• /
• v.9 no.1
• /
• pp.89-94
• /
• 2009

This paper proposes the algorithms that design the highly digital logic circuit and assign the code to each node of DCG(Directed Cyclic Graph) of length ${\zeta}$. The conventional algorithm have some problems, so this paper introduce the matrix equation from DCG of length ${\zeta}$ and proposes highly digital logic circuit design algorithms according to the DCG of length ${\zeta}$. Using the proposed circuit design algorithms in this paper, it become realized that was able to design from former algorithm. Also, making a comparison between the circuit using former algorithm and this paper's, we testify that proposed paper's algorithm is able to realize more optimized circuit design. According to proposed circuit design algorithm in this paper, it is possible to design current that DCG have natural number, so it have the following advantages, reduction of the circuit input/output digits, simplification of circuit composition, reduction of computation time and cost. And we show comparability and verification about this paper's algorithm.