偏微分方程课件
合集下载
第5章偏微分方程值解ppt课件

t
t nt , x ix , y jy , z kz
总目录
本章目录
5.1
5.2
5.3
5.4
5.2 基本离散化公式
以3对于二阶偏导,我们可以通过对泰勒展开式处 理技术得到下面离散化计算公式:
2u t 2 2u x 2 2u y 2 2u z 2
总目录
本章目录
5.1
5.2
5.3
5.4
5.3 几种常见偏微分方程的离散化计算
例下面介绍3种迭代格式: 1 u (u u u u (1)同步迭代: 4 1 u (u u u u (2)异步迭代: 4 1 u u u ) u (u 4 (3)超松弛迭代:
(5-4) (计算实例VB程序见课本)
总目录
本章目录
5.1
5.2
5.3
5.4
5.3 几种常见偏微分方程的离散化计算
2、一维流动传热传导方程的混合问题 一维流动传热传导方程的混合问题:
2 u u 2 u b f (u, t ) a 2 t x x u t 0 (x), u 0 x x l u x 0 μ1(t)
u
x0
1 (t ),u xt 2 (t )
为初值条件 为边值条件
当该波动方程只提初值条件时,称此方程为波动 方程的初值问题,二者均提时,称为波动方程的 混合问题。
总目录 本章目录
5.1
5.2
5.3
5.4
5.3 几种常见偏微分方程的离散化计算
t t
x
0
x
0
l
(a)初值问题
计算流体力学基础_P2_偏微分方程的性质 ppt课件

di( a 1 , g 2,3)
对于左边界:
条件
描述
u0 anduc u0 anduc
u0 anduc
超音速入口 亚音速入口 超音速出口
u0 anduc 亚音速出口
边界条件设定
给定3个边界条件 给定2个边界条件 无需给定边界条件 给定1个边界条件
知识点
Slide 14
5. 椭圆型方程:Laplace方程
Uu
E
0
1
0
AU f ((232)u3)u2u/2c21
(3)u c2 32u2 1 2
1
u
推导
u f(U)u2 p
u(Ep)
守恒变量:质量 密度、动量密度、 能量密度
u1 U u u2
E u3
u 1,uu 2/u 1,Eu 3
E p 1 u2 1 2
p
c/a 0
I A 0 a 2 b c 0 ( 3 )
特征方程(3)有两个互异实根 -> 矩阵A可对角化 -> 双曲型
特征方程(3) 有两个相同实根,且无法对角化 -> 抛物型
特征方程(3)无实根
-> 椭圆型
Slide 11
4. 讨论Euler方程组
一维非定常流动:
f(U)AU
x
x
U f(U) 0 t x
则有:
duaubuc ds x y
特征相容关系 (特征线上物理量的简化方程)
✓偏微方程在特征线上变成了常微分方程 Slide 7
演示: 如何利用特征线计算物理量
a(x,y)ub(x,y)uc(x,y)
x
y
特征线法是空气动力学重要的计算方 法。早期(计算机出现之前),是主 y 要的CFD手工计算方法之一。
一阶偏微分方程教程省公开课获奖课件说课比赛一等奖课件

z x ydx 得 x 2 z x y C2 ,所以得到另一种首次积分为
x2 zxy
于是原方程旳隐式通解为
2x z C1
x 2 z x y C2 7
由(3)可得
dx dy dz P(x, y, z) Q(x, y, z) R(x, y, z)
于是得到方程组(3)旳一种等价形式:
精神病药物研究需测定新药旳效果,例如治疗帕 金森症旳多巴胺旳脑部注射效果。为了精确估计药 物影响旳脑部区域,我们必须估计注射后药物在空 间旳分布形状和尺寸。
研究旳数据涉及50根圆柱组织样本中每一根所含 药物旳测量值(见表1、表2及图1)。每一圆柱旳长度 为0.76mm,直径为0.66mm。这些平行圆柱旳中心 位于1mm×0.76mm×1mm旳网格点上。所以,圆
首次积分为 y u, 2x u2
于是原方程旳隐式通解为 y u, 2x u2 0
其中 为任意二元连续可微函数。
将该解代入初始条件,得 y, 2 y 0
于是有 2x u2 2( y u) ,解得 u 1 1 2(x y)
再由初始条件得Cauchy问题旳解为
u 1 1 2(x y)
25
p(a, t ) a
p(a,t) t
(a,t,
N (t))
p(a,t)
f
(a,t),
a 0, t 0
p(a,
0)
p0 (a),
a0
(4)
p(0,
t)
(a,t, N (t)) p(a,t)da,
0
t 0
N (t) 0 p(a,t)da, t 0
26
精神病用药问题旳方程模型
• 问题旳提出
表2 前方垂直截面
163 324 432 243 166 712 4055 6098 1048 232 2137 15531 19742 4785 330 444 11431 14960 3182 301 294 2061 1036 258 188
偏微分方程课件-Ch1

《偏微分方程》第一章 绪论
第一章 绪论
1.1
《偏微分方程》第一章 绪论
在偏微分方程中最高阶偏微商的阶数叫做该偏微分方 程的阶。
在偏微分方程组中最高阶偏微商的阶数叫做该偏微分 方程组的阶。
如果一个函数在其自变量 (x1, x2 , , xn )的某变化范 围内连续并且具有方程(方程组)的一切连续偏微商 将它代入方程后使其成为恒等式,则称该函数是方程
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
有三种方程的划分(1.3.16)可知,在椭圆型区域内不存 在实的特征方向;在双曲型区域内存在两族实的特征方 向;而在抛物型的点上仅有一个实的特征方向。因此, 方程双曲型区域被两族实特征曲线网覆盖,在抛物型的 点集被一簇实特征线网覆盖
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
考虑方程
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
考虑(1.3.8)的线性主部
《偏微分方程》第一章 绪论
支覆盖了整个下半平面4,并且每个分支都与 x 轴相切。
(3)在椭圆型区域 y 0 中特征方程的其中的一个复解为
x 2i y c
去实部和虚部作变换
x 2 y
《偏微分方程》第一章 绪论
经过计算便得到方程在上半平面的标准型
u
u
1
u
0,
y
0.
《偏微分方程》第一章 绪论
第一章 绪论
1.1
《偏微分方程》第一章 绪论
在偏微分方程中最高阶偏微商的阶数叫做该偏微分方 程的阶。
在偏微分方程组中最高阶偏微商的阶数叫做该偏微分 方程组的阶。
如果一个函数在其自变量 (x1, x2 , , xn )的某变化范 围内连续并且具有方程(方程组)的一切连续偏微商 将它代入方程后使其成为恒等式,则称该函数是方程
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
有三种方程的划分(1.3.16)可知,在椭圆型区域内不存 在实的特征方向;在双曲型区域内存在两族实的特征方 向;而在抛物型的点上仅有一个实的特征方向。因此, 方程双曲型区域被两族实特征曲线网覆盖,在抛物型的 点集被一簇实特征线网覆盖
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
考虑方程
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
考虑(1.3.8)的线性主部
《偏微分方程》第一章 绪论
支覆盖了整个下半平面4,并且每个分支都与 x 轴相切。
(3)在椭圆型区域 y 0 中特征方程的其中的一个复解为
x 2i y c
去实部和虚部作变换
x 2 y
《偏微分方程》第一章 绪论
经过计算便得到方程在上半平面的标准型
u
u
1
u
0,
y
0.
《偏微分方程》第一章 绪论
偏微分方程及其求解实例ppt课件

(hn1-2.*h(k,n)+h(k,n-1))./dr.^2);
end plot(r(3:n)./ra,p(k,3:n).*theta.*2./rb)
h hi1 hi1 r i 2r
2h hi1 2hi hi1
r 2
r 2
i
P
1 rb 4
1
r
h r
2h r 2
偏微分方程的求解实例2:
2u A x2
2u B
xy
C
2u y 2
D u x
E u y
Fu
f
x,
y,u,
u x
,
u y
(1) 导热方程:
u 2u
t x2 (2) 拉普拉斯方程: 如稳态静电场和稳态温度分布模型
2u 2u 0
x2 y2
(3) 波动方程: 一维弦振动模型
2u 2 2u
t 2
x2
偏微分方程的边界条件
function PDE1Dd_CrankNicolson % 使用Crank-Nicolson有限差分方法求解一维动态传
热模型
c1 = 100; c2 = 0; a = 10; b = 8; alpha = 2; n = 6; m = 8; U = CrankNicolson(@ic,c1,c2,a,b,alpha,n,m)
h t 3 9c
9c
h3 h33
4h r 4
3
h5 4h4
6h3 4h2 r 4
h1
h t
n
V
r i 2r
2h hi1 2hi hi1
r 2
r 2
i
3h r 3
hi2
2hi1 2hi1 2r 3
计算机应用基础偏微分方程求解PPT课件

6.2 二阶偏微分方程的求解
二 抛物线型偏微分方程
第16页/共43页
6.2 二阶偏微分方程的求解
parabolic函数用于求解抛物型偏微分方程的解,调用格 式如下:
u1=parabolic(u0,tlist,b,p,e,t,c,a,f,d) b: 边界条件 u0: 初始条件 tlist;时间列表 u1:对应于tlist的解向量 p,e,t :网格数据
• 启动偏微分方程求解界面
– 在 MATLAB 下键入 pdetool
• 该界面分为四个部分
– 菜单系统 – 工具栏 – 集合编辑 – 求解区域
第20页/共43页
6.3 偏微分方程求解工具箱
菜单栏
工具栏
第21页/共43页
6.3 偏微分方程求解工具箱
第22页/共43页
5.3 偏微分方程求解工具箱
第9页/共43页
6.1 偏微分方程组求解
边界条件程序”c7mbc.m” function [pa, qa, pb, qb]=c7mpbc(xa, ua, xb, ub, t) pa=[0; ua(2)]; qa=[1; 0]; pb=[ub(1)-1; 0]; qb=[0; 1];
function u0=c7mpic(x) u0=[1; 0];
进入反应器,相当于总质量速率为G=2500kg.h-1.m2。反应管
外用速率为F 130kg h-1烟道气与反应混合物
逆流加热反应管,烟道气出口温度为620 C。其
它数据:催化剂的堆积密度=1440kg / m3,操作
压力P 1.2bar,乙苯的反应热H=140000kJ / m ol,
床层有效导热系数e 0.45w.m1.k 1,有效扩散系数
偏微分方程初步介绍公开课获奖课件百校联赛一等奖课件

0, 0
第三边值问题(Robin)
经典旳定解问题举例
热传导方程旳初、边值问题
u t
a2
2u x 2
f (x, t),
t 0,0 x L
u(x, t) (x)
t 0
u( x, t) x0
g (t), u(x, t) xL
h(t)
何为适定性?
存在性 唯一性 连续依赖性(稳定性)
自变量 未知函数
F (x, u,
u x1
,,
u xn
,
2u x12
,)
0
偏微分方程旳一般形式
某些概念
PDE旳阶 古典解
PDE旳 解
广义解
线性PDE
非线性PDE
是指这么一种函数,它本身以及它旳偏导 数在所考虑旳区域上连续,同步用满足方 程。
半线性PDE 拟线性PDE 完全非线性PDE
线性PDE: PDE中对所含未知函数及其各阶导数旳全体都是线 性旳。 线性PDE中全部具同一最高阶数旳偏导数构成旳 部分,称为线性方程旳主部。
r x2 y2
6.
u t
6u
u x
3u x3
0
KDV方程
特解都不易找到
7. ut uux eu
拟线性PDE
8.
v x v xx
v
2 y
v
yy
v2
拟线性PDE
9. a( x, y)(vxx vyy ) ev (vx vy ) 半线性PDE
10. ut ux sin u
11. ut 2 ux 2 u 2
a22
y
y
a11 x
a22
y
a11 x
a22
《偏微分方程的建立》课件

根据问题的性质和已知条件,通过数学推导和逻辑推理建立偏微分方程。
求解方程
根据建立的偏微分方程,选择适当的数学方法和计算工具进行求解。
验证解的正确性
通过对比实际数据或实验结果,验证所求偏微分方程的解是否符合实际情况。
描述物理现象的偏微分方程在力学、电磁学、光学等领域有广泛应用。
物理学
偏微分方程在金融领域的应用主要涉及资产定价和风险管理等方面。
工程问题
描述生理过程、药物动力学等。
生物医学问题
04
CHAPTER
偏微分方程的实际应用案例
总结词
描述人口随时间变化的规律
详细描述
通过偏微分方程建立人口动态模型,考虑出生率、死亡率以及迁移率等因素对人口数量的影响,预测未来人口数量变化趋势。
总结词
模拟热量在物体中的传递过程
详细描述
利用偏微分方程建立热传导模型,描述热量在材料中的扩散过程,常用于材料科学、能源工程等领域。
金融学
在机械工程、航空航天、电子工程等领域,偏微分方程被用来描述各种物理现象和工程问题。
工程学
在生态学、生理学和流行病学等领域,偏微分方程被用来描述种群动态和疾病传播等现象。
生物学Biblioteka 010302
04
03
CHAPTER
偏微分方程的求解方法
方程中的未知函数及其导数都是一次幂或常数。
线性偏微分方程
方程中的未知函数及其导数是二次幂或更高次幂。
THANKS
感谢您的观看。
VS
模拟流体运动规律和特性
详细描述
利用偏微分方程建立流体动力学模型,研究流体运动的速度场、压力场、温度场等特性,广泛应用于航空航天、船舶、能源等领域。
总结词
求解方程
根据建立的偏微分方程,选择适当的数学方法和计算工具进行求解。
验证解的正确性
通过对比实际数据或实验结果,验证所求偏微分方程的解是否符合实际情况。
描述物理现象的偏微分方程在力学、电磁学、光学等领域有广泛应用。
物理学
偏微分方程在金融领域的应用主要涉及资产定价和风险管理等方面。
工程问题
描述生理过程、药物动力学等。
生物医学问题
04
CHAPTER
偏微分方程的实际应用案例
总结词
描述人口随时间变化的规律
详细描述
通过偏微分方程建立人口动态模型,考虑出生率、死亡率以及迁移率等因素对人口数量的影响,预测未来人口数量变化趋势。
总结词
模拟热量在物体中的传递过程
详细描述
利用偏微分方程建立热传导模型,描述热量在材料中的扩散过程,常用于材料科学、能源工程等领域。
金融学
在机械工程、航空航天、电子工程等领域,偏微分方程被用来描述各种物理现象和工程问题。
工程学
在生态学、生理学和流行病学等领域,偏微分方程被用来描述种群动态和疾病传播等现象。
生物学Biblioteka 010302
04
03
CHAPTER
偏微分方程的求解方法
方程中的未知函数及其导数都是一次幂或常数。
线性偏微分方程
方程中的未知函数及其导数是二次幂或更高次幂。
THANKS
感谢您的观看。
VS
模拟流体运动规律和特性
详细描述
利用偏微分方程建立流体动力学模型,研究流体运动的速度场、压力场、温度场等特性,广泛应用于航空航天、船舶、能源等领域。
总结词
偏微分方程数值解法(抛物型方程差分法)2省名师优质课赛课获奖课件市赛课一等奖课件

j(H
sin 2
)
| j
|
1 jn
2(n 1)
| (ra2 ) |
(ra2 )
| n |
| 1 |
ra 2
极值点满足
ra2 1 2
1 4ra2 sin2 4ra2 sin2 n 1
2(n 1)
2(n 1)
(ra2 ) 1 2sin2
cos 1
2(n 1)
n1
显式差分格式稳定充分条件. h2 / 2a 2
4/17
无穷大范数定义 ||
uk
||
max
1 jn
|
ukj
|
双层差分格式
n
n
u (k ) k1 jm m
u (k ) k
jm m
f
k j
m1
m1
记矩阵
A(k )
(
) (k )
jm nn
B(k)
(
) (k )
jm nn
双层格式旳矩阵形式 A(k )uk1 B(k )uk f k
双层差分格式初值稳定概念:
2ra2 cos j )]1
j
n
1
]1
[1 4ra2 sin2 nj 1 ]1 1
2(n 1)
11/17
过渡矩阵旳谱半径
(
H
)
max
1 j N 1
|
j
(
H
)
|
1
max
1 jn
|
1
4ra
sin2 (
j
/
2(n
1))
|
1
1 4ra sin2( / 2(n 1)) 1
偏微分方程课件 云南财经大学

1.1.5. 非线性偏微分方程 我们把不是线性偏微分方程的偏微分方程统称为非线性偏 微分方程。在非线性偏微分方程中, 如果关于未知函数的所有 最高阶偏导数都是线性的, 则称它为拟线性偏微分方程。
二阶拟线性偏微分方程 二阶拟线性偏微分方程 三阶拟线性偏微分方程
在拟线性偏微分方程中, 由最高阶偏导数所组成的那一部 分, 称为方程的主部; 若主部内的系数都是常数或是自变量的 已知函数, 这时方程被称为是半线性的。
如果给定一个函数 u (x) , 将它及它对自变量的各阶偏导
数代入方程(1.1.1), 能使(1.1.1)成为恒等式, 则称函数是偏微分方 程(1.1.1)的解。
我们知道, 一个常微分方程如果有解, 就必有无穷多个解, 其表现形式是依赖于一个或几个任意常数的通解. 于是自然会 想到偏微分方程的通解也会含有任意元素.
它被称为三维Laplace方程。
利用Laplace算子
2 x2
2 y2
2 z2
,三维Laplace方程写成
u 0
对于函数 u u(x1, x2, , xn ,t) 的n维Laplace方程,利用
Laplace算子
2 x12
2 x22
2 xn2
则偏微分方程的一般形式为
实自变量 未知函数
5
机动 目录 上页 下页 返回 结束
《偏微分方程》第一章 绪论 第6页
其中是F自变量x,未知函数u及u的有限多个偏导数的已知函数. 例如关系式
等都是偏微分方程.
6
机动 目录 上页 下页 返回 结束
《偏微分方程》第一章 绪论 第7页
1.1.2. 偏微分方程的解
m
二阶拟线性偏微分方程 二阶拟线性偏微分方程 三阶拟线性偏微分方程
在拟线性偏微分方程中, 由最高阶偏导数所组成的那一部 分, 称为方程的主部; 若主部内的系数都是常数或是自变量的 已知函数, 这时方程被称为是半线性的。
如果给定一个函数 u (x) , 将它及它对自变量的各阶偏导
数代入方程(1.1.1), 能使(1.1.1)成为恒等式, 则称函数是偏微分方 程(1.1.1)的解。
我们知道, 一个常微分方程如果有解, 就必有无穷多个解, 其表现形式是依赖于一个或几个任意常数的通解. 于是自然会 想到偏微分方程的通解也会含有任意元素.
它被称为三维Laplace方程。
利用Laplace算子
2 x2
2 y2
2 z2
,三维Laplace方程写成
u 0
对于函数 u u(x1, x2, , xn ,t) 的n维Laplace方程,利用
Laplace算子
2 x12
2 x22
2 xn2
则偏微分方程的一般形式为
实自变量 未知函数
5
机动 目录 上页 下页 返回 结束
《偏微分方程》第一章 绪论 第6页
其中是F自变量x,未知函数u及u的有限多个偏导数的已知函数. 例如关系式
等都是偏微分方程.
6
机动 目录 上页 下页 返回 结束
《偏微分方程》第一章 绪论 第7页
1.1.2. 偏微分方程的解
m
《偏微分方程Ch》PPT课件

《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
分析可得上述初值问题的形式解是:
u(x,t) 1 [(x at) (x at)] 1
xat
( y)dy
2
2a xat
称此式为d’Alembert(达朗贝尔)公式
《偏微分方程》第3章 波动方程
则有达朗贝尔公式得:
1
1 xt
u(x,t) [g(x t) g(x t)]
h( y)dy
2
2a xt
从而,有 u, g, h 的定义,便得到原问题的解
u(
x,
t)
1 2 1 2
[g( [g(
x x
t) t)
g g
(x t)] (t x)]
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
的方法.它是把求解非齐次方程的问题归结为 解一个齐次方程的问题,是常微分方程中的常 数变易法在线性偏微分方程中的推广.通常称 这个方法为Duhamel 原理,又叫齐次化原理, 或从物理的角度称为 冲量原理.
《偏微分方程》第3章 波动方程
• 下面,我们以求解三维非齐次波动方程的初值问题 为例,说明这个方法的思想.考虑问题
x
x
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
分析可得上述初值问题的形式解是:
u(x,t) 1 [(x at) (x at)] 1
xat
( y)dy
2
2a xat
称此式为d’Alembert(达朗贝尔)公式
《偏微分方程》第3章 波动方程
则有达朗贝尔公式得:
1
1 xt
u(x,t) [g(x t) g(x t)]
h( y)dy
2
2a xt
从而,有 u, g, h 的定义,便得到原问题的解
u(
x,
t)
1 2 1 2
[g( [g(
x x
t) t)
g g
(x t)] (t x)]
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
《偏微分方程》第3章 波动方程
的方法.它是把求解非齐次方程的问题归结为 解一个齐次方程的问题,是常微分方程中的常 数变易法在线性偏微分方程中的推广.通常称 这个方法为Duhamel 原理,又叫齐次化原理, 或从物理的角度称为 冲量原理.
《偏微分方程》第3章 波动方程
• 下面,我们以求解三维非齐次波动方程的初值问题 为例,说明这个方法的思想.考虑问题
x
x
《偏微分方程》第3章 波动方程
偏微分方程演讲稿ppt课件

偏微分方程
PARTIAL DIFFIERENTIAL EQUATION (P.D.E)
演讲人:Marky
1
目录
• 1 偏微分方程的基本概念 • 2 有限差分方法 • 3 常系数扩散方程及初边值问题 • 4 复金兹堡-朗道方程的简单介绍
深圳大学材料学院
2
1 偏微分方程的基本概念
3
1.1 偏微分方程定义
深圳大学材料学院
17
4 复金兹堡-朗道方程的简单介绍
18
复Ginzburg-Landau方程(CGLE)形式如下:
t
A
A
(1
i
)
2 x
A
(1
i
)
A2
A
其中,A=(x,t)是关于时间t和空间x的复变量;μ是标度参数,通常
情况下,μ=1 ;实数α,β是系统参数。当α,β→∞,α/β=常数,
上方程转变为非线性薛定谔方程。当α,β→0,方程可以化为一个简
, tn1)
u(x j
,tn )
[
u t
]nj
O(
)
(1)
u(x j1, tn ) u(x j , tn ) h
[
u x
]nj
O(h)
(2)
u(x j1, tn )
2u(x j , tn ) h2
u(x j1,t n)
[
2u x 2
]nj
O(h2 )
(3)
深圳大学材料学院
11
利用(1)式和(2)有
1.2.1 偏微分方程的解
偏微分方程的解:如果给定一个函数,将它及它对自变量的各阶偏导 数代入原偏微分方程,能使方程成为恒等式,则称函数是偏微分方程的解。
PARTIAL DIFFIERENTIAL EQUATION (P.D.E)
演讲人:Marky
1
目录
• 1 偏微分方程的基本概念 • 2 有限差分方法 • 3 常系数扩散方程及初边值问题 • 4 复金兹堡-朗道方程的简单介绍
深圳大学材料学院
2
1 偏微分方程的基本概念
3
1.1 偏微分方程定义
深圳大学材料学院
17
4 复金兹堡-朗道方程的简单介绍
18
复Ginzburg-Landau方程(CGLE)形式如下:
t
A
A
(1
i
)
2 x
A
(1
i
)
A2
A
其中,A=(x,t)是关于时间t和空间x的复变量;μ是标度参数,通常
情况下,μ=1 ;实数α,β是系统参数。当α,β→∞,α/β=常数,
上方程转变为非线性薛定谔方程。当α,β→0,方程可以化为一个简
, tn1)
u(x j
,tn )
[
u t
]nj
O(
)
(1)
u(x j1, tn ) u(x j , tn ) h
[
u x
]nj
O(h)
(2)
u(x j1, tn )
2u(x j , tn ) h2
u(x j1,t n)
[
2u x 2
]nj
O(h2 )
(3)
深圳大学材料学院
11
利用(1)式和(2)有
1.2.1 偏微分方程的解
偏微分方程的解:如果给定一个函数,将它及它对自变量的各阶偏导 数代入原偏微分方程,能使方程成为恒等式,则称函数是偏微分方程的解。
偏微分方程的有限元方法市公开课一等奖省赛课获奖PPT课件

第19页
展开J
(Ju(nu)n
)
1 2
(un , un
)
(
f
,
un
)
1 2
n i 1
n
(i , j )cic j
j 1
n
( f , j )c j
j 1
令
J (un ) 0 j 1, 2, , n
c j
则c1, c2 ,, cn满足
n
(i , j )ci ( f , j ) j 1, 2, , n
第1页
偏微分方程有限元方法
一 边值问题变分原理
1 引论 (1)等周问题
在长度一定全部平面封闭曲线中,求所 围面积为最大曲线。
模型:在条件
s2
dx
2
dy
2
ds
l
下
s1 ds ds
求使得泛函 s(x, y) 1 s2 x dy y dx ds
2 s1 ds ds
到达最大函数 x(s), y(s。)
x (a,b)
J (u) 1 (Lu,u) ( f ,u)
2
1
b d p du udx
b
qu
2
dx
b
fudx
2 a dx dx
a
a
1 b ( pu2 qu2 2 fu)dx
2a
引入泛函算子
(u, v)
b
[
p
du
dv
quv]dx
a dx dx
则 J (u) 1 (u,u) ( f ,u)
x2
,,x
n
)T
ann
b (b1, b2 ,,bn )T
则J(x)可表示为:
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
偏微分方程
1.1 基本概念
数学物理方程通常是指物理学、力学、 工程技术和其他学科中出现的偏微分方 程。 反映有关的未知变量关于时间的导数和 关于空间变量的导数之间的制约关系。 连续介质力学、电磁学、量子力学等等 方面的基本方程都属于数学物理方程的 范围。
1.1 基本概念
偏微分方程是指含有未知函数以及未知 函数的某些偏导数的等式。
2.1两个自变量的二阶线性PDE的分类和标准型
定理:设φ(x,y)满足隐函数存在定理中的条件,,则φ(x,y)是方程 (2.1.10)的解的充要条件是φ(x,y)=c是一阶常微分方程
的通积分。 证明: 设φ(x,y)是方程(2.1.10)的解。
2.1两个自变量的二阶线性PDE的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
定解问题的提法是否合适? 例如:这个定解问题的解是否一定存在? 这个定解问题的解是否只有一个? 解的存在性问题 解的唯一性问题
此外,还要考虑解的稳定性问题(或称为解对定解条件或自由 项的连续依赖性问题),即当定解条件或自由项作很小的变化 时,问题的解是否也作很小的变化。 定解问题的存在性、唯一性、稳定性统称为定解问题的适定性。 如果一个定解问题的解是存在的、唯一的、稳定的,称这个问 题是适定的,即认为这样的定解问题的提法是合适的。
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
偶延拓
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
常见的定解条件,可分为初始条件与边界条件。
1.3.1 初始条件
1.3.1 初始条件
1.3.1 初始条件
1.3.1 初始条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
(1.1.1) (1.1.2) (1.1.3) (1.1.4) (1.1.5)
1.1 基本概念
偏微分方程的一般形式
注:F中可以不显含自变量和未知函数,但是, 必须含有未知函数的某个偏导数。 涉及几个未知函数及其偏导数的多个偏微分 方程构成一个偏微分方程组。 注:除非特别说明,一般假设函数u及其在 方程中的各阶偏导数连续。
2.1两个自变量的二阶线性PDE的分类和标准型
椭圆型方程的标准形
2.1两个自变量的二阶线性PDE的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
第三章 波动方程的初值(柯西)问题与行波法
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
1.2 泊松方程的导出
又由库仑定律知,静电场是有势的。即存在静电位势u=u(x,y,z),使 E=-grad u 代入上式,得静电位势u满足以下的泊松方程
即
1.2 拉普拉斯方程和泊松方程的导出
1.3 定解条件与定解问题
一个偏微分方程与定解条件一起构成对于具体问题的完整描 述,称为定解问题。
定解问题中的偏微分方程称为泛定方程。
1.2 热传导方程的导出
1.2 热传导方程的导出
1.2 热传导方程的导出
1.2 热传导方程的导出
下面考虑物体内部有热源(例如物体中通有电流,或有化学反应等情况)。 设在单位时间内单位体积中所产生的热量为F(x,y,z,t),则
则有热源的热传导方程为
1.2 热传导方程的导出
无热源的情况下得到的热传导方程:
2.1两个自变量的二阶线性PDE的分类和标准型
根据复合函数求导法则,有
2.1两个自变量的二阶线性PDE的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
记
的符号是自变量可逆 变换下的不变量
2.1两个自变量的二阶线性PDE的分类和标准型
注:混合型的
退缩的
2.1两个自变量的二阶线性PDE的分类和标准型
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
双曲型方程的第一标准形式
2.1两个自变量的二阶线性PDE的分类和标准型
双曲型方程的第二标准形式 双曲型方程的第一标准形式和第二标准形式统称为双曲型方程的标准形式
2.1两个自变量的二阶线性PDE的分类和标准型
抛物型方程的标准形式
2.1两个自变量的二阶线性PDE的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
1.2 三类经典方程的导出
1.2 热传导方程的导出
例 1.2.2 热传导方程 所谓热传导,就是物体内温度较高的点处的热量 向温度较低点处的流动。 热传导问题归结为求物体内部温度的分布规律。
1.2 热传导方程的导出
设物体在Ω内无热源。 在Ω中任取一封闭曲面S。
以函数u(x,y,z,t)表示物体在t时刻M=M(x,y,z)处的温度。
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
(1.1.1) (1.1.2) (1.1.3)
1.1 基本概念
1.1 基本概念
1.1 基本概念
1.1 基本念
1.1 基本概念
1.1 基本概念
1.2 三类经典方程的导出
例1.2.1 弦的微小横振动问题 弦振动方程是在18世纪由达朗贝尔等人首先给予系 统研究的。 设有一根长为L均匀柔软富有弹性的细弦,平衡时沿 直线拉紧,在受到初始小扰动下,作微小横振动。 试确定该弦的运动方程。
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
分析
特解 通解 通解
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
第二章 二阶线性偏微分方程的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
1.1 基本概念
对于一个非线性偏微分方程,如果它关于未知函数 的最高阶偏导数是线性的,则称它是拟线性偏微分 方程。 例
1.1 基本概念
对于线性偏微分方程而言,将方程中不含未知函数及 其偏导数的项称为自由项。 当自由项为零时,该方程称为齐次方程,否则称为非 齐次方程。 注:齐次、非齐次是对线性偏微分方程而言的。
1.2 三类经典方程的导出
假设: 1. 细弦,就是与张力相比,弦的重量可以忽略不计。 2. 有弹性,表示张力的大小可以按胡可(Hooke)定律来计算。 3. 柔软,是指弦可以弯曲,同时发生于弦中张力的方向总是沿 着弦所在曲线的切线方向。 4. 横振动,是指弦的运动只发生在一个平面上,且弦上各点的 位移与弦的平衡位置垂直。 5. 微小横振动,是指振动的幅度及弦在任意处切线的倾角都很 小。
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
1.3.2边界条件
1.3.2边界条件
1.3.3定解问题
一个偏微分方程与定解条件一起构成对于具 体问题的完整描述,称为定解问题。
1.3.3定解问题
1.3.3定解问题
1.3.3定解问题
1.3.3定解问题
1.3.3定解问题
1.3.3定解问题
1.3.3定解问题
1.3.3定解问题
1.4定解问题的适定性
1.4定解问题的适定性
1.4定解问题的适定性
1.4定解问题的适定性
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
叠加原理的适用范围非常广泛。 叠加原理对于用线性方程和线性定解条件描述的物理现象来说, 都是成立的。 例如,一维热传导方程及其定解问题的叠加原理。
称为齐次热传导方程
有热源的情况下得到的热传导方程: 称为非齐次热传导方程
1.2 热传导方程的导出
1.2 热传导方程的导出
1.2 拉普拉斯方程的导出
1.2 泊松方程的导出
设空间中有一电荷密度为ρ(x,y,z)的静电场。 在此电场内任取一由封闭曲面S包围的区域Ω, 由静电学基本原理知,通过S向外的电通量等于Ω中总电量的4π倍。 即 其中E为电场强度矢量, n为Ω上的单位外法线向量。
1.1 基本概念
(1.1.1) (1.1.2) (1.1.3) (1.1.4) (1.1.5)
1.1 基本概念
如果一个偏微分方程对未知函数及它的所有偏导数都是 线性的,且它们的系数都是仅依赖于自变量的已知函数, 则这样的偏微分方程称为线性偏微分方程。
(1.1.1) (1.1.2) (1.1.3) (1.1.4) (1.1.5)
1.1 基本概念
数学物理方程通常是指物理学、力学、 工程技术和其他学科中出现的偏微分方 程。 反映有关的未知变量关于时间的导数和 关于空间变量的导数之间的制约关系。 连续介质力学、电磁学、量子力学等等 方面的基本方程都属于数学物理方程的 范围。
1.1 基本概念
偏微分方程是指含有未知函数以及未知 函数的某些偏导数的等式。
2.1两个自变量的二阶线性PDE的分类和标准型
定理:设φ(x,y)满足隐函数存在定理中的条件,,则φ(x,y)是方程 (2.1.10)的解的充要条件是φ(x,y)=c是一阶常微分方程
的通积分。 证明: 设φ(x,y)是方程(2.1.10)的解。
2.1两个自变量的二阶线性PDE的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
定解问题的提法是否合适? 例如:这个定解问题的解是否一定存在? 这个定解问题的解是否只有一个? 解的存在性问题 解的唯一性问题
此外,还要考虑解的稳定性问题(或称为解对定解条件或自由 项的连续依赖性问题),即当定解条件或自由项作很小的变化 时,问题的解是否也作很小的变化。 定解问题的存在性、唯一性、稳定性统称为定解问题的适定性。 如果一个定解问题的解是存在的、唯一的、稳定的,称这个问 题是适定的,即认为这样的定解问题的提法是合适的。
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
偶延拓
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
常见的定解条件,可分为初始条件与边界条件。
1.3.1 初始条件
1.3.1 初始条件
1.3.1 初始条件
1.3.1 初始条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
1.3.2边界条件
(1.1.1) (1.1.2) (1.1.3) (1.1.4) (1.1.5)
1.1 基本概念
偏微分方程的一般形式
注:F中可以不显含自变量和未知函数,但是, 必须含有未知函数的某个偏导数。 涉及几个未知函数及其偏导数的多个偏微分 方程构成一个偏微分方程组。 注:除非特别说明,一般假设函数u及其在 方程中的各阶偏导数连续。
2.1两个自变量的二阶线性PDE的分类和标准型
椭圆型方程的标准形
2.1两个自变量的二阶线性PDE的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
第三章 波动方程的初值(柯西)问题与行波法
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
1.2 泊松方程的导出
又由库仑定律知,静电场是有势的。即存在静电位势u=u(x,y,z),使 E=-grad u 代入上式,得静电位势u满足以下的泊松方程
即
1.2 拉普拉斯方程和泊松方程的导出
1.3 定解条件与定解问题
一个偏微分方程与定解条件一起构成对于具体问题的完整描 述,称为定解问题。
定解问题中的偏微分方程称为泛定方程。
1.2 热传导方程的导出
1.2 热传导方程的导出
1.2 热传导方程的导出
1.2 热传导方程的导出
下面考虑物体内部有热源(例如物体中通有电流,或有化学反应等情况)。 设在单位时间内单位体积中所产生的热量为F(x,y,z,t),则
则有热源的热传导方程为
1.2 热传导方程的导出
无热源的情况下得到的热传导方程:
2.1两个自变量的二阶线性PDE的分类和标准型
根据复合函数求导法则,有
2.1两个自变量的二阶线性PDE的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
记
的符号是自变量可逆 变换下的不变量
2.1两个自变量的二阶线性PDE的分类和标准型
注:混合型的
退缩的
2.1两个自变量的二阶线性PDE的分类和标准型
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
双曲型方程的第一标准形式
2.1两个自变量的二阶线性PDE的分类和标准型
双曲型方程的第二标准形式 双曲型方程的第一标准形式和第二标准形式统称为双曲型方程的标准形式
2.1两个自变量的二阶线性PDE的分类和标准型
抛物型方程的标准形式
2.1两个自变量的二阶线性PDE的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
1.2 三类经典方程的导出
1.2 热传导方程的导出
例 1.2.2 热传导方程 所谓热传导,就是物体内温度较高的点处的热量 向温度较低点处的流动。 热传导问题归结为求物体内部温度的分布规律。
1.2 热传导方程的导出
设物体在Ω内无热源。 在Ω中任取一封闭曲面S。
以函数u(x,y,z,t)表示物体在t时刻M=M(x,y,z)处的温度。
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
(1.1.1) (1.1.2) (1.1.3)
1.1 基本概念
1.1 基本概念
1.1 基本概念
1.1 基本念
1.1 基本概念
1.1 基本概念
1.2 三类经典方程的导出
例1.2.1 弦的微小横振动问题 弦振动方程是在18世纪由达朗贝尔等人首先给予系 统研究的。 设有一根长为L均匀柔软富有弹性的细弦,平衡时沿 直线拉紧,在受到初始小扰动下,作微小横振动。 试确定该弦的运动方程。
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
分析
特解 通解 通解
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
第二章 二阶线性偏微分方程的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
2.1两个自变量的二阶线性PDE的分类和标准型
1.1 基本概念
对于一个非线性偏微分方程,如果它关于未知函数 的最高阶偏导数是线性的,则称它是拟线性偏微分 方程。 例
1.1 基本概念
对于线性偏微分方程而言,将方程中不含未知函数及 其偏导数的项称为自由项。 当自由项为零时,该方程称为齐次方程,否则称为非 齐次方程。 注:齐次、非齐次是对线性偏微分方程而言的。
1.2 三类经典方程的导出
假设: 1. 细弦,就是与张力相比,弦的重量可以忽略不计。 2. 有弹性,表示张力的大小可以按胡可(Hooke)定律来计算。 3. 柔软,是指弦可以弯曲,同时发生于弦中张力的方向总是沿 着弦所在曲线的切线方向。 4. 横振动,是指弦的运动只发生在一个平面上,且弦上各点的 位移与弦的平衡位置垂直。 5. 微小横振动,是指振动的幅度及弦在任意处切线的倾角都很 小。
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
3.1一维波动方程的初值(柯西)问题
1.3.2边界条件
1.3.2边界条件
1.3.3定解问题
一个偏微分方程与定解条件一起构成对于具 体问题的完整描述,称为定解问题。
1.3.3定解问题
1.3.3定解问题
1.3.3定解问题
1.3.3定解问题
1.3.3定解问题
1.3.3定解问题
1.3.3定解问题
1.3.3定解问题
1.4定解问题的适定性
1.4定解问题的适定性
1.4定解问题的适定性
1.4定解问题的适定性
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
1.5 线性叠加原理
叠加原理的适用范围非常广泛。 叠加原理对于用线性方程和线性定解条件描述的物理现象来说, 都是成立的。 例如,一维热传导方程及其定解问题的叠加原理。
称为齐次热传导方程
有热源的情况下得到的热传导方程: 称为非齐次热传导方程
1.2 热传导方程的导出
1.2 热传导方程的导出
1.2 拉普拉斯方程的导出
1.2 泊松方程的导出
设空间中有一电荷密度为ρ(x,y,z)的静电场。 在此电场内任取一由封闭曲面S包围的区域Ω, 由静电学基本原理知,通过S向外的电通量等于Ω中总电量的4π倍。 即 其中E为电场强度矢量, n为Ω上的单位外法线向量。
1.1 基本概念
(1.1.1) (1.1.2) (1.1.3) (1.1.4) (1.1.5)
1.1 基本概念
如果一个偏微分方程对未知函数及它的所有偏导数都是 线性的,且它们的系数都是仅依赖于自变量的已知函数, 则这样的偏微分方程称为线性偏微分方程。
(1.1.1) (1.1.2) (1.1.3) (1.1.4) (1.1.5)