• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.355-360
• /
• 2001

This study addresses a computational work of the impulsive wave which is discharged from the open end of a pipe. An initial compression wave inside the pipe is assumed to propagate toward atmosphere. The over pressure and wave-length of the initial compression wave are changed to investigate the characteristic values of the impulsive wave. The second order total variation diminishing (TVD) scheme is employed to solve the axisymmetric, compressible, unsteady Euler equations. The relationship between the initial compression wave form and impulsive wave is characterized in terms of the peak pressure of the impulsive wave and its directivity. The results obtained show that for the initial compression wave of a large wave-length the peak pressure of the impulsive wave does not depend on the over pressure of the initial compression wave, but for the initial compression wave of a very short wave-length, like a shock wave, the peak pressure of the impulsive wave is increased with an increase in the over pressure of the initial compression wave. The directivity of the impulsive wave to the pipe axis becomes significant with a decrease in the wave-length of the initial compression wave.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.7 s.250
• /
• pp.671-681
• /
• 2006

This study focuses on the development of numerical procedure to analyze the nonlinear combustion instabilities in liquid rocket engine. Nonlinear behaviors of acoustic instabilities are characterized by the existence of limit cycle in linearly unstable engines and nonlinear or triggering instability in linearly stable engines. To discretize convective fluxes with high accuracy and robustness, approximated Riemann solver based on characteristics and Euler-characteristic boundary conditions are employed. The present procedure predicts well the transition processes from initial harmonic pressure disturbance to N-like steep-fronted shock wave in a resonant pipe. Longitudinal pressure oscillations within the SSME(Space Shuttle Main Engine) engine have been analyzed using the pressure-sensitive time lag model to account for unsteady combustion response. It is observed that the pressure oscillations reach a limit cycle which is independent of the characteristics of the initial disturbances and depends only on combustion parameters and operating conditions.

• Nuclear Engineering and Technology
• /
• v.31 no.3
• /
• pp.328-334
• /
• 1999

Based on the Euler beam theory and the elastic strain energy method, the elastic stiffness formula of the holddown spring assembly consisting of several leaves was previously derived. Even though the previous formula was known to be useful to estimate the elastic stiffness of the holddown spring assembly, recently it was reported that the elastic stiffness from the previous formula deviated greatly from the test results as the number of leaves was increased. The objective of this study is to extend the previous formula in order to resolve such an increasing deviation when increasing the number of leaves. Additionally, considering the friction forces acting on the interfaces between the leaves, we obtained an extended elastic stiffness formula. The characteristic test and the elastic stiffness analysis on the various kinds of specimens of the holddown spring assembly have been carried out; the validity of the extended formula has been verified by the comparison of their results. As a result of comparisons, it is found that the extended formula is able to evaluate the elastic stiffness of the holddown spring assembly within the maximum error range of + 12%, irrespective of the number of the leaves.

• Bulletin of the Korean Mathematical Society
• /
• v.31 no.1
• /
• pp.15-23
• /
• 1994

Throughout this paper, we are working on the complex number field C. The aim of this paper is to explain the applications of Theorem 2 in .cint. 1. In the surface theory, the adjoint linear system has played important roles and many tools have been developed to understand it. In the cases of higher dimensional varieties, we don't have any useful tools so far. Theorem 2 implies that it is enough to compute the dimension of the adjoint linear system to check the birationality. We can compute, somehow, the dimension of the adjoint linear system. For example, we can get an information about $h^{0}$ (X, $O_{x}$( $K_{x}$ + D)) from Euler characteristic of vertical bar $K_{X}$ + D vertical bar and some vanishing theorems. We are going to show the applications of Theorem 2 to smooth three-folds and smooth fourfold, specially, of general type with a nef canonical divisor, smooth Fano variety, and Calabi-Yau manifold. Our main results are Theorem A and Theorem B. Most of birationality problems in Theorem A and Theorem B have been studied. (see Ando [1] and Matsuki [4] for the detail matters.) But Theorem 2 gives short and easy proofs in the cases of dimension 3 and improves the previously known results in the cases of dimension 4.4. 4.4.

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

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.66-72
• /
• 1997

Optimum nozzle design exploiting the method of characteristic(M.O.C) has been in application as an efficient design methodology targeting a less weighted and short expansion nozzle. This paper treats the optimum nozzle design and the analysis of the inviscid compressible flow inside. Based on traditional Rao's method, the optimum nozzle design is coded with minor modifications for the identification of the control surface across which the mass flux should be conserved. Internal flow field is simulated numerically by M.O.C and implicit/explicit Taylor-Galerkin finite element method(F.E.M) with the aid of adaptive remeshing to capture the shock wave, hence improve the accuracy. Designed and calculated flow fields due to the separate analyses show that the mass flux predicted by optimum nozzle design with M.O.C is not conserved across the control surface and the sonic line should be located upstream of the nozzle throat. Rao's optimum nozzle design methodology exaggerates the momentum thrust and tends to overemphasize the engine performance loss.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.23 no.11 s.170
• /
• pp.1940-1951
• /
• 1999

In this paper, the equations of motion of a Bernoulli-Euler cantilever beam fixed on a moving cart and carrying a lumped mass concentrated at an arbitrary position along the beam is derived. The motion of the beam-mass-cart system is analyzed through unconstrained modal analysis, and a unified characteristic equation for calculating the natural frequencies of the system is obtained. The changes of natural frequencies and the corresponding mode shapes with respect to the changes in mass ratios of the system and to the concentrated position of the lumped mass are investigated with the frequency equation, which can be generally applied to this kind of systems. The exact and assumed-mode solutions including the dynamics of the base cart are obtained, and the open-loop responses of the system by arbitrarily designed forcing function are given by numerical simulations. The results match well with physical phenomena even at the extreme cases where the concentrated mass is attached to the bottom and to the top of the beam.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.26 no.3
• /
• pp.281-289
• /
• 2016

In this paper, a new dynamic model for modal analysis of a rotating cantilever beam with a tip-mass is developed. The nonlinear strain such as von Karman type and the corresponding linearized stress are used to consider the geometric nonlinearity, and Euler-Bernoulli beam theory is applied in the present model. The nonlinear equations of motion and the associated boundary conditions which include the inertia of the tip-mass are derived through Hamilton's principle. In order to investigate modal characteristics of the present model, the linearized equations of motion in the neighborhood of the equilibrium position are obtained by using perturbation technique to the nonlinear equations. Since the effect of the tip-mass is considered to the boundary condition of the flexible beam, weak forms are used to discretize the linearized equations. Compared with equations related to stiffening effect due to centrifugal force of the present and the previous model, the present model predicts the dynamic characteristic more precisely than the another model. As a result, the difference of natural frequencies loci between two models become larger as the rotating speed increases. In addition, we observed that the mode veering phenomenon occurs at the certain rotating speed.

• Steel and Composite Structures
• /
• v.45 no.4
• /
• pp.547-554
• /
• 2022

This paper studies the static stability of an axially graded column with the power-law gradient varying along the axial direction. For a nonhomogeneous column with one end linked to a rotational spring and loaded by a compressive force, respectively, an Euler problem is analyzed by solving a boundary value problem of an ordinary differential equation with varying coefficients. Buckling loads through the characteristic equation with the aid of the Bessel functions are exactly given. An alternative way to approximately determine buckling loads through the integral equation method is also presented. By comparing approximate buckling loads with the exact ones, the approximate solution is simple in form and enough accurate for varying power-law gradients. The influences of the gradient index and the rotational spring stiffness on the critical forces are elucidated. The critical force and mode shapes at buckling are presented in graph. The critical force given here may be used as a benchmark to check the accuracy and effectiveness of numerical solutions. The approximate solution provides a feasible approach to calculating the buckling loads and to assessing the loss of stability of columns in engineering.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.