24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

scottxhx

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 22
帖子: 31
在线: 24.6小时
虫号: 2164937
注册: 2012-12-04

[交流] Comsol自定义场变量，只要曾经温度T>Ts，对应场变量元素为1，一旦为1，则永久保持已有3人参与

问题：comsol能否自定义一个场变量，来标记每个网格点温度T曾经是否大于Ts，一旦某个网格点温度曾经大于Ts，该网格状态量为1，即便温度降下来了，状态量不变。

这个问题非常具有代表性，其归根结底是要定义场变量f(x,y,z,t)，这个场变量依附于每一个网格，可以初始化，可被用户自定义操作。我目前发现，comsol并不能不依附任何物理场模块自定义这样的变量，退而求其次的做法是采用comsol的数学模块里的域偏微分/域常微分方程/不可逆转变/生物传热模块，这些做法逻辑是能行得通，但是求解非常不准确。采用偏微分或者常微分模块，无论是通过系数将模块数学方程u=f，还是du/dt=f这么简单的方式，其中f通过逻辑控制，即便这样，我发现comsol求解出来的u误差很大，连u=f，这里置f=1或者0，都不能严格控制计算的u=1或者0；du/dt=f的方式求解的误差更离谱，我一直定义f为正，但是u求解出来的还有一点点负值，这样就无法通过u来判断它曾经的状态。不可逆转变/生物传热一样有这样的问题，我觉得这个东西应该很简单，而且计算也不困难，不知道大家有没有什么好的方案。

还要说的是，Matlab函数我也试过，本来想把这样的状态量在matlab里定义为全局场变量，但是发现comsol送过来的网格坐标是分批的，打个比方，模型有50万网格，comsol会分好几拨送过来，也就是Matlab函数在一次时间迭代开始前会被调用好几次，而且送过来的网格批次的先后顺序也不得而知，这就无法对齐场变量元素的位置。

回复此楼