Viscosity of Asphalt Binder through Equilibrium and Non-Equilibrium Molecular Dynamics Simulations

Viscosity is a curial indicator for evaluating asphalt performance, representing its ability to resist deformation under external forces. The Green–Kubo integral in equilibrium molecular dynamics simulations and the Muller-Plathe algorithm in reverse non-equilibrium molecular dynamics simulations we...

Full description

Saved in:
Bibliographic Details
Main Authors: Xiancheng Hu, Xiaohan Huang, Yuanbin Zhou, Jiandong Zhang, Hongquan Lu
Format: Article
Language:English
Published: MDPI AG 2024-09-01
Series:Buildings
Subjects:
Online Access:https://www.mdpi.com/2075-5309/14/9/2827
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1850258303676841984
author Xiancheng Hu
Xiaohan Huang
Yuanbin Zhou
Jiandong Zhang
Hongquan Lu
author_facet Xiancheng Hu
Xiaohan Huang
Yuanbin Zhou
Jiandong Zhang
Hongquan Lu
author_sort Xiancheng Hu
collection DOAJ
description Viscosity is a curial indicator for evaluating asphalt performance, representing its ability to resist deformation under external forces. The Green–Kubo integral in equilibrium molecular dynamics simulations and the Muller-Plathe algorithm in reverse non-equilibrium molecular dynamics simulations were used to calculate the asphalt viscosity. Meanwhile, the key parameters of both methods were rationalized. The results show that in equilibrium calculations, using a 1/t weighting for the viscosity integral curve results in a well-fitted curve that closely matches the original data. The isotropy of the asphalt model improves for atomic counts exceeding 260,000, rendering viscosity calculations more reasonable. When the viscosity did not converge, it increased linearly with the number of atoms. In non-equilibrium calculations, the number of region divisions had almost no effect on the viscosity value. A momentum exchange period of 20 timesteps exhibits a favorable linear trend in velocity gradients, and an ideal momentum exchange period was found to be between 10 and 20 timesteps. As the model size increased, the linear relationship with the shear rate became more pronounced, and the isotropy of the asphalt system improved. Using an orthogonal simulation box with a side length of 75 Å effectively meets the computational requirements.
format Article
id doaj-art-4792b91fa9a54f48aeda47d6d99749ed
institution OA Journals
issn 2075-5309
language English
publishDate 2024-09-01
publisher MDPI AG
record_format Article
series Buildings
spelling doaj-art-4792b91fa9a54f48aeda47d6d99749ed2025-08-20T01:56:12ZengMDPI AGBuildings2075-53092024-09-01149282710.3390/buildings14092827Viscosity of Asphalt Binder through Equilibrium and Non-Equilibrium Molecular Dynamics SimulationsXiancheng Hu0Xiaohan Huang1Yuanbin Zhou2Jiandong Zhang3Hongquan Lu4ZCCC Hongtu Transportation Construction Co., Ltd., Hangzhou 310051, ChinaCollege of Civil Engineering and Architecture, Zhejiang University, Hangzhou 310058, ChinaCollege of Civil Engineering and Architecture, Zhejiang University, Hangzhou 310058, ChinaCollege of Civil Engineering and Architecture, Zhejiang University, Hangzhou 310058, ChinaCollege of Civil Engineering and Architecture, Quzhou University, Quzhou 324000, ChinaViscosity is a curial indicator for evaluating asphalt performance, representing its ability to resist deformation under external forces. The Green–Kubo integral in equilibrium molecular dynamics simulations and the Muller-Plathe algorithm in reverse non-equilibrium molecular dynamics simulations were used to calculate the asphalt viscosity. Meanwhile, the key parameters of both methods were rationalized. The results show that in equilibrium calculations, using a 1/t weighting for the viscosity integral curve results in a well-fitted curve that closely matches the original data. The isotropy of the asphalt model improves for atomic counts exceeding 260,000, rendering viscosity calculations more reasonable. When the viscosity did not converge, it increased linearly with the number of atoms. In non-equilibrium calculations, the number of region divisions had almost no effect on the viscosity value. A momentum exchange period of 20 timesteps exhibits a favorable linear trend in velocity gradients, and an ideal momentum exchange period was found to be between 10 and 20 timesteps. As the model size increased, the linear relationship with the shear rate became more pronounced, and the isotropy of the asphalt system improved. Using an orthogonal simulation box with a side length of 75 Å effectively meets the computational requirements.https://www.mdpi.com/2075-5309/14/9/2827highway asphaltviscositymolecular dynamics simulation
spellingShingle Xiancheng Hu
Xiaohan Huang
Yuanbin Zhou
Jiandong Zhang
Hongquan Lu
Viscosity of Asphalt Binder through Equilibrium and Non-Equilibrium Molecular Dynamics Simulations
Buildings
highway asphalt
viscosity
molecular dynamics simulation
title Viscosity of Asphalt Binder through Equilibrium and Non-Equilibrium Molecular Dynamics Simulations
title_full Viscosity of Asphalt Binder through Equilibrium and Non-Equilibrium Molecular Dynamics Simulations
title_fullStr Viscosity of Asphalt Binder through Equilibrium and Non-Equilibrium Molecular Dynamics Simulations
title_full_unstemmed Viscosity of Asphalt Binder through Equilibrium and Non-Equilibrium Molecular Dynamics Simulations
title_short Viscosity of Asphalt Binder through Equilibrium and Non-Equilibrium Molecular Dynamics Simulations
title_sort viscosity of asphalt binder through equilibrium and non equilibrium molecular dynamics simulations
topic highway asphalt
viscosity
molecular dynamics simulation
url https://www.mdpi.com/2075-5309/14/9/2827
work_keys_str_mv AT xianchenghu viscosityofasphaltbinderthroughequilibriumandnonequilibriummoleculardynamicssimulations
AT xiaohanhuang viscosityofasphaltbinderthroughequilibriumandnonequilibriummoleculardynamicssimulations
AT yuanbinzhou viscosityofasphaltbinderthroughequilibriumandnonequilibriummoleculardynamicssimulations
AT jiandongzhang viscosityofasphaltbinderthroughequilibriumandnonequilibriummoleculardynamicssimulations
AT hongquanlu viscosityofasphaltbinderthroughequilibriumandnonequilibriummoleculardynamicssimulations