横向互击型层板喷注器的喷雾特性

利用激光粒子动态分析仪 (PDA)对 180°互击型层板式喷注器的喷雾特性进行了试验研究 ,测量了喷雾典型区域的液滴速度、直径、体积通量密度等参数。分析了喷注器的雾化过程和喷雾特点 ,同时重点研究了喷口槽宽、推进剂流道面积比以及推进剂喷注速度比等参数对喷雾特性的影响。试验结果表明 :喷雾沿喷口展向可分为三个明显区域 ,中间区域流量较小但雾化较细 ,从槽两端喷出的液雾流量较大但液滴尺寸较大。对于一定的流道截面比γd 和喷注速度比γv,存在一个最佳的喷嘴槽宽w使雾化质量最优。同时 ,在喷嘴槽宽合适时 ,流道截面比γd 和喷注速度比γv 的适当增大有利于改善雾化     

磁层中重离子O＋通量密度分布的理论研究

应用求解得到的分布函数，研究了北半球子午面内极光带区起源的重离子Ｏ＋在磁层中的通量密度定态分布及其特性．结果表明：（１）重离子Ｏ＋通量密度分布沿ＧＳＭＺ轴方向呈峰状结构，其峰值随地心距离增大而减小，在近地空间减小较快；（２）磁扰动时，Ｏ＋离子占据的空间范围较磁宁静时大，通量峰较不突出；（３）理论估算的磁尾等离子体片边界层附近重离子Ｏ＋通量密度与观测结果相一致．     

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

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

基于边界涡量流的二维叶型优化建模与求解

针对二维叶片整体造型优化问题,提出一种基于边界涡量流的二维叶型优化模型,并对该模型进行求解。建立以含参数的Bezier曲线定义的二维叶型型线为边界、以边界涡量流正峰值的最小化为目标函数的优化模型,通过分析边界涡量流的生成机制,把叶型型线边界信息代入边界涡量流中,将边界涡量流用叶型型线参数和型线上的气动参数表示,进而获得新的目标函数表达式。该目标函数表达式可直接对叶型型线参数求导,仅使用简单的定步长梯度算法即可对优化模型求解。结果表明:优化后的叶型型线上的边界涡量流的正峰值显著降低。     

Unusually strong magnetic fields in Titan’s ionosphere: T42 case study

Observations of unusually large magnetic fields in the ionosphere indicate periods of maximum stress on Titan’s ionosphere and potentially of the strongest loss rates of ionospheric plasma. During Titan flyby T42, the observed magnetic field attained a maximum value of 37 nT between an altitude of 1200 and 1600 km, about 20 nT stronger than on any other Titan pass and close to five times greater in magnetic pressure. The strong fields occurred near the corotation-flow terminator rather than at the sub-flow point, suggesting that the flow which magnetized the ionosphere was from a direction far from corotation and possibly towards Saturn. Extrapolation of solar wind plasma conditions from Earth to Saturn using the University of Michigan MHD code predicts an enhanced solar wind dynamic pressure at Saturn close to this time. Cassini’s earlier exits from Saturn’s magnetosphere support this prediction because the Cassini Plasma Spectrometer instrument saw a magnetopause crossing three hours before the strong field observation. Thus it appears that Titan’s ionosphere was magnetized when the enhanced solar wind dynamic pressure compressed the Saturnian magnetosphere, and perhaps the magnetosheath magnetic field, against Titan. The solar wind pressure then decreased, leaving a strong fossil field in the ionosphere. When observed, this strong magnetic flux tube had begun to twist, further enhancing its strength.     

Interface flux reconstruction method based on optimized weight essentially non-oscillatory scheme

Aimed at the computational aeroacoustics multi-scale problem of complex configurations discretized with multi-size mesh, the flux reconstruction method based on modified Weight Essentially Non-Oscillatory (WENO) scheme is proposed at the interfaces of multi-block grids. With the idea of Dispersion-Relation-Preserving (DRP) scheme, different weight coefficients are obtained by optimization, so that it is in WENO schemes with various characteristics of dispersion and dissipation. On the basis, hybrid flux vector splitting method is utilized to intelligently judge the amplitude of the gap between grid interfaces. After the simulation and analysis of 1D convection equation with different initial conditions, modified WENO scheme is proved to be able to independently distinguish the gap amplitude and generate corresponding dissipation according to the grid resolution. Using the idea of flux reconstruction at grid interfaces, modified WENO scheme with increasing dissipation is applied at grid points, while DRP scheme with low dispersion and dissipation is applied at the inner part of grids. Moreover, Gauss impulse spread and periodic point sound source flow among three cylinders with multi-scale grids are carried out. The results show that the flux reconstruction method at grid interfaces is capable of dealing with Computational AeroAcoustics (CAA) multi-scale problems.     

On the Flux of Water and Minor Volatiles from the Surface of Comet Nuclei

Surface temperature and the available effective energy strongly influence the mass flux of H₂O and minor volatiles from the nucleus. We perform computer simulations to model the gas flux from volatile, icy components in porous ice-dust surfaces, in order to better understand results from observations of comets. Our model assumes a porous body containing dust, one major ice component (H₂O) and up to eight minor components of higher volatility (e.g. CO, CH₄, CH₃OH, HCN, C₂H₂, H₂S), The body's porous structure is modeled as a bundle of tubes with a given tortuosity and an initially constant pore diameter. Heat is conducted by the matrix and carried by the vapors. The model includes radially inward and outward flowing vapor within the body, escape of outward flowing gas from the body, complete depletion of less volatile ices in outer layers, and recondensation of vapor in deeper, cooler layers. From the calculations we obtain temperature profiles and changes in relative chemical abundances, porosity and pore size distribution as a function of depth, and the gas flux into the interior and into the atmosphere for each of the volatiles at various positions of the body in its orbit. In this paper we relate the observed relative molecular abundances in the coma of Comet C/1995 O1 (Hale-Bopp) and of Comet 46P/Wirtanen to molecular fluxes at the surface calculated from our model. This revised version was published online in June 2006 with corrections to the Cover Date.     

磁尾等离子体团型结构的数值模拟

考虑横越磁尾不同区域的数密度与离子温度的分布特点,取宁静磁尾温度、密度呈同样形态的非均匀分布,作为模拟计算的初态,对初始By为不同分布的3个算例作模拟计算．数值结果展示了3类磁结构的演化特征．考察磁尾中性片一个给定点的磁场变化,做出3个算例典型事件的磁场矢端图．它们从另一个侧面展示了磁尾通量绳型等离子体团、具有复杂闭合磁力线位形类似于“闭合环”的等离子体团,以及二维“磁岛”型等离子体团的特点．观测表明,多数磁尾等离子体团为具有强核心场的通量绳结构．通过跟踪一个通量绳型等离子体团的发展,作出等离子体团各方向速度与磁场强度随x演化的曲线．其中,vx与磁场强度|B|在近尾至中尾的计算值与Geotail卫星资料统计分析结果大致相符．此外,与Jin等取初始温度为均匀分布的计算结果相比,本文给出的通量绳型等离子体团,其内温度较高、密度较低,与资料分析结果的偏离也随之减小．     

Prediction of geosynchronous electron fluxes using an artificial neural network driven by solar wind parameters

There are hundreds of satellites operating at the geosynchronous (GEO) orbit where relativistic electrons can cause severe damage. Thus, predicting relativistic electron fluxes is significant for spacecraft safety. In this study, using GOES satellite data during 2011–2020, we propose two neural network models with two hidden layers to predict geosynchronous relativistic electron fluxes at two energy channels (>0.8 MeV and > 2 MeV). The number of input neurons of the two channels (>0.8 MeV and > 2 MeV) are determined to be 36 and 44, respectively. The > 0.8 MeV model has 22 and 9 neurons in the hidden layers, while the > 2 MeV model has 25 and 15 neurons in the hidden layers. The input parameters include the north–south component of the interplanetary magnetic field, solar wind speed, solar wind dynamic pressure and solar wind proton density. Through the analysis of different time delays, we determine that the optimal time delays of two energy channels (>0.8 MeV and > 2 MeV) are 8 days and 10 days, respectively. The training set and validation set (Jan 2011-Dec 2018) are divided by the 10-fold cross-validation method, and the remaining data (Jan 2019-Feb 2020) is used to analyze the model performance as a test set. The prediction results of both energy channels show good agreement with satellite observations indicated by low RMSE (～0.3 cm^-2sr^-1s^?1), high PE (～0.8) and CC (～0.9). These results suggest that only using solar wind parameters is capable of obtaining reasonable predictions of geosynchronous relativistic electron fluxes.     

Scientific visualization to study Flux Transfer Events at the Community Coordinated Modeling Center

In this paper we present recent additions to the visualization toolset offered by the Community Coordinated Modeling Center (CCMC). Two suites of visualization tools are available that can address different needs during the analysis of model simulations of the magnetosphere that are provided by the CCMC. The online, server-side visualization allows the user to quickly browse through simulation runs and now can create maps of magnetic field line topology in the magnetosphere. The second tool, SWX, can be used on the client computer after data have been downloaded. With this second tool the user can interact directly with the three-dimensional objects that are being rendered. We present results from a simulation of a Flux Transfer Event that was performed at the CCMC using a magnetohydrodynamic model of the Earth’s magnetosphere with a high resolution grid focused on the dayside magnetosheath and dayside magnetopause. The simulation shows that the FTE that results from localized magnetic reconnection is a complicated three-dimensional structure that requires modern visualization techniques. Visualization techniques that are presented here allow the researcher to fully appreciate the complexity contained in magnetospheric simulation results.