A Calculation Approach to Elastic Constants of Crystallines at High Pressure and Finite Temperature

Elastic constants of Na and Li metals are calculated successfully for temperatures up to 350K and pressures up to 30GPa using a scheme without involving any adjustable parameter.Elastic constants are assumed to depend only on an effective pair potential that is only determined by the average interatomic distance.Temperature has an effect on elastic constants by way of charging the equilibrium.The elastic constants can be obtained by fitting the relationship between total energy and strain tensor using the new set of lattice parameters obtained by calculating displacement of atoms at the finite temperature and at a fixed pressure.The relationship between the effective pair potential and the interatiomic distance is fitted by using a series of data of cohesive energy corresponding to lattice parameters.     

Spectrum Simulation of Li-Like Aluminium Plasma

X-ray emission spectra for L-shell of Li-like aluminium ions are simulated by using the flexible atomic code based on the collisional radiative model. Atomic processes including radiative recombination, dielectronic recombination, collisional ionization and resonance excitation from the neighbouring ion （Al^9＋ and Al^11＋ ） charge states of the target ion （Al^10＋） are considered in the model. In addition, the contributions of different atomic processes to the x-ray spectrum are analysed. The results show that dielectronic recombination, radiative recombination, collisional ionization and resonance excitation, other than direct collisional excitation, are very important processes.     

基于二次插值的非饱和土壤水流问题的特征差分方法及数值模拟

本文针对一类非饱和土壤水流问题,提出了基于二次插值的特征差分格式,得到了严谨的L2模误差估计.并作了数值试验,指明方法的有效性.     

Growth of NaBi（WO4）2 Dendrite and Mechanism

The solid-liquid interface motion of NaBi（WO4）2 （NBWO） melt crystal growth is observed in an in situ system, in which the whole processes of interface transition from fiat interface and cellular to dendrite are visualized. The spacing of the dendrite under smaller temperature gradient turns out to be larger than that under larger temperature gradient, which is found to be sensitive to the temperature distribution. The mechanism of dendrite growth of NBWO is studied based on the model of the growth units of anion coordination polyhedra. The { 001} face has two apex links, so it shows higher stability and has high growth rate and forms the arm of dendrite, whereas the {010} face has only one apex link, and thus shows relative slower growth rate and firstly forms the branches.     

Frequency Selectivity Behaviour in the Auditory Midbrain： Implications of Model Study

By numerical simulations on frequency dependence of the spiking threshold, i.e. on the critical amplitude of periodic stimulus, for a neuron to fire, we find that bushy cells in the cochlear nuclear exhibit frequency selec- tivity behaviour. However, the selective frequency band of a bushy cell is far away from that of the preferred spectral range in human and mammal auditory perception. The mechanism underlying this neural activity is also discussed. Further studies show that the ion channel densities have little impact on the selective frequency band of bushy cells. These findings suggest that the neuronal behaviour of frequency selectivity in bushy cells at both the single cell and population levels may be not functionally relevant to frequency discrimination. Our results may reveal a neural hint to the reconsideration on the busily cell functional role in auditory information processing of sound frequency.     

上海光源储存环束流托歇克寿命研究

上海同步辐射装置(SSRF)是目前在建的第3代专用同步辐射光源. 其储存环电子能量3.5GeV, 设计束团发射度3.9nm•rad. 托歇克寿命将是影响束流寿命的最主要因素, 它主要受限于储存环的能量接受度. 储存环的能量接受度不但取决于 高频电压, 同时也受动力学孔径和物理孔径的影响. 因此能量接受度的计算将是复杂的. 通过计算机模拟跟踪的方法计算储存环各点的能量接受度, 其程序是建立在Accelerator Toolbox(AT)基础上的自编程序. 通过这些能量接受度数据给出更加准确的托歇克寿命, 并且分析了感兴趣的不同运行条件下的托歇克寿命的变化情况.     

HL-2A托卡马克边界等离子体数值模拟及非日冕平衡条件下的杂质辐射问题研究

对边界层等离子体中常见的物理问题，从两点模型到二维流体描述。从原子分子物理过程到杂质的输运和辐射等进行了系统的归纳和总结，特别是对等离子体不同参数运行区如鞘层限制参数区、传导限制参数区以及脱靶参数区等的一维流体描述，在参阅相关文献的基础上使用一定的假设条件进行了简单推导。分别阐述了它们的特点。     

稀薄气体二维外部柱体绕流特性研究

本文采用直接模拟蒙特卡罗方法对稀薄气体二维外部柱体绕流问题进行了数值模拟。结果表明:外部绕流问题,在特定情况下会产生激波,激波的产生,不仅与气体的稀薄程度有关,还与来流马赫数有关。而气流与壁面之间的换热,随来流马赫数增加而增加,随气体稀薄程度增加而减小。     

混凝土动态劈裂拉伸试验的数值模拟

采用LS-DYNA有限元编码程序，对混凝土的动态劈裂拉伸试验进行了数值模拟，给出了混凝土试样在不同类型加载情况下应力分布之间的差异，证实了动态劈裂拉伸试验的有效性，提出了改进试验方法的若干结论．     

二维欧拉多介质流体力学数值模拟-TVD＋VOF方法

用TVD格式结合VOF界面处理方法编制了二维多介质高分辨欧拉程序，以解决冲击波和多介质界面处理。程序包括单介质网格高精度流体力学计算、多介质网格内界面重构、各种介质输运和压力驰豫平衡过程。其中单介质网格的计算采用Harten二阶TVD格式结合MacCormark方法计算含有源项的非齐次守恒定律方程组，通过4节点限制函数保证格式单调。多介质网格采用Youngs方法构造界面，采用x，y方向分裂格式计算体积份额输运，再根据体积份额输运计算质量、动量和能量的输运，最后利用等熵条件计算各种介质的压力驰豫平衡过程。