 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3120 | 回復(fù): 11 | ||
[求助]
求助LAMMPS兩次運算結(jié)果不一致問題！ 已有1人參與
|
|
本人打算采用LAMMPS計算高聚物的導(dǎo)熱系數(shù)，由于是新手，首先按照LAMMPS例子里面的IN文件進行了練習(xí)（文件路徑：LAMMPSExampleKAPPAin.heat），運算結(jié)果與所提供的參照結(jié)果比較吻合。而考慮到將來我的研究對象是復(fù)雜的高聚物體系，應(yīng)該要用到MS來建模，為了先探探路，于是先將上述in.heat文件內(nèi)描述的模型應(yīng)用MS進行建模，建好模型后在Discover模塊進行設(shè)置（energy里面選擇compass），隨后typing，并順利導(dǎo)出、轉(zhuǎn)換成LAMMPS所能識別的data文件（只有選擇了compass才能順利導(dǎo)出，后來在data文件中修改了pair_coeff）。 我碰到的問題是：從MS-->LAMMPS的方法計算結(jié)果僅為參照結(jié)果的一半（即按照LAMMPS例子中提供的in文件的計算結(jié)果）。我仔細(xì)對比過MS產(chǎn)生data文件與LAMMPS提供的in文件的設(shè)置，關(guān)于模型的參數(shù)（如原子數(shù)目、pair_style、pair_coeff、Mass等參數(shù)均相同），以及后續(xù)對物理過程的處理代碼都是相同的，可是就是結(jié)果差別太大。（附上兩次運算的LOG文件，其中Neighs數(shù)字有較大區(qū)別） 特請這方面經(jīng)驗豐富的朋友幫忙出出主意，或者提點寶貴意見，先感謝啦！@oxox6085@smutao |
回復(fù)此樓» 猜你喜歡
0703化學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
0703化學(xué) 305求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703化學(xué)調(diào)劑，求各位老師收留
已經(jīng)有10人回復(fù)
-
271材料工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 硅納米管在MS中的discover模塊進行不同溫度下的動力學(xué)分析 和 在lammps下的計算
已經(jīng)有3人回復(fù)
- 你會發(fā)表數(shù)據(jù)不完美的結(jié)果嗎？
已經(jīng)有20人回復(fù)
- MS生成Lammps的data文件，文件中原子的類型和順序并不是自己想要的，請問如何修改？
已經(jīng)有8人回復(fù)
- compute 命令出錯，急求
已經(jīng)有5人回復(fù)
- 用lammps做分子動力學(xué)新手求助
已經(jīng)有6人回復(fù)
- 經(jīng)過MS模型導(dǎo)入到LAMMPS，在運算時出現(xiàn)的錯誤求助
已經(jīng)有22人回復(fù)
- 求助LAMMPS關(guān)于MSD計算的問題
已經(jīng)有17人回復(fù)
- lammps計算中出現(xiàn)原子距離太近現(xiàn)象
已經(jīng)有11人回復(fù)
- lammps 離子液體，離子中相鄰隔的挺遠(yuǎn)的原子成鍵了，不知什么原因
已經(jīng)有8人回復(fù)
- LAMMPS 的坐標(biāo)文件怎么構(gòu)建呢？不可能一個一個原子的輸入吧？
已經(jīng)有8人回復(fù)
- lammps計算結(jié)果xyz文件中每個數(shù)字代表的字符長短。
已經(jīng)有10人回復(fù)
- 求助數(shù)值計算結(jié)果和仿真結(jié)果有什么不一樣的
已經(jīng)有6人回復(fù)
- LAMMPS并行計算的問題（cpu——time關(guān)系）
已經(jīng)有17人回復(fù)
- 特定分子振動的不可約表示
已經(jīng)有5人回復(fù)
- 用MS建模導(dǎo)入LAMMPS計算，是不是不能用compass力場？
已經(jīng)有9人回復(fù)
- lammps求助：計算含有不同原子體系的熱導(dǎo)率的in文件該怎么修改
已經(jīng)有7人回復(fù)
- 【討論】lammps 用sw勢計算硅的熔點怎么這么高（2400k),附in文件
已經(jīng)有13人回復(fù)
- 【求助】LAMMPS中非周期性邊界條件和Born-Mayer-Huggins勢函數(shù)不會設(shè)置
已經(jīng)有4人回復(fù)
2樓2014-04-02 17:10:32
|
您好！我重新上傳附件（其中l(wèi)og.MS為利用MS創(chuàng)建模型產(chǎn)生的數(shù)據(jù)文件的計算結(jié)果），請幫我分析分析，謝謝！ |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : Logfiles.rar
2014-04-04 09:21:30, 5.67 K
3樓2014-04-04 09:22:42
4樓2014-04-04 18:06:54
|
請看我上傳到Datafiles，其中in.heat是LAMMPS例子里面的文件，用這個文件直接運行得到 log.LAMMPS；運行in.hydrogen（其中會調(diào)用hydrogen.data）產(chǎn)生 log.MS。請幫忙檢查一下，謝謝！ |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : Datafiles.rar
2014-04-04 19:47:01, 20.23 K
5樓2014-04-04 19:49:52
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
月只藍(lán): 金幣+5, 感謝指導(dǎo)！ 2014-04-07 10:27:25
月只藍(lán): 金幣+5, 感謝指導(dǎo)！ 2014-04-07 10:27:25
|
1、data文件里#慎用，例如“Atoms # full”。 2、region hot block INF INF INF INF -0.4705175*v_rho 0.5294825*v_rho units box region cold block INF INF INF INF 9.5294825*v_rho 10.5294825*v_rho units box 這里設(shè)置可能得改改，否則報錯： Use of region with undefined lattice 3、data文件里力場參數(shù)設(shè)置改一改： Pair Coeffs 1 1.0 1.0 這樣好像可以得到類似結(jié)果。按原來的設(shè)置好像多設(shè)定了個cutoff1。 參考手冊里read_data和pair_style lj/cut部分。 |
6樓2014-04-06 17:04:11
7樓2014-04-06 19:21:17
★ ★ ★ ★ ★
月只藍(lán): 金幣+5, 感謝指導(dǎo)！ 2014-04-07 10:27:33
月只藍(lán): 金幣+5, 感謝指導(dǎo)！ 2014-04-07 10:27:33
|
試過了。你試一下這兩個文件，和你的對比一下。 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : hydrogen.data
- 附件 2 : in.hydrogen
2014-04-06 20:09:09, 353.65 K
2014-04-06 20:09:10, 1.84 K
8樓2014-04-06 20:09:12
9樓2014-04-08 08:35:51
★ ★ ★
月只藍(lán): 金幣+3, 鼓勵交流 2014-04-13 09:29:06
月只藍(lán): 金幣+3, 鼓勵交流 2014-04-13 09:29:06
|
我用的windows版本的lammps，暫時只保留了弛豫部分的run。 例子的in文件不變。自己的data文件里改一下。 你把出問題的log文件發(fā)來看看也行。 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : log.lammps
- 附件 2 : logMS.lammps
2014-04-08 09:32:28, 3.78 K
2014-04-08 09:32:36, 2.34 K
10樓2014-04-08 09:38:21
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zjnu_yc 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷 +10 | 曦熙兮 2026-03-15 | 10/500 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 機械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|