A TVD Scheme for Incompressible Flow Coupled With Different Turbulence Modles on A Ground-Mounted Square-RIB Flow

A finite-difference Total Variation Diminishing (TVD) numerical simulation model for coupling the Reynolds Averaged Navier-Stokes (RANS) equations, pressure-relative continuity equation and various k-ε turbulence models was developed to solve the incompressible flow based on the pseudo-compressibility method. The hyperbolicity of all these equations was studied and the discretization of the fully coupling equations with all the primal variables and source terms were made in this article. Numerical simulation for modeling the flow around a ground-mounted square rib was implemented and validated by comparing with the published wind tunnel experimental data. It is shown that such a numerical simulation method with a proper turbulence model has a very good accuracy to simulate the flow around a surface-mounted rib. It is concluded that the Renormalization Group (RNG) and Chen-Kim k-ε turbulence models have much better ability to predict the characteristics of the vortex structure and flow separation than the standard k-ε model.     

FINITE ELEMENT ANALYSIS OF THE FLOW INDUCED BY ROTATING BLADES IN AN INCOMPRESSIBLE VISCOUS FLUID

Aerodynamic loads on a multi-bladed helicopter rotor in hovering flight were calculated by solving the three-dimensional incompressible Navier-Stokes equations. The rotor wake effects were accounted by the correction of local geometric angle of attack according to a free-wake modeling in addition to an empirical modification for the tip flow effect. The validity and efficiency of the present method were verified by the comparisons between numerical results and experimental data.     

NUMERICAL SIMULATIONS OF 2D PERIODIC UNSTEADY CAVITATING FLOWS

1. INTRODUCTION Cavitation is a natural phenomenon especially existing in liquids. A cavitating flow generally involves a large number of vapor structures such as bubbles or vortices which are convecting downstream. When they reach high pressure zones, th…     

An Operator-integration-factor splitting method for time-dependent problems: Application to incompressible fluid flow

In this paper we present a simple, general methodology for the generation of high-order operator decomposition (splitting) techniques for the solution of time-dependent problems arising in ordinary and partial differential equations. The new approach exploits operator integration factors to reduce multiple-operator equations to an associated series of single-operator initial-value subproblems. Two illustrations of the procedure are presented: the first, a second-order method in time applied to velocity-pressure decoupling in the incompressible Stokes problem; the second, a third-order method in time applied to convection-Stokes decoupling in the incompressible Navier-Stokes equations. Critical open questions are briefly described.     

排桩式波浪能发电装置附近流场特性研究

基于雷诺时均Navier-Stokes方程和k?ω湍流模型,研究单排桩式振荡水柱式波浪能发电装置在规则波作用下的桩基附近流场特性. 通过物理模型实验验证所建数值波浪水槽的准确性. 模拟规则波作用下装置附近的流场特性,分析装置附近的涡特征以及其对装置附近泥沙冲刷情况的潜在影响. 结果表明,测试工况下,在桩式振荡水柱装置桩柱体的桩基附近观察到马蹄涡和尾涡现象,马蹄涡强度随着Keulegan–Carpenter（KC）数的增大而减小,尾涡强度随着KC数的增大而增大. 马蹄涡因强度过小,在模拟的KC数范围内不是装置桩基泥沙冲刷的主要因素,高强度的尾涡很可能成为装置桩基附近泥沙起动、输运和冲刷的重要影响因素.     

Dynamics of fluid circulation in coupled free and heterogeneous porous domains

Fluid flow in coupled free and porous domain, particularly when the porous medium is heterogeneous, is encountered in many hydro-environmental conditions, e.g., leakage from underground pipe, combined groundwater lake-subsurface interactions. One of the most difficult problems in the study of coupled flow behaviour has been the development of a universally applicable modelling scheme for combining the flow regimes. This is because the free/porous interfacial properties (e.g., shear-stress; velocity slip) that govern the coupled flow behaviour are difficult to determine experimentally under hydro-environmental conditions. On the other hand, the implications of various forms of heterogeneity in the porous media properties can be very different on the fluid-flow behaviour. Difficulties may also arise in direct coupling of the model equations that govern the fluid flow in the individual regions (e.g., Navier-Stokes for free-flow region and the Darcy's equation for the porous flow region). Consequently, models of coupled free and porous flow for hydro-environmental conditions are not very well developed at the moment. While there are some indications that fluids in coupled free and porous domains may circulate (i.e., development of flow cells), there is a lack of appropriate 3D analysis on how heterogeneities in porous media may affect such flow patterns. In this paper, we aim to analyse how porous media heterogeneity affects the dynamics of flow circulation in the porous side of a coupled free and porous domain. For this purpose, we analyse flow patterns in several model domains made up of two porous layers with differing permeabilities. The governing model equations are discretised and solved using the standard finite volume method on a staggered cell-centred mesh. The temporal discretisation is done using the explicit method. An in-house graphical user interface (GUI) has been created specifically to aid in the visualisation of otherwise complex flow patterns. The GUI contains many post-processing options and provides a comprehensive tool for the analysis of hydrodynamics and contaminant motion (not discussed in this paper) in coupled free and porous flow domains. This GUI is described in this paper briefly. The effects of altering the aspect ratio (i.e., multi-scale) of the domain on the coupled flow pattern have also been discussed.     

Shape optimization of a micromixer with staggered herringbone groove

Shape optimization of a staggered herringbone groove micromixer using three-dimensional Navier-Stokes analysis has been carried out. The analysis of the degree of mixing is performed by the calculation of spatial data statistics. The calculation of the variance of the mass fraction at various nodes on a plane in the channel is used to quantify mixing. A numerical optimization technique is applied to optimize the shape of the grooves on a single wall of the channel. Two design variables, namely, the ratio of the groove depth to channel height ratio and the angle of the groove, are selected for optimization. A mixing index is used as the objective function. The results of the optimization show that the mixing is very sensitive to the shape of the groove which can be used in controlling the mixing in microdevices. The mixing is affected by the depth of the groove much more than the angle of the groove.     

平面上不可压缩的Navier-Stokes方程五模截断方程组的标准型

把广义的类Lorenz系统推广到五维,给出了平面上不可压缩的Navier-Stokes方程五模截断方程组广义的类Lorenz系统的标准型;讨论了截断方程组标准型的形式统一性。     

基于Smoothed Particle Hydrodynamics方法的实时流体模拟

提出基于物理的、实时的技术模拟水体,根据Smoothed Particle Hydrodynamics方法求解流体动力学的Navier-Stokes方程,并运用Marching Cubes体绘制算法重建流体表面.实验表明基于粒子的方法能模拟流体所特有的飞溅、漩涡等现象,Marching Cubes在表面重建的高效性使模拟达到实时、交互的应用.     

A class of high-resolution algorithms for incompressible flows

We present a class of a high-resolution Godunov-type algorithms for solving flow problems governed by the incompressible Navier-Stokes equations. The algorithms use high-resolution finite volume methods developed in LeVeque (SIAM J Numer Anal 1996;33:627-665) for the advective terms and finite difference methods for the diffusion and the Poisson pressure equation. The high-resolution algorithm advects the cell-centered velocities using the divergence-free cell-edge velocities. The resulting cell-centered velocity is then updated by the solution of the Poisson equation. The algorithms are proven to be robust for constant-density flows at high Reynolds numbers via an example of lid-driven cavity flow. With a slight modification for the projection operator in the constant-density solvers, the algorithms also solve incompressible flows with finite-amplitude density variation. The strength of such algorithms is illustrated through problems like Rayleigh-Taylor instability and the Boussinesq equations for Rayleigh-Bénard convection. Numerical studies of the convergence and order of accuracy for the velocity field are provided. While simulations for two-dimensional regular-geometry problems are presented in this study, in principle, extension of the algorithms to three dimensions with complex geometry is feasible.