一维PDE传热的问题大家一起学习吧-111.zip

所需积分/C币:14 2019-08-13 12:30:44 16KB ZIP
39
收藏 收藏
举报

一维PDE传热的问题大家一起学习吧-111.zip 本帖最后由 remie 于 2013-3-11 23:49 编辑     最近忙一个课题,是传热传质的一个pde,现在遇到问题了,写个帖子请教请教大家,顺便把这几天对pde的理解和大家分享分享。     本人上课学pde,也只是学了些皮毛,原本是想自己写个差分方程来的,后来发现强大的matlab 似乎是可以解全部的一维PDE 的问题的,这是因为matlab 的帮助系统里面把 一维的pde问题的表示得很宽泛,如下: 1.png 形如2-2的方程,貌似都能解的。对应于doc pdepe  里面的例子(math大哥提示某个会员看的),本人的c、f、s描述的m文件为:function [c,f,s] = pdex1pde c=1;% x:mm       t:s G=1.46; rhov=0.712; R=8.314; muv=12.22*10^; rhol=954;mul=265*^;A=12400;B=610/999;mmax=100/G-66.7; f=*G*permeav11*rhov*R* 273.15)/muv ...    permeal11*rhol/mul*A*B*-20)))^/-20)))^)*DuDx; s=0;复制代码f就是本人那个非常恶心的函数,对应着方程2-2来看。其中涉及3个m文件,permeav11.m 、permeal11.m和 t10585.m,待会一起给出。 要这3个m文件是表示pde里面的变参数,tmd都是2维的变参数,我本想用拟合的方法把这参数用函数表示的,但是拟合的效果并不好,变参数的图如下: 9 空间完整的温度图像.jpg 6 permeability vapor.jpg 8 渗透性插值.jpg 试了下,f里面是可以用m文件来表达的,所以也就这样用来。m文件的调用里面的输入有一些倍数关系,如t10585,这个t/3600,只是表征物理单位不一样,其他有些地方也是同样的含义。 回来说这些m文件,其实是基于数据的插值,二维函数用的自然是interp2插值了。interp2插值,有正确的输入一般都是有输出的,这个可以参考相关的书籍和matlab的帮助系统。m.文件比较多,待会一起给出。 说到这里该把之前有些东西交代下,我有个startup.m,是在启动matlab时就做一些设置的,内容如下%% Set grids on all axis. set; set; set; %% Set the display format format compact; format short; %% And finally put down anything ye have brought forth clear all; %% project worked on open T105_85_file.m open pdex1.m open pdex1pde.m open pdex1ic.m open pdex1bc.m global darvmm darvxx Dar_v ttt2 XXX2 TTT darlmm darlxx Dar_l load myproject.mat commandwindow复制代码前面几行都是系统设置的,后面才是与这次课题相关的一些内容。load了一个mat文件,也就是我的实验数据,并且申明了全局变量,我昨天就在论坛里面问了个全局变量的问题的,这些开机就自动运行的代码其实主要是为调用变参数的m文件服务的。 mat.的文件如下mat文件好像不是论坛上传的允许文件,所以我就打包了下。。: 新建文件夹.zip 对照doc pdepe里面的描述,解pde的主程序m文件如下:function pdex1 clear; load myproject.mat m = 0; x = linspace; t = linspace; sol = pdepe; % Extract the first solution component as u. u = sol; % A surface plot is often a good way to study a solution. figure;surf hold on; plot3 title xlabel ylabel'

...展开详情
立即下载 低至0.43元/次 身份认证VIP会员低至7折
一个资源只可评论一次,评论内容不能少于5个字
您会向同学/朋友/同事推荐我们的CSDN下载吗?
谢谢参与!您的真实评价是我们改进的动力~
关注 私信
上传资源赚钱or赚积分
最新推荐
一维PDE传热的问题大家一起学习吧-111.zip 14积分/C币 立即下载
1/0