• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.11
• /
• pp.895-903
• /
• 2007

CMBT(Curved Moving Boundary Treatment) is a newly developed scheme for the treatment of a no slip condition on the curved solid wall of moving obstacle in a flow field. In our research CMBT was used to perform LBM simulation of a flow over a moving circular cylinder to determine the flow feature and aerodynamics characteristic of the cylinder. To ascertain the applicability of CMBT on the complex shape of the obstacle, it was first simulated for the case of the flow over a fixed circular cylinder in a channel and the results were compared against the solution of Navier-Stokes equation with deforming mesh technique. The simulations were performed in a moderate range of reynolds number at each moving cylinder to identify the flow feature and aerodynamic characteristics of circular cylinder in a channel. The drag coefficients of the cylinder were calculated from the simulation results. We have numerically confirmed that the critical reynolds number for vortex shedding is ar Re=250 and the result is the same as the case of fixed cylinder. As the cylinder approaching to one wall, the 2nd vortex is developed by interacting with the wall boundary-layer vorticity. As the velocity ratio increase the third vortex are generated by interacting with the 2nd vortexes developed on the upper and lower wall boundary layer. The resultant $C_d$ decrease as reynolds number increasing and the Cd approached to a value when Re>1000.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2001.04a
• /
• pp.77-80
• /
• 2001

The element-free Galerkin (EFG) method is one of meshless methods, which is an efficient method of modeling problems of fluid or solid mechanics with complex boundary shapes and large changes in boundary conditions. This paper discusses the theory of the EFG method and its applications to modeling of groundwater flow. In the EFG method, shape functions are constructed based on the moving least square (MLS) approximation, which requires only set of nodes. The EFG method can eliminate time-consuming mesh generation procedure with irregular shaped boundaries because it does not require any elements. The coupled EFG-FEM technique was introduced to treat Dirichlet boundary conditions. A computer code EFGG was developed and tested for the problems of steady-state and transient groundwater flow in homogeneous or heterogeneous aquifers. The accuracy of solutions by the EFG method was similar to that by the FEM. The EFG method has the advantages in convenient node generation and flexible boundary condition implementation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.10
• /
• pp.1337-1345
• /
• 2001

Numerical investigation is made to study the effects of thermosolutal convection on the formation of macrosegregation in a Pb-Sn alloy solidification process in a two dimensional confined rectangluar mold. The basic equations are sovled using the Contrinum Model theory with the SIMPE algorithm during the solidification process. In addition, to track the liquid-solid interface with time variations, the moving boundary condition was adopted and moving irregular interface shapes were treated with the time-dependent, boundary-fitted coordinate system. As the temperature reduces from the liquidus to the solidus, the liquid concentration of Sn, the lighter constituent, increases. Then the buoyancy-driven flow due to temperature and liquid composition gradients occurs in the mushy region and forms the complicated macrosegregation maps. belated to this phenomena, effects on the macrosegregation formation depending on the cooling condition and gravity values are examined.

• Journal of Ship and Ocean Technology
• /
• v.10 no.3
• /
• pp.1-16
• /
• 2006

Despite the high cost of memory and CPU time required to resolve the boundary layer, a viscous unstructured grid solver has many advantages over a structured grid solver such as the convenience in automated grid generation and vortex capturing by solution adaption. In present study, an unstructured Cartesian grid solver is developed on the basis of the existing Euler solver, NASCART-GT. Instead of cut-cell approach, immersed boundary approach is applied with ghost cell boundary condition, which can be easily applied to a moving grid solver. The standard $k-{\varepsilon}$ model by Launder and Spalding is employed for the turbulence modeling, and a new wall function approach is devised for the unstructured Cartesian grid solver. Developed approach is validated and the efficiency of the developed boundary condition is tested in 2-D flow field around a flat plate, NACA0012 airfoil, and axisymmetric hemispheroid.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.305-310
• /
• 2001

This paper presents the computational method for analyzing the compressible flow fields in a high voltage gas circuit breaker. There are many difficult problems in analyzing the gas flow in GCB due to complex geometry, moving boundary, shock wave and so on. In particular, the distortion problem of the grid due to the movement of moving parts can be worked out by the fixed grid technique. Numerical simulations are based on a fully implicit finite volume method of the compressible Reynolds-averaged Navier-Stokes equations to obtain the pressure, density, and velocity through the entire interruption process. The presented method is applied to the real circuit breaker model and the pressure in front of the piston is good agreement with the experimental one.

• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.16-23
• /
• 2005

The objective of this research is to estimate two dimensional ablating and charring of heat shield materials in severe aero-thermal heat transfer. This estimation requires an accurate and rapid technique for its serious heat transfer with a moving boundary. Aerodynamic heating is obtained by an explicit relation which is a function of Mach number and air condition, while a fully implicit method is used for heat transfer calculations. Moving boundary is captured by FLIAR method which is a subgroup of VOF. Thickness of ablating and charring of heat shield, temperature of the moving surface and rate of radiation heat are calculated and compared with references. The results are in good agreement with other calculations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.12
• /
• pp.1668-1678
• /
• 1997

Numerical analyses were performed to investigate the characteristics of regenerator in oscillating flow by using moving boundary method and Darcy model. In this work, periodic adiabatic boundary condition was suggested as the boundary condition of adiabatic part so that the effects of the thermal inertia of the wall could be considered. In carrying out numerical analyses, two models were applied and compared. One called isotropic model has the same thermal conductivity in radial and axial directions within a porous media. The other called aeolotropic model has different conductivity in each directions. Isotropic model could not show the advantage of energy reduction which needs to maintain constant wall temperature difference between heater and cooler. But aeolotropic model could simulate the reduction of energy consumption.

• Ocean Systems Engineering
• /
• v.6 no.3
• /
• pp.245-264
• /
• 2016

The iterative boundary element method (IBEM) developed originally before for cavitating two-dimensional (2-D) and three-dimensional (3-D) hydrofoils moving under free surface is modified and applied to the case of 2-D (two-dimensional) airfoils and 3-D (three-dimensional) wings over water. The calculation of the steady-state flow characteristics of an inviscid, incompressible fluid past 2-D airfoils and 3-D wings above free water surface is of practical importance for air-assisted marine vehicles such as some racing boats including catamarans with hydrofoils and WIG (Wing-In-Ground) effect crafts. In the present paper, the effects of free surface both on 2-D airfoils and 3-D wings moving steadily over free water surface are investigated in detail. The iterative numerical method (IBEM) based on the Green's theorem allows separating the airfoil or wing problems and the free surface problem. Both the 2-D airfoil surface (or 3-D wing surface) and the free surface are modeled with constant strength dipole and constant strength source panels. While the kinematic boundary condition is applied on the airfoil surface or on the wing surface, the linearized kinematic-dynamic combined condition is applied on the free surface. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream boundary in 2-D airfoil and 3-D wing cases and transverse boundaries in only 3-D wing case. The method is first applied to 2-D NACA0004 airfoil with angle of attack of four degrees to validate the method. The effects of height of 2-D airfoil from free surface and Froude number on lift and drag coefficients are investigated. The method is also applied to NACA0015 airfoil for another validation with experiments in case of ground effect. The lift coefficient with different clearance values are compared with those of experiments. The numerical method is then applied to NACA0012 airfoil with the angle of attack of five degrees and the effects of Froude number and clearance on the lift and drag coefficients are discussed. The method is lastly applied to a rectangular 3-D wing and the effects of Froude number on wing performance have been investigated. The numerical results for wing moving under free surface have also been compared with those of the same wing moving above free surface. It has been found that the free surface can affect the wing performance significantly.

• International Journal of Highway Engineering
• /
• v.10 no.4
• /
• pp.213-224
• /
• 2008

Flexible pavement responses to vehicular loading, such as critical stresses and strains, in each pavement layer, could be predicted by the multilayered elastic analysis. However, multilayered elastic theory suffers from major drawbacks including spatial dimension of a numerical model, material properties considered in the analysis, boundary conditions, and ill-presentation of tire-pavement contact shape and stresses. To overcome these shortcomings, three-dimensional finite element (3D FE) models are developed and numerical analyses are conducted to calculate pavement responses to moving load in this study. This paper introduces a methodology for an effective 3D FE to simulate flexible pavement structure. It also discusses the mesh development and boundary condition analysis. Sensitivity analyses of flexible pavement response to loading are conducted. The infinite boundary conditions and time-dependent history of calculated pavement responses are considered in the analysis. This study found that the outcome of 3D FE implicit dynamic analysis of flexible pavement that utilizes appropriate boundary conditions, continuous moving load, viscoelastic hot-mix asphalt model is comparable to field measurements.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.37-44
• /
• 2000

Numerical investigation is made to study the effects of natural convection on the formation of macrosegregation of a Pb-Sn alloy solidification process in a 2-D confined rectangle mold. The governing equations are calculated using previous continuum models with SIMPLE algorithm doring the solidification process. In addition. to track the solid-liquid interface with time variations. the moving boundary condition Is adopted and irregular interface shapes are treated with Boundary-Fitted Coordinate system. As the temperature reduce from the liquidus to the solidus, the liquid concentration of Sn. the lighter constituent, increases. Then the buoyancy-driven flow due to temperature and liquid composition gradients, called thermosolutal convection or double diffusion, occurs in the mushy region and forms the complicated macrosegregation maps. Related to this phnomena, effects on the macrosegregation formation depending on the cooling condition and gravity values are described.