• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.54-57
• /
• 2011

In this paper we present an algorithm about how to simulate two dimensional dam breaking with lattice Boltzmann method (LBM). LBM considers a typical volume element of fluid to be composed of a collection of particles that represented by a particle velocity distribution function for each fluid component at each grid point. We use the modified Lattice Boltzmann Method for incompressible fluid. This paper will represent detailed information on single phase flow which considers only the water instead of both air and water. Interface treatment and conservation of mass are the most important things in simulating free surface where the Interface is treated by mass exchange with the water region. We consider the surface tension on the interface and also bounce back boundary condition for the treatment of solid obstacles. We will compare the results of the simulation with some methods and experimental results.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.408-409
• /
• 2008

This paper presents a new development for topology optimization of heat-dissipating structure with forced convection. To cool down electric devices or machines, two types of convection models have been widely used: Natural convection model with a large Archimedes number and Forced convection with a small Archimedes number. Nowadays, many engineering application areas such as electrochemical conversion device or fuel cell devices adopt the forced convection to transfer generated heat. Therefore, to our knowledge, it becomes an important issue to design flow channels inside which generated heat transfer. Thus, this paper studies optimal topological designs considering fluid-heat interaction. To consider the effect of the advection in the heat transfer problem, the incompressible Navier-stokes equation is solved. This paper numerically studies the coupling phenomena and presents optimal channel design considering forced convection.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.156-161
• /
• 1998

This describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows of fluid behaviour with free-surface. The elliptic differential equations governing the flows have been linearized by means of finite-difference approximations, and the resulting equations have been solved via a fully-implicit iterative method. The free-surface is defined by the motion of a set of marker particles and interface behaviour was investigated by way of a 'Lagrangian' technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions or experimental results from the literature. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.2
• /
• pp.616-623
• /
• 1996

This paper investigates the linear stability of reacting mixing layers with special emphasis on the existence of multiple unstable modes. The governing equations for laminar flows are from two-dimensional compressible boundary-layer equations. The chemistry is a finite rate single step irreversible reaction with Arrhenius kinetics. For the incompressible reacintg mixing layer with variable density. A necessary condition for instability has been derived. The condition requires that the angular momentum, not the vorticity, to have a maximum in the flow domain. New inflectional modes of instability are found to exist in the outer part of the mixing layer. For the compressible reacting mixing layer, supersonic unstable modes may exist in the abscence of a generalized inflection point. The outer modes at high Mach numbers in the reacting mixing layer are continuations of the inflectional modes of low Mach number flows. However, the generalized inflection point is less important at supersonic flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.2
• /
• pp.741-751
• /
• 1996

The coupled conduction and convection heat transfer from the protruding heat source in a vertical channel is numerically investigated. Conjugate solution of the two-dimensional energy equation is obtained for the incompressible air flow over the rectangular block with local heat source. It was found that several recirculation zones and separation bubble near the block were related to Re and Gr. And the results show that fractions of the heat transfer through each of the block face, maximum temperature of the block and the relative effect of each parameter on the maximum temperature and heat transfer.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.798-803
• /
• 2003

The objective of the present study is to investigate the changes in the acoustic source characteristics and far-field noise propagation in an incompressible round jet at Re=10000 for single-frequency excitations using large eddy simulation and Lighthill acoustic analogy. We apply excitations at a frequency corresponding to the jet-column mode ($St_{D}=0.85$) or maximum growth rate in the shear layer ( $St_{\theta}=0.017$ ). The acoustic source derived from the Lighthill acoustic analogy is the second spatial derivative of the Reynolds stresses. In the case of $St_{D}=0.85$, vortex ring and large scale structures are dominant sources, whereas in the case of $St_{\theta}=0.017$, the main sources are located at an upstream position along the shear layer than in the uncontrolled case. Also, the far-field noise propagates along the axial direction due to excitation.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.04a
• /
• pp.636-642
• /
• 1997

Laminar vortex shedding behind a circular cylinder with and without splitter plates attached to the circular cylinder at low Reynolds numbers are simulated by solving the unsteady incompressible Navier-Strokes equations. The Strouhal number, lift and drag rapidly change with the splitter plate. Far-field noise from the vortex shedding behind the cylinder is computed using the Lighthill acoustic analogy and the Curle's solution for the Lighthill equation. The acoustic source functions are obtained from the computed near-field velocity and pressure. Numerical results show that the volume quadrupole noise is small at low Mach numbers, compared with the surface dipole noise. Also the amplitude and frequency of the acoustic density fluctuations are varied with the length of splitter plates. The scattering effects at the edge of a splitter plate are considered by using the half-plane Green's function.

• The KSFM Journal of Fluid Machinery
• /
• v.7 no.1 s.22
• /
• pp.9-16
• /
• 2004

A three-dimensional numerical analysis was carried out on incompressible flows in butterfly valves by using commercial FLUENT/UNS code. Characteristics of complex flows including cavitation effect were investigated for different valve disk angles. The butterfly-valves that had different disk angles and different disk shapes were compared with each other in detail. This study focused on the flow analysis in the conventional butterfly valve and the newly designed butterfly valve in order to prevent cavitation. The newly designed valve shows great improvement on performance and endurance.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.60-67
• /
• 1998

This paper deals with the development of a variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined elements and efficiently used for construction of a refined mesh without generating distorted elements. A modified Gaussian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1994.10a
• /
• pp.55-63
• /
• 1994

The analysis of superplastic sheet forming process is studied by the use of the finite element method using a convected coordinate system and a skew boundary condition. In the formulation, the large inelastic behavior of the superplastic material is described as incompressible, nonlinear, viscous flow. The formulation is then approximated to the finite dimensional space with the use of membrane elements, which results in algebraic linear equations. In addition to the finite element formulation, a pressure cycle control algorithm is combined in the analysis for optimization of the forming time, which deals with the maximization of the strain rate sensitivity, the protection of the thickness reduction, the consistency of the desired strain rate and improvement of formability.