• Journal of Advanced Marine Engineering and Technology
• /
• v.22 no.5
• /
• pp.670-679
• /
• 1998

The solution of hyperbolic equation with wave nature has sharp discontinuties in the medium at the wave front. Difficulties encounted in the numrtical solution of such problem in clude among oth-ers numerical oscillation and the representation of sharp discontinuities with good resolution at the wave front. In this work inviscid Burgers equation and modified heat conduction equation is intro-duced as hyperboic equation. These equations are caculated by numerical methods(explicit method MacCormack method Total Variation Diminishing(TVD) method) along various Courant numbers and numerical solutions are compared with the exact analytic solution. For inviscid Burgers equa-tion TVD method remains stable and produces high resolution at sharp wave front but for modified heat Conduction equation MacCormack method is recommmanded as numerical technique.

• Journal of Korea Game Society
• /
• v.18 no.5
• /
• pp.123-132
• /
• 2018

Various methods have been proposed to efficiently animate the motion of soft objects in realtime. In order to maintain the topology between the elements of the objects, it is required to employ constraint forces, which limit the size of the time steps for the numerical integration and reduce the efficiency. To tackle this, an implicit method with larger steps was proposed. However, the method is, in essence, a linear system with a large matrix, of which solution requires heavy computations. Several approximate methods have been proposed, but the approximation is obtained with an increased damping and the loss of accuracy. In this paper, new integration method based on harmonic oscillation with better stability was proposed, and it was further stabilized with the hybridization with approximate implicit method. GPU parallelism can be easily implemented for the method, and large-scale soft objects can be simulated in realtime.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.516-523
• /
• 2009

The artificial compressibility (AC) method for the incompressible Navier-Stokes equations in the generalized curvilinear coordinates using the primitive form is implemented. The main advantage of the AC approach is that the resulting system of equations resembles the system of compressible N-S equations and can thus be integrated in time using standard, well-established time-marching methods. The errors, which are the odd-even oscillation, for pressure field in using the artificial compressibility can be eliminated by using the $4^{th}$ order artificial dissipation term which is explicitly included. Even though this paper focuses exclusively on 2D laminar flows to validate and assess the performance of this solver, this numerical method is general enough so that it can be readily extended to carry out 3D URANS simulation of engineering flows. This algorithm yields practically identical velocity profiles and primary vortex and secondary vortices that are in excellent overall agreement with the results of the vorticity-stream function formulation (Ghia et al., 1982). However, the grid resolution have to be required to be large enough to express the various vortices.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.11
• /
• pp.948-956
• /
• 2012

This paper presents the computations of the longitudinal and directional stability derivatives for the SDM(Standard Dynamic Model). The static and dynamic derivatives are evaluated at once using forced harmonic oscillations in the pitch and yaw directions. For the numerical simulations, a 3-D Euler solver that uses a dual time stepping method for unsteady time accurate simulations is applied. This work investigates the variation of the derivatives in terms of the Mach number and the several motion parameters. Good agreement of the pitch and yaw stability derivatives with previously published numerical results and experimental results are observed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.3
• /
• pp.267-279
• /
• 1998

A numerical study was made of flow behind a circular cylinder in a uniform flow, where the cylinder was rotationally oscillated in time. The temporal behavior of vortex formation was scrutinized over broad ranges of the two externally specified parameters, i.e., the dimensionless rotary oscillating frequency (.110.leq. $S_{f}$.leq..220) and the maximum angular amplitude of rotation (.theta.$_{max}$=15 deg., 30 deg. and 60 deg.). The Reynolds number (Re= $U_{{\inf}D}$.nu.) was fixed at Re=110. A fractional-step method was utilized to solve the Navier-Stokes equations with a generalized coordinate system. The main emphasis was placed on the initial vortex formations by varying $S_{f}$ and .theta.$_{max}$. Instantaneous streamlines and pressure distributions were displayed to show the vortex formation patterns. The vortex formation modes and relevant phase changes were characterized by measuring the lift coefficient ( $C_{L}$) and the time of negative maximum $C_{L}$( $t_{-C}$$_{Lmax}$) with variable forcing conditions.s.tions.s.s.s.

• The KSFM Journal of Fluid Machinery
• /
• v.11 no.3
• /
• pp.14-21
• /
• 2008

A numerical simulation on the flow field of a valveless bidirectional piezoelectric micropump has been performed. In this type of micropump, the oscillation of the piezoelectric diaphragm generates the blowing and suction flow through the oblique channel from the pumping chamber. The angle between the oblique and main channel causes the variation of flow distribution through upstream and downstream channels in suction and blowing modes. In the suction flow mode, the working fluid flows from both the upstream and downstream of the main channel to the pumping chamber through the oblique channel. However, in the blowing flow mode, the fluid pushed out of the pumping chamber flows more toward the downstream of the main channel due to the inertia of the fluid. In the present study, the effects of geometries such as the angle of oblique channel and the shape of main channel on the flow rate of the up/downstream were investigated. The flow rate obtained from the pump and the energy required to the pump were also analyzed for various displacements and frequencies of the oscillation of the diaphragm.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.4
• /
• pp.439-451
• /
• 1994

A numerical model for harbour oscillation is presented by use of Galerkin finite element method. The governing equation is used by the modified mild slope equation derived from Chen (1986) in which bottom friction is incorporated. Since the existing absorbing boundary condition. however. is shown to be incorrect correct boundary condition and forcing term due to an incident plane wave are rederived. Computation results for a rectangular harbour are shown in comparison with both laboratory data and existing numerical results. After the values of friction factor (f) and reflection coefficient (K$_{r}$) are discussed, the set (K$_{r}$=0, 94, f=0) is found to be best fitted to the laboratory data of the rectangular harbour.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.2
• /
• pp.139-148
• /
• 2012

The MLS(Moving Least Squares) Difference Method is a numerical scheme that combines the MLS method of Meshfree method and Taylor expansion involving not numerical quadrature or mesh structure but only nodes. This paper presents an dynamic algorithm of MLS difference method for solving transient solid mechanics problems. The developed algorithm performs time integration by using Newmark method and directly discretizes strong forms. It is very convenient to increase the order of Taylor polynomial because derivative approximations are obtained by the Taylor series expanded by MLS method without real differentiation. The accuracy and efficiency of the dynamic algorithm are verified through numerical experiments. Numerical results converge very well to the closed-form solutions and show less oscillation and periodic error than FEM(Finite Element Method).

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.2
• /
• pp.21-29
• /
• 2002

The M-AUSMPW+ scheme that can capture shock waves exactly with monotonic characteristic is developed by analyzing the cause of oscillation in shock regions. Firstly shock-capturing characteristics of general FVS including the AUSM-type schemes are investigated in detail, according to the different between a cell-interface and a sonic transition position. The cause of oscillation is the improper numerical dissipation that could not represent the real physics. The M-AUSMPW+ could capture shocks exactly without oscillatory behaviors in considering the sonic transition position and an cell-interface position

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.7 no.1
• /
• pp.108-116
• /
• 1997

A laboratory experiment was made of critical transition flow modes in Czocllralski convection. Numerical computation was also made to delineate the dynamic transition. The period of temperature oscillation ($t_p$) and the interval of temperature oscillation ($\Delta\theta$) were scrutinized to capture the critical transition regime. The mixed convection parameter was varied in the range of $0.134\le Ra/PrRe^2 \le3.804$. The data from calculation were in good agreement with ones from experiment. The influence of the Prandtl number on the transition was examined for Pr ＝ 910, 4445 and 8889. To understand the transition mechanism, the detailed temperature oscillation modes, the isolines of meridional temperature and the axial velocity profiles were investigated.