| 查看: 886 | 回复: 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. 请问该怎么解决，或有效的实现分段边界条件呢？谢谢 |
回复此楼» 猜你喜欢
散金币祈福
已经有187人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有109人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有118人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
找到一些相关的精华帖子，希望有用哦~
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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
