 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2612 | 回復(fù): 3 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[求助]
lammps中關(guān)于pressure的問題 已有1人參與
|
|||
一直不理解這個問題：在lammps中，首先做結(jié)構(gòu)優(yōu)化，使用fix box/relax 和 minimize 命令，其中fix box/relax 命令設(shè)定壓強為0 bar，最后一般得到的模型的壓強接近0 bar，但是接著加入溫度后，在nvt系綜下做弛豫，壓強就會變很大，幾千bar， 遠遠大于標準氣壓，我很是不能理解，這種情況下構(gòu)型應(yīng)該不存在了吧？但查閱了很多計算，比如體積熱容和熱脹系數(shù)的計算，都有這個情況，不知道這里面的物理意義是什么？新手求解答， |
回復(fù)此樓» 猜你喜歡
277材料科學與工程080500求調(diào)劑
已經(jīng)有7人回復(fù)
-
317求調(diào)劑
已經(jīng)有18人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
280分求調(diào)劑 一志愿085802
已經(jīng)有7人回復(fù)
-
0854電子信息求調(diào)劑
已經(jīng)有3人回復(fù)
-
263求調(diào)劑
已經(jīng)有4人回復(fù)
-
石河子大學（211、雙一流）碩博研究生長期招生公告
已經(jīng)有3人回復(fù)
-
一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有13人回復(fù)
-
北科281學碩材料求調(diào)劑
已經(jīng)有6人回復(fù)
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 關(guān)于Green-Kubo方法模擬水的粘度
已經(jīng)有5人回復(fù)
- 求助lammps輸出與輸入的問題
已經(jīng)有0人回復(fù)
- 關(guān)于lammps中expected floating point parameter in input
已經(jīng)有13人回復(fù)
- 關(guān)于LAMMPS 如果定義charge
已經(jīng)有0人回復(fù)
- 【求助】關(guān)于LAMMPS建模
已經(jīng)有9人回復(fù)
- 關(guān)于LAmmps的windows下源程序編譯調(diào)試問題
已經(jīng)有7人回復(fù)
- 關(guān)于 run 0 command 重要性
已經(jīng)有0人回復(fù)
- 求助各位前輩關(guān)于在LAMMPS中angle_style的設(shè)置問題。
已經(jīng)有2人回復(fù)
- 關(guān)于towhee 回答zhang_jaj的問題
已經(jīng)有44人回復(fù)
- 一個關(guān)于非周期性邊界條件的問題
已經(jīng)有5人回復(fù)
- 【分享】mdbbs發(fā)布了三年優(yōu)質(zhì)資源匯總電子書
已經(jīng)有48人回復(fù)
- 【分享】mdbbs發(fā)布了三年優(yōu)質(zhì)資源匯總電子書
已經(jīng)有40人回復(fù)
|
你好，謝謝你的回答。您的回答很有啟發(fā)，關(guān)于壓強變動劇烈，主要來自于vriial項的說法，我很贊同，我在mail-list上看到steve也有類似的回復(fù)。其實，我的主要問題還不是壓強變動劇烈，我想換個說法可能更容易表述：比如說，我想測1bar（標準大氣壓），300k室溫下的材料熱導(dǎo)，應(yīng)該怎么做呢？我的做法是，先做minimize和fix_relax 讓pressure為1bar，然后依次做nvt ( 300 k，進入這個過程時，之前0 bar 的系統(tǒng) pressure 陡增到很大），npt (1 bar，體積變化，pressure變動劇烈)，nvt (300 k)，然后再測熱導(dǎo)，不知道這樣做對不對呢？我看到lammps的例子中并沒有考慮得這么復(fù)雜，直接nvt平衡后就開始測熱導(dǎo)了。希望能幫助解答，謝謝。 |
3樓2015-10-23 21:05:53
alundilong
木蟲 (正式寫手)
- 應(yīng)助: 129 (高中生)
- 金幣: 2019.1
- 紅花: 18
- 帖子: 309
- 在線: 56小時
- 蟲號: 595468
- 注冊: 2008-09-08
- 性別: GG
- 專業(yè): 傳熱傳質(zhì)學
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
從宏觀角度看，nvt，密度溫度確定。假設(shè)體系熱平衡，那么兩個狀態(tài)量定下來，所有其他狀態(tài)都定了，包括壓力。然而現(xiàn)在壓力變化很大，這個應(yīng)該是讓樓主困惑的地方吧？ 這個問題不好回答吧。 可能找本統(tǒng)計力學的書看看，會有收獲。 在分子動力學中，我覺得這里的壓力用應(yīng)力代替更加貼切。 這里是計算應(yīng)力的公式 http://lammps.sandia.gov/doc/compute_stress_atom.html 尤其這一段話： compute peratom all stress/atom NULL compute p all reduce sum c_peratom[1] c_peratom[2] c_peratom[3] variable press equal -(c_p[1]+c_p[2]+c_p[3])/(3*vol) thermo_style custom step temp etotal press v_press 說明了壓力就是應(yīng)力張量對角元素的均值（假設(shè)不可壓，這是成立的） 這個compute計算出的壓力就是lammps顯示的壓力。 仔細看看看公式，很大的壓力，其實來自于virial項（動能變化不大），導(dǎo)致了內(nèi)部應(yīng)力變化劇烈，表現(xiàn)出來的就是很大的壓力變化。 另外一個分析角度就是尺度了： 分子動力學計算的空間，時間尺度極小。 時間上：分子松弛時間的倍數(shù)（可以很多很多倍，然而還是不夠大） 空間上：分子自由程的倍數(shù)（可以很多很多倍，然而還是不夠） 平時所講的壓力，其實是個統(tǒng)計量（統(tǒng)計力學）。統(tǒng)計需要足夠的統(tǒng)計樣本，換言之，時間要足夠長，空間也要足夠大�，F(xiàn)在分子動力學很難達到這兩點。那么，在樣本不夠的情況下，壓力變化很大也是可以理解的。 以上為本人拙見，希望有所幫助。 |
2樓2015-10-21 01:56:14
liuxun200376
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 425.2
- 散金: 100
- 紅花: 1
- 帖子: 156
- 在線: 157.4小時
- 蟲號: 1296409
- 注冊: 2011-05-15
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學和
4樓2015-10-26 10:23:33
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +6 | chuanzhu川燭 2026-03-18 | 6/300 |
|
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
|
[考研]
 一志愿東華大學化學070300，求調(diào)劑
+6
|
2117205181 2026-03-21 | 7/350 |
|
|
|
[考研] 考研調(diào)劑 +4 | 來好運來來來 2026-03-21 | 4/200 |
|
|
|
[考研] 化學調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
|
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學大學 ，080500材料科學與工程學碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 08工學調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|