• Journal of the Korean Mathematical Society
• /
• v.40 no.6
• /
• pp.1075-1083
• /
• 2003

The purpose of this paper is to investigate the piecewise wise polynomial approximation for the curved boundary. We analyze the error of an approximated solution due to this approximation and then compare the approximation errors for the cases of polygonal and piecewise polynomial approximations for the curved boundary. Based on the results of analysis, p-version numerical methods for solving Dirichlet's problems are applied to any smooth curved domain.

• Journal of applied mathematics & informatics
• /
• v.8 no.1
• /
• pp.27-44
• /
• 2001

A fast Poisson solver on irregular domains, based on bound-ary methods, is presented. The harmonic polynomial approximation of the solution of the associated homogeneous problem provides a good practical boundary method which allows a trivial parallel processing for solution evaluation or straightfoward computations of the interface values for domain decomposition/embedding. AMS Mathematics Subject Classification : 65N35, 65N55, 65Y05.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.5 no.1
• /
• pp.47-53
• /
• 2001

A problem of approximation by irrational splines is considered. These splines have a constant curvature between interpolation nodes and need only one additional boundary condition for derivatives, which should be set only at one of two boundary nodes, that is impossible for usual polynomial splines required boundary conditions at both boundary nodal points. Some estimations for numerical differentiation and rounding error analysis are presented.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.2
• /
• pp.67-85
• /
• 2013

Reproducing Polynomial Particle Method (RPPM) is one of meshless methods that use meshes minimally or do not use meshes at all. In this paper, the RPPM is employed for free vibration analysis of shear-deformable plates of the first order shear deformation model (FSDT), called Reissner-Mindlin plate. For numerical implementation, we use flat-top partition of unity functions, introduced by Oh et al, and patchwise RPPM in which approximation functions have high order polynomial reproducing property and the Kronecker delta property. Also, we demonstrate that our method is highly effective than other existing results for various aspect ratios and boundary conditions.

• 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 Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.411-420
• /
• 2009

This study presents an extended finite difference method based on moving least squares(MLS) method for solving potential problems with interfacial boundary. The approximation constructed from the MLS Taylor polynomial is modified by inserting of wedge functions for the interface modeling. Governing equations are node-wisely discretized without involving element or grid; immersion of interfacial condition into the approximation circumvents numerical difficulties owing to geometrical modeling of interface. Interface modeling introduces no additional unknowns in the system of equations but makes the system overdetermined. So, the numbers of unknowns and equations are equalized by the symmetrization of the stiffness matrix. Increase in computational effort is the trade-off for ease of interface modeling. Numerical results clearly show that the developed numerical scheme sharply describes the wedge behavior as well as jumps and efficiently and accurately solves potential problems with interface.

• Journal of Ocean Engineering and Technology
• /
• v.36 no.2
• /
• pp.101-107
• /
• 2022

Dean (1965) proposed the use of the root mean square error (RMSE) in the dynamic free surface boundary condition (DFSBC) and kinematic free-surface boundary condition (KFSBC) as an error evaluation criterion for wave theories. There are well known wave theories with RMSE more than 1%, such as Airy theory, Stokes theory, Dean's stream function theory, Fenton's theory, and trochodial theory for deep-water waves. However, none of them can be applied for deep-water breaking waves. The purpose of this study is to provide a closed-form solution for deep-water waves with RMSE less than 1% even for breaking waves. This study is based on a previous study (Shin, 2016), and all flow fields were simplified for deep-water waves. For a closed-form solution, all Fourier series coefficients and all related parameters are presented with Newton's polynomials, which were determined by curve fitting data (Shin, 2016). For verification, a wave in Miche's limit was calculated, and, the profiles, velocities, and the accelerations were compared with those of 5^th-order Stokes theory. The results give greater velocities and acceleration than 5^th-order Stokes theory, and the wavelength depends on the wave height. The results satisfy the Laplace equation, bottom boundary condition (BBC), and KFSBC, while Stokes theory satisfies only the Laplace equation and BBC. RMSE in DFSBC less than 7.25×10^-2% was obtained. The series order of the proposed method is three, but the series order of 5^th-order Stokes theory is five. Nevertheless, this study provides less RMSE than 5^th-order Stokes theory. As a result, the method is suitable for offshore structural design.

• Coupled systems mechanics
• /
• v.1 no.2
• /
• pp.155-163
• /
• 2012

This paper presents analytical approximate solutions for the initial post-buckling deformation of the pipes in oil and gas wells. The governing differential equation with sinusoidal nonlinearity can be reduced to form a third-order-polynomial nonlinear equation, by coupling of the well-known Maclaurin series expansion and orthogonal Chebyshev polynomials. Analytical approximations to the resulting boundary condition problem are established by combining the Newton's method with the method of harmonic balance. The linearization is performed prior to proceeding with harmonic balancing thus resulting in a set of linear algebraic equations instead of one of non-linear algebraic equations, unlike the classical method of harmonic balance. We are hence able to establish analytical approximate solutions. The approximate formulae for load along axis, and periodic solution are established for derivative of the helix angle at the end of the pipe. Illustrative examples are selected and compared to "reference" solution obtained by the shooting method to substantiate the accuracy and correctness of the approximate analytical approach.

• Structural Engineering and Mechanics
• /
• v.43 no.2
• /
• pp.211-223
• /
• 2012

This paper deals with large deformation post-buckling of a linear-elastic and hygrothermal beam with axially nonmovable pinned-pinned ends and subjected to a significant increase in swelling by an alternative method. Analytical approximate solutions for the geometrically nonlinear problem are presented. The solution for the limiting case of a string is also obtained. By coupling of the well-known Maclaurin series expansion and orthogonal Chebyshev polynomials, the governing differential equation with sinusoidal nonlinearity can be reduced to form a cubic-nonlinear equation, and supplementary condition with cosinoidal nonlinearity can also be simplified to be a polynomial integral equation. Analytical approximations to the resulting boundary condition problem are established by combining the Newton's method with the method of harmonic balance. Two approximate formulae for load along axis, potential strain for free hygrothermal expansion and periodic solution are established for small as well as large angle of rotation at the end of the beam. Illustrative examples are selected and compared to "reference" solution obtained by the shooting method to substantiate the accuracy and correctness of the approximate analytical approach.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.495-500
• /
• 2007

A Meshfree is a method used to establish algebraic equations of system for the whole problem domain without the use of a predefined mesh for the domain discretization. A point interpolation method is based on combining radial and polynomial basis functions. Involvement of radial basis functions overcomes possible singularity. Furthermore, the interpolation function passes through all scattered points in an influence domain and thus shape functions are of delta function property. This makes the implementation of essential boundary conditions much easier than the meshfree methods based on the moving least-squares approximation. This study aims to investigate a stress analysis of structural element between a meshfree method and the finite element method. Examples on cantilever type plate and stress concentration problems show that the accuracy and convergence rate of the meshfree methods are high.