• The Journal of the Acoustical Society of Korea
• /
• 제25권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.

• 한국해양공학회지
• /
• 제24권6호
• /
• pp.34-40
• /
• 2010

Bragg reflection of nonlinear waves is simulated by a 2D fully nonlinear numerical wave tank (NWT). The developed NWT was based on the Boundary Element Method (BEM) with potential theory and the mixed Eulerian-Lagrangian (MEL) time marching scheme with Runge-Kutta 4th-order time integration. A spatial variation of wave elevations and their Fourier amplitudes of each component are compared to investigate the effect of sea bottom ripples and their relative heights. The incident waves over an undulated sea bottom are partially reflected and changed to partial standing waves due to Bragg reflection. The present results are verified with linear calculations and experimental data. It is found that the 1st-order wave component is mainly affected by Bragg reflection and its spatial modulation is significant in front of the bottom ripples.

• Steel and Composite Structures
• /
• 제27권3호
• /
• pp.273-283
• /
• 2018

Nonlinear low velocity impact response of sandwich beam with laminated composite face sheets and soft core is studied based on Extended High Order Sandwich Panel Theory (EHSAPT). The face sheets follow the Third order shear deformation beam theory (TSDT) that has hitherto not reported in conventional EHSAPT. Besides, the two dimensional elasticity is used for the core. The nonlinear Von Karman type relations for strains of face sheets and the core are adopted. Contact force between the impactor and the beam is obtained using the modified Hertz law. The field equations are derived via the Ritz based applied to the total energy of the system. The solution is obtained in the time domain by implementing the well-known Runge-Kutta method. The effects of boundary conditions, core-to-face sheet thickness ratio, initial velocity of the impactor, the impactor mass and position of the impactor are studied in detail. It is found that each of these parameters have significant effect on the impact characteristics which should be considered. Finally, some low velocity impact tests have been carried out by Drop Hammer Testing Machine. The contact force histories predicted by EHSAPT are in good agreement with that obtained by experimental results.

• 제어로봇시스템학회논문지
• /
• 제3권1호
• /
• pp.67-76
• /
• 1997

In this paper we developed a tactile feedback device using repulsive force of magnets. The force of the tactile feedback device was derived from the Maxwell's stress method by using the concept of magnetic charge. Magnetic repulsive force is linear function with respect to current and nonlinear to displacement. Experimental data shows these characteristics. To compensate the fact that the presented tactile feedback device can not be controlled by close loop control, we developed a simulation model which predicts output displacement and force by using Runge-Kutta method. And, this paper evaluated the presented tactile feedback device and compared it with commercial tactile feedback devices.

• 한국환경과학회지
• /
• 제6권4호
• /
• pp.399-407
• /
• 1997

The pyrolysis reactions of atomic hydrogen with chloroform were studied In a 4 cm 1.6. tubular flow reactor with low flow velocity 1518 cm/sec and a 2.6 cm 1.4. tubular flow reactor with high flow velocity (1227 cm/sec). The hydrogen atom concentration was measured by chemiluminescence titration with nitrogen dioxide, and the chloroform concentrations were determined using a gas chromatography. The chloroform conversion efficiency depended on both the chloroform flow rate and linear flow velocity, but 416 not depend on the flow rate of hydrogen atom. A computer model was employed to estimate a rate constant for the initial reaction of atomic hydrogen with chloroform. The model consisted of a scheme for chloroform-hydrogen atom reaction, Runge-Kutta 4th-order method for Integration of first-order differential equations describing the time dependence of the concentrations of various chemical species, and Rosenbrock method for optimization to match model and experimental results. The scheme for chloroform-hydrogen atom reaction Included 22 elementary reactions. The rate constant estimated using the data obtained from the 2.6 cm 1.4. reactor was to be 8.1 $\times$ $10^{-14}$ $cm^3$/molecule-sec and 3.8 $\times$ $10^{-15}$ cms/molecule-sec, and the deviations of computer model from experimental results were 9% and 12% , for the each reaction time of 0.028 sec and 0.072 sec, respectively.

• Journal of Mechanical Science and Technology
• /
• 제19권spc1호
• /
• pp.388-394
• /
• 2005

Recently, the authors have developed a method for real-time dynamics of multibody systems, which combines a semi-recursive formulation to derive the equations of motion in dependent relative coordinates, along with an augmented Lagrangian technique to impose the loop closure conditions. The following numerical integration procedures, which can be grouped into the so-called structural integrators, were tested : trapezoidal rule, Newmark dissipative schemes, HHT rule, and the Generalized-${\alpha}$ family. It was shown that, for large multi body systems, Newmark dissipative was the best election since, provided that the adequate parameters were chosen, excellent behavior was achieved in terms of efficiency and robustness with acceptable levels of accuracy. In the present paper, the performance of the described method in combination with another group of integrators, the Implicit Runge-Kutta family (IRK), is analyzed. The purpose is to clarify which kind of IRK algorithms can be more suitable for real-time applications, and to see whether they can be competitive with the already tested structural family of integrators. The final objective of the work is to provide some practical criteria for those interested in achieving real-time performance for large and complex multibody systems.

• Journal of Advanced Marine Engineering and Technology
• /
• 제20권5호
• /
• pp.82-92
• /
• 1996

A computer simulation with predict the fuel injection rates and the fuel injection pressure behaviors in diesel engine fuel injection systems would by very useful in designing or improving fuel injection systems. In this paper we developed computer program in order to predict the behaviors of the fuel injection rate and the injection pressure for Electronic Hydraulic Ultra-High Pressure Fuel Injection System. We've applied the continuity and momentum equations for the hydraulic phenomena and the dynamics of individual components of the Electronic Hydraulic Fuel Injection System. To solve all the equations numerically we've applied the Runge-kutta IV method. Water hammer equations were applied for the hydraulic pipe solution, and the method of characteristics was employed in our calculations. The simulation results were compared with the experimental results for: Accumulator pressure, Injection pressure and unjection rate. As a result, The simulation results agree very well with our experimental results. We found that a large accumulator and the high speed solenoid valve were required, and the compression volume of the fuel had to be as small as possible in order to acheive ultra-high pressure fuel injection.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.145-155
• /
• 1998

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 추계학술대회
• /
• pp.998-1003
• /
• 2003

Ultrasonic flow metering(UFM) technology is being received much attention from a variety of industrial fields to exactly measure the flow rate. The UFM has much advantage over other conventional flow meter systems, since it has no moving parts, and offers good accuracy and reliability without giving any disturbances to measure the flow rate, thereby not causing pressure losses in the flow fields. In the present study, 3-dimensional, unsteady, compressible Navier-Stokes equations are solved by a finite volume scheme, based upon the second order upwind scheme for spatial derivatives and the multi-stage Runge-Kutta integral method for time derivatives. In order to simulate multi-path ultrasonic flow meter, an excited pressure signal is applied to three different locations upstream, and the pressure signals are received at three different locations downstream. The mean flow velocities are calculated by the time difference between upstream and downstream propagating pressure signals. The obtained results show that the present CFD method simulates successfully ultrasonic meter gas flow and the mean velocity measured along the chord near the wall is considerably influenced by the boundary layers.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 추계학술대회
• /
• pp.169-174
• /
• 2003

The thermal design of the bus bar of a Gas-Insulated Switchgear(GIS) becomes important since the current-carrying capacity of the GIS is limited by maximum operating temperature. In order to predict temperature rise of the bus bar, a program has been developed. Various heat sources possibly generated in the bus bar are calculated in the program. To estimate temperature rises at the bus bar caused by the heat balance between the heat generation and heat transfer, the finite volume method as well as the $4^{th}$ order Runge-Kutta method has been employed. In the experiments, temperature rises at conductor, contact part and external tank are measured for full-scale gas-insulated bus bars. The comparisons of the predicted values of the heat balance calculation to those of the experiments are made. From the comparisons, it is concluded that the developed program can predict the temperature rise of the bus bar quite well.