当前位置: 代码迷 >> 综合 >> 定解问题(二)| 定解条件 - 初始条件、边界条件(I、II、III) | 偏微分方程(六)
  详细解决方案

定解问题(二)| 定解条件 - 初始条件、边界条件(I、II、III) | 偏微分方程(六)

热度:72   发布时间:2023-12-21 14:14:22.0

定解条件

背景:在建立方程的过程中,仅考虑了系统内部各部分间的相互作用,以及外界对系统内部的作用。而一个确定的物理过程还要受到历史情况的影响和周围环境通过边界对运动的制约。

泛定方程:反映系统内部作用导出的偏微分方程

定解条件:确定运动的制约条件。

  1. 初始条件(历史情况的影响)

  2. 边界条件(周围环境对边界的影响)

    第I类边界条件(给顶端点值): u ∣ x = x i = μ i ( t ) u|_{x=x_i}=\mu_i(t) ux=xi??=μi?(t)

    第II类边界条件(给定端点梯度): ? u ? n ∣ x = x i = f i ( t ) \frac{\partial u}{\partial n}|_{x=x_i}=f_i(t) ?n?u?x=xi??=fi?(t)

    第III类边界条件(混合I&II): [ a i u + β i ? u ? n ] x = x i = F i ( t ) [a_iu+\beta_i\frac{\partial u}{\partial n}]_{x=x_i}=F_i(t) [ai?u+βi??n?u?]x=xi??=Fi?(t)

  3. 衔接条件(系统内部边界)

定解问题:泛定方程配以适当的定解条件构成一个偏微分方程

1.初始条件

对于随着时间发展变化的物理过程,某一时刻的状态将影响该时刻以后的过程,该时刻的状态便是初始条件。如弦振动问题中,影响弦今后运动的条件有两个
初 位 移 u ∣ t = 0 = ψ ( x ) 初 速 度 ? u ? t ∣ t = 0 = ψ ( x ) 初位移 \quad u|_{t=0}=\psi(x) \\ 初速度 \quad \frac{\partial u}{\partial t}|_{t=0}=\psi(x) ut=0?=ψ(x)?t?u?t=0?=ψ(x)
在热传导问题中,影响今后温度变化的则是初始时刻的温度分布
u ∣ t = 0 = ψ ( x ) u|_{t=0}=\psi(x) ut=0?=ψ(x)
从数学上讲,初始条件是给出未知函数u及其关于某个自变量t的若干阶偏导函数在同一时刻 t = t 0 t=t_0 t=t0?的值。如果方程中关于t的最高阶导数是m阶的,则应给出 u , ? u ? t , . . . , ? m ? 1 u ? t m ? 1 u,\frac{\partial u}{\partial t},...,\frac{\partial^{m-1}u}{\partial t^{m-1}} u,?t?u?,...,?tm?1?m?1u? t = t 0 t=t_0 t=t0?的值。

2.边界条件

在弦振动问题中,对一条有限长的弦( x 1 ≤ x ≤ x 2 x_1\leq x\leq x_2 x1?xx2?),端点 x = x 1 , x = x 2 x=x_1,x=x_2 x=x1?,x=x2?的运动状态对整根弦的运动有制约。以左端点 x = x 1 x=x_1 x=x1?为例,最简单的情况是端点运动状态已知,即
u ∣ x = x 1 = u 1 ( t ) u|_{x=x_1}=u_1(t) ux=x1??=u1?(t)
称为第I类边界条件(给定端点值)。当端点固定在平衡位置时, μ 1 ( t ) ≡ 0 \mu_1(t)\equiv 0 μ1?(t)0,称为第I类齐次边界条件

如果端点负荷已知, x = x 1 x=x_1 x=x1?点受横向外力 F 1 ( t ) = F 1 ( t ) u 。 F_1(t)=F_1(t)\bold u^。 F1?(t)=F1?(t)u。如图所示,在左端点取微元 [ x 1 , x 1 + d x ] [x_1,x_1+dx] [x1?,x1?+dx],在 u 。 \bold u^。 u方向该微元满足牛顿第二定律
ρ d x ? 2 u ? t 2 ∣ x = x 1 = F 1 ( t ) + T s i n θ ∣ x 1 + d x = F 1 ( t ) + T ? u ? x ∣ x 1 + d x = F 1 ( t ) + T ? u ? x ∣ x = x 1 + T ? 2 u ? x 2 d x \rho dx\frac{\partial^2u}{\partial t^2}|_{x=x_1}=F_1(t)+Tsin\theta|_{x_1+dx} \\ =F_1(t)+ T\frac{\partial u}{\partial x}|_{x_1+dx}\\ =F_1(t)+T\frac{\partial u}{\partial x}|_{x=x_1}+T\frac{\partial^2u}{\partial x^2}dx ρdx?t2?2u?x=x1??=F1?(t)+Tsinθx1?+dx?=F1?(t)+T?x?u?x1?+dx?=F1?(t)+T?x?u?x=x1??+T?x2?2u?dx
又有波动方程得,上式方程化简为
0 = F 1 ( t ) + T ? u ? x ∣ x 1 0 = F_1(t)+T\frac{\partial u}{\partial x}|_{x_1} 0=F1?(t)+T?x?u?x1??
得边界条件
? u ? x ∣ x 1 = ? F 1 ( t ) T \frac{\partial u}{\partial x}|_{x_1}=-\frac{F_1(t)}{T} ?x?u?x1??=?TF1?(t)?
这里给出的是 ? u ? x \frac{\partial u}{\partial x} ?x?u?在端点的值,称为第II类边界条件(端点的受力已知,给定端点处的梯度)

类推得

? u ? n ∣ x i = f i ( t ) , i = 1 , 2 \frac{\partial u}{\partial n}|_{x_i}=f_i(t),i=1,2 ?n?u?xi??=fi?(t),i=1,2
在这里插入图片描述
如果端点弹性支撑,即端点除负荷 F 1 ( t ) = F 1 ( t ) u 。 \bold F_1(t)=F_1(t)\bold u^。 F1?(t)=F1?(t)u外,还有弹性力 P → = ? k u ( t , x 1 ) u 。 \overrightarrow P=-ku(t,x_1)\bold u^。 P =?ku(t,x1?)u,k为弹性系数,则端点的运动需满足。
? u ? x ∣ x 1 = ? F 1 ( t ) ? k u ( t , x 1 ) T \frac{\partial u}{\partial x}|_{x_1}=-\frac{F_1(t)-ku(t,x_1)}{T} ?x?u?x1??=?TF1?(t)?ku(t,x1?)?
即混合边界条件为
[ k u ? T ? u ? x ] x 1 = F 1 ( t ) [ku-T\frac{\partial u}{\partial x}]_{x_1}=F_1(t) [ku?T?x?u?]x1??=F1?(t)
此边界条件以 u u u ? u ? x \frac{\partial u}{\partial x} ?x?u?的线性组合给出,称为第III类边界条件(端点弹性支撑)。如果非弹性负荷 F 1 ( t ) ≡ 0 F_1(t)\equiv 0 F1?(t)0,则为第III类齐次边界条件
在这里插入图片描述
对于右端点 x = x 2 x=x_2 x=x2?,可同样导出这三类边界条件,与左端点不同的是在 ? u ? x \frac{\partial u}{\partial x} ?x?u?项前添负号。如果用 n \bold n n表示端点的外法向,则左右两端的三类边界条件可统一表示为 u u u ? u ? n \frac{\partial u}{\partial n} ?n?u?的线性组合
x [ a i u + β i ? u ? n ] x i = F i ( t ) , i = 1 , 2 x[a_iu+\beta_i\frac{\partial u}{\partial n}]_{x_i}=F_i(t),i=1,2 x[ai?u+βi??n?u?]xi??=Fi?(t),i=1,2
边界条件:三维热传导问题

第I类边界条件(当物体边界温度已知时):
u ( t , x , y , z ) ∣ ? V = μ ( t , x , y , z ) ∣ ? V u(t,x,y,z)|_{\partial V}=\mu(t,x,y,z)|_{\partial V} u(t,x,y,z)?V?=μ(t,x,y,z)?V?
当边界温度保持零度时,得第I类齐次边界条件。

当边界上沿外法向 n \bold n n的热流密度 q ( t , x , y , z ) q(t,x,y,z) q(t,x,y,z)已知时,由热传导定律导出第II类边界条件
? u ? n ∣ ? V = ? q ( t , x , y , z ) k ∣ ? V \frac{\partial u}{\partial n}|_{\partial V}=-\frac{q(t,x,y,z)}{k}|_{\partial V} ?n?u??V?=?kq(t,x,y,z)??V?
其中,k为热传导系数。常见的边界绝热情况,相应于第II类齐次边界条件。
在这里插入图片描述
如果物体通过边界与外界自由热交换,在边界面上(x,y,z)处取小面元 d s ds ds,在时间段 [ t , t + d t ] [t,t+dt] [t,t+dt]内从物体内部流入面元 d s ds ds的热量为
Q i n = ? k ? u ? n ∣ ( t , x , y , z ) d s d t Q_{in}=-k\frac{\partial u}{\partial n}|_{(t,x,y,z)}dsdt Qin?=?k?n?u?(t,x,y,z)?dsdt
根据牛顿冷却定律(冷却速率与该温度与温室的温差成正比),从外部流入面元的热量为
Q o u t = h ( T o u t s i d e ? u ) ∣ ( t , x , y , z ) d s d t Q_{out}=h(T_{outside}-u)|_{(t,x,y,z)}dsdt Qout?=h(Toutside??u)(t,x,y,z)?dsdt
h为两种物质间的热交换系数, T o u t s i d e = T o u t s i d e ( x , y , z , t ) T_{outside}=T_{outside}(x,y,z,t) Toutside?=Toutside?(x,y,z,t)为外界的温度

能量守恒定律决定了热量不能在面元上积聚,从而有
? k ? u ? n ∣ ? V d s d t + h ( T o u t s i d e ? u ) ∣ ? V d s d t = 0 -k\frac{\partial u}{\partial n}|_{\partial V}dsdt+h(T_{outside}-u)|_{\partial V}dsdt=0 ?k?n?u??V?dsdt+h(Toutside??u)?V?dsdt=0

[ h u + k ? u ? n ] ? V = h T o u t s i d e ∣ ? V [hu+k\frac{\partial u}{\partial n}]_{\partial V}=hT_{outside}|_{\partial V} [hu+k?n?u?]?V?=hToutside??V?

当外界温度为0时,为第III此齐次边界条件。

在静电场问题中,最常见的是边界接地的情况,此时电位满足第I类齐次边界条件 u ∣ ? V = 0 u|_{\partial V}=0 u?V?=0.

  相关解决方案