长江土脑子河段曲线二维动床模型研究

采用控制体积法与SIMPLE算法对一般曲线坐标下控制方程进行离散和求解,建立了长江土脑子河段曲线二维动床模型.模型控制方程离散时物理量布置采用非交错网格,避免了采用交错网格时所需的大量插值计算,提高了计算效率.模型较全面地考虑了非均匀悬移质及推移质运动,具有较好的普遍适用性.通过长江土脑子河段三峡水库蓄水以后2004年3~6月实测水沙资料对模型进行了验证,计算结果表明模型能够合理地反映土脑子河段水沙运动及河床变形规律.     

原始变量法计算封闭腔内自然对流

在文献[1，2]的基础上，利用原始变量法求解封闭腔内自然对流换热问题．对求解域采用非等距交错网格剖分，利用泰勒级数于网格点展开取二阶精度进行(u,v,θ)方程各项的离散，采用SIMPLE方法对压力及压力修正进行求解．不同瑞利数(Ra)条件下的数值试验显示，该数值方法物理概念清晰，计算稳定,可通过调节网格的疏密等方法，达到改善求解的收敛及提高求解的精度等目的．     

同位网格下平面二维水流数学模型的有限分析解法

为避免在求解平面二维水流数学模型方程时,流速、水位分别采用有限分析法和有限差分法求解而导致运动方程和连续性方程的计算精度不匹配,进而降低模型整体的计算精度和收敛速度.参照SIMPLE算法,推导出了能够在同位网格下采用有限分析法求解的一般曲线坐标系中的水位控制方程,实现在同位网格下采用有限分析方法求解全部控制方程的同精度算法,避免了交错网格,提高了计算精度,加快了收敛速度.运用于实际工程,计算结果与实测值基本吻合.     

五级旋风预热器冷模单体流场数值研究

采用雷诺应力(RSM)湍流模型对五级旋风预热器冷模单体流场进行了数值模拟，并同冷模测试结果进行了比较，取得了较好的吻合性。涡壳采用完全非结构网格，锥筒采用结构网格，有限体积法化湍流运动方程为差分方程．差分格式采用一阶迎风，压力速度耦合采用经典的SIMPLE算法数值求解湍流运动方程组，并在求解过程中采用网格自适应技术。通过对模拟结果分析，得到了旋风筒流场的基本特征。     

人体心血管系统中两相流动的数值模拟

本文用液体——颗粒两相流动的方程来描述心血管系统中的血液流动、对导出的二维两相流动基本方程用 SIMPLE 及 IPSA 算法进行了离散求解,算例的结果显示:两相系统的描述方法确实能对红细胞的一些运动特性给出恰当的解释,具有单相流动法不可比拟的优点。     

非对称突扩通道流动特性的数值模拟

针对非对称突扩通道的流动,建立了平面二维不可压稳态层流数学模型,采用守恒性能高的有限体积法对控制方程进行离散,并选用绝对稳定的混合格式,利用交错网格下的SIMPLE算法进行流场数值计算,得到了回流区形态的变化规律.研究表明,非对称突扩通道流动受雷诺数和几何形状的影响.     

粘性不可压缩三维流场分块算法研究

提出了重叠对接耦合分块方法即为隐式求解的块迭代方法．采用自主开发的三维粘性不可压缩流场计算程序对重叠对接与不重叠对接的流场分块求解方法进行了比较研究．并与他人的实验结果进行了比较．研究表明,两种算法均能获得收敛解,而利用重叠对接分块法的收敛性优于不重叠对接分块法．Schwarz迭代法用于重叠分块网格的区域解法,而Dirichlet-Neumann迭代法用于不重叠分块网格的区域分解法．像Abdallah方法一样,控制方程在非交错网格下离散．如果网格雷诺数大于2则在对流项中采用逆风格式．     

用改进的强隐式方法求解三维不可压粘性流动问题

提出了一种用改进的强隐式方法求解三维不可压缩粘性流动的交错网格有限体积法。该方法在ＳＩＭＰＥＲ算法的基础上，对Ｎ－Ｓ方程的对流项采用ＱＵＩＣＫ格式离散，利用ＭＳＩ方法求解代数方程，数值计算表明本方法具有计算稳定，收敛速度快等特点。     

基于LES模型的水轮机三维非稳定湍流及压力脉动数值模拟

本文基于N—S方程和大涡模拟（LES）模型，采用贴体坐标和四面体网格系统和SIMPLE算法，较准确地对混流式水轮机内部三维非稳定湍流流动及其尾水管内的旋涡压力脉动进行了数值模拟。     

基于两相流理论的低浓度挟沙水流运动数值模拟Ⅱ—数学模型的求解与验证

基于非结构网格上控制方程的离散思想,提出了曲线网格上低浓度挟沙水流控制方程离散的新方法.该方法采用有限体积法对控制方程进行离散,用基于同位网格的SIMPLE算法处理水流运动方程中水深和速度的耦合关系,且不需要对控制方程进行坐标转换,因而离散方程的形式简单,物理概念清晰,所得到的离散格式既适用于正交曲线网格还适用于非正交曲线网格,能够适应天然河道这类几何形状复杂的物理边界,从而使编程和计算变得更加简单有效.采用海河流域漳河的实测资料进行了验证计算,计算所得的水位、流速与河床冲淤变形与实测资料基本一致.     

深圳大运会主体育场罩棚风压分布特性的数值模拟研究

基于Reynolds时均 方程和RSM模型对深圳大运动会主体育场罩棚风压分布特性进行了数值模拟,文中采用有限体积法对控制微分方程进行离散,并采用SIMPLE压力校正算法来实现非线性方程的迭代求解。分析了风向角和风攻角对风压分布的影响,探讨了罩棚周围流场的绕流特性,并根据罩棚的结构形式及风压分布特点,给出了最不利风向下罩棚表面的区域平均风压系数以及罩棚设计风荷载的计算方法,以供工程设计的参考。     

贴体曲线坐标系中的非交错网格的SIMPLE算法

提出了在贴体坐标系中非交错网格的ＳＩＭＰＬＥ算法,通过对控制界面上的逆变速度分量作特殊加权平均处理,可避免在非交错网格中产生压力波动而连续性方程又无法检测出来的问题。     

薄壁梁桥的强非线性弯扭耦合振动解析解

采用非线性Rayleigh阻尼描述在运动初期有激发使运动加速至最大速度后开始衰减直到运动停止。本文的控制方程组是两个耦连的非线性偏微分方程,选用适合所有边界条件的两个振型函数代入控制方程组,应用Fourier级数的运算,得出两个无穷系列耦连的以时间为自变量的未知函数的二阶非线性常微分方程组,用迭代加法每次只取一个未知函数而化为一个van der Pol方程,得出其头两次的渐近解,讨论了解的跳跃特性,最后代入具体数据,得出数值结果。本文是文献（7）的深入,因（7）是弱非线性。     

基于时域边界元法的散热结构瞬态热传导分析

针对散热器结构的瞬态热传导问题,首先,在热力学理论基础上,利用问题的控制方程推导出问题的积分方程;然后针对积分方程中的域积分,采用双互易边界元法(DRBEM)进行处理,得到边界积分方程;再对其进行边界离散,获得常系数微分方程组;最后,运用精细积分法(PIM)进行方程组求解,得到内部点的温度结果.通过边界元法与有限元法计...     

圆形沉淀池内温差异重流与浮力流的数值模拟

为了了解圆形沉淀池内水流在无温差及存在温度差时的特点,建立了描述圆形沉淀池内水流运动的k-ε双层湍流数学模型,采用交错网格系统,压力与速度解耦时采用了SIMPLE方法,控制方程的求解方法为混合有限分析解法。对中心进水周边出水的圆形沉淀池的模拟结果表明所建立数学模型的模拟结果与文献中的资料吻合良好。在无温差时且低污染物浓度<150~200 mg/L时,池内存在一个很大的回流区;低温水进入高温水池时,常会产生异重流现象;高温水进入低温水池时,会产生浮力流现象。     

Static characteristics of new type externally pressurized spherical air bearings

In order to provide some theoretical guideline for the structure design of the new type externally pressurized spherical air bearings,the static characteristics and the factors affecting the static characteristics of the air bearings were analyzed.A finite volume method was adopted to discretize the three-dimensional steady-state compressible Navier-Stokes equations,and a modified SIMPLE algorithm for compressible fluid was applied to solve the discretized governing equations.The pressure field and velocity...     

深圳大运会主体育场罩棚风压分布的数值模拟

为了求得体育场罩棚的风载体型系数,基于Reynolds时均N-S方程和RSM模型对深圳大运动会主体育场罩棚风压分布特性进行了数值模拟.采用有限体积法对控制微分方程进行离散,并采用SIMPLE压力校正算法来实现非线性方程的迭代求解.将数值计算结果与风洞试验结果进行了对比,两者吻合较好,在此基础上,分析了风向角和风攻角对风压分布的影响,探讨了罩棚周围流场的绕流特性,并根据罩棚的结构形式及风压分布特点,给出了最不利风向下罩棚表面的区域平均风压系数以及罩棚设计风荷载的计算方法.结果表明:风向角对风荷载的影响较大,不同风向角下,来流的分离和漩涡脱落作用均有较大的不同.风攻角对罩棚表面风压分布的影响不大,对罩棚局部峰值的影响在15%以内,整体影响在5%以内.