• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.9-16
• /
• 2013

The heat transfer enhancement of heat sink with mist flow is studied numerically by solving the conservation equations for mass, momentum and energy in the continuous and dispersed phases. A Lagrangian method is used for tracing dispersed water droplets in the heat sink and an Eulerian species transport model for air and steam mixture. The continuous and dispersed phases are interacted with the drag and evaporation source terms. The computed results show that addition of evaporating mist droplets enhances the cooling performance of heat sink significantly.

• Journal of Korean Society for Atmospheric Environment
• /
• v.12 no.3
• /
• pp.243-254
• /
• 1996

The main purpose of this study is to investigate the collection efficiency characteristics of a cylindrical ESP. To do that, it is necessary to analyze the electric field, gas flow field, and mechanism of particle movement by numerical simulation based on EHD model. For a gas flow field, Navier-Stokes equation involving the electric source term was solved by SIMPLE algorithm. In case of the electric field, the current continuity and electric field equations were solved by S.O.R. method. The analysis of particle movement was performed on the basis of PSI-CELL model from the Lagrangian viewpoint. The results showed that the influence on the gas flow field by the electric field is almost negligible in a cylindrical ESP. The particle drift velocity $V_P$ toward the collection surface is increased continuously by the electrostatic force due to the rise of particle charge as the particle is moving to the flow direction and the particle size becomes larger. The collection efficiency is to quitely higher with the increase of applied voltage for the same particle size, while becomes smaller as the inlet velocity is increased.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.203-207
• /
• 2001

SPH is the shorthand for Smoothed Particle Hydrodynamics. This method is a Lagrangian method which means that it involves following the motion of elements of fluid. These elements have the characteristics of particles and the method is called a particle method. A useful review of SPH (Monaghan 1992) gives the basic technique and how it can be applied to numerous problems relevant to astrophysics. You can get some basic SPH programs from http) /www.maths.monash.edu.au/jjm/sphlect In the present lecture I will assume that the student has studied this review and therefore understands the basic principles. In today's lecture I plan to approach the equations from a different perspective by using a variational principle.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.369-376
• /
• 2002

The stochastic analysis of semi-infinite domain is presented using the weighted integral method, which is improved to include the higher order terms in expanding the displacement vector. To improve the weighted integral method, the Lagrangian remainder is taken into account in the expansion of the status variable with respect to the mean value of the random variables. In the resulting formulae only the 'proportionality coefficients' are introduced in the resulting equation, therefore no additional computation time and memory requirement is needed. The equations are applied in analyzing the semi-infinite domain. The results obtained by the improved weighted integral method are reasonable and are in good agreement with those of the Monte Carlo simulation. To model the semi-infinite domain, the Bettess's infinite element is adopted, where the theoretical decomposition of the strain-displacement matrix to calculate the deviatoric stiffness of the semi-infinite domains is introduced. The calculated value of mean and the covariance of the displacement are revealed to be larger than those given by the finite domain assumptions which is thought to be rational and should be considered in the design of structures on semi-infinite domains.

• Journal of Korean Society of Steel Construction
• /
• v.11 no.2 s.39
• /
• pp.153-165
• /
• 1999

In order to trace the lateral-torsional post-bucking behaviors of thin-walled space frames with non-symmetric cross sections, a geometrically non-linear finite element formulation is presented by applying incremental equilibrium equations based on the updated Lagrangian formulation and introducing Vlasov's assumption. The improved displacement field for non-symmetric thin-walled cross sections is introduced based on inclusion of second order terms of finite rotations, and the potential energy corresponding to the semitangential rotations and moments is consistently derived. For finite element analysis, tangent stiffness matrices of thin-walled space frame element are derived by using the Hermition polynomials as shape functions. A co-rotational formulation in order to evaluate the unbalanced loads is presented by separating the rigid body rotations and pure deformations from incremental displacements and evaluating the updated direction cosines and incremental member forces.

• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.35-41
• /
• 2012

This paper presents incompressible Navier-Stokes solution algorithm for 2D Free-surface flow problems on the Cartesian mesh, which was implemented to run on Graphics Processing Units(GPU). The INS solver utilizes the variable arrangement on the Cartesian mesh, Finite Volume discretization along Constrained Interpolation Profile-Conservative Semi-Lagrangian(CIP-CSL). Solution procedure of incompressible Navier-Stokes equations for free-surface flow takes considerable amount of computation time and memory space even in modern multi-core computing architecture based on Central Processing Units(CPUs). By the recent development of computer architecture technology, Graphics Processing Unit(GPU)'s scientific computing performance outperforms that of CPU's. This paper focus on the utilization of GPU's high performance computing capability, and presents an efficient solution algorithm for free surface flow simulation. The performance of the GPU implementations with double precision accuracy is compared to that of the CPU code using an representative free-surface flow problem, namely. dam-break problem.

• 제어로봇시스템학회:학술대회논문집
• /
• 2003.10a
• /
• pp.2707-2712
• /
• 2003

In this paper, with the utilization of a three-dimensional version of Leibniz’s rule, the procedure of deriving the time rate of change of an energy functional for axially moving continua is investigated. It will be shown that the method in [14], which describes the way of getting the time rate of change of an energy functional in Eulerian description, and subsequent results in [10, 11] are not complete. The key point is that the time derivatives at boundaries in the Eulerian description of axially moving continua should take into account the velocity of the moving material itself. A noble way of deriving the time rate of change of the energy functional is proposed. The correctness of the proposed method has been confirmed by other approaches. Two examples, one-dimensional axially moving string and beam equations, are provided for the purpose of demonstration. The results following the procedure proposed and the results in [14] are compared.

• The Korean Journal of Rheology
• /
• v.11 no.1
• /
• pp.18-27
• /
• 1999

Predicting the deformation behaviors of sheets in thermoforming processes has been a daunting challenge due to the strong nonlinearities arising from very large deformations, mold-polymer contact condition and hyperelasticity constitutive equations. Nonlinear numerical analysis is always required to face this challenge especially for realistic processing conditions. In this study a 3-D algorithm and the membrane approximation are developed for thermoforming processes. The constitutive equation is expressed in terms of the 2nd Piola-Kirchhoff stress tensor and the Cauchy-Green deformation tensor. The 2-term Mooney-Rivlin model is used for the material model equation. The algorithm is established by the finite element formulation employing the total Lagrangian coordinate. The deformation behavior and the stress distribution results of 3-D algorithm with various point boundary conditions are compared to those of the membrane approximation algorithm. Also, the slip boundary condition and the no-slip boundary condition are applied for the systems that have molds. Finally, the effect of sheet temperatures on the final thickness distribution is investigated for the ABS material.

• Journal of the Korea Computer Graphics Society
• /
• v.10 no.4
• /
• pp.13-21
• /
• 2004

This paper presents a physically based technique for simulating complex multiphase fluids. This work is motivated by the "stable fluids" method developed by Stam to handle gaseous fluids. We extend this technique to water, which calls for the development of methods for modeling multiphase fluids and suppressing dissipation. We construct a multiphase fluid formulation by combining the Navier-Stokes equations with the level set method. By adopting constrained interpolation profile (CIP)-based advection, we reduce the numerical dissipation and diffusion significantly. We further reduce the dissipation by converting potential1y dissipative cel1s into droplets or bubbles that undergo Lagrangian motion. Due to the multiphase formulation, the proposed method properly simulates the interaction of water with surrounding air, instead of simulating water in a void space. Moreover, the introduction of the non-dissipative technique means that, in contrast to previous methods, the simulated water does not unnecessarily lose mass and its motion is not damped to an unphysical extent. Experiments showed that the proposed method is stable and runs fast. It is demonstrated that two-dimensional simulation runs in real-time.

• Structural Engineering and Mechanics
• /
• v.26 no.6
• /
• pp.725-739
• /
• 2007

In recent years there are many plate bending elements that emerged for solving both thin and thick plates. The main features of these elements are that they are based on mix formulation interpolation with discrete collocation constraints. These elements passed the patch test for mix formulation and performed well for linear analysis of thin and thick plates. In this paper a member of this family of elements, namely, the Discrete Reissner-Mindlin (DRM) is further extended and developed to analyze both thin and thick plates with geometric nonlinearity. The Von K$\acute{a}$rm$\acute{a}$n's large displacement plate theory based on Lagrangian coordinate system is used. The Hu-Washizu variational principle is employed to formulate the stiffness matrix of the geometrically Nonlinear Discrete Reissner-Mindlin (NDRM). An iterative-incremental procedure is implemented to solve the nonlinear equations. The element is then tested for plates with simply supported and clamped edges under uniformly distributed transverse loads. The results obtained using the geometrically NDRM element is then compared with the results of available analytical solutions. It has been observed that the NDRM results agreed well with the analytical solutions results. Therefore, it is concluded that the NDRM element is both reliable and efficient in analyzing thin and thick plates with geometric non-linearity.