| 24小時(shí)熱門版塊排行榜 |
| 查看: 3121 | 回復(fù): 11 | ||
[求助]
求助LAMMPS兩次運(yùn)算結(jié)果不一致問(wèn)題! 已有1人參與
|
|
本人打算采用LAMMPS計(jì)算高聚物的導(dǎo)熱系數(shù),由于是新手,首先按照LAMMPS例子里面的IN文件進(jìn)行了練習(xí)(文件路徑:LAMMPSExampleKAPPAin.heat),運(yùn)算結(jié)果與所提供的參照結(jié)果比較吻合。而考慮到將來(lái)我的研究對(duì)象是復(fù)雜的高聚物體系,應(yīng)該要用到MS來(lái)建模,為了先探探路,于是先將上述in.heat文件內(nèi)描述的模型應(yīng)用MS進(jìn)行建模,建好模型后在Discover模塊進(jìn)行設(shè)置(energy里面選擇compass),隨后typing,并順利導(dǎo)出、轉(zhuǎn)換成LAMMPS所能識(shí)別的data文件(只有選擇了compass才能順利導(dǎo)出,后來(lái)在data文件中修改了pair_coeff)。 我碰到的問(wèn)題是:從MS-->LAMMPS的方法計(jì)算結(jié)果僅為參照結(jié)果的一半(即按照LAMMPS例子中提供的in文件的計(jì)算結(jié)果)。我仔細(xì)對(duì)比過(guò)MS產(chǎn)生data文件與LAMMPS提供的in文件的設(shè)置,關(guān)于模型的參數(shù)(如原子數(shù)目、pair_style、pair_coeff、Mass等參數(shù)均相同),以及后續(xù)對(duì)物理過(guò)程的處理代碼都是相同的,可是就是結(jié)果差別太大。(附上兩次運(yùn)算的LOG文件,其中Neighs數(shù)字有較大區(qū)別) 特請(qǐng)這方面經(jīng)驗(yàn)豐富的朋友幫忙出出主意,或者提點(diǎn)寶貴意見(jiàn),先感謝啦!@oxox6085@smutao |

|
您好!我重新上傳附件(其中l(wèi)og.MS為利用MS創(chuàng)建模型產(chǎn)生的數(shù)據(jù)文件的計(jì)算結(jié)果),請(qǐng)幫我分析分析,謝謝! |

|
請(qǐng)看我上傳到Datafiles,其中in.heat是LAMMPS例子里面的文件,用這個(gè)文件直接運(yùn)行得到 log.LAMMPS;運(yùn)行in.hydrogen(其中會(huì)調(diào)用hydrogen.data)產(chǎn)生 log.MS。請(qǐng)幫忙檢查一下,謝謝! |
|
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è)置可能得改改,否則報(bào)錯(cuò): Use of region with undefined lattice 3、data文件里力場(chǎng)參數(shù)設(shè)置改一改: Pair Coeffs 1 1.0 1.0 這樣好像可以得到類似結(jié)果。按原來(lái)的設(shè)置好像多設(shè)定了個(gè)cutoff1。 參考手冊(cè)里read_data和pair_style lj/cut部分。 |

|
試過(guò)了。你試一下這兩個(gè)文件,和你的對(duì)比一下。 |

|
我用的windows版本的lammps,暫時(shí)只保留了弛豫部分的run。 例子的in文件不變。自己的data文件里改一下。 你把出問(wèn)題的log文件發(fā)來(lái)看看也行。 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 288求調(diào)劑 +3 | 于海海海海 2026-03-19 | 3/150 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 9/450 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|