大家帮忙看过来看看我设计的耦合过程为啥会有这问题

江南财子

大家可以先去下面链接看看我发的有关耦合问题模型帖子
http://www.simwe.com/forum/viewthread.php?tid=344864

调试过程，出现如下提示：
Error in UMFPACK back substitution.

从HELP找过，但还是不得解，在这有劳大家各位大侠帮忙看看，在这先谢谢了。
下面是我的耦合设计步骤：
1.先选择模型一，PDE coefficient>Time-dependent analysis，定义变量T、u、w（代替下标为u、w的变量）
2.设置除了T、下标为w、u的变量之外的常量
3.建立.015×0.45的矩形面
4.设置子域和边界
子域：依据方程1的形式设置各系数
边界：选择第二类型，设置r。1、2、4的T、u、w各为0、0.1、0.12，边界3的T为-5，u、w为0.1、0.12

5.在选择模型二，PDE coefficient>Time-dependent analysis，定义变量T、w（代替下标为w）
6.设置子域和边界
子域：依据方程2的形式设置各系数，对于方程2右边第一项我是这么设计它的偏微分平方diff(T,x)^2/T^2（不知可否？？？？？）

边界：选择第二类型，设置r。实质同模型一，只是少设置了一个变量u的边界值

7.最后选择模型三，PDE coefficient>Time-dependent analysis，定义变量T、u、w（代替下标为u、w的变量）
8.设置子域和边界
子域：依据方程3的形式设置各系数，对于方程3右边第一项左边那个偏微分我是这么设计的diff(u,T)（不知可否？？？？？）

边界同上述模型一，大家看看我的设计步骤，是不在设计边界时有点画蛇添足之感，还有那些系数的设定是不也有点无章可询之嫌，不管怎么说，希望大家能帮忙多多提出修改建议，这是俺的毕业论文数值模拟最为关键的一部分，在这先感激大家了，学校没人学这软件，只有到这个大家庭来请教大家了。
附件附上源程序：

edodo

你的这个问题中diff(T,x)^2/T^2肯定是不可以的，最好用weak form的方式来构造这类强非线性的问题。
你的fl文件我打开后将第二个方程中所有有关T的系数中的T删除后，就可以计算。你对你的方程最好再推导一下。

江南财子

最近忙着工作的事，今天特来看看——谢谢版主了！
你说的方法我会去努力尝试的，可是weak form都是用来计算有关积分之类的方程啊，不知道版主能否详细指点一下，在这谢谢版主的帮助了。