二维浅水流动的一种普适的高性能格式——有限体积Osher格式

天然水体具有形状复杂的计算区域和水下地形。本文采用无结构的网格以适合这一情况,并能方便地根据工程应用的要求局部地和适应性地加密网格。相应地,对二维水流计算问题采用了有限体积法的数学表示。跨单元边界的法向数值通量通过用Osher格式求解黎曼问题得到。文中给出了适合于二维浅水方程组的有关单元界面和各种物理边界的法向数值通量公式。分析了这一格式所具有的优点,包括:普适性、守恒性、逆风性、单调保持性、高效性、对间断的高分辨率、边界处理和内部格式相容,以及不引入数值边界条件等。最后,通过长江口南支水流计算的实例,阐明其良好性能。

浅水流动计算中—阶有限体积法Osher格式的实现

近年来，一阶有限体积法Ｏｓｈｅｒ格式已在二维浅水明流的一批模型问题和应用实例中获得成功。本文首次讨论其算法实现的种种问题。核心是建立单元水力模型—阶梯流，在数学上可用一类特殊的黎曼问题来描述。将该问题化作气体动力学中的黎曼问题近似求解，然后对结果加以校正。还在理论分析和数值试验的基础上详细讨论了各种外部边界条件、内部边界和动边界的处理，构成完整的算法。

浅水流动有限体积法/Osher格式的二维水流-水质模拟

应用无结构网格有限体积法 /Osher格式的二维水流 -水质耦合数值模型 ,采用黎曼近似解计算模拟区域水量、动量及污染物输送通量。通过对长江南京八卦洲江段浓度场的模拟计算 ,其结果与遥感卫星图片确定的污染带分布相吻合 。

COMPUTATION OF TURBULENCE FLOW ON HYBRID GRIDS USING k-ω TURBULENCE MODEL AND OSHER SCHEME

An unstructured Reynolds-averaged NavierStokes flow solver using the finite volume method is studied. The spatial discretisation is based on the Osher approximate Riemann solvers. A two-equatlon turbulence model (k-ω model) is also developed for hybrid grids to compute the turbulence flow. The turbulence flow past NACA0012 airfoil and the double ellipsolids are computed, and the numerical results show that the above methods are very efficient.

一类带小参数微分方程组的差分格式的一致收敛性的证明

近几年来，奇异摄动的研究已取得了很大的进展，但大多数工作都是针对单个方程的，对方程组的研究还不多，笔者对这类微分方程组作了有益的尝试，证明了差分格式的一致收敛性，其收敛阶为O（h）.

河口准三维涌潮数学模型研究

本文依据近年来计算流体力学在计算格式方面的最新成果,提出了一种用于计算包括涌潮在内的近海三维流场的数学模型。模型中,控制方程采用σ坐标变换形式;底部摩阻采用对数流速率;运用k ε紊动模型计算垂向紊动粘滞系数,使模型较好地给出水平速度的垂向分布。数值方法上,采用有限体积法离散方程,并用Osher格式计算有限体积交界面通量,使得模型可高精度地捕捉间断;水平方向用无结构三角网格,较好地适应不规则形状边界;用该模型模拟了明渠均匀流和概化河口涌潮的水流运动规律。计算结果表明模型可用于河口区涌潮流场计算。

悬浮物输运的数学模型

为了研究镇江内湖悬浮物的输移运动规律,建立了二维悬浮物输运数学模型。模型考虑了风场及潮流作用下悬浮物的输运,同时还考虑了悬浮物的沉降和再悬浮机制,采用切应力概念确定悬浮物的起悬量和沉降量。在无结构网格上对水流和悬浮物输运方程进行有限体积的积分离散,应用高性能的Osher格式计算水量、动量及悬浮物输运等通量,提高了数值模拟的精度。应用2次大规模水文和悬浮物含量的野外同步监测资料对模型进行了验证。模拟值与实测值吻合较好,表明模型能合理地反映内湖水流及悬浮物的运动规律。

黄河下游花园口-夹河滩河段二维洪水模拟

黄河下游花园口至夹河滩河段系典型的游荡型河段。在该河段 ,黄河大堤内范围宽广 ,一般洪水频率年份 ,水流主要限制在主槽内 ,因此大堤内分布有不少居民点以及纵横交错的保护居民点的生产堤和不少高于地面的灌溉渠堤和公路 ,使洪水行洪范围受到了很大的限制。当洪峰流量很大时 ,洪水将造成生产堤溃决 ,极大地危害滩区居民的生活。因此 ,设计模拟模型计算网格时需要考虑大堤、生产堤、明显高于地面的道路等阻水建筑物的影响 ,使这些堤及公路成为计算格网的边。不规则四边形网格能够很好地拟合黄河这种复杂的计算域。数值模拟时采用有限体积法 ,为确保通量的单向性 ,文中使用Osher格式计算通量。通过对 1 982年洪水的模拟 。

基于GIS的黄河下游二维水沙数学模型水沙构件设计

基于无结构网格的FVM方法离散微分方程,采用Osher及LSS格式计算对流通量,通过高斯积分求扩散项中控制体形心处的梯度矢量,预测、校正二步格式求解方程组,保证了水流计算格式的解具有良好的物理基础,符合解的特征构造,具有守恒性、逆风性、单调保持性。利用实测资料对非均匀沙沉速、分组沙挟沙力、紊动黏滞性系数计算公式等进行了适应性分析,综合确定了适应于高含沙水流的计算模式和方法。

基于混合网格上的高温非平衡流的数值模拟

新一代再入飞行器及空间传输系统需要了解其飞行时的气动问题及热力学问题。数值模拟高空、低密度、高熵及非平衡流动是一具有挑战性的问题。本文针对带座舱飞船高超声速再入大气过程中存在的严重气动加热现象,利用混合网格及Osher逆风格式,数值求解了三维化学非平衡Navier-Stokes方程,其中化学模型是5组分1 7个化学反应的空气化学模型,对带座舱飞船再入高度为40 km和马赫数为2 0 ,1 0的化学非平衡流场进行了数值模拟,给出了迎角为0°和2 0°情况下的各个组分的密度分布、压力等参数,并与量完全气体的计算结果进行了比较。

回流区水流运动二维数值模拟

建立了基于无结构网格有限体积法的平面二维水流数学模型。数值方法上,采用有限体积法离散方程,并采用基于黎曼解的Osher格式计算有限体积界面通量:采用三角形单元对计算区域进行剖分;提出了一种确定紊动黏滞系数的方法。用此模型对平面上突扩引起的回流区竖轴环流进行了模拟,流速断面分布的计算结果与实测结果接近,再现了因主、副流剪切层不稳定引起的涡脱落、涡向下游输移以及被竖轴环流吞并的过程,预报了回流区的水面共振现象。

典型城市型湖泊风生流的数值模拟

湖流是湖泊中污染物扩散、悬浮物沉降及输移等的主要水动力学要素。造成湖水运动的成因机制很复杂,其中风生流是重要形式。通过应用无结构网格有限体积法sher格式的二维水流数值模型,对不同风情下的典型城市型湖泊——南京市玄武湖湖流流态进行模拟研究,探讨其吞吐流较小,沿岸无入流的情况下的风生流特征。

冲积河流枢纽上下游冲淤问题的二维数值模拟

采用有限体积法/Osher格式,通过对冲积河流平原河段枢纽上下游水沙二相流的特点和边界条件的分析,建立二维水沙数值模型。对塔里木河乌斯满水利枢纽上下游进行模拟运算,其结果与物理模型实测数据相吻合,证明了有限体积法Osher格式模型的模拟的可行性和合理性。

滨海电站温排水数值模拟

针对滨海火(核)电站温排水海洋影响预测及评价的迫切需要,建立了二阶Osher格式水流———温度模型,采用“干湿单元水力模型”处理滨海电站所处海域复杂的计算边界,并结合某大型火电站温排水的数值模拟,显示了该模型正确模拟滨海电站温排水扩散过程的能力,为海洋环境影响评价提供可靠的技术资料。

高速气流吹袭防护气动特性的数值模拟(英文)

应用CFD方法研究了配备高速气流吹袭防护装置的人椅系统的气动特性。空间离散采用Osher逆风格式,紊流模型采用基于Spalart-Allmaras一方程的DES方法,在马赫数为0.6和1.2;迎角-10～30°;侧滑角为0°和15°范围内的粘性大分离流;以及来流马赫数为0.8和1.6;迎角5～30°;侧滑角为0,-20和-50°条件下,对有无防护装置的人椅系统气动特性进行了数值模拟。计算结果与试验结果吻合较好,高速气流吹袭防护装置在高速弹射过程中有效地保护了人体。