浮动步长法在静电场数值计算的应用

静电场中的三维区域问题是比较复杂的，而蒙特卡罗（MonleCorlo）方法则是解决这类问题的一种方便且灵活的手段，浮动步长法即是一种改进的蒙特卡罗方法，文中用此方法解决拉普拉斯方程的数值计算问题，用格林函数提高计算效率，电位值和电场强度值均可直接计算，减少了误差。     

Hubbard模型中格林函数的数值计算

从Hubbard模型出发,由动态平均场论方法结合量子蒙特卡洛模拟数值求解系统的格林函数,重点讨论在顺磁序和反磁序情形下格林函数的数值结果,给出在不同温度和不同Hubbard能下格林函数的变化情形,从格林函数的变化可以分析得到在温度降低或Hubbard能升高时,系统会发生从金属到绝缘体的Mott转变.     

复变量Spence函数的数值计算

介绍了一个可计算复变量Srence 函数的Fortran 程序.该程序克服了西欧核子中心(CERN)程序库中相应程序只能计算实变量情形且仅给出其函数值实部的缺点.结果表明这个程序比CERN 库中的程序有很大改进.     

数值逼近法在酸碱计算中的应用

本文讨论了近似计算的另一种方法一数值逼近法。供从事分析化学的师生参考。     

一种新型的数值计算方法——基函数法

我们提出一种新型的数值计算方法——基函数法．此方法直接在非结构网格上离散微分算子．采用基函数展开逼近真实函数,构造出了导数的中心格式和迎风格式．取多项式为基函数并采用通量分裂法及中心格式和迎风格式相结合的技术以消除激波附近的非物理波动,我们构造出数值求解无黏可压缩流动一阶多项式的基函数格式．通过一、二、三维多个无黏可压缩流动典型算例的数值计算表明本方法是一种高精度的,对激波具有高分辨率的无波动新型数值计算方法,与网格自适应技术相结合可得到十分满意的结果．     

δ冲击模型寿命分布的积分计算及M函数的性质

冲击模型特别是δ冲击模型在可靠性理论中具有重要的作用.讨论了δ冲击模型寿命分布中多重积分的计算问题,并研究了由此积分引出的M函数的性质.     

磁轫致辐射辅助函数的特性和数值计算

磁轫致辐射函数在小宗量和大宗量区域中的展开特性。提出了一种渐近和多项式展开与数值积分相结合计算磁轫致辐射辅助函数的方法。以该算法编制的磁轫致辐射辅助函数子程序可以有效地应用于天体磁轫致辐射谱分布特性的数值计算。     

模糊数值正弦函数与余弦函数

利用扩展原理引入了模糊数值正弦函数与余弦函数，并研究了这两种模糊函数的基本性质．     

具有偶数自由度的F分布函数的数值计算

给出了关于具有偶数自由度的F分布函数在偶数自由度方向进行递推的数值计算.它与常用的交替递推算法相比,减少了计算量和累积误差.针对该方法可能出现的超字长问题,亦提出了分析和解决办法.     

激光相变硬化数值计算

根据实验的实际条件，考虑对流，辐射等复杂边界条件，给出了有限尺寸，移动热热源的激光相变硬化三维温度场数学模型；     

液力透平的数值计算与试验

设计了液力透平试验台,对一单级液力透平进行了试验,得到了外特性曲线.采用全流场和结构化网格技术对液力透平内部流动进行了数值计算.分析了液力透平在不同流量下的压力场和速度场,得到了内部流场的分布规律.应用速度三角形对液力透平叶轮和尾水管内部速度场随流量变化规律进行了研究.结果表明:离心泵反转可用作透平运行,并具有较高的效率;最高效率的数值计算与试验结果相对误差为4.85％;透平内部的压力从蜗壳进口经叶轮到尾水管逐渐减小,进出口压差随流量增加而逐渐增加;在透平叶片背面和工作面存在漩涡区域,漩涡位置和区域大小随流量而变化;在尾水管横截面上存在的圆周速度分量随流量而变化.     

有界变差函数的Müntz有理逼近

给定M＞0,设∧={λn}∞n=1是满足0≤λ1＜λ2＜…的实数序列,且对所有n≥1,有λn 1-λn≥Mn,文中得到了由Müntz系统{xλn}构成的有理函数对有界变差函数在Lp范数下逼近的一种估计.     

非完整桩扭转振动问题的数值计算

文章提出了非完整桩在瞬态激振力矩作用下扭转振动的数值计算模型,并用有限差分法计算了非完整桩的扭转振动,得到了各种典型缺陷桩对应的理论曲线。定量计算结果表明,桩周土的性质是影响扭转波衰减的一个重要因素;桩周土越软,扭转波衰减越慢,缺陷反射和桩底反射越明显,此时利用扭转波进行基桩完整性检测更有效。     

凹壁不稳定性流动的数值计算

凹壁结构在工业应用的多个领域中都是很常见的,由于离心力不稳定性,流体流过凹壁面附近时可能会产生G?rtler涡。为了研究这些涡的特性及对流场产生的影响,采用CFD技术对曲率半径为2 m的弯曲渠道的流动过程进行数值模拟。用有限体积法对模型进行空间离散,采用大涡模拟中Smagorinsky-Lilly模型进行计算。通过与实验结果进行对比发现,计算结果能较准确的反应G?rtler涡随流动的发展过程。在此基础上进行的分析表明:在速度3 m/s,湍流度0.35%的进口条件下,沿流向流过一段距离后,在靠近凹壁的边界层中产生了G?rtler涡;涡轴平行于流动方向并且相邻涡轴之间的距离保持不变;随着涡在展向和法向空间的增长,速度等值线出现"蘑菇状"结构;靠近出口位置"蘑菇状"结构破裂,流动发生转捩。     

高超音速飞行器头罩气动热流场数值模拟

对几种不同外形的高超音速飞行器头部气动热效应进行数值模拟,采用Fluent软件对各飞行器进行三维薄层N-S方程数值计算,得到了不同外形的高超音速飞行器头罩外表面热流场分布.然后将计算结果与美国宇航局的超音速飞行器X33模型风洞实验结果进行分析和比较.分析表明:在五倍音速近地飞行时,其中两种外形头部的高超音速飞行器安装中波红外光学成像探测器制导是比较适合的.     

旋流燃烧室中NO排放的数值计算

采用数值模拟的方法研究了旋流燃烧室中旋流数对NO排放的影响，其中湍流采用修正的k-ε模型，燃烧采用EBU-Arrhenius模型，NO的生成包括热力NO和瞬发NO，采用设定的PDF模型．计算结果表明，随着入口旋流数的增加，旋流燃烧室中的高温区移向入口，火焰长度缩短，宽度增大；NO排放量先增加，后减小．因此，适当调整燃烧室的旋流数可以降低NO排放。     

船舶浅水航行下沉量数值计算

为求解船舶航行中的下沉量,采用NURBS技术对船体进行曲面重构,并应用基于Rankine源格林函数的一阶面元法构造数学模型,然后结合两者得出船舶下沉量的计算方法.以标准船型KCS为例进行验证,计算值与实验值的最大和最小相对误差分别为-10.32％和-1.84％,与实验情况基本相符,说明该算法具有可行性.     

2.5MW煤粉氧燃烧旋流火焰数值计算

为了准确地数值预报煤粉旋流火焰特征和主要烟气成分,比较扩展的涡耗散模型在多湍流模型下的预测性和影响,针对IFRF(国际火焰研究基地)2.5MW煤粉氧燃烧和空气燃烧实验,利用组分输运方程,结合7步改进的总包反应,对炉内煤粉氧燃烧进行数值计算,并与实验结果和空气燃烧基况对比.结果表明:氧扩散率、发射率和比热容等物性参数修正后,湿循环、氧预混二次流的煤粉氧燃烧扩散火焰温度场与空气燃烧总体一致,火焰稳定且为Ⅱ型火焰结构;EDC模型(涡耗散概念模型)对组分预测更为准确,尤在可实现k-ε模型下的缓慢反应CO生成和IRZ(中心内回流区)高温预测;而FRED模型(有限速率涡耗散模型)在多湍流模型下对炉内温度、组分、火焰结构预测亦较准确.     

三体船附加质量的数值求解及分析

以船舶操纵水动力预报为研究背景,以无升力势流理论为基础,采用源分布面元法对半椭球体和Wigley船在水平面内运动的附加质量进行计算,将由附加质量所得孟克力矩与无升力势流理论物面压力分布所得孟克力矩理论计算结果进行对比,验证本方法计算结果的正确性及其用于孟克力矩计算的适用性.根据本计算方法求解不同侧体布局下三体船的附加质量,并分析侧体布局对三体船附加质量的影响,所得结果和结论符合实际船舶操纵运动规律.