分享到:

超声速平面钝体流动的快速解法

引言 超声速钝体绕流问题作为一个实际应用问题是一个非常老的课题,早在50一60年代,国外就对此进行了广泛的研究,提出了许多方法,如Van Dyke等人的推进解法、直线法、积分关系法等,国内60年代也有许多专家致力于钝体流动的数值解法的研究。但由于它们或出现奇点或因问题的不适定性引起误差增长,所以应用起来存在不少问题。时间相关法的出现及其各种差分格式的提出给跨声速流动的差分计算带来了很大希望,但它漫长的时间渐近过程及其算法构造的复杂往往造成人力物力的巨大浪费。 1979年,More伙i提出的求解非守恒型的非定常欧拉方程组的“几”格式。其重要思路是双曲型方程包含了信息传播的依赖域问题,因此,双曲型方程数值解的差分格式必须反映特征理论的本质。他把非守恒型欧拉方程组的导数项按照特征方向人为分为几个子项,每个子项的差分格式则根据特征值的符号分别采用后差或前差格式。使得欧拉方程的求解无论是精度还是速度都有了大大提高,因而被认为是当时最好的格...  (本文共9页) 阅读全文>>

《高等学校计算数学学报》1991年03期
高等学校计算数学学报

三维涡度方程的谱—差分方法(Ⅰ)

1 Introduction There 15 a lotofliterature eoneerning numeriealmethods for fluid flow,for instanee,fi-nite differenee methodand finite elem血tmethod(e一91,see盛koaehe[1],Riehtmyer,Morton[2],T6man[3]and Guo Ben一yu[4]).During the past ten years,many authors developedsPeetral method and PseudosPeetral method for eomPutationalfluid dynamies(e.g;,seeor、zag[5],Hald[6]and Guo Be。一yu〔71).If the Problems eonsidered are fully Peri...  (本文共12页) 阅读全文>>

《统计与决策》2017年17期
统计与决策

政策效果评估的双重差分方法

0引言如何客观评估政策和制度绩效,特别是定量考察新政策对经济影响的动态因果检验成为经济学界亟需解决的问题。20世纪80年代,国外经济学界借鉴自然科学实验效果检验方法,兴起了一种专门评估政策效果的方法——双重差分法(Differences-in-Differences Method,简称DID)[1,2],由于DID方法思路简洁,模型简单易用,估计方法成熟,在西方学界被广泛应用于政策效果、制度绩效和项目评价等方面,目前国内DID方法应用也呈现快速增长态势。本文对DID方法运用中的问题结合国内外一些代表性论文,围绕如何在自然实验条件下进行规范的实证研究展开。1双重差分方法理论基础与政策效果评估双重差分法源于自然科学比较样本组在实验前后差异的科学传统,具有数百年历史,而西方经济学界引入双重差分法却不到半个世纪,其中原因在于双重差分模型建立条件及其严苛[3]。自然科学实验是建立在控制实验基础上得到实验效果,而社会科学研究几乎不能进行控制...  (本文共4页) 阅读全文>>

《力学学报》1985年06期
力学学报

二维拟线性双曲型方程组的差分方法和检验

一、问题的提出 一阶双曲型方程组常常出现在流体力学里,它描述非定常流动过程和定常超音速流动,而且几乎毫无例外地归结为拟线性方程组和初边值问题,对这类问题的解析求解通常是困难的.有限差分法却是一种广泛使用又行之有效的数值解法比2,,但迄今只对常系数和变系数偏微分方程发展了对数值解法的分析检验方法〔3].对于非线性情形的检验,理论上是不成熟的,必须借助于数值检验.这就需要合适的检验问题.由于人们不仅要了解确切的数值误差,而且要认识误差产生的机理,这就不仅要求检验问题的解具有足够高精度的可运算性和充分的非奇异性,还希望解有显式的表示形式. 本文讨论在x一t平面有限矩形区域a(x镇夕,t。成r《t‘上的拟线性双曲型方程组的初边值问题: u,一a(“,,)u:+b(u,,)。二,。,一c(u,。)u二+d(u,,),二或以矩阵的形式, 。,一、(。)v二,v一{:」和初值 U(x,以及根据矩阵一A(U)的特征值,A(U)=。(u,。)c(...  (本文共10页) 阅读全文>>

《计算机工程与科学》1982年04期
计算机工程与科学

钝体超音速区定常差分方法及其在超大型机上的实施方案

讨1渗~Jl仁1 由于机器容量和速度的限制,不定常差分方法只能用于计算钝休前缘一个不大的混合型区域。要用不定常方法计算整个典型外形的流踢是困难的,而且是不必要的。因为后体的流踢均处于超音速区:其数学模型变为求解双曲型方程的初值问题。它有很多较简便又快速的方法。 本文介绍的定常差分方法是与前绿的非定常差分方法相配套,用于轴对称物体动力计算。它虽没有非定常方法精度好,但程序简单,计算速度快,与非定常方法相比较,计算时间几乎可忽略不计,这样就使这套完整的方法具有很好的实用价值,而被广泛地采用。 本文除了完整地给出这个方法的计算公式和计算思想外,还要给出在超大型机上的实施方案与具有的特点,最后给出初步的效能分析,证明本方法在超大型机上是能高效地运行。二、方程和边界条件 由于方法是把激波和物面当做间断面来处理,面点的方程。 1.内点 考虑光滑凸体定常可压流,其柱坐标系(/,所以必须分别给出内点,激波点和物甲,:)中的基本方程为一知一)~‘...  (本文共7页) 阅读全文>>

《水动力学研究与进展(A辑)》1990年04期
水动力学研究与进展(A辑)

计算管道液压动特性的差分方法

方程组与条件液压动态方程组为1 0.:,,厂.T,万厂‘,犷“二,‘,__二__1,,,.^,、一夕~“一万J、丁宁。不尹1)为(砚1 工T/ 0.。、.,二。欢二~、犷厂曰申才~公少十a一犷苗二U尸其中尸(二,约为流体压强,犷(戈,f)为流速,p为流体密度,a为音速,夕为重力加逮度,管轴与水平面的夹角,f(t+△约为单位长度管内动态摩擦压强损失,良 厂(才+么才、=,旦岑丁即_汽犷(t、+IRr,,.(t+△*、+,,。(,+At、、,,。(,、A,、,飞 J、‘宁。‘,~,几币厂飞。I·。犷、‘少个上OL刀i气‘一。‘少一y忿、不十。i)一夕吕气丁十。T)J少 D““一”一、一‘’一’‘口孟、一’一’‘’汀‘、”一产’汀.、”一’了‘J y:(t+△t)=夕:(t)e一“,’o,,‘,IDZ+扔〔犷(*+△r)一犷(r)〕 夕2(t+△t)=夕:(*)e一“。。,J,尸:+5 .1〔厂(r+△t)一厂(t)〕 夕s(t+A*...  (本文共9页) 阅读全文>>