• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.289-294
• /
• 2010

Numerical stud of an oscillating body in incompressible fluid is performed. Stabilized finite element method comprising of Streamline-Upwind/Petrov-Galerkin (SUPG) and Pressure-Stabilizing/Petrov-Galerkin (PSPG) formulations of linear triangular elements was employed to solve 2D incompressible Navier-Stokes equations whereas the motion of the body was considered by incorporating the arbitrary Langrangian-Eulerian(ALE) formulation. An algebraic moving mesh strategy is utilized for obtaining body conforming mesh deformation at each time step. Two tests cases, namely motion of a circular cylinder and of an airfoil in incompressible flow were analyzed. The model is first validated against the stationary cases and then the capability to handle moving boundaries is demonstrated.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.9
• /
• pp.2322-2330
• /
• 1994

The torsion problem in a cylindrical rod is usually formulated in terms of either the warping function or the Prandtl stress function. In a rod whose cross-section is bounded by circles and rectangles, we develop an analytic solution approach based on the warping function, which satisfies Laplace's equation. The present formulation employs polynomials and The Fourier series-type solutions, both of which satisfy exactly the governing differential equation. Using the present method, the maximum shear stress and torsional rigidity are efficiently and accurately calculated and the present results are compared with those by other methods. The specific numerical examples include the case with eccentric holes which was investigated earlier. The finite element results are also compared with the present results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.1
• /
• pp.72-80
• /
• 2004

Steady flow of an ER (electro-rheological) fluid in a two-dimensional electrode channel is studied by using FEM. Hydrodynamic interactions between the particles and the fluid are calculated by solving the Navier-Stokes equation combined with the equation of motion for each particle, where the multi-body electrostatic interaction is described by using point-dipole model. Motion of the particles in the ER fluid is elucidated in conjunction with the mechanisms of the flow resistance and the increase of viscosity. The ER effects have been studied by varying the Mason number and volume fraction of particles. These parameters have an influence on the formation of the chains resulting in the changes of the fluid velocity and the effective viscosity of ER fluids.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.5
• /
• pp.850-858
• /
• 1997

A new a posteriori error estimator is introduced and applied to variational inequalities occurring in problems of flow through porous media. In order to construct element-wise a posteriori error estimates the global error is localized by a special mixed formulation in which continuity conditions at interfaces are treated as constraints. This approach leads to error indicators which provide rigorous upper bounds of the element errors. A discussion of a compatibility condition for the well-posedness of the local error analysis problem is given. Two numerical examples are solved to check the compatibility of the local problems and convergence of the effectivity index both in a local and a global sense with respect to local refinements.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.336-340
• /
• 2001

We propose an implicit numerical implementation for the Leblond's transformation plasticity constitutive equations, which are widely used in welded steel structure. We apply generalized trapezoidal rule to integrate the equations and determine the consistent tangent moduli. The implementation may be used with updated Lagrangian formulation. We test a simple butt-welding process to compare with SYSWELD and discuss the accuracy.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.3
• /
• pp.194-200
• /
• 2009

Suspension of a single solid particle in a channel flow with a constant pressure gradient is studied numerically. The interaction of a circular particle with a surrounding Newtonian fluid is formulated using a combined formulation. Numerical results are presented using two dimensionless variables: the sedimentation Reynolds number and the generalized Froude number. From the present results, it has been shown that a solid particle is suspended at a smaller generalized Froude number as the viscosity of the surrounding fluid increases. The time taken for equilibrium position is found to be smaller as fluid viscosity increases when both : the sedimentation Reynolds number and the generalized Froude number are the same while, at the same situation, the dimensionless time taken for equilibrium position is to be nearly the same regardless of fluid viscosity when a dimensionless time variable is introduced

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.15-18
• /
• 2004

A numerical method for simulating two-phase flows including surface force is presented. The method is based on fractional step method of finite volume formulation and the interface is tracked with PLIC VOF method. In the CSF model, as color function, f, representing the location of interface varies steeply in the interface region, we need to use smoothed function f to get accurate unit normal and the curvature. Peskin kernel is used to get smoothed function f. A spherical drop in static equilibrium and three-dimensional merging of gas bubble are tested, resulting in the validation of this method

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.215-220
• /
• 2004

A mathematical formulation based on two-phase, two-fluid hyperbolic conservation laws is developed to investigate propagation of shock waves in one- and two-dimensions. We used a high resolution upwind scheme called the split-coefficient matrix method. Two extreme cases are computed for validation of the computer code: the states of a pure gas and a pure liquid. Computed results agreed well with the previous experimental and numerical results. It is studied how the shock wave propagation pattern is affected by the void fraction in the two-phase flow. The shock structure in a two-phase flow turned out, in fact, much deviated from the shape well known in the gas only phase.

• Journal of the Korean Society for Precision Engineering
• /
• v.15 no.5
• /
• pp.120-127
• /
• 1998

Differential inequalities occurring in problems of obstacle contact problems are recast into variational inequalities and analyzed by finite element methods. A new a posteriori error estimator, which is essential in adaptive finite element method, is introduced to capture the errors in finite element approximations of these variational inequalities. In order to construct a posteriori error estimates, saddle point problems are introduced using Lagrange parameters and upper bounds are provided. The global upper bound is localized by a special mixed formulation, which leads to upper bounds of the element errors. A numerical experiment is performed on an obstacle contact problem to check the effectivity index both in a local and a global sense.

• Journal of the Korean Society for Precision Engineering
• /
• v.8 no.4
• /
• pp.95-106
• /
• 1991

A numerical computation of turbulent wirling flow in a stationary pipe is presented in this work. Major concerns of this study are: 1) To approve similarity laws which were verified experimentally. 2) To investigate the effects of curvature modification for the K- .epsilon. model. To account for effects of swirl, Rodi's curvature correction and Kim & Chung's are applied. The governing differential equations for eliptic flow are discretized by control volume formulation method, and the discretized equations are calculated ay line by line TDMA and SIMPLE algorithm. The computational results also satisfy similarity laws which are based on swirl angle as in experiments. And the curvature modification of Rodi improves compuational accuacy than the standard K- .epsilon. model. But such lower order closure models are not adequate for the prediction of this complex flow.