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

• KIPS Transactions on Software and Data Engineering
• /
• v.6 no.11
• /
• pp.527-536
• /
• 2017

With the recent development of smart grid industry, the necessity for efficient EMS(Energy Management System) has been increased. In particular, in order to reduce electric load and energy cost, sophisticated electric load forecasting and efficient smart grid operation strategy are required. In this paper, for more accurate electric load forecasting, we extend the data collected at demand time into high time resolution and construct an artificial neural network-based forecasting model appropriate for the high time resolution data. Furthermore, to improve the accuracy of electric load forecasting, time series data of sequence form are transformed into continuous data of two-dimensional space to solve that problem that machine learning methods cannot reflect the periodicity of time series data. In addition, to consider external factors such as temperature and humidity in accordance with the time resolution, we estimate their value at the time resolution using linear interpolation method. Finally, we apply the PCA(Principal Component Analysis) algorithm to the feature vector composed of external factors to remove data which have little correlation with the power data. Finally, we perform the evaluation of our model through 5-fold cross-validation. The results show that forecasting based on higher time resolution improve the accuracy and the best error rate of 3.71% was achieved at the 3-min resolution.

• Journal of Korean Society for Geospatial Information Science
• /
• v.23 no.1
• /
• pp.89-99
• /
• 2015

Recently, our construction industry is actively participating in numerous city planning projects in the third world countries. Considering the current depression of domestic real estate market, the emerging foreign demands could certainly provide substantial opportunities for the domestic industry to overcome the trough. For the field planners dealing with such foreign projects, though, the immediate problem is the lack of public statistics and geographic information to perform spatial analyses and/or prepare master plans. This study, in this context, tries to simulate a process to construct a digitized basemap of the case area, 'Bab Ezzouar,' in Algeria of Northern Africa. The area is a typical municipality that lacks the IT databases. To overcome the data shortage, the study uses the satellite map tiles so as to digitize the roads and building structures. It then estimates the block-wise populations based on the building image interpolation as well as the supplementary field survey data. The topographic TINs are also built by the SRTM (Shuttle Radar Topography Mission) digital elevation maps so that the three-dimensional configuration of the structures and terrains are rendered to check the urban scenery and skylines.

• Geophysics and Geophysical Exploration
• /
• v.16 no.1
• /
• pp.27-35
• /
• 2013

In this study, we developed the trial product on 3D magnetic inversion tentatively named 'KMag3D'. Also, we briefly introduced its own function and graphic user interface on which especially focused through the development in the form of user manual. KMag3D is consisted of two fundamental frame for the 3D magnetic inversion. First, algebraic reconstruction technique was selected as a 3D inversion algorithm instead of least square method conventionally used in various magnetic inversion. By comparison, it was turned out that algebraic reconstruction algorithm was more effective and economic than that of least squares in aspect of both computation time and memory. Second, for the effective determination of the 3D initial and a-priori information model required in the execution of our algorithm, we proposed the practical technique based on the assemblage of 2D forward modeling and inversion results for individual user-selected 2D profiles. And in succession, initial and a-priori information model were constructed by appropriate interpolation along the strke direction. From this, we concluded that our technique is both suitable and very practical for the application of 3D magentic inversion problem.

• Journal of Navigation and Port Research
• /
• v.46 no.6
• /
• pp.471-482
• /
• 2022

This study aims to investigate the characteristics of wave generation in a deep ocean engineering basin and to develop a meta-model of the transfer function of the wavemaker that reflects the geometric characteristics of the deep ocean engineering basin. To this end, the two-dimensional frequency domain boundary element method was applied to achieve an efficient analysis that reflects the geometric characteristics of the deep ocean engineering basin. The developed numerical method was validated through comparison with the analytical solution. Numerical analyses were conducted for the boundary value problem of the wavemaker according to various periods and the positions of the movable bottom. The numerical results were used to investigate the effect of the geometric characteristics of the deep ocean engineering basin on the transfer function of the wavemaker, and the effect of depth on wave generation was checked by changing the position of the movable bottom. To efficiently utilize the various results of the boundary element method, a meta-model, an approximate model of the transfer function of the wave maker, was developed using a thin plate spline interpolation model. The validity of the developed meta-model was confirmed through a comparison of the results of the model tests.

• The Journal of the Acoustical Society of Korea
• /
• v.13 no.6
• /
• pp.75-82
• /
• 1994

The ray paths and travel times of sound wave in the ocean depend on the physical properties of the propagating media. Ocean Acoustic Tomography(OAT), which is inversely estimate the travel time variations between fixed sources and receivers the physical properties of the corresponding media can he understood. To apply ocean survey technology by using the OAT, the tomographic procedure requires forward problem that variation of the travel times be identified with the variability of the medium. Also, received signals must be satisfied the necessary conditions of ray path stability, identification and resolution in order for OAT to work. The canonical ocean has been determined based on the historical data and its travel time and ray path are used as reference values. The sound speed of canonical ocean in the East Sea is about 1523 m/s at the surface and 1458 m/s at the sound channel axis(400m). Sound speeds in the East Sea are perturbed by warm eddy whose horizontal extension is more than 100 km with deeper than 200 m in depth scale. In this study, an acoustic source and receiver are placed at the depth above the sound channel axis, 350 m, and are separated by 200 km range. Ray paths are identified by the ray theory methed in a range dependent medium whose sound speeds are functions of a range and depth. The eigenray information obtained from interpolation between the rays bracketing the receiver are used to simulate the received signal by convolution of source signal with the eigenray informations. The source signal is taken as a 400 Hz rectangular pulse signal, bandwidth is 16 Hz and pulse length is 64 ms. According to the analysis of the received signal and identified ray path by using numerical model of underwater sound propagation, simulated signals satisfy the necessary conditions of OAT, applied in the East Sea.