加热圆管流动热阻力及熵产的数值分析

本文利用数值计算方法,对圆管加热流动的热阻力和熵产进行了综合分析研究。以空气加热圆管为例,通过联立求解流体力学基本方程和熵产方程,获得了热阻力、熵产及其随加热量、管长而变化的关系。数值分析结果表明,加热流动中的热阻力与熵产二者之间存在着对应相关的关系。     

Measuring atmospheric density using GPS–LEO tracking data

We present a method to estimate the total neutral atmospheric density from precise orbit determination of Low Earth Orbit (LEO) satellites. We derive the total atmospheric density by determining the drag force acting on the LEOs through centimeter-level reduced-dynamic precise orbit determination (POD) using onboard Global Positioning System (GPS) tracking data. The precision of the estimated drag accelerations is assessed using various metrics, including differences between estimated along-track accelerations from consecutive 30-h POD solutions which overlap by 6 h, comparison of the resulting accelerations with accelerometer measurements, and comparison against an existing atmospheric density model, DTM-2000. We apply the method to GPS tracking data from CHAMP, GRACE, SAC-C, Jason-2, TerraSAR-X and COSMIC satellites, spanning 12 years (2001–2012) and covering orbital heights from 400 km to 1300 km. Errors in the estimates, including those introduced by deficiencies in other modeled forces (such as solar radiation pressure and Earth radiation pressure), are evaluated and the signal and noise levels for each satellite are analyzed. The estimated density data from CHAMP, GRACE, SAC-C and TerraSAR-X are identified as having high signal and low noise levels. These data all have high correlations with anominal atmospheric density model and show common features in relative residuals with respect to the nominal model in related parameter space. On the contrary, the estimated density data from COSMIC and Jason-2 show errors larger than the actual signal at corresponding altitudes thus having little practical value for this study. The results demonstrate that this method is applicable to data from a variety of missions and can provide useful total neutral density measurements for atmospheric study up to altitude as high as 715 km, with precision and resolution between those derived from traditional special orbital perturbation analysis and those obtained from onboard accelerometers.     

Debris flux comparisons from the Goldstone Radar,Haystack Radar,and Hax Radar prior,during, and after the last solar maximum

The continual monitoring of the low Earth orbit (LEO) debris environment using highly sensitive radars is essential for an accurate characterization of these dynamic populations. Debris populations are continually evolving since there are new debris sources, previously unrecognized debris sources, and debris loss mechanisms that are dependent on the dynamic space environment. Such radar data are used to supplement, update, and validate existing orbital debris models. NASA has been utilizing radar observations of the debris environment for over a decade from three complementary radars: the NASA JPL Goldstone radar, the MIT Lincoln Laboratory (MIT/LL) Long Range Imaging Radar (known as the Haystack radar), and the MIT/LL Haystack Auxiliary radar (HAX). All of these systems are highly sensitive radars that operate in a fixed staring mode to statistically sample orbital debris in the LEO environment. Each of these radars is ideally suited to measure debris within a specific size region. The Goldstone radar generally observes objects with sizes from 2 mm to 1 cm. The Haystack radar generally measures from 5 mm to several meters. The HAX radar generally measures from 2 cm to several meters. These overlapping size regions allow a continuous measurement of cumulative debris flux versus diameter from 2 mm to several meters for a given altitude window. This is demonstrated for all three radars by comparing the debris flux versus diameter over 200 km altitude windows for 3 nonconsecutive years from 1998 to 2003. These years correspond to periods before, during, and after the peak of the last solar cycle. Comparing the year to year flux from Haystack for each of these altitude regions indicate statistically significant changes in subsets of the debris populations. Potential causes of these changes are discussed. These analysis results include error bars that represent statistical sampling errors.     

TRIP2.0软件的确认：DPWⅡ复杂组合体的数值模拟

 采用“亚跨超声速计算流体力学软件平台”（TRIP2.0）数值模拟了阻力预测小组（AIAA CFD Drag Prediction Workshop Ⅱ,DPWⅡ）翼/身/架/舱复杂组合体运输机构型,数值模拟采用的多块对接网格、测压和测力的试验结果均来自DPWⅡ,对比计算采用了CFL3D的结果。重点针对DLR-F6翼/身/架/舱复杂组合体构型,详细研究了网格密度和湍流模型对总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致。采用SA一方程和SST两方程模型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布有一定的影响。     

Use of two-line element data for thermosphere neutral density model calibration

Traditional empirical thermospheric density models are widely used in orbit determination and prediction of low-Earth satellites. Unfortunately, these models often exhibit large density errors of up to around 30% RMS. Density errors translate into orbit errors, adversely affecting applications such as re-entry operations, manoeuvre planning, collision avoidance and precise orbit determination for geodetic missions. The extensive database of two-line element (TLE) orbit data contains a wealth of information on satellite drag, at a sufficiently high spatial and temporal resolution to allow a calibration of existing neutral density models with a latency of one to two days. In our calibration software, new TLE data for selected objects is converted to satellite drag data on a daily basis. The resulting drag data is then used in a daily adjustment of density model calibration parameters, which modify the output of an existing empirical density model with the aim of increasing its accuracy. Two different calibration schemes have been tested using TLE data for about 50 objects during the year 2000. The schemes involve either height-dependent scale factors to the density or corrections to CIRA-72 model temperatures, which affect the density output based on a physical model. Both schemes have been applied with different spherical harmonic expansions of the parameters in latitude and local solar time. Five TLE objects, varying in perigee altitude between 280 and 530 km, were deliberately not used during calibration, in order to provide independent validation. Even with a single daily parameter, the RMS density model error along their tracks can already be reduced from the 30% to the 15% level. Adding additional parameters results in RMS errors lower than 12%.     

对称面圆周角对轴对称降落伞流场特性的影响

根据降落伞的特点,通过3点假设(伞衣薄膜、伞的轴对称和流场定常),将三维复杂流动问题转化成二维轴对称问题,以节约计算时间.定义对称面圆周角,保持伞衣幅底部直径和顶孔直径不变,选取对称面圆周角在80°~140°范围内变化,建立一系列轴对称降落伞的计算模型.利用数值模拟手段,求解RNG (Renormalization Group)k-epsilon湍流模型下的N-S方程组,获得与有关单位试验相吻合的计算结果.分析发现对称面圆周角和伞衣幅高度对降落伞阻力影响很小.算例中阻力随对称面圆周角的变化在±0.28%以内.对称面圆周角的变化对轴对称降落伞尾流区流场的拓扑结构没有影响.对称轴上存在2个鞍点,随着对称面圆周角的增大,第1个鞍点的位置几乎不变,第2个鞍点的位置向尾流方向推移.      

高升阻比乘波体外形设计及气动特性计算研究

采用数值模拟方法研究在设计点(马赫数6,迎角0°)锥导乘波体气动外形的设计方法及其基本气动特性,以及在非设计点时该乘波体的气动特性,即各个气动系数随迎角和马赫数的变化特性.研究表明:基于无粘锥形流的乘波体气动外形的反设计方法是成功的;在设计点附近选取合适的半基准圆锥激波角并考虑粘性影响时,可得到乘波体最大升阻比为3.36;给出了采用这一布局的单级入轨运载器的可行的飞行控制方案;分析提出了进一步提高该乘波体气动布局升阻比的有效途径.      

航空发动机试车台附加阻力修正方法

附加阻力(进气冲量阻力和外部冲量阻力)是航空发动机试车台架推力的重要组成部分,准确地确定附加阻力对提高发动机试车台推力测量的准确性有十分重要的意义。本文对发动机总推力和试车台的对比标定进行了介绍,对发动机在露天试车台、室内地面台和高空台上附加阻力的确定方法和修正方法进行了详细地研究,给出了高空台与高空台、高空台与地面台的对比标定试验结果。      

改装双三角翼对歼击机升阻特性的影响

本文通过理论分析和风洞试验，研究改装双三角翼对飞机动性能的影响。所得结果对飞机设计，特别是改型设计有一定实用价值。     

较低雷诺数零升阻力估算方法研究

根据层流和湍流的摩擦阻力估算公式,建立较低雷诺数条件下当量摩擦阻力系数的计算模型来估算零升阻力系数。通过对估算方法进行验证和分析,研究较低雷诺数时零升阻力系数随雷诺数的变化趋势。结果表明,建立的估算方法合理,并且准确地反映了零升阻力系数在较大雷诺数范围内的变化规律。