udf据所求流场中某一点的温度来修正边界条件

superongst

虚心请教：

问题背景：想将传热反问题(已知温度场求边界条件)的共轭梯度法用udf实现。

遇到的问题：在假定边界条件的情况下求解温度场，再根据所求流场中某一点的温度来修正边界条件，直到得到合适的解为止。

具体如图：
AB,CD,DA的边界条件已知（绝热），现在要根据O点（中心）的温度值来逐步修正BC的边界条件，来求出跟实验值一致的O点温度：

1.        假设BC边界：h=5000，（环境温度300K，求解区域初始温度400K），计算出0.2s后的温度场（SI）；
2.        判断中心温度To：if（To>=395）h=h+1000；else停止计算，接受h；
3.        If成立时，重新初始化进行计算（放弃已经计算 0.2s数据，用新的h=h+1000重新计算0.2s后的温度场）
4.        回带步骤2

之前方向是有限元方法，最近接触udf，希望得到大家指导，越细节越好。

[ 本帖最后由 superongst 于 2011-10-14 11:00 编辑 ]