• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.4 s.247
• /
• pp.320-327
• /
• 2006

The computations of chemically reacting laminar and turbulent flows are performed using the preconditioned Navier-Stokes solver coupled with turbulent transport and multi-species equations. A low-Reynolds number $k-\varepsilon$ turbulence model proposed by Chien is used. The presence of the turbulent kinetic energy tenn in the momentum equation can materially affect the overall stability of the fluids-turbulence system. Because of this coupling effect, a fully coupled formulation is desirable and this approach is taken in the present study. Choi and Merkle's preconditioning technique is used to overcome the convergence difficulties occurred at low speed flows. The numerical scheme used for the present study is based on the implicit upwind ADI algorithm and is validated through the comparisons of computational and experimental results for laminar methane-air diffusion flame and $ H_2/O_2$ reacting turbulent shear flow. Preconditioning formulation shows better convergence characteristics than that of non-preconditioned system by approximately five times as much.

• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.13-21
• /
• 2013

Cavitating flow is usually formed on the surface of a high speed underwater object. When a object moves near a free surface at very high speed, the cavity signature becomes one of the major factors to be overcome by sensors of military satellite. The present work was to study the free surface effect on the ventilated cavitation. The governing equations were Navier-Stokes equations based on a homogeneous mixture model. The multiphase flow solver used an implicit preconditioning method in the curvilinear coordinate system. The cavitation model used here was the one first presented by Merkle et al.(2006) and redeveloped by Park & Ha(2009). Computations considered the free surface effects were carried out with a NACA0012 hydrofoil and the corresponding results were compared with the experimental data to have a good agreement. Calculations were then performed considering the ventilated cavitation, including the effect of non-condensable gas under the free surface effects.

• International Journal of Aeronautical and Space Sciences
• /
• v.12 no.4
• /
• pp.318-330
• /
• 2011

This paper investigates the convergence characteristics of the modified artificial compressibility method proposed by Turkel. In particular, a focus is mode on the convergence characteristics due to variation of the preconditioning factor (${\alpha}_u$) and the artificial compressibility (${\beta}$) in conjunction with an upwind method. For the investigations, a code using the modified artificial compressibility is developed. The code solves the axisymmetric incompressible Reynolds averaged Navier-Stokes equations. The cell-centered finite volume method is used in conjunction with Roe's approximate Riemann solver for the inviscid flux, and the central difference discretization is used for the viscous flux. Time marching is accomplished by the approximated factorization-alternate direction implicit method. In addition, Menter's k-${\omega}$ shear stress transport turbulence model is adopted for analysis of turbulent flows. Inviscid, laminar, and turbulent flows are solved to investigate the accuracy of solutions and convergence behavior in the modified artificial compressibility method. The possible reason for loss of robustness of the modified artificial compressibility method with ${\alpha}_u$ >1.0 is given.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.2 s.32
• /
• pp.138-146
• /
• 1999

The viscous flow around a ship hull is calculated by the use of RANS(Reynolds-averaged Navier-Stokes) solver. Reynolds stresses are midelled by using the k-${epsilon}$ turbulence model and the law is applied near the body. Body fitted corrdinates are introduced for the treatment of the complex boundary of the ship hull form and the governing equations in the physical domain transformed into ones in the computational domain. The transformed equations are numerically solved by an employment of FVM(Finite Volume Method). SIMPLE(Semi-Implicit Pressure Linked Equation) method is adopted in the calculation of pressure and the solution of the sidcretized equation is obtained by the line-by-line method with the use of TDMA(Tri-Diagonal Matrix Algorithme). To assure the proprietty of this computing method, HSVA tanker and Dyne hull are calculated ar both model and ship scale Reynolds number. Their reaults of pressure distributions on fore and aft body, axial velocity contours and transverse velocity velocity vectors and viscous resistance coefficients are compared with other's experiments and calculations.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2540-2543
• /
• 2008

Liquid crystal displays (LCD's) are continuously coated with some chemicals in the clean room of a factory. Spreading of these chemicals is causing serious problems both in controlling clean room quality as well as to the workers inside the factory. It is required to alleviate or properly control the offensive odor which is mainly composed of propylene glycol mono ethyl acetate, novolak resin and photo active compound. The control strategy employed is to bleed the offensive odor gas out the clean room. A full scale 3D CFD model was created with anisotropic porous media, chemical species transport with no volumetric reaction, and thermal diffusion with propane gas (tracer gas) to simulate the odor spreading. A segregated implicit solver with standard k-$\varepsilon$ model is employed. The detailed CFD analysis made it possible to develop an effective method of ventilating the coater room and optimizing their capacities.

• Journal of the Korean GEO-environmental Society
• /
• v.18 no.2
• /
• pp.53-58
• /
• 2017

In this paper, we propose a new method to debond between mixed finite elements for porous media in ABAQUS (2014). ABAQUS just provides debonding algorithm for the u-p model using cohesive elements in standard version. However, this approach has a drawback that it is hard to simulate complex debonding problems like element separation, rigid body motion, and contact between separated elements in standard version. ABAQUS-explicit can resolve these complex problems, but cohesive elements for the u-p model cannot be applied. We introduce a new algorithm for debonding for porous media instead of using cohesive elements. In this method, subroutines VUMAT to apply constitutive models and VDISP to separate elements in ABAQUS are used to simulate debonding problems. In addition, a simple 2-D example is demonstrated in the ABAQUS-explicit solver.

• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• Composites Research
• /
• v.26 no.6
• /
• pp.363-372
• /
• 2013

Numerical simulation for Resin Transfer Molding manufacturing process is attempted by using the eXtended Finite Element Method (XFEM) combined with the level set method. XFEM allows to obtaining a good numerical precision of the pressure near the resin flow front, where its gradient is discontinuous. The enriched shape functions of XFEM are derived by using the level set values so as to correctly describe the interpolation with the resin flow front. In addition, the level set method is used to transport the resin flow front at each time step during the mold filling. The level set values are calculated by an implicit characteristic Galerkin FEM. The multi-frontal solver of IPSAP is adopted to solve the system. This work is validated by comparing the obtained results with analytic solutions. Moreover, a localization method of XFEM and level set method is proposed to increase the computing efficiency. The computation domain is reduced to the small region near the resin flow front. Therefore, the total computing time is strongly reduced by it. The efficiency test is made with a simple channel flow model. Several application examples are analyzed to demonstrate ability of this method.