• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.1 s.94
• /
• pp.29-38
• /
• 2005

The track pull-in behavior analysis in an optical disk drive (ODD) using plane phase and the evaluations for effecting parameters of it are discussed. Track pull-in, track capture procedure to do track following control, is a key factor to increase data transfer rate. First, the relative velocity between the beam spot of an optical pick-up and the target track of an optical disk is analyzed during the track pull-in procedure. In this process, it is showed that the track error signal has nonlinear characteristics which are depending on the time. Second, Runge-Kutta method to solve the nonlinear equation is applied to find the track pull-in behavior, and some optimal parameters to get stable and fast pull-in condition are obtained. Then, the phase plane analysis for track pull-in procedure is presented. Finally, some comments for the simulated results are discussed briefly.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.12
• /
• pp.3287-3296
• /
• 1994

For the numerical computation of three-dimensional compressible flow field within a blade row in a centrifugal compressor, a quasi 3-dimensional solver which combines a reversible B-B flow and an irreversible H-S flow using finite element methods was developed. In a reversible B-B flow, the governing coordinates are modified in order to be applied to any type of turbomachinery, and two kinds of stream functions are introduced in order to make the Kutta condition exactly satisfied. In an irreversible H-S flow, the changes of entropy in the irreversible governing equations are determined not by empirical source but by the theoretical treatment of dissipation forces. The dissipation forces are obtained from the distribution of shear stresses in the flow passage which are given from the wall shear stresses using the exponential functions. A more accurate quasi-3-dimensional solver is established where the effect of body forces is involved in the non-axisymmetric H-S flow. Some numerical results obtained from authors' previous studies for axial flow machines assure that the present method is able to predict well as long as the flow is subsonic and not under strong viscous effect.

• 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 Mechanical Science and Technology
• /
• v.15 no.8
• /
• pp.1165-1173
• /
• 2001

This paper deals with dynamic analysis of Pipeline Inspection Gauge (PIG) flow control in natural gas pipelines. The dynamic behaviour of PIG depends on the pressure differential generated by injected gas flow behind the tail of the PIG and expelled gas flow in front of its nose. To analyze dynamic behaviour characteristics (e.g. gas flow, the PIG position and velocity) mathematical models are derived. Tow types of nonlinear hyperbolic partial differential equations are developed for unsteady flow analysis of the PIG driving and expelled gas. Also, a non-homogeneous differential equation for dynamic analysis of the PIG is given. The nonlinear equations are solved by method of characteristics (MOC) with a regular rectangular grid under appropriate initial and boundary conditions. Runge-Kutta method is used for solving the steady flow equations to get the initial flow values and for solving the dynamic equation of the PIG. The upstream and downstream regions are divided into a number of elements of equal length. The sampling time and distance are chosen under Courant-Friedrich-Lewy (CFL) restriction. Simulation is performed with a pipeline segment in the Korea gas corporation (KOGAS) low pressure system. Ueijungboo-Sangye line. The simulation results show that the derived mathematical models and the proposed computational scheme are effective for estimating the position and velocity of the PIG with a given operational condition of pipeline.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.3 no.4
• /
• pp.241-249
• /
• 1991

A quasi-three dimensional calculation method is presented on the basis of Wu's idea using finite element methods. In B-B flow the governing equations are cast into a single equation to overcome the restriction of the type of turbomachinery, and Kutta condition is exactly assured by introducing a combination of two kinds of stream functions. In H-S flow a dissipative force which is assumed to be opposed to the relative velocity is added to the governing equation for a consistent loss model. The entropy change along each streamline is then calculated by assuming that the dissipative force may be a force coming from laminar viscous stresses with inviscid velocity distributions. Both the flow solvers are combined to build a three-dimensional flow field through a few iterations. For an effect of the distortion of H-S flow surface the body forces are computed after each B-B flow calculation is finished. Mizuki's centrifugal impellers are tested numerically. The reliability of the numerical solution compared with experimental data is guaranteed.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.3
• /
• pp.425-436
• /
• 1987

Accurate predictions of heat transfer coefficient of vertical laminar film-boiling are very important in many engineering applications. There are many predictions, however they are not exact as yet, since they have used the assumption of constant thermodynamic properties in the analysis. In this paper, heat transfer of vertical film boiling was analysized by Runnge Kutta method using veriable thermodynamic properties. 1/4 interval method was exployed for the prediction of unknown wall boundary condition. Numerical computations were performed with varying the wall temperature and the free stream velocity of liquid. Results show that assumption of constant thermodynamic properties induced considerable error in predicting the heat transfer coefficient, friction factor, film thickness, and critical length for transition to turbulent flow. Comparision of the predicted heat transfer coefficient of present analysis with that from Bromley's correlation shows that the use of general latent heat in Bromely equation instead of modified latent heat is more desireable since it makes the coefficient of Bromley equation into constant.

• Journal of the Korean Institute of Navigation
• /
• v.4 no.1
• /
• pp.31-50
• /
• 1980

It does not seem necessarily practicable to keep the system always in optimal condition, athough the control system of the follow-up mechanism on the most marine gyro compasses is to be adjusted by the operator through the gain adjustment. Sometimes a sustained oscillation or an incorrect gyro reading occurs to the system. For such a system any systematical research or theoretical basis of the guide for the optimal gain adjustment has not been reported yet. As a basic investigation of the theoretical system analysis to solve the problems concerned, the author attempts in this paper to express the system in a mathematical model deduced from the results of the theoretical approach and the experimental observation of each element contained in the follow-up mechanism of Hokshin D-1 gyro compass, and to constitute an over-all closed loop transfer function. This funciton being reverted to a fourth orderlinear differential equation, the first order simultaneous differential equations are obtained by means of the state-variables. The latter equations are solved by the Runge-Kutta method with digital computer. By comparing the characteristic of the simulated over-all output with that of the experimental result, it is shown that both outputs are nearly consistent with each other. It is also expected that the system representation proposed by this paper is valid and will be a prospective means in a further study on the design and optimal adjustment of the system.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.2 s.140
• /
• pp.113-123
• /
• 2005

A higher order panel method based on B-spline representation for both the geometry and the solution is developed for the analysis of steady flow around marine propellers. The self-influence functions due to the normal dipole and the source are desingularized through the quadratic transformation, and then shown to be evaluated using conventional numerical quadrature. By selecting a proper order for numerical quadrature, the accuracy of the present method can be increased to the machine limit. The far- and near-field influences are shown to be evaluated based on the same far-field approximation, but the near-field solution requires subdividing the panels into smaller subpanels continuously, which can be effectively implemented due to the B-spline representation of the geometry. A null pressure jump Kutta condition at the trailing edge is found to be effective in stabilizing the solution process and in predicting the correct solution. Numerical experiments indicate that the present method is robust and predicts the pressure distribution on the blade surface, including very close to the tip and trailing edge regions, with far fewer panels than existing low order panel methods.

• Structural Engineering and Mechanics
• /
• v.90 no.4
• /
• pp.357-370
• /
• 2024

Due to the fact that the mechanism of the effects of temperature and initial geometric imperfection on low-velocity impact problem of axially moving plates is not yet clear, the present paper is to fill the gap. In the present paper, the nonlinear dynamic behavior of axially moving imperfect graphene platelet reinforced metal foams (GPLRMF) plates subjected to lowvelocity impact in thermal environment is analyzed. The equivalent physical parameters of GPLRMF plates are estimated based on the Halpin-Tsai equation and the mixing rule. Combining Kirchhoff plate theory and the modified nonlinear Hertz contact theory, the nonlinear governing equations of GPLRMF plates are derived. Under the condition of simply supported boundary, the nonlinear control equation is discretized with the help of Gallekin method. The correctness of the proposed model is verified by comparison with the existing results. Finally, the time history curves of contact force and transverse center displacement are obtained by using the fourth order Runge-Kutta method. Through detailed parameter research, the effects of graphene platelet (GPL) distribution mode, foam distribution mode, GPL weight fraction, foam coefficient, axial moving speed, prestressing force, temperature changes, damping coefficient, initial geometric defect, radius and initial velocity of the impactor on the nonlinear impact problem are explored. The results indicate that temperature changes and initial geometric imperfections have significant impacts.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.31 no.4
• /
• pp.372-378
• /
• 1995

This study describes an analysis on the sinking characteristics in purse seine. The experiment was carried out using three simplified model seines in a flume tank under still water condition. The densities($\rho$) of netting materials were 0.91g/$cm^3$ for PP seine, 1.14g/$cm^3$ for PA seine and 1.38g/$cm^3$ for PES seine. Differential equations were derived from the conservation of momenta of the model seines and used to determine the sinking speeds of the depths of leadline and the other portions of seines. An analysis carried out by simultaneous differential equations for numerical method by subroutine Runge-Kutta-Gill. The results obtained were as follows: 1) Sinking speed of net margin was fastest for PP seine, followed by PA and PES seines. 2) The coefficient of resistance for netting of seine was estimated to be $K_D=0.061({\frac{\rho}{{\rho}_w}})^4$. 3) The coefficient of resistance for netting bundle of seine was estimated to be $C_R=0.91({\frac{\rho}{{\rho}_w}}$. 4) In all seines, the calculated depths of leadline closely agreed with the measured ones(meas.=0.99cal.).