曲线坐标系下平面二维浅水模型的修正与应用

借鉴有关弯道水流流速分布的研究成果,计入弯道环流引起的横向的动量交换,对曲线坐标系平面二维浅水方程做了修正;采用边界处正交的曲线网格生成技术,处理复杂的计算区域边界。采用修正后的模型对90°弯道水流进行了数值模拟,并与原模型的计算结果及实测资料进行了比较,结果表明该模型能够有效地模拟流线弯曲的复杂水流的水力特性。     

弯道水流数值模拟研究

采用边界拟合坐标技术,生成边界处正交的曲线网格,在曲线坐标系中进行弯道水流数值计算。借鉴现有关于弯道水流流速分布的研究成果,对曲线坐标系平面二维浅水方程作了修正,修正后的方程组能够计入弯道环流引起的横向的动量交换,即考虑了二次流对流线弯曲的复杂水力特性的影响。通过连续弯道实例计算,验证了该模型的可靠性和实用性。     

常曲率弯道二维流动稳定性与非线性演化规律分析

国内外对弯曲河流的研究从线性特性转移到不稳定性与非线性特性上来, 弯曲河流作为一个不稳定的系统, 主要是受到水流和边界的不稳定性的影响.针对弯曲河流中存在的水流动力的不稳定性与非线性特征,利用弱非线性理论与摄动法进行展开, 在微弯条件下建立了时间模式下的常曲率 弯道二维水流扰动幅值与扰动幅角的非线性演化方程.研究了在不同弯曲度下的扰动发展特征,探讨了扰动流场在弯曲度影响下的时空分布的规律,分析了扰动流场中出现的扰动漩涡的运动过程, 阐述了弯曲度对弯曲边界内水流不稳定性的影响, 具体表现为在微弯的情况下, 弯道的弯曲度增大会提高河湾内水流的稳定性, 扰动振幅与扰动幅角会随弯曲度的增大其衰减速度更快, 并且扰动流速的对称性在弯曲度较大时减弱, 逐渐向弯道的凸岸偏斜, 在中性状态附近 的弯道水流动力具有对流不稳定性和非线性衰减的特征.本文的研究成果为构建水流动力非线性与床面形态几何非线性相互作用的模型提供了思路.      

弯道水流二维浅水模型的修正研究

借鉴有关弯道水流流速分布的研究成果,计入深度平均流速与真实流流速分布差值引起的扩散效应,在正交曲线坐标系下建立了平面二维浅水模型.采用以标准κ-ε模型为基础的曲率效应修正紊流模型模拟紊动应力项,在一定程度上考虑了流线弯曲水流紊动应力的各向异性.应用控制体积法和交错网格法离散方程,并用SIMPLEC算法求解离散方程;同时采用修正后的模型对90°弯道水流进行了数值模拟,并与原模型的计算结果及实测资料进行了比较,结果表明该模型能够有效地模拟流线弯曲水流的水力特性.     

摆动河槽水动力稳定性特征分析

河流形态与水动力结构息息相关,形态约束水动力结构,水动力结构则通过泥沙运动进一步塑造形态,在自然界河流中形成一对辩证互馈关系.天然河流形态形式多样,大致可以分为顺直、微弯、分叉和散乱游荡几种类型,其中微弯及多个弯曲构成的河型为河流动力演化中最重要的一环.多个弯曲构成的河型可用正弦派生曲线来描述,它也是天然河流主槽与水动力结构复杂相互作用的结果.作为探讨这一过程的力学作用机理,构建摆动槽道并研究槽道摆动与其内部流动结构的互馈关系,既是流体力学研究的热点内容,也是目前河流动力过程研究的基础内容.在此重点讨论这一互馈关系前一部分,即水流对摆动边界的响应.文中建立了随体坐标系下摆动河槽与内部水流动力响应理论模型,通过给定摆动弯曲槽道的不同特征参数,研究讨论了正弦派生型摆动边界下的槽道水流动力稳定性特征,明确了弯曲槽道摆动对其内部主流及扰动水流结构的影响,确定弯曲槽道摆动波数、摆动频率对扰动流发展影响的相应参数定量关系,得到了槽道弯曲度和摆动特征对其内部水流不同尺度扰动影响的阈值选择性范围.     

混凝土面板砂砾石坝漫顶溃坝数值模拟

通过溃坝案例的现场调查和水槽模型试验，建立了面板砂砾石坝漫顶溃坝数学模型，以模拟水流特性、砂砾石料以及混凝土面板三者之间的相互作用。该模型主要由3部分组成：(1)考虑漩涡水流和砂砾石料之间的冲蚀特性，引入RNG k-ε湍流模型模拟漩涡水流的反向侵蚀；(2)通过泥沙输移公式和流体体积法(VOF法)追踪砂粒与水流侵蚀交界面，并考虑了砂砾石料的孔隙特性；(3)基于弯矩平衡法判定混凝土面板在自重和水荷载作用下的破坏过程。通过水槽模型试验对数学模型进行验证，结果表明，建立的数学模型能够准确地模拟溃决过程中漩涡水流对大坝的反向侵蚀和面板溃口的详细演化过程；该模型计算的溃口流量、堆石体溃口的发展、溃坝历时以及面板的折断长度，与实测数据的相对误差均小于15%,验证了模型的合理性。与参数化模型相比，本研究提出的溃坝模拟数学模型具有更详细的结果。     

均匀流中双色波传播特性的模拟研究

针对双色波浪与均匀流相互作用问题,采用时域高阶边界元方法建立自由水面满足完全非线性边界条件的数学模型。求解中采用混合欧拉-拉格朗日方法追踪流体瞬时水面,运用四阶龙格库塔方法更新下一时间步的波面和速度势。通过与已发表试验结果对比,验证了本模型的准确性。通过数值计算研究了水流参数对各组成波及衍生的高阶波幅值、波浪和水流间能量交换的影响规律。     

锥阀附近水力特性的三维数值模拟研究

采用可压缩的两相流模型,辅以Realizable k-ε湍流模型来模拟锥阀附近的水流水力特性,对于自由水面的处理采用了VOF法。通过对计算的压力场以及流线的分析,指出锥阀在水力消能中的作用。数值模拟结果与试验资料的比较分析表明两者吻合较好,从而验证了该模型在锥阀湍流模拟中的可靠性。     

海域水流及水质之综合智慧型系统

本文研究之目的为将近年来先进的人工智能技术与现有的数学模型互相结合，开发一套有关海域水流及水质之综合智能型系统。其应用之技术包括智能型系统、人工神经网络及模糊推理系统。此系统可作辅助设计或训练之工具，能帮助水力及环境工程师熟习最先进的海域水流及水质模型，以期使科研与实际执行者之间的距离缩小。其结果可导致对各种不同方式数学模型之好处、应用或限制，有更透彻的了解，甚至从而达成突破及贡献。随着环境可持续性愈来愈被重视，对开发此类综合系统之需求亦逐渐增加，正好可辅助决策者迅速地达成决定，并同时可向公众提供方便及公开的水质信息。     

页岩气储层压裂数值模拟技术研究进展

页岩气储层水力压裂数值模拟既要考虑页岩储层岩石的特性,又要兼顾水平外分段压裂施工工艺,是一个非常棘手的力学难题.本文简述了页岩气储层岩石具有的地质力学特征和页岩气储层开发常用的水平井分段压裂技术;详述了扩展有限元、边界元、离散元在水力压裂裂缝模拟上的应用现状,指出了它们在处理裂缝问题的局限性和优越性,总结出边界元三维位移不连续法是模拟多裂缝扩展的有效方法.     

Effect of streamline curvature on flow field of a turbulent plane jet in cross-flow

The effects of the streamline curvature and finite edge velocity on the flow field of a turbulent plane jet in cross-flow are studied numerically by incorporating the curvature effect in the k–ε turbulence model. The improvement in the predictions by the streamline curvature model is assessed by comparing its prediction with those by the standard k–ε model. The predictions by both the models are compared with available experimental data. It has been observed that the performance of the k–ε model with streamline curvature modification is superior to the standard k–ε model.     

Effect of Streamline Curvature on Three-Dimensionality of Transitional Near-Wall Flow in a Linear Hydrofoil Cascade: A DNS Investigation

Flow, Turbulence and Combustion - The laminar flow on a curved surface transits to turbulent induced by streamline curvature which generates pressure gradient field and separated shear layer flow....     

Nonlinear turbulence models for predicting strong curvature effects

Prediction of the characteristics of turbulent flows with strong streamline curvature, such as flows in turbomachines, curved channel flows, flows around airfoils and buildings, is of great importance in engineering applications and poses a very practical challenge for turbulence modeling. In this paper, we analyze qualitatively the curvature effects on the structure of turbulence and conduct numerical simulations of a turbulent Uduct flow with a number of turbulence models in order to assess their overall performance. The models evaluated in this work are some typical linear eddy viscosity turbulence models, nonlinear eddy viscosity turbulence models （NLEVM） （quadratic and cubic）, a quadratic explicit algebraic stress model （EASM） and a Reynolds stress model （RSM） developed based on the second-moment closure. Our numerical results show that a cubic NLEVM that performs considerably well in other benchmark turbulent flows, such as the Craft, Launder and Suga model and the Huang and Ma model, is able to capture the major features of the highly curved turbulent U-duct flow, including the damping of turbulence near the convex wall, the enhancement of turbulence near the concave wall, and the subsequent turbulent flow separation. The predictions of the cubic models are quite close to that of the RSM, in relatively good agreement with the experimental data, which suggests that these models may be employed to simulate the turbulent curved flows in engineering applications.     

A hybrid k-ε turbulence model of recirculating flow

This investigation deals with the modification of streamline curvature effects in the k-ε turbulence model for the case of recirculating flows. Based upon an idea that the modification of curvature effects in C₂ should not be made in regions where the streamline curvature is small, a hybrid k-ε model extended from the modification originally proposed by Srinivasan and Mongia is developed. A satisfactory agreement of model predictions with experimental data reveals that the hybrid k-ε model can perform better simulation of recirculating turbulent flows.     

A study of the effects of curvature and compression on the behavior of a supersonic turbulent boundary layer

The influences of concave streamline curvature and bulk compression on a turbulent boundary layer in a supersonic flow were studied. The effects of curvature and compression were identified by comparing two flows with identical pressure distributions but different amounts of curvature. The effect of the pressure rise was to increase the wall stress and amplify the Reynolds stresses. The influence of curvature complemented the destabilizing effect of the compression, and the larger the curvature the greater the increase in the wall friction and Reynolds stresses in the boundary layer.This work was supported by AFOSR grant 90-0261, monitored by Dr. L. Sakell     

Finite volume computation of the turbulent flow over a hill employing 2D or 3D non-orthogonal collocated grid systems

A general numerical method for the solution of the complete Reynolds-averaged Navier-Stokes equations for 2D or 3D flows is described. The method uses non-orthogonal co-ordinates, Cartesian velocity components and a pressure-velocity-coupling algorithm adequate for non-staggered grid systems. The capability of the method and the overall performance of the κ–? eddy viscosity model are demonstrated by calculations of 2D and 3D flow over a hill. Solution error estimations based on fine grids, e.g. 320 × 192 control volumes, together with comparisons with standard turbulence model modifications, low-Reynoldsnumber or streamline curvature effects, have allowed the investigation of model drawbacks in predicting turbulent flows over surface-mounted hills.     

Computation of turbulent asymmetric wake

The development of asymmetric wake behind an aerofoil in turbulent incompressible flow has been computed using finite volume scheme for solving two-dimensional Navier-Stokes equations along with the k-ε model of turbulence. The results are compared with available experimental data. It is observed that the computed shift of the point of minimum velocity with distance is sensitive to the prescribed value of the normal component of velocity at the trailing edge of the aerofoil. Making the model constant C_u as a function of streamline curvature and changing the production term in the equation for ε, has only marginal influence on the results.     

Modeling the Turbulent Vortex Wake Behind a High-Lift Wing

A mathematical model and the results of a numerical simulation of the multiple vortex structure behind the high-lift wing (with deflected flaps) of a civil aircraft model are presented. The calculations are performed within the framework of the simplified (with allowance for the specific features of the flow under consideration) three-dimensional Reynolds equations using three turbulence models. The results obtained are compared with the experimental data; a considerable influence of the streamline curvature, typical of vortex flows, not only on the turbulence parameters but also on the parameters of the average flow is demonstrated. As a result, only one of the turbulence models considered, namely, that taking the streamline curvature and rotation effects on the turbulence into account, ensures acceptable accuracy of the main wake characteristics.     

Modeling rotation and curvature effects within scalar eddy viscosity model framework

Two approaches to incorporate the effects of rotation and curvature in scalar eddy viscosity models are explored. One is the “Modified coefficients approach” – to parameterize the model coefficients such that the growth rate of turbulent kinetic energy is suppressed or enhanced. The other is the “Bifurcation approach” – to parameterize the eddy viscosity coefficient such that the equilibrium solution bifurcates from healthy to decaying solution branches. Simple, yet, predictive models in each of these two approaches are proposed and validated on some benchmark test cases characterized by profound effects of system rotation and/or streamline curvature. The results obtained with both the models are encouraging.