24小时热门版块排行榜    

查看: 1338  |  回复: 6

bluethky

木虫 (著名写手)

[求助] 有界的扩散方程如何通过调和函数分析法来求解 已有1人参与

求解的是一个有界(-L,L)的fick扩散方程,作者说可以通过harmonic analysis ,就得到扩散方程的解,用一个级数的形式表示的,没有看出来怎么来的,貌似是离散的fourier变换,不过我也没有算出来
求大侠指点,谢谢!


附件是图片,已经相关的文献有界的扩散方程如何通过调和函数分析法来求解
01.jpg


有界的扩散方程如何通过调和函数分析法来求解-1
02.jpg
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : supplementary.pdf
  • 2014-05-31 08:57:45, 112.63 K

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xcenxcen

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
这是一个两点边值问题
P=X(x)*T(t)
XT'=DX''T
T'/(D*T)=X''/X=-lamda^2
X''+lamda^2*X=0
X=C1*cos(lamda*x)+C2*sin(lamda*x)
利用边值条件: X'=C1*lamda*(-sin(lamda*x))+C2*cos(lamda*x)*lamda
x=0, X'=0 =====> C2=0
x=L, X'=0 =====> lamda*sin(lamda*L)=0
lamda*L=k*pi
lamda=k*pi/L
Xk(x)=C1*cos(k*pi/L*x)
Tk(t)=C3*exp(-k^2*pi^2*D/L^2*t)
Pk(x,t)=Ck*cos(k*pi/L*x)*exp(-k^2*pi^2*D/L^2*t)
P=sigma(Ck*cos(k*pi/L*x)*exp(-k^2*pi^2*D/L^2*t))
P(x,ts;xs,ts)=delta(x-xs)
sigma(Ck*cos(k*pi/L*x)*exp(-k^2*pi^2*D/L^2*ts))=delta(x-xs)
sigma(Ck*<cos(k*pi/L*x),cos(n*pi/L*x)>*exp(-k^2*pi^2*D/L^2*ts)=cos(n*pi/L*xs)
Cn*exp(-k^2*pi^2*D/L^2*ts)*integrate(0,pi)dx(1+cos(2*n*x))/2*L/pi=cos(n*pi/L*xs)
Cn*exp(-k^2*pi^2*D/L^2*ts)*L/2=cos(n*pi/L*xs)
Cn=2/L*cos(n*pi/L*xs)*exp(k^2*pi^2*D/L^2*ts)
P(x,t;xs,ts)=sigma(k=0,+infinity)2/L*cos(n*pi/L*xs)*cos(n*pi/L*x)*exp(-k^2*pi^2*D*(t-ts)/L^2).......
待我长发及腰,遮住一身肥膘。纵然虎背熊腰,也要高冷傲娇。
2楼2014-05-31 17:08:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluethky

木虫 (著名写手)

引用回帖:
2楼: Originally posted by xcenxcen at 2014-05-31 17:08:32
这是一个两点边值问题
P=X(x)*T(t)
XT'=DX''T
T'/(D*T)=X''/X=-lamda^2
X''+lamda^2*X=0
X=C1*cos(lamda*x)+C2*sin(lamda*x)
利用边值条件: X'=C1*lamda*(-sin(lamda*x))+C2*cos(lamda*x)*lamda
x=0, X'=0 = ...

你用的分离变量法吗?先谢谢再说我好好看看的你如何算的
3楼2014-05-31 20:47:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluethky

木虫 (著名写手)

送红花一朵
首先先谢谢你的提示,我有几点不太懂,
1)T'/(D*T)=X''/X=-lamda^2     为什么他们的值要设为 -lamda^2,是从题中隐含的吗?

2)sigma(Ck*<cos(k*pi/L*x),cos(n*pi/L*x)>*exp(-k^2*pi^2*D/L^2*ts)=cos(n*pi/L*xs),这个式子的含义是什么?  不太理解

3)Cn*exp(-k^2*pi^2*D/L^2*ts)*integrate(0,pi)dx(1+cos(2*n*x))/2*L/pi=cos(n*pi/L*xs) ,这个我也不懂,式子中引入了一个(0,pi)的积分,没有明白您的思路
4楼2014-06-01 13:16:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluethky

木虫 (著名写手)

昨天晚上我按您的思路计算了一下,自己不懂的地方,再向您请教一下,先谢谢啦
5楼2014-06-01 13:17:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xcenxcen

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
bluethky: 金币+10, ★★★★★最佳答案 2014-06-01 23:44:46
引用回帖:
4楼: Originally posted by bluethky at 2014-06-01 13:16:55
首先先谢谢你的提示,我有几点不太懂,
1)T'/(D*T)=X''/X=-lamda^2     为什么他们的值要设为 -lamda^2,是从题中隐含的吗?

2)sigma(Ck*<cos(k*pi/L*x),cos(n*pi/L*x)>*exp(-k^2*pi^2*D/L^2*ts)=cos(n ...

1 设为-lamda^2比较好算,左边是纯t的函数 右边是纯x的函数 所以应该和t,x都没有关系,所以就得是常数啊。
2 这个就是两边同乘以cos(n*pi/L*x),然后在[0,L]上积分。
3 左边的级数,除了k=n的时候不是零以外,其他项都是零,右边呃呃呃。
待我长发及腰,遮住一身肥膘。纵然虎背熊腰,也要高冷傲娇。
6楼2014-06-01 17:45:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluethky

木虫 (著名写手)

引用回帖:
6楼: Originally posted by xcenxcen at 2014-06-01 17:45:46
1 设为-lamda^2比较好算,左边是纯t的函数 右边是纯x的函数 所以应该和t,x都没有关系,所以就得是常数啊。
2 这个就是两边同乘以cos(n*pi/L*x),然后在上积分。
3 左边的级数,除了k=n的时候不是零以外,其他项 ...

好的,谢谢!我在好好看看,谢谢指点!
7楼2014-06-01 23:44:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 逃课的小屁孩 的主题更新
信息提示
请填处理意见