• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.23-28
• /
• 1999

A mutigrid convergence acceleration technique is presented for computing hypersonic inviscid and viscous flows in equilibrium state. The governing equations are solved using an explicit Runge-Kutta method. Curve fitting data in NASA Reference Publication 1181, 1260 are used to calculate equilibrium properties. In order to ensure stability, damped prolongation and modified implicit residual smoothing are proposed. Blunt body test cases are presented to demonstrate the robustness and the efficiency in performance of the proposed methods

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

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

• Journal of applied mathematics & informatics
• /
• v.31 no.1_2
• /
• pp.131-145
• /
• 2013

In this paper, we present an initial value technique for solving singularly perturbed differential difference equations with a boundary layer at one end point. Taylor's series is used to tackle the terms containing shift provided the shift is of small order of singular perturbation parameter and obtained a singularly perturbed boundary value problem. This singularly perturbed boundary value problem is replaced by a pair of initial value problems. Classical fourth order Runge-Kutta method is used to solve these initial value problems. The effect of small shift on the boundary layer solution in both the cases, i.e., the boundary layer on the left side as well as the right side is discussed by considering numerical experiments. Several numerical examples are solved to demonstate the applicability of the method.

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

• Journal of Convergence for Information Technology
• /
• v.10 no.7
• /
• pp.131-137
• /
• 2020

The characteristics of forced vibration with excited sinusoidal force was introduced. Also, numerical analyses and FRF in frequency domain were performed in detail. In this regard, the responses of displacement, velocity and acceleration were investigated in a forced vibration model. The FRF characteristics in real and imaginary part around natural frequency are also discussed. This response approach of forced vibration in time domain is used for the identification and monitoring of sinusoidal forced vibration. For acquiring a displacement, velocity and acceleration, a numerical technique of Runge-Kutta-Gill method was performed. For the FRF(frequency response function), These responses are used. Also, the FRF can represent the intrinsic characteristics of the forced vibration. These performed results and analysis are successful in each damped condition for the forced vibration model. After numerical analysis of the different mass, damping and stiffness, the forced vibration response characteristics with sinusoidal force was discriminated considering its amplitude and frequency simultaneously.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.41 no.1
• /
• pp.39-48
• /
• 2021

In this paper, an analytical model is proposed for assessing the natural frequency of barrettes subjected to vertical loading. The differential equation governing the free vibration of rectangular friction piles embedded in inhomogeneous soil is derived. The governing equation is numerically integrated by Runge-Kutta technique and the eigenvalue of natural frequency is computed by Regula-Falsi method. The numerical solutions for the natural frequency of barrettes compare well with those obtained from finite element analysis. Illustrated examples show that the natural frequencies increase with an increase of the cross-sectional aspect ratio, the friction resistance ratio and the soil stiffness ratio, and decrease with an increase of the friction aspect ratio, the slenderness ratio and the load factor, respectively.

• Journal of Convergence for Information Technology
• /
• v.11 no.3
• /
• pp.132-139
• /
• 2021

The response characteristics of the forced vibration generated when the high-damped vehicle pass the bumped barrier was studied, and in particular, the response behavior of displacement, velocity and acceleration was analyzed for the forced vibration model. In addition, in order to obtain responses such as displacement, velocity, and acceleration, a numerical analysis technique of the Runge-Kutta-Gill method was performed in time domain. The response was successfully obtained in detail under several high damping conditions. As a numerical analysis result, the response of the vehicle was obtained by considering the vehicle body to which the impulse impact was applied. Also, the analysis result was compared with the experimental result in order to verify the validity of vehicle model. The amplitude and natural frequency of the vehicle were considered and analyzed. The Nyquist diagram of the vehicle model was also obtained and the relationship could be analyzed. And the vibration response was analyzed on different mass, damping and stiffness.

• Journal of the Korean institute of surface engineering
• /
• v.41 no.5
• /
• pp.205-213
• /
• 2008

A 3D particle code is developed to analyze electron behavior in a planar magnetron sputtering cathode either in balanced or unbalanced configuration. Three types of collisions are included; electron - neutral elastic, excitation to a metastable state and ionization. Flight path is calculated by a 4-th order Runge-Kutta method with a time step of 10 ps. Effects of electron starting position, magnetic field intensity and configuration were analyzed. For a more efficient and accurate modeling, multithreading technique is considered for multicore CPU computers. Under an assumption of cold ion approach, target erosion profiles are predicted for a flat target surface.

• Journal of Computational Design and Engineering
• /
• v.3 no.4
• /
• pp.330-336
• /
• 2016

The present work is concerned with heat and mass transfer of an electrically conducting second grade MHD fluid past a semi-infinite stretching sheet with convective surface heat flux. The analysis accounts for thermophoresis and thermal radiation. A similarity transformations is used to reduce the governing equations into a dimensionless form. The local similarity equations are derived and solved using Nachtsheim-Swigert shooting iteration technique together with Runge-Kutta sixth order integration scheme. Results for various flow characteristics are presented through graphs and tables delineating the effect of various parameters characterizing the flow. Our analysis explores that the rate of heat transfer enhances with increasing the values of the surface convection parameter. Also the fluid velocity and temperature in the boundary layer region rise significantly for increasing the values of thermal radiation parameter.