一种简单的局部离散速度统一气体动理论格式

统一气体动理论格式UGKS(Unified Gas-Kinetic Scheme)是一种适用于从连续流到自由分子流的全流域计算格式。在该格式中一般使用统一的离散速度空间。而在高速流动中,不同节点的分布函数往往差异很大。为了保证计算的精度,离散速度空间必须满足所有节点的需要,占用了大量的内存。采用局部的均匀离散速度空间,离散速度的范围随节点状态的变化而变化,从而降低了内存的需要,并通过引入背景网格避免了不同节点离散速度的插值。最后,通过两个一维算例对该方法进行了测试。测试结果显示,采用局部离散速度空间能够得到可靠的结果,并且在模拟高速流动时计算效率明显提高。     

水槽涡的变化规律数值模拟

建立了不可压缩Navier-Stokes方程的Crank-Nicolson有限差分方法,数值模拟水槽晃动中流场及其涡流的数值变化规律。将数值解与解析解和前人的数值解进行比较,数值验证了不可压缩Naver-Stokes方程有限差分方法的有效性。通过数值模拟得到水槽在不同程度的倾斜激励晃动下流场及涡流的数值变换规律,当倾斜激励晃动的频率接近或远离共振频率时,水槽涡场的变化逐步由双涡变成单涡,再到不规则的涡场。当倾斜激励晃动的频率靠近共振频率ω_p=0.95ω_1附近时,水槽流场上部形成一个小涡,然后小涡扩大成整个水槽中的大涡,大涡下沉分裂成两个单涡,最后在底部消失;当倾斜激励晃动的频率在ω_p=0.75ω_1附近时,水槽底部形成一个小涡,然后扩大成大的单涡,最后在自由面消失;当倾斜激励晃动的频率在ω_p=0.55ω_1附近时,水槽底部出现小涡,然后扩大成大的单涡,大涡在自由面消失,继而出现不规则的大涡和不规则的小涡。     

Two-color laser wavelength effect on intense terahertz generation in air

Near-IR femtosecond lasers have been proposed to produce high-field terahertz radiation in the air via the laser-plasma interaction, but the physical mechanism still needs to be further explored. In this work, we theoretically investigate the effect of the two-color laser wavelength on the terahertz generation in the air based on a transient photocurrent model.We show that the long wavelength laser excitation can greatly enhance the terahertz amplitude for a given total laser intensity. Furthermore, we utilize a local current model to illustrate the enhancement mechanism. Our analysis shows that the terahertz amplitude is determined by the superposition of contributions from individual ionization events, and for the long wavelength laser excitation, the electron production concentrates in a few ionization events and acquires the larger drift velocities, which results in the stronger terahertz radiation generation. These results will be very helpful for understanding the terahertz generation process and optimizing the terahertz output.     

α decay properties of ~(297)Og within the two-potential approach

The α decay half-life of the unknown nucleus ~(297)Og is predicted within the two-potential approach, andα preformation probabilities of 64 odd-A nuclei in the region of proton numbers 82 Z 126 and neutron numbers 152 N 184, from ~(251)Cf to ~(295)Og, are extracted. In addition, based on the latest experimental data, a new set of parameters for α preformation probabilities considering the shell effect and proton-neutron interaction are obtained.The predicted α decay half-life of ~(297)Og is 0.16 ms within a factor of 4.97. The predicted spin and parity of the ground states for ~(269)Sg,~(285)Fl and ~(293)Lv are 3/2~+, 3/2~+ and 5/2~+, respectively.     

Physics potential of searching for 0vββ decays in JUNO

In the past few decades, numerous searches have been made for the neutrinoless double-beta decay (0vββ) process, aiming to establish whether neutrinos are their own antiparticles (Majorana neutrinos), but no 0vββ decay signal has yet been observed. A number of new experiments are proposed but they ultimately suffer from a common problem: the sensitivity may not increase indefinitely with the target mass. We have performed a detailed analysis of the physics potential by using the Jiangmen Underground Neutrino Observatory (JUNO) to improve the sensitivity to 0vββ up to a few meV, a major step forward with respect to the experiments currently being planned. JUNO is a 20 kton low-background liquid scintillator (LS) detector with 3%E(MeV) energy resolution, now under construction. It is feasible to build a balloon filled with enriched xenon gas (with ¹³⁶Xe up to 80%) dissolved in LS, inserted into the central region of the JUNO LS. The energy resolution is ～1.9% at the Q-value of ¹³⁶Xe 0vββ decay. Ultra-low background is the key for 0vββ decay searches. Detailed studies of background rates from intrinsic 2vββ and ⁸B solar neutrinos, natural radioactivity, and cosmogenic radionuclides (including light isotopes and ¹³⁷Xe) were performed and several muon veto schemes were developed. We find that JUNO has the potential to reach a sensitivity (at 90% C. L.) to T_1/2^0vββ of 1.8×10²⁸ yr (5.6×10²⁷ yr) with ～50 tons (5 tons) of fiducial ¹³⁶Xe and 5 years exposure, while in the 50-ton case the corresponding sensitivity to the effective neutrino mass, m_ββ, could reach (5-12) meV, covering completely the allowed region of inverted neutrino mass ordering.     

改进差分进化算法求解装载率凹费用装箱问题

研究了广泛存在于物流作业中一类新型的装箱问题,主要特征体现在箱子使用费用是关于装载率的凹函数。为求解问题,提出了一种基于分组编码策略的改进差分进化算法,以避免常规实数和整数编码方法存在放大搜索空间的不足。针对分组编码策略,定制化设计了以促进优秀基因传播为导向的新型变异和交叉操作,另外还嵌入了以物品置换为邻域的自适应局部搜索操作以增强局部搜索能力。对以往文献给出算例在不同凹费用函数下进行测试,实验结果显示所提出的算法明显优于BFD启发式算法,并且较遗传算法也有显著性改进。     

晶硅掺镓抑制太阳电池光致衰减效应的研究

采用掺镓晶硅和掺硼晶硅制备的寿命片和太阳电池片分别进行光致衰减实验,用WT-2000少子寿命测试仪和Halm电池电学性能测试仪等研究了它们受光照前后少子寿命和电学性能的变化.发现掺镓单晶寿命片的少子寿命衰减率比掺硼单晶寿命片低50;左右,掺镓单晶PERC电池和掺镓多晶常规电池转换效率的衰减率比掺硼单晶PERC电池和掺硼多晶常规电池分别降低3.41;和0.92;.这些结果表明晶硅太阳电池的光致衰减效应主要是晶硅中少子寿命降低导致的,晶硅掺镓后能有效抑制太阳电池的光致衰减现象.     

爆炸与冲击问题的大规模高精度计算

爆炸与冲击问题的数值模拟在国防和民用安全领域具有重要的工程实用价值.由于爆炸与冲击问题是一个多物质在高应变率、高温及高压条件下的强非线性的瞬态动力学问题,给数值模拟带来了很多的困难,为此,针对爆炸与冲击问题数值模拟中的一些关键和难点问题开展了研究.提出了三维非线性双曲守恒系统的伪弧长自适应网格算法,分析了算法的实现过程,数值结果表明该算法有效地提高了冲击波强间断处的分辨率.发展了针对气相爆轰数值模拟的附加龙格-库塔方法,对非线性对流项进行显示计算,化学反应源项进行半隐式计算,有效地解决了源项引起的刚性问题,计算结果表明该算法可以准确地捕捉和描述爆轰波的复杂结构和典型特征.针对三维工程实际物理问题中的大规模计算需求,给出了三维多物质流体动力学欧拉数值方法的并行化方法,开发了三维爆炸与冲击问题并行计算程序,并给出了针对该并行程序的测试方法.上述工作有利地解决了爆炸与冲击问题大规模、高精度计算中的一些难题.最后,开展了大口径聚能射流侵彻混凝土靶问题的数值模拟和实验研究,通过典型爆炸与冲击工程问题的计算验证了所研究数值方法的有效性.      

All-speed Roe scheme for the large eddy simulation of homogeneous decaying turbulence

As a type of shock-capturing scheme, the traditional Roe scheme fails in large eddy simulation (LES) because it cannot reproduce important turbulent characteristics, such as the famous k^?5/3 spectral law, as a consequence of the large numerical dissipation. In this work, the Roe scheme is divided into five parts, namely, ξ, δU_p, δp_p, δU_u, and δp_u, which denote basic upwind dissipation, pressure difference-driven modification of interface fluxes, pressure difference-driven modification of pressure, velocity difference-driven modification of interface fluxes, and velocity difference-driven modification of pressure, respectively. Then, the role of each part in the LES of homogeneous decaying turbulence with a low Mach number is investigated. Results show that the parts δU_u, δp_p, and δU_p have little effect on LES. Such minimal effect is integral to computational stability, especially for δU_p. The large numerical dissipation is due to ξ and δp_u, each of which features a larger dissipation than the sub-grid scale model. On the basis of these conditions, an improved all-speed Roe scheme for LES is proposed. This scheme can provide satisfactory LES results even for coarse grid resolutions with usually adopted second-order reconstructions for the finite volume method.     

Skin formation in drying a film of soft matter solutions:Application of solute based Lagrangian scheme

When a film of soft matter solutions is being dried, a skin layer often forms at its surface, which is a gel-like elastic phase made of concentrated soft matter solutions. We study the dynamics of this process by using the solute based Lagrangian scheme which was proposed by us recently. In this scheme, the process of the gelation(i.e., the change from sol to gel) can be naturally incorporated in the diffusion equation. Effects of the elasticity of the skin phase, the evaporation rate of the solvents, and the initial concentration of the solutions are discussed. Moreover, the condition for the skin formation is provided.