• Proceedings of the Korean Institute of IIIuminating and Electrical Installation Engineers Conference
• /
• 2003.11a
• /
• pp.233-235
• /
• 2003

The voltage and current of fluoescent lamp has been successfully calculated in the case of frequency variation. Energy states of mercury atom in the discharge process are regarded as six levels. These calculations have been accurately solved by numerically employing mixed the FDM and the 2nd Runge-Kutta method. The theoretical calculation results and experimental results were presented to verify the feasibility of the modeling. Calculation and experimental results were presented to verify the feasibility of the calculation.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.22 no.5
• /
• pp.143-150
• /
• 2021

To develop the technique for predicting the horizontal displacement of a retaining wall induced by an earthquake, an equation of motion that depicts the retaining wall-soil vibrating system was derived. The resulting differential equation was solved using the Runge-Kutta-Nystr?m method. Considering the pre-mentioned derivation process, the analysis procedures for obtaining horizontal displacement induced by an earthquake were programmed. The core algorithm of the displacement-force relationship, which is the main engine of the developed program, was suggested. Considering the results obtained by adopting the developed program to the assumed retaining wall under an earthquake, the relationships between the time-displacement, time-force, and displacement-force were reasonable. According to the results computed by the program, the displacements to the front direction of the wall occurred, and the displacement per cycle converged after some cycles elapsed. Displacements with a natural period were calculated, which showed that the maximum displacement was observed when the natural frequency was slightly different from the excitation frequency rather than the same values of the two frequencies. This happens because the vibrating system was modeled by two springs with different stiffness.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.5
• /
• pp.1383-1393
• /
• 2014

Recently, with rapid improvement in computer hardware and theoretical development in the field of computational fluid dynamics, high-order accurate schemes also have been applied in the realm of computational hydraulics. In this study, numerical solutions of 1D shallow water equations are presented with TVD Runge-Kutta discontinuous Galerkin (RKDG) finite element method. The transcritical flows such as dam-break flows due to instant dam failure and transcritical flow with bottom elevation change were studied. As a formulation of approximate Riemann solver, the local Lax-Friedrichs (LLF), Roe, HLL flux schemes were employed and MUSCL slope limiter was used to eliminate unnecessary numerical oscillations. The developed model was applied to 1D dam break and transcritical flow. The results were compared to the exact solutions and experimental data.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.4
• /
• pp.251-258
• /
• 1993

The main purpose of this paper is to present both the natural frequencies and the buckling loads of beam-columns on Winkler-type foundations. The ordinary differential equations governing the free vibrations and the buckling loads of beam-columns on Winkler-type foundation are derived as nondimensional forms. The Runge-Kutta method and Determinant Search method are used to perform the integration of the differential equations and to determine the eigenvalues(natural frequencies and buckling loads), respectively. Hinged-hinged and damped-clamped end constraints are applied in numerical examples. The relation between frequency parameter and elastic foundation parameter is presented in figure. The effects of axial loads on the natural frequencies of beam-columns on elastic foundations are investigated and the relation between buckling load parameter and elastic foundation parameter is also analyzed. The relation between foundation rested ratio and frequency parameter, buckling load parameter are investigated. The beam-columns on non-homogeneous elastic foundation are analyzed and typical mode shapes are also presented.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.3
• /
• pp.101-109
• /
• 1987

The governing differential equations for the free vibration of general arch are derived including the effect of rotary inertia in addition to the usual actions. These differential equations are applied to the sinusoidal arch and the numerical methods are developed to analyze these equations. A trial eigenvalue method and the Runge-Kutta method are used to determine the natural frequencies and to perform the integration of the differential equations, respectively. A detailed studies are made of the lowest three vibration frequencies for hinged arches with the span length equal to 10 m. The effect of the rotary inertia is analyzed. And as the numerical results the frequency versus the rise of arch and the radius of gyration are presented in figures.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.1
• /
• pp.1-8
• /
• 2001

The computations of the turbulent flow around the ship models with the free-surface effects were carried out. Incompressible Reynolds-Averaged Navier-Stokes equations were solved by using an explicit finite-difference method with the nonstaggered grid system. The method employed second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration. For the turbulence closure, a modified Baldwin-Lomax model was exploited. The location of the free surface was determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and a free-surface conforming grid was generated at each time step so that one of the grid boundary surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition was applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method developed in the present study, the computations were carried out for beth Wigley and Series 60 $C_B=0.6$ ship model and the computational results showed good agreements with the experimental data.

• Steel and Composite Structures
• /
• v.36 no.5
• /
• pp.603-615
• /
• 2020

The present manuscript focuses on the flow and heat transfer of the dusty fluid along exponentially stretching cylinder. Enormous attempts are made for fluid flow along cylinder but the study of fluid behavior along exponentially stretching cylinder is discussed lately. Using appropriate transformations, the governing partial differential equations are converted to non-dimensional ordinary differential equations. The transformed equations are solved numerically using Shooting technique with Runge-Kutta method. The influence of the physical parameters on the velocity and temperature profiles as well as the skin fraction coefficient and the local Nusselt number are examined in detail. The essential observations are as the fluid velocity decreases but temperature grows with rise in particle interaction parameter, and both the fluid velocity and temperature fall with increase in mass concentration parameter, Reynold number, Particle interaction parameter for temperature and the Prandtl number.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Journal of KSNVE
• /
• v.10 no.6
• /
• pp.1067-1074
• /
• 2000

Numerical methods for calculating both the natural frequencies and buckling loads of columns with the multiple elastic springs are developed. In order to derive the governing equations of such columns, each elastic spring is modeled as a discrete elastic foundation with the finite longitudinal length. By using this model, the differential equations governing both the free vibrations and buckled shapes, respectively, of such columns are derided. These differential equations are solved numerically. The Runge- Kutta method is used to integrate the differential equations, and the determinant search method combined with Regula-Falsi method is used to determine the eingenvalues. namely natural frequencies and buckling loads. In the numerical examples, the clamped-clamped. clamped-hinged, hinged-clamped and hinged-hinged end constraints are considered. Extensive numerical results including the frequency parameters, mode shapes of free vibrations and buckling load parameters are presented in the non-dimensional forms.

• Journal of Mechanical Science and Technology
• /
• v.20 no.8
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.