Simulation Techniques for Cosmological Simulations

Modern cosmological observations allow us to study in great detail the evolution and history of the large scale structure hierarchy. The fundamental problem of accurate constraints on the cosmological parameters, within a given cosmological model, requires precise modelling of the observed structure. In this paper we briefly review the current most effective techniques of large scale structure simulations, emphasising both their advantages and shortcomings. Starting with basics of the direct N-body simulations appropriate to modelling cold dark matter evolution, we then discuss the direct-sum technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and the tree algorithms. Simulations of baryonic matter in the Universe often use hydrodynamic codes based on both particle methods that discretise mass, and grid-based methods. We briefly describe Eulerian grid methods, and also some variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.     

Particle Acceleration Mechanisms

In this paper we review the possible mechanisms for production of non-thermal electrons which are responsible for the observed non-thermal radiation in clusters of galaxies. Our primary focus is on non-thermal Bremsstrahlung and inverse Compton scattering, that produce hard X-ray emission. We first give a brief review of acceleration mechanisms and point out that in most astrophysical situations, and in particular for the intracluster medium, shocks, turbulence and plasma waves play a crucial role. We also outline how the effects of the turbulence can be accounted for. Using a generic model for turbulence and acceleration, we then consider two scenarios for production of non-thermal radiation. The first is motivated by the possibility that hard X-ray emission is due to non-thermal Bremsstrahlung by nonrelativistic particles and attempts to produce non-thermal tails by accelerating the electrons from the background plasma with an initial Maxwellian distribution. For acceleration rates smaller than the Coulomb energy loss rate, the effect of energising the plasma is to primarily heat the plasma with little sign of a distinct non-thermal tail. Such tails are discernible only for acceleration rates comparable or larger than the Coulomb loss rate. However, these tails are accompanied by significant heating and they are present for a short time of <10⁶ years, which is also the time that the tail will be thermalised. A longer period of acceleration at such rates will result in a runaway situation with most particles being accelerated to very high energies. These more exact treatments confirm the difficulty with this model, first pointed out by Petrosian (Astrophys. J. 557:560, 2001). Such non-thermal tails, even if possible, can only explain the hard X-ray but not the radio emission which needs GeV or higher energy electrons. For these and for production of hard X-rays by the inverse Compton model, we need the second scenario where there is injection and subsequent acceleration of relativistic electrons. It is shown that a steady state situation, for example arising from secondary electrons produced from cosmic ray proton scattering by background protons, will most likely lead to flatter than required electron spectra or it requires a short escape time of the electrons from the cluster. An episodic injection of relativistic electrons, presumably from galaxies or AGN, and/or episodic generation of turbulence and shocks by mergers can result in an electron spectrum consistent with observations but for only a short period of less than one billion years.     

几种航空铝材动态力学性能实验

利用分离式Hopkinson压杆实验设备测定了4种航空铝合金材料的动态力学性能,根据材料的应力应变曲线特性,选择Cowper-Symonds模型来拟合4种材料的动态本构.在此基础上从实验结果中分析和比较了这4种铝合金在动态冲击下的破坏形式和动态力学特性.实验结果表明LF6R和LF21M具有良好的塑性,而LY12CZ和LC4CS具有一定的脆性;随着4种材料屈服应力的提高,材料切线模量也增加,二者呈线性关系;实验的4种航空铝材中LC4CS应变率敏感程度最小,LF6R敏感程度最大,LF21M,LY12CZ和LC4CS的应变率因子曲线形状相同.      

Sharp boundaries of solar wind plasma structures and their relationship to solar wind turbulence

Sharp (<10 min) and large (>20%) solar wind ion flux changes are common phenomena in turbulent solar wind plasma. These changes are the boundaries of small- and middle-scale solar wind plasma structures which can have a significant influence on Earth’s magnetosphere. These solar wind ion flux changes are typically accompanied by only a small change in the bulk solar wind velocity, hence, the flux changes are driven mainly by plasma density variations. We show that these events occur more frequently in high-density solar wind. A characteristic of solar wind turbulence, intermittency, is determined for time periods with and without these flux changes. The probability distribution functions (PDF) of solar wind ion flux variations for different time scales are calculated for each of these periods and compared. For large time scales, the PDFs are Gaussian for both data sets. For small time scales, the PDFs from both data set are more flat than Gaussian, but the degree of flatness is much larger for the data near the sharp flux change boundaries.     

Great solar bursts of October 19, 22 and 23, 1989

Itapetinga measurements at 48 GHz with the multibeam technique are used to determine the relative position of solar burst centroid of emission with high spatial accuracy and time resolution. For the Great Bursts of October 19,22, 1989, with a large production of relativistic particles, and October 23, it is suggested that, at 48 GHz, the bursts might have originated in more then one source in space and time. Additionally the October 19 and 22 Ground Level Events exhibited very unusual intensity-time profiles including double component structures for the onset phase. The Bern observatory spectral radio emission data show a strong spectral flattening typical for large source inhomogeneties. The interpretation for this is that large solar flares are a superposition of a few strong bursts (separated both in space and time) in the same flaring region.     

Cyclic variations in the solar lower atmosphere

The Ca K line has been measured regularly nearly every month since 1974 at Kitt Peak. It is well known that the K₁ component of the Ca K line is formed in the temperature minimum region (TMR) of the solar atmosphere. Our study of the data of CaII K profiles over two solar cycles indicates that both in full disc integrated spectra and in center disc spectra, the distance between the red K₁ and the blue K₁ of the profiles and its average intensity show periodic variations. But the variation for the full disc integrated spectra fluctuates in the same way as the sunspot number does, while that for the center disc spectra has a time delay with respect to sunspot number. Non-LTE computations yield a cyclic temperature variation of about 17 K of the TMR in the quiet-Sun atmosphere and a cyclic variation of about 15–20 km in the height position of the TMR.     

悬停状态下直升机旋翼水滴撞击特性研究

介绍了直升机在悬停状态下,求解旋翼水滴撞击特性的数值方法.该方法根据直升机悬停状态下的准定常特性,对旋翼运动模型进行简化,并在对绕桨叶运动的气流场计算的基础上,采用有限差分方法求解水滴运动方程,得到水滴运动轨迹.进而,确定了水滴对旋翼表面的撞击特性参数,为旋翼防冰系统设计提供条件.最后,对NACA0012翼型的水滴撞击特性进行分析.      

对流层顶变化对上对流层/下平流层臭氧分布的影响

上对流层和下平流层(UT／LS)，位于8-25km高度之间，是大气中一个很特殊的区域．大部分的臭氧分布在下平流层，在下平流层臭氧的含量发生一个很小的变化，就会对气候和地面的紫外辐射产生很大的影响．而作为气象参数的对流层顶，是充分混合、缺乏臭氧的上对流层和层结稳定、臭氧丰富的下平流层之间的边界或过渡层，其变化对臭氧总量和分布有直接和明显的影响．本文使用二维模式模拟研究对流层顶变化对臭氧在UT／LS分布的影响．模拟结果表明对流层顶的季节变化对UT/LS的臭氧分布有明显的影响，臭氧的局地变化可以超过10％在冬季北半球中纬度对流层顶高度升高1km时，模式结果表明对臭氧分布的影响比较显著，局地变化可超过6％，但是对臭氧总量的影响较小，变化不超过5DU，小于观测资料统计分析的结果。     

Models of P/Borrelly: Activity and dust mantle formation

Comet 19P/Borrelly was observed by Deep Space One spacecraft on September 22, 2001 (Soderblom et al., 2002).The DS1 images show a very dark and elongate nucleus with a complex topography; the IR spectra show a strong red-ward slope consistent with a very hot and dry surface (345K to 300K). During DS1 encounter the comet coma was dominated by a prominent jet but most of the comet was inactive, confirming the Earth-based observations that <10% of the surface is actively sublimating. We have developed a thermal evolution model of comet PBorrelly, using a numerical code that is able to solve the heat conduction and gas diffusion equations at the same time across an idealized spherical nucleus ( De Sanctis et al., 1999, 2000; Capria et al., 2000; Coradini et al., 1997a,b). The comet nucleus is composed by water, volatiles ices and dust in different proportions. The refractory component is made by grains that are embedded in the icy matrix. The code is able to account for the dust release, contributing to the dust flux, and the formation of dust mantles on the comet surface. The model was applied to a cometary nucleus with the estimated physical and dynamical characteristics of P/Borrelly in order to infer the status and activity level of a body on such an orbit during the DS1 observation. The comet gas flux, differentiation and thermal behavior were simulated and reproduced. The model results are in good agreement with the DS1 flyby results and the ground based observations, in terms of activity, dust coverage and temperatures of the surface.     

Measurements of bidirectional reflectance of ordinary chondrite for muses-C in-situ detection

MUSES-C, a Japanese sample return mission, is targeting a small near Earth asteroid, 1998SF36, which is considered an S-type asteroid and is similar in spectroscopy to LL class ordinary chondrite meteorite ([Binzel et al., 2001]). Although this mission will bring us detailed photometric data, that is, disk-resolved bidirectional reflectance data of the asteroid, there were few bidirectional reflectance data of ordinary chondrite meteorites. For the purpose of comparison with the data obtained by the in-situ observation, we have performed measurements of bidirectional reflectance of ordinary chondrite samples.

Here we summarize the result of our laboratory measurements of the bidirectional reflectance and compare them with the scattering property of 1998SF36. Although the geometric albedo of 1998SF36 is higher than the typical value of S-type asteroids, however, the laboratory data of ordinary chondrite are similar to or brighter than the model disk-resolved reflectance of 1998SF36 derived from disk-integrated ground-based data. We found in our laboratory data that there is a positive correlation between the surface roughness and the strength of the opposition effect. A stronger opposition effect in the reflectance of 1998SF36 than the laboratory data may therefore indicate that the surface of the asteroid has rougher structure than our laboratory samples.