扩散问题地偏微分方程模型,数学建模.docxVIP

标准实用 第七节 扩散问题的偏微分方程模型 物质的扩散问题，在石油开采、环境污染、疾病流行、化学反应、新闻传播、煤矿瓦斯 爆炸、 农田墒情、 水利工程、 生态问题、 房屋基建、 神经传导、 药物在人体内分布以及超导、 液晶、 燃烧等诸多自然科学与工程技术领域， 十分普遍地存在着 . 显然， 对这些问题的研究 是十分必要的， 其中的数学含量极大 . 事实上， 凡与反应扩散有关的现象， 大都能由线性或 非线性抛物型偏微分方程作为数学模型来定量或定性地加以解决 . ＭＣＭ的试题来自实际，是“真问题 数学建模 计算机处理”的“三合一”准科研 性质的一种竞赛， 对上述这种有普遍意义和数学含量高， 必须用计算机处理才能得到数值解 的扩散问题，当然成为试题的重要来源，例如，ＡＭＣＭ -90 Ａ，就是这类试题；ＡＭＣＭ -90 Ａ要研究治疗帕金森症的多巴胺（ dopamine ）在人脑中的分布，此药液注射后在脑子里 经历的是扩散衰减过程，可以由线性抛物型方程这一数学模型来刻划 . ＡＭＣＭ -90 Ａ要研 究单层住宅混凝土地板中的温度变化， 也属扩散（热传导） 问题，其数学模型与ＡＭＣＭ -90 Ａ一样，也是线性抛物型方程 . 本文交代扩散问题建模的思路以及如何推导出相应的抛物型方程， 如何利用积分变换求 解、如何确定方程与解的表达式中的参数等关键数学过程，且以ＡＭＣＭ -90 Ａ题为例，显 示一个较细致的分析、建模、求解过程 . §1 抛物型方程的导出 设 u(x, y, z, t) 是 t 时刻点 ( x, y, z) 处一种物质的浓度 . 任取一个闭曲面 S ，它所围的区 域是 ，由于扩散，从 t 到 t t 时刻这段时间内，通过 S 流入 的质量为 M 1 t t u cos b2 u cos c2 u cos )dSdt . ( a2 t S x y z 由高斯公式得 M 1 t t 2 2 u b 2 2u c 2 2u ( a x 2 y 2 z 2 )d xdydzdt （ 1） t . 其中， a2 , b2 , c2 分别是沿 x, y, z 方向的扩散系数 . 由于衰减（例如吸收、代谢等） ， 内的质量减少为 M 2 t t k2udxdydzdt ， （ 2） t 其中 k 2 是衰减系数 . 内的质量为 M1 M2. 由物质不灭定律，在 内由于扩散与衰减的合作用，积存于 换一种角度看， 内由于深度之变化引起的质量增加为 M 3 [ u(x, y, z, t t) u( x, y, z, t)]dxdydz t t t  u dxdydzdt. (3) 显然 M3 M1 M2，即 文案大全 标准实用 t t udxdydzdt t t t t (a 2 2u b 2 2u c 2 2u 2 u)dxdydzdt. t 2 2 2 k x y z 由 t, t, 之任意性得 u a2 2 u b2 2u c2 2 u k2u （ 4） t x2 y2 z2 方程（４）是常系数线性抛物型方程， 它就是有衰减的扩散过程的数学模型， 对于具体问题， 尚需与相应的定解条件（初始条件与边界条件等）匹配才能求得确定情况下的解 . 2 Ｄirac 函数 物理学家 Dirac 为了物理模型之需要， 硬是引入了一个当时颇遭微词的， 使得数学与物理学传统密切关系出现裂痕的“怪”函数： ( x) 0, x 0, 1. （ 5） , x (x)dx 0, 它的背景是清晰的，以一条无穷长的杆子为例，沿杆建立了一维坐标系，点的坐标为 x ，杆 的线密度是 ( x) ，在 ( , x] 段，杆子质量为 m( x) ，则有 dm( x) x dx ( x), ( x)dx m( x). . （ 6） 设此无穷长的杆子总质量为１，质量集中在x x0 点，则应有 m( x) 1,x x0 , H (x x0 ) ， 0, x 或写成 m(x) x0, 其中 H (x) 为 1,x 0, H (x) 0, x 0, 如果沿用（ 6）中的算法，则在质量集中分布的这种情形有 0, x x0 , ( x) 0. , x x 且( x)d x H (x x0 ) ， 于是得 (x)dx 1. . （ 7） 但是，从传统数学观点看，若一个函数除某点处处为零，则不论哪种意义下的积分，都必定 为零，（ 7）式岂能成立！但是， 函数对于物理学而言是如此之有用，以致物理学家正当地 拒绝放弃它 . 尽管当时数学家们大都嘲笑这种函数， 但Ｐ . Ａ . Ｍ . Ｄ irac 及其追随者们在物 理领域却收获颇丰， Dirac 于 1933 年获诺贝尔物理奖 . 当然 Dirac 也意识到 ( x) 不是一个 通常的函数，