• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.2
• /
• pp.1-13
• /
• 1986

The nonlinear hydrodynamic forces acting on a two-dimensional circular cylinder, oscillating with large amplitude in the free surface, are calculated by using the Semi-Lagrangian Time-Step-ping Method used by O.M. Faltinsen. In present calculation the position and the potential value of free surface are calculated using the exact kinematic and dynamic free surface boundary condition. At each time step an integral equation is solved to obtain the value of potential and normal velocity along the boundaries, consisting of both the body surface and the free surface. Some effort was devoted to the elimination of instability arising in the range of high frequency. Numerical simulations were performed up to the 3rd or 4th period which seems to be enough for the transient effect to die out. Each harmonic component and time-mean force are obtained by the Fourier transform of forces in time domain. The results are compared with others' experimental and theoretical results. Particularly, the calculation shows the tendency that the acceleration-phase 1st-harmonic component(added mass) increases as the motion amplitude increases and a reverse tendency in the velocity-phase 1st-harmonic component(damping coefficient). The Yamashita's experimental result also shows the same tendency. In general, the present result show relatively good agreement with the Yamashita's experimental result except for the time-mean force.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.249-252
• /
• 2003

Acoustic behavior in combustion chamber with acoustic cavity is numerically investigated by adopting linear acoustic analysis. Helmholtz-type resonator is employed as a cavity model to suppress acoustic instability. The tuning frequency of acoustic cavity is adjusted by varying the sound speed in acoustic cavity. Acoustic pressure responses of chamber to acoustic oscillating excitation are shown md acoustic damping effect of acoustic cavity is quantified by damping factor. As the tuning frequency approaches the target frequency of the resonant mode, mode split from the original resonant mode to lower and upper modes appears and thereby damping effect is degraded. Considering mode split and damping effect as a function of tuning frequency, it is desirable to make acoustic cavity tuned to maximum frequency of those of the possible splitted upper modes.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.1 s.71
• /
• pp.39-46
• /
• 2006

This paper analyzes the continuous Galerkin method for the space-time discretization of wave equation. The method of space-time finite elements enables the simple solution than the usual finite element analysis with discretization in space only. We present a discretization technique in which finite element approximations are used in time and space simultaneously for a relatively large time period called a time slab. The weighted residual process is used to formulate a finite element method for a space-time domain. Instability is caused by a too large time step in successive time steps. A stability problem is described and some investigations for chosen types of rectangular space-time finite elements are carried out. Some numerical examples prove the efficiency of the described method under determined limitations.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.2
• /
• pp.21-28
• /
• 1987

The paper describes the analysis of large displacement problems including instability phenomena. The element used in this is a degenerated isoparametric shell element with eight nodes. Total Lagrangian formulation has been adopted in this study using Newton-Raphson iteration method with incremental load. The linear stability analyses performed usually for the initial position can be repeated at several advanced fundamental states on the non-linear buckling path. Thus a current estimate of the failure load is given. The numerical examples of a cylindrical panel under uniform load, simply supported plate under axial load, and clamped plate under uniform load are carried out. The examples applying degenerated isoparametric elements to bifurcation buckling and nonlinear collapse problems are also performed.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.369-374
• /
• 2010

This paper presents a numerical investigation of the Deflagration to Detonation Transition (DDT) of flame acceleration by a shock wave in combustible gas mixture. A model consisting of the reactive compressible Navier-Stokes equations is used. The effects of viscosity, thermal conduction, species diffusion, and chemical reactions are included. Using this model, the generation of hot spots by repeated shock and flame interaction in front and back of flame and the change of detonation occurrence by various shock intensities (Ms=1.1, 1.2, 1.3) are studied. The simulations show that as the incident shock intensity increases, the Richtmyer-Meshkov (RM) instability becomes stronger and DDT occurrence time is reduced.

• Structural Engineering and Mechanics
• /
• v.52 no.5
• /
• pp.903-918
• /
• 2014

Large space structures may have resonant low eigenvalues and often these appear with closely-spaced natural frequencies. Owing to the coupling among modes with closely-spaced natural frequencies, each eigenvector corresponding to closely-spaced eigenvalues is ill-conditioned that may cause structural instability. The subspace to an invariant subspace corresponding to closely-spaced eigenvalues is well-conditioned, so a method is presented to design the feedback control law of intelligent structures with closely-spaced eigenvalues in this paper. The main steps are as follows: firstly, the system with closely-spaced eigenvalues is transformed into that with repeated eigenvalues by the spectral decomposition method; secondly, the computation for the linear combination of eigenvectors corresponding to repeated eigenvalues is obtained; thirdly, the feedback control law is designed on the basis of the system with repeated eigenvalues; fourthly, the system with closely-spaced eigenvalues is regarded as perturbed system on the basis of the system with repeated eigenvalues; finally, the feedback control law is applied to the original system, the first order perturbations of eigenvalues are discussed when the parameter modifications of the system are introduced. Numerical examples are given to demonstrate the application of the present method.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1915-1920
• /
• 2004

An edge tone is the discrete tone or narrow-band sound produced by an oscillating free shear layer impinging on a rigid surface. In this paper we present a two-dimensional edge tone to predict the frequency characteristics of the discrete oscillations of a jet-edge feedback cycle by the finite difference lattice Boltzmann method. We use a new lattice BGK compressible fluid model that has an additional term and allow larger time increment comparing a conventional FDLB model, and also use a boundary fitted coordinates. The jet is chosen long enough in order to guarantee the parabolic velocity profile of the jet at the outlet, and the edge consists of a wedge with an angle of ${\alpha}=23^{\circ}$ . At a stand-off distance ${\omega}$ , the edge is inserted along the centreline of the jet, and a sinuous instability wave with real frequency f is assumed to be created in the vicinity of the nozzle and to propagate towards the downstream. We have succeeded in capturing very small pressure fluctuations result from periodically oscillation of jet around the edge. That pressure fluctuations propagate with the sound speed. Its interaction with the wedge produces an irrotational feedback field which, near the nozzle exit, is a periodic transverse flow producing the singularities at the nozzle lips.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.12 no.8
• /
• pp.771-781
• /
• 2000

The Karman boundary-layer, has been numerically investigated for the disturbance wave number, wave velocity, azimuth angle and radius (Reynolds number, Re). The disturbed flow over rotating disk can lead to transition at a much lower Re than that of the well-known Type 1 mode of instability. This early transition is due to the excitation of the Type II mode. Presented are the neutral stability results concerning these modes by solving new formulated vorticity equations with consideration of whole convective terms. When the present numerical results are compared with the previously known results, the value of critical Re corresponding to Type I is moved from Rec,! =285.3 to 270.2 and the value corresponding to Type II is from $Re_{c,2}$=69.4 to 36.9, respectively. Also, the corresponding wave number is moved from $k_1$ =0.378 to $k_1$ =0.389 for Type I; from $k_2$ =0.279 to $k_2$=0.385 for Type II. For Type II, the upper limit of wave number and azimuth angle is $k_U$=0.5872,$varepsilon_U=-18^{\circ}$ , while its lower limit is$k_L$ =0.05, $varepsilon_L=-27^{\circ}$ This implies that the disturbances will be relatively fast amplified at small Re and within narrow bands of wave number compared with the previous results.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.62-68
• /
• 2004

A comprehensive numerical analysis has been carried out for both non-reacting and reacting flows in a scramjet engine combustor with and without a cavity. The theoretical formulation treats the complete conservation equations of chemically reacting flows with finite-rate chemistry of hydrogen-air. Turbulence closure is achieved by means of a k-$\omega$ two-equation model. The governing equations are discretized using a MUSCL-type TVD scheme, and temporally integrated by a second-order accurate implicit scheme. Transverse injection of hydrogen is considered over a broad range of injection pressure. The corresponding equivalence ratio of the overall fuel/air mixture ranges from 0.167 to 0.50. The work features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous studies. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the .underlying physical mechanisms. Much of the flow unsteadiness is related not only to the cavity, but also to the intrinsic unsteadiness in the flow-field. The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The roles of the cavity, injection pressure, and heat release in determining the flow dynamics are examined systematically.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.251-253
• /
• 2017

Numerical study of air-blast type injector for low emission aircraft engines was conducted. Volume-of-fluids approach was used to track interface of fuel and air. Primary atomization of fuel stream was visualized, and thickness and mean velocity at the injector exit was calculated. Liquid fuel injected from fuel slots joined together as a thin film on preflimer surface, and interacted with swirling air. As instability on the fuel surface increased, separation of fuel as ligaments and droplets occured. The film thickness and velocity were used to as fuel injection boundary condition for reactive flow simulation. Primary reaction zone was formed in vicinity of the fuel nozzle, creating a stable flame inside the combustor.