• Korean Journal of Hydrosciences
• /
• v.5
• /
• pp.85-97
• /
• 1994

A hybrid finite difference method for the longitudinal dispersion equation, which is based on combining the Holly-Preissmann scheme with fifth-degree Hermite interpolating polynomial and the generalized Crank-Nicholson scheme, is described and comparatively evaluated with other characteristics-based numerical methods. Longitudinal dispersion of an instantaneously-loaded pollutant source is simulated, and computational results are compared with the exact solution. The present method is free from wiggles regardless of the Courant number, and exactly reproduces the location of the peak concentration. Overall accuracy of the computation increases for smaller value of the weighting factor, $\theta$of the model. Larger values of $\theta$ overestimates the peak concentration. Smaller Courant number yields better accuracy, in general, but the sensitivity is very low, especially when the value of $\theta$ is small. From comparisons with the hybrid method using cubic interpolating polynomial and with splitoperator methods, the present method shows the best performance in reproducing the exact solution as the advection becomes more dominant.

• Journal for History of Mathematics
• /
• v.32 no.4
• /
• pp.159-173
• /
• 2019

In this paper, we examine the brief history of the ring of four almonds regarding Mesopotamian mathematics, and present reasons why the Omar Khayyam's triangle, a special right triangle in a ring of four almonds, was essential for artisans due to its unique pattern. We presume that the ring of four almonds originated from a point symmetry figure given two concentric squares used in the proto-Sumerian Jemdet Nasr period (approximately 3000 B.C.) and a square halfway between two given concentric squares used during the time of the Old Akkadian period (2340-2200 B.C.) and the Old Babylonian age (2000-1600 B.C.). Artisans tried to create a new intricate pattern as almonds and 6-pointed stars by subdividing right triangles in the pattern of the popular altered Old Akkadian square band at the time. Therefore, artisans needed the Omar Khayyam's triangle, whose hypotenuse equals the sum of the short side and the perpendicular to the hypotenuse. We presume that artisans asked mathematicians how to construct the Omar Khayyam's triangle at a meeting between artisans and mathematicians in Isfahan. The construction of Omar Khayyam's triangle requires solving an irreducible cubic polynomial. Omar Khayyam was the first to classify equations of integer polynomials of degree up to three and then proceeded to solve all types of cubic equations by means of intersections of conic sections. Omar Khayyam's triangle gave practical meaning to the type of cubic equation $x^3+bx=cx^2+a$. The work of Omar Khayyam was completed by Descartes in the 17th century.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.20 no.3
• /
• pp.265-272
• /
• 2002

By this time, many methods have been developed for computing the pit excavation volumes, ranging from a simple formula to more complicated numerical methods. Earlier the standard methods for pit excavation volume computation requires that the considered area be divided the boundary ranges of x and y directions into a rectangular grid. whereas these methods may not calculate the estimation of pit excavation volume that is often required in many surveying situation exactly. In Easa methods(1998), the rectangular grid is divided into the same linear in the range x and y directions respectively. This method employs a cubic Hermite polynomial for individual intervals in both directions of the grid. Because the height data over the same boundary of x and y interval ranges have to be exist, it is not possible to choose the governing points of the terrain boundary such as points of maximum and minimum height. In this study, a method of volume computation, that combines the advantages of Easa methods(1998) and avoids the drawbacks of it, is presented. The proposed method employs a cubic Hermite polynomial for individual intervals in both directions of the non-grid, the all over intervals of it may be unequal grid x in width and y in length y, partially. The new proposed method should produce better accuracy than the other conventional methods.

• Journal of the Korean Geosynthetics Society
• /
• v.15 no.1
• /
• pp.1-10
• /
• 2016

Recently, offshore wind farms are increasingly expected, because there are huge resource and large site in offshore. Jeju island has optimum condition for constructing a wind energy farm. Unlike the mainland, Jeju island has stratified structure distribution between rock layers sediments due to volcanic activation. In these case, it can be occur engineering problems in whole structures as well as the safety of foundation as the thickness and distribution of sediment under top rock layer can not support sufficiently the structure. In this study, field lateral load test of the pile for analyzing lateral behavior of the offshore wind turbine which is embedded in basalt. After calculating the subgrade resistance and the horizontal deflection from the measured strain to derive p-y curve from the lateral load test results, the subgrade resistance amplifies the error in the process of differentiation and the error of piecewise polynomial curve fitting is the smallest. In order to calculate the horizontal deflection from the measured strain, the six-order polynomial was used.

• Structural Engineering and Mechanics
• /
• v.38 no.2
• /
• pp.211-229
• /
• 2011

In this paper, we present a new explicit procedure using periodic cubic B-spline interpolation polynomials to solve linear and nonlinear dynamic equation of motion governing single degree of freedom (SDOF) systems. In the proposed approach, a straightforward formulation was derived from the approximation of displacement with B-spline basis in a fluent manner. In this way, there is no need to use a special pre-starting procedure to commence solving the problem. Actually, this method lies in the case of conditionally stable methods. A simple step-by-step algorithm is implemented and presented to calculate dynamic response of SDOF systems. The validity and effectiveness of the proposed method is demonstrated with four examples. The results were compared with those from the numerical methods such as Duhamel integration, Linear Acceleration and also Exact method. The comparison shows that the proposed method is a fast and simple procedure with trivial computational effort and acceptable accuracy exactly like the Linear Acceleration method. But its power point is that its time consumption is notably less than the Linear Acceleration method especially in the nonlinear analysis.

• Current Optics and Photonics
• /
• v.7 no.5
• /
• pp.518-528
• /
• 2023

As an effective means of remotely detecting the spectral information of the object, the spectral calibration for the Savart polarization interference imaging spectrometer (SPIIS) is a basis and prerequisite of information quantification, and its experimental calibration scheme is firstly proposed in this paper. In order to evaluate the accuracy of the spectral information acquisition, the linear interpolation, cubic spline interpolation, and piecewise cubic interpolation algorithms are adopted, and the precision of the quadratic polynomial fitting is the highest, whose fitting error is better than 5.8642 nm in the wavelength range of [500 nm, 820 nm]. Besides, the inversed value of the spectral resolution for the monochromatic light is greater than the theoretical value, and the deviation between them becomes larger with the wavelength increasing, which is mainly caused by the structural design of the SPIIS, together with the rationality of the spectral restoration algorithm and the selection of the maximum optical path difference (OPD). This work demonstrates that the SPIIS has achieved high performance assuring the feasibility of its practical use in various fields.

• Journal of Navigation and Port Research
• /
• v.26 no.4
• /
• pp.473-479
• /
• 2002

The estimation of the volume of a pit excavation is often required in many surveying, soil mechanics, highway applications and transportation engineering situations. The calculation of earthwork plays a major role in plan or design of many civil engineering projects such as seashore reclamation, and thus it has become very important to improve the accuracy of earthwork calculation. In this paper the spot height method, proposed formulas(A, B, C), and chen and Line method are compared with the volumes of the pits in these examples. And we proposed an algorithm of finding a terrain surface with the free boundary conditions and both direction spline method drawback, i.e., the modeling curves form peak points at the joints. To avoid this drawback, the cubic spline polynomial was chosen as the methematical model of the new method. From the characteristics of the cubic spline polynomial, the modeling curve of the new method was smooth and matched the ground profile well. As a result of this study, algorithm of proposed three methods to estimate pit excavation volume provided a better accuracy than spot height, chamber, chen and Lin method. And the mathematical model mentioned makes is thought to give a maximum acccuracy in estimating the volume of a pit excavation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.3
• /
• pp.271-280
• /
• 2004

Assessment of image registration for Pressure Sensitive Paint (PSP) was performed. A 16 bit camera and LED lamp were used with Uni-FIB paint (ISSI). Because of model displacement and deformation at 'wind-on' condition, a large error of the intensity ratio was induced between 'wind-on' and' wind-off images. To correct the error, many kinds of image registrations were tested. At first, control points were marked on the model surface to find the coefficients of polynomial transform functions between the 'wind-off' 'wind-on' images. The 2nd-order polynomial function was sufficient for representing the model displacement and deformation. An automatic detection scheme was introduced to find the exact coordinates of the control points. The present automatic detection algorithm showed more accurate and user-friendly than the manual detection algorithm. Since the coordinates of transformed pixel were not integer, five interpolation methods were applied to get the exact pixel intensity after transforming the 'wind-on' image. Among these methods, the cubic convolution interpolation scheme gave the best result.

• Journal of the Korean Society for Precision Engineering
• /
• v.23 no.9 s.186
• /
• pp.200-207
• /
• 2006

In this paper a new approach which combines implicit surface scheme and NURBS surface interpolation method is proposed in order to generate a complete surface model from unorganized point cloud data. In the method a base surface was generated by creating smooth implicit surface from the input point cloud data through which the actual surface would pass. The implicit surface was defined by a combination of shape functions including quadratic polynomial function, cubic polynomial functions and radial basis function using adaptive domain decomposition method. In this paper voxel data which can be extracted easily from the base implicit surface were used in order to generate rectangular net with good quality using the normal projection and smoothing scheme. After generating the interior points and tangential vectors in each rectangular region considering the required accuracy, the NURBS surface were constructed by interpolating the rectangular array of points using boundary tangential vectors which assure C$^1$ continuity between rectangular patches. The validity and effectiveness of this new approach was demonstrated by performing numerical experiments for the various types of point cloud data.

• Journal of Information Processing Systems
• /
• v.10 no.3
• /
• pp.395-411
• /
• 2014

Temporal medical data is often collected during patient treatments that require personal analysis. Each observation recorded in the temporal medical data is associated with measurements and time treatments. A major problem in the analysis of temporal medical data are the missing values that are caused, for example, by patients dropping out of a study before completion. Therefore, the imputation of missing data is an important step during pre-processing and can provide useful information before the data is mined. For each patient and each variable, this imputation replaces the missing data with a value drawn from an estimated distribution of that variable. In this paper, we propose a new method, called Newton's finite divided difference polynomial interpolation with condition order degree, for dealing with missing values in temporal medical data related to obesity. We compared the new imputation method with three existing subspace estimation techniques, including the k-nearest neighbor, local least squares, and natural cubic spline approaches. The performance of each approach was then evaluated by using the normalized root mean square error and the statistically significant test results. The experimental results have demonstrated that the proposed method provides the best fit with the smallest error and is more accurate than the other methods.