excおよびmixの計算でini,endをaveと別建てで計算している

[9rnsr]

FlexID/FlexID.Calc/SubRoutine.cs

Lines 97 to 100 in ff6601e

	// 放射能[Bq/day] = 流入元の放射能[Bq/day] * 流入元の生物学的崩壊定数[/day] * 流入割合[-]
	ini += Act.IterNow[inflow.Organ.Index].ini * beforeBio * inflow.Rate;
	ave += Act.IterNow[inflow.Organ.Index].ave * beforeBio * inflow.Rate;
	end += Act.IterNow[inflow.Organ.Index].end * beforeBio * inflow.Rate;

FlexID/FlexID.Calc/SubRoutine.cs

Lines 124 to 127 in ff6601e

	// 放射能[Bq/day] = 流入元の放射能[Bq/day] * 流入元の生物学的崩壊定数[/day] * 流入割合[-]
	ini += Act.IterNow[inflow.Organ.Index].ini * beforeBio * inflow.Rate;
	ave += Act.IterNow[inflow.Organ.Index].ave * beforeBio * inflow.Rate;
	end += Act.IterNow[inflow.Organ.Index].end * beforeBio * inflow.Rate;

FlexIDでは、コンパートメント間の放射能の移動については計算時間メッシュ内で一定であるものとして計算するため、移動元の平均放射能のみが重要なはずである。蓄積コンパートメントはこの例外であり、計算時間メッシュ内での崩壊定数を考慮した減衰を計算するために初期放射能iniと末期放射能endの計算が行われる。

このことを念頭に上に示す実装箇所を見ると、流入する放射能についてiniとendをaveとは別建てで計算しており、何か間違ったことをしているように見える。

[9rnsr]

いや、間違ってはないのか？　例えばn個のacc → mix という経路を考えたとき、mixに流入する放射能のiniとendは、n個のaccのiniとaveに連動するわけだから、現行のiniとendの計算方法自体は間違っているとまでは言い切れない。

mixがexcになってもこれは同じだ。

[9rnsr]

もう少し考えたが、少なくともmixについては、iniとaveはゼロになるべきだと考えられる。以下に自分の考えを示す。

• まず、mixは流入元と流出先を持ち、時間メッシュ内にコンパートメントを通過する放射能はゼロではない。言い換えると、コンパートメントに一時的に存在する放射能を時間メッシュの幅で平均した値=aveは非ゼロになりえる。
• ただし、mixは生物学的崩壊定数が1.0、すなわち流入した放射能を全量流出させるものとして動作する。これはつまり、時間メッシュの末期にコンパートメントに残留する放射能=endはゼロとなるべきであり、すなわち次の時間メッシュの初期放射能=iniもゼロになるべきだと考えられる。
• mixコンパートメントに残留する放射能がどの時点でもゼロならば、これらを積算した放射能=totalもゼロになるべきである。