• The KSFM Journal of Fluid Machinery
• /
• v.4 no.3 s.12
• /
• pp.53-61
• /
• 2001

Vortex type Fluidic Device(FD) which is installed at the bottom of Safety Injection Tank(SIT) controls the discharge flow rate from the tank. In case of loss of coolant accident the injection water flows into primary system in two steps; initial high flow rate for certain period of time and subsequent low flow rate. By two-step control of the discharge flow rate, FD can ensure the effective use of water in the tank. A small-scale FD has been tested to obtain a required flow characteristics maintaining full pressure and height of prototype, which are the major contributing parameters. Through the testing of many different arrangements of internal geometry of FD, most appropriate one was selected and its performance data was obtained. As characteristics of FD, time dependent Euler number, flow rate and pressure are presented and discussed. Also a method to predict the full size FD is presented.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.175-181
• /
• 1995

Three-Dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for the various contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreements.

• 유체기계공업학회:학술대회논문집
• /
• 2001.11a
• /
• pp.76-80
• /
• 2001

Unsteady compressible Euler equation is solved and the high-order, high-resolution numerical solver, physical boundary condition, adaptive nonlinear artificial dissipation model and conformal mapping are applied to computation of steady transonic flow and unsteady acoustics. The acoustic characteristics of axi-symmetric duct and two dimensional straight/S channel are studied and the computation results shows good agreements with linear analysis. In transonic case, local time stepping and canceling-the-residual techniques are used for convergence acceleration. The aspect of flow and acoustics in S-channel and the Pattern of noise radiation is changed by inflow Mach no. and static pressure at fan-face.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.153-159
• /
• 1998

The Newton-Krylov method on the unstructured grid flow solver using the cell-centered spatial discretization oi compressible Euler equations is presented. This flow solver uses the reconstructed primitive variables to get the higher order solutions. To get the quadratic convergence of Newton method with this solver, the careful linearization of face flux is performed with the reconstructed flow variables. The GMRES method is used to solve large sparse matrix and to improve the performance ILU preconditioner is adopted and vectorized with level scheduling algorithm. To get the quadratic convergence with the higher order schemes and to reduce the memory storage. the matrix-free implementation and Barth's matrix-vector method are implemented and compared with the traditional matrix-vector method. The convergence and computing times are compared with each other.

• International Journal of Fluid Machinery and Systems
• /
• v.10 no.2
• /
• pp.154-163
• /
• 2017

Performance characteristics in pump mode of pump-turbines are vital for the safe and effective operation of pumped storage power plants. However, the head characteristics are different under different guide vane openings. In this paper, 3-D steady simulations were performed under 13mm, 19mm and 25mm guide vane openings. Three groups of operating points under the three GVOs were chosen based on experimental validation to investigate the influence of guide vane setting on flow patterns upstream and downstream. The results reveal that, the guide vane setting will obviously change the flow pattern downstream, which in turn influences the flow upstream. It shows a strong effect on hydraulic loss (power dissipation) in the guide and stay vanes. It is also found that the hydraulic loss mainly comes from the flow separation and vortices. In addition, in some operating conditions, the change of guide vane opening will change the flow angle at the runner inlet and outlet, which will change the Euler momentum (power input). The joint action of Euler momentum and hydraulic loss results in the change of the head characteristics.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.3
• /
• pp.181-189
• /
• 2016

An analytical model was developed by Sim to estimate the two-phase damping ratio for upward two-phase flow perpendicular to horizontal tube bundles. The parameters of two-phase flow, such as void fraction and pressure loss evaluated in the model, were calculated based on existing experimental formulations. However, it is necessary to implement a few improvements in the formulations for the case of tube bundles. For the purpose of the improved formulation, we need more information about the two-phase parameters, which can be found through experimental test. An experiment is performed with a typical normal square array of cylinders subjected to the two-phase flow of air-water in the tube bundles, to calculate the two-phase Euler number and the two-phase friction multiplier. The pitch-to-diameter ratio is 1.35 and the diameter of cylinder is 18mm. Pressure loss along the flow direction in the tube bundles is measured with a pressure transducer and data acquisition system to calculate the two-phase Euler number and the two-phase friction multiplier. The void fraction model by Feenstra et al. is used to estimate the void fraction of the two-phase flow in tube bundles. The experimental results of the two phase friction multiplier and two-phase Euler number for homogeneous and non-homogeneous two-phase flows are compared and evaluated against the analytical results given by Sim's model.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.3
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.3
• /
• pp.27-36
• /
• 2002

Numerical stability is studied based on numerics and mathematics. It is frequently observed in the region where velocity is zero. In that region, the Euler equation have numerous solutions and, thus, it is impossible to determine a unique solution with only governing equations. However, a unique solution can be determined by additional outer flow conditions or outer numerical discontinuity calculation since the information or a unique solution under undisturbed conditions is lost by disturbances. In this reason, the numerical scheme comsistent with Euler equations cannot remove shock instability completely.

• Nuclear Engineering and Technology
• /
• v.44 no.7
• /
• pp.735-744
• /
• 2012

In order to quantify the flow distribution characteristics of APR+ reactor, a test was performed on a test facility, ACOP ($\underline{A}$PR+ $\underline{C}$ore Flow & $\underline{P}$ressure Test Facility), having a length scale of 1/5 referring to the prototype plant. The major parameters are core inlet flow and outlet pressure distribution and sectional pressure drops along the major flow path inside reactor vessel. To preserve the flow characteristics of prototype plant, the test facility was designed based on a preservation of major flow path geometry. An Euler number is considered as primary dimensionless parameter, which is conserved with a 1/40.9 of Reynolds number scaling ratio. ACOP simplifies each fuel assembly into a hydraulic simulator having the same axial flow resistance and lateral cross flow characteristics. In order to supply boundary condition to estimate thermal margins of the reactor, the distribution of inlet core flow and core exit pressure were measured in each of 257 fuel assembly simulators. In total, 584 points of static pressure and differential pressures were measured with a limited number of differential pressure transmitters by developing a sequential operation system of valves. In the current study, reactor flow characteristics under the balanced four-cold leg flow conditions at each of the cold legs were quantified, which is a part of the test matrix composing the APR+ flow distribution test program. The final identification of the reactor flow distribution was obtained by ensemble averaging 15 independent test data. The details of the design of the test facility, experiment, and data analysis are included in the current paper.

• Journal of the Chungcheong Mathematical Society
• /
• v.31 no.4
• /
• pp.363-368
• /
• 2018

The aim of this work is to derive a partial differential equation that explains the movement of vortex lines on a vortex trajectory surface in a three dimensional incompressible inviscid flow.