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

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.1
• /
• pp.1-8
• /
• 2001

The computations of the turbulent flow around the ship models with the free-surface effects were carried out. Incompressible Reynolds-Averaged Navier-Stokes equations were solved by using an explicit finite-difference method with the nonstaggered grid system. The method employed second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration. For the turbulence closure, a modified Baldwin-Lomax model was exploited. The location of the free surface was determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and a free-surface conforming grid was generated at each time step so that one of the grid boundary surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition was applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method developed in the present study, the computations were carried out for beth Wigley and Series 60 $C_B=0.6$ ship model and the computational results showed good agreements with the experimental data.

• Fire Science and Engineering
• /
• v.28 no.1
• /
• pp.52-57
• /
• 2014

In this paper, the unsteady incompressible Navier-Stokes equations coupled with energy equation were solved to find out the optimal location of electrical heat trace for anti-freeze of water inside the pipe for fire protection. Since the conduction equation of pipe was coupled with the natural convection of water, the analysis of conjugate heat transfer was conducted. A commercial code (ANSYS-FLUENT) based on SIMPLE-type algorithm was used for investigating the unsteady flows and temperature distributions in water region. From the numerical experiments, the isotherms and the vector fields in water region were obtained. Furthermore, it was found that the lowest part of the pipe cross-section was an optimal position of electrical heat trace assuming the constant thermal expansion coefficient of water since the minimum temperature of the water with the position is higher than those with the other positions.

• Fire Science and Engineering
• /
• v.29 no.6
• /
• pp.33-39
• /
• 2015

In this paper, a numerical experiment was conducted to find out the optimal location of electrical heat trace for anti-freeze of water inside the CPVC pipe for fire protection. The unsteady incompressible Navier-Stokes equations coupled with energy equation were solved. Since the conduction equation of pipe was coupled with the natural convection of water, the analysis of conjugate heat transfer was conducted. A commercial code (ANSYS-FLUENT) based on SIMPLE-type algorithm was used for investigating the unsteady flows and temperature distributions in water region. From the present numerical experiment, it has been found that the vector field of water inside the PVC pipe is opposite to the case of steel because of the huge difference of material properties of the two pipes. Furthermore, it was found that the lowest part of the pipe was an optimal position for electrical heat trace since the minimum water temperature of the case was higher than those of the other cases.

• Fire Science and Engineering
• /
• v.30 no.1
• /
• pp.49-56
• /
• 2016

This paper numerically examines, straight and twisted electrical heat trace installations for their anti-freezing effects on water inside a pipe. The unsteady incompressible Navier-Stokes equations coupled with an energy equation were solved to compare the two installation methods. The heat conduction of the pipe with a heat source interacts with the natural convection of the water, and the conjugate heat transfer was considered using a commercial code (ANSYS-FLUENT) based on a SIMPLE-type algorithm. Numerical experiments, were done to investigate the isotherms and the vector fields in the water region to extract the evolutions of the minimum and maximum temperatures of the water inside the pipe. There was no substantial difference in the anti-freezing effects between the straight and twisted. Therefore, the straight installation is recommended after considering the damage and short circuit behavior of the electrical heat trace.

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

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.1-10
• /
• 1998

This paper presents a vorticity-based numerical method for analyzing an incompressible Newtonian viscous flow around an impulsively started cylinder. The Navier-Stockes equations have a natural Helmholtz decomposition. The vorticity transport equation and the pressure equation are derived from this decoupled form. The associated boundary conditions are dynamic for the vorticity and pressure variables representing the coupling relation between them and the force balance on the wall. The various numerical treatments for solving the governing equations are introduced. According to Wu et al.(1994), the boundary conditions are decoupled, keeping the dynamic relation between vorticity and pressure. The vorticity transport equation is formulated by FVM and TVD(Total Variation Diminishing) scheme is used for the convection term. An integral approach similar to the panel method is used to obtain the velocity field for a given vorticity field and the pressure field, instead of the conventional differential approaches. In the numerical process, the structured grid is generated. The results are compared to existing numerical and analytic results for the validity of the present method.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.60-68
• /
• 1991

The updated Lagrangian Finite Element Method is introduced to analyse rigid body-fluid impact problem which is characterized by incompressible Navier-Stokes equations and impact-contact conditions between free surface and rigid body. For the convenience of numerical computation, velocity fields are splinted into vicous and pressure parts, and then the governing equations and boundary conditions are decomposed in accordance with the decomposition. However, Viscous stresses acting an the solid boundaries are neglected on the assumption that very small velocity gradients may occur during extremely small time interval of the impact. Four coded quadrilateral elements are used to discretize the space domain and the fully explicit time-marching algorithm is employed with a reasonably small time step. At the beginning of each time step, contact velocity of the rigid body is computed from the momentum balance between the body and the fluid. The velocity field is then computed to satisfy the discretized equations of motions and incompressibility and contact constraints as well as an exact free surface boundary condition. At the end of each time step, the fluid domain is updated from the velocity field. In the present time stepping numerical analysis, behaviour of the free surface near the body can be observed without any difficulty which is very important in the water impact problem. The applicability of the algorithm is illustrated by a wedge type falling body problem. The numerical solutions for time-varying pressure distributions and impact loadings acting ion the surface are obtained.

• The Journal of the Acoustical Society of Korea
• /
• v.36 no.5
• /
• pp.314-320
• /
• 2017

A pressure relief valve is generally used to prevent piping systems from being broken due to high pressure gas flows. However, the sudden pressure drop caused by the pressure relief valve produces high acoustic energy which propagates in the form of compressible acoustic waves in the pipe and sometimes causes severe vibration of the pipe structure, thereby resulting in its failure. In this study, internal aerodynamic noise due to valve flow is estimated for a simple contraction-expansion pipe by combining the LES (Large-Eddy Simulation) technique with the wavenumber-frequency analysis, which allows the decomposition of fluctuating pressure into incompressible hydrodynamic pressure and compressible acoustic pressure. In order to increase the convergence, the steady Reynolds-Averaged Navier-Stokes equations are numerically solved. And then, for the unsteady flow analysis with high accuracy, the unsteady LES is performed with the steady result as the initial value. The wavenumber-frequency analysis is finally performed using the unsteady flow simulation results. The wavenumber-frequency analysis is shown to separate the compressible pressure fluctuation in the flow field from the incompressible one. This result can provide the accurate information for the source causing so-called acoustic-induced-vibration of a piping system.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2294-2303
• /
• 1993

This paper investigated the characteristics of the turbulent incompressible flow past the orifice ring in an axi-symmetric pipe. The flow field was the turbulent pulsatile flow for Reynolds number of $2{\times}10^{5}$ which was defined based on the maximum velocity and the pipe diameter at the inlet, with oscillating frequence $(f_{os})=1/4{\pi}$ which was considered as quasi-steady state frequence. In the present investigation, finite analytic method was used to solve the governing equations in Navier Stokes and turbulent transport formulations. Particularly at high Reynolds number and low oscillation frequency, the effects of orifice ring on the flow were numerically investigated. The separation zone behind the orifice ring during the acceleration phase was found to be decreased. However, during the deceleration phase, the separation behind the orifice ring for pulsatile flow continuously grow to a size even larger than that in steady flow. The pressure drop in steady flow was found to be constant and always positive while for pulsatile flow the pressure drop change with time. And large turbulent kinetic energy, dissipation rate were found to be located in the region where the flow passes through the orifics ring. The maximum turbulent kinetic energy, generally occurs along the shear layer where the velocity gradient is large.