• 한국유체기계학회 논문집
• /
• 제4권3호
• /
• 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년도 추계 학술대회논문집
• /
• 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년도 유체기계 연구개발 발표회 논문집
• /
• 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년도 추계 학술대회논문집
• /
• 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
• /
• 제10권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.

• 대한기계학회논문집B
• /
• 제40권3호
• /
• pp.181-189
• /
• 2016

수평 관군에 대하여 수직이고 상향으로 흐르는 2상 유동에 의한 감쇠비를 예측하기 위한 해석모델이 Sim에 의하여 개발되었다. 이 모델에서 평가된 2상 유동의 기공률, 압력손실 등의 유동변수는 기존의 실험식을 사용하여 계산하였다. 그러나 관군의 경우에 사용하기에는 약간의 개선이 요구된다. 따라서 관군 내에 흐르는 2상 유동의 유동 변수에 대한 더 많은 정보를 획득하기 위하여 실험적으로 연구할 필요가 있다. 실험은 공기 - 물의 2상 유동이 흐르는 정사각형 배열 관군에서의 압력계수와 2상 유동 마찰승수를 계산하기 위하여 수행되었다. 피치 직경 비는 1.35이었고, 실린더의 직경은 18 mm이다. 압력센서와 신호처리 장치를 이용하여 관군에서의 압력차를 측정하였다. 2상 유동 마찰승수와 오일러수를 계산하기 위하여 관군에 적용되는 비균질 유동의 기공률은 Feenstra 등의 실험식을 사용하여 계산하였다. 균질과 비균질 2상 유동의 마찰승수와 오일러의 수를 실험적으로 구하고 Sim의 어림적 모델에 근거한 이론적 해석 결과와 비교 분석하였다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제20권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.

• 한국항공우주학회지
• /
• 제30권3호
• /
• pp.27-36
• /
• 2002

AUSM계열 수치기법의 수치적 불안정성에 대한 원인과 해결방안에 대한 연구를 수행하였다. Euler 유동에서 수치적 불안정성은 제어면에 수직한 방향의 유동속도가 영인 영역에서 발생하며 이 영역에서 Eule r 방정식은 근본적으로 부정해를 가지게 되어 무수히 많은 해를 가지게 된다. 지배방정식 자체로는 유일해를 찾는 것이 불가능하고 주위의 유동조건이나 외부교란에 의해 유일해를 결정하게 된다. 이러한 특징은 충격파 영역에서 교란이 존재할 경우 초기 상태에 대한 정보를 상실하게 되어 충격파 불안정성을 유발하게 된다. slip유동을 정확히 계산할 수 있는, 즉 유일해를 결정할 수 없는, 수치기법은 충격파 불안정성을 근본적으로 제거할 수 없다.

• Nuclear Engineering and Technology
• /
• 제44권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.

• 충청수학회지
• /
• 제31권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.