| 查看: 877 | 回复: 0 | ||
[求助]
matlab 解传质方程时出现的问题（分段边界条件）
|
|
大家好，最近在学习着用MATLAB处理传质问题，但在边界条件的处理中遇到一些问题，想问问大家。所要讨论的模型如下： ∂C/∂t+F ∂q/∂t+u ∂C/∂z=D (∂^2 C)/(∂z^2 ) 其中：F、u、D 是常数；q=f(C) 初始条件：C ( z,t = 0)= 0 边界条件：C ( z=0, t) =C0 , 0<=t<=tp C ( z = 0, t)= 0, tp<t 我目前用的代码是这样的 m = 0; x = linspace(0,1,20); t = linspace(0,20,80); sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,x,t); % Extract the first solution component as u. u = sol(:,:,1); % A surface plot is often a good way to study a solution. surf(x,t,u) title('Numerical solution computed with 20 mesh points.') xlabel('Distance x') ylabel('Time t') % A solution profile can also be illuminating. figure plot(t,u(:,1)) %-------------------------------------------------------- function [c,f,s] = pdex1pde(x,t,u,DuDx) a0=1; b0=0.3; k0=1; c0=30; d0=0; c=1; f=-d0*DuDx; s=-(a0/b0)*DuDx-((1-b0)/b0)*k0*(u-c0); %-------------------------------------------------------- function [pl,ql,pr,qr] = pdex1bc(xl,ul,xr,ur,t) pl=ul-((t-7)>=0)*100; ql = 0; pr = 0; qr = 0.1; %-------------------------------------------------------- function u0=pdex1ic(x) u0=0; %-------------------------------------------------------- 但是模拟时出现了错误 Warning: Failure at t=7.000000e+000. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.421085e-014) at time t. > In ode15s at 819 In pdepe at 320 Warning: Time integration has failed. Solution is available at requested time points up to t=6.835443e+000. > In pdepe at 326 ??? Error using ==> surf at 70 Data dimensions must agree. 请问该怎么解决，或有效的实现分段边界条件呢？谢谢 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有167人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
找到一些相关的精华帖子，希望有用哦~
matlab解偏微分方程求助
已经有12人回复
- 求Matlab解非线性方程代码
已经有12人回复
- 用matlab如何解方程？
已经有6人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- 请问matlab拟合二元方程的问题
已经有3人回复
- 关于matlab线性规划 背包问题的 解法 求救呀~~方程知道 但不知道如何解答····
已经有4人回复
- matlab解微分方程
已经有10人回复
- 求matlab 解方程的问题
已经有5人回复
- 高金求助matlab解微分方程组
已经有12人回复
- matlab解微分方程组
已经有15人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- 求助Matlab解方程问题
已经有6人回复
- MATLAB解方程出现??? Subscripted assignment dimension mismatch.
已经有3人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】matlab 解矩阵方程 重金悬赏解决方案
已经有10人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
30