• Korea-Australia Rheology Journal
• /
• v.18 no.4
• /
• pp.171-181
• /
• 2006

We present a short review for authors' previous work on direct numerical simulations for inertialess hard particle suspensions formulated either with a Newtonian fluid or with viscoelastic polymeric fluids to understand the microstructural evolution and the bulk material behavior. We employ two well-defined bi-periodic domain concepts such that a single cell problem with a small number of particles may represent a large number of repeated structures: one is the sliding bi-periodic frame for simple shear flow and the other is the extensional bi-periodic frame for planar elongational flow. For implicit treatment of hydrodynamic interaction between particle and fluid, we use the finite-element/fictitious-domain method similar to the distributed Lagrangian multiplier (DLM) method together with the rigid ring description. The bi-periodic boundary conditions can be effectively incorportated as constraint equations and implemented by Lagrangian multipliers. The bulk stress can be evaluated by simple boundary integrals of stresslets on the particle boundary in such formulations. Some 2-D example results are presented to show effects of the solid fraction and the particle configuration on the shear and elongational viscosity along with the micro-structural evolution for both particles and fluid. Effects of the fluid elasticity has been also presented.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.299-306
• /
• 2006

Using a level set method and topological derivatives, a topological shape optimization method that is independent of an initial design is developed for linearly elastic structures. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. The 'Hamilton-Jacobi (H-J)' equation and computationally robust numerical technique of 'up-wind scheme' lead the initial implicit boundary to an optimal one according to the normal velocity field while minimizing the objective function of compliance and satisfying the constraint of allowable volume. Based on the asymptotic regularization concept, the topological derivative is considered as the limit of shape derivative as the radius of hole approaches to zero. The required velocity field to update the H -J equation is determined from the descent direction of Lagrangian derived from optimality conditions. It turns out that the initial holes is not required to get the optimal result since the developed method can create holes whenever and wherever necessary using indicators obtained from the topological derivatives. It is demonstrated that the proper choice of control parameters for nucleation is crucial for efficient optimization process.

• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.361-379
• /
• 2012

This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.22 no.1
• /
• pp.111-120
• /
• 1998

This paper is the first one of two-parted research efforts focusing on the modeling of rubber pad forming process. The rubber pad, driven by the pressurized fluid during the forming process, pushes the sheet metal to solid tool half and forms a part to final shape. In this part of the paper, a numerical procedure for the FE analysis of the rubber pad deformation is presented. The developed three-dimensional FE model is based on the total Lagrangian description of rubber maerial characterized by nearly incompressible hyper-elastic behavior under a large deformation assumption. Validity of the model as well as effects of different algorithms corresponding to incompresibility constraints and time integration methods on numerical solution responses are also demonstrated.

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.8-13
• /
• 2008

We present a numerical simulation technique and some preliminary results of the impact and spreading of a droplet containing particles on the solid substrate in 2D. We used the 2nd-order Adams-Bashforth / Crank-Nicholson method to solve the Navier-Stokes equation and employed the level-set method with the continuous surface stress for description of droplet spreading with interfacial tension. The impact velocity has been generated by the instantaneous gravity. The distributed Lagrangian-multipliers method has been combined for the implicit treatment of rigid particles and the discontinuous Galerkin method has been used for the stabilization of the interface advection equation. We investigated the droplet spreading by the inertial force and discussed effects of the presence of particles on the spreading behavior using an example problem. We observed reduced oscillation and spread for the particulate droplet.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.05a
• /
• pp.662-667
• /
• 2005

Nonlinear dynamics of active piezolaminated plates are investigated with respect to the thermopiezoelastic behaviors. For largely deformed structures with small strain, the incremental total Lagrangian formulation is presented based on the virtual work principles. A multi field layer wise finite shell element is proposed for assuring high accuracy and non-linearity of displacement, electric and thermal fields. For dynamic consideration of thermopiezoelastic snap through phenomena, the implicit Newmark's scheme with the Newton-Raphson iteration is implemented for the transient response of various piezolaminated models with symmetric or eccentric active layers. The bifurcate thermal buckling of symmetric structural models is first investigated and the characteristics of piezoelectric active responses are studied for finding snap through piezoelectric potentials and the load path tracking map. The thermoelastic stable and unstable postbuckling, thermopiezoelastic snap through phenomena with several attractors are proved using the nonlinear time responses for various initial conditions and damping loss factors. Present results show that thermopiezoelastic snap through phenomena can result in the difficulty of buckling and postbuckling control of intelligent structures.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1992.03a
• /
• pp.201-212
• /
• 1992

The plane strain analysis for simulating the stretch/draw forming operation of arbitrarily-shaped tool profiles and arbitrarily draw-in conditions is introduced. An implicit, incremental, updated Lagrangian formulation is employed, introducing a rigid-viscoplastic constitutive equation. Contact and friction are considered through the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshe without depending on the explicit spatial derivatives of tool surfaces. The FEM formulation is tested in the sections automotive inner panel and two-side draw-in. Not only the excellent agreement between measured and computed strains in the stretched section is obtained, but also the numerical stability of current formulation is verified in the two-side draw-in section.

• 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.

• Coupled systems mechanics
• /
• v.7 no.4
• /
• pp.421-440
• /
• 2018

In this study, a computational method for wall shear stress combined with an implicit direct-forcing immersed boundary method is presented. Near the immersed boundaries, the sub-grid stress is determined by a wall model in which the wall shear stress is directly calculated from the Lagrangian force on the immersed boundary. A coupling mathematical model of the transition process for a model Francis turbine comprising turbulent flow and rotating rigid guide vanes is established. The spatiotemporal distributions of pressure, velocity, vorticity and turbulent quantity are gained with the transient process; the drag and lift coefficients as well as other forces (moments) are also obtained as functions of the attack angle. At the same time, analysis is conducted of the characteristics of pressure pulsation, velocity stripes and vortex structure at some key parts of flowing passage. The coupling relations among the turbulent flow, the dynamical force (moment) response of blade and the rotating of guide vane are also obtained.

• Journal of the Korean Society for Precision Engineering
• /
• v.13 no.11
• /
• pp.143-151
• /
• 1996

Structural reaction injection molding을 이용한 복합재료의 제품 셩형을 위한 모델링 전략을 설명하였다. 사용된 모델은 두 평행한 원판형 mold에 있는 불 균일한 온도조건의 fiber preform을 통과하는 reactive srsin의 방사형 유동을 시뮬레이션 한다. 이러한 모델은 중요한 작동인자와 공정인자(주입온도, mold의 온도, 유량, cavity의 두께와 섬유의 조밀도)등이 유동속도, 변화(monomer, radical, inhibitor) 및 온도분포 등에 미치는 영향을 예견한다. 열전달과 질량전달 및 화학반응을 고려하여 모델을 개발하였다. 중요한 공정인자를 평가하기 위한 효울적인 공정창( process window)을 제공하는데 본 연구의 목적을 두었다. 2차원적인 Lagrangian 방식에 1차원적인 유동과 제한적인 2차원 열전달을 가정하여 모델을 유도하였고, 방정식은 implicit difference scheme에 의해 전개하였다. 이 모델은 Gonzalez-Romero의 실험 결과와 비교함에 의해 확인되었고, 실험결과가 잘 일치함을 보였다.