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

• 한국연소학회:학술대회논문집
• /
• 2012.11a
• /
• pp.61-62
• /
• 2012

The present study has numerically investigated the effects of fuel-side dilution and pressure on flame structure and extinction scalar dissipation rate of turbulent syngas nonpremixedd flames. Numerical results indicate that for highly diluted case, peak temperature is decreased and stoichiometric mixture fraction is increased. By decreasing the pressure and the nitrgen dilution levelcreased, the extinction scalar dissipation rate is increased.

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

• Journal of Mechanical Science and Technology
• /
• v.19 no.7
• /
• pp.1517-1525
• /
• 2005

Liquid sloshing displays the highly nonlinear free surface fluctuation when either the external excitation is of large amplitude or its frequency approaches natural sloshing frequencies. Naturally, the accurate tracking of time-varying free surface configuration becomes a key task for the reliable prediction of the sloshing time-history response. However, the numerical instability and dissipation may occur in the nonlinear sloshing analysis, particularly in the long-time beating simulation, when two simulation parameters, the relative time-increment parameter a and the fluid mesh pattern, are not elaborately chosen. This paper intends to examine the effects of these two parameters on the potential-based nonlinear finite element method introduced for the large amplitude sloshing flow.

• Advances in aircraft and spacecraft science
• /
• v.9 no.5
• /
• pp.463-477
• /
• 2022

This study compares various ways of calculating flows for the problems with the presence of shock waves by first-order schemes and higher-order DG method on the tests from the Quirk list, namely: Quirk's problem and its modifications, shock wave diffraction at a 90 degree corner, the problem of double Mach reflection. It is shown that the use of HLLC and Godunov's numerical schemes flows in calculations can lead to instability, the Rusanov-Lax-Friedrichs scheme flow can lead to high dissipation of the solution. The most universal in heavy production calculations are hybrid schemes flows, which allow the suppression of the development of instability and conserve the accuracy of the method.

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

• Geomechanics and Engineering
• /
• v.8 no.4
• /
• pp.495-510
• /
• 2015

The applications of water jet cutting (WJC) in coal mine have progressed slowly. In this paper, we analyzed the possibility and reasonableness of WJC application to pressure relief in hard coal seam, simulated the distributive characteristics of stress and energy fields suffered by hard coal roadway wallrock and the internal relationships of the fields to the instability due to WJC (including horizontal radial slot and vertical annular slot) on roadway wallrock. The results showed that: (1) WJC can unload hard coal seam effectively by inducing stress release and energy dissipation in coal mass near its slots; its annular slots also can block or weaken stress and energy transfer in coal mass; (2) the two slots may cause "the beam structure" and "the small pillar skeleton", and "the layered energy reservoir structure", respectively, which lead to the increase in stress concentration and energy accumulation in coal element mass near the slots; (3) the reasonable design and optimization of slots' positions and their combination not only can significantly reduce the scope of stress concentration and energy accumulation, but also destroy coal mass structure on a larger scale to force stress to transfer deeper coal mass.

• Earthquakes and Structures
• /
• v.13 no.3
• /
• pp.279-288
• /
• 2017

Having the ability to keep on yielding stable solutions in problems involving high potential of instability, composite time integration methods have become very popular among scientists. These methods try to split a time step into multiple sub-steps so that each sub-step can be solved using different time integration methods with different behaviors. This paper proposes a new composite time integration in which a time step is divided into two sub-steps; the first sub-step is solved using the well-known Newmark method and the second sub-step is solved using Simpson's Rule of integration. An unconditional stability region is determined for the constant parameters to be chosen from. Also accuracy analysis is perform on the proposed method and proved that minor period elongation as well as a reasonable amount of numerical dissipation is produced in the responses obtained by the proposed method. Finally, in order to provide a practical assessment of the method, several benchmark problems are solved using the proposed method.

• The Journal of the Acoustical Society of Korea
• /
• v.25 no.1E
• /
• pp.32-35
• /
• 2006

The screech tone of underexpanded jet is numerically calculated without any specific modeling for the screech tone itself. Fourth-order optimized compact scheme and fourth-order Runge-Kutta method are used to solve the 2D axisymmetric Euler equation. Adaptive nonlinear artificial dissipation model and generalized characteristic boundary condition are also used. The screech tone, generated by a closed loop between instability waves and quasi-periodic shock cells at the near field, is reasonably analyzed with present numerical methods for the underexpanded jet having Mach number 1.13. First of all, the centerline mean pressure distribution is calculated and compared with experimental and other numerical results. The instantaneous density contour plot shows Mach waves due to mixing layer convecting supersonically, which propagate downstream. The pressure signal and its Fourier transform at upstream and downstream shows the directivity pattern of screech tone very clearly. Most of all, we can simulate the axisymmetric mode change of screech tone very precisely with present method. It can be concluded that the basic phenomenon of screech tone including the frequency can be calculated by using high-order and high-resolution schemes without any specific numerical modeling for screech tone feedback loop.