扩散问题的偏微分方程模型

扩散问题的偏微分方程模型 第七节 扩散问题的偏微分方程模型 物质的扩散问题，在石油开采、环境污染、疾病流行、化学反应、新闻传播、煤矿瓦斯爆炸、农田墒情、水利 工程 路基工程安全技术交底工程项目施工成本控制工程量增项单年度零星工程技术标正投影法基本原理 、生态问题、房屋基建、神经传导、药物在人体内分布以及超导、液晶、燃烧等诸多自然科学与工程技术领域，十分普遍地存在着. 显然，对这些问题的研究是十分必要的，其中的数学含量极大. 事实上，凡与反应扩散有关的现象，大都能由线性或非线性抛物型偏微分方程作为数学模型来定量或定性地加以解决. ,,,的 试题 中考模拟试题doc幼小衔接 数学试题 下载云南高中历年会考数学试题下载N4真题下载党史题库下载 来自实际，是“真问题数学建模计算机处理”的“三合一”准科研,, 性质的一种竞赛，对上述这种有普遍意义和数学含量高，必须用计算机处理才能得到数值解的扩散问题，当然成为试题的重要来源，例如，,,,,-90,，就是这类试题;,,,,-90,要研究治疗帕金森症的多巴胺(dopamine)在人脑中的分布，此药液注射后在脑子里经历的是扩散衰减过程，可以由线性抛物型方程这一数学模型来刻划. ,,,,-90,要研究单层住宅混凝土地板中的温度变化，也属扩散(热传导)问题，其数学模型与,,,,-90,一样，也是线性抛物型方程. 本文交代扩散问题建模的思路以及如何推导出相应的抛物型方程，如何利用积分变换求解、如何确定方程与解的表达式中的参数等关键数学过程，且以,,,,-90,题为例，显示一个较细致的分析、建模、求解过程. ?1 抛物型方程的导出 设是时刻点处一种物质的浓度. 任取一个闭曲面，它所围的区uxyzt(,,,)(,,)xyztS域是，由于扩散，从到时刻这段时间内，通过流入的质量为 ,,ttt，,S tt，,,,,uuu222. (coscoscos)dSd,,,Mabct,，，1,,,t,,,xyzS 由高斯 公式 小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载 得 222tt，,,,,uuu222()ddddMabcxyzt,，，1,,,,222t . (1) ,,,xyz, 222其中，分别是沿方向的扩散系数. xyz,,abc,, 由于衰减(例如吸收、代谢等)，,内的质量减少为 tt，,2， (2) Mkuxyzt,dddd2,,,,t,2其中是衰减系数. k ,,由物质不灭定律，在内由于扩散与衰减的合作用，积存于内的质量为. MM,12 ,换一种角度看，内由于深度之变化引起的质量增加为 Muxyzttuxyztxyz,，,,[(,,,)(,,,)]ddd3,,,, tt，,,u,dddd. (3)xyzt,,,,t,t, 1 显然，即 MMM,,312 tt，,,uddddxyzt,,,,t,t, 222tt，,,,,uuu2222()dddd.,，，,abckuxyzt222,,,,t,,,xyz, 由之任意性得 ,,tt,, 222,,,,uuuu2222 (4) ,，，,abcku222,,,,txyz 方程(,)是常系数线性抛物型方程，它就是有衰减的扩散过程的数学模型，对于具体问题，尚需与相应的定解条件(初始条件与边界条件等)匹配才能求得确定情况下的解. ?2 ,irac函数 物理学家Dirac为了物理模型之需要，硬是引入了一个当时颇遭微词的，使得数学与物理学传统密切关系出现裂痕的“怪”函数: ，,0,0,x,, (5) ,,() ()1.xxdx,,,,,,,,,0,x, 它的背景是清晰的，以一条无穷长的杆子为例，沿杆建立了一维坐标系，点的坐标为，杆x的线密度是，在段，杆子质量为，则有 ,()x(,],,xmx() xd()mx. (6) ,,,,(), ()d().xxxmx,,,dx 设此无穷长的杆子总质量为,，质量集中在点，则应有 xx,0 1,,xx,,0 或写成 ， mxHxx()(),,mx(),,00,,xx,0, 其中为 Hx() 1,0,x,, Hx(),,0,0,x,, 如果沿用(6)中的算法，则在质量集中分布的这种情形有 0,,xx,,0 ,()x,,,,,0.x, x且， ,()d()xxHxx,,0,,, 于是得 ，,. (7) ,()d1.xx,,,, 但是，从传统数学观点看，若一个函数除某点处处为零，则不论哪种意义下的积分，都必定为零，(7)式岂能成立～但是，函数对于物理学而言是如此之有用，以致物理学家正当地, 拒绝放弃它. 尽管当时数学家们大都嘲笑这种函数，但,.,.,.,irac及其追随者们在物理领域却收获颇丰，Dirac于1933年获诺贝尔物理奖. 当然Dirac也意识到不是一个通常,()x的函数，至于找一种什么办法来阐明这一符号的合法性，那就是数学家的任务了. 1940,()x 2 年，法国数学家许瓦兹(L.Schwartz)严格 证明 住所证明下载场所使用证明下载诊断证明下载住所证明下载爱问住所证明下载爱问 了应用的正确性，把函数置于坚实的,()x, 数学基础上;1950年，L. Schwartz获数学界最高奖Fields奖. 函数的重要性质有: , ，,1). (8) ,()d1xxx,,0,,,，,2). (9) ,()()d()xxfxxfx,,00,,, 其中，即摘出了在的值. fxC()(,),,,，,,()xx,fx()xx,00 dHxx(),0. (10) 3),,,()xx0dx 4)的导数是存在的，不过要到积分号下去理解: ,()x ，,,, (11) ,()()(),xxfxdxfx,,,00,,,，,()()nnn (12) ,()()(1)().xxfxdxfx,,,00,,, 事实上，由于在处为零，则形式地用分部积分公式 ,()xx,，,,,,0，,，,,,,()()()()dxxfxxxfxx,,,00,,,,, ，,, ()()d,,,xxfxx,0,,,n其中，，于是有(11)与(12)公式. fxC()(,),,,，, 5)对于，有 ,()(,)xC,,,，, . (13) ,,,,()()()()xxxxxx,,,000 16),,,,. (14) ()() (0)bxxb||b 7). (15) ,,,,(,,)()()()xxyyzzxxyyzz,,,,,,,000000 8)付立叶变换 ,ix,0 (16) F[()].,yye,,0 (17) F[()]1.,x, (18) FFF[()()][()][()].CxxCxxCxxCxx,,,,,，,,,，,11221122 9)拉普拉斯变换 ,,x0 (19) FF[(),[()1.,,xxex,,,0 (20) FFF[()]()][()[()].CxxCxxCxxCxx,,,,,，,,,，,11221122 从上面的定义与性质看出，Delta函数与一般可微函数还是有重大区别的，我们说,()x 它是“广义函数. ” ?3 Cauchy问题的解 设扩散源在点(,,)xyz处，则此扩散问题满足Cauchy问题 000 222,uuuu,,,,2222abcku, (21),，，,,222 txyz,,,,, ,,,,uxyzMxxyyzz(,,,0)()()(). (22),,,,000, 对(21)(22)进行付立叶变换，且令 3 ˆ， ,,,,,,,(,,), (,)[(,,,)]utuxyztF123 由于 222uuu,,,222ˆˆˆ ,,,[], [], [],uuuFFF,,,,,,123222xyz,,, F[(,,,0)]uxyz ,,,,FFF,,, [()][()][()]Mxxyyzz000 ,，，ixyz,,,()102030, ,Me 故得常微分方程Cauchy问题 du,2222222ˆ()0,,,,，，，，,abcku,123 dt,,，，ixyz,,,()102030,ˆ(0,).uMe,,, 得唯一解 2222222,，，，,，，()()abcktixyz,,,,,,123102030ˆ. (23) utMe(,),, ,1对(23)求逆变换，由于 F2,1,22,21ax,,4a， F[]ee, ,2,,ix,110,e2,212,,a4a， F[]()eexx,,0,故得 ,1ˆuxyztu(,,,)[],F 222 ,,()()()xxyyzz,,,M2000,,,,,expkt,,2223222444atbtct,(2)tabc,, 222,,()()()xxyyzz,,,M2000 (24) exp.kt,,,,,,,222444atbtct,,8tabct,, 如果认为经过了相当长时间后，扩散已经终止，物质分布处于平衡状态，则方程(4) ,u中的，于是有线性椭圆型方程的边值问题 ,0,t 222,,,,uuu22220, (,,)abckuxyzD，，,,,,222,,,xyz , ,,(,,)(,,). uxyzxyz,,D, ,u也可以用付立叶变换求解. 当然，根据实际情况，还可以考虑第二边条件(,,),,xyz,n,D ,uD或第三边条件等，其中是区域的边界，是外法线方向，[](,,),,,n,D，,uxyz,D,n 是实常数. ,,, ?4 参数估计 4 在Cauchy问题(21)(22)的解(23)中，有四个未知的参数，它们分别是扩abck,,,散与衰减过程中的扩散系数与衰减系数的算术平方根. 至于点源的质量与位置 是已知的.设观测取样为: Mxyz,(,,)(,,,), (,,,),,(,,,),xyzmxyzmxyzm00011112222nnnn 22取样时刻为(不然设是取样时间，则(21)变成 ttt,,, UtaUtbU,，，t,100txxyy0022，对而言，取样时间为,,而方程形状与(21)一致)，把在点tcUtkU,(,,)xyz,iii00zz 观测到的物质密度与公式(24)都取对数，令，则 mt,1i 222()()()xxyyzz,,,M2000. (25) ln(,,,1)lnln[]uxyzabck,,,，，，2223444abc,(2) 222()()()xxyyzz,,,111000 令,,,XYZ,,,,,,,,,,,,,,,222444abc M2，则(25)写成 ,,,,lnlnabck3(2), ， (26) WuxyzXYZ,,，，，ln(,,,1),,,,而我们已观测得的数据，用三元回归分析方法求出的(,,,)1,2,,XYZWin,,,,,,,,iiii 估计值如下: ˆˆˆˆ， (27) ,,,,,,，，WXYZ() 其中 nnnn1111WWXXYYZZ,,,, , , , ,,,,,kiiinnnn,,,,1111kkkk ˆˆˆ满足方程组 ,,,,, ˆ,ˆˆ,,,，，,llll,,11121310,,ˆˆˆ ，，,llll,,,,,,21222320 ,ˆˆˆllll,.，，,,,,31323330,, 其中 nn lWXWWlYYWW,,,,,,()(), ()(),,,1020kkkk,,kk11 n lZZWW,,,()(),,30kk,1k nnn222lXXlYYlZZ,,,,,,(), (), (), ,,,112233kkk,,,111kkk nn lXXYYlXXZZ,,,,,,()(), ()(),,,1213kkkk,,kk11 n lYYZZl,,,()(), ,,,lllll, , .,12311332232321kk,k1 111222222ˆˆˆˆˆˆ,,,,,,由可求得的估计值，即. 又由于 abc,,abc, , ,,,,,ˆ,,, 5 M2， (28) kabc,，,,lnln3(2), ˆˆˆˆ由(27)式可得，再把代入(28)得 ,abc,, M2ˆˆˆˆˆ. (29) ,，,,kabclnln3(2), 222222222222ˆˆˆˆ至此得到参数的估计值，把它们代入(24)分别替代，abck,,,abck,,,abck,,, 则得不含未知参数的解的近似表达式. uxyzt(,,,) ?5 竞赛试题分析 AMCM-90A不可用本文的思路与方法加以解决;该试题由东华盛顿大学数学系,ves Nievergelt提供，要求研究药物在脑中的分布，题文称: “研究脑功能失调的人员欲测试新的药物的效果，例如治疗帕金森症往脑部注射多巴胺(Dopamine)的效果，为了精确估计药物影响到的脑部区域，它们必须估计注射后药物在脑内空间分布区域的大小和形状. “研究数据包括50个圆柱体组织样本的每个样本药物含量的测定值(如图6-1)，每个圆柱体长0.76mm，直径0.66mm，这些互相平行的圆柱体样本的中心位于网格距为1mm×0.76×mm×1mm的格点上，所以圆柱体互相间在底面上接触，侧面互不接触. 注射是在最高计数的那个圆柱体的中心附近进行的. 自然在圆柱体之间以及由圆柱体样本的覆盖的区域外也有药物. “试估计受到药物影响的区域中药物的分布. ” -18“一个单位表示一个闪烁微粒的计数，或多巴胺的4.753×10克分子量，例如表6-1指出位于后排当中那个圆柱体的含药量是28353个单位. ” 后方垂直截面 164 442 1320 414 188 480 7022 14411 5158 352 2091 23027 28353 13138 681 789 21260 20921 11731 727 213 1303 3765 1715 453 前方垂直截面 163 324 432 243 166 712 1055 6098 1048 232 2137 15531 19742 4785 330 444 11431 14960 3182 301 294 2061 1036 258 188 图6-1 数学模型只是实际问题的近似，要建立数学模型，一般首先要对所研究的实际问题进行必要和允许的简化与假设，而且，不同的简化与假设，又可能导致不同的数学模型，例如[2]是抛物型方程模型，而[3]则是椭圆方程模型. 假设: (1)注射前大脑中的多巴胺含量可以忽略不计. 6 (2)大脑中多巴胺注射液经历着扩散与衰减的过程，且沿三个方向的扩散系数xyz,, 分别是常数，衰减使质量之减少与深度成正比. (3)注射点在后排中央那个圆柱中心，即注射点的坐标已知，注射量有医(,,)xyz000 疗记录可查，是已知的. (4)注射瞬间完成，可视为点源delta函数. (5)取样也是瞬间完成，取样时间已知为. t,1 (6)样本区域与整个大脑相比可以忽略，样本组织远离脑之边界，不受大脑边界面的影响. 在以上假设之下，显然可以用本文前面讲过的思路来建模，于是得AMCM-90A的数学模型为Cauchy问题(21)(22)，解的表达式为(24)，且用三元回归分析来估出参数，abck,,,于是可以求得任意位置任意时刻药物的深度. 如果所给数据认为是在平衡状态测得的，药物注射进脑后，从高深度处向低深度处扩散，与扩散同时，一部分药物进入脑细胞被吸收固定，扩散系数与吸收系数都是常数，但过一段时间，所有药物都被脑细胞所固定，达到了平衡态. 在这种假设下，[3]给出了下述的分析、建模、求解过程. 设是时刻在点处游离的药物浓度，是时刻vxyzt(,,,)(,,)xyzwxyzt(,,,)(,,)xyztt点处吸收固定的药物浓度，是达到平衡态时点处吸收固定的药物浓度. uxyz(,,)(,,)xyz 又设游离药物在各方向上有相同的扩散系数，吸收系数为，于是有 kh ,v. (30) ,,,kvhv,t ,w又，即吸收速度与游离的浓度成正比，代入(30)得 ,hv,t ,,,vkw. (31) ,,,()w,,,thtt 对(31)关于从,到，,积分得 t k，,，,，,. (32) ,,,vwwttt000,,,h 由于最后无游离药物，故，又开始时无被吸收的药物，故vxyz(,,,)0，,,(0)t, ;平衡状态在时达到，这时 wxyzwxyz(,,,0)0, (,,,0)0,,,uxyz(,,),t,，, ，于是由(32)得 wxyz(,,,)，, k， (33) ,,，,uuvxyz(,,,0)h 其中vxyz(,,,0)是开始时的浓度分布，近似于注射点的点源脉冲函数. 把此注射点取为坐标 k,,2,原点，则是注射量，于是记, (0,0,0)vxyzLxyzL(,,,0)(,,),,,,,h,, 2， (34) ,,，,,,uuLxyz(,,) 作付立叶变换得 2222ˆˆ,,,(),，，，,suuL Lˆu,,2222，，，1()s,,, 再作反变换得 7 1,1222， (35) ,,，，,uLCxyzexp{}222，，xyz 其中是可计算常数. C 222如果考虑各向不同性，设方向上扩散系数分别为，注射点在xyz,,abc,, ，则 (,,)xyz000 222,,uuu,,,222， abcuLxxyyzz,,,()()(),，，，,,,,,,000222xyz,,,,, 于是解为 DL uxyz(,,), 111222()()()xxyyzz,，,，,000222abc ,,111,,222 ，(36) exp1()()(),,，,，,xxyyzz,,000222abc,,,, (36)中的D可计算常数. 用前面类似的方法可以进行参数估计. 在建模过程中，点源函数的使用显然与实况有差别;尤其是认为扩散系数与吸收系数都是常数，对于人脑这种有复杂结构的区域，这种假设与实际不会完全符合;夜间与白天(睡与醒)对这些系数有无影响,脑中各点这些系数是否有变,除时间位置应考虑外，可能还与药液浓度有关. 如此看来，脑内药液分布的数学模型很可能不是常系数线性偏微分方程，而是函数系数的线性微分方程甚至是非线性偏微分方程. 这时，其解不再能用封闭公式来表达，求解过程会变得极为复杂，所以也可以考虑是否试用其他数学模型来解，例如在平衡态的假设下，用回归分析方法建立药液的模拟分布. ufxyz,(,,) 对一个实际问题，其数学模型未必唯一，各模型间孰优孰劣，没有一般的判别法，须经实践来检验. 参 考 文 献 [1]叶其孝，大学生数学建模竞赛辅导教材，湖南教育出版社，1993. [2]Christopher, R. Malone, Gian Pauletto, James, I. Zoellick, Distribution of Dopamine in the Brain, The Journal of Under graduate Mathematics, and its Applications, vol. 12(1991), Special Issue: The 1991 Mathematical Contest in Modeling, pp. 211-223. [3]孙晓东，荆秦，梁俊，脑中药物分布的数学模型，数学的实践与认识，1991年No. 4，63-69. [4]中国科学院数理统计组，常用数理统计方法，科学出版社，19784. 8