一维守恒的Lagrangian ADER格式

基于欧拉框架下ADER格式,构造一维守恒只有一个时间步的、高精度中心型拉格朗日ADER(LADER)格式.构造r阶LADER格式包括:从欧拉方程出发推导拉格朗日框架下积分形式的方程、采用WENO方法高精度重构节点处守恒量和从1阶到r-1阶的空间导数、求拉氏框架下这些变量的Godunov值,并计算1阶到r-1阶的时间全导数,最后高精度离散积分形式的流通量函数.对光滑流场的模拟表明,LADER格式达到设计的精度;对含强间断的流场模拟表明,数值解在间断附近基本无振荡.     

Fully HOC Scheme for Mixed Convection Flow in a Lid-Driven Cavity Filled with a Nanofluid

A fully higher-order compact (HOC) finite difference scheme on the 9-point two-dimensional (2D) stencil is formulated for solving the steady-state laminar mixed convection flow in a lid-driven inclined square enclosure filled with water-$Al_2O_3$ nanofluid. Two cases are considered depending on the direction of temperature gradient imposed (Case I, top and bottom; Case II, left and right). The developed equations are given in terms of the stream function-vorticity formulation and are non-dimensionalized and then solved numerically by a fourth-order accurate compact finite difference method. Unlike other compact solution procedure in literature for this physical configuration, the present method is fully compact and fully higher-order accurate. The fluid flow, heat transfer and heat transport characteristics were illustrated by streamlines, isotherms and averaged Nusselt number. Comparisons with previously published work are performed and found to be in excellent agreement. A parametric study is conducted and a set of graphical results is presented and discussed to elucidate that significant heat transfer enhancement can be obtained due to the presence of nanoparticles and that this is accentuated by inclination of the enclosure at moderate and large Richardson numbers.     

分辨率优化的混合WENO格式

为提高有限差分格式的分辨率,利用傅里叶分析对WENO格式进行色散及耗散优化,并给出优化的线性权重.用优化后的WENO格式与保单调格式（MP）进行加权混合,得到新的加权混合WENO格式（H-WENO）.通过一维激波管问题、Shu-Osher问题及二维双Mach反射问题及R-T不稳定性问题对格式进行数值测试.结果显示,新格式具有强健的激波捕捉能力和对小尺度波结构的高分辨率,与原WENO格式相比改进明显.     

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.     

轴对称射流气动声场的数值模拟

构造了Kirchhoff积分和CFD计算相结合的算法,采用高精度差分格式对不同喷口马赫数的轴对称射流进行数值模拟,并以此作为近场声源,运用Kirchhoff积分方法研究远场气动噪声.数值结果表明,远场噪声具有方向性,噪声声压在离开对称轴20°处达到最大值.随着传播距离增大,噪声方向性逐渐减弱.     

Simulating high Reynolds number flow in two-dimensional lid-driven cavity by multi-relaxation-time lattice Boltzmann method

By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number (Re) flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50,000 in the literature; however, this paper can readily extend the maximum to 1000,000 with the above combination.     

计算机仿真环境下电子双缝衍射的实验方案

在一般的近代物理实验教材中，关于电子衍射实验采用的是晶体表面反射衍射或金属箔透射衍射。本文提出一套基于双缝衍射的、用于研究微观粒子波粒二象性的电子衍射实验方案，并且鉴于在一般的实验中不具备进行电子双缝衍射实验的条件，本实验方案是在电子双缝衍射计算机模拟的仿真环境下进行的。     

Evidence for shape coexistence at medium spin in Rb

Four previously known rotational bands in ⁷⁶Rb have been extended to moderate spins using the Gammasphere and Microball γ ray and charged particle detector arrays and the ⁴⁰Ca(⁴⁰Ca, 3pn) reaction at a beam energy of 165 MeV. The properties of two of the negative-parity bands can only readily be interpreted in terms of the highly successful Cranked Nilsson-Strutinsky model calculations if they have the same configuration in terms of the number of g9/2 particles, but they result from different nuclear shapes (one near-oblate and the other near-prolate). These data appear to constitute a unique example of shape coexisting structures at medium spins.     

The fine structure levels and spin-singlet contributions to zero-field-splitting parameters of Cr2+ ion in CdGa2S4

A theoretical investigation is reported of the fine structure levels and the spin-singlet contributions to zero-field-splitting (ZFS) parameters for Cr²⁺ ion in CdGa₂S₄ crystals. Firstly, the complete energy matrix including all spin states for a 3d ⁴ ion in tetrahedral D _2d symmetry is constructed according the double-group chain in the strong-field scheme. Then, by diagonalizing the complete energy of electron–electron interactions, the crystal field and the spin–orbit coupling for the Cr²⁺ (3d ⁴) ion in CdGa₂S₄ crystal, the fine structure levels and the spin-singlet contributions to ZFS parameters a, D and F are calculated. The results show that the spin-singlet contribution to D is negligible, but the contributions to a and F are very important. So, to obtain more accurate ZFS parameters for 3d ⁴ ions in the tetrahedral crystals, all spin states should be considered.     

A numerical study of the motion of a wave running up a beach

The main difficulty for the numerical calculation of the wave running up a beach is the treatment of its moving water boundary. In this paper a scheme of turning the free boundary problem into a fixed boundary problem is designed. The calculated run-up height is consistent with the experiments. Some interesting wave phenomena are also found.