• 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.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• Geophysics and Geophysical Exploration
• /
• v.6 no.4
• /
• pp.165-170
• /
• 2003

Recently, resistivity surveys have been frequently carried out over the irregular terrain such as mountainous area. Such an irregular terrain itself can produce significant anomalies which may lead to misinterpretations. In this study, topographic effects in resistivity survey were studied using the physical scale modeling as well as the numerical one adopting finite element method. The scale modeling was conducted at a pond, so that we could avoid the edge effect, the inherent problem of the scale modeling conducted in a water tank in laboratory. The modeling experiments for two topographic features, a ridge and a valley with various slope angles, confirmed that the results by the two different modeling techniques coincide with each other fairly well for all the terrain models. These experiments adopting dipole-dipole array showed the distinctive terrain effects, such that a ridge produces a high apparent resistivity anomaly at the ridge center flanked by zones of lower apparent resistivity. On the other hand, a valley produces the opposite anomaly pattern, a central low flanked by highs. As the slope of a terrain model becomes steeper, the terrain-induced anomalies become stronger, and moreover, apparent resistivity can become even negative for the model with extremely high slope angle. All the modeling results led us to the conclusion that terrain effects should be included in the numerical modeling and/or the inversion process to interpret data acquired at the rugged terrain area.