找回密码
 注册
查看: 3063|回复: 2

怎样把以分段函数形式存在的定压比热容Cp做成UDF或者放到物性里头?

[复制链接]
发表于 2012-8-26 18:26:56 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

x
各位大侠:
    我在fluent计算中定义物性的时候,定压比热容是随温度剧烈变化的,如图11:
    我打算用UDF编程嵌入到material面板里头,发现密度,黏度和导热系数都可以嵌入udf程序但比热容不能。
    后来又考虑放弃udf,使用piecewise-polynomial功能做分段的多项式函数,但是发现piecewise-polynomial功能中的分段函数最多只能分三段(右上角的range值最大为3),而我的比热容分段函数拟合了15段。如下:
   ( y为比热容,x为开尔文温度)
    300K-303K    y = 0.266x2 - 159.446x +23898.34
       303K-304K      y = 3.050x2 - 1846.8x + 279571
       304K-304.4K    y = 21.59x2 - 13119.1x + 1992970
       304.4K-304.75K    y = 261.3x2 - 159095x + 24216588
       304.75K-304.92K  y = 1122.x - 41901
       304.92K-305.05K     y = -1131.x + 345073   
       305.05K-305.3K    y = 574.1x2 - 350556x + 53513892
       305.3K-305.8K    y = 42.49x2 - 25992.5x + 3975141   
       305.8K-307K    y = 3.754x2 - 2306.28x + 354227
       307K-313K      y = 0.139x2 - 86.904x + 13587.26
       313K-323K      y = 0.008x2 - 5.228x +856.082
       323K-333K      y = 0.001x2 - 0.761x + 143.897
       333K-343K       y = -0.020x + 3.511
       343K-350K          y = -0.013x + 6.021
       遇到这种情况我该如何将这个函数关系嵌进去?感谢各位大侠。各位大侠也可以回复我的邮箱forger0819@163.com.

二氧化碳定压比热随温度变化

二氧化碳定压比热随温度变化
 楼主| 发表于 2012-8-26 18:28:11 | 显示全部楼层

图片虽然是叉号,但单击后仍然可以查看,发现实际上里头有个很尖锐的峰值

图片虽然是叉号,但单击后仍然可以查看,发现实际上里头有个很尖锐的峰值,因此不得不分段拟合。
发表于 2012-8-27 09:09:41 | 显示全部楼层
用最新的fluent14就有自定义Cp。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表