• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.8-15
• /
• 2012

A code developed using the flux-difference splitting scheme on the hybrid Cartesian/immersed boundary method is applied to simulate three-dimensional internal waves. The material interface is regarded as a moving contact discontinuity and is captured on the basis of mass conservation without any additional treatment across the interface. Inviscid fluxes are estimated using the flux-difference splitting scheme for incompressible fluids of different density. The hybrid Cartesian/immersed boundary method is used to enforce the boundary condition for a moving three-dimensional body. Immersed boundary nodes are identified within an instantaneous fluid domain on the basis of edges crossing a boundary. The dependent variables are reconstructed at the immersed boundary nodes along local normal lines to provide the boundary condition for a discretized flow problem. The internal waves are simulated, which are generated by an pitching ellipsoid near an material interface. The effects of density ratio and location of the ellipsoid on internal waves are compared.

• Journal of Mechanical Science and Technology
• /
• v.16 no.5
• /
• pp.683-695
• /
• 2002

The objective of the present study is to analyze the fluid flow with moving boundary using a finite element method. The algorithm uses a fractional step approach that can be used to solve low-speed flow with large density changes due to intense temperature gradients. The explicit Lax-Wendroff scheme is applied to nonlinear convective terms in the momentum equations to prevent checkerboard pressure oscillations. The ALE (Arbitrary Lagrangian Eulerian) method is adopted for moving grids. The numerical algorithm in the present study is validated for two-dimensional unsteady flow in a driven cavity and a natural convection problem. To extend the present numerical method to engine simulations, a piston-driven intake flow with moving boundary is also simulated. The density, temperature and axial velocity profiles are calculated for the three-dimensional unsteady piston-driven intake flow with density changes due to high inlet fluid temperatures using the present algorithm. The calculated results are in good agreement with other numerical and experimental ones.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.760-763
• /
• 2010

The expansion of 1D numerical code to 2D or 3D is needed in order to improve the analysis accuracy of the interior ballistics. The cut cell method has been imposed for the code expansion to multi dimensions. The MUSCL-Hancock scheme as a high resolution method has been selected. A feasibility of the cut cell method has been verified by analyzing the free piston problem.

• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.24-30
• /
• 2005

Unstructured grid system is suitable for flows of complex geometries. For problems with moving boundary walls, the grid system must be time-dependently changing and deforming according to the movement of the boundaries when we use a body fitted grid system. In this paper, a new moving-grid finite-volume method on unstructured grid system is proposed and developed for unsteady compressible flows with shock waves. To assure geometric conservation laws on moving grid system, a control volume on the space-time unified domain is adopted for estimating numerical flux. The method is described and applied for two-dimensional flows.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.659-665
• /
• 2008

The 1-D numerical study of the interior ballistics has been conducted. The unsteady compressible 1-D CFD code using SIMPLER algorithm and QUICK scheme has been developed. The mathematical model of the two-phase flow has been established for the behavior of the interior ballistics. The moving boundary due to the projectile motion as the physical phenomena of the interior ballistics results in the varied control volume. In order to analyze the moving boundary, the numerical codes, which apply the ghost-cell extrapolation method and the Lagrangian method respectively, have been developed. The ghost-cell extrapolation method has been used in the Eulerian coordinate system. The Lagrangian method has been used in Non-Eulerian coordinate system. These codes have been verified through the analysis of the free piston motion problem in the tube. Through this study, the basic techniques of the numerical code for the multi-dimensional two-phase flow of the interior ballistics have been obtained.

• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.93-100
• /
• 2009

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved successively in the moving least square sense. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as a FVM.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.165-169
• /
• 2009

The effect of the Dirichlet boundary condition for the redistance equation of level set method on the solutionof sloshing problem is investigated by adopting four Dirichlet boundary conditions. For the solution of the incompressible Navier-Stokes equations, P1P1 four-step fractional finite element method is employed and a least-square finite element method is used for the solutions of the two hyperbolic type equations of level set method; advection and redistance equation. ALE (Arbitrary Lagrangian Eulerian) method is used to deal with a moving computational domain. It has been shown that the free surface motion in a sloshing tank is strongly dependent on the type of the Dirichlet boundary condition and the results of broken dam and sloshing problems using various Dirichlet boundary conditions are discussed and compared with the existing experimental results.

• Selected Papers of The Society of Naval Architects of Korea
• /
• v.2 no.1
• /
• pp.140-155
• /
• 1994

The hydrodynamic problem is treated here when a pressurized bag is submerged partially in the water and the end points of it oscillate. SES(Surface Effect Ship) has a bag filled with pressurized air at the stern in order to prevent the air leakage, and the pitch motion of SES is largely affected by the hydrodynamic force of the bag. The shape of a bag can be determined with the pressure difference between inside and outside. Once the hydrodynamic pressure is given, the shape of a bag can be obtained, however in order to calculate the hydrodynamic pressure we should know the shape change of the bag, and vice versa. Therefore the type of boundary condition on the surface of a bag is a moving boundary like a free surface boundary. The present paper describes the formulation of this problem and treats a linearized problem. The computations of the radiation problem for an oscillating bag are shown in comparison with the case that the bag is treated as a rigid body. The hydrodynamic forces are calculated for various values of the pressure inside the bag and the submerged depth.

• Journal of the Korean Society for Precision Engineering
• /
• v.16 no.7
• /
• pp.158-167
• /
• 1999

Heat and moisture transfer associated with porous materials are investigated. The heat and moisture transfer in porous materials caused by the interaction of moisture gradient, temperature gradient, conduction, and evaporation are considered. The variations of temperature and moisture not only change the volume but also induce the hygro-thermal stress. The finite element formulation for solving the temperature and moisture transfer as well as the associated hygro-thermal stresses is developed. In order to verify the finite element formulation, the heat and moisture moving boundary problem in a half space and the hygro-thermo-mechanical problem in an infinite plate with a circular hole are analyzed. Temperature profile, moisture profile, and hygro-thermal stresses are compared with those of analytic solution and other investigator. Good agreements are examined

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.11
• /
• pp.1049-1055
• /
• 2008

Modern multidisciplinary computational fluid dynamics often incorporates moving boundaries, as would be required in the applications such as design optimization, aeroelasticity, or forced boundary motion. It is challenging to develop robust, efficient grid deformation algorithms when large displacement of the moving boundaries is required. In this paper, a volume grid deformation code is developed based on the finite macro-element and the transfinite Interpolation, and then interfaces to a structured multi-block Navier-Stokes in-house code. As demonstrated by an airfoil with pitching motion, the hysteresis loops of lift, drag and moment coefficients of the developed method are shown to be in good agreement with those of experimental data.