• International Journal of Naval Architecture and Ocean Engineering
• /
• 제7권6호
• /
• pp.1034-1043
• /
• 2015

This is the continuation of a numerical study on vortex shedding from a blunt trailing-edge of a hydrofoil. In our previous work (Lee et al., 2015), numerical schemes for efficient computations were successfully implemented; i.e. multiple domains, the approximation of domain boundary conditions using cubic spline functions, and particle-based domain decomposition for better load balancing. In this study, numerical results through a hybrid particle-mesh method which adopts the Vortex-In-Cell (VIC) method and the Brinkman penalization model are further rigorously validated through comparison to experimental data at the Reynolds number of $2{\times}10^6$. The effects of changes in numerical parameters are also explored herein. We find that the present numerical method enables us to reasonably simulate vortex shedding phenomenon, as well as turbulent wakes of a hydrofoil.

• Journal of computational fluids engineering
• /
• 제12권3호
• /
• pp.29-40
• /
• 2007

An implicit discontinuous Galerkin method for the two-dimensional Euler equations was developed on unstructured triangular meshes. The method can achieve high-order spatial accuracy by using hierachical basis functions based on Legendre polynomials. Numerical tests were conducted to estimate the convergence order of numerical solutions to the Ringleb flow and the supersonic vortex flow for which analytic solutions are available. Also, the flows around a 2-D circular cylinder and an NACA0012 airfoil were numerically simulated. The numerical results showed that the implicit discontinuous Galerkin methods couples with a high-order representation of curved solid boundaries can be an efficient method to obtain very accurate numerical solutions on unstructured meshes.

• Honam Mathematical Journal
• /
• 제43권3호
• /
• pp.373-383
• /
• 2021

In this study, a numerical method deals with the Chebyshev wavelet collocation and Adomian decomposition methods are proposed for solving Korteweg-de Vries equation. Integration of the Chebyshev wavelets operational matrices is derived. This problem is reduced to a system of non-linear algebraic equations by using their operational matrix. Thus, it becomes easier to solve KdV problem. The error estimation for the Chebyshev wavelet collocation method and ADM is investigated. The proposed method's validity and accuracy are demonstrated by numerical results. When the exact and approximate solutions are compared, for non-linear or linear partial differential equations, the Chebyshev wavelet collocation method is shown to be acceptable, efficient and accurate.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• 제19권7호
• /
• pp.719-725
• /
• 2009

An efficient time-domain numerical method for the analysis of broadband noise generation and propagation due to turbulence-cascade interaction is developed. The core algorithm of the present method is based on the B-periodicity of the acoustic response function of the flat-airfoil cascade to the ingesting gust(B denotes the number of airfoils in the cascade). To confirm this periodicity, gust-cascade interaction problem are solved by using the time-domain method, which shows that the incident gust with the circumferential mode number having the same remainders when divided by the airfoil number excites the same acoustic response of the cascade. Using the proposed fast algorithm with this periodicity, we show that the total computation time for the model broadband problem using the total 525 incident gust modes can be reduced to about 1/4 of that taken in using the previous time-domain program.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 한국소음진동공학회 2009년도 춘계학술대회 논문집
• /
• pp.477-482
• /
• 2009

An efficient time-domain numerical method for the analysis of broadband noise generation and propagation due to turbulence-cascade interaction is developed. The core algorithm of the present method is based on the B-periodicity of the acoustic response function of the flat-airfoil cascade to the ingesting gust (B denotes the number of airfoils in the cascade). To confirm this periodicity, gust-cascade interaction problem are solved by using the time-domain method, which shows that the incident gust with the circumferential mode number having the same remainders when divided by the airfoil number excites the same acoustic response of the cascade. Using the proposed fast algorithm with this periodicity, we show that the total computation time for the model broadband problem using the total 525 incident gust modes can be reduced to about 1/4 of that taken in using the previous time-domain program.

• Proceedings of the Korean Operations and Management Science Society Conference
• /
• 한국경영과학회 1998년도 추계학술대회 논문집
• /
• pp.67-70
• /
• 1998

The computational speed of interior point method of linear programming depends on the speed of Cholesky factorization to solve AΘA$^{T}$ $\Delta$y=b. If the coefficient matrix A has dense columns then the matrix AΘA$^{T}$ becomes a dense matrix. This causes Cholesky factorization to be slow. The Schur complement method is applied to treat dense columns in many implementations but suffers from its numerical unstability. We study efficient implementation of Schur complement method. We achieve improvements in computational speed and numerical stability.rical stability.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제20권4호
• /
• pp.277-293
• /
• 2016

We develop an efficient numerical method for pricing the Derivative Linked Securities (DLS). The payoff structure of the hybrid DLS consists with a standard 2-Star step-down type ELS and the range accrual product which depends on the number of days in the coupon period that the index stay within the pre-determined range. We assume that the 2-dimensional Geometric Brownian Motion (GBM) as the model of two equities and a no-arbitrage interest model (One-factor Hull and White interest rate model) as a model for the interest rate. In this study, we employ the Monte Carlo simulation method with the Compute Unified Device Architecture (CUDA) parallel computing as the General Purpose computing on Graphic Processing Unit (GPGPU) technology for fast and efficient numerical valuation of DLS. Comparing the Monte Carlo method with single CPU computation or MPI implementation, the result of Monte Carlo simulation with CUDA parallel computing produces higher performance.

• Journal of applied mathematics & informatics
• /
• 제40권1_2호
• /
• pp.267-281
• /
• 2022

Herein, an algorithm for efficient evaluation of oscillatory Fourier-integrals with Jacobi-Cauchy type singularities is suggested. This method is based on the use of the traditional Clenshaw-Curtis (CC) algorithms in which the given function is approximated by the truncated Chebyshev series, term by term, and the oscillatory factor is approximated by using Bessel function of the first kind. Subsequently, the modified moments are computed efficiently using the numerical steepest descent method or special functions. Furthermore, Algorithm and programming code in MATHEMATICA^® 9.0 are provided for the implementation of the method for automatic computation on a computer. Finally, selected numerical examples are given in support of our theoretical analysis.

• Journal of computational fluids engineering
• /
• 제19권3호
• /
• pp.52-61
• /
• 2014

The present paper deals with the efficient computation of higher-order CFD methods for compressible flow using graphics processing units (GPU). The higher-order CFD methods, such as discontinuous Galerkin (DG) methods and correction procedure via reconstruction (CPR) methods, can realize arbitrary higher-order accuracy with compact stencil on unstructured mesh. However, they require much more computational costs compared to the widely used finite volume methods (FVM). Graphics processing unit, consisting of hundreds or thousands small cores, is apt to massive parallel computations of compressible flow based on the higher-order CFD methods and can reduce computational time greatly. Higher-order multi-dimensional limiting process (MLP) is applied for the robust control of numerical oscillations around shock discontinuity and implemented efficiently on GPU. The program is written and optimized in CUDA library offered from NVIDIA. The whole algorithms are implemented to guarantee accurate and efficient computations for parallel programming on shared-memory model of GPU. The extensive numerical experiments validates that the GPU successfully accelerates computing compressible flow using higher-order method.

• Proceedings of the Korea Concrete Institute Conference
• /
• 한국콘크리트학회 2006년도 춘계학술발표회 논문집(I)
• /
• pp.462-465
• /
• 2006

This paper is aimed at proposing the sampling method to reduce variance of statistical parameters in uncertainty analysis of concrete structures. The proposed method is a modification of Latin Hypercube sampling method. This uses specially modified tables of random permutations of rank number. Also, the Spearman coefficient is used to make modified tables. Numerical analysis is carried out to predict the uncertainty of axial shortening in prestressed concrete bridge. The numerical results show that the method is efficient for uncertainty analysis of complex structural system such as prestressed concrete bridges.