三维复杂流动的间断有限元方法模拟

本文基于三维可压缩Euler方程,采用基于Runge-Kutta时间离散的间断有限元方法(RKDG方法),对三维前台阶、三维Riemann问题和球Riemann等问题进行了模拟。结果表明,本文的RKDG方法能够在很少的网格内清晰地捕捉到三维复杂流场中的激波和接触间断;同时,将球Riemann问题中z=0.4平面压强沿到对称轴距离的分布与文献中的近似精确解相比,吻合较好,这也验证了本文的RKDG方法不仅能够进行三维复杂流场的定性描述,也能够应用于三维复杂流场的定量计算。     

Theoretical analysis of non-uniform extensional flows in relation to processing rheology and predictions of some constitutive equations

In this paper, a characteristic equation involving the stream function, already given by one of the authors in a previous work for classifying axisymmetric incompressible flows, is re-considered. Non-uniform nearly extensional flows are derived as particular solutions from this equation. Using experimental data in the literature for polymer solutions and melts, it is proved that particular solutions of the characteristic equation lead to kinematics very close to those encountered in the fiber-spinning process. The kinematic equations satisfactorily correlating the fiber-spinning data are used in order to determine the ability of constitutive equations to predict realistic stresses in the flow domain. The rheological parameters of the fluids, obtained from experiments, are used for computation of differential and integral constitutive equations in the spinning conditions. Comparisons with the stress response of adequate constitutive equations are given and discussed.Also affiliated to: Université Joseph Fourier Grenoble I and Institut National Polytechnique de Grenoble, Associé au CNRS (URA 1510)     

Analysis of plane and axisymmetric flows of incompressible fluids with the stream tube method: Numerical simulation by trust-region optimization algorithm

New concepts for the study of incompressible plane or axisymmetric flows are analysed by the stream tube method. Flows without eddies and pure vortex flows are considered in a transformed domain where the mapped streamlines are rectilinear or circular. The transformation between the physical domain and the computational domain is an unknown of the problem. In order to solve the non-linear set of relevant equations, we present a new algorithm based on a trust region technique which is effective for non-convex optimization problems. Experimental results show that the new algorithm is more robust compared to the Newton-Raphson method.     

Numerical simulation of the effect of void fraction and inlet velocity on two-phase turbulence in bubble-liquid flows

There are contradicted opinions on whether bubbles enhance or reduce the liquid turbulence. In this paper, the effect of void fraction and inlet velocity on the bubble–liquid two-phase turbulence of the multiple bubble–liquid jets in a two-dimensional channel is studied by using the two-phase second-order moment turbulence model. The results confirm the phenomena observed in experiments and reported in references that at a low void fraction and low inlet velocities the bubbles enhance the liquid turbulence, whereas at a high void fraction and high inlet velocities the bubbles reduce the liquid turbulence.The project supported by the China Special Funds for Major State Basic Research (G-1999-0222-08) and the Innovation and Technology Commission of Hong Kong and Aoyagi (H.K.) Ltd, Hong Kong, under the Grant No. UIM/122. The English text was polished by Keren Wang.     

Numerical simulation of viscous flows with free surface around realistic hull forms with transom

This paper describes a method for simulation of viscous flows with a free surface around realistic hull forms with a transom, which has been developed based on a FINFLO RANS solver with a moving mesh. A dry‐transom model is proposed and implemented for the treatment of flows off the transom. The bulk RANS flow with the artificial compressibility is solved by a cell‐centred finite volume multigrid scheme and the free surface deformed by wave motions is tracked by satisfying the kinematic and dynamic free‐surface boundary conditions on the actual location of the surface. The effects of turbulence on flows are evaluated with the Baldwin–Lomax turbulence model without a wall function. A test case is modern container ship model with a transom, the Hamburg Test Case. The calculated results are validated and they agree well with the measured results in terms of the free‐surface waves and the total resistance coefficient. Furthermore, the numerical solutions successfully captured many important features of the complicated interaction of the free surface with viscous flows around transom stern ships. In addition, the convergence performance and the grid refinement studies are also investigated. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical modelling of 3D stratified free surface flows: a case study of sediment dumping

A three‐dimensional numerical model has been developed to simulate stratified flows with free surfaces. The model is based on the Reynolds‐averaged Navier–Stokes (RANS) equations with variable fluid density. The equations are solved in a transformed σ‐coordinate system with the use of operator‐splitting method (Int. J. Numer. Meth. Fluids 2002; 38 :1045–1068). The numerical model is validated against the one‐dimensional diffusion problem and the two‐dimensional density‐gradient flow. Excellent agreements are obtained between numerical results and analytical solutions. The model is then used to study transport phenomena of dumped sediments into a water body, which has been modelled as a strongly stratified flow. For the two‐dimensional problem, the numerical results compare well with experimental data in terms of mean particle falling velocity and spreading rate of the sediment cloud for both coarse and medium‐size sediments. The model is also employed to study the dumping of sediments in a three‐dimensional environment with the presence of free surface. It is found that during the descending process an annulus‐like cloud is formed for fine sediments whereas a plate‐like cloud for medium‐size sediments. The model is proven to be a good tool to simulate strongly stratified free surface flows. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical simulation of complex flows of non-Newtonian fluids using the stream tube method and memory integral constitutive equations

In this paper a memory integral viscoelastic equation is considered for simulating complex flows of non-Newtonian fluids by stream tube analysis. A formalism is developed to take into account co-deformational memory equations in a mapped computational domain where the transformed streamlines are parallel and straight. The particle-tracking problem is avoided. Evolution in time and related kinematic quantities involved with a K-BKZ integral constitutive model are easily taken into account in evaluating the stresses. Successive subdomains, the stream tubes, may be considered for computing the main flow in abrupt axisymmetric contractions from the wall to the central flow region. The ‘peripheral stream tube’ close to the duct wall is determined by developing a non-conventional modified Hermite element. A mixed formulation is adopted and the relevant non-linear equations are solved numerically by the Levenberg-Marquardt algorithm. Although the singularity at the section of contraction is not involved explicitly, the results obtained for the peripheral stream tube clearly show the singularity effects and the extent of the recirculating zone near the salient corner. The algorithm is stable even at high flow rates and provides satisfactory solutions when compared with similar calculations in the literature.     

Assessment of nonlinear strain measures for extensional and shearing flows of polymer melts

From stress-strain experiments in extensional and shearing flows, nonlinear strain measures and effective damping functions are derived for a polyisobutylene melt. The strain measures determined in planar extensional flow and in simple shear flow coincide. Experimental results are compared with predictions of two molecular theories, the Doi-Edwards model and the molecular stress function approach of Wagner and Schaeffer. Discrepancies between theories and experiment lead to a reconsideration of the classification of extensional flows. The symmetry of the flow field is identified and quantified as an important parameter influencing the strain measure, and a unifying strain measure for general extensional and shearing flows of polymer melts is presented.     

Numerical simulation and hydrodynamic visualization of transient viscous flow around an oscillating aerofoil

Unsteady viscous flow around a large-amplitude and high-frequency oscillating aerofoil is examined in this paper by numerical simulation and experimental visualization. The numerical method is based on the combination of a fourth-order Hermitian finite difference scheme for the stream function equation and a classical second-order scheme to solve the vorticity transport equation. Experiments are carried out by a traditional visualization method using solid tracers suspended in water. The comparison between numerical and experimental results is found to be satisfactory. Time evolutions of the flow structure are presented for Reynolds numbers of 3 × 10³ and 10⁴. The influence of the amplitude and frequency of the oscillating motion on the dynamic stall is analysed.     

The strain-hardening behaviour of linear and long-chain-branched polyolefin melts in extensional flows

By generalizing the Doi-Edwards model to the Molecular Stress Function theory of Wagner and Schaeffer, the extensional viscosities of polyolefin melts in uniaxial, equibiaxial and planar constant strain-rate experiments starting from the isotropic state can be described quantitatively. While the strain hardening of four linear polymer melts (two high-density polyethylenes, a polystyrene and a polypropylene) can be accounted for by a tube diameter that decreases affinely with the average stretch, the two long-chain-branched polymer melts considered (a low-density polyethylene and a long-chain branched polypropylene) show enhanced strain hardening in extensional flows due to the presence of long-chain branches. This can be quantified by a molecular stress function, the square of which is quadratic in the average stretch and which follows from the junction fluctuation theory of Flory. The ultimate magnitude of the strain-hardening effect is governed by a maximum value of the molecular stress, which is specific to the polymer melt considered and which is the only free non-linear parameter of the theory. Received: 1 June 1999/Accepted: 24 November 1999     

Numerical simulation of expansion/compression effect on shock induced turbulent flow

The interaction of homogeneous and isotropic turbulence with a shock wave is observed by solving the Reynolds-averaged Navier–Stokes equations with the k–? turbulence model. All turbulent fluctuations are measured at the period of expansion in the turbulent field and during compression by the reflected shock on turbulent field, and it is observed that the longitudinal turbulent velocity fluctuation is enhanced more at the period of expansion due to incident shock wave movement far from the turbulent field. The amplification of the turbulent kinetic energy (TKE) level in the shock/turbulence interaction depends on the shock wave strength and the longitudinal velocity difference across the shock wave. On decreasing the longitudinal velocity difference across the shock, the turbulent kinetic energy (TKE) level is less amplified. The TKE level is amplified by the factor of 1.5–1.8 in the shock/turbulence interaction where the dissipation rate of TKE decreases in all cases of shock/turbulence interaction. After the shock/turbulence interaction, the turbulent dissipative-length scale is amplified slightly and the amplification of the length scales decreases when increasing the shock strength. To cite this article: M.A. Jinnah, K. Takayama, C. R. Mecanique 333 (2005).     

Numerical simulation of thermally and chemically nonequilibrium flows and heat transfer in underexpanded induction plasmatron jets

The results of numerical simulation are presented for thermally and chemically nonequilibrium air plasma flows in a plasmatron discharge channel and underexpanded dissociated and partially ionized air jets flowing past a cylindrical model with a blunt leading edge and cooled copper surface under the experimental conditions realized in a VGU-4 100 kW induction plasmatron (Institute for Problems in Mechanics of the Russian Academy of Sciences) (see, for example, [1, 2]). The nonequilibrium excitation of the vibrational degrees of freedom of the molecules in the modal approximation and the difference between the electron and translational heavy-particle temperatures are taken into account in the calculations. The calculated data on the heat transfer and pressure at the stagnation point are compared with the results obtained within the framework of the thermally equilibrium model. Comparison with the experimental data obtained in the Institute for Problem in Mechanics of the Russian Academy of Sciences (Laboratory for interaction between plasma and radiation and materials) and kindly provided for comparison purposes gives satisfactory agreement.     

用CE/SE法对弯曲与分叉河道的溃坝洪水波的数值研究

对作者的二维溃坝洪水波的数学模型进一步推广,得到了一般形式的基于任意四边形网格的时空守恒元和解元方法(简称CE/SE法)的新的格式.CE/SE法从守恒积分型浅水方程出发,设立守恒元和解元,严格保证其物理意义上的守恒律,并且构造思想简单,格式通用性好.首先采用CE/SE法计算等宽矩形河道的溃坝洪水波,并与Stoker解析解进行比较,在此基础上,数值模拟了180度强弯曲河道、45度三支分叉河道的二维溃坝洪水波的演进过程,揭示了溃坝洪水波在弯曲河道中内外两岸速度与水位的变化,在分叉河道中自动进行流量与动量的再分配,在分叉点处形成旋涡,水位变化剧烈等复杂的运动特征,算例结果表明基于任意四边形网格的CE/SE法精度高,稳定性好,该格式对各种不规则几何区域内的溃坝问题具有较强的适应性,对溃坝洪水波的间断具有较高的分辨率.     

人群行走激励的数值模拟及动力耦合作用研究

大跨轻柔结构需要考虑人群行走激励下的结构振动响应以及人体的舒适度.根据单步落足曲线,结合影响人群行走效应的因素如人群分布、行走频率及个体生理特征等,采用随机方法模拟人群行走激励.对不同的人群倍增因子计算方法进行了比较并提出改进建议.将人群等效为具有主体部分和耦合均布质量的时变动力模型,建立了梁板结构与移动人体相耦合的动...     

Numerical simulation of axisymmetric turbulent flow in combustors and diffusers

Numerical studies of turbulent flow in an axisymmetric 45° expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k–? model were solved in a non-orthogonal curvillinear co-ordinate system. A zonal grid method, wherein the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. However, treatment of the zonal boundaries required special handling. Boundary overlap and interpolating techniques were used and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite differencing methods—hybrid, quadratic upwind and skew upwind—were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.     

Numerical simulation of turbulent free surface flow with two‐equation k–ε eddy‐viscosity models

This paper presents a finite difference technique for solving incompressible turbulent free surface fluid flow problems. The closure of the time‐averaged Navier–Stokes equations is achieved by using the two‐equation eddy‐viscosity model: the high‐Reynolds k–ε (standard) model, with a time scale proposed by Durbin; and a low‐Reynolds number form of the standard k–ε model, similar to that proposed by Yang and Shih. In order to achieve an accurate discretization of the non‐linear terms, a second/third‐order upwinding technique is adopted. The computational method is validated by applying it to the flat plate boundary layer problem and to impinging jet flows. The method is then applied to a turbulent planar jet flow beneath and parallel to a free surface. Computations show that the high‐Reynolds k–ε model yields favourable predictions both of the zero‐pressure‐gradient turbulent boundary layer on a flat plate and jet impingement flows. However, the results using the low‐Reynolds number form of the k–ε model are somewhat unsatisfactory. Copyright © 2004 John Wiley & Sons, Ltd.