数值解法的另一种欧拉改进法

在求常微分方程初值问题的数值解时,本文在欧拉方法计算结果上再补充一个拉格朗日余项的近似值,以期达到提高精确度的效果,并由此而构造了一个新的计算格式.     

城市大气中污染物扩散的施密特数研究

借助反问题方法对风洞实验工况选取不同施密特数进行数值模拟试验,分析比较了不同施密特数时街道峡谷内各个观测点的数值模拟数据和风洞实验数据.研究结果表明,城市街道内大气污染物扩散的施密特数为0.7～0.8.     

低渗气藏多次应力敏感测试及应用

低渗透气田在我国已开发气藏中占有相当大比例，应力敏感性对气田开发的影响已逐渐引起人们重视。针对目前大多数低渗透气田应力敏感测试仍然采用定内压变围压测试方式所存在的与实际气田开发过程中内压变化上覆压力恒定的变化规律不一致的问题，研究提出了变内压定围压测试方法，采用该方式测试低渗透岩芯在净上覆压力六升六降过程中应力敏感性,并应用测试的应力敏感曲线对低渗透气井试井特征及单井生产动态进行研究。结果表明：变内压定围压测试方法更接近于气田实际压力变化过程；（不）考虑应力敏感试井分析两者在试井特征曲线第II阶段相差较大，五开五关试井分析表明渗透率应力敏感效应的存在对气藏的开发是不利的；应力敏感对单井模拟井底流压影响很大，考虑应力敏感时，气井生产压差明显增大。因此，在低渗透气田开发过程中，必须注意应力敏感对气田开发的影响。     

顶部送风方腔内混合对流换热的数值研究

采用非稳态数学模型对具有对称和非对称结构的顶部送风二维方腔内混合对流流场与温度场进行了数值求解.数值计算中控制参数Ra取为固定值106,Re在1 000~3 000范围内变化.数值结果显示,随Re的增大,具有对称结构问题的数值解会分别出现定常解、周期性振荡解和非周期性振荡解;而对于非对称结构问题,数值解的最大网格Pe虽然大于对称结构问题的最大网格Pe,数值解是定常的,并未发生解的振荡.因此,判明具有对称结构的顶部送风二维方腔内混合对流问题数值解的振荡是客观存在的物理振荡,而非数值方法不稳定所引起的数值振荡.     

数值模式匹配法解波导混合模问题

针对任意非均匀介质填充波导问题，用数值模式匹配法进行了全波分析，给出了波导中电磁场的完整表达式和色散方程。所导出的公式皆为不含任何积分的代数运算，从而具有精度好效率高的优点。最后以部分介质填充波导为例进行了计算，所得结果与文献报道吻合一致。     

误差函数计算方法的研究

误差函数已有多种计算方法，其中按e^-t^2的幂级数展开式为基础的算法，数学上是收敛的．且在科技应用范围内，数值上也是收敛的．数值积分法，如梯形法是计算误差函数更好的方法，文中给出了控制积分变量等分数目的计算公式，并得到了很好的计算结果．     

环保税制下政府规制对企业减排的影响——基于演化博弈的分析

以新环保税制为背景,将环境保护税与政府其他环境规制行为结合起来,运用演化博弈理论,研究地方政府与排污企业的策略互动行为,并将公众举报纳入政府收益函数当中。基于地方政府及排污企业不同策略间净收益的大小,得到12种演化路径,以此为基础将企业减排创新的过程分为发展、成长和成熟期进行数值仿真。结果表明:政府监管强度对企业减排的影响程度视后者减排净收益大小而定,而前者强弱与否主要取决于监管成本的高低;企业减排发展初期过低的环保税率并不会引致企业进行减排技术创新;减排成长期的企业进行减排创新的决策则视政府策略而定;处于成熟期的减排企业则会主动进行减排技术创新。因此,在多种规制工具并行的前提下,政府应视企业所处发展周期的不同制定环保税率,同时辅以一定的税收减免和补贴力度,以降低企业治污成本,促使排污企业在生产过程中向减排技术创新演化。     

调节阀内三维流动与启闭过程的数值模拟及分析

运用计算流体力学方法（CFD）对调节阀的定常与非定常水力特性进行了数值模拟．定常模拟得到各种开度下的三维流场分布以及压降与流量的特性曲线，特性曲线与试验结果相吻合．非定常计算模拟了阀门开启和关闭过程中流动特性的变化，得到流量和阀芯轴向力随时间变化曲线，并对启闭速度的影响进行了分析．     

建筑物绕流的环境效应

本文建立了一个适合建筑物绕流等环境流动的可压缩湍流粘性数值模型，并用该模型模拟了单个和多个建筑物的绕流，得出了建筑物问的滞留旋涡、建筑物“热岛”等对环境有很大影响的流谱．     

移动最小二乘微分求积法解Helmholtz方程

用一种新型的数值方法--移动最小二乘微分求积法（MLSDQ）求解二维Helmholtz方程。MLSDQ方法是一种直接将微分方程离散的方法，它是将未知函数的各阶偏导数在离散点处的值用域内各配点的函数值加权组合来表示，权系数则直接用移动最小二乘Galerkin法中的形函数求导得到，通过MLSDQ技术将Helmholtz方程和相应的边界条件转化成为一组关于各配点位势的线性代数方程组，求解这组代数方程，便可得到各配点的位势，通过求解几个具有精确解的算例，讨论了方法的收敛性和数值精度，结果表明：该方法较适合于求解小波数的Helmholtz方程,对高波数的方程，需要设置大量的域内配点才能有较好的数值结果。