扩散问题的偏微分模型.ppt_第1页
扩散问题的偏微分模型.ppt_第2页
扩散问题的偏微分模型.ppt_第3页
扩散问题的偏微分模型.ppt_第4页
扩散问题的偏微分模型.ppt_第5页
已阅读5页,还剩9页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

扩散问题的偏微分模型,物质的扩散问题,在石油开采、环境污染、疾病流行、化学反应、新闻传播、煤矿瓦斯爆炸、农田墒情、水利工程、生态问题、房屋基建、神经传导、药物在人体内分布以及超导、液晶、燃烧等诸多自然科学与工程技术领域,十分普遍地存在着。 凡与反映扩散有关的现象,大都能由线性或非线性抛物型偏微分方程作为数学模型来定量或定性地加以解决。,有衰减的扩散问题:设有一扩散源,某物质从此扩散源向四周扩散,沿 x,y,z 三个方向的扩散系数分别为常数,衰减(例如吸收、代谢等)使质量的减少与浓度成正比,扩散前周围空间此物质的浓度为零,估计物质的分布。,设 u(x, y, z, t) 是 t 时刻点 (x, y, z) 处某物质的浓度。任取一个闭曲面 S,它所围的区域是 ,由于扩散,从 t 到 t + t 时刻这段时间内,通过 S流入 的质量为,其中 a2,b2,c2 分别是沿 x,y,z 方向的扩散系数。,由高斯公式,由于衰减, 内的质量减少为,其中 k2 为衰减系数。 由物质不灭定律,在 t 到 t + t 时刻间 内由于扩散与衰减的合作用,积存于 内的质量为M1 M2。,换一个角度看,在 t 到 t + t 时刻间 内由于浓度的变化引起的质量增加为,显然,M3 = M1 M2,即,由 t,t, 的任意性得:,上述方程是常系数线性抛物型方程,它就是有衰减的扩散过程的数学模型。,设扩散源在点 (x0, y0, z0) 处,则此扩散问题满足 Cauchy 问题:,其中 M 为扩散源的质量。用傅立叶变换可求得Cauchy 问题的解析解为,但值得注意的是,在实际应用中,参数 a,b,c,k 往往是很难获得的,通常都是利用观测取样值进行估计,从而得出 u(x, y, z, t) 的近似表达式。 参数估计:目的是对上式中出现的参数 a,b,c,k 进行估计。 已知条件: 点源(扩散源)的质量 M; 点源(扩散源)的位置:(x0, y0, z0); t0 时刻的观测取样值 (xi, yi, zi, mi),mi 为t0 时刻 (xi, yi, zi) 处物质的浓度,i = 1, , n。,首先考虑取样时刻。事实上,取样时刻是未知的,但若设取样时刻为 t0,作变量替换 t = t0,则有 = t/t0,从而,即,上式仍然是常系数线性抛物型方程,与有衰减的扩散过程的数学模型形状完全一致,故可令观测取样值的取样时刻为 t0 = 1。于是,(xi, yi, zi, mi) 满足,其次考虑参数估计。对上式两端取对数,有,令,则有关系式: W = lnu(x, y, z, 1) = X + Y + Z + 由于我们获得的观测取样值 (xi, yi, zi, mi) 可以转化为相应的观测取样值 (Xi, Yi, Zi, Wi),于是利用多元回归分析可以求出 、 的估计值,从而得到参数

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论