• Journal of Astronomy and Space Sciences
• /
• v.3 no.2
• /
• pp.117-122
• /
• 1986

The broadband electrostatic noise has been observed in the boundary layer region of the earth's magnetosphers. These electrostatic waves believed to be generated by drifting ion beams in the matnetotail. We have shown the numerical simulation result of ion beam acoustic instability in the magnetotail. This instability heats both background and beam ion in the boundary layer of neutral sheet observed by satellite.

• Structural Engineering and Mechanics
• /
• v.71 no.1
• /
• pp.65-75
• /
• 2019

In this article, the dynamic fracture instability characteristics, including dynamic crack propagation and crack branching, in PMMA brittle solids under dynamic loading are investigated using the discrete element method (DEM) simulations. The microscopic parameters in DEM are first calibrated using the comparison with the previous experimental results not only in the field of qualitative analysis, but also in the field of quantitative analysis. The calibrating process illustrates that the selected microscopic parameters in DEM are suitable to effectively and accurately simulate dynamic fracture process in PMMA brittle solids subjected to dynamic loads. The typical dynamic fracture behaviors of solids under dynamic loading are then reproduced by DEM. Compared with the previous experimental and numerical results, the present numerical results are in good agreement with the existing ones not only in the field of qualitative analysis, but also in the field of quantitative analysis. Furthermore, effects of dynamic loading magnitude, offset distance of the initial crack and initial crack length on dynamic fracture behaviors are numerically discussed.

• Journal of the Korean Mathematical Society
• /
• v.57 no.1
• /
• pp.249-281
• /
• 2020

A reaction-diffusion model with spatiotemporal delay modeling the dynamical behavior of a single species is investigated. The parameter regions for the local stability, global stability and instability of the unique positive constant steady state solution are derived. The conditions of the occurrence of Turing (diffusion-driven) instability are obtained. The existence of time-periodic solutions, the existence and nonexistence of nonconstant positive steady state solutions are proved by bifurcation method and energy method. Numerical simulations are presented to verify and illustrate the theoretical results.

• 한국연소학회:학술대회논문집
• /
• 2015.12a
• /
• pp.115-116
• /
• 2015

In the present approach, the flame-acoustics interaction is represented by FDF (Flame Describing Function) which is a important source term in the Helmholtz' equation. In this study, the combustion instability is analyzed by the forced mode strategy with the measured FDF. Numerical results indicate that the present approach reasonably well predicts the essential features of the combustion instability characteristics in the lean premixed combustor under the gas-turbine like environment.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.2
• /
• pp.139-155
• /
• 2019

In this paper, we study the asymptotic behavior and Hopf bifurcation of the modified Holling-Tanner models for the predator-prey interactions in the absence of diffusion. Further the direction of Hopf bifurcation and stability of bifurcating periodic solutions are investigated. Diffusion driven instability of the positive equilibrium solutions and Turing instability region regarding the parameters are established. Finally we illustrate the theoretical results with some numerical examples.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.281-283
• /
• 2001

As a companion to an adiabatic version developed by Ryu and his coworkers, we have built an isothermal magnetohydrodynamic code for astrophysical flows. It is suited for the dynamical simulations of flows where cooling timescale is much shorter than dynamical timescale, as well as for turbulence and dynamo simulations in which detailed energetics are unimportant. Since a simple isothermal equation of state substitutes the energy conservation equation, the numerical schemes for isothermal flows are simpler (no contact discontinuity) than those for adiabatic flows and the resulting code is faster. Tests for shock tubes and Alfven wave decay have shown that our isothermal code has not only a good shock capturing ability, but also numerical dissipation smaller than its adiabatic analogue. As a real astrophysical application of the code, we have simulated the nonlinear three-dimensional evolution of the Parker instability. A factor of two enhancement in vertical column density has been achieved at most, and the main structures formed are sheet-like and aligned with the mean field direction. We conclude that the Parker instability alone is not a viable formation mechanism of the giant molecular clouds.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.8
• /
• pp.702-710
• /
• 2007

Direct numerical simulation (DNS) is performed to investigate suppressed wake transition and occurrence of lock-on in the wake of a circular cylinder disturbed by sinusoidal perturbation at the Reynolds number of 220 (A-mode instability regime). The sinusoidal perturbation, of which the frequency is near twice the natural shedding frequency, is superimposed on the free stream velocity. It is shown that the wake transition behind the circular cylinder can be suppressed due to the perturbation of the free stream velocity. This change causes a jump in the Strouhal number from the value corresponding to A-mode instability regime to the value corresponding to retarded wake transition regime (extrapolated from laminar shedding regime) in the Strouhal-Reynolds number relationship. As a result, vortex shedding frequency is locked on the perturbation frequency depending not on the natural shedding frequency but on the modified shedding frequency.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.84-90
• /
• 2001

Secondary instability of flow past a circular cylinder is examined using direct numerical simulation at Reynolds number 220 and 250. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme, and the pressure and viscosity terms are applied by the 4th order central difference. In spanwise, Navier-Stokes equation is distributed using Spectral Method. The critical Reynolds number for this instability is found to be about Re=190. The secondary instability leads re three-dimensionality with a spanwise wavelength about 4 cylinder diameters at onset (A-mode). Results of three-dimensional effect in wake of a circular cylinder are represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.4
• /
• pp.415-423
• /
• 2003

The hydrodynamic instability of the three-dimensional boundary-layer over a rotating disk has been numerically investigated for these flows; Ro = -1, -0.5, and 0, using linear stability theory. Detailed numerical values of the disturbance wave number. wave frequency. azimuth angle. radius (Reynolds number, Re) and other characteristics have been calculated for the pre-swirl flows. On the basis of Ekman and Karman boundary layer theory, the instability of the pre-swirl flows have been investigated for the unstable criteria. The disturbance will be relatively fast amplified at small fe and within wide bands of wave number compared with previously known Karman boundary-layer results. The flow (Ro =-0.5) is found to be always stable for a disturbance whose dimensionless wave number is greater than 0.9. It has a larger range of unstable interval than Karman boundary layer and can be unstable at smaller Re.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.7
• /
• pp.695-702
• /
• 2008

Flow instability in the rocket turbo pump systems can be caused by various elements such as valve, orifice and venturi and etc. The formation of cavitation specially in the propellant feeding system can trigger the mass flow and pressure oscillation due to cyclic formation and depletion of cavitations. If the cryogenic propellant are used, which is very sensitive to temperature variation, the change of propellant properties due to thermodynamic effect should be accounted for in the flow analysis. This study focuses on the formation of cryogenic cavitation adopting MUSHY IDM model suggested by Shyy and coworkers. Also, the flow instability is investigated with developed numerical code in the downstream of orifice flow. To this end, three different orifices are selected and investigated by the numerical calculation.