• Nuclear Engineering and Technology
• /
• 제33권4호
• /
• pp.409-422
• /
• 2001

The dynamic character of a system of the governing differential equations for the one- dimensional two-fluid model, where the momentum flux parameters are employed to consider the velocity and void fraction distribution in a flow channel, is investigated. In response to a perturbation in the form of a＇traveling wave, a linear stability analysis is peformed for the governing differential equations. The expression for the growth factor as a function of wave number and various flow parameters is analytically derived. It provides the necessary and sufficient conditions for the stability of the one-dimensional two-fluid model in terms of momentum flux parameters. It is demonstrated that the one-dimensional two-fluid model employing the physical momentum flux parameters for the whole range of dispersed flow regime, which are determined from the simplified velocity and void fraction profiles constructed from the available experimental data and $C_{o}$ correlation, is stable to the linear perturbations in all wave-lengths. As the basic form of the governing differential equations for the conventional one-dimensional two-fluid model is mathematically ill posed, it is suggested that the velocity and void distributions should be properly accounted for in the one-dimensional two-fluid model by use of momentum flux parameters.s.

• 한국전산유체공학회지
• /
• 제12권4호
• /
• pp.20-27
• /
• 2007

Two-dimensional stagnation flow toward a plane wall coated with magnetic fluid of uniform thickness is investigated. The flow field is represented as a similarity solution of the Navier-Stokes equation for this incompressible laminar flow. The resulting third order ordinary differential equation is solved numerically by using the shooting method and by determining two shooting parameters so as to satisfy the boundary and interface conditions. Features of the flow including streamline patterns are investigated for the varying values of density ratio, viscosity ratio, and Reynolds number. An adverse flow with double eddy pair in magnetic fluid region is found to emerge as the Reynolds number becomes higher than a threshold value. The results for the interface velocity, interface and wall shear stress, and boundary layer and displacement thickness are also presented.

• Nuclear Engineering and Technology
• /
• 제54권4호
• /
• pp.1495-1507
• /
• 2022

The improvement of thermal-hydraulic analysis techniques is essential to ensure the safety and reliability of nuclear power plants. The one-dimensional two-fluid model has been adopted in state-of-the-art thermal-hydraulic system codes. Current constitutive equations used in the system codes reach a mature level. Some exceptions are the partition method of wall friction in the momentum equation of the two-fluid model and the interfacial drag force model for a horizontal two-phase flow. This study is focused on deriving the partition method of wall friction in the momentum equation of the two-fluid model and modeling the interfacial drag force model for a horizontal bubbly flow. The one-dimensional momentum equation in the two-fluid model is derived from the local momentum equation. The derived one-dimensional momentum equation demonstrates that total wall friction should be apportioned to gas and liquid phases based on the phasic volume fraction, which is the same as that used in the SPACE code. The constitutive equations for the interfacial drag force are also identified. Based on the assessments, the Rassame-Hibiki correlation, Hibiki-Ishii correlation, Ishii-Zuber correlation, and Rassame-Hibiki correlation are recommended for computing the distribution parameter, interfacial area concentration, drag coefficient, and relative velocity covariance of a horizontal bubbly flow, respectively.

• Journal of Advanced Marine Engineering and Technology
• /
• 제29권1호
• /
• pp.25-33
• /
• 2005

When a circular tube with uniform heat generation within the wall was placed in a cross flow, heat flows by conduction in the circumferential direction due to the asymmetric nature of the fluid flow around the perimeter of the circular tube The circumferential heat flow affects the wall temperature distribution to such an extent that. in some cases, significantly different results may be obtained for geometrically similar surfaces. In the present investigation, the effect of circumferential wall heat conduction is investigated for forced convection around circular tube in cross flow of air and water Two-dimensional temperature distribution $T_w(r,{\theta})$ is calculated through the numerical analysis. The difference between one-dimensional and two-dimensional solutions is demonstrated on the graph of local heat transfer coefficients. It is observed that the effect of working fluid is very remarkable.

• Journal of Mechanical Science and Technology
• /
• 제17권7호
• /
• pp.1073-1082
• /
• 2003

Many researchers have investigated the blood flow characteristics through bileaflet mechanical heart valves using computational fluid dynamics (CFD) models. Their numerical approach methods can be classified into three types; steady flow analysis, pulsatile flow analysis with fixed leaflets, and pulsatile flow analysis with moving leaflets. The first and second methods have been generally employed for two-dimensional and three-dimensional calculations. The pulsatile flow analysis interacted with moving leaflets has been recently introduced and tried only in two-dimensional analysis because this approach method has difficulty in considering simultaneously two physics of blood flow and leaflet behavior interacted with blood flow. In this publication, numerical calculation for pulsatile flow with moving leaflets using a fluid-structure interaction method has been performed in a three-dimensional geometry. Also, pulsatile flow with fixed leaflets has been analyzed for comparison with the case with moving leaflets. The calculated results using the fluid-structure interaction model have shown good agreements with results visualized by previous experiments. In peak systole. calculations with the two approach methods have predicted similar flow fields. However, the model with fixed leaflets has not been able to predict the flow fields during opening and closing phases. Therefore, the model with moving leaflets is rigorously required for advanced analysis of flow fields.

• 대한조선학회지
• /
• 제24권3호
• /
• pp.9-16
• /
• 1987

In this paper, the fluid flow in the two-dimensional tank is analyzed by the Finite Difference Method. The Navier-Stokes equation is modified for the tank fixed coordinate system. For the treatment of the free surface, the Volume of Fluid Method by Hirt and Nichols is adopted. The continuity equation and the Poisson equation which is derived from the Navier-Stokes equation to find the pressure are solved by the Successive-Line-Overrelaxation Method. The comparison of the calculated results with experimental data show a favorable agreement. The fluid flow in the two-dimensional tank can be predicted reasonably before the free surface reaches breaking by this numerical method.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2005년도 추계 학술대회논문집
• /
• pp.183-186
• /
• 2005

We developed an upwind numerical formulation based on the eigenvalues of the approximate Jacobian matrix in order to solve the hyperbolic conservation laws governing the two-fluid two-phase flow models. We obtained eight analytic eigenvalues in the two dimensions that can be used for estimate of the wave speeds essential in constructing an upwind numerical method. Two-dimensional underwater cavitation in a flow past structural shapes or by underwater explosion can be solved using this method. We present quantitative prediction of cavitation for the water tunnel wall and airfoils that has both experimental data as well as numerical results by other numerical methods and models.

• International Journal of Fluid Machinery and Systems
• /
• 제3권1호
• /
• pp.58-66
• /
• 2010

A steam control valve is used to control the flow from the steam generator to the steam turbine in thermal and nuclear power plants. During startup and shutdown of the plant, the steam control valve is operated under a partial flow conditions. In such conditions, the valve opening is small and the pressure deference across the valve is large. As a result, the flow downstream of the valve is composed of separated unsteady transonic jets. Such flow patterns often cause undesirable large unsteady fluid force on the valve head and downstream pipe system. In the present study, various flow patterns are investigated in order to understand the characteristics of the unsteady flow around the valve. Experiments are carried out with simplified two-dimensional valve models. Two-dimensional unsteady flow simulations are conducted in order to understand the experimental results in detail. Scale effects on the flow characteristics are also examined. Results show three types of oscillating flow pattern and three types of static flow patterns.

• Journal of Biosystems Engineering
• /
• 제20권1호
• /
• pp.66-72
• /
• 1995

A numerical study for a two dimensional multi-impingement jet with crossflow of the spent fluid has been carried out. To study the flow characteristics especially in the jet flow region, three different distributions of mass flow rate at 5-jet exits were assumed. For each distribution, various Reynolds numbers ranging from laminar to turbulent flows were considered. Calculations drew the following items as conclusion. 1) A periodical fully developed flow was observed from the third protrusion. This was also observed from previous experimentally by Whidden at al. The Nessult number at the protrusion surface increased mildly as going downstream. 2) The low Reynolds number turbulence model of Launder and Sharma was found to be adequate for the prediction of fluid flow and heat transfer characteristics of two dimensional multi-jet configuration. 3) The Nusselt number at the protrusion surface was nearly proportional to the square root of the Reynolds number.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 추계학술대회논문집B
• /
• pp.591-596
• /
• 2000

Flow through compliant tubes with linear taper in wall thickness is numerically simulated by finite element analysis. Two models are examined: a planar two-dimensional channel, and an axisymmetric tube. For verification of the numerical method, flow through a compliant stenotic vessel is simulated and compared to existing experimental data. Computational results for an axisymmetric tube show that as cross-sectional area falls with a reduction in downstream pressure, flow rate increases and reaches a maximum when the speed index (mean velocity divided by wave speed) is near unity at the point of minimum cross-section area, indicative of wave speed flow limitation or "choking" (flow speed equals wave speed) in previous one-dimensional studies. For further reductions in downstream pressure, flow rate decreases. Cross-sectional narrowing is significant but localized. When the ratio of downstream-to-upstream wall thickness is ${\le}$ 2 the area throat is located near the downstream end; as wall taper is increased to ${\ge}$ 3 the constriction moves to the upstream end of the tube. In the planar two-dimensional channel, area reduction and flow limitation are also observed when outlet pressure is decreased. In contrast to the axisymmetric case, however, the elastic wall in the two-dimensional channel forms a smooth concave surface with the area throat located near the mid-point of the elastic wall. Though flow rate reaches a maximum and then falls, the flow does not appear to be choked.