• 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.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.7
• /
• pp.1-8
• /
• 2003

Numerical characteristics of air chemistry associated with hypersonic flows are described and are compared with those of hydrogen oxygen combustion, applying the partially implicit time integration method to air chemistry. This paper reveals that the time integration of air chemistry needs a chemical Jacobian for stable calculations. However the positive real eigenvalues in air chemistry are relatively smaller than those of hydrogen combustion, and the numerical integration is less sensitive than that with combustion. lt is also found that the application of the partia1ly irnplicit method reduces the computing time without numerical instabilities.

• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.76-84
• /
• 2010

In order to analyze the periodic unsteady flow problem efficiently the partially implicit harmonic balance (PIHB) method was developed. Contrary to the existing harmonic balance method, this method handles the harmonic source term explicitly and deals with flux terms implicitly. This method has a good convergence in comparison with the full explicit harmonic method and it is easy to apply this method because there is no need to calculate the complicated flux Jacobian term by comparing with the full implicit harmonic method. With the multigrid method about the each harmonic it turns out that this method has a good convergence regardless of the number of harmonics. The oscillating flows over NACA0012 airfoil is considered to verify this method then the result correponsed to both the result of dual time stepping and explicit Runge-Kutta method.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.52-53
• /
• 2003

• Journal of Ocean Engineering and Technology
• /
• v.24 no.5
• /
• pp.31-38
• /
• 2010

The Floating storage and re-gasification unit (FSRU), which has large cargo storage tanks, is a floating liquefied natural gas (LNG) import terminal. The sloshing motion in tanks that are partially filled with LNG can cause impact pressure on the containment system and affect the global motion of the FSRU. Therefore, the accurate prediction of sloshing motion has been a significant issue in the offshore gas production industry. In this paper, a particle method based on the moving particle semi-implicit (MPS) method proposed by Koshizuka and Oka (1996) has been modified to predict sloshing motion accurately in a rectangular tank with the filling ratio of water. The simulation results, including the violent sloshing of the fluid, were validated by comparison with the original MPS method.

• Journal of Aerospace System Engineering
• /
• v.14 no.3
• /
• pp.17-23
• /
• 2020

In this paper, the efficient periodic unsteady flow analysis is developed by using a Chebyshev collocation operator applied to the time differential term of the governing equations. The partial implicit time integration method was also applied in the governing equation for a fluid, which means flux terms were implicitly processed for a time integration and the time derivative terms were applied explicitly in the form of the source term by applying the Chebyshev collocation operator. To verify this method, we applied the 1D unsteady Burgers equation and the 2D oscillating airfoil. The results were compared with the existing unsteady flow frequency analysis technique, the Harmonic Balance Method, and the experimental data. The Chebyshev collocation operator can manage time derivatives for periodic and non-periodic problems, so it can be applied to non-periodic problems later.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.4
• /
• pp.609-620
• /
• 1992

A numerical solution of temperature and thermally induced stress in a wafer during rapid thermal processing (R.T.P) is obtained, and an analysis of onset and propagation of slip is performed and compared with experiment. In order to calculate temperature distribution of a wafer in R.T.P system, heat conduction equation that incorporated with radiative and convective heat transfer model is solved, and the solution of the equation is calculated numerically using alternating direction implicit (A.D.I) method. In dealing with radiative heat transfer, a partially transparent body that absorbs the radiation energy is assumed and this transparent body undergoes multiple internal reflections and absorptions. Two dimensional (assuming plane stress) thermoelastic constitutive equation is used to calculate thermal stress induced in a wafer and finite element method is employed to solve the equation numerically. The stress resolved in the slip directions on the slip planes of silicon is compared with the yield stress of silicon in order to predict the slip. The result of the analysis shows that the wafer temperature at which slip occurs is affected by the heating rate of the R.T.P system. It is observed that once slip occurs in the wafer, the slip grows.