核燃料循环软件：ORIGEN二次开发_（6）.活度与衰变热计算.docx

PAGE1

PAGE1

活度与衰变热计算

活度计算的基本公式

活度计算的基本公式如下：

A

其中：

At是在时间t

Nit是在时间t时第i

λi是第i

N是核素的总数。

核素衰变链

在多核素系统中，核素之间的衰变链关系复杂，需要使用衰变链矩阵来表示。衰变链矩阵B中的元素Bij表示第i种核素通过衰变生成第j种核素的速率。衰变链矩阵的对角线元素B

活度计算的数值方法

活度计算通常采用数值方法，如欧拉法、龙格-库塔法等。这些方法可以有效地解决复杂的衰变链问题。以下是一个简单的欧拉法示例：

importnumpyasnp

#定义时间步长和总时间

dt=0.1#时间步长（秒）

total_time=100#总时间（秒）

#定义初始核素数量

N0=np.array([100,0,0])#初始核素数量[N1,N2,N3]

#定义衰变常数

lambda1=0.1#第一种核素的衰变常数（1/秒）

lambda2=0.2#第二种核素的衰变常数（1/秒）

lambda3=0.3#第三种核素的衰变常数（1/秒）

#定义衰变链矩阵

B=np.array([

[-lambda1,0,0],

[lambda1,-lambda2,0],

[0,lambda2,-lambda3]

])

#定义时间步长函数

defeuler_step(N,B,dt):

returnN+dt*B@N

#计算活度

N=N0

A=np.zeros(total_time//dt)

fortinrange(0,total_time,int(dt)):

A[t//dt]=np.sum(N*np.array([lambda1,lambda2,lambda3]))

N=euler_step(N,B,dt)

#输出活度随时间的变化

print(A)

代码解释

时间步长和总时间：定义了时间步长dt和总时间total_time。

初始核素数量：定义了初始核素数量N0。

衰变常数：定义了每种核素的衰变常数lambda1,lambda2,lambda3。

衰变链矩阵：定义了衰变链矩阵B，其中对角线元素为负值，表示衰变速率。

时间步长函数：定义了欧拉法的时间步长函数euler_step，用于计算下一时间步的核素数量。

活度计算：在时间步长循环中，计算每一步的活度A并更新核素数量N。

输出结果：输出活度随时间的变化。

衰变热计算原理

衰变热是指核素在衰变过程中释放的热量，是评估核燃料和放射性废物热效应的重要参数。衰变热的计算通常基于以下公式：

P

其中：

Pt是在时间t

Nit是在时间t时第i

λi是第i

Ei是第i

N是核素的总数。

衰变热计算的数值方法

衰变热计算同样可以采用数值方法。以下是一个使用龙格-库塔法（四阶）的示例：

importnumpyasnp

#定义时间步长和总时间

dt=0.1#时间步长（秒）

total_time=100#总时间（秒）

#定义初始核素数量

N0=np.array([100,0,0])#初始核素数量[N1,N2,N3]

#定义衰变常数

lambda1=0.1#第一种核素的衰变常数（1/秒）

lambda2=0.2#第二种核素的衰变常数（1/秒）

lambda3=0.3#第三种核素的衰变常数（1/秒）

#定义衰变链矩阵

B=np.array([

[-lambda1,0,0],

[lambda1,-lambda2,0],

[0,lambda2,-lambda3]

])

#定义衰变热能量

E1=1.0#第一种核素释放的能量（焦耳）

E2=2.0#第二种核素释放的能量（焦耳）

E3=3.0#第三种核素释放的能量（焦耳）

#定义龙格-库塔法（四阶）函数

defrunge_kutta_step(N,B,dt):

k1=B@N

k2=B@(N+dt*k1/2)

k3=B@(N+dt*k2/2)

k4=B@(N+dt*k3)

returnN