共查询到20条相似文献,搜索用时 15 毫秒
1.
We propose a new model and a solution method for two‐phase two‐fluid compressible flows. The model involves six equations obtained from conservation principles applied to a one‐dimensional flow of gas and liquid mixture completed by additional closure governing equations. The model is valid for pure fluids as well as for fluid mixtures. The system of partial differential equations with source terms is hyperbolic and has conservative form. Hyperbolicity is obtained using the principles of extended thermodynamics. Features of the model include the existence of real eigenvalues and a complete set of independent eigenvectors. Its numerical solution poses several difficulties. The model possesses a large number of acoustic and convective waves and it is not easy to upwind all of these accurately and simply. In this paper we use relatively modern shock‐capturing methods of a centred‐type such as the total variation diminishing (TVD) slope limiter centre (SLIC) scheme which solve these problems in a simple way and with good accuracy. Several numerical test problems are displayed in order to highlight the efficiency of the study we propose. The scheme provides reliable results, is able to compute strong shock waves and deals with complex equations of state. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
2.
Consideration is given in this paper to the numerical solution of the transient two‐phase flow in rigid pipelines. The governing equations for such flows are two coupled, non‐linear, hyperbolic, partial differential equations with pressure dependent coefficients. The fluid pressure and velocity are considered as two principle dependent variables. The fluid is a homogeneous gas–liquid mixture for which the density is defined by an expression averaging the two‐component densities where a polytropic process of the gaseous phase is admitted. Instead of the void fraction, which varies with the pressure, the gas–fluid mass ratio (or the quality) is assumed to be constant, and is used in the mathematical formulation. The problem has been solved by the method of non‐linear characteristics and the finite difference conservative scheme. To verify their validity, the computed results of the two numerical techniques are compared for different values of the quality, in the case where the liquid compressibility and the pipe wall elasticity are neglected. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
3.
Yang‐Yao Niu 《国际流体数值方法杂志》2001,36(3):351-371
In this study, the advection upwinding splitting method (AUSM) is modified for the resolution of two‐phase mixtures with interfaces. The compressible two‐fluid model proposed by Saurel and Abgrall is chosen as the model equations. Dense and dilute phases are described in terms of the volume fraction and equations of state to represent multi‐phase mixtures. Test cases involving an air–water shock tube, water faucet, and dilute particulate turbulent flows through a 90° bend are used to verify the current work. It is shown that the AUSM based on flux differences (AUSMD) contains the mechanism to correctly capture the contact discontinuity and interfaces between phases. In addition, a successful application to dilute particulate turbulence flows by the AUSMD is demonstrated. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
4.
This paper describes the development of numerical method for the solution of condensing steam flow in internal aerodynamics problems. The numerical method is based on the fractional step method, where the resulting set of ODEs is solved by the two‐stage Runge–Kutta method and the homogeneous set of PDEs by a finite volume method. The flow does not contain both phases (gas and liquid) in the whole domain, therefore we discuss properties of used finite volume methods in several cases of single‐phase transonic flow in a channel and a turbine cascade. We present numerical results of two‐phase flow of condensing steam in 2D nozzle achieved by several numerical methods and show the differences in results caused by numerical diffusion. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
5.
A conservative local interface sharpening scheme has been developed for the constrained interpolation profile method with the conservative semi‐Lagrangian scheme, because the conservative semi‐Lagrangian scheme does not feature a mechanism to control the interface thickness, thus causing an increase of numerical error with the advance of the time step. The proposed sharpening scheme is based on the conservative level set method proposed by Olsson and Kreiss. However, because their method can cause excessive deformation of the free‐surface in certain circumstances, we propose an improvement of the method by developing a local sharpening technique. Several advection tests are presented to assess the correctness of the advection and the improved interface sharpening scheme. This is followed by the validations of dam‐breaking flow and the rising bubble flows. The mass of the fluid is exactly conserved and the computed terminal velocity of the rising bubble agrees well with the experiments compared with other numerical methods such as the volume of fluid method (VOF), the front tracking method, and the level set method. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
6.
This paper investigates the performance of preconditioned Krylov subspace methods used in a previously presented two‐fluid model developed for the simulation of separated and intermittent gas–liquid flows. The two‐fluid model has momentum and mass balances for each phase. The equations comprising this model are solved numerically by applying a two‐step semi‐implicit time integration procedure. A finite difference numerical scheme with a staggered mesh is used. Previously, the resulting linear algebraic equations were solved by a Gaussian band solver. In this study, these algebraic equations are also solved using the generalized minimum residual (GMRES) and the biconjugate gradient stabilized (Bi‐CGSTAB) Krylov subspace iterative methods preconditioned with incomplete LU factorization using the ILUT(p, τ) algorithm. The decrease in the computational time using the iterative solvers instead of the Gaussian band solver is shown to be considerable. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
7.
Current existing main nuclear thermal‐hydraulics (T‐H) system analysis codes, such as RALAP5, TRACE, and CATHARE, play a crucial role in the nuclear engineering field for the design and safety analysis of nuclear reactor systems. However, two‐fluid model used in these T‐H system analysis codes is ill posed, easily leading to numerical oscillations, and the classical first‐order methods for temporal and special discretization are widely employed for numerical simulations, yielding excessive numerical diffusion. Two‐fluid seven‐equation two‐pressure model is of particular interest due to the inherent well‐posed advantage. Moreover, high‐order accuracy schemes have also attracted great attention to overcome the challenge of serious numerical diffusion induced by low‐order time and space schemes for accurately simulating nuclear T‐H problems. In this paper, the semi‐implicit solution algorithm with high‐order accuracy in space and time is developed for this well‐posed two‐fluid model and the robustness and accuracy are verified and assessed against several important two‐phase flow benchmark tests in the nuclear engineering T‐H field, which include two linear advection problems, the oscillation problem of the liquid column, the Ransom water faucet problem, the reversed water faucet problem, and the two‐phase shock tube problem. The following conclusions are achieved. (1) The proposed semi‐implicit solution algorithm is robust in solving two‐phase flows, even for fast transients and discontinuous solutions. (2) High‐order schemes in both time and space could prevent excessive numerical diffusion effectively and the numerical simulation results are more accurate than those of first‐order time and space schemes, which demonstrates the advantage of using high‐order schemes. 相似文献
8.
9.
Theoretical studies have been made to determine the pressure drops caused by abrupt flow area expansion/contraction in small circular pipes for two‐phase flow of air and water mixtures at room temperature and near atmospheric pressure. Two‐phase computational fluid dynamics (CFD) calculations, using Eulerian–Eulerian model (with the air phase being compressible for pipe contraction case) are employed to calculate the pressure drop across sudden expansion and contraction. The pressure drop is determined by extrapolating the computed pressure profiles upstream and downstream of the expansion/contraction. The larger and smaller tube diameters are 1.6 and 0.84 mm, respectively. Computations have been performed with single‐phase water and air, and two‐phase mixtures in a range of Reynolds number (considering all‐liquid flow) from 1000 to 12 000 and flow quality from 1.2 × 10?3 to 1.6 × 10?2. The numerical results are validated against experimental data from the literature and are found to be in good agreement. The expansion and contraction loss coefficients are found to be different for single‐phase flow of air and water, and they agreed reasonably well with the commonly used theoretical predictions. Based on the numerical results as well as experimental data, correlations are developed for two‐phase flow pressure drops caused by the flow area contraction as well as expansion. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
10.
This paper presents the extension of a high‐resolution conservative scheme to the one‐dimensional one‐pressure six‐equation two‐fluid flow model. Only mixtures of water and air have been considered in this study, both fluids have been characterized using simple equations of state, namely stiffened gas for the liquid phase and perfect gas for the gas phase. The resulting scheme is explicit and first‐order accurate in space and time. A second‐order version of the scheme has also been derived using the MUSCL strategy and slope limiters. Some numerical results show the good capabilities of this type of schemes in the solution of discontinuities in two‐fluid flow problems, all of them are based on water/air numerical benchmarks widely used in the two‐phase flow literature. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
11.
A method is introduced to discretize irregular and complex two‐dimensional fractured media. The geometry of the fractured media is first analysed by searching and treating the complex configurations. Based on that, the method generated a good mesh quality and allows for including finer grids. An incompressible two‐phase flow problem is solved to compare the developed method and a public method based on the approximation of a 1D fracture by the edges of a 2D finite element grid of the porous media. The comparison showed that the developed method (i) represents better the fractured domain by maintaining the geometric integrity of input surfaces and geologic data, (ii) provides, for sample and complex fractured domains, excellent and more accurate results, and (iii) is much less sensitive to the grid sizes. Furthermore, the method has to be more efficient than the other methods for transport problems and has to provide better predictable results; this is mainly based on point (ii) and because the method produces optimal triangular grids. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
12.
O. Coutier‐Delgosha R. Fortes‐Patella J. L. Reboud N. Hakimi C. Hirsch 《国际流体数值方法杂志》2005,48(2):135-167
A computational method is proposed to simulate 3D unsteady cavitating flows in spatial turbopump inducers. It is based on the code FineTurbo, adapted to take into account two‐phase flow phenomena. The initial model is a time‐marching algorithm devoted to compressible flow, associated with a low‐speed preconditioner to treat low Mach number flows. The presented work covers the 3D implementation of a physical model developed in LEGI for several years to simulate 2D unsteady cavitating flows. It is based on a barotropic state law that relates the fluid density to the pressure variations. A modification of the preconditioner is proposed to treat efficiently as well highly compressible two‐phase flow areas as weakly compressible single‐phase flow conditions. The numerical model is applied to time‐accurate simulations of cavitating flow in spatial turbopump inducers. The first geometry is a 2D Venturi type section designed to simulate an inducer blade suction side. Results obtained with this simple test case, including the study of its general cavitating behaviour, numerical tests, and precise comparisons with previous experimental measurements inside the cavity, lead to a satisfactory validation of the model. A complete three‐dimensional rotating inducer geometry is then considered, and its quasi‐static behaviour in cavitating conditions is investigated. Numerical results are compared to experimental measurements and visualizations, and a promising agreement is obtained. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
13.
We derive a partially analytical Roe scheme with wave limiters for the compressible six‐equation two‐fluid model. Specifically, we derive the Roe averages for the relevant variables. First, the fluxes are split into convective and pressure parts. Then, independent Roe conditions are stated for these two parts. These conditions are successively reduced while defining acceptable Roe averages. For the convective part, all the averages are analytical. For the pressure part, most of the averages are analytical, whereas the remaining averages are dependent on the thermodynamic equation of state. This gives a large flexibility to the scheme with respect to the choice of equation of state. Furthermore, this model contains nonconservative terms. They are a challenge to handle right, and it is not the object of this paper to discuss this issue. However, the Roe averages presented in this paper are fully independent from how those terms are handled, which makes this framework compatible with any treatment of nonconservative terms. Finally, we point out that the eigenspace of this model may collapse, making the Roe scheme inapplicable. This is called resonance. We propose a fix to handle this particular case. Numerical tests show that the scheme performs well. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
14.
15.
16.
We develop and analyse an improved version of the multi‐stage (MUSTA) approach to the construction of upwind Godunov‐type fluxes whereby the solution of the Riemann problem, approximate or exact, is not required. The new MUSTA schemes improve upon the original schemes in terms of monotonicity properties, accuracy and stability in multiple space dimensions. We incorporate the MUSTA technology into the framework of finite‐volume weighted essentially nonoscillatory schemes as applied to the Euler equations of compressible gas dynamics. The results demonstrate that our new schemes are good alternatives to current centred methods and to conventional upwind methods as applied to complicated hyperbolic systems for which the solution of the Riemann problem is costly or unknown. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
17.
The effect of wavelength and relative velocity on the disturbed interface of two‐phase stratified regime is modeled and discussed. To analyze the stability, a small perturbation is imposed on the interface. Growth or decline of the disturbed wave, relative velocity, and surface tension with respect to time will be discussed numerically. Newly developed scheme applied to a two‐dimensional flow field and the governing Navier–Stokes equations in laminar regime are solved. Finite volume method together with non‐staggered curvilinear grid is a very effective approach to capture interface shape with time. Because of the interface shape, for any time advancement, a new grid is performed separately on each stratified field, liquid, and gas regime. The results are compared with the analytical characteristics method and one‐dimensional modeling. This comparison shows that solving the momentum equation including viscosity term leads to physically more realistic results. In addition, the newly developed method is capable of predicting two‐phase stratified flow behavior more precisely than one‐dimensional modeling. It was perceived that the surface tension has an inevitable role in dissipation of interface instability and convergence of the two‐phase flow model. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
18.
Eric Goncalvès 《国际流体数值方法杂志》2012,68(8):1053-1072
A compressible, multiphase, one‐fluid Reynolds‐averaged Navier–Stokes solver has been developed to study turbulent cavitating flows. The interplay between turbulence and cavitation regarding the unsteadiness and structure of the flow is complex and not well understood. This constitutes a critical point to accurately simulate the dynamic behavior of sheet cavities. In the present study, different formulations based on a k ? ? transport‐equation model are investigated and a scale‐adaptive formulation is proposed. Numerical results are given for a Venturi geometry and comparisons are made with experimental data. The scale‐adaptive model shows several improvements compared with standard turbulence models. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
19.
Charalampos Kouris Yannis Dimakopoulos Georgios Georgiou John Tsamopoulos 《国际流体数值方法杂志》2002,39(1):41-73
A Galerkin/finite element and a pseudo‐spectral method, in conjunction with the primitive (velocity‐pressure) and streamfunction‐vorticity formulations, are tested for solving the two‐phase flow in a tube, which has a periodically varying, circular cross section. Two immiscible, incompressible, Newtonian fluids are arranged so that one of them is around the axis of the tube (core fluid) and the other one surrounds it (annular fluid). The physical and flow parameters are such that the interface between the two fluids remains continuous and single‐valued. This arrangement is usually referred to as Core‐Annular flow. A non‐orthogonal mapping is used to transform the uneven tube shape and the unknown, time dependent interface to fixed, cylindrical surfaces. With both methods and formulations, steady states are calculated first using the Newton–Raphson method. The most dangerous eigenvalues of the related linear stability problem are calculated using the Arnoldi method, and dynamic simulations are carried out using the implicit Euler method. It is shown that with a smooth tube shape the pseudo‐spectral method exhibits exponential convergence, whereas the finite element method exhibits algebraic convergence, albeit of higher order than expected from the relevant theory. Thus the former method, especially when coupled with the streamfunction‐vorticity formulation, is much more efficient. The finite element method becomes more advantageous when the tube shape contains a cusp, in which case the convergence rate of the pseudo‐spectral method deteriorates exhibiting algebraic convergence with the number of the axial spectral modes, whereas the convergence rate of the finite element method remains unaffected. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
20.
In this paper, several numerical schemes are extended to obtain approximate solutions to the system of equations encountered in the analysis of multiphase mixtures of gas and particles. Both dense and dilute mixtures are studied, the gas is modelled as a perfect gas and the solid is considered incompressible. Although the tests employed throughout this work for studying the behaviour of the schemes are essentially one dimensional, the finite volume method developed permits its application to multidimensional problems in unstructured grids. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献