• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.4
• /
• pp.101-109
• /
• 1993

On the basis of Nikuradse laboratory experiments. two transition flow regimes are defined with respect to the characteristics of boundary layer. One is the transitional turbulent flow which has a transitional characteristics between smooth turbulent flow and rough turbulent flow, and the other may be called as transitional laminar flow which has transitional features between laminar flow and turbulent flow. The laboratory results of Nikuradse are carefully re-examined, and the flow regions are clearly defined. The velocity profile of the transitional turbulent flow is described by newly formulated equation, and the Darcy-Weisbach friction coefficient for the transitional turbulent flow is determined based on the theoretical form of profile equation, which is far better accurate than Colebrook-White equation. Difficulties still arise for the description of velocity profile when the flow undergoes from laminar to turbulent. In this case a linear interpolation procedure is proposed for the estimation of friction coefficient.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.474-481
• /
• 2001

This paper addresses a numerical method for predicting transient temperature distributions in the wall of a curved pipe subjected to internally thermal stratification flow. A simple and convenient numerical method of treating the unsteady conjugate heat transfer in the non-orthogonal coordinate systems is presented. The proposed method is implemented in a finite volume thermal-hydraulic computer code based on a cell-centered, non-staggered grid arrangement, the SIMPLEC algorithm, a higher-order bounded convection scheme, and the modified version of momentum interpolation method. Calculations are performed for the transient evolution of thermal stratification in two curved pipes, where the one has thick wall and the other has so thin wall that its presence can be negligible in the heat transfer analysis. The predicted results show that the thermally stratified flow and transient conjugate heat transfer in a curved pipe with a finite wall thickness can be satisfactorily analyzed by the present numerical method, and that the neglect of wall thickness in the prediction of pipe wall temperature distributions can provide unacceptably distorted results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.6
• /
• pp.619-627
• /
• 2004

In most industrial applications, the geometrical complexity is combined with the moving boundaries. These problems considerably increase the computational difficulties since they require, respectively, regeneration and deformation of the grid. As a result, engineering flow simulation is restricted. In order to solve this kind of problems the immersed boundary method was developed. In this study, the immersed boundary method is applied to the numerical simulation of stationary, rotating and oscillating cylinders in the 2-dimensional square cavity. No-slip velocity boundary conditions are given by imposing feedback forcing term to the momentum equation. Besides, this technique is used with a second-order accurate interpolation scheme in order to improve the accuracy of flow near the immersed boundaries. The governing equations for the mass and momentum using the immersed boundary method are discretized on the non-staggered grid by using the finite volume method. The results agree well with previous numerical and experimental results. This study presents the possibility of the immersed boundary method to apply to the complex flow experienced in the industrial applications. The usefulness of this method will be confirmed when we solve the complex geometries and moving bodies.

• The KSFM Journal of Fluid Machinery
• /
• v.14 no.5
• /
• pp.5-11
• /
• 2011

In this study, the effect of cone angle of seven-hole probe on the accuracy of measured flow angle and velocity has been investigated. The seven-hole probe consisted of seven 1mm OD stainless inner tubes and one 3mm ID stainless tube. Six cone angles of $45^{\circ}$, $60^{\circ}$, $75^{\circ}$, $90^{\circ}$, $105^{\circ}$ and $120^{\circ}$ were tested. Calibrations of the seven-hole probes were conducted within ${\pm}60^{\circ}$ range with the interval of $5^{\circ}$. Analysis results show that the effect on the cone angle was not significant on the accuracy of the measured flow angle, pressure, and velocity. However, the data reduction method had more effect on the measurement accuracy.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.5
• /
• pp.1661-1670
• /
• 1996

Incompressible flow over a backward-facing step is computed by low Reynolds number turbulence models in order to compare with direct simulation results. In this study, selected low Reynolds number 1st and 2nd (Algebraic Stress Model : ASM) moment closure turbulence models are adopted and compared with each other. Each turbulence model predicts different flow characteristics, different re-attachment point, velocity profiles and Reynolds stress distribution etc. Results by .kappa.-.epsilon. turbulence models indicate that predicted re-attachment lengths are shorter than those by standard model. Turbulent intensity and eddy viscosity by low Reynolds number .kappa.-.epsilon. models are still greater than DNS results. The results by algebraic stress model (ASM) are more reasonable than those by .kappa.-.epsilon. models. The convective scheme is QUICK (Quadratic Upstream Interpolation for Convective Kinematics) and SIMPLE algorithm is adopted. Reynolds number based on step height and inlet free stream velocity is 5100.

• Water Engineering Research
• /
• v.2 no.3
• /
• pp.179-185
• /
• 2001

The time consuming and expensive nature of experimental research on scouring processes caused by flowing water makes it attractive to develop numerical tools for the predication of the interaction of the fluid flow and the movable bed. In this paper the numerical simulation of scour by a wall jet is presented. The flow is assumed to be two-dimensional, and the alluvium is cohesionless. The solution process, repeated at each time step, involves simulation of a turbulent wall jet flow, solution of the convection-diffusion of sand concentration, and prediction of the bed deformation. For simulation of the jet flow, the governing equations for momentum, mass balance and turbulent parameters are solved by the finite volume method. The SIMPLE scheme with momentum interpolation is used for pressure correction. The convection-diffusion equation is solved for sediment concentration. A boundary condition for concentration at the bed, which takes into account the effect of bed-load, is implemented. The time rate of deposition and scour at the bed is obtained by solving the continuity equation for sediment. The shape and position of the scour hole and deposition of the bed material downstream of the hole appear realistic.

• Journal of Mechanical Science and Technology
• /
• v.20 no.10
• /
• pp.1741-1752
• /
• 2006

A combined Streamline Upwind Petrov-Galerkin method (SUPG) and segregated finite element algorithm for solving conjugate heat transfer problems where heat conduction in a solid is coupled with heat convection in viscous fluid flow is presented. The Streamline Upwind Petrov-Galerkin method is used for the analysis of viscous thermal flow in the fluid region, while the analysis of heat conduction in solid region is performed by the Galerkin method. The method uses the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the presented method is to consistently couple heat transfer along the fluid-solid interface. Four test cases, which are the conjugate Couette flow problem in parallel plate channel, the counter-flow in heat exchanger, the conjugate natural convection in a square cavity with a conducting wall, and the conjugate natural convection and conduction from heated cylinder in square cavity, are selected to evaluate efficiency of the presented method.

• Journal of Power System Engineering
• /
• v.10 no.4
• /
• pp.56-64
• /
• 2006

This paper shows the development process of a straight-type five-hole pressure probe for measuring three-dimensional flow velocity components. The data reduction method using a bi-cubic curve-fitting program in a new calibration map was introduced in this study. This new calibration map can be applied up to the application angle, ${\pm}55^{\circ}$ of a probe. As a result, for the application angle of ${\pm}45^{\circ}$, an error for yaw and pitch angles appeared from $-1.76^{\circ}\;to\;1.83^{\circ}$ and from $-1.91^{\circ}\;to\;1.75^{\circ}$, respectively. Moreover, an error for a vector magnitude and a static pressure compared with a dynamic one showed from -7.83% to 4.87% and from -0.73 to 0.77, respectively. Even though this data reduction method showed unsatisfactory errors in a vector magnitude, it resulted in an easy and simple application method. Especially, when it was applied to an actual flow field including a swirling flow, a good result came out on the whole. However, in order to obtain a better result, it is thought that a more sophisticated interpolation method needs to be introduced.

• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.1
• /
• pp.47-57
• /
• 2019

Recent CFD solvers in engineering have to treat geometrically complex domains and moving body problems. In ship hydrodynamics, flow around the stern and ship motions in waves are examples of such cases mentioned before. The unstructured grid scheme is successfully applied for these problems, but it has weakness of inefficient memory usage and intensive computational time as compared to the structured grid method. Overset grid scheme is one of the alternatives for structured grid system taking advantage of fast and memory efficiency. Overset grid scheme is especially useful for moving body problem because there is no need to re-mesh around the body. In this paper, we adopted the Suggar++, the grid connectivity and interpolation utility for the overlapping grid, to WAVIS which is the in-house flow solver of KRISO. Then we introduced some applications using the overset grid method for flow analysis around the ships. The computed results show that WAVIS with Suggar++ is practically feasible and has an advantages for moving geometry cases.