中子扩散方程高阶谐波可用于重构堆芯中子注量率分布,但传统源迭代与源修正迭代法求解时的收敛速度慢,计算耗时长。采用隐式重启Arnoldi方法(Implicitly Restarted Arnoldi Method,IRAM)求解本征值问题的中子扩散方程获得谐波数据,通过本征正交分解(Proper Orthogonal Decomposition,POD)与伽辽金(Galerkin)投影相结合的方法构建POD-Galerkin低阶模型,并重构二维稳态TWIGL基准题中子注量率分布。研究结果表明:IRAM方法在求解中子扩散方程的高阶本征值和谐波问题上具有较高的精度;基于POD-Galerkin低阶模型重构中子注量率分布具有较高的保真性与计算效率,有效增值系数与参考解的误差为8.7×10-5,对角线上快群和热群中子注量率最大相对误差为2.56%,且低阶模型计算用时仅为全阶模型的10.18%。本研究为堆芯中子注量率重构提供了一种可靠且高效的方法,该方法不仅可用于重构稳态时堆芯中子注量率分布,还具有在瞬态情况下预测中子注量率分布的潜力,有望在未来的应用中进一步拓展。
中子扩散方程 隐式重启Arnoldi方法 本征正交分解 伽辽金投影 中子注量率重构 Neutron diffusion equation Implicitly restarted Arnoldi method Proper orthogonal decomposition Galerkin projection Neutron flux reconstruction
烟台大学物理与电子信息学院,山东 烟台 264005
生物表面组织是多层介质,目前只有光在单层介质中传输的P3方程,本文给出了光在两层平板介质中传播的P3方程的时域解,编写蒙特卡罗模拟进行验证。结果表明P3时域方程与蒙特卡罗模拟符合很好,说明双层平板介质的P3方程正确地反映了光子在介质中的迁移。与扩散时域方程相比,在峰值位置附近,P3时域方程测量的反射率和透射率更加准确,在远离峰值时,P3方程与扩散方程结果一致。双层平板介质的P3时域方程扩展了P3方程的应用范围,可以替代扩散方程,也可以代替单层平板介质的P3方程,可以更好地研究光在生物组织中的传播过程。
光子迁移 P3近似 扩散方程 蒙特卡罗模拟 光学学报
2023, 43(13): 1320001
天津大学精密测试技术及仪器国家重点实验室, 天津 300072
近红外漫反射光谱法可用于物质成分定性、 定量分析, 具有快速、 无损检测等优点, 在食品、 制药、 环境监测、 生命科学等领域有广泛应用。 采用近红外漫反射光谱对浑浊散射介质进行成分定量分析时, 由于光散射作用, 介质吸光度随测量距离不再呈线性变化, 不同测量距离下测量模型也很难转换。 从扩散方程出发, 提出了从漫射光谱中提取线性吸光度的方法, 获得了与测量距离无关的介质光学参数——有效衰减系数的光谱, 该光谱可反映介质吸收变化, 并用于物质成分测量。 具体采用了双位置差分法, 可在线性吸光度的范围内任意选取两个光源-探测器距离测量, 通过差分获得与测量位置无关的有效衰减系数光谱, 其中测量距离可依据波段灵活地进行选择。 另外, 经理论推导可知, 差分运算的同时还大大削减了光子扩散对测量的影响, 有利于不同散射特性的介质间光谱模型的相互借鉴。 经理论计算、 蒙特卡洛模拟和实验验证, 在1 000~1 360 nm波段, 该方法适用于强散射的浑浊水溶液。 对散射系数处于较大范围(28.53~87.47 cm-1)的intralipid水溶液进行了测试, 获得了线性吸光度的光源-探测器距离范围。 进一步以葡萄糖为被测成分为例, 获得了葡萄糖在intralipid水溶液中的有效衰减系数光谱。 采用测量距离连续可调的漫反射光谱测量系统进行了实验, 测量范围位于0~0.3 cm, 测试了3种典型的介质(3%, 5%, 10%Intralipid水溶液)。 结果可知, 3种介质中葡萄糖光谱呈现一致性, 它们之间的测量模型可以方便地进行线性转换。 研究对象intralipid水溶液常作为生物组织的仿体, 其光学参数也覆盖了常见的牛奶、 果汁等被测浑浊溶液, 因此测量结果具有广泛的适应性。 综上, 该研究实现了对浑浊水溶液中物质成分线性吸光度的提取, 该方法不仅使得测量位置的选择更加灵活方便, 还有利于实现不同散射介质间模型的线性转换, 特别适用于介质成分复杂多变的场合下的测量应用, 如对人体组织、 牛奶或食品中的成分测量等。
近红外光谱法 散射介质 扩散方程 蒙特卡洛模拟 线性吸光度 有效衰减系数 Near infrared spectroscopy Scattering medium Diffusion equation Monte Carlo simulation Linear attenuance Effective attenuation coefficient 光谱学与光谱分析
2019, 39(12): 3711
杭州电子科技大学 信息工程学院, 杭州 310018
引入核函数法对随机扩散方程(SDE)样本的密度分布进行统计, 希望用核函数来减少统计涨落。由于SDE样本的密度随时间发展, 越来越稀疏, 所以核函数也应该越来越大, 也就是说核函数应该随时间在变化。通过一个瞬时释放的二维扩散问题(具有解析解), 从定性和定量两个角度比较了变带宽核函数法和传统统计方法在密度分布统计中的性能差别, 论述了变带宽核函数法的优缺点, 变带宽核函数法在牺牲部分峰值的前提下可以很好地解决SDE样本密度分布统计涨落问题, 在工程应用中值得推广。
随机扩散方程 扩散输运方程 拉格朗日模型 核函数法 stochastic diffusion equation advection-diffusion transport equation Lagrangian model kernel function method 强激光与粒子束
2017, 29(12): 126009
中国核动力研究设计院 核反应堆系统设计技术重点实验室, 成都 610041
三维多群中子扩散方程的精确、高效求解是核动力堆芯设计及燃料管理的基础。应用有限差分方法求解该方程具有简便、精确、成熟的优点; 然而,该方法的计算量和存储量均较大,极大地限制了它的计算规模和应用范围。本文基于大规模并行计算,研究三维多群中子扩散方程有限差分方法:采用中心有限差分格式离散中子扩散方程; 基于MPI并行编程模型,采用空间区域分解的方式实现大规模并行计算; 采用多群多区域耦合PGMRES算法进行并行加速。在集群服务器上开发了ParaFiDi程序,并采用IAEA3D,PHWR等多个基准题对该程序进行验证。数值结果表明,ParaFiDi程序具有较高的计算精度和计算效率。
三维多群中子扩散方程 有限差分方法 大规模并行计算 区域分解算法 多群多区域耦合PGMRES算法 3D multi-group neutron diffusion equation finite difference method large-scale parallel computation domain decomposition algorithm multi-group multi-domain coupled PGMRES algorithm 强激光与粒子束
2017, 29(8): 086001
清华大学 核能与新能源技术研究院, 先进核能技术协同创新中心, 先进反应堆工程与安全教育部重点实验室, 北京 100084
研究了JFNK框架下高温堆中子扩散问题的求解方法。研究结果表明,JFNK方法在求解与源迭代相同形式中子扩散方程时,相对残差下降趋势为逐渐加快并趋于稳定,有利于更高求解精度的实现。使用通量归一化附加方程可以获得更好的JFNK非线性迭代特性,但在算例中其部分牛顿修正方程求解时间偏多,总计算时间高于显式有效增殖系数附加方程法,需要研究更高效的JFNK预处理方法对线性求解环节进行改善。
高温堆 中子扩散方程 附加方程 源迭代 high temperature reactor neutron diffusion equation additional equation source iteration JFNK JFNK 强激光与粒子束
2017, 29(3): 036024
1 安徽理工大学理学院, 淮南 232001
2 中国科学技术大学材料科学与工程系, 合肥 230026
我们利用相干态表象将经典扩散方程自然地推广到量子力学形式, 并讨论了相干态在扩散通道中的演化问题, 发现相干纯态在扩散通道中会演化为混态, 并给出任意时刻的终态表达式。同时, 我们得到光子数在扩散通道中会增加。
相干态 扩散方程 光子数 Coherent State Diffusion Equation Photon Number
北京应用物理与计算数学研究所, 北京 100094
针对惯性约束聚变领域需要求解的含强间断导热系数的扩散方程,借鉴各向异性导热的思想,给出了一种导热系数定义在网格节点的支撑算子改进格式。改进后的方法保留了支撑算子格式的算子相容性优点。数值算例表明,新方法可以很好地处理含强间断和强非线性的热扩散问题,计算结果保持较高的精度,适合于辐射流体力学问题的模拟。with strongly discontinuous coefficients
惯性约束聚变 扩散方程 支撑算子格式 强间断 强非线性 inertial confinement fusion diffusion equation support operator method strongly discontinuous coefficients strong nonlinearity 强激光与粒子束
2015, 27(9): 092014
1 西安电子科技大学 物理与光电工程学院, 西安 710071
2 桂林电子科技大学 生命与环境科学学院, 广西 桂林 541004
3 上海卫星工程研究所, 上海 200240
针对复杂背景下红外弱小目标检测效果不佳的问题, 结合多尺度分析法和各向异性扩散方程, 利用图像尺度和方向信息, 提出一种新的红外弱小目标检测算法.首先, 采用Surfacelet变换对图像进行多尺度、多方向分解, 得到不同尺度下高频子带系数和低频子带系数;其次根据不同频带的特点, 分别采用改进的各向异性扩散方程差分滤波和局部去均值滤波对高频子带系数和低频子带系数进行处理;最后, 采用逆Surfacelet变换重构处理后的子带系数, 并采用自适应阈值分割对重构的图像进行分割, 以实现目标检测.采用多组红外图像进行试验, 并与小波滤波以及各向异性扩散滤波进行比较, 实验结果显示, 该算法能有效抑制背景及其边缘, 可以获得比另外两种算法更好的红外弱小目标检测效果.
红外与夜视技术 目标检测 多尺度分析 图像处理 弱小目标 图像滤波 各向异性扩散方程 阈值分割 Infrared night vision technique Target detection Multiscale analysis Image processing Dim and small target Image filtering Anisotropic diffusion equation Threshold segmentation