非正交曲线坐标系平面二维电厂温排水模拟

结合河道水流及温排水的运动特性，从非正交曲线坐标系下温排水基本方程出发，采用有限体积法及SIMPLE算法离散求解方程，建立了基于非正交曲线坐标系下非交错网格的平面二维温排水数学模型．用该模型模拟天然河流电厂温排水运动的计算结果表明，该数学模型能合理反映河段电厂温排水运动。     

贴体网格中几种调节因子比较

基于有限差分求解椭圆型方程的方法，对假定的扩展弯曲河段和明渠扩展段，采用几种常见的调节因子分别生成了正交贴体坐标系，并对计算结果进行了比较。     

伶仃洋三维流场数值模拟

使用σ坐标变换技术，将（ｘ，ｙ，ｚ，ｔ）坐标系中的三维非线性Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程转换成（ｘ，ｙ，ａ，ｔ）坐标系中的方程，采用有限差分方法进行离散求解，用垂向（σ向）积分的连续方程求解水位，用动量方程求解流速垂直结构，在水平方向上使用非交错正方网格，σ向等间距分层，从而建立三维流场数学模型。应用本模型对多口门、多岛屿、地形及岸线复杂的伶洋的枯季、洪季流场进行了数值模拟。将模拟结果与多个水位、流速同步测站的实测资料进行比较，无论是水位、垂线平均流速、流向过程线还是流速的垂向分布形式与量值，二者均吻合较好。根据模拟计算结果，对伶洋流场的平面及空间结构特征进行了分析。     

平面二维河床冲淤计算有限体积法

给出了边界层坐标系下平面二维河床冲淤计算数学模型方程 ,并用河势贴体正交曲线二维网格形式来划分计算区域 ,模型控制方程及定解条件采用有限体积法进行离散 ,数值求解基于SIMPLE算法。该模型用武汉河段天然实测资料进行了验证 ,并在工程实践中得到了大量的应用和检验     

应用RNG k-ε湍流模式数值模拟90°弯曲槽道内的湍流流动

本文将Yakhot与Orszag提出的RNGk-ε湍流模式与壁面律相结合,应用于90°弯曲槽道湍流流动的数值模拟.计算在任意非正交曲线坐标系下进行.并采用非交错网格有限体积方法求解三维不可压N-S方程.文中给出了详细的计算结果,并与实验结果进行了比较,结果表明,RNG k-ε湍流模式能够有效地模拟有曲率影响的湍流流动.     

非交错曲线网络的水流数值计算模型

在非正交边界拟合曲线坐标系下，建立非交锚网格变量布置的水流计算 模式。采用逆变速度分量作为独立计算变量。为避免出现棋盘格式压力场物理上的非真实解，采用动量插 值方法计算交界面处的流速插值，运用压力 校正法（SLMPLE、SIMPLEC或 SIMPLER）求解离散方程，并对S型弯道水流运动 进行数值计算。     

长江口边界拟合坐标的三维潮流数学模型

在正交曲线坐标系下，采用坐标变换，有限差分、交替隐式求解、动边界处理等方法，建立了边界拟合网格的三维潮流数学模型。应用在长江口潮流场的计算中，得到了与长江口实测潮流，流速和流向相一致的结果，能较好地反映长江口复杂边界和地形下的三维潮流变化。     

通量差发裂上风差分格式求解船舶粘性流动

本文阐述了采用人工可压缩方法计算不主缩粘性流体的流动问题，对一般曲线坐标系下的流动控制方程采用了通量差分分裂上风差分格式的二阶中心差分格式离散。在隐式求解格式下，代数方程组用线松弛方法求解。     

弯曲河流三维数值模型

将平面正交曲线坐标系和垂向σ伸缩坐标系相结合,建立适合弯曲河流的三维水流数值模型。采用模式分裂技术将控制方程分解为外、内模式联合求解,方程的离散采用有限差分法。采用嘉陵江广元城区河段5个弯段内7个典型横断面地形、水位和流速资料对本文模型进行验证,模型成功地模拟出了水面横比降和弯道二次流,模拟结果与实测数据吻合较好。     

船舶螺旋桨周围粘流场数值预报与流场分析

本文叙述了通过直接求解雷诺平均应力方程（ＲＡＮＳ）来获得船舶螺旋桨粘流场数值解的方法与数值求解步骤，该方法采用非交错网格系统，利用幂函数格式离散动量方程，预报－校正方法求解速度－压力耦合问题，应用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式求解雷诺应力项使方程组封闭，以此来获得螺旋桨粘流场，为了验证数值方法的稳定性和可靠性，文中以ＤＴＲＣ４１１９桨为算例，对螺旋桨粘流场进行了数值求解，通过对计算结果的分     

HIGH ACCURATE SOLUTION OF INCOMPRESSIBLE VISCOUS FLOW WITH PRIMATIVE VARIABLE

This paper presents a higher order difference scheme for the computationof the incompressible viscous flows.The discretization of the two-dimensional incompress-ible viscous Navier-Stokes equations,in generalized curvilinear coordinates and tensor for-mulation,is based on a non-ataggered grid.The momentum equations are integrated intime using the four-stage explicit Runge-Kutta algorithm [1]and discretized in space us-ing the fourth-order accurate compact scheme[2]The pressure-Poisson equation is dis-cretized using the nine-point compact scheme.In order to satisfy the continuity constraintand ensure the smoothness of pressure field,an optimum procedure to derive a discretepressure equation is proposed [9][3]The method is applied to calculate the driven cavityflow on a stretched grid with the Reynolds numbers from 100 to 10000.The numerical re-sults are in very good agreement with the results obtained by Ghia et al [7]and includethe periodic solutions for high Reynolds numbers.     

在低中高三种Re数下圆管入口流的有限元分析

圆管入口流问题的研究已有很长历史,但理论结果与实验结果不能很好地符合。本文用改进的Newton-Raphson方法线化由Navier-Stokes方程导出的有限元方程,用波阵法求解线化后的方程组。最后在0～800的Re数范围内得到了速度场及入口长度,并给出了描述入口长度随Re数变化关系的拟合公式。数值解与实验结果的比较是令人满意的。本文还对简化Navier-Stokes方程(SNSE)与完整Navier-Stokes方程(CNSE)在求解本问题过程中的一致性作了数值研究,得出了当Re≥80时CNSE和SNSE结果一致的结论。     

三维方腔拖曳流的数值模拟

本文推导出一种适用于定常和不定常粘性不可压缩Navier-Stokes方程的分裂步方法。采用Taylor-Galerkin有限元格式进行求解，对有限元等式中关于速度的时间项进行三点向后差分，深入考虑粘性不可压缩流Navier-Stokes方程中对流项的作用，利用二阶Taylor展开完成时间项向空间项的转化，采用张量分析的方法推导了N-S方程分裂步方法的有限元离散格式，并采用低Reynolds数三维方腔拖曳粘性流^[23][24]作为基本算例，检验了这种分裂步方法的稳定性和有效性，同时与大涡模拟相结合对Reynolds数为10000的三维方腔拖曳湍流流场进行了相关的分析，进一步揭示了方腔回流运动的非定常非对称性、流动结构表现为竖轴环流与立面环流相叠加、流速沿垂线分布相对均匀等流动规律，显示了该方法与大涡模拟相结合能够有效地捕捉涡系及其时变过程。     

VOF方法模拟水面振荡流场

应用VOF(Volume of Fluid)方法结合不可压缩流体的N-S方程和连续性方程对矩形水池中的水面振荡流场进行数值模拟.采用共轭剩余法结合差分法求解N-S方程和连续性方程,采用施主-受主法求解流体体积函数控制方程.计算得到一个振荡周期内典型时刻矩形水槽水面振荡流场分布情况并与理论分析结果进行比较.并将计算得到的...     

WAVE-BODY INTERACTIONS FOR A SURFACE-PIERCING BODY IN WATER OF FINITE DEPTH

Nonlinear wave-body interactions for a stationary surface-piercing body in water of finite depth with flat and sloping bottoms are simulated in a two-dimensional numerical wave tank,which is constructed mainly based on the spatially averaged Navier-Stokes equations with the model for simulating the turbulence.The equations are discretized based on the finite volume method and the scheme of the pressure implicit splitting of operators is employed to solve the Navier-Stokes equations.By using the force time h...     

Large eddy simulation of water flow over series of dunes

Large eddy simulation was used to investigate the spatial development of open channel flow over a series of dunes. The three-dimensional filtered Navier-Stokes (N-S) equations were numerically solved with the fractional-step method in sigma coordinates. The subgrid-scale turbulent stress was modeled with a dynamic coherent eddy viscosity model proposed by the authors. The computed velocity profiles are in good agreement with the available experimental results. The mean velocity and the turbulent Reynolds st...     

A STUDY OF OSCILLATING FLOWS OVER PROPAGATING RIPPLES:PART I. STEADY STREAMING

1.　ＩＮＴＲＯＤＵＣＴＩＯＮＯｓｃｉｌｌａｔｉｎｇｆｌｏｗｓｏｖｅｒｐｅｒｉｏｄｉｃｒｉｐｐｌｅｓａｒｅｏｆｐｒａｃｔｉｃａｌａｓｗｅｌｌａｓｓｃｉｅｎｔｉｆｉｃｉｎｔｅｒｅｓｔｂｅｃａｕｓｅｏｆｔｈｅｉｒｒｅｌｅｖａｎｃｅｔｏｂｅａｃｈｐｒｏｃｅｓ...     

应用等位函数法模拟溃坝流动

 采用等位函数法配合三维不可压缩纳维叶-斯托克斯方程式的计算方法求解溃坝流场。控制方程式采用有限差分法求解,并利用投影法导出压力的波松方程,使得能自动满足连续方程式。此外,利用了WENO法求解等位函数方程的空间离散项。最后以三维溃坝问题做计算,发现数值计算结果与试验数据相当吻合,显示出数值模式的成功。     

Irreversibility analysis of unsteady couette flow with variable viscosity

This paper investigates numerically the inherent irreversibility in unsteady generalized Couette flow between two parallel plates with variable viscosity. The nonlinear governing equations are derived from the Navier-Stokes equations and solved numerically using a semi-discretization finite difference method together with the Runge-Kutta-Fehlberg integration scheme. The profiles of velocity and the temperature obtained are used to compute the entropy generation number, Bejan number, skin friction and Nusselt number. The effects of embedded parameters on entire flow structure are presented graphically and discussed quantitatively.     

A NEW FINITE DIFFERENCE SCHEME FOR 2-D FLOWS AT LARGE Re

Based on the conventional ADI method and SOR method for solving the Navier-Stokes equations, averaging those linkage terms which also are the terms of vary-ing coefficients of the equations, a new finite difference scheme-Averaging Finite Diffe-rence (AFD) scheme for 2-D flows was proposed. A 2-D driven cavity flow was calculated numerically as an example with the presented scheme at Re= 100, 1000, 2000, 3200, 10000. The results were discussed and compared to those obtained with the conventional methods as well as experimental data. It showed that a slight change of the approximation pattern of the conventional scheme in the terms of varying-coefficients of the governing equations seems to have an obvious influence on the solutions at high Re which will be erro- neous if the conventional schemes was employed.