• Journal of applied mathematics & informatics
• /
• v.30 no.1_2
• /
• pp.235-244
• /
• 2012

In this paper, we solve the three-dimensional biharmonic equation with Dirichlet boundary conditions of second kind using the full multigrid (FMG) algorithm. We derive a finite difference approximations for the biharmonic equation on a 18 point compact stencil. The unknown solution and its second derivatives are carried as unknowns at grid points. In the multigrid methods, we use a fourth order interpolation to producing a new intermediate unknown functions values on a finer grid, and the full weighting restriction operators to calculating the residuals at coarse grid points. A set of test problems gives excellent results.

• Journal of the Chungcheong Mathematical Society
• /
• v.7 no.1
• /
• pp.237-242
• /
• 1994

An quintic spline interpolate to a function in $C^{10}$[a, b] and its O($h^6$) error behavior are presented when its fourth derivative satisfies some kind of end conditions. The O($h^6$) relations between its derivatives up to fourth order and the m-th derivatives of the given function are also given at the nodes.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.33 no.5
• /
• pp.279-286
• /
• 2020

In this study, we investigate the accuracy of higher order derivatives in the moving least square (MLS) difference method. An interpolation function is constructed by employing a Taylor series expansion via MLS approximation. The function is then applied to the mixed variational theorem in which the displacement and stress resultants are treated as independent variables. The higher order derivatives are evaluated by solving simply supported beams and cantilevers. The results are compared with the analytical solutions in terms of the order of polynomials, support size of the weighting function, and number of nodes. The accuracy of the higher order derivatives improves with the employment of the mean value theorem, especially for very high-order derivatives (e.g., above fourth-order derivatives), which are important in a classical asymptotic analysis.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.8
• /
• pp.693-698
• /
• 2007

A numerical method for accurate simulations of rotor aerodynamics is proposed. The numerical diffusion in the typically coarse grids away from the rotor blades is improved by implying a fourth-order of interpolation of local characteristic variables of the flow in the reconstruction stage of MUSCL approach in the framework of a finite volume formulation. In addition, different slope limiters are applied to the different characteristic fields, such as compressive limiters to linear characteristic fields to reduce the numerical dissipation whereas, diffusive limiters to nonlinear characteristic fields to increase numerical stability. Various exemplary problems related to the rotor aerodynamics applications are tested and the numerical results show a significant improvement in wake capturing capability. However, rotor aeroacoustic calculations show no meaningful difference over traditional MUSCL approach.

• Structural Engineering and Mechanics
• /
• v.68 no.4
• /
• pp.385-398
• /
• 2018

In the present paper, we offer a new flat shell finite element. It is the result of the combination of a membrane element and a bending element, both based on the strain-based formulation. It is known that $C^{\circ}$ plane membrane elements provide poor deflection and stress for problems where bending is dominant. In addition, they encounter continuity and compliance problems when they connect to C1 class plate elements. The reach of the present work is to surmount these problems when a membrane element is coupled with a thin plate element in order to construct a shell element. The membrane element used is a triangular element with four nodes, three nodes at the vertices of the triangle and the fourth one at its barycenter. Each node has three degrees of freedom, two translations and one rotation around the normal. The coefficients related to the degrees of freedom at the internal node are subsequently removed from the element stiffness matrix by using the static condensation technique. The interpolation functions of strain, displacements and stresses fields are developed from equilibrium conditions. The plate element used for the construction of the present shell element is a triangular four-node thin plate element based on Kirchhoff plate theory, the strain approach, the four fictitious node, the static condensation and the analytic integration. The shell element result of this combination is robust, competitive and efficient.