适用于气动弹性的小波非定常气动力降阶方法

发展非定常气动力模型降阶技术旨在缩减计算耗费并且使得计算流体力学信息能够应用于气动伺服弹性以及设计优化当中。采用小波方法建立基于Volterra级数的非定常气动力降阶模型。在模型识别过程中,激励是光滑连续的且以零阶保持方式进行离散,采样频率选择二的整数幂次。近似的一阶Volterra核基于Haar尺度函数族展开,鉴于Volterra核在系统响应中的衰减特性,可在合适的有限时长截断一阶核。为了获得一阶核的各项展开系数,需要求解由输入/输出数据组成的超定方程,其中涉及到奇异值分解算法。为了验证小波方法的有效性,算例选取了二维的NACA64a010翼型。数值仿真结果表明该方法能够比较准确地预测结构小扰动引起的非定常气动力响应且能描述一定的非线性现象。     

基于奇偶方程的FADS传感器故障检测方法

重复使用运载器（RLV）的嵌入式大气数据传感系统（FADS）中传感器的高可靠性是RLV飞行控制系统高可靠性的保障。结合FADS采用多个测压点冗余配置的特点,利用各传感器测量值之间存在的解析冗余关系,设计奇偶方程,实现对各个测压点故障传感器的有效检测。     

外涵偏轴分开排气喷管的流场和声场数值计算

以某型涡扇发动机为原型,建立了1/10缩比的分开排气式涡扇发动机排气喷管物理模型,采用大涡模拟(LES)计算了不同外涵偏轴距离的排气喷管非稳态流场,利用FW-H(Ffows Williams-Hawkings)声学模型和傅里叶变换得到了排气喷流噪声声压级空间分布特征.计算结果表明:外涵偏轴式分开排气喷管有抑制排气喷流噪声的作用.喷流噪声整体抑制效果与偏轴距离有关,在偏轴距离为0.15D时,排气喷管下方的噪声降低达到最大,且排气噪声峰值由原型的146.7dB降低到139.5dB,降低了7.2dB.研究的排气结构和数值计算方法可作为分开排气式涡扇发动机降噪设计的参考.      

单圆盘转子非线性振动模型及机理

根据转子非线性振动的力学性质及大变形假设,结合工程实际构建了转子非线性振动具有Duffing方程形式的力学模型,推导了非线性系统的等效频率和频响方程,研究了刚度元件、阻尼元件与非线性元件对系统振动的影响.通过单圆盘转子非线性振动实验,验证了非线性振动模型结构与Duffing方程求解及频响方程的合理性与有效性,揭示了共振、突变、滞后等非线性现象的机理,为转子振动特性分析提供了一种新的研究方法.      

含裂纹固体推进剂装药的内流场数值模拟

利用“N－S方程＋粒子仿真”方法来计算含裂纹翼柱型装药内流场。计算结果表明：环形裂纹内的流动对主流的影响大于纵向裂纹，环形裂纹向上游位置流动对主流的影响较下游位置的大，并可保持到喷管喉部，翼槽附近环形裂纹内的流动对主流的影响相对较弱。     

高超声速内收缩进气道分步优化设计方法

提出了基准流场与唇口平面形状分步优化的高超声速内收缩进气道设计方法。基准流场以反射激波不均匀性最小和总压恢复最大进行多目标优化设计,使用结合Tayler-Maccoll方程的有旋特征线方法(MOC)进行流场计算,获得双拐点母线内收缩锥基准流场。进气道唇口形状以沿流线积分(Streamline Integral Method, SIM)获得的进气道无黏阻力最小为目标进行优化设计,获得类椭圆形唇口平面形状。针对优化设计结果进行数值模拟,与传统直母线基准流场相比,双拐点母线基准流场反射激波后流动不均匀性下降40%左右,总压损失减少35%左右,总体性能提升明显。类椭圆唇口进气道在设计点的单位质量流量无黏阻力相较于圆形唇口降低6%,具有良好的压缩特性和气动效率,能够减弱进气系统对飞行器气动性能的不利影响。研究结果表明该方法是一种高效且实用的高超声速内收缩进气道设计方法。     

Evolution and History in a new “Mathematical SETI” model

In a recent paper (Maccone, 2011 [15]) and in a recent book (Maccone, 2012 [17]), this author proposed a new mathematical model capable of merging SETI and Darwinian Evolution into a single mathematical scheme. This model is based on exponentials and lognormal probability distributions, called “b-lognormals” if they start at any positive time b (“birth”) larger than zero. Indeed:

• 1.Darwinian evolution theory may be regarded as a part of SETI theory in that the factor f_l in the Drake equation represents the fraction of planets suitable for life on which life actually arose, as it happened on Earth.
• 2.In 2008 (Maccone, 2008 [9]) this author firstly provided a statistical generalization of the Drake equation where the number N of communicating ET civilizations in the Galaxy was shown to follow the lognormal probability distribution. This fact is a consequence of the Central Limit Theorem (CLT) of Statistics, stating that the product of a number of independent random variables whose probability densities are unknown and independent of each other approached the lognormal distribution if the number of factors is increased at will, i.e. it approaches infinity.
• 3.Also, in Maccone (2011 [15]), it was shown that the exponential growth of the number of species typical of Darwinian Evolution may be regarded as the geometric locus of the peaks of a one-parameter family of b-lognormal distributions constrained between the time axis and the exponential growth curve. This was a brand-new result. And one more new and far-reaching idea was to define Darwinian Evolution as a particular realization of a stochastic process called Geometric Brownian Motion (GBM) having the above exponential as its own mean value curve.
• 4.The b-lognormals may be also be interpreted as the lifespan of any living being, let it be a cell, or an animal, a plant, a human, or even the historic lifetime of any civilization. In Maccone, (2012 [17, Chapters 6, 7, 8 and 11]), as well as in the present paper, we give important exact equations yielding the b-lognormal when its birth time, senility-time (descending inflexion point) and death time (where the tangent at senility intercepts the time axis) are known. These also are brand-new results. In particular, the σ=1 b-lognormals are shown to be related to the golden ratio, so famous in the arts and in architecture, and these special b-lognormals we call “golden b-lognormals”.
• 5.Applying this new mathematical apparatus to Human History leads to the discovery of the exponential trend of progress between Ancient Greece and the current USA Empire as the envelope of the b-lognormals of all Western Civilizations over a period of 2500 years.
• 6.We then invoke Shannon's Information Theory. The entropy of the obtained b-lognormals turns out to be the index of “development level” reached by each historic civilization. As a consequence, we get a numerical estimate of the entropy difference (i.e. the difference in the evolution levels) between any two civilizations. In particular, this was the case when Spaniards first met with Aztecs in 1519, and we find the relevant entropy difference between Spaniards an Aztecs to be 3.84 bits/individual over a period of about 50 centuries of technological difference. In a similar calculation, the entropy difference between the first living organism on Earth (RNA?) and Humans turns out to equal 25.57 bits/individual over a period of 3.5 billion years of Darwinian Evolution.
• 7.Finally, we extrapolate our exponentials into the future, which is of course arbitrary, but is the best Humans can do before they get in touch with any alien civilization. The results are appalling: the entropy difference between aliens 1 million years more advanced than Humans is of the order of 1000 bits/individual, while 10,000 bits/individual would be requested to any Civilization wishing to colonize the whole Galaxy (Fermi Paradox).
• 8.In conclusion, we have derived a mathematical model capable of estimating how much more advanced than humans an alien civilization will be when SETI succeeds.

     

剪切气流中无黏液体横向射流破碎机理

航空发动机燃烧室中,液体燃料的雾化过程、特别是初始雾化过程非常复杂,至今无法建立准确的初始雾化模型。忽略液体射流黏性,采用线性不稳定性分析法对无黏液体（工质为水）在二维剪切横向气流中的破碎机理进行研究。通过建立射流的色散方程,根据其表面波的增长率及波数的发展情况对射流破碎进行预测。当气体韦伯数或液体韦伯数增大时,射流表面波增长率显著增加,最佳波长明显减小。液气动量比大于临界值时,Kelvin-Helmholtz （K-H）不稳定性占主导作用;反之,Rayleigh-Taylor （R-T）不稳定性占主导作用。二维剪切气流在射流方向上具有速度梯度,只改变横向气动力对射流表面波的作用。气体韦伯数与液体韦伯数对射流破碎的作用与均匀气流相似;保持气流流量相同时,负梯度剪切气流可以加速射流的破碎。     

奇异积分方程方法在接触压力分析中的应用

总结回顾了奇异积分方程的数值解法,对于第二类奇异积分方程使用分段连续函数法进行了求解;对于两个弹性体接触问题,通过接触体之间的滑移函数和间隙函数,建立求解接触压力的奇异积分方程.分别针对圆柱体与弹性半空间体、抛物线型压头与弹性半空间体、榫头与榫槽3类接触问题,确立奇异积分方程的具体表达式,而后使用分段连续函数方法进行求解,获得接触面上的接触压力.最后将计算所得的接触压力分别与理论解和有限元解进行了对比.对于圆柱体与弹性半空间体接触问题,奇异积分方程法的最大接触压力与理论解和有限元解的相对误差分别为0.3%和0.5%;对于榫头与榫槽接触问题,奇异积分方程方法计算所得的最大接触压力与有限元解的相对误差为1.8%,验证了奇异积分方程方法的有效性.      

微纳卫星共面伴飞相对运动椭圆短半轴最省燃料控制

针对伴随微纳卫星资源受限、轨控需要尽可能节省燃料的现实问题,基于希尔（Hill）方程,研究推导了共面编队伴飞卫星的轨控时机和轨控方向对相对运动椭圆短半轴控制效率的影响。理论推导和仿真均表明：当控制量大小｜ΔV｜与相对运动椭圆短半轴b满足｜ΔV｜≤nb/2关系时（n为参考星平均轨道角速度）,在相对运动椭圆上下点进行横向或反横向控制,最大效率地将相对运动椭圆短半轴改变了｜Δb｜=2｜ΔV｜/n。其中,在上点反横向或下点横向进行控制,可以最大效率地增大椭圆短半轴;在上点横向或下点反横向进行控制,可以最大效率地减小椭圆短半轴。