• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.72-78
• /
• 2014

A hybrid particle-mesh method as the combination between the Vortex-In-Cell (VIC) method and penalization method has been achieved in recent years. The VIC method, which is based on the vorticity-velocity formulation, offers particle-mesh algorithms to numerically simulate flows past a solid body. The penalization method is used to enforce boundary conditions at a body surface with a decoupling between body boundaries and computational grids. The main advantage of the hybrid particle-mesh method is an efficient implementation for solid boundaries of arbitrary complexity on Cartesian grids. However, a numerical simulation of flows in large domains is still not too easy. In this study, a multi-domain approach is thus proposed to further reduce computation cost and easily implement it. We validate the implementation by numerical simulations of an incompressible viscous flow around an impulsively started circular cylinder.

• Journal of Mechanical Science and Technology
• /
• v.19 no.spc1
• /
• pp.388-394
• /
• 2005

Recently, the authors have developed a method for real-time dynamics of multibody systems, which combines a semi-recursive formulation to derive the equations of motion in dependent relative coordinates, along with an augmented Lagrangian technique to impose the loop closure conditions. The following numerical integration procedures, which can be grouped into the so-called structural integrators, were tested : trapezoidal rule, Newmark dissipative schemes, HHT rule, and the Generalized-${\alpha}$ family. It was shown that, for large multi body systems, Newmark dissipative was the best election since, provided that the adequate parameters were chosen, excellent behavior was achieved in terms of efficiency and robustness with acceptable levels of accuracy. In the present paper, the performance of the described method in combination with another group of integrators, the Implicit Runge-Kutta family (IRK), is analyzed. The purpose is to clarify which kind of IRK algorithms can be more suitable for real-time applications, and to see whether they can be competitive with the already tested structural family of integrators. The final objective of the work is to provide some practical criteria for those interested in achieving real-time performance for large and complex multibody systems.

• Journal of Advanced Marine Engineering and Technology
• /
• v.26 no.1
• /
• pp.116-124
• /
• 2002

The bubble model by Keller and Prosperetti is adapted to solve the nonlinear oscillation of a gas bubble. This formulation leads to accurate results since it introduces the energy equation instead of the polytropic assumption for the bubble interior. The numerical method used in this study is stable enough to handle large amplitude of bubble oscillation. The numerical results show some interesting nonlinear phenomena fur the bubble oscillator. The excitation changes the natural frequency of the bubble and makes some harmonic resonances at $f/f_0=1/2, 1/3$ and so on. The natural frequency of a bubble oscillator decreases compared with the linear case result, which means that the nonlinear bubble oscillation system is a "softening"system. In addition, the frequency response curve jumps up or down at a certain frequency. It is also found that there exist multi-valued regions in the frequency response curve depending on the initial conditions of bubble. The dependency of the bubble motion on the initial condition can generate extremely large pressure and temperature which might be the cause of the acoustic cavitation and the sonoluminescence.inescence.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.7-14
• /
• 2005

In this study, the shallow water equations on the sphere is simulated by the proposed method which has high spatial resolution and is based on the CIP(Cubic Interpolated Pseudoparticle) method. The governing equations are approximated on the longitude-latitudinal coordinate system. To avoid the problems resulting from the convergence of the meridians toward high-latitude and singularities on the poles, the semi-Lagrangian and finite volume method are employed. in addition, the CIP method is employed to solve the advection equations and is extended to apply on the equations in the polar coordinate system. The mathematical formulation and numerical results are also described. To verify of the efficiency, accuracy and capability of proposed algorithm, the standard test cases proposed bv Williamson et. al are simulated and the results are compared with other results. As a result, it is found that the present scheme gives a good properties in preserving shapes of solution and settles the pole problems in solving the shallow water equations on the sphere.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1509-1517
• /
• 1993

The purpose of this study is to develop a program for the automatic derivation of the symbolic equations necessary for the kinematic and dynamic analyses of the spatial mechanism. For this purpose, a symbolic manipulation package called MCSYMA is used. Every symbolic equation is formulated using relative joint coordinate to obtain the numerically efficient system equations. These equations are produced in FORTRAN statements and linked to a FORTRAN program for numerical analysis. Several examples are taken for comparison with the commercial package called DADS which is using Cartesian coordinate approach. Also, this symbolic formulation approach is compared with a conventional numerical approach for an example. The results show that this symbolic approach with relative joint coordinate system is most efficient in computational time among three and is recommended for the derivation of macro elements frequently used.

• Korea-Australia Rheology Journal
• /
• v.13 no.1
• /
• pp.37-45
• /
• 2001

Numerical simulation of coextrusion process of viscoelastic fluids within a die has been carried out. In the coextrusion process velocity profile at the outflow boundary is not known a priori, which makes it difficult to impose the proper boundary condition at the outflow boundary. This difficulty has been avoided by using the open boundary condition (OBC) method. In this study, elastic viscous stress splitting (EVSS) formulation with streamline upwind (SU) method has been used in the finite element method. In order to test the validity of the OBC method, comparison between the results of fully developed condition at the outlet and those of OBC has been made for a Newtonian fluid. In the case of upper convected Maxwell (UCM) fluid, the effect of outflow boundary condition on the interface position has been investigated by using two meshes having different downstream lengths. In both cases, the results with the OBC method showed reasonable interface shape. In particular, for the UCM fluid the interface shape calculated with OBC was independent of the downstream length, while the results with the zero traction condition showed oscillation of interface position close to the outlet. Viscosity difference was found to be more important than elasticity difference in determining the final interface position. However, the overshoot of interface position near the con-fluent point increased with elasticity.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2003.11a
• /
• pp.465-470
• /
• 2003

In this study, the state-space Newmark method based on average velocity is presented to analyse the transient dynamic response for general dynamic system. The conventional Newmark method based on average acceleration cannot he directly to the first-order state-space differential equations introducing the state-space vector. To overcome this problem, the time-step integration algorithm, based on average velocity concept, suitable for the first-order state-space differential equations is proposed In results, the proposed method has %he numerical stability and order of accuracy, which is proved analytically, equal to those of the conventional Newmark method based on average acceleration. Also, the formulation for numerical solution is very simple and the calculation time Is nearly equal to that of the conventional Newmark method based on average acceleration in spite of an increase of two times over matrix size. This method will be look forward to applying the general dynamic system to calculate the transient dynamic response.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.2
• /
• pp.185-191
• /
• 2008

Optimal damping layout of the constrained viscoelastic damping layer on beam is identified with temperatures by using a gradient-based numerical search algorithm. An optimal design problem is defined in order to determine the constrained damping layer configuration. A finite element formulation is introduced to model the constrained layer damping beam. The four-parameter fractional derivative model and the Arrhenius shift factor are used to describe dynamic characteristics of viscoelastic material with respect to frequency and temperature. Frequency-dependent complex-valued eigenvalue problems are solved by using a simple re-substitution algorithm in order to obtain the loss factor of each mode and responses of the structure. The results of the numerical example show that the proposed method can reduce frequency responses of beam at peaks only by reconfiguring the layout of constrained damping layer within a limited weight constraint.

• International Journal of Precision Engineering and Manufacturing
• /
• v.5 no.4
• /
• pp.50-60
• /
• 2004

A very flexible beam can be used to model various types of continuous mechanical parts such as cables and wires. In this paper, the dynamic properties of a very flexible beam, included in a multibody system, are analyzed using absolute nodal coordinates formulation, which is based on finite element procedures, and the general continuum mechanics theory to represent the elastic forces. In order to consider the dynamic interaction between a continuous large deformable beam and a rigid multibody system, a combined system equations of motion is derived by adopting absolute nodal coordinates and rigid body coordinates. Using the derived system equation, a computation method for the dynamic stress during flexible multibody simulation is presented based on Euler-Bernoulli beam theory, and its reliability is verified by a commercial program NASTRAN. This method is significant in that the structural and multibody dynamics models can be unified into one numerical system. In addition, to analyze a multibody system including a very flexible beam, formulations for the sliding joint between a very deformable beam and a rigid body are derived using a non-generalized coordinate, which has no inertia or forces associated with it. In particular, a very flexible catenary cable on which a multibody system moves along its length is presented as a numerical example.

• Journal of the Korean Society for Precision Engineering
• /
• v.29 no.6
• /
• pp.632-639
• /
• 2012

A fully coupled thermo-mechanical model is adopted to study the temperature distribution and the material deformation in friction stir welding(FSW) process. Rotational speed is most important parameters in this research. Three dimension results under different process parameters were presented. Result indicate that the maximum temperature is lower than the melting point of the welding material. The higher temperature gradient occurs in the leading side of the workpiece. The maximum temperature can be increased with increasing the tool angular velocity, rpm in the current numerical modeling. In this research ABAQUS Ver.6.7 is to analyze a fully coupled thermo-mechanical model. ALE(Arbitrary Lagrangian-Eulerian) finite element formulation is used for the large deformation in FSW process and using the Mass scaling for the analysis time efficiency.