• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.26 no.3
• /
• pp.281-289
• /
• 2016

In this paper, a new dynamic model for modal analysis of a rotating cantilever beam with a tip-mass is developed. The nonlinear strain such as von Karman type and the corresponding linearized stress are used to consider the geometric nonlinearity, and Euler-Bernoulli beam theory is applied in the present model. The nonlinear equations of motion and the associated boundary conditions which include the inertia of the tip-mass are derived through Hamilton's principle. In order to investigate modal characteristics of the present model, the linearized equations of motion in the neighborhood of the equilibrium position are obtained by using perturbation technique to the nonlinear equations. Since the effect of the tip-mass is considered to the boundary condition of the flexible beam, weak forms are used to discretize the linearized equations. Compared with equations related to stiffening effect due to centrifugal force of the present and the previous model, the present model predicts the dynamic characteristic more precisely than the another model. As a result, the difference of natural frequencies loci between two models become larger as the rotating speed increases. In addition, we observed that the mode veering phenomenon occurs at the certain rotating speed.

• Journal of Energy Engineering
• /
• v.17 no.4
• /
• pp.233-240
• /
• 2008

Conventionally, the second-order self-adjoint neutron transport equations have been studied using the even parity and the odd parity equations. Recently, however, the SAAF(self-adjoint angular flux) form of neutron transport equation has been introduced as a new option for the second-order self-adjoint equations. In this paper we validated the SAAF equation mathematically and clarified how it relates with the existing even and odd parity equations. We also developed a second-order SAAF differencing formula including DSA(diffusion synthetic acceleration) from the first-order difference equations. Numerical result is attached to show that the proposed methods increases accuracy with effective computational effort.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.13 no.2
• /
• pp.149-157
• /
• 2001

From the weakly nonlinear mild-slope wave equations introduced by Nadaoka et al.(1994, 1997), a set of weakly nonlinear wave equations for rapidly varying topography are derived by including the bottom curvature and slope-squared tenns ignored in the original equations ofNadaoka et al. To solve the linear version of extended wave equations derived in this study one-dimensional finite difference numerical model is con¬structed. The perfonnance of the model is tested for the case of wave reflection from a plane slope with various inclination. The numerical results are compared with the results calculated using other numerical models reported earlier. The comparison shows that the accuracy of the numerical model is improved significantly in comparison with that of the original equations ofNadaoka et al. by including a complete set of bottom curva1w'e and slope¬squared terms for a rapidly varying topography.

• Steel and Composite Structures
• /
• v.18 no.6
• /
• pp.1465-1478
• /
• 2015

The effect of cutting off fibers on transient load in a polymeric matrix composite lamina was studied in this paper. The behavior of fibers was considered to be linear elastic and the matrix behavior was considered to be linear viscoelastic. To model the viscoelastic behavior of matrix, a three parameter solid model was employed. To conduct this research, finite difference method was used. The governing equations were obtained using Shear-lag theory and were solved using boundary and initial conditions before and after the development of break. Using finite difference method, the governing integro-differential equations were developed and normal stress in the fibers is obtained. Particular attention is paid the dynamic overshoot resulting when the fibers are suddenly broken. Results show that considering viscoelastic properties of matrix causes a decrease in dynamic load concentration factor and an increase in static load concentration factor. Also with increases the number of broken fibers, trend of increasing load concentration factor decreases gradually. Furthermore, the overshoot of load in fibers adjacent to the break in a polymeric matrix with high transient time is lower than a matrix with lower transient time, but the load concentration factor in the matrix with high transient time is lower.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.55-63
• /
• 1995

The flow field around a high-speed train including cross-wind effects has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations in the inertial frame using the iterative time marching scheme. The governing equations are differenced with 1st-order accurate backward difference scheme for the time derivatives, 3th-order accurate QUICK scheme for the convective terms and 2nd-order accurate central difference scheme for the viscous terms. The Marker-and-Cell concept was applied to efficiently solve continuity equation, which is differenced with 2nd-order accurate central difference scheme. The 4th-order artificial damping is added to the continuity equation for numerical stability. A C-H type of elliptic grid system is generated around a high-speed train including ground. The Baldwin-Lomax turbulent model was implemented to simulate the turbulent flows. To validate the present procedure, the flow around a high speed train at constant yaw angle of $45^{\circ}\;and\;90^{\circ}$ has been simulated. The simulation shows 3-D vortex generation in the lee corner. The flow separation is also observed around the rear of the train. It has concluded that the results of present study properly agree with physical flow phenomena.

• Nuclear Engineering and Technology
• /
• v.16 no.1
• /
• pp.21-28
• /
• 1984

The reduction of the truncation error including numerical diffusion, has been one of the most important tasks in the development of numerical methods. The stream line method is used to cancel cross numerical diffusion and some of the non-diffusion type truncation error. The two-step stream line method which is the combination of the stream line method and finite difference methods is developed in this work for the solution of the govern ing equations of incompressible buoyant turbulent flow. This method is compared with the finite difference method. The predictions of both classes of numerical methods are compared with experimental findings. Truncation error analysis also has been performed in order to the compare truncation error of the stream line method with that of finite difference methods.

• Solar Energy
• /
• v.5 no.2
• /
• pp.77-83
• /
• 1985

An analysis is performed to investigate the influence of the buoyancy force and the thickness variation of melting layer in the containment that is filled with phase-change Material surrounding a cylindrical heating tube during melting process. The phase-change material is assumed to be initially solid at its phase-change temperature and the remaining solid at any given time is still at the phase-change temperature and neglecting the effect of heat transfer occuring within the solid. At the start of melting process, the thickness of melting layer is assumed to be a stefan-problem and after the starting process, the change of temperature and velocity is calculated using a two dimensional finite difference method. The governing equations for velocity and temperature are solved by a finite difference method which used SIMPLE (Semi Implicit Method Pressure linked Equations) algorithm. Results are presented for a wide range of Granshof number and in accordance with the time increment and it is founded that two dimensional fluid flow occurred by natural convection decreases the velocity of melting process at the bottom of container. The larger the radius of heating tube, the higher heat transfer is occurred in the melting layer.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.411-420
• /
• 2009

This study presents an extended finite difference method based on moving least squares(MLS) method for solving potential problems with interfacial boundary. The approximation constructed from the MLS Taylor polynomial is modified by inserting of wedge functions for the interface modeling. Governing equations are node-wisely discretized without involving element or grid; immersion of interfacial condition into the approximation circumvents numerical difficulties owing to geometrical modeling of interface. Interface modeling introduces no additional unknowns in the system of equations but makes the system overdetermined. So, the numbers of unknowns and equations are equalized by the symmetrization of the stiffness matrix. Increase in computational effort is the trade-off for ease of interface modeling. Numerical results clearly show that the developed numerical scheme sharply describes the wedge behavior as well as jumps and efficiently and accurately solves potential problems with interface.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.21 no.5
• /
• pp.577-582
• /
• 2015

This paper is the fundamental study for the application of cold storage system to the transportation equipment by sea and land. This numerical study presents the solid-liquid phase change phenomenon of calcium chloride solution of 30wt %. The governing equations are 1-dimensional unsteady state heat transfer equations of $1^{st}$ order partial differential equations. This type of latent heat storage material is often usable in fishery vessel for controlling the temperature of container with constant condition. The governing equation was discretized with finite difference method and the program was composed with Mathcad program. The main parameters of this solution were the initial temperature of heat storage material, ambient temperature of cold air and the velocity of cold air. The data of boundary layer thickness becomes thin with the increasing of cold air flowing velocity and also the heat storage completion time become shorten.

• Journal of applied mathematics & informatics
• /
• v.39 no.3_4
• /
• pp.381-403
• /
• 2021

Consider the three-dimensional system of difference equations $x_{n+1}=\frac{{\prod_{j=0}^{k}}z_n-3j}{{\prod_{j=1}^{k}}x_n-(3j-1)\;\(a_n+b_n{\prod_{j=0}^{k}}z_n-3j\)}$, $y_{n+1}=\frac{{\prod_{j=0}^{k}}x_n-3j}{{\prod_{j=1}^{k}}y_n-(3j-1)\;\(c_n+d_n{\prod_{j=0}^{k}}x_n-3j\)}$, $z_{n+1}=\frac{{\prod_{j=0}^{k}}y_n-3j}{{\prod_{j=1}^{k}}z_n-(3j-1)\;\(e_n+f_n{\prod_{j=0}^{k}}y_n-3j\)}$, n ∈ ℕ₀, where k ∈ ℕ₀, the sequences $(a_n)_{n{\in}{\mathbb{N}}_0$, $(b_n)_{n{\in}{\mathbb{N}}_0$, $(c_n)_{n{\in}{\mathbb{N}}_0$, $(d_n)_{n{\in}{\mathbb{N}}_0$, $(e_n)_{n{\in}{\mathbb{N}}_0$, $(f_n)_{n{\in}{\mathbb{N}}_0$ and the initial values x_-3k, x_-3k+1, …, x₀, y_-3k, y_-3k+1, …, y₀, z_-3k, z_-3k+1, …, z₀ are real numbers. In this work, we give explicit formulas for the well defined solutions of the above system. Also, the forbidden set of solution of the system is found. For the constant case, a result on the existence of periodic solutions is provided and the asymptotic behavior of the solutions is investigated in detail.