实现应变率为105~106 s-1的阻抗梯度飞片复杂加载波形计算分析

为了在气炮上实现应变率为105~106 s-1的复杂加载技术研究,采用自行研制的拉格朗日程序MLEP(multi-material Lagrangian elastic-plastic)对Al-Cu-W材料体系的阻抗梯度飞片复杂加载不锈钢靶板进行数值模拟,计算设计并分析了阻抗梯度飞片的厚度和密度分布指数对靶板压力、速度和应变率峰值等波形的影响。结果表明:密度指数分布越大,加载时间越短,加载后期的压力、速度和应变率峰值曲线更陡峭;同时, 为了避免靶板/LiF窗口界面反射的稀疏波早于阻抗梯度飞片后界面反射的稀疏波达到碰撞面位置,计算设计中还考虑了飞片厚度的影响。此外,对基于理论设计的阻抗梯度飞片进行了动态考核实验,实验结果基本反映了预期的设计,为材料强度的测量奠定了基础。     

稳定分层流密度界面处湍流混合与分形结构

采用室内混合箱研究稳定分层流(上层淡水、下层盐水) 无剪切密度界面处的湍流混合与分形结构. 湍流通过浸没在盐水层中的振动格栅产生, 密度界面结构通过在盐水层中添加荧光剂或染料可视化, 共进行12 组实验. 实验观测并记录了:(1) 淡、盐水密度界面距混合箱底部的平均高程(h);(2) 淡、盐水层的密度(ρ₀,ρ), (3) 淡、盐水密度界面. 其中, 淡、盐水密度界面通过照片、录像进行记录. 观测结果用于计算:(1) 盐水层密度;(2) 卷挟速度, (3) 整体理查孙数(Ri_o), (4) 二维、三维密度界面, (5) 二维、三维密度界面的分形维度. 结果分析发现:(1) 湍流卷挟率随Rio 增大而减小, 并且满足Ri_o的-3=2 或-7=4 幂律;表明随着湍流强度的减弱, 混合的速度也越来越缓慢;(2) 二维密度界面分形维度大于1, 三维密度界面分形维度大于2;表明二维、三维密度界面存在分形结构;(3) 分形维度随Ri_o的增大而减小;表明随着湍流强度的减弱, 密度界面也越来越趋于光滑.      

Interfacial area concentration in gas–liquid bubbly to churn flow regimes in large diameter pipes

This study performed a survey on existing correlations for interfacial area concentration (IAC) prediction and collected an IAC experimental database of two-phase flows taken under various flow conditions in large diameter pipes. Although some of these existing correlations were developed by partly using the IAC databases taken in the low-void-fraction two-phase flows in large diameter pipes, no correlation can satisfactorily predict the IAC in the two-phase flows changing from bubbly, cap bubbly to churn flow in the collected database of large diameter pipes. So this study presented a systematic way to predict the IAC for the bubbly-to-churn flows in large diameter pipes by categorizing bubbles into two groups (group 1: spherical or distorted bubble, group 2: cap bubble). A correlation was developed to predict the group 1 void fraction by using the void fraction for all bubble. The group 1 bubble IAC and bubble diameter were modeled by using the key parameters such as group 1 void fraction and bubble Reynolds number based on the analysis of Hibiki and Ishii (2001, 2002) using one-dimensional bubble number density and interfacial area transport equations. The correlations of IAC and bubble diameter for group 2 cap bubbles were developed by taking into account the characteristics of the representative bubbles among the group 2 bubbles and the comparison between a newly-derived drift velocity correlation for large diameter pipes and the existing drift velocity correlation of Kataoka and Ishii (1987) for large diameter pipes. The predictions from the newly-developed two-group IAC correlation were compared with the collected experimental data in gas–liquid bubbly to churn flow regimes in large diameter pipes and their mean absolute relative deviations were obtained to be 28.1%, 54.4% and 29.6% for group 1, group 2 and all bubbles respectively.     

二维联合概率分布函数构造方法对结构串联系统可靠度影响分析

简要介绍了两种构造联合分布函数近似方法：基于Pearson相关系数的近似方法P和基于Spearman相关系数的近似方法S。推导了基于直接积分方法的串联系统失效概率计算公式,提出了两构件功能函数间负相关时串联结构系统失效概率上限值的计算公式。以理论联合概率分布函数是二维极值分布为例研究了两种近似方法在串联系统可靠度计算中的精度。结果表明,两种近似方法能够有效地计算串联结构系统可靠度,且精度很高,为不完备概率信息条件下串联结构系统可靠度分析提供了一条有效的途径。当两构件的功能函数正相关时,两种近似方法误差随串联系统失效概率的减小而增加,但近似方法与精确方法系统失效概率的差别最大也不会超过2倍;当两构件的功能函数负相关时,两种近似方法误差随系统失效概率的减小而减小,但两种近似方法的失效概率几乎与精确解一样。此外,两种近似方法误差不是随构件间相关性的增加而单调增加,而是存在一极大值。     

STOCHASTIC P-BIFURCATION OF TRI-STABLE VAN DER POL-DUFFING OSCILLATOR

分析了乘性和加性噪声作用下三稳态Van der Pol-Duffing振子的随机P分岔. 首先用随机平均法得到系统的随机微分方程,求得系统响应幅值的稳态概率密度函数. 然后应用分岔分析的奇异性理论,求得随机P分岔发生的临界参数条件,得到多种定性不同的稳态概率密度曲线. 讨论了2种激励噪声强度和系统阻尼对响应稳态概率密度曲线峰的个数、各峰值相对大小的影响. 通过Monte-Carlo数值模拟对理论计算结果进行了验证.该方法可用于其他系统的随机P分岔分析.     

Monte Carlo simulations of down-scattered neutron and knock-on deuteron spectra in deuterium-tritium capsule implosions

A Monte Carlo particle tracking (MCPT) code has been developed and coupled to hydrodynamic simulations to generate and track primary and scattered neutrons in addition to scattered charged particles by post-processing. This code uses one dimensional (1-D) deuterium-tritium implosion profiles, and the computed down-scattered neutron and knock-on deuteron spectra are analyzed for different areal densities. The mixing effects on the spectra of down-scattered neutron and knock-on deuteron are also investigated. The implementation of the numerical scheme is analyzed, and the particle splitting technique is adopted, which is proven to efficiently reduce computational effort.     

邻二甲苯分子结构和红外光谱的DFT研究

用 Onsager 模型和密度泛函B3LYP方法在6-31G^＊基组水平上对邻二甲苯分子的几何构型进行了全优化，计算了最优构型的红外光谱．将计算所得的邻二甲苯的红外频率与其Sadtler标准红外光谱进行了比较，确定了邻二甲苯的分子结构，并对红外光谱进行了归属．同时，对 Onsager 模型的适用性进行了讨论，进而得出了 Onsager 模型适于研究纯液态的邻二甲苯分子结构和红外光谱的结论．     

Zero density of L-functions related to Maass forms

Let f(z) be a Hecke-Maass cusp form for SL ₂(?), and let L(s, f) be the corresponding automorphic L-function associated to f. For sufficiently large T, let N(σ, T) be the number of zeros ρ = β +iγ of L(s, f) with |γ| ? T, β ? σ, the zeros being counted according to multiplicity. In this paper, we get that for 3/4 ? σ ? 1 ? ?, there exists a constant C = C(?) such that N(σ,T) ? T ^2(1?σ)/σ(logT) ^C , which improves the previous results.     

Sums of primes and quadratic linear recurrence sequences

Let u be a sequence of positive integers which grows essentially as a geometric progression. We give a criterion on u in terms of its distribution modulo d, d = 1, 2,..., under which the set of positive integers expressible by the sum of a prime number and an element of u has a positive lower density. This criterion is then checked for some second order linear recurrence sequences. It follows, for instance, that the set of positive integers of the form p ＋ [（2 ＋ √3）n], where p is a prime number and n is a positive integer, has a positive lower density. This generalizes a recent result of Enoch Lee. In passing, we show that the periods of linear recurrence sequences of order m modulo a prime number p cannot be ＂too small＂ for most prime numbers p.     

有限域上长为2mpn 的负循环码

设F_q 是奇数阶有限域. 本文主要借助X^2mpn+1 在F_q 上的不可约因式分解来确定有限域F_q上所有长为2^mpⁿ 的负循环码和自对偶的负循环码的生成多项式, 这里p 是q-1 的奇素因子, m 和n是正整数.