| 查看: 2342 | 回复: 13 | |||||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 yakamozyk 的 20 个金币 | |||||
[交流]
comsol livelink with matlab 使用交流
|
|||||
|
最近在使用comsol&MATLAB时候想要尝试通过类似于tutorials里，案例“pseudoperiodicity llmatlab”那样用计算区域的边界条件替换来计算更大的计算区域，案例中计算的是一个在圆形管道中流动的流体受到底部加热的情况，计算总区域有6m长，为了提高效率，使用了一个1m长的计算区域，在计算完成后将出口条件赋值给入口条件，通过这样的方式计算下一段，一直计算到最后一段，计算的变量是流体的温度，速度使用的是管内均匀流动。我想使用类似的方法来计算流动，想通对过一段较短的计算区域反复计算来求出整个空间分布，但是在边界条件传递的时候会有问题，为了简单，我使用圆形流道，壁面光滑无滑移，法相速度入口边界，压力出口条件，计算可以进行但是在MATLAB中想要通过边界条件传递来反复计算时出现问题，报错提示是我给新的入口边界赋值是出现错误，我给入口边界赋值的语句是 inlet1=model.physics('spf').feature('inlet1'); inlet1.set('U0in',3,'inletVelo(x,y)'); 其中inletVelo(x,y)是我仿照案例中对出口温度数据进行提取时用的插值函数 案例中的语句是这样的 %添加一个名为inletTemp的插值函数。通过读取后面建立的包含出口边界温度分布的文本文件。 int1=model.func.create('int1','Interpolation'); int1.model('comp1'); int1.set('source','file'); filename=fullfile(tempdir,'pp_data.txt'); int1.set('filename',filename); int1.set('nargs','2'); model.func('int1').setIndex('funcs','inletTemp',0,0); %% %为inlteTemp(x,y)插值函数设置入口边界温度 temp1=model.physics('ht').feature('temp1'); temp1.set('T0',1,'inletTemp(y,z)'); %% 我仿照案例写的插值函数定义和使用是这样的 %添加一个名为inletVelo的插值函数。通过读取后面建立的包含出口边界速度分布的文本文件。 int1=model.func.creat('int1','Interpolation'); int1.model('comp1'); int1.set('source','file'); filename=fullfile(tempdir,'Untitled1_datav.txt'); int1.set('filename',filename); int1.set('nargs','2'); model.func('int1').setIndex('funcs','inletVelo',0,0); %% %% %为inlteVelo(x,y)插值函数设置入口边界速度 inlet1=model.physics('spf').feature('inlet1'); inlet1.set('U0in',3,'inletVelo(x,y)'); 在运行时这里报错 提示这句话出错，还有A scalar value expected 我在考虑是不是要被赋值的变量应该被赋予一个标量而不是这样的一个分布，如果是这样那我想用那种方法进行计算该怎么做呢，出口边界条件或者进口条件该选哪一种 这个问题也搞了有几天了一直没搞出来，希望有会的朋友能帮帮忙 |
回复此楼» 收录本帖的淘贴专辑推荐
 科研软件及应用 | comsol with matlab |
» 猜你喜欢
工材口青年基金大概什么样能上会？
已经有13人回复
-
想要用氢化钠拔掉吲哚N上的氢取代酰氯
已经有4人回复
-
刚刚收到科研之友邮件
已经有27人回复
-
江南大学到瑞士招聘，称取消非升即走，改预聘+长聘
已经有23人回复
-
怎么带研究生？
已经有8人回复
-
国产期刊影响因子大于12的有20多个个了
已经有6人回复
-
求机械类四区sci推荐
已经有7人回复
-
投稿求助
已经有4人回复
-
初始合伙人来啦！（生物试剂耗材标准品）
已经有26人回复
-
ACS AMI 返回审稿意见，一个大修，两个据稿，编辑给的修改重投
已经有6人回复
头文字d
铁虫 (正式写手)
- 应助: 9 (幼儿园)
- 金币: 1373.8
- 散金: 435
- 红花: 6
- 沙发: 1
- 帖子: 717
- 在线: 70.4小时
- 虫号: 3581083
- 注册: 2014-12-07
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2016-02-25 09:39:35
2楼2016-02-24 09:15:00
 本帖仅楼主可见 |
3楼2016-02-24 09:27:20
relive012429
铁杆木虫 (职业作家)
- 应助: 236 (大学生)
- 贵宾: 0.283
- 金币: 8609.8
- 散金: 4004
- 红花: 77
- 沙发: 1
- 帖子: 3207
- 在线: 1304.4小时
- 虫号: 2210325
- 注册: 2012-12-27
- 性别: GG
- 专业: 传热传质学
4楼2016-02-24 12:08:56
纳米之心
木虫 (文坛精英)
- 应助: 10 (幼儿园)
- 金币: 3978.7
- 散金: 4216
- 红花: 25
- 帖子: 12436
- 在线: 339.5小时
- 虫号: 2445904
- 注册: 2013-05-03
- 专业: 细胞信号转导
6楼2016-02-25 10:01:57
7楼2016-03-01 09:49:26
relive012429
铁杆木虫 (职业作家)
- 应助: 236 (大学生)
- 贵宾: 0.283
- 金币: 8609.8
- 散金: 4004
- 红花: 77
- 沙发: 1
- 帖子: 3207
- 在线: 1304.4小时
- 虫号: 2210325
- 注册: 2012-12-27
- 性别: GG
- 专业: 传热传质学
8楼2016-03-01 15:52:12
9楼2016-03-03 19:13:38
relive012429
铁杆木虫 (职业作家)
- 应助: 236 (大学生)
- 贵宾: 0.283
- 金币: 8609.8
- 散金: 4004
- 红花: 77
- 沙发: 1
- 帖子: 3207
- 在线: 1304.4小时
- 虫号: 2210325
- 注册: 2012-12-27
- 性别: GG
- 专业: 传热传质学
吾心需心吾10楼2016-03-03 19:44:11
50