• Kyungpook Mathematical Journal
• /
• v.61 no.2
• /
• pp.383-393
• /
• 2021

In this paper, numerical techniques are presented for solving initial value problems of fractional differential equations with variable coefficients. The method is derived by applying a Taylor vector approximation. Moreover, the operational matrix of fractional integration of a Taylor vector is provided in order to transform the continuous equations into a system of algebraic equations. Furthermore, numerical examples demonstrate that this method is applicable and accurate.

• Structural Engineering and Mechanics
• /
• v.14 no.3
• /
• pp.263-285
• /
• 2002

This paper presents a numerical procedure for solving initial-value problems using the special functions which belong to a class of Rvachev's basis functions $R_{bf}$ based on algebraic and trigonometric polynomials. Because of infinite derivability of these functions, derivatives of all orders, required by differential equation of the problem and initial conditions, are used directly in the numerical procedure. The accuracy and stability of the proposed numerical procedure are proved on an example of a single degree of freedom system. Critical time step was also determined. An algorithm for solving multiple degree of freedom systems by the collocation method was developed. Numerical results obtained by $R_{bf}$ functions are compared with exact solutions and results obtained by the most commonly used numerical procedures for solving initial-value problems.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.2
• /
• pp.196-202
• /
• 2012

In this study, numerical analysis of viscous flows is carried out based on the unstructured grid. There exist some difficulties in expressing and computing numerical derivatives on the unstructured grid due to lack of the structured characteristics. The general computer algorithms are developed to perform numerical derivatives easily and extended to be applicable to various geometries composed of hybrid meshes. And the optimal method of strongly implicit procedure is newly contrived to accelerate the rate of convergence in solving the pressure Poisson equation. To verify numerical schemes, the driven cavity problems of 2 and 3 dimension are simulated. The numerical results are compared with others and our numerical schemes are shown to be valid.

• Structural Engineering and Mechanics
• /
• v.49 no.5
• /
• pp.555-568
• /
• 2014

In this article, stability of composite conical shells subjected to dynamic external pressure is investigated by numerical and experimental methods. In experimental tests, cross-ply glass woven fabrics were selected for manufacturing of specimens. Hand-layup method was employed for fabricating the glass-epoxy composite shells. A test-setup that includes pressure vessel and data acquisition system was designed. Also, numerical analyses are performed. In these analyses, effect of actual geometrical imperfections of experimental specimens on the numerical results is investigated. For introducing the imperfections to the numerical models, linear eigen-value buckling analyses were employed. The buckling modes are multiplied by very small numbers that are derived from measurement of actual specimens. Finally, results are compared together while a good agreement between results of imperfect numerical analyses and experimental tests is observed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.341-349
• /
• 1997

The characteristics of defect correction method are discussed in a sample heat conduction problem showing the numerical solution of the error correction equation can predict the error of the numerical solution of the original governing equation. A way of using defect correction method combined with the existing algorithm for the incompressible fluid flow, is proposed and subsequently tested for the driven square cavity problem. The error correction equations for the continuity equation and the momentum equations are considered to estimate the errors of the numerical solutions of the original governing equations. With this new approach, better velocity and pressure fields can be obtained by correcting the original numerical solutions using the estimated errors. These calculated errors also can be used to estimate the orders of magnitude of the errors of the original numerical solutions.

• Journal of electromagnetic engineering and science
• /
• v.3 no.1
• /
• pp.39-44
• /
• 2003

An unconditionally stable compact 2D Alternating-Direction Implicit (ADI) FDTD method for calculating dispersion characteristics of waveguide structures is proposed. The numerical stability and numerical dispersion relation of the proposed method are also presented and discussed. Numerical wavelengths for the dominant and higher order modes in a hollow waveguide are obtained from numerical simulations and compared with those from the analytical dispersion relation. The numerical results show that the proposed scheme has the potential to successfully analyze a class of waveguides having locally fine geometry with reduced numerical costs.

• Journal of Industrial Technology
• /
• v.28 no.B
• /
• pp.185-191
• /
• 2008

This study tried the 1st, 2nd dimensional numerical analysis according to the pier's shape, size and installing method in order to compare the depth of scour calculation method using the variables calculated by using the 2nd dimensional numerical analysis with the calculated depth of scour value by using the calculated variables by using the 1st dimensional numerical analysis. And then verified the problems occurring when the depth of scour is calculated by using the calculated values by using the 1st dimensional numerical analysis, as calculating the deviation depending on it.

• Communications of the Korean Mathematical Society
• /
• v.13 no.3
• /
• pp.655-681
• /
• 1998

Discontinuous solutions or interfaces are common in nature, for examples, shock waves or material interfaces. However, their numerical computation is difficult by the feature of discontinuities. In this paper, we summarize the numerical approaches for discontinuities and interfaces appearing mostly in the system of hyperbolic conservation laws, and explain various numerical methods for them. We explain two numerical approaches to handle discontinuities in the solution: shock capturing and shock tracking, and illustrate their underlying algorithms and mathematical problems. The front tracking method is explained in details and the level set method is outlined briefly. The several applications of front tracking are illustrated, and the research issues in this field are discussed.

• Journal of Mechanical Science and Technology
• /
• v.17 no.7
• /
• pp.986-998
• /
• 2003

In this paper, an adaptive numerical integration scheme, which does not need non-overlapping and contiguous integration meshes, is proposed for the MLPG (Meshless Local Petrov-Galerkin) method. In the proposed algorithm, the integration points are located between the neighboring nodes to properly consider the irregular nodal distribution, and the nodal points are also included as integration points. For numerical integration without well-defined meshes, the Shepard shape function is adopted to approximate the integrand in the local symmetric weak form, by the values of the integrand at the integration points. This procedure makes it possible to integrate the local symmetric weak form without any integration meshes (non-overlapping and contiguous integration domains). The convergence tests are performed, to investigate the present scheme and several numerical examples are analyzed by using the proposed scheme.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.103-106
• /
• 2002

This paper presents an effective numerical method for analyzing three-dimensional unsteady conjugate heat transfer problems of a curved pipe subjected to infernally thermal stratification. In the present numerical analyses, the thermally stratified flows in the pipe are simulated using the standard $k-{\varepsilon}$turbulent model and the unsteady conjugate heat transfer is treated numerically with a simple and convenient numerical technique. The unsteady conjugate heat transfer analysis method is implemented in a finite volume thermal-hydraulic computer code based on a non-staggered grid arrangement, SIMPLEC algorithm and higher-order bounded convection scheme. Numerical calculations have been performed far the two cases of thermally stratified pipe flows where the surging directions are opposite each other i.e. In-surge and out-surge. The results show that the present numerical analysis method is effective to solve the unsteady flow and conjugate heat transfer in a curved pipe subjected to infernally thermal stratification.