核工程仿真软件：COMSOL二次开发_（2）.核反应堆物理与热工水力学仿真.docx

PAGE1

PAGE1

核反应堆物理与热工水力学仿真

1.核反应堆物理仿真基础

核反应堆物理仿真是核工程领域中一个重要的分支，它涉及到反应堆中中子输运、裂变、吸收等物理过程的模拟。COMSOLMultiphysics软件虽然不是专门用于核反应堆物理仿真的工具，但通过二次开发，可以实现对这些过程的精确模拟。本节将介绍如何在COMSOL中进行核反应堆物理仿真的基础内容。

1.1中子输运方程

中子输运方程是描述中子在核反应堆中行为的基本方程。它考虑了中子的产生、吸收、散射和泄漏等过程。方程的形式如下：

$$

+(r,E,t)+_t(r,E)(r,E,t)=S(r,E,t)+_s(r,E’E)(r,E’,t)dE’+(r,E)_f(r,E’)(r,E’,t)dE’

$$

其中：

?r

jr

Σt

Σs

Σf

ν是每个裂变中子的平均数。

χr

Sr

1.1.1中子输运方程在COMSOL中的实现

在COMSOL中，可以通过PDE接口来实现中子输运方程。以下是一个简单的例子，展示如何在COMSOL中设置和求解中子输运方程。

例子：一维中子输运方程

假设我们有一个一维的简单反应堆模型，中子输运方程可以简化为：

$$

++_t(x,E)(x,E,t)=S(x,E,t)+_s(x,E’E)(x,E’,t)dE’+(x,E)_f(x,E’)(x,E’,t)dE’

$$

在COMSOL中，我们可以通过以下步骤来实现：

创建几何模型：

定义一个一维的几何模型，例如一个长度为1米的反应堆。

定义材料属性：

设置材料的中子宏观截面Σt,Σs,Σf和裂变中子能谱

设置边界条件：

定义反应堆的边界条件，例如反射边界或吸收边界。

编写PDE：

使用COMSOL的PDE接口来编写中子输运方程。

代码示例

#导入COMSOL接口

importcomsol

#创建模型

model=comsol.Model(NeutronTransport1D)

#定义几何模型

geometry=model.geometry()

geometry.add_interval(0,1)#定义一个长度为1米的一维区间

#定义材料属性

materials=model.materials()

fuel=materials.add(Fuel)

fuel.property(Sigma_t,10/meter)

fuel.property(Sigma_s,8/meter)

fuel.property(Sigma_f,2/meter)

fuel.property(nu,2.5)

fuel.property(chi,1)

#定义PDE接口

pde=model.physics.add(PDE,NeutronTransport)

pde.depending_variable(phi)

#设置PDE方程

pde.equation(1/v*d(phi(x,E,t),t)+d(j(x,E,t),x)+Sigma_t*phi(x,E,t)=S+intop(Sigma_s*phi(x,E1,t),E1)+chi*intop(nu*Sigma_f*phi(x,E1,t),E1))

#定义边界条件

boundary=model.physics.boundary()

boundary.add(Reflective,left)

boundary.add(Absorbing,right)

#设置求解器

solver=model.solver()

solver.add(TimeDependent)

solver.add(Stationary)

#求解

model.solve()

#输出结果

results=model.results()

results.plot(phi,NeutronFlux)

1.2裂变和吸收过程

裂变和吸收过程是核反应堆物理仿真的核心内容。裂变过程涉及中子引发的核裂变反应，而吸收过程则涉及中子被材料吸收的物理过程。在COMSOL中，可以通过定义相应的源项和吸收项来模拟这些过程。

1.2.1裂变过程的模拟

裂变过程可以通过定义裂变源项来实现。裂变源项的形式为：

$$

S_f=_f

$$

其中：

ν是每个裂变中子的平均数。

Σf

?