• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.3
• /
• pp.166-177
• /
• 2005

This study developed the jet integral model to analyze the behavior of the wastewater discharge in the near field using the fourth order Runge-Kutta method in order to numerically solve the problems of six ordinary differential equations and six unknowns. This jet integral model used the entrainment hypothesis and the manipulation of sonle shape constant. This study also conducted the hydraulic experiments fnr single horizontal buoyant Jet using LIF through the calibration procedure. The results calculated by the previous models, CORMIX 1 and VISJET, and the proposed jet integral model were compared to the hydraulic experimental results. The centerline trajectories predicted by the proposed model were in good agreements with the experimental results in the transition region whereas the trajectories calculated by the VISJET model agreed well with the measured data in the momentum and buoyancy-dominated regions. The centerline dilution calculated by the proposed model agreed generally with the measured dilution in the intial and transition regions while the centerline dilution predicted by the CORMIX 1 was in good agreements with the experimental results in the momentum and buoyancy-dominated regions.

• Structural Engineering and Mechanics
• /
• v.12 no.6
• /
• pp.657-668
• /
• 2001

Fiber reinforced cementitious composites are nowadays widely applied in civil engineering. The postcracking performance of this material depends on the interaction between a steel fiber, which is obliquely across a crack, and its surrounding matrix. While the partly debonded steel fiber is subjected to pulling out from the matrix and simultaneously subjected to transverse force, it may be modelled as a Bernoulli-Euler beam partly supported on an elastic foundation with non-linearly varying modulus. The fiber bridging the crack may be cut into two parts to simplify the problem (Leung and Li 1992). To obtain the transverse displacement at the cut end of the fiber (Fig. 1), it is convenient to directly solve the corresponding differential equation. At the first glance, it is a classical beam on foundation problem. However, the differential equation is not analytically solvable due to the non-linear distribution of the foundation stiffness. Moreover, since the second order deformation effect is included, the boundary conditions become complex and hence conventional numerical tools such as the spline or difference methods may not be sufficient. In this study, moment equilibrium is the basis for formulation of the fundamental differential equation for the beam (Timoshenko 1956). For the cantilever part of the beam, direct integration is performed. For the non-linearly supported part, a transformation is carried out to reduce the higher order differential equation into one order simultaneous equations. The Runge-Kutta technique is employed for the solution within the boundary domain. Finally, multi-dimensional optimization approaches are carefully tested and applied to find the boundary values that are of interest. The numerical solution procedure is demonstrated to be stable and convergent.

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

• Wind and Structures
• /
• v.23 no.5
• /
• pp.421-447
• /
• 2016

Wind-induced and earthquake-induced excitations on tall structures can be effectively controlled by Tuned Liquid Damper (TLD). This work presents a numerical simulation procedure to study the performance of tuned liquid tank- structure system through ${\sigma}$-transformation based fluid-structure coupled solver. For this, a 'C' based computational code is developed. Structural equations are coupled with fluid equations in order to achieve the transfer of sloshing forces to structure for damping. Structural equations are solved by fourth order Runge-Kutta method while fluid equations are solved using finite difference based sigma transformed algorithm. Code is validated with previously published results. The minimum displacement of structure is observed when the resonance condition of the coupled system is satisfied through proper tuning of TLD. Since real-time excitations are random in nature, the performance study of TLD under random excitation is also carried out in which the Bretschneider spectrum is used to generate the random input wave.

• Proceedings of the KIEE Conference
• /
• 1996.07a
• /
• pp.104-106
• /
• 1996

In this paper, the dynamic motion driven by electromagnetic force of transformer windings is modeled and its characteristics are numerically analyzed. The electromagnetic field is obtained using the 2D finite element method taking account of anisotropic property of iron core, and the electromagnetic force on the transformer winding is calculated from Lorenz's force formula using the field distribution result. The system motion equation driven by electromagnetic force and gravitational force is numerically analyzed using the 4-order Runge-Kutta algorithm. Above analyses procedure is applied to a single-phase core-type transformer to validate its algorithm.

• Structural Engineering and Mechanics
• /
• v.61 no.5
• /
• pp.657-661
• /
• 2017

In this paper, it has been tried to propose a new semi analytical approach for solving nonlinear vibration of conservative systems. Hamiltonian approach is presented and applied to high nonlinear vibration systems. Hamiltonian approach leads us to high accurate solution using only one iteration. The method doesn't need any small perturbation and sufficiently accurate to both linear and nonlinear problems in engineering. The results are compared with numerical solution using Runge-Kutta-algorithm. The procedure of numerical solution are presented in detail. Hamiltonian approach could be simply apply to other powerfully non-natural oscillations and it could be found widely feasible in engineering and science.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2004.05a
• /
• pp.894-899
• /
• 2004

In this paper, the track pull-in behavior analyses in an optical disk drive (ODD) using plane phase is treated and the parameters affecting it are discussed. Track pull-in is the track capture procedure to do track following control and it is key factor to increase data transfer rate. Simulation method, Runge-Kutta method to solve nonlinear equation, is used to evaluate the track pull-in conditions, and the real servo loop parameters are applied in this process to get the more real condition. Finally, the comments for the acquired results are discussed briefly.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.5
• /
• pp.621-626
• /
• 1985

Unsteady groundwater flow in a three-layer unconfined aquifer has been studied theoretically and experimentally. Two different methods have been used in solving the governing equations of the flow, the nonlinear partial differential equations; (1) The governing equations are linearized for each layer and approximate solutions are obtained. (2) The governing equations are transformed to nonlinear ordinary differential equations, which are solved numerically by Runge-Kutta procedure. Fine, middle sized and coarse sands are used in the experiments. It is found that the solutions from the method(2) ( the reduction of partial differential equations to ordinary differential equations) give better agreement with the experimental results than the solution from the method(1).

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.11
• /
• pp.3031-3040
• /
• 1995

Two-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with the unstructured triangular meshes. The pressure-velocity coupling is handled by the artificial compressibility algorithm due to its computational efficiency associated with the hyperbolic nature of the resulting equations. The convective fluxes are obtained by the Roe's flux difference splitting scheme using edge-based connectivities and higher-order differences are achieved by a reconstruction procedure. The time integration is based on an explicit four-stage Runge-Kutta scheme. Numerical procedures with local time stepping and implicit residual smoothing have been implemented to accelerate the convergence for the steady-state solutions. Comparisons with experimental data and other numerical results have proven accuracy and efficiency of the present unstructured approach.

• Structural Engineering and Mechanics
• /
• v.42 no.4
• /
• pp.449-469
• /
• 2012

For generally damped linear systems with repeated eigenvalues and defective eigenvectors, this study provides a decomposition method based on residue matrix, which is suitable for engineering applications. Based on this method, a hybrid approach is presented, incorporating the merits of the modal superposition method and the residue matrix decomposition method, which does not need to consider the defective characteristics of the eigenvectors corresponding to repeated eigenvalues. The method derived in this study has clear physical concepts and is easily to be understood and mastered by engineering designers. Furthermore, this study analyzes the applicability of step-by-step methods, including the Newmark beta and Runge-Kutta methods for dynamic response calculation of defective systems. Finally, the implementation procedure of the proposed hybrid approach is illustrated by analyzing numerical examples, and the correctness and the effectiveness of the formula are judged by comparing the results obtained from the different methods.