7个方程7个未知数，求助，有详图

homogenization · 发表于 2010-2-3 15:37:26

本帖最后由 homogenization 于 2010-2-3 15:41 编辑

7个方程7个未知数，其中4个偏微分，3个不是，方程比较怪，详见图。

我用weak form解，对7个全部写出域和边界的weak form，求解时报错：initial guess leads to undefined function value。

在submains setting中任意给定initial value，能计算。

但是我的问题是稳态的，与给不给初值无关，请大家给我看看，是不是在写weak form有问题，还是可以用另外的方法来做？

kinggun · 发表于 2010-2-5 11:13:40

疑问：
你推导的弱解形式正确吗？

homogenization · 发表于 2010-2-5 14:40:30

谢谢king版。
请king版给我检查一下。
我的弱形式推导1.pdf，以及模型文件pphh.mph. (用的3.2a老版本)

homogenization · 发表于 2010-2-7 12:28:27

自己顶一下，期待达人指正

homogenization · 发表于 2010-2-9 21:14:32

～～TBE_Legen版，你天天在的，帮小弟看看啊

homogenization · 发表于 2010-2-18 19:22:10

自己在顶一下

caoer · 发表于 2010-2-23 06:47:16

后三个不是Pde怎么写弱解？只需要4个弱解公式

shanwent · 发表于 2010-2-23 09:32:25

我好好的研究了你的问题，你的模型有两个错误。
1、模型应该用一维对称建模，而不是2D；
2、把弱形式改成X求导。
由于你的方程是非线性的，只能用迭代法，可能需指定迭代初值。

homogenization · 发表于 2010-2-23 19:20:20

谢谢楼上的兄弟！
1.你说的没有错，我的问题是y方向的一维问题，用一维建模是对的，但是，二维我觉得同样适用。如果采用二维，将左右边界的weak、dweak以及constr设置为零应该是等价的，不知你是不是这个意思？
2.问题是y一维方向的传输，当然，用一维的话，也可以用x的符号，是不是这样理解？
3.是非线性的，迭代初值应该如何确定？有可能要结合问题本身的物理意义来确定？
大家有兴趣的话，帮忙改动一下模型，原模型文件（pphh.mph）我已经贴出来了.

homogenization · 发表于 2010-2-23 19:25:44

7# caoer
谢谢兄弟。
又有PDE又有其他方程的方程组我是同一次遇到。如果弱解形式只是针对PDE的？那么其他三个方程方程，如果利用？有点糊涂了。

wxch · 发表于 2012-1-17 22:32:00

学习一下

wangli_nju · 发表于 2012-3-13 20:49:49

在建模时设定7个变量，然后PDE变成矢量方程，系数是张量形式，这样应该就能解了。

maxiao0 · 发表于 2012-4-21 09:33:44

您为什么选择弱形式求解，这几个方程用系数型PDE方程也能解啊

www8851111 · 发表于 2012-4-22 16:47:58

如果迭代方法的话，初值可以加快收敛，不会影响求解结果吧

maxiao0 · 发表于 2012-4-22 17:15:09

1. 后3个方程把分母移到分子上
2. 最后一个方程还有变量u7 因此你的弱形式推导符号错了
3. 第一个方程原项符号一致当然可能您在设置参数的时候把符号改过了
以上还不算大碍，主要是下面：
4. 边界条件给的不对

王小城 · 发表于 2012-6-9 09:42:40

学习了，以后注意

xxoo还是xxoo · 发表于 2012-8-19 18:35:02

楼主，问题解决了吗？我也在寻求偏微分方程组的解法，也是遇到这种报错情况。我个人认为是方程存在奇点和边界条件不当所致，也在寻求解答中

ydangliu · 发表于 2012-9-7 09:27:05

学习中

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[01.PDE] 7个方程7个未知数，求助，有详图

本帖子中包含更多资源

评分

相关帖子

本帖子中包含更多资源

评分

评分

评分