EA
E.I. Assaf Martinez-Streignard
info
Please Note
<p>This page displays the records of the person named above and is not linked to a unique person identifier. This record may need to be merged to a profile.</p>
13 records found
1
Heavy petroleum streams such as heavy crudes, residues, and bitumen support modern energy and infrastructure, yet their properties vary strongly with composition, temperature, and chemical history. Because these materials consist of chemically diverse and evolving molecular populations rather than well-defined compounds, conventional structure–property correlations often lack robustness.
This dissertation develops a composition-based structure–property framework that combines Molecular Dynamics simulations with neural-network surrogate models to predict thermophysical and mechanical properties directly from molecular descriptors and temperature.
By enabling near-instantaneous property screening, comparison, and optimization of candidate formulations, the resulting draw-and-predict workflow enhances engineering design, accelerating the development of advanced hydrocarbon materials for next-generation energy and infrastructure systems. ...
This dissertation develops a composition-based structure–property framework that combines Molecular Dynamics simulations with neural-network surrogate models to predict thermophysical and mechanical properties directly from molecular descriptors and temperature.
By enabling near-instantaneous property screening, comparison, and optimization of candidate formulations, the resulting draw-and-predict workflow enhances engineering design, accelerating the development of advanced hydrocarbon materials for next-generation energy and infrastructure systems. ...
Heavy petroleum streams such as heavy crudes, residues, and bitumen support modern energy and infrastructure, yet their properties vary strongly with composition, temperature, and chemical history. Because these materials consist of chemically diverse and evolving molecular populations rather than well-defined compounds, conventional structure–property correlations often lack robustness.
This dissertation develops a composition-based structure–property framework that combines Molecular Dynamics simulations with neural-network surrogate models to predict thermophysical and mechanical properties directly from molecular descriptors and temperature.
By enabling near-instantaneous property screening, comparison, and optimization of candidate formulations, the resulting draw-and-predict workflow enhances engineering design, accelerating the development of advanced hydrocarbon materials for next-generation energy and infrastructure systems.
This dissertation develops a composition-based structure–property framework that combines Molecular Dynamics simulations with neural-network surrogate models to predict thermophysical and mechanical properties directly from molecular descriptors and temperature.
By enabling near-instantaneous property screening, comparison, and optimization of candidate formulations, the resulting draw-and-predict workflow enhances engineering design, accelerating the development of advanced hydrocarbon materials for next-generation energy and infrastructure systems.
Scymol
A python-based software package for initializing and running molecular dynamics simulations using LAMMPS
Scymol is a Python-based software package specifically designed to facilitate the setup and execution of molecular simulations in LAMMPS. It comes equipped with a user-friendly interface, which simplifies the process of initializing molecular systems and defining simulation parameters. Moreover, the software generates and executes LAMMPS simulation sequences, enabling researchers to establish comprehensive simulation schemes, such as heating or deformation cycles, in a single run. Through its successful application in diverse research projects and its modular design, Scymol demonstrates considerable promise as an indispensable tool for researchers aiming to carry out molecular dynamics simulations without sacrificing complexity or high-throughput capabilities in their methodologies.
...
Scymol is a Python-based software package specifically designed to facilitate the setup and execution of molecular simulations in LAMMPS. It comes equipped with a user-friendly interface, which simplifies the process of initializing molecular systems and defining simulation parameters. Moreover, the software generates and executes LAMMPS simulation sequences, enabling researchers to establish comprehensive simulation schemes, such as heating or deformation cycles, in a single run. Through its successful application in diverse research projects and its modular design, Scymol demonstrates considerable promise as an indispensable tool for researchers aiming to carry out molecular dynamics simulations without sacrificing complexity or high-throughput capabilities in their methodologies.
The chemo-mechanical properties of bitumen undergo significant alternations during aging and rejuvenation, posing challenges for accurately evaluating and enhancing rejuvenation efficiency in asphalt recycling. This study investigates how bitumen source, aging degree, rejuvenator type and dosage influence the chemical and rheological performance of rejuvenated bitumen. Comprehensive characterizations are performed using saturate, aromatic, resin, and asphaltene (SARA) fractionation, elemental analysis, gel permeation chromatography (GPC), and dynamic shear rheometer (DSR) tests. To elucidate chemo-rheological correlations, statistical techniques (Pearson correlation, analysis of variance (ANOVA), and Chi-square tests) are combined with artificial neural networks (ANN). Results indicate that the NB bitumen with more colloidal stability and less sulfur content exhibits the highest resistance to long-term aging. FB bitumen with 4.3 % sulfur achieves the best high-temperature deformation resistance with rutting failure temperature (RFT) higher than 80 °C, and TB bitumen exhibits the longest fatigue life. Rejuvenation using bio-oil is most effective on reducing relaxation time by up to 60 % and increasing creep compliance (Jnr3.2) by 1.7–2.5 times, depending on bitumen type. Rejuvenator dosage sensitivity for relaxation stress follows the trend: bio-oil < engine-oil < naphthenic-oil, while aromatic-oil shows variability depending on its source. Among the tested rejuvenators, bio-oil proves most effective, particularly for rejuvenating TB and FB bitumen. The ANN model demonstrates strong predictive performance for rheological properties, achieving R2 values between 0.90 and 0.98, with the highest accuracy observed for relaxation indices, followed by fatigue and rutting properties.
...
The chemo-mechanical properties of bitumen undergo significant alternations during aging and rejuvenation, posing challenges for accurately evaluating and enhancing rejuvenation efficiency in asphalt recycling. This study investigates how bitumen source, aging degree, rejuvenator type and dosage influence the chemical and rheological performance of rejuvenated bitumen. Comprehensive characterizations are performed using saturate, aromatic, resin, and asphaltene (SARA) fractionation, elemental analysis, gel permeation chromatography (GPC), and dynamic shear rheometer (DSR) tests. To elucidate chemo-rheological correlations, statistical techniques (Pearson correlation, analysis of variance (ANOVA), and Chi-square tests) are combined with artificial neural networks (ANN). Results indicate that the NB bitumen with more colloidal stability and less sulfur content exhibits the highest resistance to long-term aging. FB bitumen with 4.3 % sulfur achieves the best high-temperature deformation resistance with rutting failure temperature (RFT) higher than 80 °C, and TB bitumen exhibits the longest fatigue life. Rejuvenation using bio-oil is most effective on reducing relaxation time by up to 60 % and increasing creep compliance (Jnr3.2) by 1.7–2.5 times, depending on bitumen type. Rejuvenator dosage sensitivity for relaxation stress follows the trend: bio-oil < engine-oil < naphthenic-oil, while aromatic-oil shows variability depending on its source. Among the tested rejuvenators, bio-oil proves most effective, particularly for rejuvenating TB and FB bitumen. The ANN model demonstrates strong predictive performance for rheological properties, achieving R2 values between 0.90 and 0.98, with the highest accuracy observed for relaxation indices, followed by fatigue and rutting properties.
Previous work demonstrated that Random Forest Regressors (RFRs) could estimate the physical properties of bitumen using molecular descriptors derived from Molecular Dynamics (MD) simulations, thereby reducing the need for computationally intensive simulations. However, due to their decision-tree structure, RFRs lack true predictive capabilities, particularly for interpolation and extrapolation beyond the training data.
This study advances that foundation by employing Artificial Neural Networks (ANNs), which—when properly trained—can capture complex relationships with greater continuity and generalizability. Beyond simply replacing RFRs, we develop a fully automated framework for constructing Machine Learning Models (MLMs) to predict density and thermal expansion coefficients of bitumen. Using Optuna for hyperparameter optimization, we ensure that the information extracted from MD simulations is utilized effectively.
The resulting ANN models accurately reproduce MD-predicted densities, achieving R2>0.99, MSEs below 0.1 %, and maximum absolute errors below 5 % on test data. In addition to reducing computational cost, the models exhibit improved interpolation and extrapolation capabilities, enabling reliable predictions for properties, ranges, and compositions not explicitly simulated.
Key aspects of our approach include:
• Transitioning from RFRs to ANNs, improving generalization, interpolation, and predictive accuracy.
• Automated hyperparameter optimization, leveraging Optuna to maximize model efficiency.
• Expanding applicability, enabling property prediction for unseen compositions without additional MD simulations. ...
This study advances that foundation by employing Artificial Neural Networks (ANNs), which—when properly trained—can capture complex relationships with greater continuity and generalizability. Beyond simply replacing RFRs, we develop a fully automated framework for constructing Machine Learning Models (MLMs) to predict density and thermal expansion coefficients of bitumen. Using Optuna for hyperparameter optimization, we ensure that the information extracted from MD simulations is utilized effectively.
The resulting ANN models accurately reproduce MD-predicted densities, achieving R2>0.99, MSEs below 0.1 %, and maximum absolute errors below 5 % on test data. In addition to reducing computational cost, the models exhibit improved interpolation and extrapolation capabilities, enabling reliable predictions for properties, ranges, and compositions not explicitly simulated.
Key aspects of our approach include:
• Transitioning from RFRs to ANNs, improving generalization, interpolation, and predictive accuracy.
• Automated hyperparameter optimization, leveraging Optuna to maximize model efficiency.
• Expanding applicability, enabling property prediction for unseen compositions without additional MD simulations. ...
Previous work demonstrated that Random Forest Regressors (RFRs) could estimate the physical properties of bitumen using molecular descriptors derived from Molecular Dynamics (MD) simulations, thereby reducing the need for computationally intensive simulations. However, due to their decision-tree structure, RFRs lack true predictive capabilities, particularly for interpolation and extrapolation beyond the training data.
This study advances that foundation by employing Artificial Neural Networks (ANNs), which—when properly trained—can capture complex relationships with greater continuity and generalizability. Beyond simply replacing RFRs, we develop a fully automated framework for constructing Machine Learning Models (MLMs) to predict density and thermal expansion coefficients of bitumen. Using Optuna for hyperparameter optimization, we ensure that the information extracted from MD simulations is utilized effectively.
The resulting ANN models accurately reproduce MD-predicted densities, achieving R2>0.99, MSEs below 0.1 %, and maximum absolute errors below 5 % on test data. In addition to reducing computational cost, the models exhibit improved interpolation and extrapolation capabilities, enabling reliable predictions for properties, ranges, and compositions not explicitly simulated.
Key aspects of our approach include:
• Transitioning from RFRs to ANNs, improving generalization, interpolation, and predictive accuracy.
• Automated hyperparameter optimization, leveraging Optuna to maximize model efficiency.
• Expanding applicability, enabling property prediction for unseen compositions without additional MD simulations.
This study advances that foundation by employing Artificial Neural Networks (ANNs), which—when properly trained—can capture complex relationships with greater continuity and generalizability. Beyond simply replacing RFRs, we develop a fully automated framework for constructing Machine Learning Models (MLMs) to predict density and thermal expansion coefficients of bitumen. Using Optuna for hyperparameter optimization, we ensure that the information extracted from MD simulations is utilized effectively.
The resulting ANN models accurately reproduce MD-predicted densities, achieving R2>0.99, MSEs below 0.1 %, and maximum absolute errors below 5 % on test data. In addition to reducing computational cost, the models exhibit improved interpolation and extrapolation capabilities, enabling reliable predictions for properties, ranges, and compositions not explicitly simulated.
Key aspects of our approach include:
• Transitioning from RFRs to ANNs, improving generalization, interpolation, and predictive accuracy.
• Automated hyperparameter optimization, leveraging Optuna to maximize model efficiency.
• Expanding applicability, enabling property prediction for unseen compositions without additional MD simulations.
PDB2DAT
Automating LAMMPS data file generation from PDB molecular systems using Python, Rdkit, and Pysimm
Pdb2dat, developed in Python, is an open-source, self-contained utility that facilitates the conversion of PDB files into LAMMPS data files, catering to the need of initializing atomistic simulation from initial atomic configurations. It extracts molecular details from PDB files, uses Rdkit and Xyz2mol for bonding analysis and 3D conformer generation, and uses Pysimm for assigning force field types and charges. Designed to be lightweight and fully Pythonic, pdb2dat is suitable for use in privilege-limited high-throughput environments. The output details system topologies for use in MD simulations, significantly simplifying the preparatory steps needed by researchers to explore materials phenomena through LAMMPS.
...
Pdb2dat, developed in Python, is an open-source, self-contained utility that facilitates the conversion of PDB files into LAMMPS data files, catering to the need of initializing atomistic simulation from initial atomic configurations. It extracts molecular details from PDB files, uses Rdkit and Xyz2mol for bonding analysis and 3D conformer generation, and uses Pysimm for assigning force field types and charges. Designed to be lightweight and fully Pythonic, pdb2dat is suitable for use in privilege-limited high-throughput environments. The output details system topologies for use in MD simulations, significantly simplifying the preparatory steps needed by researchers to explore materials phenomena through LAMMPS.
This study enhances the molecular analysis of bitumen by transitioning from traditional chemical descriptors, such as SARA (Saturates, Aromatics, Resins, and Asphaltenes) fractions and elemental compositions, to specific force field atom types in Molecular Dynamics (MD) models. This shift improves the precision in predicting material properties critical for bituminous material characterization. Machine Learning Models (MLMs) were developed to use these atom types as input features, inherently reflecting fundamental chemical characteristics. Trained on data from over 1,770 LAMMPS simulations of diverse bitumen types and conditions, these MLMs enable the prediction of properties like density, heat capacity, solubility parameters, and thermal expansion coefficients without the need for additional MD simulations. The models utilize 30 chemical descriptors corresponding to specific atom types in the PCFF force field, which collectively account for over 95% of the influence on these properties. By accurately predicting fundamental, thermodynamic, and kinetic properties, the use of MLMs and force field atom types allows researchers to efficiently tweak the chemical nature of organic molecules and mixtures to achieve desired properties. With near-instantaneous prediction times, these MLMs offer valuable insights for advancing bitumen research in the construction and petroleum industries, reducing the need for more intensive simulation techniques.
...
This study enhances the molecular analysis of bitumen by transitioning from traditional chemical descriptors, such as SARA (Saturates, Aromatics, Resins, and Asphaltenes) fractions and elemental compositions, to specific force field atom types in Molecular Dynamics (MD) models. This shift improves the precision in predicting material properties critical for bituminous material characterization. Machine Learning Models (MLMs) were developed to use these atom types as input features, inherently reflecting fundamental chemical characteristics. Trained on data from over 1,770 LAMMPS simulations of diverse bitumen types and conditions, these MLMs enable the prediction of properties like density, heat capacity, solubility parameters, and thermal expansion coefficients without the need for additional MD simulations. The models utilize 30 chemical descriptors corresponding to specific atom types in the PCFF force field, which collectively account for over 95% of the influence on these properties. By accurately predicting fundamental, thermodynamic, and kinetic properties, the use of MLMs and force field atom types allows researchers to efficiently tweak the chemical nature of organic molecules and mixtures to achieve desired properties. With near-instantaneous prediction times, these MLMs offer valuable insights for advancing bitumen research in the construction and petroleum industries, reducing the need for more intensive simulation techniques.
SMI2PDB
A self-contained Python tool to generate atomistic systems of organic molecules using their SMILES notations
The advent of computational techniques, particularly atomistic simulations, has lessened the dependency on physical experiments in various scientific fields. Yet, the preparation complexity for simulations using platforms like LAMMPS and GROMACS persists. We introduce SMI2PDB, a Python tool that automates molecular systems assembly from SMILES to PDB format, easing molecular dynamics simulation setups. SMI2PDB manages molecule configuration and quantification effortlessly, establishes stable conformers, applies random rotations, and positions them in a simulation box with a Sobol sequence to reduce overlaps. This script facilitates the rapid preparation of complex organic mixtures for use in simulations, enhancing the exploration of novel materials.
...
The advent of computational techniques, particularly atomistic simulations, has lessened the dependency on physical experiments in various scientific fields. Yet, the preparation complexity for simulations using platforms like LAMMPS and GROMACS persists. We introduce SMI2PDB, a Python tool that automates molecular systems assembly from SMILES to PDB format, easing molecular dynamics simulation setups. SMI2PDB manages molecule configuration and quantification effortlessly, establishes stable conformers, applies random rotations, and positions them in a simulation box with a Sobol sequence to reduce overlaps. This script facilitates the rapid preparation of complex organic mixtures for use in simulations, enhancing the exploration of novel materials.
AA2UA
Converting all-atom models into their united atom coarse grained counterparts for use in LAMMPS
Atomistic simulations are crucial for understanding material properties at the molecular level but are limited by high computational costs, especially for large, complex systems like bituminous materials. Our team developed a Force-matched United Atom (UA) Coarse Graining (CG) force field to enhance computational efficiency while retaining atomic detail. However, converting all-atom models to CG models is complex, requiring detailed atom-to-bead mapping and compatibility with molecular dynamics (MD) engines like LAMMPS. To address this, we introduce AA2UA, an open-source software that simplifies the conversion of PDB files into LAMMPS-readable structure topology files, facilitating broader use of the developed UA force field.
...
Atomistic simulations are crucial for understanding material properties at the molecular level but are limited by high computational costs, especially for large, complex systems like bituminous materials. Our team developed a Force-matched United Atom (UA) Coarse Graining (CG) force field to enhance computational efficiency while retaining atomic detail. However, converting all-atom models to CG models is complex, requiring detailed atom-to-bead mapping and compatibility with molecular dynamics (MD) engines like LAMMPS. To address this, we introduce AA2UA, an open-source software that simplifies the conversion of PDB files into LAMMPS-readable structure topology files, facilitating broader use of the developed UA force field.
This paper presents a United Atom (UA) force field for simulating hydrocarbon molecules in bituminous materials, integrating explicit hydrogens into beads with their parent atom. This method simplifies all-atom molecular models, significantly accelerating Molecular Dynamics (MD) simulations of bitumen by 10 to 100 times. Key advantages include halving the particle count, eliminating complex hydrogen interactions, and decreasing the degrees of freedom of the molecules. Developed by mapping forces from an all-atom model to the centers of mass of UA model beads, the force field ensures accurate replication of energies, forces, and molecular conformations, mirroring properties like pressure and density. It features 17 bead types and 287 interaction types, encompassing various hydrocarbon molecules. The UA force field's stability, surpassing all-atom models, is a notable achievement. This stability, stemming from smoother potential energy surfaces, leads to consistent property measurements and improved stress tensor accuracy. It enables the extension of MD simulations to larger spatiotemporal scales, crucial for understanding complex phenomena such as phase separation in bituminous materials. This foundational work sets the stage for future developments, including refining parameters and introducing new bead types, to enhance the modeling capabilities of the force field, thereby advancing the application and understanding of bituminous materials.
...
This paper presents a United Atom (UA) force field for simulating hydrocarbon molecules in bituminous materials, integrating explicit hydrogens into beads with their parent atom. This method simplifies all-atom molecular models, significantly accelerating Molecular Dynamics (MD) simulations of bitumen by 10 to 100 times. Key advantages include halving the particle count, eliminating complex hydrogen interactions, and decreasing the degrees of freedom of the molecules. Developed by mapping forces from an all-atom model to the centers of mass of UA model beads, the force field ensures accurate replication of energies, forces, and molecular conformations, mirroring properties like pressure and density. It features 17 bead types and 287 interaction types, encompassing various hydrocarbon molecules. The UA force field's stability, surpassing all-atom models, is a notable achievement. This stability, stemming from smoother potential energy surfaces, leads to consistent property measurements and improved stress tensor accuracy. It enables the extension of MD simulations to larger spatiotemporal scales, crucial for understanding complex phenomena such as phase separation in bituminous materials. This foundational work sets the stage for future developments, including refining parameters and introducing new bead types, to enhance the modeling capabilities of the force field, thereby advancing the application and understanding of bituminous materials.
This study explores the use of chemical descriptors derived from force field atom types to predict Fickian diffusion coefficients of rejuvenators in bitumen, utilizing machine learning models trained on data from 240 non-equilibrium molecular dynamics simulations. The simulations cover three bitumen types (NO, TO, FO), five aging degrees, and four temperatures (60 °C, 120 °C, 160 °C, 200 °C), capturing diffusion coefficients ranging from 0.0068e-10 m2/s in highly aged bitumens at 60 °C to 4.35e-10 m2/s in fresher samples at 200 °C. The MLM, built with 18 chemical descriptors for bitumen and rejuvenator sides, achieves an R2 of 0.97, accurately predicting diffusion across varied conditions. This approach abstracts away from the need for repeated MD simulations, enabling diffusion predictions even for systems outside the original dataset. The manuscript presents three case studies to illustrate how the model can be used for the iterative design of rejuvenators by optimizing molecular structures based on critical chemical features, such as rejuvenator oxygen content, bitumen sulfur content, and molecular weights. It also demonstrates how the model offers a practical framework for understanding the diffusion and performance of rejuvenators by linking time-dependent factors—such as concentration, depth, and rejuvenation time—with the bulk properties of bitumen-rejuvenator systems, facilitating industrial applications.
...
This study explores the use of chemical descriptors derived from force field atom types to predict Fickian diffusion coefficients of rejuvenators in bitumen, utilizing machine learning models trained on data from 240 non-equilibrium molecular dynamics simulations. The simulations cover three bitumen types (NO, TO, FO), five aging degrees, and four temperatures (60 °C, 120 °C, 160 °C, 200 °C), capturing diffusion coefficients ranging from 0.0068e-10 m2/s in highly aged bitumens at 60 °C to 4.35e-10 m2/s in fresher samples at 200 °C. The MLM, built with 18 chemical descriptors for bitumen and rejuvenator sides, achieves an R2 of 0.97, accurately predicting diffusion across varied conditions. This approach abstracts away from the need for repeated MD simulations, enabling diffusion predictions even for systems outside the original dataset. The manuscript presents three case studies to illustrate how the model can be used for the iterative design of rejuvenators by optimizing molecular structures based on critical chemical features, such as rejuvenator oxygen content, bitumen sulfur content, and molecular weights. It also demonstrates how the model offers a practical framework for understanding the diffusion and performance of rejuvenators by linking time-dependent factors—such as concentration, depth, and rejuvenation time—with the bulk properties of bitumen-rejuvenator systems, facilitating industrial applications.
This study employs strain-controlled oscillatory deformations in Molecular Dynamics (MD) simulations to evaluate the dynamic properties of all-atom molecular systems, specifically targeting the SARA fractions of bitumen. Twelve molecular systems representing these fractions were modeled using the PCFF force field. The simulations effectively captured their viscoelastic properties across multiple frequency domains, including Elastic, Glassy, Rubbery, and Viscous responses. Reported storage and loss moduli range from thousands to tens of megapascals, with viscosities from tens to near-zero Pascal-seconds across various frequencies and temperatures, aligning well with experimental observations. Saturates and Aromatics were identified as the softest and most thermally susceptible fractions, while Resins and Asphaltenes were the stiffest and least susceptible. The study reveals that the relaxation time of all-atom molecular systems is significantly shorter than in experimental setups, necessitating careful comparison of stress-related phenomena across equivalent relaxation times. Although this allows for the exploration of response profiles in computationally tractable simulations, the nature of all-atom force fields and simulation algorithms introduces spatiotemporal scale discrepancies that must be addressed in future simulations involving the study of stress-related phenomena using MD.
...
This study employs strain-controlled oscillatory deformations in Molecular Dynamics (MD) simulations to evaluate the dynamic properties of all-atom molecular systems, specifically targeting the SARA fractions of bitumen. Twelve molecular systems representing these fractions were modeled using the PCFF force field. The simulations effectively captured their viscoelastic properties across multiple frequency domains, including Elastic, Glassy, Rubbery, and Viscous responses. Reported storage and loss moduli range from thousands to tens of megapascals, with viscosities from tens to near-zero Pascal-seconds across various frequencies and temperatures, aligning well with experimental observations. Saturates and Aromatics were identified as the softest and most thermally susceptible fractions, while Resins and Asphaltenes were the stiffest and least susceptible. The study reveals that the relaxation time of all-atom molecular systems is significantly shorter than in experimental setups, necessitating careful comparison of stress-related phenomena across equivalent relaxation times. Although this allows for the exploration of response profiles in computationally tractable simulations, the nature of all-atom force fields and simulation algorithms introduces spatiotemporal scale discrepancies that must be addressed in future simulations involving the study of stress-related phenomena using MD.
Journal article
(2024)
-
Yangming Gao, Xueyan Liu, Shisong Ren, Eli I. Assaf Martinez-Streignard, Pengfei Liu, Yuqing Zhang
Bitumen fatigue resistance is critical to determine the overall fatigue performance and service life of asphalt pavements. However, the mechanisms responsible for fatigue damage of bitumen have previously not been well understood. Molecular dynamics (MD) simulation has recently emerged as a powerful computer-aided numerical technique to model the microscopic failure behaviours in materials. This study aims to use the MD method to investigate the molecular origin of bitumen fatigue damage. The molecular models of the virgin and aged PEN70/100 bitumen were firstly constructed based on their saturate, aromatic, resin and asphaltene (SARA) four fractions. An MD equilibrium was run on the developed bitumen models with the assigned interatomic potentials. Following an MD-based tensile simulation, a strain-controlled fatigue simulation was performed to study the nanostructure and damage behaviours of the virgin and aged bitumen under fatigue loading by calculating the stress-strain response, potential energy, molecular structure and nanovoid volumes. Furthermore, a rheometer measurement was also conducted to characterise the fatigue damage of the bitumen directly by a crack length at the macroscale. Results indicate that the bitumen molecules become unfolded and tend to align along the loading direction when fatigue loading was applied. The change in the molecular configuration helped the molecular chains move closer together and thus contributed to the reduction of the intermolecular interactions including the van der Waals and Coulombic energies. With the increasing load cycles, nanovoids were formed and grew in the bitumen through molecular rearrangement and movement, leading to microscopic fatigue damage of the bitumen. It was found that the aged bitumen produced more severe fatigue damage than the virgin bitumen, which was indicated by the MD-based nanovoid volume at the molecular scale and the DSR-based crack length at the macroscale. The findings from MD simulation provide a fundamental understanding of the molecular origin of fatigue damage, that cannot be experimentally detected for bitumen materials.
...
Bitumen fatigue resistance is critical to determine the overall fatigue performance and service life of asphalt pavements. However, the mechanisms responsible for fatigue damage of bitumen have previously not been well understood. Molecular dynamics (MD) simulation has recently emerged as a powerful computer-aided numerical technique to model the microscopic failure behaviours in materials. This study aims to use the MD method to investigate the molecular origin of bitumen fatigue damage. The molecular models of the virgin and aged PEN70/100 bitumen were firstly constructed based on their saturate, aromatic, resin and asphaltene (SARA) four fractions. An MD equilibrium was run on the developed bitumen models with the assigned interatomic potentials. Following an MD-based tensile simulation, a strain-controlled fatigue simulation was performed to study the nanostructure and damage behaviours of the virgin and aged bitumen under fatigue loading by calculating the stress-strain response, potential energy, molecular structure and nanovoid volumes. Furthermore, a rheometer measurement was also conducted to characterise the fatigue damage of the bitumen directly by a crack length at the macroscale. Results indicate that the bitumen molecules become unfolded and tend to align along the loading direction when fatigue loading was applied. The change in the molecular configuration helped the molecular chains move closer together and thus contributed to the reduction of the intermolecular interactions including the van der Waals and Coulombic energies. With the increasing load cycles, nanovoids were formed and grew in the bitumen through molecular rearrangement and movement, leading to microscopic fatigue damage of the bitumen. It was found that the aged bitumen produced more severe fatigue damage than the virgin bitumen, which was indicated by the MD-based nanovoid volume at the molecular scale and the DSR-based crack length at the macroscale. The findings from MD simulation provide a fundamental understanding of the molecular origin of fatigue damage, that cannot be experimentally detected for bitumen materials.
Journal article
(2023)
-
Eli I. Assaf, Xueyan Liu, Peng Lin, Sandra Erkens, Sayeda Nahar, Liz I.S. Mensink
Conventional Molecular Dynamics (MD) models of bitumen are built by homogeneously mixing molecules in a volume without considering that the molecules in bitumen are known to exhibit phase behavior and form distinctive molecular arrangements. These are known to have a significant impact in the behavior of bitumen, and considering their existence is paramount in producing improved representations of bitumen using computational models. This study explores whether MD models of bitumen that are conventionally assumed to be in equilibrium can still undergo significant phase separation over considerably long simulation times. It also aims to establish a more formal pathway to build and study models with highly heterogeneous arrangements of their molecules. Moreover, it aims to evaluate whether the presence of distinct morphologies have a significant impact in numerous physical properties of bitumen. The study shows that conventional and widely used models of bitumen exhibit significant molecular rearrangements over long times (>360 ns). It also shows that building heterogeneous morphologies is possible and result in energetically favorable conformations. Moreover, it proves that studying properties regularly used to validate MD models of bitumen (e.g., density) are insufficient in assessing the impact of different morphologies; more thorough methods are required to evaluate them.
...
Conventional Molecular Dynamics (MD) models of bitumen are built by homogeneously mixing molecules in a volume without considering that the molecules in bitumen are known to exhibit phase behavior and form distinctive molecular arrangements. These are known to have a significant impact in the behavior of bitumen, and considering their existence is paramount in producing improved representations of bitumen using computational models. This study explores whether MD models of bitumen that are conventionally assumed to be in equilibrium can still undergo significant phase separation over considerably long simulation times. It also aims to establish a more formal pathway to build and study models with highly heterogeneous arrangements of their molecules. Moreover, it aims to evaluate whether the presence of distinct morphologies have a significant impact in numerous physical properties of bitumen. The study shows that conventional and widely used models of bitumen exhibit significant molecular rearrangements over long times (>360 ns). It also shows that building heterogeneous morphologies is possible and result in energetically favorable conformations. Moreover, it proves that studying properties regularly used to validate MD models of bitumen (e.g., density) are insufficient in assessing the impact of different morphologies; more thorough methods are required to evaluate them.