液体燃料熔盐堆三维瞬态耦合分析

根据熔盐堆的流动特性,采用多群扩散理论,建立液体燃料熔盐堆的三维中子动力学模型和流动传热模型,开发针对液体燃料熔盐堆的三维核热耦合程序,并分析瞬态情况下MOSART堆的物理热工特性.结果表明：堆芯热工参数的瞬态变化直接影响物理参数的变化,液体熔盐堆具有较明显的温度负反馈特性.     

叶栅内定常及非定常粘性流动数值模拟

本文给出了一个模拟叶栅内准三维定常和非定常粘性流动的数值方法。对于定常流动,采用ＴＶＤ Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式和代数湍流模型求解雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,使用当地时间步长和多网格技术使计算加速收敛到定常状态;对于非定常流动,使用双时间步长和全隐式离散,采用与求解定常流动相似的多网格方法求解隐式离散方程。文中给出了ＶＫＩ透平叶栅内的定常流结果和１．５级透平叶栅内的非定常数值结果。     

新型钍基熔盐堆多物理计算模型及分析

针对熔盐堆系统特点,提出了包含堆芯及其他主回路系统在内的多物理紧密耦合计算模型,并在此基础上自主开发了多物理分析程序TANG-MSR。利用该程序进行了新型钍基熔盐堆(TMSR)的设计,并对设计方案进行了稳态及瞬态分析。相关计算结果表明,TANG-MSR所采用的多物理模型能够很好地捕捉熔盐堆的主要物理现象,提出的新型熔盐堆设计在安全性和可持续性方面表现优异。关键词:多物理模型;新型钍基熔盐堆;稳态;瞬态 Abstract: Key words:     

双曲型守恒律的一种高精度TVD差分格式

构造了一维双曲型守恒律方程的一个高精度高分辨率的守恒型TVD差分格式.其主要思想是:首先将计算区域划分为互不重叠的小单元,且每个小单元再根据希望的精度阶数分为细小单元;其次,根据流动方向将通量分裂为正、负通量,并通过小单元上的高阶插值逼近得到了细小单元边界上的正、负数值通量,为避免由高阶插值产生的数值振荡,进一步根据流向对其进行TVD校正;再利用高阶Runge KuttaTVD离散方法对时间进行离散,得到了高阶全离散方法.进一步推广到一维方程组情形.最后对一维欧拉方程组计算了几个算例.     

A central Rankine–Hugoniot solver for hyperbolic conservation laws

A numerical method in which the Rankine–Hugoniot condition is enforced at the discrete level is developed. The simple format of central discretization in a finite volume method is used together with the jump condition to develop a simple and yet accurate numerical method free of Riemann solvers and complicated flux splittings. The steady discontinuities are captured accurately by this numerical method. The basic idea is to fix the coefficient of numerical dissipation based on the Rankine–Hugoniot (jump) condition. Several numerical examples for scalar and vector hyperbolic conservation laws representing the inviscid Burgers equation, the Euler equations of gas dynamics, shallow water equations and ideal MHD equations in one and two dimensions are presented which demonstrate the efficiency and accuracy of this numerical method in capturing the flow features.     

Czochralski晶体生长中非定常热对流过程的数值模拟

本文应用有限控制体算法,采用中心差分格式处理具有Boussinesq近似的N-S方程的对流和扩散通量项,计算了无坩埚和籽晶旋转情况下CZ晶体生长中熔体的非定常自然对流。数值计算的结果得到了某些典型时刻的流场结构和温度分布情况。结合计算得到的CZ熔体总动能和总热能随时间变化的规律,分析了流场和温度场的不稳定变化过程。此外,计算结果还给出了在非稳定热对流影响下的晶-液界面热通量的变化。非定常计算的结果对于研究自然对流和CZ单晶生长,具有理论和实际意义。 关键词：     

CAA broadband noise prediction for aeroacoustic design

The current status of a computational aeroacoustics (CAA) approach to simulate broadband noise is reviewed. The method rests on the use of steady Reynolds averaged Navier-Stokes (RANS) simulation to describe the time-averaged motion of turbulent flow. By means of synthetic turbulence the steady one-point statistics (e.g. turbulence kinetic energy) and turbulent length- and time-scales of RANS are translated into fluctuations having statistics that very accurately reproduce the initial RANS target-setting. The synthetic fluctuations are used to prescribe sound sources which drive linear perturbation equations. The whole approach represents a methodology to solve statistical noise theory with state-of-the-art CAA tools in the time-domain. A brief overview of the synthetic turbulence model and its numerical discretization in terms of the random particle-mesh (RPM) and fast random particle-mesh (FRPM) method is given. Results are presented for trailing-edge noise, slat noise, and jet noise. Some problems related to the formulation of vortex sound sources are discussed.     

金属型脉冲堆惯性效应计算

金属型脉冲堆的反应性负反馈主要来自燃料的热膨胀。当脉冲堆爆发短周期脉冲时,由于燃料的热膨胀不能及时响应温度的变化,致使负反馈滞后,产生惯性效应。采用ANSYS程序计算了在突然热加载条件下的振荡特性,同时采用随机中子输运程序计算了平均关闭系数,然后把带有振荡特性的热反馈带入点堆动力学方程,采用数值方法求解得到了考虑惯性效应的裂变脉冲波形。计算结果与Godiva实验测量结果一致。采用的数值方法中,脉冲周期与振荡效应存在定量关系。由计算结果可知,脉冲周期短,惯性效应明显;脉冲周期越长,惯性效应的影响越小,金属型脉冲堆爆发脉冲时应尽量避免短周期脉冲。     

LDA侧面泵浦薄片激光器的泵浦光和温度分布

建立了激光二极管阵列(LDA)环绕侧面泵浦复合薄片激光器的数值模型,LDA的快轴垂直于薄片表面并被压缩。依据LDA的输出光束特性,考虑到介质与空气的对流热交换和材料热导率的温度相关性,根据经典热传导方程,运用有限单元法,得出了薄片内的泵浦光和温度分布。分析了LDA个数、泵浦距离、吸收系数和光束发散角对薄片内泵浦光分布和吸收效率的影响规律,讨论了温度与泵浦功率、换热系数、冷却液温度和时间的变化规律,所得的有关规律与相关实验相符合。计算结果可为LDA泵浦固体激光器的结构优化设计和实验研究提供理论参考。     

启明星Ⅱ次临界非稳态热中子注量率测量

中子注量率及分布是反应堆的重要参数,本工作通过核数守恒在非稳态情况下的推导和求解,从理论上论证了次临界反应堆非稳态情况下中子注量率测量的可行性。将活化法和固体径迹法有机结合,利用固体径迹探测器标定活化片的测量数据,测量了启明星Ⅱ零功率装置的He-3管实验孔道内及反应堆外壁的中子注量率的分布,并与模拟计算结果进行了比较,利用MCNPX程序得到的模拟计算结果与实验结果的趋势一致,证明了该测量方法可以测量低通量的中子注量率,可实现反应堆不同时刻、不同位置的中子注量率测量,为CiADS技术的研发提供了实验数据与技术支撑。Neutron flux measurements were carried out at VENUS-Ⅱ lead-based zero power reactor by neutron activation method combined with solid-state nuclear track detectors (SSNTD). This experimental method was proposed based on the principle of nuclear number conservation when a foil was irradiated in an unsteady-state neutron field. By this method, thermal neutron flux distributions inside the He-3 duct were measured when VENUS-Ⅱ was operated under unsteady-state. The neutron flux distributions were also calculated with MCNPX code and were consistent with the experimental data. In addition, the neutron fluxes in the outer layer of VENUSⅡ were measured under steady-state. These results would benefit the further study of experimental methods for neutron flux measurement and provide important support for the design of CiADS.     

Efficient stochastic Galerkin methods for random diffusion equations

We discuss in this paper efficient solvers for stochastic diffusion equations in random media. We employ generalized polynomial chaos (gPC) expansion to express the solution in a convergent series and obtain a set of deterministic equations for the expansion coefficients by Galerkin projection. Although the resulting system of diffusion equations are coupled, we show that one can construct fast numerical methods to solve them in a decoupled fashion. The methods are based on separation of the diagonal terms and off-diagonal terms in the matrix of the Galerkin system. We examine properties of this matrix and show that the proposed method is unconditionally stable for unsteady problems and convergent for steady problems with a convergent rate independent of discretization parameters. Numerical examples are provided, for both steady and unsteady random diffusions, to support the analysis.     

Chemical composition of the lunar surface from neutron leakage fluxes

The neutron leakage fluxes from the lunar surface are calculated by Monte Carlo transport code based on Geant4. The integral fluxes of fast neutrons, epi-thermal neutrons and thermal neutrons are analyzed. Numerical results for 20 kinds of lunar soils and 7 minerals show that the fast neutron fluxes are linearly related to the average atomic mass numbers of the lunar materials used in simulations. Meanwhile, the average atomic mass numbers are strongly modulated by the abundances of iron (Fe) and titanium (Ti), and a linear relationship between the average atomic mass numbers and the abundances of Fe and Ti is found. Furthermore, the results show that the ratios of epi-thermal to thermal fluxes for lunar soils are linearly related to the macroscopic absorption cross sections of lunar materials, and that the macroscopic absorption cross sections monotonically increase with the abundances of Fe and Ti by a simple function. Then we reach the conclusion that the neutron fluxes can provide the information about the Fe and Ti contents.     

Mathematical Development and Verification of a Finite Volume Model for Morphodynamic Flow Applications

The accuracy and efficiency of a class of finite volume methods are investigated for numerical solution of morphodynamic problems in one space dimension. The governing equations consist of two components, namely a hydraulic part described by the shallow water equations and a sediment part described by the Exner equation. Based on different formulations of the morphodynamic equations, we propose a family of three finite volume methods. The numerical fluxes are reconstructed using a modified Roe's scheme that incorporates, in its reconstruction, the sign of the Jacobian matrix in the morphodynamic system. A well-balanced discretization is used for the treatment of the source terms. The method is well-balanced, non-oscillatory and suitable for both slow and rapid interactions between hydraulic flow and sediment transport. The obtained results for several morphodynamic problems are considered to be representative, and might be helpful for a fair rating of finite volume solution schemes, particularly in long time computations.     

Prediction of the thermophysical properties of molten salt fast reactor fuel from first-principles

Molten fluorides are known to show favourable thermophysical properties which make them good candidate coolants for nuclear fission reactors. Here we investigate the special case of mixtures of lithium fluoride and thorium fluoride, which act both as coolant and as fuel in the molten salt fast reactor concept. By using ab initio parameterised polarisable force fields, we show that it is possible to calculate the whole set of properties (density, thermal expansion, heat capacity, viscosity and thermal conductivity) which are necessary for assessing the heat transfer performance of the melt over the whole range of compositions and temperatures. We then deduce from our calculations several figures of merit which are important in helping the optimisation of the design of molten salt fast reactors.     

扩散方程的有限方向差分方法

研究二维散乱点集上数值求解非线性扩散方程的有限方向差分方法。利用五个邻点信息构造具有最小模板的离散格式,并且离散系数具有显式表达式。另外,利用五点公式获得了间断问题物质界面的离散格式,该格式对界面流的计算具有近似二阶精度。不同计算区域及不同类型的离散点集上的计算结果验证了方法的有效性。     

铅冷快堆M2LFR-1000堆芯燃料管理方案设计

堆芯燃料管理是反应堆设计中极为重要而且复杂的工作,直接影响着堆芯的经济性。目前国内外对于压水堆等传统热堆已有了较为丰富和成熟的燃料管理计算方法,但对于快堆,由于其中子能谱硬,与传统热堆相比有着不同的控制方式和功率分布,快堆的堆芯燃料管理缺乏系统研究。针对中国科学技术大学自主研发的强迫循环冷却的铅基快堆M2LFR-1000,应用SRAC/COREBN软件包进行堆芯燃耗计算,根据燃耗深度提取核素核子密度,计算伪平衡循环参数进行燃料管理预估,然后进行首循环装料、过渡循环和平衡循环燃料管理方案设计。结果表明:对M2LFR-1000堆芯外区燃料换料组件Pu的富集度进行优化,可以延长换料周期到540 d,提高平均卸料燃耗深度;伪平衡循环结果与平衡循环基本一致,伪平衡循环可以用于燃料管理预估。     

A reconstructed discontinuous Galerkin method for the compressible Navier–Stokes equations on arbitrary grids

A reconstruction-based discontinuous Galerkin (RDG) method is presented for the solution of the compressible Navier–Stokes equations on arbitrary grids. The RDG method, originally developed for the compressible Euler equations, is extended to discretize viscous and heat fluxes in the Navier–Stokes equations using a so-called inter-cell reconstruction, where a smooth solution is locally reconstructed using a least-squares method from the underlying discontinuous DG solution. Similar to the recovery-based DG (rDG) methods, this reconstructed DG method eliminates the introduction of ad hoc penalty or coupling terms commonly found in traditional DG methods. Unlike rDG methods, this RDG method does not need to judiciously choose a proper form of a recovered polynomial, thus is simple, flexible, and robust, and can be used on arbitrary grids. The developed RDG method is used to compute a variety of flow problems on arbitrary meshes to demonstrate its accuracy, efficiency, robustness, and versatility. The numerical results indicate that this RDG method is able to deliver the same accuracy as the well-known Bassi–Rebay II scheme, at a half of its computing costs for the discretization of the viscous fluxes in the Navier–Stokes equations, clearly demonstrating its superior performance over the existing DG methods for solving the compressible Navier–Stokes equations.     

LDA侧面抽运棒状激光器热透镜效应的有限单元法分析

建立了激光二极管阵列端面抽运棒状激光介质的数值模型.考虑到介质与空气的对流换热和介质的热力学参量随温度的变化,根据经典热传导方程和热弹性方程,运用有限单元法,得出了棒状介质内瞬态温度、热应力和应变的分布.分析了热透镜焦距随抽运功率的变化规律,所得的规律与有关文献相符合.理论分析结果可为激光二极管阵列抽运固体激光器的结构优化设计和实验提供理论参考.     

基于CFD/CSD的舵面控制气动弹性数值模拟方法

在动态网格上通过耦合求解流动控制方程和结构动力学方程, 发展了一种舵面控制下飞行器运动响应过程中气动弹性数值模拟研究方法.流动控制方程采用N-S方程, 结构动力学采用线性模态叠加方法, 其中流动控制方程空间离散采用基于非结构网格的有限体积方法, 对流通量采用计算HLLC格式, 非定常时间离散采用基于LU-SGS的双时间步长方法.模拟中, 气动运动和结构变形在双时间步长方法推进过程中采用改进松耦合方法, 气动网格与结构网格之间信息交换采用无限平板样条法实现, 飞行器的运动和变形采用基于重叠网格和Delaunay图映射变形网格相结合的方法进行处理.采用多个考核算例对发展的数值方法进行考核验证, 结果表明该方法可以高效精确模拟舵面开环控制下飞行器运动响应过程中的气动弹性特性.