• Nuclear Engineering and Technology
• /
• v.32 no.2
• /
• pp.169-179
• /
• 2000

A detailed numerical analysis of the evolution of thermal stratification in a curved piping system in a nuclear power plant is performed. A finite volume based thermal-hydraulic computer code has been developed employing a body-fitted, non-orthogonal curvilinear coordinate for this purpose. The cell-centered, non-staggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is prevented by the application of momentum interpolation method. The SIMPLE algorithm is employed for the pressure and velocity coupling, and the convection terms are approximated by a higher-order bounded scheme. The thermal-hydraulic computer code developed in the present study has been applied to the analysis of thermal stratification in a curved duct and some of the predicted results are compared with the available experimental data. It is shown that the predicted results agree fairly well with the experimental measurements and the transient formation of thermal stratification in a curved duct is also well predicted.

• Korean Journal of Computational Design and Engineering
• /
• v.10 no.5
• /
• pp.314-327
• /
• 2005

This paper describes the volumetric data modeling and analysis methods that employ volumetric NURBS or VNURBS that represents heterogeneous objects or fields in multidimensional space. For volumetric data modeling, we formulate the construction algorithms involving the scattered data approximation and the curvilinear grid data interpolation. And then the computational algorithms are presented for the geometric and mathematical analysis of the volume data set with the VNURBS model. Finally, we apply the modeling and analysis methods to various field applications including grid generation, flow visualization, implicit surface modeling, and image morphing. Those application examples verify the usefulness and extensibility of our VNUBRS representation in the context of volume modeling and analysis.

• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• pp.163-169
• /
• 1999

The numerical procedure has been developed for simulating incompressible viscous flow around a turbine stage with rotor-stator interaction. This study solves 2-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system. The Marker-and-Cell concept is applied to efficiently solve continuity equation. To impose an accurate boundary condition, O-H multiblocked grid system is generated. O-type grid and H-type grid is generated near and outer rotor-stator The cubic-spline interpolation is applied to handle a relative motion of a rotor to the stator. Turbulent flows have been modeled by the Baldwin- Lomax turbulent model. To validate present procedure, the time averaged pressure coefficients around the rotor and stator are compared with experiment and a good agreement obtained.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.900-907
• /
• 2000

The three-dimensional turbulent flow in curved pipes susceptible to flow-accelerated corrosion has been analyzed numerically to predict the pressure and shear stress distributions on the inner surface of the pipes. The analysis employs the body-fitted non-orthogonal curvilinear coordinate system and a standard $ {\kappa}-{\varepsilon}$ turbulence model with wall function method. The finite volume method is used to discretize the governing equations. The convection term is approximated by a high-resolution and bounded discretization scheme. The cell-centered, non-staggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is prevented by the application of a modified version of momentum interpolation scheme. The SIMPLE algorithm is employed for the pressure and velocity coupling. The numerical calculations have been performed for two curved pipes with different bend angles and curvature radii, and discussions have been made on the distributions of the primary and secondary flow velocities, pressure and shear stress on the inner surface of the pipe to examine applicability of the present analysis method. As the result it is seen that the method is effective to predict the susceptible systems or their local areas where the fluid velocity or local turbulence is so high that the structural integrity can be threatened by wall thinning degradation due to flow-accelerated corrosion.

• Journal of computational fluids engineering
• /
• v.7 no.1
• /
• pp.10-19
• /
• 2002

The non-staggered(collocated) grid approach in which all the solution variables are located at the centers of control volumes is very popular for incompressible flow analyses because of its numerical efficiency on the curvilinear or unstructured grids. Rhie and Chow's paper is the first in using non-staggered grid method for SIMPLE algorithm, where pressure weighted interpolation was used to prevent decoupling of pressure and velocity. But it has been known that this non-staggered grid method has stability problems when pressure fields are nonlinear like in natural convection flows. Also Rhie-Chow scheme generates large numerical diffusion near curved walls. The cause of these unwanted problems is too large pressure damping term compared to the magnitude of face velocity. In this study the magnitude of pressure damping term of Rhie-Chow's method is limited to 1∼10% of face velocity to prevent physically unreasonable solutions. The wall pressure extrapolation which is necessary for cell-centered FVM is another source of numerical errors. Some methods are applied in a unstructured FV solver and analyzed in view of numerical accuracy. Here, two natural convection problems are solved to check the effect of the Rhie-Chow's method on numerical stability. And numerical diffusion from Rhie-Chow's method is studied by solving the inviscid flow around a circular cylinder.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.10
• /
• pp.2535-2542
• /
• 1993

Two-noded curved beam elements, CMLC (field-consistent membrane and linear curvature) and IMLC(field-inconsistent membrane and linear curvature) are developed on the basis of Timoshenko's beam theory and curvilinear coordinate. The curved beam element is developed by the separation of the radial deflection into the bending deflection. In the CMLC element, field-consistent axial strain interpolation is adapted for removing the membrane locking. The CMLC element shows the rapid and stable convergence on the wide range of curved beam radius to thickness. The field-consistent axial strain and the separation of radial deformation produces the most efficient linear element possible.

• Proceedings of the KIEE Conference
• /
• 1998.07a
• /
• pp.146-148
• /
• 1998

For the shape design of electromagnetic devices using the FEM, the choice of design parameters influence to the success of the optimization process. If the design parameter distribution has a one to one corespondence with finite element model, we can encounter not only serious accuracy problem but also obtain a zigzag shape along the interface. The nodes between those design parameters can be parameterized by interpolating using one among many interpolation methods. The conventional parameterization of design parameters has a limit of application for shape, because design parameters and movable nodes are linearly intepolated. In this paper, using the B-spline curve that use to present any interfaces in computer graphics, the curvilinear parameterization between design parameters and node points is compared with the linear parameterization.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.3
• /
• pp.59-70
• /
• 1992

A Finite volume method for the computation of the two-dimensional, incompressible, steady, laminar Navier-Stokes equation is developed using a non-staggered grid system in a general curvilinear coordinate. The numerical pressure fluctuations, usually encountered when the non-staggered grid system is used, is suppressed by the momentum interpolation method. Flows around a NACA0012 foil section have been computed by the present method and the results show good agreements with other experimental and numerical ones.

• 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.

• Journal of Korea Water Resources Association
• /
• v.50 no.7
• /
• pp.455-465
• /
• 2017

Understanding of the two-dimensional velocity field is crucial in terms of analyzing various hydrodynamic and fluvial processes in the riverine environments. Until recently, many numerical models have played major roles of providing such velocity field instead of in-situ flow measurements, because there were limitations in instruments and methodologies suitable for efficiently measuring in the broad range of river reaches. In the last decades, however, the advent of modernized instrumentations started to revolutionize the flow measurements. Among others, acoustic Doppler current profilers (ADCPs) became very promising especially for accurately assessing streamflow discharge, and they are also able to provide the detailed velocity field very efficiently. Thus it became possible to capture the velocity field only with field observations. Since most of ADCPs measurements have been mostly conducted in the cross-sectional lines despite their capabilities, it is still required to apply appropriate interpolation methods to obtain dense velocity field as likely as results from numerical simulations. However, anisotropic nature of the meandering river channel could have brought in the difficulties for applying simple spatial interpolation methods for handling dynamic flow velocity vector, since the flow direction continuously changes over the curvature of the channel shape. Without considering anisotropic characteristics in terms of the meandering, therefore, conventional interpolation methods such as IDW and Kriging possibly lead to erroneous results, when they dealt with velocity vectors in the meandering channel. Based on the consecutive ADCP cross-sectional measurements in the meandering river channel. For this purpose, the geographic coordinate with the measured ADCP velocity was converted from the conventional Cartesian coordinate (x, y) to a curvilinear coordinate (s, n). The results from application of A-VIM showed significant improvement in accuracy as much as 41.5% in RMSE.