Higher-order approximate solutions of fractional stochastic point kinetics equations in nuclear reactor dynamics

Stochastic point kinetics equations(SPKEs) are a system of Ito? stochastic differential equations whose solution has been obtained by higher-order approximation.In this study, a fractional model of SPKEs has been analyzed. The efficiency of the proposed higher-order approximation scheme has been discussed in the results section. The solutions of SPKEs in the presence of Newtonian temperature feedback have also been provided to further discuss the physical behavior of the fractional model.

Developed mathematical technique for fractional stochastic point kinetics model in nuclear reactor dynamics

Fractional stochastic kinetics equations have proven to be valuable tools for the point reactor kinetics model, where the nuclear reactions are not fully described by deterministic relations. A fractional stochastic model for the point kinetics system with multi-group of precursors,including the effect of temperature feedback, has been developed and analyzed. A major mathematical and inflexible scheme to the point kinetics model is obtained by merging the fractional and stochastic technique. A novel split-step method including mathematical tools of the Laplace transforms, Mittage–Leffler function, eigenvalues of the coefficient matrix, and its corresponding eigenvectors have been used for the fractional stochastic matrix differential equation. The validity of the proposed technique has been demonstrated via calculations of the mean and standard deviation of neutrons and precursor populations for various reactivities: step, ramp, sinusoidal, and temperature reactivity feedback. The results of the proposed method agree well with the conventional one of the deterministic point kinetics equations.

Analytical solution of modified point kinetics equations for linear reactivity variation in subcritical nuclear reactors adopting an incomplete gamma function approximation

The present work aims to achieve a fast and accurate analytical solution of the point kinetics equations applied to subcritical reactors such as ADS (Accelerator-Driven System), assuming a linear reactivity and external source variation. It was used a new set of point kinetics equations for subcritical systems based on the model proposed by Gandini & Salvatores. In this work it was employed the integrating factor method. The analytical solution for the case of interest was obtained by using only an approximation which consists of disregarding the term of the second derivative for neutron density in relation to time when compared with the other terms of the equation. And also, it is proposed an approximation for the upper incomplete gamma function found in the solution in order to make the computational processing faster. In addition, for purposes of validation and comparison a numerical solution was obtained by the finite differences method. Finally, it can be concluded that the obtained solution is accurate and has fast numerical processing time, especially when compared with the results of numerical solution by finite difference. One can also observe that the gamma approximation used achieve a high accuracy for the usual parameters. Thus we got satisfactory results when the solution is applied to practical situations, such as a reactor startup.

一类分数阶随机微分方程的均方渐近概周期解

关于分数阶随机微分方程解的性质研究是近几年数学界的热门方向之一。针对Hilbert空间上一类线性分数阶随机微分方程,研究其均方渐近概周期温和解的存在性和唯一性,然后将这类线性分数阶随机微分方程的结论推广到对应的半线性分数阶随机微分方程中,利用Banach不动点定理讨论这类半线性分数阶随机微分方程均方渐近概周期温和解的存在唯一性,再利用Schauder不动点定理讨论这类方程在非Lipschitz条件下均方渐近概周期温和解的存在性。

Stochastic Differential Equations Driven by Multi-fractional Brownian Motion and Poisson Point Process

In this paper,we study a class of stochastic differential equations with additive noise that contains a non-stationary multifractional Brownian motion(mBm)with a Hurst parameter as a function of time and a Poisson point process of class(QL).The differential equation of this kind is motivated by the reserve processes in a general insurance model,in which between the claim payment and the past history of liability present the long term dependence.By using the variable order fractional calculus on the fractional Wiener-Poisson space and a multifractional derivative operator,and employing Girsanov theorem for multifractional Brownian motion,we prove the existence of weak solutions to the SDEs under consideration,As a consequence,we deduce the uniqueness in law and the pathwise uniqueness.

APPROXIMATE CONTROLLABILITY OF FRACTIONAL IMPULSIVE NEUTRAL STOCHASTIC INTEGRO-DIFFERENTIAL EQUATIONS WITH NONLOCAL CONDITIONS AND INFINITE DELAY

This paper is concerned with the approximate controllability of nonlinear fractional impulsive neutral stochastic integro-differential equations with nonlocal conditions and infinite delay in Hilbert spaces under the assumptions that the corresponding linear system is approximately controllable. By the Krasnoselskii-Schaefer-type fixed point theorem and stochastic analysis theory, some sufficient conditions are given for the approximate controllability of the system. At the end, an example is given to illustrate the application of our result.

临界反应堆阶跃正反应性输入时中子密度响应的近似修正解

通过修正单组缓发中子先驱核衰变常量A值，使点堆中子动力学方程单组缓发中子模型在正反应性阶跃输入时的数值计算结果趋近于六组缓发中子模型数值计算结果。在此基础上，用修正后的单组模型解析方法进行计算。结果表明：采用修正后的单组解析方法计算阶跃正反应性输入的中子密度响应，计算结果与六组的接近，满足工程计算精度要求，同时计算简便，避免了刚性问题，可以实现快速计算。

高次端点浮动法—解点堆中子动力学方程

本文介绍了解点堆中子动力学方程的高次多项式近似的“端点浮动法”，该方法不但具有克服方程刚性和灵活应用的优点，而且能够推广应用于刚性不显著或消失时的大反应性情况。由于近似多项式次数的提高、截断误差的减少、计算的精度有进一步的提高。

ADS次临界反应堆的点堆中子动力学方程

加速器驱动的次临界系统(ADS)中的次临界反应堆与临界反应堆相比,中子注量率的空间分布具有严重的不均匀性,同时中子平均能量较高且中子能量变化复杂,中子价值变化大,因而传统的点堆动力学方程不能较为真实地模拟ADS次临界反应堆。本文从含多群中子多组缓发中子先驱核的动力学方程出发,给出其共轭方程。然后利用稳态扩散方程及其共轭方程的共轭关系,推导得出含有归一化功率的动力学方程表达式。进而定义多个特征算子,导出了含有源中子价值的点堆中子动力学方程,并对几种简单情况进行了初步验证,为进一步分析ADS次临界反应堆的动态过程奠定了基础。

单组点堆中子动力学方程的解析解

针对有外加中子源单组点堆中子动力学方程组的快速求解问题,提出了常数变易法和同伦分析方法,通过理论推导获得了单组点堆中子动力学方程组解析解,然后对3种不同刚性情况下的阶跃反应性引入的例子进行了求解和数值分析,验证了所得解析解的正确性。

零功率反应堆传递函数的实验测量

为研究零功率反应堆的传递函数,先以点堆动态方程为出发点,推导零功率反应堆对外加反应性扰动输入响应的传递函数理论表达式。然后,对外加反应性扰动为三角波时进行傅里叶分析,得到传递函数与输入频率及反应性之间的关系。最后,选取DF-Ⅵ快中子临界装置进行验证实验。在测量了扰动的反应性后,给出确定频率的外加扰动为三角波时的传递函数实验值,实验值与传递函数的理论计算值相符合,证实了零功率反应堆传递函数的理论表达式是正确的,反之,也可从已知的传递函数求出扰动的反应性。

点堆随机动力学方程伊藤解与中子数概率分布解析解的对比分析

为了研究反应堆弱中子源启动过程中的中子数密度和缓发中子先驱核随机涨落现象,我们推导和建立了点堆随机动力学方程组,在传统的点堆动力学方程组中引入了伊藤随机项。为了验证方程组的伊藤解方法和计算精度,我们在简化物理条件和方程形式下,对定态系统的中子数密度分别用随机动力学方程伊藤解和中子数概率分布函数解析解进行了对比分析。结果表明,伊藤解是一种有效、具有较高计算精度的方法,计算精度满足 sigma 的标准,置信水平在 95%以上。

一类中立型随机偏微分方程概周期解的存在性和唯一性

最近,P.Bezandry和T.Diagana(P.Bezandry,T.Diagana.Appl.Anal.,2007,117:1-10.)引入了均值概周期的概念,研究了一类随机时滞演化方程,并获得了其均值概周期存在和唯一的充分条件.我们将应用不动点理论和分数幂算子方法,获得一类中立型随机偏微分方程在均方意义下的概周期解的存在性和唯一性的充分条件.

分数阶脉冲中立型随机微积分方程的适定性

利用Sadovskii不动点定理以及α-预解算子理论讨论了一类在Hilbert空间中带无限时滞的分数阶脉冲中立型随机微积分方程温和解的适定性,并通过举例说明了结果的有效性.

基于可视化方法的点堆动力学方程刚性分析

为了更清晰地表现方程的刚性,采用可视化方法对点堆动态方程进行仿真,研究和分析该方程的刚性性质.结果表明,在可视化条件下,能够较为直观地了解点堆动力学方程的刚性性质,有助于更加全面认识和掌握该方程.

分数阶混合随机泛函微分方程的能控性

主要研究了分数阶混合随机泛函微分方程的能控性.在无限维空间下,假设所考虑方程线性部分生成半群不是紧的,使用非紧性测度技术和Mnch不动点定理,给出了方程能控性充分条件,并通过一个例子说明了结论的有效性.

优化时间步长的数值方法解核反应堆点动态学方程

通过分析步长与数值误差的关系,论证了调整和优化步长的必要性。从步长与误差的关系出发,介绍了数值求解初值微分方程问题的两种优化步长的方法:两步计算法和嵌入式龙格-库塔法。在用隐式欧拉法对反应堆点动态学方程进行数值求解时,采用两步计算法对步长进行优化,计算的结果表明控制误差对步长及数值解的精确性有决定性的影响;通过求解反应堆点动态学方程,对MATLAB自带的嵌入式龙格-库塔法的微分方程求解函数ode23及ode45的精确性及使用情况进行了探讨。

Existence of p-mean Almost Periodic Mild Solution for Fractional Stochastic Neutral Functional Differential Equation

A class of fractional stochastic neutral functional differential equation is analyzed in this paper.With the utilization of the fractional calculations,semigroup theory,fixed point technique and stochastic analysis theory,a sufficient condition of the existence for p-mean almost periodic solution is obtained,which are supported by two examples.

由分数布朗运动驱动的带跳的随机微分方程的解

论述了由分数布朗运动驱动的带跳的随机微分方程的系数满足Lipschitz条件时解的存在性和唯一性。

系统瞬态分析程序中子学模块点堆求解验证

系统瞬态分析程序是用于示范快堆设计的工具软件之一,由中国原子能科学研究院自主开发,中子学模块是此程序中用于计算反应堆功率的模块。为了确保软件计算的正确性,需要对其进行验证与确认(V&V),本文对中子学模块求解点堆计算功率的结果进行了验证,是整个程序V&V工作中的一个部分。结果显示,在不同的条件下,中子学模块计算的结果与多种数值方法符合较好,满足程序计算的需求。