任意形状长轴对称体的Oseen流动

研究了任意形状长轴对称体Oseen流动的尾流效应,选取Oseen流子作为基本奇点沿用文[6]中的方法得到了流场的封闭形式的级数解,对不同细长的长比球进行了数值计算并求出了长球的阻力系数、研究了在不同细长比和雷诺数下尾流效应对长球表面上压力分布的影响。     

弱可压缩流体流动模型及数值求解

 介绍了小马赫数下弱可压缩流体流动的模型,对该模型的数值求解方法也作了较为详尽的叙 述。讨论了紊流模型模拟涡流的方法。简要介绍了水力机械内部流场的计算实例。     

弱可压缩流体流动模型及数值求解

介绍了小马赫数下弱可压缩流体流动的模型 ,对该模型的数值求解方法也作了较为详尽的叙述。讨论了紊流模型模拟涡流的方法。简要介绍了水力机械内部流场的计算实例。     

微通道流体物性对流动与传热影响的数值模拟和实验研究

针对微矩形通道内单相层流流动和传热特性进行实验研究和三维数值模拟研究.微矩形通道宽为50μ m,高为200μ m,工质采用去离子水.讨论流体热物性随温度变化对流动传热特性的影响.并把实验数据、模拟结果和文献参考值比较.结果显示,恒热流边界条件下,由于热量加入导致流体温度沿流动方向升高,摩擦系数沿流动方向有递减趋鼢势.进口处的Nu数有最大值,在L*=0.2处降到最小值,之后又攀升到充分状态时的Nu数值.在充分发展段,Nu数保持不变.黏性耗散对热起始段的作用不能忽略.     

用边界元法计算任意形状混合边界的油藏压力

建立了考虑源（汇）影响、任意形状、混合边界条件、含有不渗透区油藏稳定渗流的数学模型,采用边界元方法对其进行了求解,获得了油藏内的压力分布,绘制了相应的压力剖面图。通过对比压力剖面图和平均压力发现,由于油藏的封闭边界部分不能向生产井提供能量,致使地层能量被消耗较多,油藏平均压力降低;由于油藏内不渗透区的存在,能量在绕过其向生产井传播时要有部分能量损失,所以油藏的平均压力会降低。     

任意拉格朗日—欧拉（ALE）法求解三维湍流流动

本文推导了旋转叶轮内流湍流Ｎ－Ｓ方程计算的ＡＬＥ格式。以网格的虚拟运动计算对流作用；采用粘性子循环时步、拉格朗日子循环时步加速收敛；改进上风格式为部分施主单元格式来减小扩散同时实现良好的计算性能；由通量限制控制耗散截断误差；计算中引入短波长阻尼以消除不规则性；以亚网格（ＳＧＳ）湍流模式模拟湍流运动。     

基于任意形状网格和精确几何边界的数值计算

有限元网格形状要尽可能规则,网格之间必须通过结点连接,这些要求给复杂形状求解域的数值计算带来很大的前处理工作负担,而且实际的曲线边界一般要离散成有限单元能够描述的形式,难以模拟CAD模型的精确几何。针对这些问题,基于独立覆盖流形法提出任意形状且任意连接的覆盖网格,在CAE分析中模拟CAD模型的精确几何边界及其边界条件:将求解域划分为可包含曲线边的任意形状的块体网格,可以采用单纯形解析积分和数值积分2种方式进行块体积分;仅需在积分过程中考虑块体之间的窄条形(包括曲线条)的覆盖重叠区域,而不必在计算模型中生成这些条形;通过边界条实现本质边界条件的严格施加,包括曲线上的边界条件;给出2个数值算例验证了方法的有效性。任意形状的覆盖网格将为实现基于精确几何模型的数值计算及其完全自动化的前处理开辟新的路径。     

任意形状港湾波浪组合绕射和折射的数值分析

本文基于Berk hoff的缓坡方程,用杂交元方法对波浪组合绕射和折射进行了数值分析。针对大区域中带有复杂边界条件的建筑物编制了通用程序,命名为“CRDWCE”。本文对杂交元方法(HybridElement Method,HEM)进行了改进,改进的HEM法可通过增加级数解项数来提高解的精度,而不增加联立方程组的数目。本文采用同类问题中较简单的泛函公式,简化了计算。本文采用不完全Choleski分解法(Incomplete Choleski-Conjugate Gradient,ICCG)解大型复代数方程组以节省CPU,用索引存贮法以节省内存,使得小机算大题成为可能。通用程序发展了功能很强的自动剖分程序。本文结果与理论解吻合很好,其算法可供实际工程规划和设计中应用。     

含中心孔的椭圆区域内位势流动分析

本文利用相似流动替换方法，解决了中心有圆孔的椭圆形区域内的位势流动问题；采用分区域解法，给出了中心有椭圆孔的椭圆形区域内位势流动问题的解析通解。该问题在油气田开发工程等许多领域中有重要意义。     

任意形状混凝土骨料的数值模拟及其应用

提出了平面二维、空间三维任意形状混凝土骨料的模拟方法。通过一次性形成所有同级骨料的基骨料，再对所有的基骨料进行随机延拓，生成任意形状的随机骨料。提出的算法避免了逐个投放效率低、生成的骨料含量低的缺点。三维骨料任意形状的模拟中，假定同一级配中的骨料粒径呈正态分布，建立了以体积为标度的凸型骨料的生长模式和骨料凸性条件，根据多面体的每个面的面积是否大于限定面积，确定该面是否需要延拓，同时以延生点与延伸面以外的所有边界面组成的体积都应该大于零作为凸性限定条件。采用体积和测试法进行骨料的侵入判别。在此基础上研制了相应的二维骨料、三维骨料随机投放软件2D－RAS和3D－RAS。生成的数值混凝土二维骨料含量可达65％以上，三维骨料含量可达55％左右，基本可以满足全级配混凝土的生成需要。最后还对三维混凝土任意骨料形状的弯拉梁的破坏过程进行了初步研究，得出任意多面体骨料试件的极限承载力要大于圆球骨料试件，任意多面体骨料试件所对应的位移相对较大等结论。     

任意形状二维物体上的二阶波浪力

用局部变分法计算了二维物体上的二阶波浪力。结果表明,在若干个频率处出现了“不规则频率”,这是一个新现象。但与实验结果的比较表明:用简单的内插通过这些不规则频率,不影响结果的精度。     

任意形状复合油藏压力动态的边界元分析

油气藏形状对井底压力动态响应有着显著的影响,常规的解析方法只能计算规则形状油气藏的井底压力响应。该文建立了任意形状复合油藏不稳定渗流的数学模型,并将区域分割法与边界元方法相结合,研究了任意形状复合油藏不稳定压力响应的计算方法,应用杜哈美原理考虑了表皮效应与井储系数对井底压力动态响应的影响。绘制并分析了复合油藏的典型曲线,通过分析指出油藏的边界性质会影响典型曲线后期形态,而内区形状的也可能掩盖复合油藏典型曲线的特点。     

高粘幂律非牛顿流体在管内流动和传热的计算

本文的数值研究考虑了变物性,自然对流和粘性耗散的影响,采用SIMPLE方法进行数值求解,计算结果与实验值吻合很好,流场和温度场耦合,使进口段大大延长,计算结果还显示了高粘物系加热时径向温差很大。     

边坡稳定分析的任意形状滑动面的简化Bishop法

放弃了传统简化Bishop法的整体力矩平衡条件的假设。根据土体沿滑动面极限平衡的充分必要条件定义滑动稳定安全系数,使简化Bishop法适用于任意形状的潜在滑动面。给出了适用于任意形状滑动面的简化Bishop法的安全系数计算方法以及相应的最危险滑动面的搜索方法。计算结果表明本文方法具有较高精度,且滑动面位置符合实际,该方法可以应用于实际工程的边坡稳定分析。改进的方法概念简明,计算简单,适用性强。     

求解浅水流动的分步有限元方法

本文将分步有限元的计算方法引入到浅水方程组的求解中。该方法起源于Taylor-Ga!erkin(T-G)方法，但数值稳定性优于T-G法并具有三阶精度。由于计算中没有引入高阶的空间导数项，实现起来比Taylor-Galerkin方法简单，适用于非线性和多维问题的求解。计算模型中包含了零方程和双方程的紊流模型，可以根据需要选择。文中详细介绍了初始和边界条件的取法，并通过五个算例验证了计算模型的可靠性。     

水准网中任意点权倒数的求解方法

在水利水电测绘外业生产中，水准网的布设常用于高程控制。但水准网如何布设，网中各点精度分布情况如何，是在外业实施以前就应该知道的，因此必须具备计算网中任意点权倒数的能力。本文以双结点水准网为例，给出网中任意点权倒数的数学模型。1等权代替法图1是双结点水准网图形，图中A、B、C、D是已知点，K1、K2是结点，S1、S2、S3、S4是已知点至结点的路线长度，S5是结点之间的长度，S′为结点K2至任意点G之间的长度，①、②、③、④、⑤为线号，Ⅰ、Ⅱ、Ⅲ为环号。若以每公里观测高差权重为单位权，则：由A至几算得K1的高程，其权…     

波纹板内流体流动特性的数值模拟

以室温下的水为介质，以涡量一流函数为理论基础，运用代数变换法生成适体坐标，通过有限差分对波纹板板内流场进行了数值模拟；并以Re=500为例，求出其板内流体的速度场和压力场，进而对同一周期内和不同周期的速度场和压力场进行分析，给出了速度和压力分布变化曲线图，得出了不同雷诺数下流动阻力损失变化曲线与阻力关联式，为波纹板的强化传热研究提供了有利的参考。     

双极流体在狭窄管道内的流动

本文研究了双极流体通过带有轴对称柔和狭窄的刚圆管时的流动特性。运用Bleustein和Green的双极流体理论,导出了速度剖面、流量、阻抗和壁面最大剪应力的公式,并做了数值计算,分析了狭窄度与双极参数对上述有关量的影响。