24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3223  |  回复: 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的回帖

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的回帖

onnmel

新虫 (小有名气)

【答案】应助回帖

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

chenboxin

新虫 (初入文坛)

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

你好你现在还做相场模拟吗
6楼2018-05-02 20:50:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenboxin

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by qinbao0401 at 2017-03-23 12:52:43
你好,你会三维的相场吗?...

你好你还在做相场模拟吗
7楼2018-05-02 20:51:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenboxin

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by onnmel at 2018-03-21 09:17:44
三维的需要并行 目前matal不知道能不能实现 大部分都是用Fortran

你好现在还做相场模拟吗
8楼2018-05-02 20:51:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

繁星半夏

金虫 (小有名气)

我既用了有限差分法也用了半隐式傅里叶,发现用傅里叶时参数只在很小的一个取值范围内稳定,略取大一点便震荡十分强烈,请问这是什么原因?你也有相似的体验吗?
9楼2018-07-28 01:40:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

itmtwlig

新虫 (初入文坛)

你好,这个附件可以再分享一下吗,现在下载不了了

发自小木虫Android客户端
10楼2021-02-27 16:01:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 husthsc 的主题更新
信息提示
请填处理意见