• Journal of applied mathematics & informatics
• /
• v.7 no.3
• /
• pp.941-958
• /
• 2000

Euler method is generalized to solve the system of nonlinear differential equations. The generalization is carried out by taking a special constant matrix S so that exp(tS) can be exactly computed. Such a matrix S is extracted from the Jacobian matrix of the given problem. Stability of the generalized Euler process is discussed. It is shown that the generalized Euler process is comparable to the fourth order Runge-Kutta method. We also exemplify that the important qualitative and geometric features of the underlying dynamical system can be recovered by the generalized Euler process.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.1
• /
• pp.3-12
• /
• 1986

This paper investigates the static configurations and the dynamic behaviors of a single point mooring line. To obtain the static configuration and static tension distribution along the mooring line, including the effect of fluid nonlinear drag and the elasticity of the line, the Runge-Kutta fourth order numerical method was used. The relationship between the horizontal excursion and the horizontal restoring force component of the mooring line, which is very important to a mooring line design, and the effect of a subsurface buoy on the static configuration are presented. In nonlinear dynamic analysis including nonlinear fluid drag acting on the line and geometrical nonlinearity for large deflections, finite element method using updated Lagrangian was used to obtain the solution. In the case of upper end harmonic excitation of the mooring line, the dynamic motion and the tension were also presented.

• Journal of Korea Multimedia Society
• /
• v.14 no.5
• /
• pp.642-649
• /
• 2011

Among the various aspects of offline handwritten character patterns, it is the great variety of writing styles and variations that renders the task of computer recognition very hard. The immense variety of character shape has been recognized but rarely studied during the past decades of numerous research efforts. This paper tries to address the problem of measuring image distortions and handwritten character patterns with respect to reference patterns. This work is based on mass-spring mesh model with the introduction of simulated electric charge as a source of the external force that can aid decoding the shape distortion. Given an input image and a reference image, the charge is defined, and then the relaxation procedure goes to find the optimum configuration of shape or patterns of least potential. The relaxation process is based on the fourth order Runge-Kutta algorithm, well-known for numerical integration. The proposed method of modeling is rigorous mathematically and leads to interesting results. Additional feature of the method is the global affine transformation that helps analyzing distortion and finding a good match by removing a large scale linear disparity between two images.

• Computers and Concrete
• /
• v.27 no.4
• /
• pp.369-383
• /
• 2021

This article deals with the frequency analysis of viscoelastic sandwich disk with graphene nano-platelets (GPLs) reinforced viscoelastic concrete (GPLRVC) face sheets and honeycomb core. The honeycomb core is made of aluminum due to its low weight and high stiffness. The rule of the mixture and modified Halpin-Tsai model are engaged to provide the effective material constant of the concrete. By employing Hamilton's principle, the governing equations of the structure are derived and solved with the aid of the Generalize Differential Quadrature Method (GDQM). In this paper, viscoelastic properties are modeled according to Kelvin-Voigt viscoelasticity. The deflection as the function of time can be solved by the fourth-order Runge-Kutta numerical method. Afterward, a parametric study is carried out to investigate the effects of the outer to inner radius ratio, hexagonal core angle, thickness to length ratio of the concrete, the weight fraction of GPLs into concrete, and the thickness of honeycomb core to inner radius ratio on the frequency of the viscoelastic sandwich disk with honeycomb core and FG-GPLRVC face sheet.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.7
• /
• pp.9-17
• /
• 2005

The wake patterns and thrust characteristics of dual flapping airfoils in a biplane configuration are investigated using an unsteady panel method. To trace complicated wake shapes behind airfoils, a core addition scheme, a vortex core model, and the fourth order Runge-Kutta convection scheme are employed. Present results are verified by comparing them with flow visualization, exact solution and published computed results. The thickness and camber of thick airfoils has an effect of decreasing thrust. The airfoils produce maximum thrust when the phase angles between plunging and pitching motions are both 90 and 120 degrees. Thrust increases as the plunge velocity is increased, which is also found as the pitch amplitude is stepped up. Thrust decreases when the distance between the airfoils is less than 0.6c.

• Earthquakes and Structures
• /
• v.24 no.3
• /
• pp.217-235
• /
• 2023

Orthogonal pairs of rollers on concave beds (OPRCB) are a low-cost, low-tech rolling-based isolating system, whose high efficiency has been shown in a previous study. However, seismic performance of OPRCB isolators has only been studied in the two-dimensional (2D) state so far. This is while their performance in the three-dimensional (3D) state differs from that of the 2D state, mainly since the vertical accelerations due to rollers' motion in their beds, simultaneously in two orthogonal horizontal directions, are added up and resulting in bigger vertical inertia forces and higher rolling resistance. In this study, first, Lagrange equations were used to derive the governing equations of motion of the OPRCB-isolated buildings in 3D. Then, some regular shear-type OPRCB-isolated buildings were considered subjected to three-component excitations of far- and near-source earthquakes, and their responses were compared to those of their fixed-base counterparts. Finally, the effects of more realistic modeling and analysis were examined by comparing the responses of isolated buildings in 2D and 3D states. Response histories were obtained by the fourth-order Runge-Kutta-Nystrom method, considering the geometrical nonlinearity of isolators. Results reveal that utilizing the OPRCB isolators effectively reduces the acceleration response, however, depending on the system specifications and earthquake characteristics, the maximum responses of isolated buildings in the 3D state can be up to 40% higher than those in the 2D state.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.8
• /
• pp.969-977
• /
• 1999

A similarity solution of a steady laminar flow of micropolar fluids past wedges has been studied. The similarity variables found by Falkner and Skan are employed to reduce the streamwise-dependence in the coupled nonlinear boundary layer equations. Numerical solutions of the equations are then obtained using the fourth-order Runge-Kutta method and the distribution of velocity, micro-rotation, shear and couple stress across the boundary layer are obtained. These results are compared with the corresponding flow problems for Newtonian fluid past wedges with various wedge angles. Numerical results show that, keeping ${\beta}$ constant, the skin friction coefficient is lower for a micropolar fluid, as compared to a Newtonian fluid. For the case of constant material parameter K, however, the velocity distribution for a micropolar fluid is higher than that of a Newtonian fluid.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.1
• /
• pp.83-102
• /
• 2015

The present paper analyses the radiation effects of mass transfer on steady nonlinear MHD boundary layer flow of a viscous incompressible fluid over a nonlinear porous stretching surface in a porous medium in presence of heat generation. The liquid metal is assumed to be gray, emitting, and absorbing but non-scattering medium. Governing nonlinear partial differential equations are transformed to nonlinear ordinary differential equations by utilizing suitable similarity transformation. The resulting nonlinear ordinary differential equations are solved numerically using Runge-Kutta fourth order method along with shooting technique. Comparison with previously published work is obtained and good agreement is found. The effects of various governing parameters on the liquid metal fluid dimensionless velocity, dimensionless temperature, dimensionless concentration, skin-friction coefficient, Nusselt number and Sherwood number are discussed with the aid of graphs.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.3
• /
• pp.245-256
• /
• 2014

This work is focused on the boundary layer and heat transfer characteristics of hydromagnetic flow over a stretching sheet with variable thickness. Steady, two dimensional, nonlinear, laminar flow of an incompressible, viscous and electrically conducting fluid over a stretching sheet with variable thickness and power law velocity in the presence of variable magnetic field and variable temperature is considered. Governing equations of the problem are converted into ordinary differential equations utilizing similarity transformations. The resulting non-linear differential equations are solved numerically by utilizing Nachtsheim-Swigert shooting iterative scheme for satisfaction of asymptotic boundary conditions along with fourth order Runge-Kutta integration method. Numerical computations are carried out for various values of the physical parameters and the effects over the velocity and temperature are analyzed. Numerical values of dimensionless skin friction coefficient and non-dimensional rate of heat transfer are also obtained.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.2
• /
• pp.336-347
• /
• 1992

The three dimensional inviscid transonic cascade flow was investigated numerically, incorporation a four stage Runge-Kutta integration method proposed by Jameson. Time marching to the steady state was accelerated by using optimum time step and enthalpy damping. In describing the boundary conditions at inlet and outlet, Riemann invariants are considered. By adding a second and a fourth order artificial viscocities, the numerical instability due to the propagation of undamped disturbance or the rapid change of state near the shock has been prevented. The numerical results for are bump cascade, cambered two dimensional turbine cascade and three dimensional stator cascade agreed reasonably well with previous results. It has been known that the accuracy of the solution depended a lot on the modeling of the leading or trailing edge.