两种方式推导人口平衡方程(PBE)

前言

之前玩了玩雷诺运输定理(Reynolds transport theorem)，当时用的是人口平衡方程(Population Balance Equation)，后来一想不对啊，这引入没头没尾的。其实PBE在各个领域上都有运用，最最原始的版本是反映某一时期人口数量及其自然变动，到了后来就运用到物化生各个学科上面，描述的是物质的量能变化情况。现在主要落地的模型是经济学那一块，和机器学习倒是关联不大。不过现在MLDL的黑箱情况那么严重，以后发展方向会不会重新回到模型依靠型也不好说。
总之嘛，PBE也算是常见常用的偏微分模型了，一维形式的推导还是值得掌握一下的，就是我的偏微分忘得差不多了，推出差池是可以预见的……

化学推导

考虑物质M分布在截面积为A的长管中，对于∀区间[xa,Xb],物质M的浓度变化为：
M总量变化率=M流入区间速率-M流出区间速率
再考虑区间中发生的化学反应（M的生成和消耗），有：
M浓度变化率=M流入区间速率-M流出区间速率+M生成速率-M消耗速率
定义：
1、c(𝑥,𝑡):物质M在𝑥处𝑡时刻的单位面积浓度
2、J(𝑥,𝑡):物质M在边界𝑥处𝑡时刻从左向右的净流动线速度
3、𝑓(𝑥,𝑡,𝑐): 物质M在𝑥处𝑡时刻浓度为𝑐时，单位体积单位时间的净产生量，如上定义𝑐为关于𝑥,𝑡的函数
那么在如图所示𝑓的小区间中：
物质M的总量为：${\int }_{x_b}^{x_a}c(x,t)Adx$
物质M的净通过量：$AJ\left(x_a,t\right)-AJ\left(x_b,t\right)$
物质M净产生率为：${\int }_{x_b}^{x_a}f(x,t,c)Adx$

联立三式写出区间内M浓度变化率，然后化简一下，大致像这样
$$\begin{gathered} \frac{d}{dt}{\int }_{x_b}^{x_a}c(x,t)Adx=({\int }_{x_b}^{x_a}f(x,t,c(x,t))Adx+A\left(J\left(x_a,t\right)-J\left(x_b,t\right)\right) \\ \Rightarrow \frac{d}{dt}{\int }_{} \end{gathered}$$