基于支板燃烧室的喷管化学非平衡效应

采用有限体积法全隐式格式和代用燃料C₁₂H₂₃的10组元13步化学反应Arrhenius有限速率模型研究煤油燃料超燃冲压发动机单边膨胀喷管(SERN, Single Expansion Ramp Nozzle)内的化学非平衡流动,通过建立支板燃烧室——喷管模型有效解决了单边膨胀喷管模型的"入口薄层"问题.计算结果表明,整个单边膨胀喷管内,流动呈现化学非平衡效应,喷管入口附近区域尤其显著;非平衡流动喷管性能明显高于冻结流动,随发动机当量混合比 ε增加,非平衡流动的喷管推力系数和升力系数相对冻结流动的百分比增量 δ不断升高,当 ε=0.8时,推力系数百分比增量 δ_CF达到9.41%,升力系数百分比增量 δ_CY达到16.39%,化学非平衡效应对煤油燃料超燃冲压发动机尾喷管性能的影响不可忽略.     

热气防冰系统内表面弦向传热性能衰减规律

利用反向热流法实验研究了弦向3排射流喷口,且中央喷口正对前缘驻点结构下射流冲击热气防冰系统前缘内表面的局部传热特性。重点研究了喷孔射流出口雷诺数Re_j、笛形管相对于前缘的距离与孔径比Z_n/d、射流冲击驻点区弦向弧长与孔径比r/d对冲击前缘靶面的局部传热性能分布的影响。实验中Re_j范围为2.5×104~1×105,Z_n/d范围为1.736~27.5,r/d范围为13.21~61.5。结果表明,局部传热性能分布曲线为从驻点开始向两侧衰减的钟形曲线,包括稳定段、下降段和结束段3部分,其中稳定段只受参数r/d的影响,而参数Re_j、Z_n/d和r/d均能够明显影响下降段的下降速率,参数r/d则几乎不会改变下降段的总下降幅度。总结了此类喷口结构下射流冲击前缘换热性能衰减分布曲线的通用实验关联式,以指导防冰系统的设计和热性能的评估。      

基于Transition SST模型的高雷诺数圆柱绕流数值研究

为了研究高雷诺数下圆柱绕流边界层的转捩现象和圆柱尾迹近壁区的流动特征,首先通过在典型雷诺数下采用Transition SST四方程转捩模型模拟圆柱绕流得到的结果与实验结果及采用SST k-ω两方程湍流模型模拟结果的对比分析,验证了Transition SST模型在模拟高雷诺数下圆柱绕流的优越性,并较为准确地预测出了圆柱绕流边界层的转捩现象及尾迹近壁区的流动特征。然后分别对亚临界区、临界区、超临界区和过临界区的圆柱绕流问题进行了数值模拟,分析了不同雷诺数下圆柱绕流的流场结构及圆柱表面压力系数、摩擦力系数的变化规律,研究了圆柱绕流近壁区的流动特征、边界层转捩的流动机理、转捩位置及其随雷诺数的变化规律。结果表明,亚临界区二维圆柱绕流边界层发生层流分离,无分离泡和转捩现象;临界区和超临界区二维圆柱绕流边界层先产生了分离泡现象,之后流动发生了转捩并在转捩后发生湍流分离;过临界区二维圆柱绕流边界层流动在转捩之后发生湍流分离,无分离泡现象;在临界区、超临界区和过临界区,二维圆柱绕流边界层转捩位置随雷诺数增大向前驻点移动。     

激发跃迁速率对热力学非平衡氮气紫外辐射的影响

基于氮气的碰撞-辐射（CR）模型,计算了速度为6.2 km/s、初始压力为133 Pa的高超声速流动激波中N₂和N₂+分子电子能级的分布情况,分析了不同激发跃迁速率模型对电子能级分布及辐射光谱模拟的影响。针对流动中热力学非平衡区域和平衡区域,在300~440 nm的辐射光谱开展了逐线法的数值模拟,并与激波管实验测量光谱进行了对比。结果表明,目前的激发跃迁速率均存在偏差,综合Park模型的爱因斯坦系数和Johnston模型的碰撞激发速率可以得到与实验结果最为接近的辐射光谱。      

空间站乘员睡眠区二氧化碳聚集现象

为了分析乘员在没有良好通风的压力舱内停留时,呼出CO₂的聚集过程并评估其可能导致的医学危害,建立了乘员在特定空间站睡眠舱的计算流体动力学模型,对乘员在睡眠区的呼吸过程进行了非稳态模拟.模拟中乘员的CO₂呼出浓度和肺通气量根据生理实验随CO₂吸入浓度变化.通过模拟给出了睡眠区典型位置处CO₂分压随时间的变化规律.最后提出了无通风条件下乘员在睡眠区停留的最大容许时间,睡眠区CO₂浓度监测仪的布置原则和类似情况下为保障乘员生命安全可采取的一些措施.      

改进SA模型对翼型分离流动的数值模拟

采用(CFD,Computational Fluid Dynamics)软件Fluent,选取(SA,Spalart-Allmaras)模型,对风力机翼型S825进行了二维数值模拟,并与实验数据进行比较.针对SA模型不能准确预测翼型尾部流动分离的情况,分析了分离区域内湍流的强非平衡输运特性.研究表明,SA模型中系数C_b1与湍流的非平衡输运特性密切相关,进而提出了修正模型系数C_b1改进SA模型的方法,结果表明该方法能够更准确地模拟风力机翼型分离流动.     

含泡沫镍中间层陶瓷/不锈钢接头微观组织和性能

连接温度850℃,保温时间60 min,用泡沫镍金属作为中间层真空钎焊Al₂O₃陶瓷与1Cr18Ni9Ti不锈钢,利用扫描电子显微镜（SEM）、能谱（EDS）分析钎焊接头的微观组织,利用剪切实验检测接头的力学性能。实验结果表明,不添加中间层时,接头平均剪切强度只有7.7 MPa,断裂位置发生在陶瓷侧;添加泡沫镍中间层时,接头平均剪切强度达到101.7 MPa,断裂位置发生在陶瓷与泡沫镍金属连接界面处。不添加中间层时,Ti元素主要分布在钎料与陶瓷以及钎料与不锈钢反应界面处;添加中间层后钎焊接头中Ti元素主要分布在中间层,与Ni元素形成TiNi₃,Ag、Cu和Ti元素分布更加均匀。      

用时域法EIS技术研究LY12CZ 铝合金的局部腐蚀

通过周期浸泡腐蚀实验和时域法电化学阻抗谱(EIS)技术研究了LY12CZ铝合金的腐蚀损伤,用金相法测量最大腐蚀深度,用统计分析方法处理实验数据.结果表明,LY12铝合金经周期浸泡腐蚀试验后,低频阻抗倒数1/R_d符合Gumbel分布,1/R_d随腐蚀时间变化的动力学曲线呈S形,该曲线分成点蚀和剥蚀两个阶段,点蚀阶段符合Sigmoidal曲线关系,剥蚀阶段符合线性关系.1/R_d与最大腐蚀深度在空间统计分布和腐蚀动力学规律两个方面都有一致性,因此可以用1/R_d表征腐蚀损伤程度.由此可见,时域法EIS技术可作为一种快速、无损、定量的腐蚀评估手段.      

提高k-ωSST模型对翼型失速特性的模拟能力

采用CFD软件fluent,k-ω SST模型对S825翼型进行了二维数值模拟.针对k-ω SST模型对翼型分离失速特性模拟不准的情况,对分离区域内湍流的强非平衡输运特性进行了分析.研究表明:原始模型由于不能准确模拟分离区内湍流的强非平衡输运特性,导致其对翼型失速特性模拟失效;提出了通过修正模型系数a₁及β_∞^*,提高k-ω SST模型预测非平衡湍流输运特性的方法,从而提高对翼型失速特性的模拟精度.     

基于伪寿命的加速退化机理一致性边界检验

根据寿命型随机变量的加速机理一致性条件,提出一种基于伪寿命的加速退化机理一致性边界检验方法,通过由低到高逐一对各个加速应力水平下的退化机理一致性进行检验,来确定加速退化机理一致性的边界应力水平.针对工程常用的两类寿命型分布——对数正态分布和Weibull分布,分别以对数标准差和形状参数为退化机理表征量,在产品伪寿命估计的基础上建立了退化机理表征量齐性检验的F统计量,导出了退化机理一致性检验拒绝域的解析式,给出了退化机理一致性边界检验的具体步骤.该方法能够充分利用多个应力水平下表征产品退化机理一致的横向信息,有效提高了检验的精度.以某镀膜平板玻璃热退化机理一致性边界温度的确定为例,分析说明了该方法的合理性与有效性.     

Geomagnetic effects on mid-latitude railways: A statistical study of anomalies in the operation of signaling and train control equipment on the East-Siberian Railway

Intense geomagnetically induced currents (GIC) can hamper rail traffic by disturbing signaling and train control systems. GIC threats have been a concern for technological systems at high-latitude locations due to geomagnetic disturbances driven by substorm expansion electrojet or convection electrojet intensifications. However, other geomagnetic storm processes such as storm sudden commencement (SSC) and geomagnetic pulsations can also cause GIC concerns for technological systems. We present in this paper the first evidence based on statistical data for links between geomagnetic disturbances and faulty operations (anomalies) in the functioning of railway automatics and telemetry. We analyze anomalies of automatic signaling and train control equipment which occurred in 2004 on the East-Siberian Railway (corrected geomagnetic latitude _m = 46–51°N and longitude λ_m = 168–187°E). Our results reveal a seasonal effect in the number of anomalies per train similar to the one observed in geomagnetic activity (Kp, Ap, Dst indices). We also found an increase by a factor of 3 in the total duration of daily anomalies during intense geomagnetic storms (local geomagnetic index specific to Siberian Observatory A_max > 30), with a significant correlation between the daily sum of durations of anomalies with geomagnetic activity. Special attention was paid to failures not related to recognized technical malfunctions. We found that the probability of these failures occurring in geomagnetically disturbed periods was 5–7 times higher than the average anomaly occurrence.     

发动机进气道水滴撞击特性分析

针对某发动机进气道的防冰系统设计,对进气道的水滴撞击特性进行了研究.在对进气道流场进行计算的基础上,采用差分法对水滴运动方程进行了数值求解,得到了水滴运动轨迹,从而确定了水滴的撞击极限、总收集系数和局部水收集系数等水滴撞击特性参数;此外,还研究了飞行高度、飞行速度及水滴半径对水滴撞击特性的影响;研究发现水滴撞击极限、总收集系数和局部水收集系数随飞行速度的减小、飞行高度的增加或者水滴半径的增加而增大.这些为进气道的防冰系统设计奠定了基础.     

PANS方法在双圆柱绕流数值模拟中的性能分析

PANS方法以传统的雷诺平均Navier-Stokes(RANS)方程为模板,引入控制参数对原方程进行修改,是一类能够有效模拟分离湍流流动的混合方法.采用基于Menter 剪切应力输运(SST)湍流模型的PANS方法,选取双圆柱绕流作为研究对象, 分析双圆柱流动结构,验证方法的准确性与可靠性,考察PANS方法中模型参数选取对于计算结果的影响,评估不同类型PANS方法的性能.研究表明:模化湍动能比例这一模型参数对计算结果有显著影响,且计算域内统一的模型参数取值难以处理复杂流动问题.在此基础上使方法中模化湍动能比例取值随流场物理信息和网格尺度变化,所得结果与SST DES计算结果及实验结果吻合较好,表明该方法适用于复杂湍流数值模拟.      

Human adaptation to simulated gravitational fields

We present the results of manned studies in which test subjects were exposed to simulated zero $\text{g}$ (water immersion or head-down tilt at ?6°) and head-to-feet acceleration. The findings give evidence that humans have different individual tolerances to an acceleration of +3 G_z after exposure to zero $\text{g}$, whether simulated by immersion or by head-down tilt. The paper discusses the functional relationship between water balance and cardiac output in the establishment of adaptive reactions to simulated zero $\text{g}$.     

Hierarchical global search for fuel-optimal impulsive transfers between lunar libration orbits

The attention to the periodic orbit in the Earth-Moon restricted three-body system continues to grow due to its special environment and locations. This research investigates the feasibility of constructing fuel-optimal single and multiple impulse transfers between unstable periodic orbits at L₁ and L₂ points. Invariant manifolds, which could provide the appropriate initial trajectories for optimization, are analyzed deeply to enable previously unknown orbit options and potentially to reduce mission cost. A global search strategy based on comparing the orbital state of the unstable and stable manifolds, incorporated with low-thrust techniques, is performed to seek a suitable matching point for maneuver application. Then the sequential quadratic programming (SQP) is adopted to further optimize the velocity increment and obtain the single/multiple impulse optimal transfers. The associated constraint gradients are derived to achieve higher accuracy and rapidity of the algorithm. To highlight the effectivity of the transfer scheme, three-dimensional low-energy transfers between different types and spatial regions of performing single and multiple impulses are explored. The total Delta-V required varies between a few meters per second and tens of meters per second, and the related flight time is about several weeks, mainly depending on the energy of periodic orbits and the invariant manifold structure. The results obtained in this paper can provide a useful reference for the selection of escape and capture site along the manifolds, maneuver magnitude and transfer time.     

Some effects in the upper atmosphere during geomagnetic storms

This paper examines the response of the high latitude ionosphere–thermosphere system during two intense geomagnetic storms. For that, data taken by instruments on board Dynamic Explorer 2 taken at heights of the F2-layer are used. These results represent a comparison of simultaneous measurements of storm disturbances in gas composition, electron density and temperature in common local time sectors. Documented are an increase in electron temperature and a decrease in electron density; increases both in electron temperature and electron density; and the correlation between electron density decreases and increases in the ratio N₂/O. It is noticed that the decrease in electron density is sometimes due to an increase in the molecular nitrogen density N₂ and not always is attributed both to the increase in N₂ density and the simultaneous decrease in the atomic oxygen density.     

A linear stability analysis of large-Prandtl-number thermocapillary liquid bridges

A linear stability analysis is applied to determine the onset of oscillatory thermocapillary convection in cylindrical liquid bridges of large Prandtl numbers (4  Pr  50). We focus on the relationships between the critical Reynolds number Re_c, the azimuthal wave number m, the aspect ratio Γ and the Prandtl number Pr. A detailed Re_c–Pr stability diagram is given for liquid bridges with various Γ. In the region of Pr > 1, which has been less studied previously and where Re_c has been usually believed to decrease with the increase of Pr, we found Re_c exhibits an early increase for liquid bridges with Γ around one. From the computed surface temperature gradient, it is concluded that the boundary layers developed at both solid ends of liquid bridges strengthen the stability of basic axisymmetric thermocapillary convection at large Prandtl number, and that the stability property of the basic flow is determined by the “effective” part of liquid bridge.     

Stellar β-decay rates of iron isotopes and its implications in astrophysics

β-Decay and positron decay are believed to play a consequential role during the late phases of stellar evolution of a massive star culminating in a supernova explosion. The β-decay contributes in maintaining a respectable lepton-to-baryon ratio, Y_e, of the core prior to collapse which results in a larger shock energy to produce the explosion. The positron decay acts in the opposite direction and tends to decrease the ratio. The structure of the presupernova star is altered both by the changes in Y_e and the entropy of the core material. Recently the microscopic calculation of weak interaction mediated rates on key isotopes of iron was introduced using the proton–neutron quasiparticle random phase approximation (pn-QRPA) theory with improved model parameters. Here I discuss in detail the improved calculation of β^±-decay rates for iron isotopes (^54,55,56Fe) in stellar environment. The pn-QRPA theory allows a microscopic “state-by-state” calculation of stellar rates as explained later in text. Excited state Gamow–Teller distributions are much different from ground state and a microscopic calculation of decay rates from these excited states greatly increases the reliability of the total decay rate calculation specially during the late stages of stellar evolution. The reported decay rates are also compared with earlier calculations. The positron decay rates are in reasonable agreement with the large-scale shell model calculation. The main finding of this work includes that the stellar β-decay rates of ^54,55,56Fe are around 3–5 orders of magnitude smaller than previously assumed and hence irrelevant for the determination of the evolution of Y_e during the presupernova phase of massive stars. The current work discourages the inclusion of ^55,56Fe in the list of key stellar β-decay nuclei as suggested by former simulation results.