水库异重流一维水沙耦合模型

建立了基于库区不规则断面的一维非恒定异重流数学模型,并采用明流与异重流水沙输移模型交替运算的两步模式,即用潜入条件动态判别异重流计算的上游边界位置,将潜入点上游的明流浑水段与下游异重流段计算连接起来。水流运动、泥沙输移与河床变形过程完全耦合,采用TVD(Total Variation Diminishing)形式的MUSCL-Hancock格式进行数值求解。将该模型应用于恒定流量与释放定量悬沙两种条件下的异重流水槽实验模拟,比较了有无水面梯度项对模拟精度的影响,计算结果表明该模型能较为准确地预测异重流的厚度、含沙量分布及传播过程。     

拟焦沙模拟低含沙量异重流运动初步分析

低含沙量异重流带来大量细颗粒浑水,由于颗粒沉速极小,在山区、高原地区的引调水、水环境、水处理等方面带来难以估量的困难,引起了社会和设计部门重视。以金沙江支流牛栏江上德泽水库回水变动区上下游河段为原型,利用概化水槽试验研究拟焦沙模拟低含沙量异重流形成与运动,探究低含沙浑水异重流运行时头部流速及含沙量沿程变化特征。结果表明:头部流速及含沙量沿程逐渐减小;相同流量不同含沙量所形成的中层、底层异重流头部流速,含沙量越大,同一位置头部流速越大;相同流量和含沙量所形成的表层、中层、底层异重流头部流速,表层异重流头部流速最大,底层异重流次之,中层异重流最小。     

各向异性密度流模型及其验证

针对密度流流动各向异性的特性,采用非线性k-ε模式求取雷诺应力,结合RNG模型及浮力修正建立了三维密度流模型.应用该三维模型计算了水库温差异重流,分析了温差异重流演进过程,并就流速分布与水温下泄过程与实测资料进行了比较.结果表明,所建模型由于引入非线性RNG模式,能有效反映密度流的各向异性特性,较好地模拟了流场和温度场强烈耦合的水库温差异重流,能准确预测温度分层过程与异重流潜行时间,具有较高精度.     

泥沙异重流与环境物质交换经验式对比

采用异重流层平均水沙耦合数学模型,模拟开闸式和恒定入流式泥沙异重流水槽实验,对比分析异重流与环境之间物质交换经验式的适应性和不确定性。考虑4个水卷吸经验式:e_w59、e_w86、e_w87和e_w01,5个泥沙侵蚀经验式:E_s77、E_s86、E_s87、E_s93和E_s04。数值研究表明:水卷吸对于水槽异重流影响较小,应用综合考虑底床摩擦和剪切不稳定的e_w经验式时模拟结果较好;开闸式异重流对床面侵蚀能力有限。对于恒定入流式异重流,应用E_s87和E_s93侵蚀经验式计算所得淤积厚度与实测值吻合较好,可能是率定时综合考虑了异重流实验数据。     

水库浑水异重流潜入点判别条件

泥沙淤积是影响多沙河流水库寿命的一大难题,而异重流排沙是减少库区淤积的重要措施之一。异重流的潜入现象是异重流开始形成的直观标志,研究异重流潜入条件的判别方法有助于掌握异重流在库区内的演进规律。总结了水库异重流潜入条件的定性描述及定量计算方法,指出已有的潜入点判别公式的优缺点及适用范围,改进了描述异重流运动的动量方程,同时分析了异重流流速与含沙量沿垂线不均匀分布对动量传递的影响;在此基础上提出新的异重流潜入条件判别式,并用多组室内及野外实测资料对该判别条件进行率定与验证。分析结果表明,新的计算公式可用于判别小浪底库区异重流的潜入条件。     

辐流式二沉池中异重流的计算流体力学模型研究

异重流的形成是影响活性污泥系统二沉池水力学性能的主要因素,形成异重流的主要原因包括沉淀池内不均匀的污泥分布和温度分布。采用多相流欧拉模型、标准k-ε湍流模型和Boussinesq假设研究了二沉池中活性污泥和表面散热引起的异重流。表面散热的计算采用经验方程,数值计算结果同实验结果作了比较,结果表明模型能较好地预测辐流式二沉池的流场、温度场和污泥分布。研究还表明,Fr_p2Gr/Re2表示了不均匀温度分布引起的浮升力与不均匀污泥分布引起的作用力的比值,用它可以判断两种作用力在异重流的形成当中的相对大小。     

粘性土壤溶质运移新模型的应用

将Moldrup等提出的模拟水流和溶质运移的运动平均斜率模型和运动浓度斜率模型联合,模拟了重粘土在入渗和蒸发条件下的水盐动态并与特征-有限差分数值模型的计算结果进行了比较。在模拟过程中新模型显示出更加高效的优点,这对粘性土壤地区盐渍化趋势的长期预报具有重要的实际意义。     

冲绳海槽西南端1.3ka以来异重流沉积记录及其古气候响应

以冲绳海槽西南端的HOBAB4?S1岩芯为研究对象,综合利用粒度、沉积构造、AMS14C测年等资料,对研究区异重流沉积特征进行研究,对异重流发育的时间段与晚全新世气候事件进行对应。结果表明,岩芯中发育17段异重流沉积,其沉积物组分以砂质粉砂为主,粒度频率分布曲线多呈以70~130 μm为中心的单峰,C?M图上样品点集中分布区间大致平行于C=M基线,且位于PQ段以下,表明沉积物搬运方式为重力流悬浮搬运。异重流内部发育平行层理、爬升沙纹层理和粒序层理等沉积构造。异重流沉积类型主要有两类,一类是厚层异重流沉积,底部侵蚀面发育,内部发育多组逆—正粒序组合,指示了水动力较强,侵蚀作用和沉积作用均较明显的异重流近端沉积;另一类是薄层异重流沉积,底部侵蚀面不发育,内部不发育或仅发育一组逆—正粒序组合,指示了水动力较弱的异重流边部沉积。HOBAB4?S1岩芯的异重流层段集中发育在800~1 300 A.D.之间,指示了当时气候条件为高温高湿,台风、洪水频发,降雨量较大,验证了“中世纪暖期”在东亚地区的存在。     

涠西南凹陷陡坡带流一段上亚段异重流沉积新发现

以涠西南凹陷北部陡坡带中段流一段下亚段为研究对象,综合利用三维地震、测井、录井、钻井岩心及分析测试等资料,对研究区异重流成因深水扇的沉积特征、控制因素及沉积模式进行了研究。结果表明,断陷盆地陡坡带至湖盆底部可发育洪水型异重流沉积,依次为顺直水道、弯曲水道和朵叶体,其延伸距离约为15 km,宽度为50~150 m。其中,水道地震反射明显强于周围深水泥岩,其内部沉积以砂岩和含砾砂岩为主;朵叶体沉积主要为细砂岩或粉砂岩,富含陆源植物碎屑,其粒度概率曲线以上拱弧式及宽缓上拱式为主,样品点分布大致平行于C=M基线,指示重力流沉积动力特征。垂向上,单层砂岩一般表现为逆粒序与正粒序成对出现,中部粒度最粗,可见泥质碎屑呈叠瓦状排列,且发育层内微侵蚀面。此外,异重流沉积内部发育块状层理、平行层理、波状交错层理及爬升砂纹层理等沉积构造。“源汇”系统决定了异重流的沉积特征与展布,万山隆起东部岩浆岩等物源输入是异重流形成的物质基础,强构造运动、陡坡折带、湿润半湿润气候及低湖水密度是控制异重流发育程度的主要因素。研究区异重流的发现不仅丰富了断陷湖盆陡坡带沉积相类型,也为陡坡带进一步油气勘探提供了理论依据。     

倾倒式岩质崩塌运动过程数值模拟分析

在延长县崩塌详细调查的基础上,采用数值模拟方法对倾倒式岩质崩塌运动过程进行了模拟,分析了其运动规律。根据崩塌的破坏形式和运动特点,将倾倒式崩塌的整个运动过程分为四个阶段：岩体原有节理的开裂变形阶段、岩体沿基座支点的倾倒运动阶段、崩塌体在一定初速度下的落体运动阶段、崩塌体的碰撞与堆积阶段。通过对倾倒式崩塌影响范围的理论计算结果、数值模拟结果和实际调查结果对比得出,数值模拟结果和实际调查结果基本一致,两者远大于理论计算结果,因此数值模拟结果可作为倾倒式崩塌的影响范围,可为崩塌灾害的防治提供科学依据。     

河道复杂采砂坑附近流场的数值模拟

数值模拟研究紊流问题的关键是充分考虑雷诺应力的差异,应用各向异性三维代数应力紊流数学模型(ASM)模拟研究了天然河道的次生流问题。以方形管道为算例与试验的比较表明,所应用的ASM模型精度较高;将ASM与非结构网格有机结合模拟分析河道复杂采砂坑,计算成果与前人结果类比,令人满意。ASM模型为研究河床稳定与变形提供了有力工具。分析表明,采砂坑使原稳定的水流形态发生变化,沿主流方向形成纵向涡旋导致采砂坑上游缘口冲刷,坑内的横向次生流则造成横向侵蚀;紊动特性研究表明采砂坑严重影响河床稳定。     

A corrected 3-D SPH method for breaking tsunami wave modelling

A 3-D large eddy simulation model that was first transformed to smoothed particle hydrodynamics (LES-SPH)-based model was employed to study breaking tsunami waves in this paper. LES-SPH is a gridless (or mesh-free), purely Lagrangian particle approach which is capable of tracking the free surface of violent deformation with fragmentation in an easy and accurate way. The Smagorinsky closure model is used to simulate the turbulence due to its simplicity and effectiveness. The Sub-Particle Scale scheme, plus the link-list algorithm, is applied to reduce the demand of computational power. The computational results show that the 3-D LES-SPH model can capture well the breaking wave characteristics. Spatial evolution features of breaking wave are presented and visualized. The detailed mechanisms of turbulence transport and vorticity dynamics are demonstrated as well. This application also presents an example to validate the SPH model.     

振荡流底层悬沙运动的数值研究

建立了平底振荡流底层立面二维水沙数值模型,利用Smagrionsky(SGS)格子涡模型封闭二维Navier Storkes方程水流运动方程组,控制方程采用SMAC法求解。该模型能较精确地模拟振荡流底层水流流动特性,以及含沙量沿垂线分布和随相位变化的情况,且与水槽实验的实测资料基本吻合。     

Fractal and geostatistical methods for modeling of a fracture network

The modeling of fracture networks is useful for fluid flow and rock mechanics studies. About 6600 fracture traces were recorded on drifts of a uranium mine in a granite massif. The traces have an extension of 0.20–20 m. The network was studied by fractal and by geostatistical methods but can be considered neither as a fractal with a constant dimension nor a set of purely randomly located fractures. Two kinds of generalization of conventional models can still provide more flexibility for the characterization of the network: (a) a nonscaling fractal model with variable similarity dimension (for a 2-D network of traces, the dimension varying from 2 for the 10-m scale to 1 for the centimeter scale, (b) a parent-daughter model with a regionalized density; the geostatistical study allows a 3-D model to be established where: fractures are assumed to be discs; fractures are grouped in clusters or swarms; and fracturation density is regionalized (with two ranges at about 30 and 300 m). The fractal model is easy to fit and to simulate along a line, but 2-D and 3-D simulations are more difficult. The geostatistical model is more complex, but easy to simulate, even in 3-D.This paper was presented at Emerging Concepts, MGUS-87 Conference, Redwood City, California, 13–15 April 1987.     

波流边界层水沙运动数值模拟——I:水动力模拟

波流边界层水动力模拟对研究波流相互作用和泥沙运动具有重要的理论意义和实践价值。开发了波流边界层1DV垂向一维水动力数值模型,可用于模拟漩涡沙波床面和平底床面水动力特征。模型的构建基于边界层控制方程,平底床面采用k-ε模型,沙波床面采用双层模型,提出了漩涡层和紊动扩散层交界面紊动动能和紊动耗散率表达式。试验资料验证表明,模型较好地模拟了波浪-水流-床面共同作用下的边界层水动力特征,包括波周期内不同相位流速分布、紊动动能、剪切应力等以及波致时均流速分布和波流相互作用下的时均流速分布等。根据所建模型,讨论了不同床面和波流组合条件下的水动力特征。该模型可为研究波流边界层内水动力特征提供工具。     

A modeling study of the Satilla River estuary,Georgia. I: Flooding-drying process and water exchange over the salt marsh-estuary-shelf complex

The flooding-drying process over the intertidal zone of the Satilla River estuary of Georgia was examined using a three-dimensional (3-D) primitive equations numerical model with Mellor and Yamada's (1982) level 2.5 turbulent closure scheme. The model was forced by the semi-diurnal M₂, S₂, and N₂ tides and freshwater discharge at the upstream end of the estuary. The intertidal salt marsh was treated using a 3-D wet-dry point treatment technique that was developed for the σ-coordinate transformation estuary model. Good agreement was found between model-data comparison at anchor monitoring sites and also along the estuary that suggested that the model provided a reasonable simulation of the temporal and spatial distribution of the 3-D tidal current and salinity in the Satilla River estuary. Numerical experiments have shown that the flooding-drying process plays a key role in the simulation of tidal currents in the main river channel and in water transport over the estuarine-salt marsh complex. Ignoring this process could lead to a 50% under-estimation of the amplitude of tidal currents. The model results also revealed a complex spatial structure of the residual flow in the main channel of the river, with characteristics of multiple eddy-like cell circulations. These complicated residual currents are formed due to tidal rectification over variable topography with superimposition of inertial effects, asymmetry of tidal currents, and baroclinic pressure gradients. Water exchanges over the estuary-intertidal salt marsh complex are asymmetric across the estuary, and tend to vary periodically on the northern side while quickly washing out of the marsh zone on the southern side. Strong Stokes’ drifting velocity was predicted in the estuary, so that the Lagrangian trajectories of particles were characterized by strong nonlinear processes that differ significantly from those estimated by the Eulerian residual currents.     

长江口斜压诊断模式三维流场数值模拟

建立了σ坐标系下长江口斜压诊断模式三维流场数学模型,采用k-kl二方程紊流闭合模型求解垂向涡粘系数,计算域内恒定、非均匀盐度场反映了计算域密度斜压效应。验证结果表明,斜压诊断模式的模拟成果较切合实际。模型成功复演了长江口水域由往复流向旋转流过渡的流场特性;潮波从四条汊道传入后在分汊口附近相遇形成多个会潮点,其位置随着径流的大小、潮汐的强弱、汊道的形态等因素而变化。     

河型转化机理及其数值模拟——Ⅰ.模型建立

为研究河型转化过程机理,建立了考虑弯道二次流影响与边岸崩塌过程的平面二维河流数学模型,包括水流模型、泥沙模型和边岸崩塌模型。通过在水流动量守恒方程中增加弥散应力项以考虑弯道二次流的影响,并采用室内水槽实验结果对水流模型进行了验证;利用上荆江沙市至石首天然长河段的水沙过程和河道演变资料,对泥沙模型进行了验证;结果表明本模型数值计算量合适,有较好的适应范围。模型中提出了边岸崩塌过程的模拟技术,相对于传统平面二维水沙模型而言,可以更好地模拟天然河道的横向摆动以及洲滩消长过程。     

Dynamics of settling-driven convection beneath a sediment-laden buoyant overflow: Implications for the length-scale of deposition in lakes and the coastal ocean

A series of laboratory experiments was conducted in order to determine how settling-driven convection influences the length-scale over which the majority of particles settle beneath a buoyant sediment-laden plume spreading over a denser saline layer. This system is analogous to sediment-laden river water spreading into a lake or the coastal ocean. The key dimensionless parameter that controls the settling dynamics of such flows is the density ratio, defined as the ratio of density differences due to the added salt and sediment. For a buoyant plume, this ratio has to be greater than unity, so that the experiments in the current study were performed for density ratios between one and five. When density ratio was close to one, settling-driven convection was vigorous and the length-scale over which sedimentation occurred was very small. A strong secondary turbidity current was generated in this case. On the other hand, for larger values of density ratio, the predicted length-scale over which a secondary plume was generated increased in proportion to the density ratio. A complete mathematical expression for this length-scale was developed using recent theory that described the timescale over which settling-driven convection evolved. The theoretically predicted propagation length-scale showed very good quantitative agreement with laboratory experiments. The use of the dimensionless density ratio allows the expression to predict which sediment-laden river plumes in lakes and the coastal ocean could quickly form secondary turbidity currents.     

近岸波生流运动三维数值模拟及验证

开发建立了近岸波生流运动三维数值计算模式。模式中,引入了三维时均剩余动量、破波表面水滚、波浪水平与垂向紊动作为主要驱动力,同时考虑了波流共同作用的底部剪切力。推导了可综合反映底坡、能量传递率和密度影响的水滚能量传输方程;将Larson-Kraus的二维波浪水平紊动系数表达式拓展至三维。采用大量实测数据和文献资料测试验证了所建模式,表明所建模式可有效模拟波浪增减水、底部离岸流、沿岸流、裂流、堤后环流等不同维度的波生流现象。此外,研究也表明破波水滚效应可解释波生流峰值向岸推移的物理现象,从而在模拟中不能忽略;破波带内沿岸流速垂向较为均匀的现象与波浪附加垂向紊动有关。