今天把之前关于黑油模型的混合元形式整理了一下,由于以前把所有的参数之间的依赖关系都搞明白了,这次重新回头看以前写的形式,才发现了之前老板一直强调的非线性的重要性。黑油模型真是繁琐,看着那么将近一页A4纸的方程和眼花缭乱的参数就有些头晕,现在真是佩服自己当初把这些理顺了.当然难解的程度可想而知,先不说过泡点问题,就是单独的强非线性就搞得人头晕眼花.
放假前写了代数多重网格(AMG)以及Uzawa算法的程序,对于线性情况来说,计算到几百万上千万的时候迭代步数也只有20几步,而且这与网格无关.可以说,对于线性情况来说,这已经是最优的不能最优的了.今天在思考能不能把uzawa也用到非线性的情况上来.但是重写方程之后发现,我的形式和一般的非线性情况有些出入,好像是简单一点,但看起来又有些奇怪,自己也拿不准.
周五要开始上讨论班,看看老师的意见,这两天再好好看看这方面的论文和算法.
Comments