• Journal of computational fluids engineering
• /
• v.16 no.1
• /
• pp.73-82
• /
• 2011

In this paper, the vortex-induced vibration of circular cylinders is studied using the immersed boundary method on the Cartesian mesh. The Reynolds numbers considered is from 100 to 200. Using the configuration of tendemly arranged multiple circular cylinders, the vortex shedding behind of the cylinders and their flow-induced motion are investigated. The staggered MAC grid arrangement, which is the typical grid system for the incompressible flow on the Cartesian meshes, is utilized. Pressure correction method is applied for solving the divergence-free incompressible velocity field. The body motion is described by immersed boundary technique that has advantages for moving object on the fixed computational domain. It is also discussed for the computational noise in hydrodynamic forces when body motion is represented by the immersed boundary method. The Predictor/Corrector method is used for simulating the nonlinear response of the elastically mounted cylinder excited by vortex-shedding.

• International Journal of Precision Engineering and Manufacturing
• /
• v.3 no.4
• /
• pp.54-63
• /
• 2002

In recent years, mechanical systems such as high speed vehicles and railway trains moving on elastic beam structures have become a very important issue to consider. In this paper, a general approach, which can predict the dynamic behavior of a constrained mechanical system moving on a flexible beam structure, is proposed. Various supporting conditions for the foundation support are considered for the elastic beam structure. The elastic structure is assumed to be a non-uniform and linear Bernoulli-Euler beam with a proportional damping effect. Combined differential-algebraic equation of motion is derived using the multi-body dynamics theory and the finite element method. The proposed equations of motion can be solved numerically using the generalized coordinate partitioning method and predictor-corrector algorithm, which is an implicit multi-step integration method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.11 s.242
• /
• pp.1455-1462
• /
• 2005

In this paper, we intend to introduce a nonlinear finite element method based on the fully nonlinear potential flow theory in order to simulate the large amplitude sloshing flow in two-dimensional baffled tank subject to horizontally forced excitation. The free surface is tracked by a direct time differentiation scheme with the four-step predictor-corrector time integration method. The flow velocity is accurately recovered from the velocity potential by second-order least square method. In order to maintain the finite element mesh regularity and total mass, the semi-Lagrangian surface tracking method with area conservation is applied. According to the numerical formulae, we perform the parametric experiments by varying the installation height and the opening width of baffles, in order to examine the effects of baffle on the nonlinear liquid sloshing. From the numerical results, the hydrodynamic characteristics of the large amplitude sloshing are investigated.

• Structural Engineering and Mechanics
• /
• v.45 no.4
• /
• pp.521-542
• /
• 2013

In this study a truss model is used for the geometrically nonlinear static and dynamic analysis of a thin shallow arch subject to snap-through. Thanks to the very simple geometry of a truss, the equilibrium conditions can be easily written and the global stiffness matrix can be easily updated with respect to the deformed structure, within each step of the analysis. A very coarse discretization is applied; so, in a very simple way, the high frequency modes are suppressed from the beginning and there is no need to develop a complicated reduced-order technique. Two short computer programs have been developed for the geometrically nonlinear static analysis by displacement control of a plane truss model of a structure as well as for its dynamic analysis by the step-by-step time integration algorithm of trapezoidal rule, combined with a predictor-corrector technique. These two short, fully documented computer programs are applied on the geometrically nonlinear static and dynamic analysis of a specific thin shallow arch subject to snap-through.

• 한국연소학회:학술대회논문집
• /
• 2002.06a
• /
• pp.68-75
• /
• 2002

The dynamic behaviors of the single vortex and flame-vortex interaction in a $CH_4/Air$ diffusion flame with addition of $CO_2$ were investigated numerically. The numerical method was based on a predictor-corrector for low Mach number flow. A two-step global reaction mechanism was adopted as a combustion model. Through comparison of results by effect of $CO_2$ added either on the fuel or oxidizer side, it was found that the growth of single vortex and entrainment of surrounding fluid by $CO_2$addition on the fuel side are larger than those by $CO_2$ addition on oxidizer side. Also, when $CO_2$ is added on fuel side, flame-vortex interaction becomes more significant than on air side.

• Structural Engineering and Mechanics
• /
• v.5 no.5
• /
• pp.491-505
• /
• 1997

This paper deals with the formulation and solution of a wide class of structures, in the presence of both geometric and material nonlinearities, as a particular mathematical programming problem. We first present key ideas for the nonholonomic (path dependent) rate formulation for a suitably discretized structural model before we develop its computationally advantageous stepwise holonomic (path independent) counterpart. A feature of the final mathematical programming problem, known as a nonlinear complementarity problem, is that the governing relations exhibit symmetry as a result of the introduction of so-called nonlinear "residuals". One advantage of this form is that it facilitates application of a particular iterative algorithm, in essence a predictor-corrector method, for the solution process. As an illustrative example, we specifically consider the simplest case of plane trusses and detail in particular the general methodology for establishing the static-kinematic relations in a dual format. Extension to other skeletal structures is conceptually transparent. Some numerical examples are presented to illustrate applicability of the procedure.

• Proceedings of the Korean Operations and Management Science Society Conference
• /
• 1996.04a
• /
• pp.530-533
• /
• 1996

This paper deals with finding an initial solution and modifying search direction by the centrering force in the predictor-corrector method which is a variant of the primal-dual barrier method. These methods were tested with NETLIB problems. Initial solutions which are located close to the center of the feasible set lower the number of iterations, as they enlarge the step length. Three heuristic methods to find such initial solution are suggested. The new methods reduce the average number of iterations by 52% to at most, compared with the old method assigning 1 to initial valurs. Solutions can move closer to the central path fast by enlarging the centering force in early steps. It enlarge the step length, so reduces the number of iterations. The more effective this method is the closer the initial solution is to the boundary of the feasible set.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.11
• /
• pp.2804-2811
• /
• 2000

A Nordsick form opf the multirate integration scheme has been proposed for flexible multibody dynamic systems. It is assumed that vibrational modal coordinates in the equations of motion are treated as fast variables, whereas the relative joint coordinates are treated as slow variables. In the multirate integration, the fast variables are integrated with small step-size, and the slow variables are integrated with larger step-size. The proposed multirate integration method is based on the Adams-Bashforth-Moulton predictor-corrector method and implemented in the Nordsieck vector form. The Nordsieck form of multrate integration method provides effective step-size control and at the same time, inherits the efficiency from the Adams integration method. Simulations of a flexible gun and turret system of the military tank have been carried out to show the effectiveness and efficiency of the proposed method.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.551-556
• /
• 2001

The filling pattern technique based on the finite element method and Eulerian mesh advancement approach has been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The penalty and predictor-corrector methods are used effectively for finite element formulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. Using the proposed numerical technique, the collapse of a dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions versus time have been compared with the reported experimental result.

• Interaction and multiscale mechanics
• /
• v.2 no.3
• /
• pp.263-276
• /
• 2009

This paper is intended to investigate interaction response of a train running over a suspension bridge undergoing support settlements. The suspension bridge is modeled as a single-span suspended beam with hinged ends and the train as successive moving oscillators with identical properties. To conduct this dynamic problem with non-homogeneous boundary conditions, this study first divides the total response of the suspended beam into two parts: the static and dynamic responses. Then, the coupled equations of motion for the suspended beam carrying multiple moving oscillators are transformed into a set of nonlinearly coupled generalized equations by Galerkin's method, and solved using the Newmark method with an incremental-iterative procedure including the three phases: predictor, corrector, and equilibrium-checking. Numerical investigations demonstrate that the present iterative technique is available in dealing with the dynamic interaction problem of vehicle/bridge coupling system and that the differential movements of bridge supports will significantly affect the dynamic response of the running vehicles but insignificant influence on the bridge response.