应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » [求助：fick第二扩散方程求解]
51/1返回列表
查看: 485  |  回复: 4
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 mytopot 的 8 个金币

mytopot

银虫 (小有名气)

[交流] [求助：fick第二扩散方程求解]

fick第二扩散方程：
dC(x,t)/dt=D*d^2C(x,t)/dx^2+a*C(x,t)/[1+b*C(x,t)]
x=0时，C=Cs。
求C(x,t)

小弟看了好几天了，一直出不来~~望大家帮帮忙~~~~谢谢
回复此楼

» 猜你喜欢
1楼2008-08-26 19:56:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chemofish

木虫 (正式写手)
★ ★ ★ ★
mytopot(金币+2,VIP+0):谢谢
zzgyb(金币+2,VIP+0):谢谢你的参与，欢迎再次光临计算模拟版！
最简单的办法是差分法,一般的Numerical methods for PDEs 的书上都有介绍
而且MATLAB软件可以编个很简单的程序就可以搞定的.
赞一下(12人) 回复此楼
2楼2008-08-26 20:35:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mytopot

银虫 (小有名气)
能不能详细点呢？我以前没用过matlab，麻烦了哈~~~
赞一下 回复此楼
3楼2008-08-26 21:39:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mytopot

银虫 (小有名气)
我用psolve做了一下：
>> y=dsolve('D2y=a*y/(b*y+1)','y(0)=1','y(1)=1','x')
Warning: Explicit solution could not be found.
> In dsolve at 333

y =

[ empty sym ]

感觉还是不对，常数a，b是不是要设置什么？help了dsolve的用法，里面没有关于常数设置的描述~~~~~~~~
赞一下 回复此楼
4楼2008-09-03 19:40:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mytopot

银虫 (小有名气)
在GUI界面中求带常数的微分方程，运行后也显示错误，请问常数如何设置呢？
先选solve中parameters选项，设置use nonlinear solver 中 nonlinear tolerance为1E-1；
在PDE TOOL中设置区域：正方形，bottle=0，left=0，right=1，height=1；
边界条件：左边界=1，右边界=1；
方程类型：椭圆，c=1，a=0，f=α*u/(1+β*u)；
网格加密；
求解时弹出对话框：undifined function or variable 'a'.

在GUI界面中如何设置常数呢？
赞一下 回复此楼
5楼2008-09-03 19:44:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mytopot 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭