格子演化方法求解非线性PoissonBoltzmann方程.pdf

中国T程热物理学会 传热传质学 2005学术会泌论文 编号：053380 格子演化方法求解非线性 王金库1，王沫然2，李志信1 ‘清华丈学航天航空学院，传热与能源利用北京市重点实验室北京，100084 2D印a岫肌tofMech肌icalEngine嘶lg，me 电话 Em砸l：l地l@蜓D出!自出!翊 界条件的实现方法。本方法不需要对方程进行线性化处理，直接求解原方程，并且通过 数值模拟，验证了该方法的正确性与有效性。 关键词：非线性Pojsson方程，格子演化方法，数值模拟 1 引 言 电动现象(EJe曲_okjnedcpb曲omena)在化学、化工及生物等领域有着广泛的应用，特 别是在单分子DNA或蛋白质检测方面有非常重要的应用前景，近年来是研究的一个热 点。当固壁f或大分子)带电，它会在周围的电解质溶液中诱导出双电层结构，根据双电 层扩散模型，可以得到描述稀溶液内电势P分布的非线性Poisson-B01tmann方程， 1 厂 ．。 、 俨妒一÷∑即‰cXp卜筹y (1) 。60 i ＼“p ／ 很多学者用一维或者二维线性简化或多重网格法求解该方程，多重网格法在求解非 面已经成功应f【{j，该方法也被用来求解其他偏微分方程。chen等人【11用格子气自动机求 解Poisson方程，由于格子气自动机的数值稳定性问题，应用很少。chen等人心提出了 数没有关系，只与位置有关。 且通过数值模拟，验证了其有效性和准确性。 2 数值方法 2．1演化方程 常规二维离散的格子Boltzmann演化方程如方稗(3)所示 1291 (2) 六(r+气巧∥f+4，。)一厶(r，f)=一当[厶(r，，)一口(r，f)] ：兰+pu．vu=一vP+∥v2u (3)、’ 西 詈卸2y+‰(r彤) (4) 其中‰=击≯‰e印(_嚣矿]，代表搠t懒的右飙 g。(r+气4，。，f+4，。)一g名(r，，)=一当[吕。(r，f)一g≯(r，f)] ‘5’ 。+(1一笋4鼎g。 其中平衡分布g?的形式如下， 茎主巍≥： ca， g?=％y，w曲％=f≥乏 4矿鲁 (7) 。=警+os ㈣ y：∑乳+去谚弗。 (9， 通过这样的处理，我们得到了一种格子演化方法求解非线性Poisson方程的方}}一 1292 2+2边界条件 边界条件的实现对于格子演化方法的准确性十分重要‘3”。为了便于说明，以幽1