• Structural Engineering and Mechanics
• /
• v.65 no.1
• /
• pp.91-98
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.

• Structural Engineering and Mechanics
• /
• v.65 no.1
• /
• pp.121-128
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.

• Journal of Mechanical Science and Technology
• /
• v.19 no.9
• /
• pp.1811-1820
• /
• 2005

Thermo-acoustic instability of methane/ air flames in an industrial gas-turbine combustor is numerically investigated adopting CFD analysis. The combustor has 37 EV burners through which methane and air are mixed and then injected into the chamber. First, steady fuel! air mixing and flow characteristics established by the burner are investigated by numerical analysis with single burner. And then, based on information on the flow data, the burners are modeled numerically via equivalent swirlers, which facilitates the numerical analysis with the whole combustion system including the chamber and numerous burners. Finally, reactive flow fields within the chamber are investigated numerically by unsteady analysis and thereby, spontaneous instability is simulated. Based on the numerical results, scaling analysis is conducted to find out the instability mechanism in the combustor and the passive control method to suppress the instability is proposed and verified numerically.

• The Bulletin of The Korean Astronomical Society
• /
• v.39 no.1
• /
• pp.37.2-37.2
• /
• 2014

Gas in disk galaxies interacts nonlinearly with a underlying stellar spiral potential to form galactic spiral shocks. Numerical simulations typically show that these shocks are unstable to the wiggle instability, forming non-axisymmetric structures with high vorticity. While previous studies suggested that the wiggle instability may arise from the Kelvin-Helmholtz instability or orbit crowding of gas elements near the shock, its physical nature remains uncertain. It was even argued that the wiggle instability is of numerical origin, caused by the inability of a numerical code to resolve a shock that is inclined to numerical grids. In this work, we perform a normal-mode linear stability analysis of galactic spiral shocks as a boundary-value problem. We find that the wiggle instability originates physically from the potential vorticity generation at a distorted shock front. As the gas follows galaxy rotation, it periodically passes through multiple shocks, successively increasing its potential vorticity. This sets up a normal-mode that grows exponentially, with a growth rate comparable to the orbital angular frequency. We show that the results of our linear stability analysis are in good agreement with the those of local hydrodynamic simulations of the wiggle instability.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.9
• /
• pp.659-668
• /
• 2008

Combustion instability has been considered as very important issue for developing gas turbine and rocket engine. There is a need for fundamental understanding of combustion instability. In this study, combustion instability was numerically and experimentally investigated in a dump combustor with bluff body. The fuel and air mixture had overall equivalence ratio of 0.9 and was injected toward dump combustor. The pressure oscillation with approximately 256Hz was experimentally obtained. For numerical simulation, the standard k-$\varepsilon$ model was used for turbulence and the hybrid combustion model (eddy dissipation model and kinetically controlled model) was applied. After calculating steady solution, unsteady calculation was performed with forcing small perturbation on initial that solution. Pressure amplitude and frequency measured by pressure sensor is nearly the same as those predicted by numerical simulation. Furthermore, it is clear that a combustion instability involving vortex shedding is affected by acoustic-vortex-combustion interaction. The phase difference between the pressure and velocity is $\pi$/2, and that between the pressure and heat release rate is in excitation range described by Rayleigh, which is obvious that combustion instability for the bluff body combustor meets thermoacoustic instability criterion.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.141-150
• /
• 2002

Numerical simulations based on the ALE finite element method are carried out to examine the aerodynamics of an oscillating circular cylinder when the separated shear flows around the cylinder are stimulated by periodic jet excitation with a shear layer instability frequency. The excitation is applied to the flows from two points on the cylinder surface. The numerical results showed that the excitation with a shear layer instability frequency can reduce the negative damping and thereby stabilize the aerodynamics of the oscillating cylinder. The change of the lift phase seems important in stabilizing the cylinder aerodynamics. The change of lift phase is caused by the merger of the vortices induced by the periodic excitation with a shear layer instability frequency, and the vortex merging comes from the high growth rate, the rapid increase of wave number and decrease of phase velocity for the periodic excitation in the separated shear flows.

• Coupled systems mechanics
• /
• v.10 no.1
• /
• pp.79-102
• /
• 2021

In this paper we deal with instability problems of structures under nonconservative loading. It is shown that such class of problems should be analyzed in dynamics framework. Next to analytic solutions, provided for several simple problems, we show how to obtain the numerical solutions to more complex problems in efficient manner by using the finite element method. In particular, the numerical solution is obtained by using a modified Euler-Bernoulli beam finite element that includes the von Karman (virtual) strain in order to capture linearized instabilities (or Euler buckling). We next generalize the numerical solution to instability problems that include shear deformation by using the Timoshenko beam finite element. The proposed numerical beam models are validated against the corresponding analytic solutions.

• Journal of the Korean Society of Combustion
• /
• v.8 no.1
• /
• pp.46-56
• /
• 2003

The present paper is concerned with the development of the computational biology in the past half century and its relationship with combustion. The modem computational biology is considered to be initiated by the work of Alan Turing on the morphogenesis in 1952. This paper first touches the life and scientific achievement of Alan Turing and his theory on the morphogenesis based on the reactive-diffusive instability, called the Turing instability. The theory of Turing instability was later extended to the nonlinear realm of the reactive-diffusive systems, which is discussed in the framework of the excitable media by using the Oregonator model. Then, combustion analogies of the Turing instability and excitable media are discussed for the cellular instability, pattern forming combustion phenomena and flame edge. Finally, the recent efforts on numerical simulations of biological systems, employing the detailed bio-chemical knietic mechanism is discussed along with the possibility of applying the numerical combustion techniques to the computational cell biology.

• Structural Engineering and Mechanics
• /
• v.62 no.2
• /
• pp.151-162
• /
• 2017

Two $Chang-{\alpha}$ dissipative family methods and two $KR-{\alpha}$ family methods were developed for time integration recently. Although the four family methods are in the category of the dissipative structure-dependent integration methods, their performances may be drastically different due to the detrimental property of weak instability or overshoot for the two $KR-{\alpha}$ family methods. This weak instability or overshoot will result in an adverse overshooting behavior or even numerical instability. In general, the four family methods can possess very similar numerical properties, such as unconditional stability, second-order accuracy, explicit formulation and controllable numerical damping. However, the two $KR-{\alpha}$ family methods are found to possess a weak instability property or overshoot in the high frequency responses to any nonzero initial conditions and thus this property will hinder them from practical applications. Whereas, the two $Chang-{\alpha}$ dissipative family methods have no such an adverse property. As a result, the performances of the two $Chang-{\alpha}$ dissipative family methods are much better than for the two $KR-{\alpha}$ family methods. Analytical assessments of all the four family methods are conducted in this work and numerical examples are used to confirm the analytical predictions.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1650-1657
• /
• 2019

Fluidelastic instability and nonlinear dynamics of tube bundles is a key issue in a steam generator. Especially, once the post-instability motion of the tube becomes larger than the clearance gap to other tubes, effective contact or impact between the tubes under consideration and the other tube inevitable. There is seldom theoretical analysis to the nonlinear dynamic characteristics of a tube array in two-phase flow. In this paper, experimental and numerical studies were utilized to obtain the critical velocity of the flow-induced instability of a rotated triangular tube array. The calculation results agreed well with the experimental data. To explore the post-instability dynamics of the tube array system, a Runge-Kutta scheme was used to solve the nonlinear governing equations of tube motion. The numerical results indicated that, when the flow pitch velocity is larger than the critical velocity, the tube array system is undergoing a limit cycle motion, and the dynamic characteristics of the tube array are almost similar for different void fractions.