Universal scaling laws for fully-developed magnetic field turbulence near and far upstream of the Earth’s bow shock

We analyze the multifractal scaling of the modulus of the interplanetary magnetic field near and far upstream of the Earth’s bow shock, measured by Cluster and ACE, respectively, from 1 to 3 February 2002. The maximum order of the structure function is carefully estimated for each time series using two different techniques, to ensure the validity of our high-order statistics. The first technique consists of plotting the integrand of the pth order structure function, and the second technique is a quantitative method which relies on the power-law scaling of the extreme events. We compare the scaling exponents computed from the structure functions of magnetic field differences with the predictions obtained by the She–Lévêque model of intermittency in anisotropic magnetohydrodynamic turbulence. Our results show a good agreement between the model and the observations near and far upstream of the Earth’s bow shock, rendering support for the modelling of universal scaling laws based on the Kolmogorov phenomenology in the presence of sheet-like dissipative structures.     

Numerical modeling of propagation of breaking nonlinear acoustic-gravity waves from the lower to the upper atmosphere

Acoustic-gravity waves (AGWs) observed in the upper atmosphere may be generated near the Earth’s surface due to a variety of meteorological sources. Two-dimensional simulations of vertical propagation and breaking of nonlinear AGWs in the atmosphere are performed. Forcing near the Earth’s surface is used as the AGW source in the model. We use a numerical method based on finite-difference analogues of fundamental conservation laws for solving atmospheric hydrodynamic equations. This approach selects physically correct generalized solutions of the wave hydrodynamic equations. Numerical simulations are performed in a representative region of the Earth’s atmosphere up to altitude 500 km. Vertical profiles of temperature, density, molecular viscosity and heat conductivity were taken from the standard atmosphere model MSIS-90 for January. Calculations were made for different amplitudes and frequencies of lower boundary wave forcing. It is shown that after activating the tropospheric wave forcing, the initial pulse of AGWs may very quickly propagate to altitudes of 100 km and above and relatively slowly dissipate due to molecular viscosity and heat conduction. This may increase the role of transient nonstationary waves in effective energy transport and variations of atmospheric parameters and gas admixtures in a broad altitude range.     

锥柱形药柱环形槽内传热现象的数值研究

采用标准k-ε湍流模型和离散传递热辐射模型对不同来流雷诺数下锥柱形药柱前向环形槽内的流场与传热进行了数值研究。对槽内各表面与热气流之间的换热量进行了比较,并对努塞尔数沿环形槽各表面的变化进行了分析。结果表明来流雷诺数的大小影响环形槽内的旋涡数和旋涡区域的分布,进而影响槽内温度的分布,最终影响槽内装药表面与点火热气流之间的换热量。整个计算结果可为装药设计和分析推进剂的热解燃烧提供参考。     

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

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

用基于混合网格的多重网格法数值模拟绕物体的紊流流动

本文在混合网格基础上用多重网格法求解了紊流N-S方程。混合网格是由在物面附近采用的柱状网格与其他区域则采用的非结构网格组成。方程的求解采用Jamson的有限体积法，紊流模型采用两层Baldwin-Lomax代数紊流模型。用多重网格法来加速解的收敛。数值算例表明，用混合网格及多重网格法来求解紊流Navier-Stokes是非常有效的。     

On the effects of dynamical turbulence on the perpendicular diffusion of low-energy cosmic ray electrons

Diffusion perpendicular to the heliospheric magnetic field plays an integral role in the transport of charged particles in the heliosphere. In this study the perpendicular diffusion coefficient of low-energy cosmic ray electrons is calculated, using an equation derived from the random ballistic decorrelation interpretation of nonlinear guiding centre theory. An observationally motivated 2D turbulence power spectrum is assumed and the effects of various turbulence inputs on the resulting perpendicular diffusion coefficient are investigated. The perpendicular diffusion coefficients are first determined at 1 AU, for both magnetostatic and dynamical turbulence conditions. These solutions are also evaluated for radial distances of 0.1 AU to 10 AU to further investigate the values of the perpendicular diffusion coefficients in the very inner heliosphere. The results of this study show that the dissipation range of the turbulence power spectrum provides a negligible contribution towards the perpendicular diffusion coefficient, and that solutions derived using only the energy containing range serve as good approximations for solutions derived assuming the full 2D turbulence power spectrum. Finally, it is shown that the effects of dynamical turbulence, as considered in the present study, do not affect the perpendicular diffusion coefficients derived from the scattering theory considered here.     

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.     

Relation between different theories for cosmic ray cross field diffusion

In the past few year several theories have been presented for describing cosmic ray scattering across the mean magnetic field. It is the purpose of the present article to discuss the relation between these different theories in order to improve our understanding of cosmic ray perpendicular scattering and to explore the parameter regimes for which these different theories are valid.     

Random walk of magnetic field lines: Subdiffusive,diffusive, and superdiffusive regimes

One- and two-dimensional models of magnetic field fluctuations and turbulence are widely used in space-, astrophysical, and laboratory contexts. In the present article we use a generalized form of the turbulence wave spectrum to calculate field line diffusion coefficients analytically and numerically. General conditions are derived for which field line wandering behaves subdiffusively, diffusively, and superdiffusively.     

Study on dependence of hydrodynamic characteristics on gravity in a vertical upward channel bubbly flow

Since the bubbly flow has extensive applications in the space field, it is very necessary to comprehend the dependence of hydrodynamic characteristics on gravity. In this paper, the dependence of the microbubble distribution and the liquid turbulence modulation by the microbubbles on gravity was investigated in detail with the Euler–Lagrange method. The liquid Navier–Stokes equation was solved using direct numerical simulations (DNS), and the microbubble motion was tracked with Newtonian motion equation considering drag force, shear lift force, added mass force, pressure gradient force, and wall lift force. The coupling between the gas and the liquid phases regarded the interphase force as a momentum source term in the momentum equation of the liquid. The results showed that the phase profile and the turbulence modulation by the microbubbles strongly depend on the magnitude of gravity. When the influence of gravity is relatively weak, the microbubbles approximately uniformly disperse in the wide central region of the channel, and the average statistics of the liquid turbulence are almost not influenced due to the addition of the microbubbles. However, when the influence of gravity is comparatively important, the majority of the microbubbles accumulate near the wall of the channel, and the injection of the microbubbles modifies the profiles of the liquid average statistics.