Flux difference splitting for open-channel flows

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow-water equations in open channels, together with an extension to two-dimensional flows. A linearized problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearized problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a one-dimensional dam-break problem, and to a problem of flow in a river whose geometry induces a region of supercritical flow. The scheme is also applied to a two-dimensional dam-break problem. The numerical results are compared with the exact solution, or other numerical results, where available.     

A novel non‐upwind,interconnected, multi‐grid,overlapping numerical procedure for problems involving fluid flow

A novel numerical procedure for heat, mass and momentum transfer in fluid flow is presented. The new scheme is passed on a non‐upwind, interconnected, multi‐grid, overlapping (NIMO) finite‐difference algorithm. In 2D flows, the NIMO algorithm solves finite‐difference equations for each dependent variable on four overlapping grids. The finite‐difference equations are formulated using the control‐volume approach, such that no interpolations are needed for computing the convective fluxes. For a particular dependent variable, four fields of values are produced. The NIMO numerical procedure is tested against the exact solution of two test problems. The first test problem is an oblique laminar 2D flow with a double step abrupt change in a passive scalar variable for infinite Peclet number. The second test problem is a rotating radial flow in an annular sector with a single step abrupt change in a passive scalar variable for infinite Peclet number. The NIMO scheme produced essentially the exact solution using different uniform and non‐uniform square and rectangular grids for 45 and 30° angle of inclination. All other schemes were unable to capture the exact solution, especially for the rectangular and non‐uniform grids. The NIMO scheme was also successful in predicting the exact solution for the rotating radial flow, using a uniform cylindrical‐polar coordinate grid. Copyright © 2008 John Wiley & Sons, Ltd.     

A coupled multigrid-domain-splitting technique for simulating incompressible flows in geometrically complex domains

This paper describes a domain decomposition numerical procedure for solving the Navier-Stokes equations in regions with complex geometries. The numerical method includes a modified version of QUICK (quadratic upstream interpolation convective kinematics) for the formulation of convective terms and a central difference scheme for the diffusion terms. A second-order-accurate predictor-corrector scheme is employed for the explicit time stepping. Although the momentum equations are solved independently on each subdomain, the pressure field is computed simultaneously on the entire flow field. A multigrid technique coupled with a Schwarz-like iteration method is devised to solve the pressure equation over the composite domains. The success of this strategy depends crucially on appropriate methods for specifying intergrid pressure boundary conditions on subdomains. A proper method for exchanging information among subdomains during the Schwarz sweep is equally important to the success of the multigrid solution for the overall pressure field. These methods are described and subsequently applied to two forced convection flow problems involving complex geometries to demonstrate the power and versatility of the technique. The resulting pressure and velocity fields exhibit excellent global consistency. The ability to simulate complex flow fields with this method provides a powerful tool for analysis and prediction of mixing and transport phenomenon.     

An alternate approach for modelling of transient vaporous cavitation

Simulation of cavitating flow has been a thrust area of research for long period due to its practical and economic importance. The major hurdle in developing a numerical model for such flows is the difficulty in representing the quick phase changes, in general, and the alternate change of flow from single phase to two phase and back, in particular. In this case, instability due to sharp variation of flow characteristics also restricts the development of numerical models. The present study demonstrates the use of a relatively simple formulation for the analysis of flow characteristics in a quasi‐rigid pipeline under abrupt phase changes due to cavitation. A popular scheme—MacCormack scheme—was used for developing a numerical solution for this problem. It uses the conservative form of the governing equations, viz. conservation of mass and momentum, the transport equation and the constitutive relationship. The model can handle variable properties of the water–vapour mixture, which is highly compressible. A newly introduced pressure under‐relaxation method overcomes the numerical instability due to sharp variation of flow characteristics during phase change. The model could predict the instant of occurrence of vapour pressure, duration of persistence of vapour pressure and the rise of pressure due to vapour collapse to satisfactory levels with published data and experimental results. Copyright © 2009 John Wiley & Sons, Ltd.     

Extension of a high‐resolution scheme to 1D liquid–gas flow

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.     

Numerical simulation of particle migration in asymmetric bifurcation channel

In this work we provide numerical validation of the particle migration during flow of concentrated suspension in asymmetric T-junction bifurcation channel observed in a recent experiment [1]. The mathematical models developed to explain particle migration phenomenon basically fall into two categories, namely, suspension balance model and diffusive flux model. These models have been successfully applied to explain migration behavior in several two-dimensional flows. However, many processes often involve flow in complex 3D geometries. In this work we have carried out numerical simulation of concentrated suspension flow in 3D bifurcation geometry using the diffusive flux model. The simulation method was validated with available experimental and theoretical results for channel flow. After validation of the method we have applied the simulation technique to study the flow of concentrated suspensions through an asymmetric T-junction bifurcation composed of rectangular channels. It is observed that in the span-wise direction inhomogeneous concentration distribution that develops upstream persists throughout the inlet and downstream channels. Due to the migration of particles near the bifurcation section there is almost equal partitioning of flow in the two downstream branches. The detailed comparison of numerical simulation results is made with the experimental data.     

Numerical simulation of flow‐induced noise using LES/SAS and Lighthill's acoustic analogy

We present a finite element (FE) formulation of Lighthill's acoustic analogy for the hybrid computation of noise generated by turbulent flows. In the present approach, the flow field is computed using large eddy simulation and scale adaptive simulation turbulence models. The acoustic propagation is obtained by solving the variational formulation of Lighthill's acoustic analogy with the FE method. In order to preserve the acoustic energy, we compute the inhomogeneous part of Lighthill's wave equation by applying the FE formulation on the fine flow grid. The resulting acoustic nodal loads are then conservatively interpolated to the coarser acoustic grid. Subsequently, the radiated acoustic field can be solved in both time and frequency domains. In the latter case, an enhanced perfectly matched layer technique is employed, allowing one to truncate the computational domain in the acoustic near field, without compromising the numerical solution. Our hybrid approach is validated by comparing the numerical results of the acoustic field induced by a corotating vortex pair with the corresponding analytical solution. To demonstrate the applicability of our scheme, we present full 3D numerical results for the computed acoustic field generated by the turbulent flow around square cylinder geometries. The sound pressure levels obtained compare well with measured values. Copyright © 2009 John Wiley & Sons, Ltd.     

Single Phase Flow in Partially Fissured Media

Totally fissured media in which the individual cells are isolated by the fissure system are effectively described by double porosity models with microstructure. Such models contain the geometry of the individual cells in the medium and the flux across their interface with the fissure system which surrounds them. We extend these results to a dual-permeability model which accounts for the secondary flux arising from direct cell-to-cell diffusion within the solid matrix. Homogenization techniques are used to construct a new macroscopic model for the flow of a single phase compressible fluid through a partially fissured medium from an exact but highly singular microscopic model, and it is shown that this macroscopic model is mathematically well posed. Preliminary numerical experiments illustrate differences in the behaviour of solutions to the partially fissured from that of the totally fissured case.     

Exit-flow velocity survey of two single-tangential-inlet vaneless turbine volutes

The most cited analytical technique for designing turbine volutes is to assume the throughflow is free from torque, although for this assumption to hold, the volute walls must lie near what would be streamlines in an unbounded free vortex-plus-sink flow. The single tangential inlet design, with inlet offset decreased and diameter increased to attain the weak exit swirl required by high specific speed turbines, deviates from such a shape, and the volute’s internal geometry is no longer torque-free. It is desired to know the actual time-averaged flow leaving such a volute, so that a rotor can be designed to compliment it.For two existing single tangential inlet volutes, time-averaged radial and tangential velocity and static pressure measurements of exit flow have been obtained on a cylindrical cut plane through the radial-inflow section using a three-port yawmeter in air. The Reynolds numbers based on inlet pipe mean conditions, around 10⁵, are well into the fully-turbulent regime and on the order of comparable water turbines.A comprehensive map of time-averaged exit flow of both volutes is presented. The integrated values of gross angular momentum flux change and total pressure loss coefficient are tabulated. Circumferential variation of flowrate and swirl strength highlight unexpected differences in outlet flow between the two volute designs. Results are presented alongside corresponding numerical results from the commercial package Fluent (Fluent, Inc., Lebanon, NH, USA) using Reynolds stress, k-ω, and inviscid flow models.In both volutes, measured gross exit angular momentum flux was more than 1.7 times what the zero-torque assumption would predict when blindly applied to the volute as a whole. This discrepancy is attributed to significant turning near the volute’s inlet region leading to an updated view of what an appropriate control volume is when applying the zero-torque assumption. Additionally, variation of both radial and tangential velocity in both the circumferential and axial directions on the order of 15% of the mean value reveal that volute swirl characterization by a single measurement would have a significant associated uncertainty.     

Prediction of Local Losses of Low <Emphasis Type="Italic">Re</Emphasis> Flows in Non-uniform Media Composed of Parrallelpiped Structures

A method is presented to predict the local losses of low Re flow through a porous matrix composed of layers of orthogonally oriented parallelepipeds for which the local geometry varies discreetly in the direction of bulk flow. In each layer, the variations in the pore lengths perpendicular to and parallel to the direction of bulk flow are restricted to be proportional to one another so that the variation in the geometry of each layer may be characterized by a single parameter, \(\beta \). The solutions to the Navier–Stokes equations are determined for flows through geometries that vary in a forward expansion about this parameter. These provide the data used in the development of a correlation that is able to directly relate local hydraulic permeability to the variation in local pore geometry. In this way, the local pressure losses (as well as the relationship between the volumetric flow rate and the total pressure drop) may be determined without requiring the explicit solution of the entire flow field. Test cases are presented showing that the correlation predicts the local pressure losses to be within 0.5% of the losses determined from the numerical solution to the Navier–Stokes equations. When the magnitude of the variation to the geometry is such that the change in the parameter \(\beta \) between layers is constant throughout the medium, a reduced form of the correlation (requiring the evaluation of only three constants) is able to provide predictions of flow rate and interface pressures that agree to within about 1% with the results of the numerical solutions to the Navier–Stokes equations.     

Optimal numerical flux of power-law fluids in some partially full pipes

Consider the steady state pressure driven flow of a power-law fluid in a partially filled straight pipe. It is known that an increase in flux can be achieved for a fixed pressure by partially filling the pipe and having the remaining volume either void or filled with a less viscous, lubricating fluid. If the pipe has circular cross section, the fluid level which maximizes flux is the level which avoids contact with exactly 25% of the boundary. This result can be proved analytically for Newtonian fluids and has been verified numerically for certain non-Newtonian models.

This paper provides a generalization of this work numerically to pipes with non-circular cross sections which are partially full with a power-law fluid. A simple and physically plausible geometric condition is presented which can be used to approximate the fluid level that maximizes flux in a wide range of pipe geometries. Additional increases in flux for a given pressure can be obtained by changing the shape of the pipe but leaving the perimeter fixed. This computational analysis of flux as a function of both fluid level and pipe geometry has not been considered to our knowledge.

Fluxes are computed using a special discretization scheme, designed to uncover general properties which are only dependent on fluid level and/or pipe cross-sectional geometry. Computations use finite elements and take advantage of the variational structure inherent in the power-law model. A minimization technique for approximating the critical points of the associated non-linear energy functional is used. In particular, the numerical scheme for the non-linear partial differential equation has been proved to be convergent with known error estimates. The numerical results obtained in this work can be useful for designing pipes and canals for transportation of non-Newtonian fluids, such as those in chemical engineering and food processing engineering.     

Series solutions for steady unsaturated flow in irregular porous domains

For the most part, analytical solutions for steady unsaturated infiltration have been restricted to infinite and semi-infinite seepage geometries, using the quasi-linear approximation for the hydraulic conductivity. We provide analytical series methods to solve the steady quasi-linear flow equations, in finite irregular seepage geometries. Unlike the classical approach, the series method has been modified, to cater for arbitrary boundary geometry and surface recharge distributions. The matrix flux potential and the stream function both satisfy the same governing partial differential equation, and the stream function formulation is used to estimate the series coefficients. For a finite vadose zone, the stream function solution does not uniquely determine the matrix flux potential, when flux boundary conditions are used. Consequently, the stream function solution applies to a range of moisture distributions, for given infiltration and evapotranspiration rates through the surface.     

On high‐resolution schemes for solving unsteady compressible two‐phase dilute viscous flows

A high‐resolution numerical scheme based on the MUSCL–Hancock approach is developed to solve unsteady compressible two‐phase dilute viscous flow. Numerical considerations for the development of the scheme are provided. Several solvers for the Godunov fluxes are tested and the results lead to the choice of an exact Riemann solver adapted for both gaseous and dispersed phases. The accuracy of the scheme is proven step by step through specific test cases. These simulations are for one‐phase viscous flows over a flat plate in subsonic and supersonic regimes, unsteady flows in a low‐pressure shock tube, two‐phase dilute viscous flows over a flat plate and, finally, two‐phase unsteady viscous flows in a shock tube. The results are compared with well‐established analytical and numerical solutions and very good agreement is achieved. Copyright © 1999 John Wiley & Sons, Ltd.     

Simulation of multiple shock–shock interference using implicit anti‐diffusive WENO schemes

Accurate computations of two‐dimensional turbulent hypersonic shock–shock interactions that arise when single and dual shocks impinge on the bow shock in front of a cylinder are presented. The simulation methods used are a class of lower–upper symmetric‐Gauss–Seidel implicit anti‐diffusive weighted essentially non‐oscillatory (WENO) schemes for solving the compressible Navier–Stokes equations with Spalart–Allmaras one‐equation turbulence model. A numerical flux of WENO scheme with anti‐diffusive flux correction is adopted, which consists of first‐order and high‐order fluxes and allows for a more flexible choice of first‐order dissipative methods. Experimental flow fields of type IV shock–shock interactions with single and dual incident shocks by Wieting are computed. By using the WENO scheme with anti‐diffusive flux corrections, the present solution indicates that good accuracy is maintained and contact discontinuities are sharpened markedly as compared with the original WENO schemes on the same meshes. Computed surface pressure distribution and heat transfer rate are also compared with experimental data and other computational results and good agreement is found. Copyright © 2009 John Wiley & Sons, Ltd.     

A NEW FLUX-CONSERVING NEWTON SCHEME FOR THE NAVIER-STOKES EQUATIONS

SUMMARY

A new numerical method is developed for the two-dimensional, steady Navier-Stokes equations. Using local polynomial expansions to represent the discrete primitive variables on each cell, we construct a scheme which has the following properties: First, the local discrete primitive variables are functional solutions of both the integral and differential forms of the Navier-Stokes equations. Second, fluxes are balanced across cell interfaces using exact functional expressions (to the order of accuracy of the local expansions). No interpolation, flux models, or flux limiters are required. Third, local and global conservation of mass, momentum, and energy are explicitly provided for. Finally, the discrete primitive variables and their derivatives are treated in a unified and consistent manner. All are treated as unknowns to be solved together for simulating the local and global flux conservation.

A general third-order formulation for the steady, compressible Navier-Stokes equations is presented. As a special case, the formulation is applied to incompressible flow, and a Newton's method scheme is developed for the solution of laminar channel flow. H is shown that, at Reynolds numbers of 100, 1000, and 2000, the developing channel flow boundary layer can be accurately resolved using as few as six to ten cells per channel width.     

气体动理学格式研究进展

介绍了近年来气体动理学格式(gas-kinetic scheme, GKS, 亦简称BGK 格式) 的主要研究进展, 重点是高阶精度动理学格式及适合从连续流到稀薄流全流域的统一动理学格式. 通过对速度分布函数的高阶展开和对初值的高阶重构, 构造了时间和空间均为三阶精度的气体动理学格式. 研究表明, 相比于传统的基于Riemann 解的高阶格式, 新格式不仅考虑了网格单元界面上物理量的高阶重构, 而且在初始场的演化阶段耦合了流体的对流和黏性扩散, 也能够保证解的高阶精度. 该研究为高精度计算流体力学(computatial uiddymamics, CFD) 格式的建立提供了一条新的途径. 通过分子离散速度空间直接求解Boltzmann 模型方程,在每个时间步长内将宏观量的更新和微观气体分布函数的更新紧密地耦合在一起, 建立了适合任意Knudsen(kn) 数的统一格式, 相比于已有的直接离散格式具有更高的求解效率. 最后, 本文还讨论了合理的物理模型对数值方法的重要性. 气体动理学方法的良好性能来自于Boltzmann 模型方程对计算网格单元界面上初始间断的时间演化的准确描述. 气体自由运动与碰撞过程的耦合是十分必要的. 通过分析数值激波层内的耗散机制,我们认识到采用Euler 方程的精确Riemann 解作为现代可压缩CFD 方法的基础具有根本的缺陷, 高马赫数下的激波失稳现象不可避免. 气体动理学格式为构造数值激波结构提供了一个重要的可供参考的物理机制.      

A shock-capturing scheme for body-fitted meshes

A finite difference scheme based on flux difference splitting is presented for the solution of the two-dimensional Euler equations of gas dynamics in a generalized co-ordinate system. The scheme is based on numerical characteristic decomposition and solves locally linearized Riemann problems using upwind differencing. The decomposition is for a generalized co-ordinate system and a convex equation of state. This ensures good shock-capturing properties when incorporated with operator splitting and the advantage of using body-fitted co-ordinates. The resulting scheme is applied to supersonic flow of real air' past a circular cylinder.     

Research of influence of reduced-order boundary on accuracy and solution of interior points

The flow field with a high order scheme is usually calculated so as to solve complex flow problems and describe the flow structure accurately.However,there are two problems,i.e.,the reduced-order boundary is inevitable and the order of the scheme at the discontinuous shock wave contained in the flow field as the supersonic flow field is low.It is questionable whether the reduced-order boundary and the low-order scheme at the shock wave have an effect on the numerical solution and accuracy of the flow field inside.In this paper,according to the actual situation of the direct numerical simulation of the flow field,two model equations with the exact solutions are solved,which are steady and unsteady,respectively,to study the question with a high order scheme at the interior of the domain and the reduced-order method at the boundary and center of the domain.Comparing with the exact solutions,it is found that the effect of reduced-order exists and cannot be ignored.In addition,the other two model equations with the exact solutions,which are often used in fluid mechanics,are also studied with the same process for the reduced-order problem.     

Thermosolutal mixed convection and flow-reversal in an inclined parallel-plate channel

The steady-state laminar mixed convection of a binary gas mixture in a parallel-plate channel is investigated. The channel walls are subjected to different combinations of first-type thermal and solutal boundary conditions and different wall inclinations have been considered. A second-order accurate control-volume based numerical scheme is used for the resolution. In parallel with the numerical investigation, the governing conservation equations are also simplified for fully developed conditions and are shown to be controlled by a single parameter. An exact analytical solution is obtained for the main flow variables and transfer rates and serves as a validation tool for the numerical model. In addition, it establishes a criterion based on the two Grashof numbers, the Reynolds number and the channel inclination for the existence of a reversed flow.     

Flux difference splitting for the Euler equations in one spatial co-ordinate with area variation

An approximate (linearized) Riemann solver is presented for the solution of the Euler equations of gas dynamics in one spatial co-ordinate. This includes cylindrically and spherically symmetric geometries and also applies to narrow ducts with area variation. The method is Roe's flux difference splitting with a technique for dealing with source terms. The results of two problems, a spherically divergent infinite shock and a converging cylindrical shock, are presented. The numerical results compare favourably with those of Noh's recent survey and also with those of Ben-Artzi and Falcovitz using a more complicated Riemann solver.