24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3230  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

husthsc

金虫 (初入文坛)

[求助] 半隐式傅里叶谱法求解相场方程已有2人参与

参考附件中文章利用半隐式傅里叶谱法求解TDGL相场方程,用的matlab。

---------------------------------------(1)
其中f(n)为:
----------------------------------------(2)

先用有限差分法求解看起来没什么问题,但是用半隐式傅里叶谱方法出问题,计算结果发散了。
其方法是将(1)式进行傅里叶变换变成(3),在傅里叶空间中用explicit Euler (4)求解出n,最后再讲n逆傅里叶变换回来。
-------------------------(3)
----------------------(4)

附件中代码是我写的有限差分和半隐式傅里叶谱法代码,求大神帮解决问题。半隐式傅里叶谱法求解相场方程
1.png


半隐式傅里叶谱法求解相场方程-1
2.png


半隐式傅里叶谱法求解相场方程-2
3.png


半隐式傅里叶谱法求解相场方程-3
4.png
回复此楼

» 本帖附件资源列表

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onnmel

新虫 (小有名气)

【答案】应助回帖

三维的需要并行 目前matal不知道能不能实现 大部分都是用Fortran
5楼2018-03-21 09:17:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

shunshierfa

新虫 (小有名气)

【答案】应助回帖

你好,请问你的问题解决了吗?我把你的程序的时间步长缩小了10倍,程序可以运行了,但结果与差分的有所不同。我目前也在做这方面的东西,可以交流交流吗?
2楼2016-04-23 10:20:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qinbao0401

新虫 (小有名气)

引用回帖:
2楼: Originally posted by shunshierfa at 2016-04-23 10:20:33
你好,请问你的问题解决了吗?我把你的程序的时间步长缩小了10倍,程序可以运行了,但结果与差分的有所不同。我目前也在做这方面的东西,可以交流交流吗?

你好,你会三维的相场吗?
3楼2017-03-23 12:52:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qinbao0401

新虫 (小有名气)

引用回帖:
2楼: Originally posted by shunshierfa at 2016-04-23 10:20:33
你好,请问你的问题解决了吗?我把你的程序的时间步长缩小了10倍,程序可以运行了,但结果与差分的有所不同。我目前也在做这方面的东西,可以交流交流吗?

你好,你会三维的相场吗?
4楼2017-03-23 13:58:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见