The University of Southampton
University of Southampton Institutional Repository

Simulation of Lubricant Properties and their Interactions with Surfaces

Simulation of Lubricant Properties and their Interactions with Surfaces
Simulation of Lubricant Properties and their Interactions with Surfaces
The behaviour of lubricants at operational conditions, such as at high temperatures and pressures, is a topic of great industrial interest. In particular, viscosity and the viscosity-pressure relation are especially important for applications and their determination by computational simulations is very desirable.

In this thesis we evaluate methods to compute these quantities based on fully atomistic molecular dynamics simulations which are computationally demanding but also have the potential to be most accurate. We tested several molecules that are used as lubricants, such as 9,10-dimethyloctadecane, main component of PAO-2 base oil, which was used as the main lubricant for our simulations. The methods used for the viscosity simulations are the Green-Kubo equilibrium molecular dynamics (EMD-GK), the direct computation of viscosity from shear during non-equilibrium MD (NEMD) and the use of confined NEMD, where the fluid is confined within explicitly defined iron oxide wall surfaces, at pressures of up to 1.0 GPa and various temperatures (40-150 degrees Celsius). We present the theory behind these methods and investigate how the simulation parameters affect the results obtained, to ensure viscosity convergence with respect to the simulation intervals and all other parameters. We show that by using each method in its regime of applicability, we can achieve good agreement between simulated and measured values. NEMD simulations at high pressures captured zero shear viscosity successfully, while at 40 degrees Celsius EMD-GK is only applicable to pressures up to 0.3 GPa, where the viscosity is lower. In NEMD, longer and multiply repeated simulations reduce the standard deviation of viscosity, which is essential at lower pressures.

Additionally, by using confined NEMD simulations, it was demonstrated that the film thickness of the fluid affects viscosity, and as we increase the number of lubricant molecules, we approach the viscosity value of the bulk fluid derived from NEMD simulations.

Another aspect of these methods is the choice of the utilised force field for the atomic interactions. This was investigated by selecting three different commonly used force fields. We have explored several methods for calculating viscosity and we obtained results of particular industrial interest.
University of Southampton
Mathas, Dimitrios
aee09891-6e68-489b-982a-cf73908c16d2
Mathas, Dimitrios
aee09891-6e68-489b-982a-cf73908c16d2
Skylaris, Chris-Kriton
8f593d13-3ace-4558-ba08-04e48211af61
Wang, Ling
c50767b1-7474-4094-9b06-4fe64e9fe362

Mathas, Dimitrios (2022) Simulation of Lubricant Properties and their Interactions with Surfaces. University of Southampton, Doctoral Thesis, 158pp.

Record type: Thesis (Doctoral)

Abstract

The behaviour of lubricants at operational conditions, such as at high temperatures and pressures, is a topic of great industrial interest. In particular, viscosity and the viscosity-pressure relation are especially important for applications and their determination by computational simulations is very desirable.

In this thesis we evaluate methods to compute these quantities based on fully atomistic molecular dynamics simulations which are computationally demanding but also have the potential to be most accurate. We tested several molecules that are used as lubricants, such as 9,10-dimethyloctadecane, main component of PAO-2 base oil, which was used as the main lubricant for our simulations. The methods used for the viscosity simulations are the Green-Kubo equilibrium molecular dynamics (EMD-GK), the direct computation of viscosity from shear during non-equilibrium MD (NEMD) and the use of confined NEMD, where the fluid is confined within explicitly defined iron oxide wall surfaces, at pressures of up to 1.0 GPa and various temperatures (40-150 degrees Celsius). We present the theory behind these methods and investigate how the simulation parameters affect the results obtained, to ensure viscosity convergence with respect to the simulation intervals and all other parameters. We show that by using each method in its regime of applicability, we can achieve good agreement between simulated and measured values. NEMD simulations at high pressures captured zero shear viscosity successfully, while at 40 degrees Celsius EMD-GK is only applicable to pressures up to 0.3 GPa, where the viscosity is lower. In NEMD, longer and multiply repeated simulations reduce the standard deviation of viscosity, which is essential at lower pressures.

Additionally, by using confined NEMD simulations, it was demonstrated that the film thickness of the fluid affects viscosity, and as we increase the number of lubricant molecules, we approach the viscosity value of the bulk fluid derived from NEMD simulations.

Another aspect of these methods is the choice of the utilised force field for the atomic interactions. This was investigated by selecting three different commonly used force fields. We have explored several methods for calculating viscosity and we obtained results of particular industrial interest.

Text
PhD_thesis_mathas_for_award - Version of Record
Available under License University of Southampton Thesis Licence.
Download (41MB)
Text
Final-thesis-submission-Examination-Mr-Dimitrios-Mathas
Restricted to Repository staff only

More information

Published date: November 2022

Identifiers

Local EPrints ID: 471557
URI: http://eprints.soton.ac.uk/id/eprint/471557
PURE UUID: 86d53aa7-7192-4d5e-a681-1848d3fa183d
ORCID for Dimitrios Mathas: ORCID iD orcid.org/0000-0003-0116-0156
ORCID for Chris-Kriton Skylaris: ORCID iD orcid.org/0000-0003-0258-3433
ORCID for Ling Wang: ORCID iD orcid.org/0000-0002-2894-6784

Catalogue record

Date deposited: 11 Nov 2022 17:33
Last modified: 17 Mar 2024 03:07

Export record

Contributors

Author: Dimitrios Mathas ORCID iD
Thesis advisor: Chris-Kriton Skylaris ORCID iD
Thesis advisor: Ling Wang ORCID iD

Download statistics

Downloads from ePrints over the past year. Other digital versions may also be available to download e.g. from the publisher's website.

View more statistics

Atom RSS 1.0 RSS 2.0

Contact ePrints Soton: eprints@soton.ac.uk

ePrints Soton supports OAI 2.0 with a base URL of http://eprints.soton.ac.uk/cgi/oai2

This repository has been built using EPrints software, developed at the University of Southampton, but available to everyone to use.

We use cookies to ensure that we give you the best experience on our website. If you continue without changing your settings, we will assume that you are happy to receive cookies on the University of Southampton website.

×