还是关于电弧模拟的

fangdian

大家好！
    还是发一个电弧模拟的问题，我用ansys做电弧模拟，想先模拟电场，再将电场计算的电流密度导入磁场进行计算，得出焦耳热，再将此焦耳热以生热率的方式加载到fluid142单元的flotran计算中，得出电弧的温度场，有关建模的帖子见我以前发的：
http://forum.simwe.com/thread-994870-1-1.html
  我电场计算的结果感觉是对的，如下：

局部放大图：


磁场计算结果似乎也没有太大的问题：

磁场计算结果俯视图


但是流场计算后，温度场差距很大，最低温度和最高温度都太差，如下：


流场计算的速度场为：

流场速度场中的一个面局部放大为：


感觉差距太大，我看过人家文献中计算的温度场为：

  我的问题可能在哪些方面？
  感觉无从查起。

fangdian

在ansys中，我感觉我没有好办法处理电弧的辐射问题，我查到了我所用其它的辐射损失（单位W/M3，随温度变化）
，

所以辐射我采用表格加载生热率的方式（将辐射损失值变为负的生热率,以bf方式加载到点上），如下：
*dim,radiat,table,73,1,1,temp,,     !等离子体辐射项

radiat(1,0,1)=5000,8098.290598,8397.435897,8696.581197,8995.726496,9294.871795,9604.700855,9903.846154,10202.99145,10502.13675,
radiat(11,0,1)=10801.28205,11100.42735,11399.57265,11698.71795,11997.86325,12297.00855,12596.15385,12895.29915,13205.12821,13504.2735,
radiat(21,0,1)=13803.4188,14102.5641,14401.7094,14700.8547,15000,15299.1453,15598.2906,15897.4359,16196.5812,16495.7265,
radiat(31,0,1)=16794.87179,17104.70085,17403.84615,17702.99145,18002.13675,18301.28205,18600.42735,18899.57265,19198.71795,19497.86325,
radiat(41,0,1)=19797.00855,20096.15385,20395.29915,20705.12821,21004.2735,21303.4188,21602.5641,21901.7094,22200.8547,22500,
radiat(51,0,1)=22799.1453,23098.2906,23397.4359,23696.5812,23995.7265,24294.87179,24604.70085,24903.84615,25202.99145,25502.13675,
radiat(61,0,1)=25801.28205,26100.42735,26399.57265,26698.71795,26997.86325,27297.00855,27596.15385,27895.29915,28205.12821,28504.2735,
radiat(71,0,1)=28803.4188,29102.5641,29401.7094

radiat(1,1,1)=0,-1335566.76,-2268284.633,-4756815.357,-12623258.01,-26864575.27,-54170060.-22,-110848146.2,-191051371.5,-311991459,
radiat(11,1,1)=-517041052.1,-784456182.4,-1196030567,-1677672305,-2296266893,-2877391359,-3587943375,-4745165140,-5417006022,-6526742603,
radiat(21,1,1)=-8299708372,-9290757363,-10148217728,-10976618837,-12227197960,-12287311518,-12968389593,-13553622272,-13487313507,-13487313507,
radiat(31,1,1)=-13822133392,-13687219403,-13487313507,-13553622272,-13620257037,-13620257037,-12968389593,-12841808759,-12904943977,-12904943977,
radiat(41,1,1)=-12968389593,-13225306655,-13553622272,-13421329145,-13890088260,-14804509950,-15472602389,-16170844255,-16817912863,-18190762563,
radiat(51,1,1)=-19293455187,-20462991134,-22461445394,-24776286181,-27062933282,-29851996081,-32607089391,-36144354123,-39869333689,-43123875632,
radiat(61,1,1)=-48749023556,-52987649376,-60787259812,-66397423113,-72170540429,-78831286752,-86106765088,-96862457770,-1.04257E+11,-1.12767E+11,
radiat(71,1,1)=-1.24389E+11,-1.37208E+11,-1.51349E+11

esel,s,mat,,1   
nsle,s
bf,all,hgen,%radiat%

不知道这样加载能否合适？
  我采用的是bf加载到点上，也试过bfe加载到单元上，只是也是有问题（由于我先将磁场计算的生热率加载到单元上，如果再用bfe将辐射率加载到单元上会将原来的值替换掉，所以最终选择了bf）
  唉，搞了好久了，就是没有出来，谢谢啊

wgq521

你的阴极没有尖端，电弧在阳极与阴极之间，你这个是不是没有起弧？

wgq521

而且你的思路，弧是由电极间的压差产生，产生电弧以后，温度升高，阳极与阴极之间导通，你单独计算电流的话，那温度怎么办呢，这是一个电磁场与流场耦合计算的问题，不能分别计算。

fangdian

我初次计算的时候给氩气一个初始高温6000℃，就能引弧了，至于温度场，我打算通过电磁场计算出生热率导入流场计算温度，再把温度返回到电磁场，这样的话就能考虑温度了。

huangwei30

楼主  我也是做焊接电弧模拟的，只是现在懂的东西太少，希望能与楼主交流

taojunhust

有些复杂，先学习了先。

fangdian

huangwei30 发表于 2011-9-8 14:18
楼主  我也是做焊接电弧模拟的，只是现在懂的东西太少，希望能与楼主交流 ...

可以啊，我qq59158681

andylove55

楼主我是个新手，能把你提到的文献名跟我说下吗，谢谢楼主

小七工作室

我们计算过不少电弧 包括ansys高斯热源  fluent激光焊接  vof焊接等

simonsnow

有些复杂，先学习了先。

fishgo

到现在对电类还不怎么了解，好难的啊...楼主人才.

chenqian92

支持一个，楼主我最近也在做电弧模拟，新手一个，以后希望多多指教

chenqian92

楼主，我现在用热电耦合来仿电弧产生的整个过程，我想在这个热电场中加入外电路，请教一下如何实现？
或者说因为我在边界处加载的电流跟我计算过程中的场变量有关，我如何将这个场变量即时反馈过来？

神愤业炎

好复杂的样子哦