微分方程模型——偏微分方程

数学建模 专栏收录该内容
4 篇文章 3 订阅

1. 简介

微分方程:描述自然界中存在的物理现象和普遍规律。

  • 常微分方程(ODE)
  • 偏微分方程(PDE)

偏微分方程理论:

  • 物理/工程问题————翻译(建模)/物理工程规律————》数学问题(PDE)
  • 物理/工程问题————求解/数学理论————》数学结果
  • 物理/工程问题————分析————》数学公式/物理意义

偏微分方程的基本概念:

  • 定义:未知函数及其偏导数所满足的方程;F(x, u(x), Du, D2u,…, Dnu) = 0;
  • 阶数:偏微分方程中偏导数的最高阶数,有n阶就为n。
  • 线性偏微分方程:方程中的未知函数或其偏导数的项都是一次项
  • 齐次方程:每项都含有未知函数或未知函数(定义的u = u(x, y))的偏导数。

导出微分方程的基本步骤:

  1. 明确研究的物理量
  2. 明确物理量遵守的物理规律
  3. 按物理规律写出微分方程(泛定方程)
  4. 微元法:划出一个微元,分析临近部分和ta的相互作用

微分方程解决的主要问题:

  • 描述对象特征随时间变化的过程
  • 分析对象特征的变化规律
  • 预报对象特征的未来形态
  • 研究控制对象特征的手段

常用的物理定律概述

  1. 牛顿第二定律:F = ma;
  2. 胡克定律,在弹性限度内,弹性体的张应力和弹性体的形变量成正比:张应力 = 杨氏模量 * 相对伸长,
  3. 热传导的傅立叶定律:在dt时间内,通过面积元dS流入一个小体积元的热量dQ,与沿着面积元法线方向的温度梯度成正比,也与dS和dt成正比:(k是导热系数,由材料决定)
  4. 牛顿冷却定律:单位时间内从周围介质,传到边界上单位面积的热量,与表面和外界的温度差成正比:(这里u1是外界媒质的温度,h为比例系数)
  5. 扩散定律(斐克定律):单位时间流过某横截面的杂质量dm与该横截面积S和浓度梯度成正比:(式中D为扩散系数,负号是指杂质浓度在减少)

2. 栗子——香烟过滤嘴的作用

问题:

  • 过滤嘴的作用和ta的材料与长度有什么关系;
  • 人体吸入的毒品量和哪些因素有关,因素的影响程度怎么样

模型分析:

  • 分析吸烟时毒物进入人体的模型,建立吸烟过程模型;
  • 假设一个机器人持续吸烟,并且环境不变

模型假设(假设的变量都在这儿哦):

  • 烟草长 l1,过滤嘴长 l2,总长度 l = l1 + l2,假设毒品均匀分布;
  • 毒物进入空气和沿香烟穿行的数量比:a’ : a, a’ + a = 1;
  • 未点燃烟草和过滤嘴的吸收率(过滤程度)分别为b和ß;
  • 烟雾沿着香烟的穿行速度为常速v,香烟的燃烧速度为u,v >> u;
  • Q为毒物进入人体的总量
  • q(x, t)为毒物的流量,w(x, t)为毒物密度

模型建立:

  • 建立毒物进入身体的总量:Q = ƒ0T q(l, t)dt, T = l1 / u;

  • 求q(x, 0) = q(x)——流量守恒

    q(x) - q(x + ∆x) = bq(x)∆∂, 0 ≤ x ≤ l1;

    q(x) - q(x + ∆x) = ßq(x)∆∂, l1 ≤ x ≤ l2;

    (∆∂ = ∆x / v)

    同时假设:

    • q(0) = aH0;(香烟毒物被吸入的量)
    • H0 = uw0;(香烟毒物的减少量)

    可得到

    dq/dx = -b/v * q(x), 0 ≤ x ≤ l1;

    dq/dx = -ß/v * q(x), l1 ≤ x ≤ l2;

    &

    q(x) = aH0e-bx/v, 0 ≤ x ≤ l1;

    q(x) = aH0e-bx/ve-ß(x-l1)/v, l1 ≤ x ≤ l2;

  • 目标是求q(l, t),假设t时刻,香烟燃烧到了x = ut;

    将q(x)中的H0拓展为H(t),x随时间会变成x-ut,l1会变成l1 - ut

    可以求得:

    q(x, t) = aH(t)e-b(x-ut)/v,ut ≤ x ≤ l1;

    q(x, t) = aH(t)e-b(l1-ut)/ve-ß(x-l1)/v,l1 ≤ x ≤ l;

  • 求w(ut, t)求∆t内毒物密度的增量。

    w(x, t + ∆t) - w(x, t) = b(q(x, t))/v * ∆t(单位长度烟雾被吸收的部分)

    我们已知:

    • q(x, t) = aH(t)e-b(x-ut)/v;
    • H(t) = uw(ut, t);(H0的拓展公式)

    结果可得到:

    w(ut, t) = w0/a’ * (1 - ae-a’but/v), a’ = 1 - a;

  • 计算Q,吸入一致烟毒物进入人体的总量;

    根据求出来的 w(x, t) 和 q(l, t) 可以代入Q = ƒ0T q(l, t)dt, T = l1 / u;求出现在的Q = aM e-ßl2/v µ®, r = a’bl1/v, µ® = 1 - e-r / r;

结果分析:

  1. Q和a,M成正比,aM是毒物集中在x = l处的吸入量;

  2. e-ßl2/v是过滤嘴的因素;

  3. µ®是烟草的吸收作用;

  4. 判断与不带过滤嘴的香烟比较

    Q1和Q2的差别为b和ß;

    我们已知:

    • 带滤嘴:Q1 = (aw0v / a’b) e-ßl2/v(1 - e-a’bl1/v);
    • 不带滤嘴:Q~12 = (aw0v / a’b) e-bl2/v(1 - e-a’bl1/v);

    所以明显是滤嘴提高了吸收能力

特点:

  • 引入两个基本函数:流量q(x,t)和密度w(x,t),运用物理学的守恒定律建立微分方程,构造动态模型。

3. 偏微方程的导出

3.1. 波动方程的导出

  • 传输线方程(电报方程 )
  • 均匀薄膜的微小横振动方程
  • 流体力学与声学方程
  • 电磁波方程

ta们的物理本质根本不同,但这些方程的数学形式与弦振动方程和杆纵振动方程完全一样:

3.2. 扩散方程的导出

3.2.1. 细杆热传导

现象描述:导热细杆各点的温度是不均匀的,热量由高到低传导。

目的:求出细杆中温度的变化方程。

定律:

  • 热传导的傅立叶定律:q = -k∆u(单位时间内流过单位时间的热量q与温度的梯度成正比,k为热传导系数);
  • 热量守恒定律:热量变化 = 边界流入 + 热源放出;

第一种情况(系统无热源):

(热传导仅由物体内部温度不均引发)

u(x, t)为x处t时刻的温度。

一维情况下热传导的傅立叶定律有: q = -k∂u/∂x;(q为热流强度)

在x方向上(微元法):

  • dt时间流入左表面的热量为:q|x Adt;
  • dt时间流出右表面的热量为:q|x+dx Adt;
  • 所以净流入为:q|x Adt - q|x+dx Adt = -∂q/∂xAdxdt =(然后把q换掉再整合)= kuxxAdxdt;

因为热量Q = c(øAdx)[u|t - u|t+dt] = c(øAdx)utdt;(假设ø为密度)

所以更具热量守恒定律得到kuxxAdxdt = c(øAdx)utdt;

结果:ut - a2uxx = 0, a2 = k/cp(这就得到了均匀物体内部无热源时的热传导方程)

第二种情况(系统内有热源):

比第一问不过就是多个热源,也就是加一个热量呗(设F(x, t)为单位时间,单位体积内产生的热量:kuxxAdxdt + F(x, t)Adxdt = c(øAdx)utdt;

得到:ut - a2uxx = F(x, t) / cp = f(x, t);(ƒ(x, t)为热源强度)

三维情况下就是:ut - a2∆u = f(x, t);

3.2.2. 扩散问题

扩散现象:系统浓度不均用时,物质从高浓度转移向低浓度的现象。

目的:建立空间各点浓度u(x, y, z, t)的方程。

物理规律:扩散定律和粒子数守恒定律

  • 扩散定律(裴克定律):单位时间内流过单位面积的分子数或质量,与浓度的梯度成正比:(D为扩散系数)

    沿着n方向的大小:

  • 粒子数守恒定律:单位时间流入一定体积粒子数流出同一体积粒子数的差,等于该体积单位时间内粒子数的增加量。

    同样是微元法,划出一个小立方体v分析浓度变化规律。

    • 计算单位时间沿x-方向的净流入量(负号是表示和浓度梯度的方向相反):(利用两个平面的流入流出的积分)
    • 同理可以求出y方向和z方向的净流入量:
    • 所以总的净流入量为
    • 单位时间的粒子增加数为:
    • 根据粒子数守恒定律就可以得到:

    我们在前面得到过这样的结果:

    所以我们将扩散定律代入就可以得到三维扩散方程

    我们484看到了吗有很多D,如果扩散是均匀的,D就是一个常数,我们可以令D = a2,则有

    那如果这个空间里面有扩散源,也就是扩散是从这儿来的,那么我们可以将ta们相等得到公式:

3.3. 稳定场方程

3.3.1. 热传导方程

如果边界条件与热源内不随时间变化,一定时间后,内部温度会达到稳态。

u = u(x, y, z), ∂u/∂t = 0;

so

ut - a2∆u = f(x, y, z, t) -> ∆u = -f(x, y, z)泊松方程

3.3.2. 静电场下的泊松方程和Laplace方程

从高斯定理出发,我们可以结合扩散定律简单的得到:

我们可以得到

还因为

我们可以得到两个方程:

  • 泊松方程:
  • Laplace方程(当的时候):

4. PDE导出总结

系统的物理状态除了取决于自己状态,还取决于系统环境的状态。

  • 物理、工程在t>0时刻的系统环境,在数学中称为边界条件;
  • 物理、工程在t=0时刻的系统状态,在数学中称为初始条件;

定解条件:是边界条件和初始条件的总体,反映了问题的个性。

泛定方程:不带边界和初始条件的方程,反映问题的个性。

4.1. 初始条件 ——描述系统的初始状态

  • 波动方程:含有时间的二阶导数,所以需要两个初始条件
    1. u|t=0 = u(x) 系统各点的初始位移
    2. ∂u/∂t|t=0 = v(x) 系统各点的初始速度
  • 热传导方程:含有时间的一阶导数
    1. u(x, t)|t=0 = u(x)初始时刻的温度分布
  • 泊松方程和拉普拉斯方程:不含时间的倒数,不需要初始条件

注意:初始条件给的是初始状态下物理量的分布,而不是指某一位置

4.2. 边界条件

未知函数在边界满足条件:

  1. 第一类Dirchlet边界条件:已知未知函数在边界上的函数值:

    基本形式:u(x, y, z, t)| = ƒ(M, t);

  2. 第二类Neumann:已知未知函数在边界上法线方向的导数值;
    基本形式:∂u(x, t)/∂n|x0 = ƒ(x0, t);

  3. 第三类Robin混合边界条件:混合牛顿冷却定律、傅立叶实验定律(一维);

    牛顿冷却定律:

    q = h(u - ø)n,u为物体表面温度,ø是周围介质温度,h是热交换系数。一维条件下简化为q = h(u - ø)。

    牛顿冷却定律可以得出流出热量与外界温差成正比。


    傅立叶实验定律:

    q|x=a = -k∂u/∂x|x=a

    基本形式:(u + H∂u/∂n)| = ƒ(x0, y0, z0, t);

4.3. 小结

5. 偏微分方程的求解

5.1. 达朗贝尔公式:行波解

基本思想:

  1. 求PDE的通解;
  2. 用定解求特解

关键步骤:通过变量变换,将波动方程化为齐次二阶偏微分方程。(所以也叫做波动方程的初始问题,或者柯西问题)

求解问题:

我们根据式1可以拆分得到,变换变量也就可以转换为

通过复合求导我们可以得到:
——》

一些不要的东西删了:

求通解:

对两个创建的变量进行积分:

积分得到

再积分得到

我们的通解就是

其中ƒ1, ƒ2 是连续可微的一元函数。

ƒ1(x - at)和ƒ2(x + at)的意义?

  • x - at为正方向运动的行波
  • x + at为反方向运动的行波

确定待定函数形式,求特解:

我们已经有初始条件:中的后两条。

我们可以根据之前的条件可以进行两个公式的化简:

  1. image-20210714150847740
  2. image-20210714150921073

所以可以得到我们求出的达朗贝尔公式:
image-20210714151011495

5.2. 分离变量法

基本思想:将偏微分方程通过分离变量变成一个常微分方程!

关键步骤:求特征值问题

适用问题:有界域上的波动方程、热传导方程、稳定场方程等

求解问题:

特征值问题:含有未知常数的常微分方程,求非零解的问题;

特征值:是方程有非零解的常数值;

特征函数:和特征值相对于的非零解

5.3. Fourier变换法(傅立叶变换法)

基本思想:将偏微分方程通过Fourier变换变成一个常微分方程~

关键步骤:求常微分方程定解问题和ta的解的方法,Fourier逆变换

适用问题:无界域上的波动方程、热传导方程等

求解问题:

Fourier变换法基本步骤:

  1. 对偏微分方程与初始条件中实行傅立叶变换,转化为常微分方程;
  2. 解常微分方程的定解问题,得到相应的傅立叶变换式;
  3. 对该式进行逆傅立叶变换,求的原问题的解
  • 1
    点赞
  • 0
    评论
  • 14
    收藏
  • 打赏
    打赏
  • 扫一扫,分享海报

参与评论 您还未登录,请先 登录 后发表或查看评论
©️2022 CSDN 皮肤主题:书香水墨 设计师:CSDN官方博客 返回首页

打赏作者

19-hanhan

你的鼓励将是我创作的最大动力

¥2 ¥4 ¥6 ¥10 ¥20
输入1-500的整数
余额支付 (余额:-- )
扫码支付
扫码支付:¥2
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值