• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.12
• /
• pp.1153-1161
• /
• 2010

An efficient solution method for harmonic balance techniques with Fourier transform is presented for periodic unsteady flow problems. The present partially-implicit harmonic balance treats the flux terms implicitly and the harmonic source term is solved explicitly. The convergence of the partially Implicit method is much faster than the explicit Runge-Kutta harmonic balance method. The method does not need to compute the additional flux Jacobian matrix from the implicit harmonic source term. Compared with fully implicit harmonic balance method, this partial approach turns out to have good convergence property. Oscillating flows over NACA0012 airfoil are considered to verify the method and to compare with results of explicit R-K(Runge-Kutta) and dual time stepping methods.

• Communications of the Korean Mathematical Society
• /
• v.14 no.4
• /
• pp.833-842
• /
• 1999

In this paper, we construct a fully discrete solution of the incompressible Navier Stokes equations using implicit Runge kutta method. We prove the existence of the fully discrete solution.

• Bulletin of the Korean Mathematical Society
• /
• v.29 no.1
• /
• pp.31-40
• /
• 1992

We introduce the generalized Runge-Kutta methods with the exponentially dominant order .omega. in [3], and the convergence theorems of the generalized explicit Euler method are derived in [4]. In this paper we will study the convergence of the generalized implicit Euler method.

• Journal of applied mathematics & informatics
• /
• v.15 no.1_2
• /
• pp.1-28
• /
• 2004

We first apply a first order splitting to a semilinear reaction-diffusion equation and then discretize the resulting system by an $H^1$-Galerkin mixed finite element method in space. This semidiscrete method yields a system of differential algebraic equations (DAEs) of index one. A priori error estimates for semidiscrete scheme are derived for both differ-ential as well as algebraic components. For fully discretization, an implicit Runge-Kutta (IRK) methods is applied to the temporal direction and the error estimates are discussed for both components. Finally, we conclude the paper with a numerical example.

• Journal of applied mathematics & informatics
• /
• v.4 no.1
• /
• pp.17-30
• /
• 1997

In this work we approximate the solution of initialboun-dary value problem using a Galerkin-finite element method for the spatial discretization and Implicit Runge-Kutta method for the spatial discretization and implicit Runge-Kutta methods for the time stepping. To deal with the nonlinear term f(x, t, u), we introduce the well-known extrapolation sheme which was used widely to prove the convergence in $L_2$-norm. We present computational results showing that the optimal order of convergence arising under $L_2$-norm will be preserved in $L_{\infty}$-norm.

• Tribology and Lubricants
• /
• v.33 no.6
• /
• pp.309-314
• /
• 2017

This paper presents a theoretical investigation of the dynamic characteristics of externally pressurized air pad bearings with closed loop grooves. These grooves are made on the surface of bearings to reduce the number of supply holes so that manufacturing costs can be reduced. The semi-implicit method is applied to calculate the time varying pressure profile on the air bearing surface owing to the advantages of numerical stability and fast time tracing characteristics. The static pressure of the groove bearings is much higher than that without grooves, so the groove bearings can provide high load carrying capacity. The equation of motion considering vertical motion and tilting motion are also solved using the Runge-Kutta 4th order method. By combining the semi-implicit method and the Runge-Kutta method, fast calculations of the dynamic behavior of the air bearing can be achieved. The variations of bearing reaction force, air film reaction moment, height, and tilting angle are investigated for the step force input, which is 20% higher than the bearing reaction, when the nominal clearance is 6 mm. The effect of the groove width and the groove depth are investigated by calculating the dynamic behavior. The possibility of the air hammering with the depth of the groove is found and discussed.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.8-12
• /
• 1997

The three-dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for three contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreement.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.113-117
• /
• 1995

Viscous analysis on incompressible flows is performed using unstructured triangular meshes. A two-dimensional and axisymmetric incompressible Navier-Stokes equations are solved in time-marching form by artificial compressibility method. The governing equations are discretized by a cell-centered based finite-volume method. and a centered scheme is used for inviscid and viscous fluxes with fourth order artificial dissipation. An explicit multi-stage Runge-Kutta method is used for the time integration with local time stepping and implicit residual smoothing. Convergence properties are examined and solution accuracies are also validated with benchmark solution and experiment.

• Korean Journal of Mathematics
• /
• v.13 no.2
• /
• pp.161-176
• /
• 2005

We consider a model of population dynamics whose mortality function is unbounded. We note that the regularity of the solution depends on the growth rate of the mortality near the maximum age. We propose Gauss-Legendre methods along the characteristics to approximate the solution when the solution is smooth enough. It is proven that the scheme is convergent at fourth-order rate in the maximum norm. We also propose discontinuous Galerkin finite element methods to approximate the solution which is not smooth enough. The stability of the method is discussed. Several numerical examples are presented.

• Journal of the Korean Society of Combustion
• /
• v.13 no.1
• /
• pp.17-22
• /
• 2008

A partly implicit/quasi-explicit method is introduced for the solution of detailed chemical kinetics with stiff source terms based on the standard fourth-order Runge-Kutta scheme. Present method solves implicitly only the stiff reaction rate equations, whereas the others explicitly. The stiff equations are selected based on the survey of the chemical Jaconian matrix and its Eigenvalues. As an application of the present method constant pressure combustion was analyzed by a detailed mechanism of hydrogen-air combustion with NOx chemistry. The sensitivity analysis reveals that only the 4 species in NOx chemistry has strong stiffness and should be solved implicitly among the 13 species. The implicit solution of the 4 species successfully predicts the entire process with same accuracy and efficiency at half the price.