虚构解方法程序验证

对求解偏微分方程问题的程序的严格验证,一直以来由于其精确解有限,较难进行,针对这种情况,本文给出了一种新的程序验证方法,虚构解方法,该方法旨在解决大型科学计算程序如何较严格地进行程序验证这一问题;该方法通过构造虚构解,修改原控制方程,然后通过对计算结果的分析处理达到对程序的精确验证;文中给出了一组通用的虚构解的选取办法,解决了虚构解方法验证过程中确定虚构解的问题;同时还给出了运用该方法对非结构二维Euler计算程序进行验证的具体例子,网格收敛分析结果表明该方法是一种有效的程序验证方法,并且具有较好的通用性。     

大迎角非线性飞行动力学（二）

五、大迎角飞行时动态特性研究方法飞机在大迎角飞行时,是一个复杂的非线性系统。非线性系统的研究与线性系统有很大的区别,它不像线性系统那样有一般性的有效的求解方法。由于系统非线性的形式是多种多样的,因此研究的方法也有许多种,其理论研究方法大致可以分为二大类:解析的方法和几何的(或拓扑的)方法。对由非线性常微分方程所描述的非线性动力学系统,其解析法通常可以分成三种类型:“精确的”、“时间连续的”和“平均的”。存在有解析的、封闭形式的精确解的问题是极少的。时间连续法,众所皆知的Runge-Kutta数值积分法就是代表方法之一,其计算程序     

解不可压N—S方程的PCIMG法

本文提出了一种解不可压N-S方程的预测、修正、迭代、多重网格法。这种算法程序简单,结果精确及收敛快速。     

后掠飞机全机自由自由振动的矩阵初参数解法

本文研究如何应用矩阵初参数法求解后掠飞机全机自由自由振动问题,指出了几种解法。对其中一种解法导出全部矩阵公式,设计出电子计算机的程序。附有计算例题, 与精确解及试验数据作了比较。悬臂机翼及导弹弹身的振动问题只是本问题的特例。     

亚临界情况下翼型的非线性压缩性修正

本文从有限基本解方法出发,介绍一种计算亚音速流中非线性压缩性效应的方法。采用逐步逼近的办法,可以把有限基本解方法推广到求解非线性的无粘可压流体,从而为非线性方程的求解提供了新的途径。 编制了FORTRAN语言程序。与精确解和其他理论值的比较表明,结果比传统的有限差分法好,尤其是在翼型头部。     

三维跨声速多级压气机内确定应力耦合问题研究

为开展确定应力耦合效果的研究,基于三维定常Denton程序发展了三维非定常计算程序,以多级跨声速压气机近失速点为算例,进行确定应力耦合精度的研究。除定常计算和非定常计算外,还采用了三种计算方式,分别是多级环境下非定常时均方程耦合确定应力的解(COU2),提取非定常精确解时均的值作为S2进口边界条件进行单独S2计算域的定常计算(ST1);在ST1基础上耦合确定应力的解(COU1)。结果表明:三维压气机算例中出口总压和马赫数的径向曲线对比显示COU1的误差水平分别为ST1的53%和34%;验证了确定应力项的耦合可以使定常计算的结果更接近非定常计算的时均解;对比发现ST1计算精度较COU2高,误差水平分别降低至37%和32%。说明在该算例中精确的边界条件比精确的确定应力发挥更大的作用。非定常时均结果中近压力面的熵值比近吸力面的熵值要高,反映出非定常相互作用中上游尾迹向压力面迁移的效应,Adamczyk时均方程耦合确定应力后得到了类似的结果。     

某型飞机燃油计算软件算法研究

实时准确地测量油箱中的剩余油量可以精确计算飞机续航时间,保证飞机安全飞行.本文对某型飞机燃油计算软件算法进行研究,提出了一种稳定、可靠、准确的燃油计算方法,该算法包含油面高度和油面角的计算、燃油质量特性数据库查询、差值解算和信息后处理,可精确得出各油箱内的燃油油量.     

太阳辐射对ER3BP拉格朗日点的影响

为了研究Halo轨道中继通信以及行星际轨道转移等问题,需要首先对拉格朗日点(平动点)的理论进行分析。使用考虑太阳辐射的椭圆型限制性三体建模,代替原有的圆型限制性三体建模,提高了建模的精度。提出了在L1,L2,L3点的更为精确的新解析解,经过Matlab仿真,其相对于精确数值解的误差,分别为Ammar的解析解的13.13%,29.51%,0.46%。     

相对温升法测量导温系数的热损效应

本文从 Cowan 方程出发,分析了热损对温升法测量导温系数的影响。结果发现,当按本文所述原则选取无量纲时间ω时,可得到不受热损影响的导温系数精确解。     

航天飞机再入大气轨道分析

本文根据航天飞机再入大气轨道的运动方程,求出了升阻比为常值条件下的数值解。着重研究了改变升阻比和再入角对再入轨道的影响。并将精确数值解与近似解进行比较。此外,文中还估算了作用于航天飞机的最大热流。     

跨音速压气机动静叶排相干的三维非定常流动数值分析

采用时空二阶精度ＮＮＤ显式差分格式求解三维Ｅｕｌｅｒ方程，对跨音速压气机级中动静叶排相干形成的非定常流动进行了数值分析。除动静叶排流动的控制方程体系、数值方法外，着重讨论了动静叶排计算域、网格以及边界条件的处理。对某压气机第１级动静叶排相干的非定常流场的计算，给出了叶片表面的马赫数分布、压力波动幅值，以及叶片槽道中瞬态密度等值线等，并分析了结果的合理性。     

二维轴流压气机转—静非定常Euler方程数值计算

将Euler方程写在包括转子和静子的绝对坐标系下,通过在转—静之间设计一个可以随时间扭动的平行四边形网格计算了二维轴流压气机转—静非定常流场。由于整个方程统一写在动网格下,转—静之间没有使用插值。      

用多重网格方法计算旋翼跨声速无粘流场

发展了一种加快悬停旋翼无粘流场计算收敛速度的多重网格方法。由于悬停旋翼流场中存在不可压区域，同时旋翼尾涡系统的发展需要较长的时间，使得旋翼流场的收敛速度远低于固定翼流场，因此研究旋翼流场的多重网格算法具有重要意义。空间离散格式采用了中心有限体积方法，时间推进应用了五步龙格-库塔法。采用3层网格的V循环，对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明：尽管多重网格方法对旋翼流场的加速收敛作用不如对固定翼流场的加速收敛效果，但是多重网格方法仍然可以显著地加快旋翼流场收敛。     

跨音速压气机级三维紊流流场数值模拟

在对跨音速压气机级的动／静双叶排中的三维定常紊流流场数值的模拟中，利用三阶高分辨率ＮＮＤ格式和ＬＵ－ＳＧＳ隐式推进迭代法，既保证了流场中激波的模拟质量和粘性流动特征的正确预估，也实现了求解过程的高效率。对某单级跨音速压气机的计算及其与实验数据的对比证实了本文方法的有效性。将计算推广至多级轴流叶轮机械的情况将会是十分简便的     

总压畸变对风扇流场影响的全流道数值研究

采用全流道非定常数值模拟方法对进口稳态总压畸变条件下跨声速风扇流场进行求解.详细分析了总压畸变造成风扇性能下降的原因.分析表明,周向不均匀流场使得动静叶进口气流角变得不均匀,从而使进口攻角发生很大变化;动叶在不同周向和径向位置流动状况完全不同;畸变区内静叶根部至中径附近存在大范围的分离区.      

A New Single-blade Based Hybrid CFD Method for Hovering and Forward-flight Rotor Computation

A hybrid Euler/full potential/Lagrangian wake method,based on single-blade simulation,for predicting unsteady aerodynamic flow around helicopter rotors in hover and forward flight has been developed.In this method,an Euler solver is used to model the near wake evolution and transonic flow phenomena in the vicinity of the blade,and a full potential equation(FPE) is used to model the isentropic potential flow region far away from the rotor,while the wake effects of other blades and the far wake are incorporated into the flow solution as an induced inflow distribution using a Lagrangian based wake analysis.To further reduce the execution time,the computational fluid dynamics(CFD) solution and rotor wake analysis(including induced velocity up-date) are conducted parallelly,and a load balancing strategy is employed to account for the information exchange between two solvers.By the developed method,several hover and forward-flight cases on Caradonna-Tung and Helishape 7A rotors are per-formed.Good agreements of the loadings on blade surface with available measured data demonstrate the validation of the method.Also,the CPU time required for different computation runs is compared in the paper,and the results show that the pre-sent hybrid method is superior to conventional CFD method in time cost,and will be more efficient with the number of blades increasing.     

高负荷涡轮端区非定常流动相互作用研究

采用三维黏性非定常数值模拟方法研究了某型高压涡轮端区非定常流动相互作用,着重研究了上游静子尾迹与转子二次流的非定常作用机制,同时还分析了负荷分布、激波等对端区非定常流动的影响。结果表明,静子尾迹的非定常作用一定程度减小了转子轮毂二次流的径向涡量;尾迹对流向涡量的影响取决于尾迹沿叶高的分布,当吸力面一侧的尾迹具有与二次流方向相反的流向涡量时,二次流的流向涡量减小;非定常效应还使得转子叶片根部负荷略为减小,也一定程度上抑制了转子轮毂二次流的发展。此外,受静子尾缘激波的影响,转子叶片表面负荷分布发生明显的周期性变化,导致叶片表面承受较强的非定常力,在涡轮设计中必须考虑。另外,通过计算涡轮级中的熵增和熵产,定量地分析了端区非定常相互作用产生的损失,并得到了一些有意义的结论。     

压气机叶片排中三维流场计算的多网格法

以三维有限差分形式的MacCormack格式为基础, 对相对圆柱坐标系下弱守恒型控制方程, 导出了残差在粗、细网格上的传播与回插公式, 形成了一种多网格算法。计算表明, 这种算法大幅度地提高了计算收敛速度, 而且方法本身简便易行。      

先进直升机旋翼悬停状态气动性能计算

 为更好地模拟先进直升机旋翼流场,同时准确计算其悬停状态的气动性能,建立了一套基于Navier-Stokes/Euler方程的混合CFD方法。该方法的求解域由两部分组成：一是围绕旋翼桨叶周围的黏性区域,采用可压缩Navier-Stokes方程来模拟旋翼附近的黏性流动和近场尾涡的捕捉;二是离桨叶较远、黏性可以忽略的远场区域,用Euler方程来描述其流动。在该方法中,将三阶逆风格式（MUSCL）与通量差分分裂方法相结合,无需添加人工黏性,因而可有效地减少旋翼尾迹数值耗散。?阌诹鞒》智蠼庖约爸芷谛员呓缣跫氖凳?采用了嵌套网格方法,并给出旋翼网格与背景网格交界面的信息传递方式。应用所建立的计算方法,首先对二维翼型、三维M6机翼的流场进行了数值模拟,以验证计算方法;然后,着重计算了有实验结果可供对比的具有先进气动外形的HELISHAPE 7A模型旋翼和UH-60A直升机旋翼,通过计算得到旋翼表面压力分布、桨叶展向拉力系数分布、桨叶表面细节流动以及气动性能等,进一步验证了该方法的有效性。     

轴流涡轮基元级动叶稠度特性的数值研究

利用商用软件数值模拟了5个不同动叶稠度的轴流涡轮基元级的非定常流动情况,以研究动叶稠度对轴流涡轮基元级性能和流动情况的影响.通过对动叶稠度对基元级反力度、叶片进出口气流角、转子和静子中的流场及损失影响情况的考察研究,发现动叶稠度的改变对涡轮基元性能和流动情况的影响与静叶稠度存在重要关系.静叶稠度不变时,动叶稠度的改变通过影响流过涡轮基元级的流量来使基元级的反力度发生变化.当动叶稠度过大时,气流在转子中会过度膨胀加速而产生激波损失及其与附面层干涉形成的流动分离损失.动叶稠度过小时,转子进口会出现极大的正攻角致使动叶吸力面发生大范围的流动分离.静叶稠度一定时,存在一个最佳的动叶稠度,使涡轮基元级呈现最好的性能.