• Korean Journal of Optics and Photonics
• /
• v.12 no.6
• /
• pp.472-478
• /
• 2001

We determined the optical constants and thickness of organic electroluminescence thin films using variable angle spectroscopic ellipsometry. Using the measured transmittance spectra and the spectroscopic ellipsomeoy data of the organic films on glass substrates in the optically transparent region, we determined the effective thickness and the refractive indices of organic thin films. Then by applying a numerical inversion method to variable angle spectro-ellipsometry data, we determined the complex refractive index at each wavelength including the optically absorbing region, as well as the thickness and surface micro-roughness of the organic thin films. The calculated transmittance spectra showed a tight agreement with the measured ones, confining the validity of the present model analysis.

• Computers and Concrete
• /
• v.11 no.5
• /
• pp.441-462
• /
• 2013

When using the conventional finite element method, a great number of grid nodes are necessary to describe the large and uneven temperature gradients in the concrete around cooling pipes when calculating the temperature field of mass concrete with cooling pipes. In this paper, the temperature gradient properties of the concrete around a pipe were studied. A new calculation method was developed based on these properties and an explicit iterative algorithm. With a small number of grid nodes, both the temperature distribution along the cooling pipe and the temperature field of the concrete around the water pipe can be correctly calculated with this new method. In conventional computing models, the cooling pipes are regarded as the third boundary condition when solving a model of concrete with plastic pipes, which is an approximate way. At the same time, the corresponding parameters have to be got by expensive experiments and inversion. But in the proposed method, the boundary condition is described strictly, and thus is more reliable and economical. And numerical examples were used to illustrate that this method is accurate, efficient and applicable to the actual engineering.

• The Journal of the Acoustical Society of Korea
• /
• v.38 no.3
• /
• pp.283-289
• /
• 2019

RLS (Recursive-least-squares) algorithm is known to have good convergence and excellent error level after convergence. However, there is a disadvantage that numerical instability is included in the algorithm due to inverse matrix calculation. In this paper, we propose an algorithm with no matrix inversion to avoid the instability aforementioned. The proposed algorithm still keeps the same convergence performance. In the proposed algorithm, we adopt an averaged gradient-based step size as a self-adjusted step size. In addition, a variable forgetting factor is introduced to provide superior performance for time-varying channel estimation. Through simulations, we compare performance with conventional RLS and show its equivalency. It also shows the merit of the variable forgetting factor in time-varying channels.

• Journal of the Korean earth science society
• /
• v.36 no.4
• /
• pp.315-329
• /
• 2015

We investigated the characteristics of meteorological conditions related to the strong downslope wind over the leeward side of the Taebaek Mountains during the period 2005~2010. The days showing the strong wind exceeding $14ms^{-1}$ in Gangwon province were selected as study cases. A total of 15 days of strong wind were observed at Sokcho, Gangneung, Donghae, and Taebaek located over the Yeongdong region. Seven cases related to tropical cyclone (3 cases) and heavy snowfall (2 cases) and heavy rainfall (2 cases) over the Yeongdong region were excluded. To investigate the characteristics of the remaining 8 cases, we used synoptic weather chart, Sokcho radiosonde, Gangneung wind profiler and numerical model. The cases showed no precipitation (or ${\leq}1mm\;day^{-1}$). From the surface and upper level weather chart, we found the pressure distribution of southern high and northern low pattern over the Korean peninsula and warm ridge over the Yeongdong region. Inversion layer (or stable layer) and warm ridge with strong wind were located in about 1~3 km (925~700 hPa) over mountains. The Regional Data Assimilation and Prediction System (RDAPS) indicated that warm core and temperature ridge with horizontal temperature gradient were $0.10{\sim}0.23^{\circ}C\;km^{-1}$ which were located on 850 hPa pressure level above mountaintop. These results were summarized as a forecasting guidance of downslope windstorm in the Yeongdong region.

• Journal of Environmental Science International
• /
• v.13 no.8
• /
• pp.721-731
• /
• 2004

The dispersion of suspended particulates in the coastal complex terrain of mountain-inland basin (city)-sea, considering their recycling was investigated using three-dimensional non-hydrostatic numerical model and lagrangian particle model (or random walk model). Convective boundary layer under synoptic scale westerly wind is developed with a thickness of about I km over the ground in the west of the mountain, while a thickness of thermal internal boundary layer (TIBL) is only confined to less than 200m along the eastern slope of the mountain, below an easterly sea breeze circulation. At the mid of the eastern slop of the mountain, westerly wind confronts easterly sea breeze, which goes to the height of 1700 m above sea level and is finally eastward return flow toward the sea. At this time, particulates floated from the ground surface of the city to the top of TIBL go along the eastern slope of the mountain in the passage of sea breeze, being away the TIBL and reach near the top of the mountain. Then those particulates disperse eastward below the height of sea-breeze circulation and widely spread out over the coastal sea. Total suspended particulate concentration near the ground surface of the city is very low. On the other hand, nighttime radiative cooling produces a shallow nocturnal surface inversion layer (NSIL) of 200 m thickness over the inland surface, but relatively thin thickness less than 100m is found near the mountain surface. As synoptic scale westerly wind should be intensified under the association of mountain wind along the eastern slope of mountain to inland plain and further combine with land-breeze from inland plain toward sea, resulting in strong wind as internal gravity waves with a hydraulic jump motion bounding up to about 1km upper level in the atmosphere in the west of the city and becoming a eastward return flow. Simultaneously, wind near the eastern coastal side of the city was moderate. Since the downward strong wind penetrated into the city, the particulate matters floated near the top of the mountain in the day also moved down along the eastern slope of the mountain, reaching the. downtown and merging in the ground surface inside the NSIL with a maximum ground level concentration of total suspended particulates (TSP) at 0300 LST. Some of them were bounded up from the ground surface to the 1km upper level and the others were forward to the coastal sea surface, showing their dispersions from the coastal NSIL toward the propagation area of internal gravity waves. On the next day at 0600 LST and 0900 LST, the dispersed particulates into the coastal sea could return to the coastal inland area under the influence of sea breeze and the recycled particulates combine with emitted ones from the ground surface, resulting in relatively high TSP concentration. Later, they float again up to the thermal internal boundary layer, following sea breeze circulation.

• 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 Computational Structural Engineering Institute of Korea
• /
• v.32 no.4
• /
• pp.249-256
• /
• 2019

This paper presents an optimization framework of electrical impedance tomography for characterizing electrical conductivity profiles of concrete structures in two dimensions. The framework utilizes a partial-differential-equation(PDE)-constrained optimization approach that can obtain the spatial distribution of electrical conductivity using measured electrical potentials from several electrodes located on the boundary of the concrete domain. The forward problem is formulated based on a complete electrode model(CEM) for the electrical potential of a medium due to current input. The CEM consists of a Laplace equation for electrical potential and boundary conditions to represent the current inputs to the electrodes on the surface. To validate the forward solution, electrical potential calculated by the finite element method is compared with that obtained using TCAD software. The PDE-constrained optimization approach seeks the optimal values of electrical conductivity on the domain of investigation while minimizing the Lagrangian function. The Lagrangian consists of least-squares objective functional and regularization terms augmented by the weak imposition of the governing equation and boundary conditions via Lagrange multipliers. Enforcing the stationarity of the Lagrangian leads to the Karush-Kuhn-Tucker condition to obtain an optimal solution for electrical conductivity within the target medium. Numerical inversion results are reported showing the reconstruction of the electrical conductivity profile of a concrete specimen in two dimensions.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.12
• /
• pp.15-24
• /
• 2006

The evaluation of shear modulus is very important in various fields of civil engineering. Non-destructive seismic methods can be used to determine shear wave velocity ($V_s$) profile. Non-destructive seismic methods geneally consist of three steps: field testing, evaluation of dispersion curve, and determination of Vs profile by inversion process. Non-destructive seismic methods can be divided into two categories according to the number of receivers used for data reduction: two-channel tests and multi-channel tests. Two channel tests use apparent velocity dispersion curve and multi-channel tests use mode dispersion curves. Multi-channel tests using mode dispersion curve can reduce calculation time to determine soil profile and uncertainties in inversion process. So far, only multi-channel tests can determine mode dispersion curves but multi-channel test needs many receivers to determine reasonable mode dispersion curves. In this paper, HWAW (Harmonic Wavelet Analysis of Wave) method is applied to determine mode dispersion curves. HWAW method uses short test setup which consists of two receivers with a spacing of 1 to 3 m. Through numerical simulations and field application, it is shown that HWAW can determine resonable mode disperson curves.

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1027-1036
• /
• 2023

Coastal topographic information is crucial in coastal management, but point measurment based approeaches, which are labor intensive, are generally applied to land and underwater, separately. This study introduces an efficient method enabling land and undetwater surveys using an unmanned aerial vehicle (UAV). This method involves applying two different algorithms to measure the topography on land and water depth, respectively, using UAV imagery and merge them to reconstruct whole coastal digital elevation model. Acquisition of the landside terrain is achieved using the Structure-from-Motion Multi-View Stereo technique with spatial scan imagery. Independently, underwater bathymetry is retrieved by employing a depth inversion technique with a drone-acquired wave field video. After merging the two digital elevation models into a local coordinate, interpolation is performed for areas where terrain measurement is not feasible, ultimately obtaining a continuous nearshore terrain. We applied the proposed survey technique to Jangsa Beach, South Korea, and verified that detailed terrain characteristics, such as berm, can be measured. The proposed UAV-based survey method has significant efficiency in terms of time, cost, and safety compared to existing methods.

• Geophysics and Geophysical Exploration
• /
• v.27 no.3
• /
• pp.181-195
• /
• 2024

This tutorial explains that the static effect in the magnetotelluric (MT) survey is a physical phenomenon caused by charges accumulated on the boundaries of subsurface inhomogeneities. To facilitate understanding of the physical phenomenon, differences between static induction and charge accumulation on the boundary are explained and analyzed with help of schematic illustrations. Subsequently, from the electromagnetic (EM) integral equation formulation, it is clearly shown that the secondary electric field due to charges accumulated on the interface in the presence of the primary field appears as the static effect. Therefore, except in the cases of the layered earth or a two-dimensional earth with transverse magnetic (TM) mode excitation, the static effect always exists in MT responses and further, it is not 'static' but rather frequency dependent. Despite the fact that the static effect is a secondary electric field due to inhomogeneity, inevitable under-sampling in the frequency and spatial domains prevent the effect from being handled properly in numerical inversion. Therefore, considering the practical aspects of the MT survey, which cannot be a continuous measurement covering the entire survey area over a wide frequency band, a three-dimensional (3-D) inversion incorporating the static shift as a constraint with the Gaussian distribution is introduced. To enhance understanding of the integral equation EM modeling, the formulation of the 3-D integral equation and mathematical analyses of the Green tensor and scattering current are described in detail in the Appendix.