方程变分形式的妙处在于很容易得到物理上的守恒关系,守恒关系包括全局守恒和局部守恒。前者表示在整个求解区域上产生的粒子与消失的粒子相等,后者表示在一个单元上产生和消失的粒子相等。守恒关系为:消失=产生,消失包括:泄露,反应,产生由源项产生。以下对输运方程的各种变分形式的守恒形式进行说明
CG 方程
中子输运方程连续有限元的变分形式为
$$<φ,Ω⋅nϕ> - (Ω⋅∇φ,ϕ)+(φ,σϕ) = (φ,s)$$
令检验函数 $φ$在整个区域中恒为 1,那么
$$<1,ω⋅nϕ> +(1,σϕ) = (1,s)$$1,ω⋅nϕ>
该式第一项为泄漏率,第二项为反应消失率,右端项为产生率
SAAF 方程
SAAF方程的变分形式为
$$<φ,Ω⋅nϕ> + (Ω⋅∇φ, \frac{1}{σ}Ω⋅∇ϕ)+(φ,σϕ) = (Ω⋅∇φ,\frac{1}{σ}s)+(φ,s)$$
令检验函数 $φ$在整个区域中恒为 1,那么
$$<1,ω⋅nϕ> +(1,σϕ) = (1,s)$$1,ω⋅nϕ>
满足守恒关系
LS 方程
LS 方程的变分形式为
$$(Ω⋅∇φ,Ω⋅∇ϕ) + (σφ,Ω⋅∇ϕ) + (Ω⋅∇φ,σϕ) + (σφ,σϕ) -
令检验函数 $φ$在整个区域中恒为 1,那么
$$<σ,Ω⋅nϕ> + (σ,σϕ) -
如果在整个区域上 $σ$ 为大于零的常数,并且 $c=σ$;方程才是全局守恒的,否则不是。
- 对于真空区域
$$(Ω⋅∇φ,Ω⋅∇ϕ) -^{-} = 0$$
所以
$$(Ω⋅∇φ,Ω⋅∇ϕ)+ (cφ,Ω⋅∇ϕ) + (Ω⋅∇φ,cϕ)-^{+} + ^{-} = 0$$
令检验函数 $φ$在整个区域中恒为 1,那么
$$ (c,Ω⋅∇ϕ)-^{+} + ^{-} = 0$$
只有当 $Ω⋅∇ϕ $ 恒等于零时,才满足守恒关系。即在真空区域要恒满足输运方程 $Ω⋅∇ϕ=0$,数值误差为0。
SUPG 方程
SUPG 方程的变分形式为
$$<φ,Ω⋅nϕ> -(Ω⋅∇φ,ϕ)+(φ,σϕ) + (τΩ⋅∇φ,Ω⋅∇ϕ)+(τΩ⋅∇φ,σϕ) = (φ,s)+ (τΩ⋅∇φ,s)$$
令检验函数 $φ$在整个区域中恒为 1,那么
$$<1,ω⋅nϕ> +(1,σϕ) = (1,s)$$1,ω⋅nϕ>
满足守恒关系
SGS 方程
SGS 方程的变分形式为
$$<φ_c,Ω⋅nϕ> - (Ω⋅∇φ_c,ϕ)+(φ_c,σϕ) = (φ_c,s)$$
$$<φ_d,Ω⋅nϕ_d> - (Ω⋅∇φ_d,ϕ_d)+(φ_d,Ω⋅∇ϕ_c)+(φ_d,σϕ) = (φ_d,s)$$
令检验函数 $φ_c$在整个区域中恒为 1,那么
$$<1,ω⋅nϕ> +(1,σϕ) = (1,s)$$1,ω⋅nϕ>
满足全局守恒关系
令检验函数 $φ_d$在单元中恒为 1,那么
$$<1,ω⋅nϕ_d> +(1,Ω⋅∇ϕ_c)+(1,σϕ) = (1,s)$$1,ω⋅nϕ_d>
$$<1,ω⋅nϕ> +(1,σϕ) = (1,s)$$
满足局部守恒关系1,ω⋅nϕ>