Normalized ionization yield function for various nuclei obtained with full Monte Carlo simulations

Several recent results important for production of ion pairs in the Earth atmosphere by various primary cosmic ray nuclei are presented. The direct ionization by various primary cosmic ray nuclei is explicitly obtained. The longitudinal profile of atmospheric cascades is sensitive to the energy and mass (charge) of the primary particle. In this study different cosmic ray nuclei are considered as primaries, namely Helium, Oxygen and Iron nuclei. The cosmic ray induced ionization is obtained on the basis of CORSIKA 6.52 code simulations using FLUKA 2006 and QGSJET II hadronic interaction models. The energy of the primary particles is normalized to GeV per nucleon. In addition, the ionization yield function Y is normalized as ion pair production per nucleon. The obtained ionization yield functions Y for various primaries are compared. The presented results and their application are discussed.     

Streamwise-body-force-model for rapid simulation combining internal and external flow fields

A streamwise-body-force-model (SBFM) is developed and applied in the overall flow simulation for the distributed propulsion system, combining internal and external flow fields. In view of axial stage effects, fan or compressor effects could be simplified as body forces along the streamline. These body forces which are functions of local parameters could be added as source terms in Navier-Stokes equations to replace solid boundary conditions of blades and hubs. The val-idation of SBFM with uniform inlet and distortion inlet of compressors shows that pressure perfor-mance characteristics agree well with experimental data. A three-dimensional simulation of the integration configuration, via a blended wing body aircraft with a distributed propulsion system using the SBFM, has been completed. Lift coefficient and drag coefficient agree well with wind tun-nel test results. Results show that to reach the goal of rapid integrated simulation combining inter-nal and external flow fields, the computational fluid dynamics method based on SBFM is reasonable.     

Photolytic behaviour of methane at Lyman-α and 248 nm: Studies in the frame of a simulation program of Titan’s atmosphere (S.E.T.U.P.)

This paper first describes briefly some of the forefront global simulations of Titan’s atmosphere that have been carried out up to now. In these experiments, an initial gaseous mixture of N₂/CH₄ is submitted to a single energy source and the retrieved gas and/or solid phase(s) is/are analyzed by different techniques.     

Propagation and interaction of interplanetary transient disturbances. Numerical simulations

We study the heliocentric evolution of ICME-like disturbances and their associated transient forward shocks (TFSs) propagating in the interplanetary (IP) medium comparing the solutions of a hydrodynamic (HD) and magnetohydrodynamic (MHD) models using the ZEUS-3D code [Stone, J.M., Norman, M.L., 1992. Zeus-2d: a radiation magnetohydrodynamics code for astrophysical flows in two space dimensions. i – the hydrodynamic algorithms and tests. Astrophysical Journal Supplement Series 80, 753–790]. The simulations show that when a fast ICME and its associated IP shock propagate in the inner heliosphere they have an initial phase of about quasi-constant propagation speed (small deceleration) followed, after a critical distance (deflection point), by an exponential deceleration. By combining white light coronograph and interplanetary scintillation (IPS) measurements of ICMEs propagating within 1 AU [Manoharan, P.K., 2005. Evolution of coronal mass ejections in the inner heliosphere: a study using white-light and scintillation images. Solar Physics 235 (1–2), 345–368], such a critical distance and deceleration has already been inferred observationally. In addition, we also address the interaction between two ICME-like disturbances: a fast ICME 2 overtaking a previously launched slower ICME 1. After interaction, the leading ICME 1 accelerates and the tracking ICME 2 decelerates and both ICMEs tend to arrive at 1 AU having similar speeds. The 2-D HD and MHD models show similar qualitative results for the evolution and interaction of these disturbances in the IP medium.     

有剪切速度的Sweet-Parker薄电流片在非线性阶段的不稳定性研究

以Harris Sheet作为初始条件, 使用数值模拟的方法, 研究了二级磁岛不稳定重联的一些性质. 在模拟中随着初始扰动的加入, Harris Sheet将演化到非线性阶段, 形成更薄的有剪切速度的电流片, 并伴有一级磁岛产生. 当Lundquist数大于或等于10⁵时, 非均匀剪切速度的Sweet-Parker电流片开始不稳定, 并有二级磁岛出现. 不稳定Sweet-Parker电流片对应的临界长宽比为65. Lundquist数越大, 演化形成的Sweet-Parker电流片越薄, 更多的二级磁岛将出现, 且沿电流片两边向外喷出, 随时间变大, 相互合并. 二级磁岛的出现使重联率增大, 并与Lundquist数之间不再满足S^-1/2的关系, 而似乎对它的依赖关系不明显.     

A Class of TVD Type Combined Numerical Scheme for MHD Equations With a Survey About Numerical Methods in Solar Wind Simulations

It has been believed that three-dimensional, numerical, magnetohydrodynamic (MHD) modelling must play a crucial role in a seamless forecasting system. This system refers to space weather originating on the sun; propagation of disturbances through the solar wind and interplanetary magnetic field (IMF), and thence, transmission into the magnetosphere, ionosphere, and thermosphere. This role comes as no surprise to numerical modelers that participate in the numerical modelling of atmospheric environments as well as the meteorological conditions at Earth. Space scientists have paid great attention to operational numerical space weather prediction models. To this purpose practical progress has been made in the past years. Here first is reviewed the progress of the numerical methods in solar wind modelling. Then, based on our discussion, a new numerical scheme of total variation diminishing (TVD) type for magnetohydrodynamic equations in spherical coordinates is proposed by taking into account convergence, stability and resolution. This new MHD model is established by solving the fluid equations of MHD system with a modified Lax-Friedrichs scheme and the magnetic induction equations with MacCormack II scheme for the purpose of developing a combined scheme of quick convergence as well as of TVD property. To verify the validation of the scheme, the propagation of one-dimensional MHD fast and slow shock problem is discussed with the numerical results conforming to the existing results obtained by the piece-wise parabolic method (PPM). Finally, some conclusions are made. This revised version was published online in August 2006 with corrections to the Cover Date.     

Ionization effect of solar protons in the Earth atmosphere – Case study of the 20 January 2005 SEP event

The cosmic ray ground level enhancement on January 20, 2005 is among the largest recorded events in the history of cosmic ray measurements. The solar protons of MeV energies cause an excess of ionization in the atmosphere, specifically over polar caps following major solar disturbances. The ionization effect in the Earth atmosphere is obtained for various latitudes on the basis of solar proton energy spectra, reconstructed from GOES 11 measurements and subsequent full Monte Carlo simulation of cosmic ray induced atmospheric cascade. The estimation of ionization rates is based on a numerical model for cosmic ray induced ionization. The evolution of atmospheric cascade is performed with the CORSIKA 6.52 code using FLUKA 2006b and QGSJET II hadron interaction models. The atmospheric ion rate ionization is explicitly obtained for 40°N, 60°N and 80°N latitudes. The time evolution of obtained ion rates is presented. It is demonstrated that ionization effect is negative for 40°N and small for 60°N, because of accompanying Forbush decrease. The ionization effect is significant only in sub-polar and polar atmosphere during the major ground level enhancement of 20 January 2005.     

Closed-loop distance-keeping for long-term satellite cluster flight

Cluster flight is a term used for describing multiple satellites that are being held within pre-defined minimum and maximum distances for long time intervals, possibly the entire mission. This technology is required for a myriad of space architectures and missions, including disaggregated space architectures. Whereas the literature is abundant with works on control laws for satellite formation flying, there are only a handful of works on control of cluster flight. The purpose of the current work is to develop a cluster flight control algorithm, which is able to keep the satellites of the cluster within pre-specified minimum and maximum distances, while utilizing small amounts of propellant. The newly developed algorithm relies on the natural inter-satellite distance dynamics. The algorithm incorporates realistic mission constraints, such as constant-magnitude thrust, and is implemented in feedback form, steering the mean elements to judiciously selected reference values. Simulations indicate that a few tens of grams of propellent are sufficient for operating a cluster flight mission in excess of 1 year, using low specific-impulse thrusters.     

直升机主桨叶的鸟撞有限元数值模拟

以某型直升机主桨叶为例,使用瞬态动力学分析软件(MSC.Dytran),采用先进的任意拉格朗日-欧拉(ALE)流固耦合算法,进行鸟撞过程的数值模拟.模拟结果表明:撞击过程中的最大应力达到199MPa,没有超出桨叶材料的屈服极限.研究结果为桨叶的抗鸟撞损伤设计提供了参考.