欧洲计算流体力学——个人观点

介绍了70年代初到20世纪末欧洲科学家对CFD发展的主要贡献,并重点关注了对国际流体力学数值方法会议(ICNMFD)和国际计算力学学术会议(ISCFD)的成功举办具有重大贡献的科学家.ICNMFD始于1969年(新西伯利亚),由参与了美国和苏联太空竞赛的两国科学家联合创办,而ISCFD则作为ICNMFD在亚太地区的同类会议初创于1995年(东京).由Oshima教授和本人提议,得益于上述两个会议全体科学委员会成员的共同努力,ICNMFD和ISCFD自2000年开始,合并成为每两年举办一次的国际计算流体力学会议(ICCFD).绝大部分本文讨论所涉及的学者通常都在受邀进行大会报告后,成为了前述三大国际会议中的一个或多个科学委员会的成员.作者对文中讨论所涉及的研究人员的圈定负有全权责任,上述人员的选择更多地基于应用而非理论计算流体力学,因此尚存在不完善之处.此外由于文中所讨论的学者大都是作者在40多年计算流体力学研究生涯中的私交好友,因此并没有将年轻一代的科学家考虑进来,但这并非由于他们的贡献不重要.另外值得注意的是欧洲的计算流体力学同美国的计算流体力学有非常紧密的联系,本文讨论所涉及的研究学者有些最早从欧洲开始涉足计算流体力学而现在则生活和工作于美国(例如Glowinski,Roe,Van Leer等人),而有些则是在美国取得了博士学位(例如Launder,Napolitano),他们都同美国的科学家进行过合作,因此很难清晰地对欧洲和美国的计算流体力学进行界定.对那些现已移居美国的欧洲科学家们,本文优先考虑了他们移居之前的贡献.最后但并非是最不重要的,作者通过电子邮件一一联系了文中所涉及的科学家,希望他们本人能够给出其自身贡献的评价草稿,从而帮助作者以最权威的视角在文中相应章节对其工作进行评述.本文向提供自述评价的诸位致以谢意,同时也希望其他同仁能够对本人不恰当的擅自评述给予谅解.     

可重复使用运载器机翼外形优化

研究了有翼可重复使用运载器机翼平面外形参数优化设计问题。在高超声速配平和低速水平着陆性能等约束条件下，利用混合遗传算法和粒子群优化算法来选择最优机翼外形参数，使机翼重量最小化。在已知的不同质心位置下求得相应的最小重量机翼外形。亚声速气动力采用面源法计算，高超声速气动力采用修正的牛顿流理论计算。为了降低计算耗费，使用均匀设计和逐步回归方法建立了低速气动力回归模型。整个优化过程在Matlab集成环境下完成。     

A direct position determination method with combined TDOA and FDOA based on particle filter

The localization of a stationary transmitter using moving receivers is considered. The original Direct Position Determination (DPD) methods, with combined Time Difference of Arrival (TDOA) and Frequency Difference of Arrival (FDOA), do not perform well under low Signal-to-Noise Ratio (SNR), and worse still, the computation cost is difficult to accept when the computational capabilities are limited. To get better positioning performance, we present a new DPD algorithm that proves to be more computationally efficient and more precise for weak signals than the conventional approach. The algorithm partitions the signal received with the same receiver into multiple non-overlapping short-time signal segments, and then uses the TDOA, the FDOA and the coherency among the short-time signals to locate the target. The fast maximum likelihood estimation, one iterative method based on particle filter, is designed to solve the problem of high computation load. A secondary but important result is a derivation of closed-form expressions of the Cramer-Rao Lower Bound (CRLB). The simulation results show that the algorithm proposed in this paper outperforms the traditional DPD algorithms with more accurate results and higher computational efficiency, and especially at low SNR, it is more close to the CRLB.     

Saw-tooth substorms: Inconsistency of repetitive bay-like magnetic disturbances with behavior of aurora

The relationships between the magnetic disturbance onsets, aurora dynamics and particles injections at the geostationary orbit have been analyzed in detail for 25 sawtooth substorms. It is shown that inconsistency between the above signatures of the substorms onset is typical of the powerful sawtooth substorms, unlike the isolated (“classical”) magnetospheric substorms. The distinguishing feature of the aurora in case of saw-tooth substorms is permanently high level of auroral activity irrespective of the magnetic disturbance onsets and the double oval structure of the aurora display. The close relationship between the aurora behavior and the particle injections at geostationary orbit is also broken. The conclusion is made, that the classical concept of the substorm development, put forward by Akasofu (1964) for isolated substorms, is not workable in cases of the sawtooth disturbances, when the powerful solar wind energy pumping into the magnetosphere provides a permanent powerful aurora particle precipitation into the auroral zone.     

Diffusive shock acceleration at interplanetary perpendicular shock waves: Influence of the large scale structure of turbulence on the maximum particle energy

We calculate the maximum energy that a particle can obtain at perpendicular interplanetary shock waves by the mechanism of diffusive shock acceleration. The influence of the energy range spectral index of the two-dimensional modes of the interplanetary turbulence is explored. We show that changes in this parameter lead to energies that differ in at least one order of magnitude. Therefore, the large scale structure of the turbulence is a key input if the maximum particle energy is calculated.     

Study of the 28 October 2003 and 20 January 2005 solar flares by means of 2.223 MeV gamma-emission line

By the data on intensity-time profiles of the neutron capture line of 2.223 MeV we have studied some characteristics of two solar flares, 28 October 2003 and 20 January 2005 (INTEGRAL and CORONAS-F observations, respectively). The SINP code was applied making allowance for the main processes of neutron interactions and deceleration in the solar plasma, character of neutron source, losses of neutrons and density model of the solar atmosphere. Comparison of the computed time profiles of 2.223 MeV line with observed ones for the flare of 28 October 2003 confirms the results obtained earlier for three other flares. Namely, the effect of density enhancement (EDE) in the sub-flare region, as well as the variations (hardening) of accelerated particle spectrum in the course of the event have been confirmed. The usual modeling procedure by the SINP code, however, seems to be inapplicable to the event of 20 January 2005. Possible causes of density enhancements during some flares and peculiarities of the 20 January 2005 flare are discussed.     

Cosmic Ray research in Armenia

Cosmic Ray research on Mt. Aragats began in 1934 with the measurements of East–West anisotropy by the group from Leningrad Physics-Technical Institute and Norair Kocharian from Yerevan State University. Stimulated by the results of their experiments in 1942 Artem and Abraham Alikhanyan brothers organized a scientific expedition to Aragats. Since that time physicists were studying Cosmic Ray fluxes on Mt. Aragats with various particle detectors: mass spectrometers, calorimeters, transition radiation detectors, and huge particle detector arrays detecting protons and nuclei accelerated in most violent explosions in Galaxy. Latest activities at Mt. Aragats include Space Weather research with networks of particle detectors located in Armenia and abroad, and detectors of Space Education center in Yerevan.     

Simultaneous observations of Pc 1 micropulsation activity and stratospheric electrodynamic perturbations on 27 January 2003

The 2nd Polar Patrol Balloon campaign (2nd-PPB) was carried out at Syowa Station in Antarctica during 2002–2003. Identical stratospheric balloon payloads were launched as close together in time as allowed by weather conditions to constitute a cluster of balloons during their flights. A very pronounced negative ion conductivity enhancement was observed at 32 km in the stratosphere below the auroral zone on 27 January 2003 from 1500 to 2200 UT. During this event, the conductivity doubled for an interval of about 7 h. This perturbation was associated with an extensive Pc 1 or Pi 1 wave event that was observed by several Antarctic ground stations, balloon PPB 10, and the Polar spacecraft. No appreciable X-ray precipitation was observed in association with this event, which would point to >60 Mev proton precipitation as a possible magnetosphere–stratosphere coupling mechanism responsible for the conductivity enhancement. Such precipitation is consistent with the wave data. During the latter half of the event, E_z was briefly positive. There was a tropospheric Southern Ocean storm system underneath the balloon during this interval. If the event was associated with this storm system and not energetic proton precipitation, the observations imply an electrified Southern Ocean storm and major perturbations in stratospheric conductivity driven by a tropospheric disturbance. This event represents a poorly understood source for global circuit current. Precipitating energetic proton data from Akebono and NOAA POES spacecraft show significant >16 MeV precipitation was occurring at the location of PPB 8 but not PPB 10, suggesting that proton precipitation was, in fact, the responsible coupling mechanism.     

The high-latitude ionosphere structure on 22 March, 1979 magnetic storm from multi-satellite and ground-based observation

Variations in the high-latitude ionosphere structure during March 22, 1979 geomagnetic storm are examined. Electron density Ne and temperature Te from the Cosmos-900 satellite, NmF2, Ne and He⁺ from the ISS-b satellite, precipitation of soft electrons from the Intercosmos-19 satellite, and the global picture of the auroral electron precipitation from the DMSP, TIROS and P78 satellites are used. These multi-satellite databases allow us to investigate the storm-time variations in the locations of the following ionospheric structures: the day-time cusp, the equatorial boundary of the diffuse auroral precipitation (DPB), the main ionospheric trough (MIT), the day-time trough, the ring ionospheric trough (RIT) and the light ions trough (LIT). The variations in NmF2, Ne, He⁺ and Te in the high-latitude ionosphere for the different local time sectors are analyzed also. The features of the high-latitude ionospheric response to a strong magnetic storm are described.