Selecting the Optimal Calculation Method and Chemical Reagents in Surface Energy Tests of Asphalt Materials

TL;DR

This paper improves the accuracy of surface energy tests for asphalt by optimizing calculation methods and reagent combinations.

Contribution

The study introduces total least squares and a systematic reagent selection method to enhance surface energy test reliability.

Findings

Total least squares reduces fitting errors compared to classical least squares in calculating asphalt surface energy parameters.

The reagent combination WFSD (distilled water + formamide + dimethyl sulfoxide + diiodomethane) yields the most accurate and stable results.

WFSD reduces error rates by 17.71% and 64.80% compared to WFEG for two types of asphalt.

Abstract

In surface energy tests of asphalt materials, the inaccuracy of the calculation method (e.g., least squares (LS)) and the arbitrary selection of chemical reagent combinations lead to unstable results, threatening the quantitative evaluation of asphalt–aggregate adhesion durability. This study addresses these two scientific deficiencies with the following findings: (1) when simultaneous equations are used to calculate the asphalt surface energy parameters, the total least squares method should be used instead of the classical least squares method to reduce the fitting error; (2) the selection of the reagent combination should be based on which one is the most rational in terms of the physical characterization, leap degree, abnormal values, and other requirements, and the reagent combination with the fewest abnormal values should be chosen as the best scheme. The results show that (1) compared with the classical least squares method, the total least squares method reduces the fitting error between the calculated and real values of asphalt surface energy parameters and improves the accuracy and stability of the calculation results; (2) the best reagent combination scheme is WFSD (distilled water + formamide + dimethyl sulfoxide + diiodomethane). The calculated values of asphalt surface energy parameters were more accurate and reasonable, and the calculation results had no abnormal values. Compared with WFEG (distilled water + formamide + ethylene glycol + glycerol), the error rate of the reagent combination scheme WFSD in calculating the total surface energy of two kinds of asphalt was reduced by 17.71% and 64.80%, respectively. These findings establish a reliable framework for the accurate quantification of surface energy, addressing the critical issue of reagent-dependent variability in the results and strengthening the scientific basis for evaluating the durability of asphalt pavement.

Figures

Full text

1. Introduction

In the design of surface asphalt mixtures, asphalt–aggregate adhesion is a key concern, directly impacting road performance parameters such as fatigue life, self-healing ability, and water stability. Current domestic/foreign specifications commonly use the boiling test to evaluate adhesion; however, this test only provides qualitative grades [1] without quantitative results. The surface energy method determines the binding energy of interface adhesion as a quantitative index, which is closely linked to interface cracking [2,3], and is widely applied to study mixture aging, self-healing, water damage, and other factors during service [4,5,6,7].

The primary prerequisite for the quantitative evaluation of the adhesion mentioned above is to accurately measure the surface energy parameters of the asphalt and aggregate. For this reason, the Good–van Oss–Chaudhury (GvOC) theoretical system has been widely adopted in the road industry; in this system, the asphalt and aggregate have three basic surface energy parameters—the nonpolar component $[eqn]$ , the polar acid component $[eqn]$ , and the polar base component $[eqn]$ —which together constitute the total surface energy γ [8,9,10,11]. For asphalt materials, the usual test methods for measuring surface energy parameters are the Wilhelmy plate method and the sessile drop method. Compared to the sessile drop method, the Wilhelmy plate method provides more stable contact angle data [12,13]. Thus, the Wilhelmy plate method is the preferred method for testing asphalt surface energy and will be applied later in this study.

To evaluate the asphalt surface energy using the Wilhelmy plate method, it is necessary to select three or more kinds of chemical reagents with known surface energy parameters as a combination of reagents and measure the contact angles θ between various chemical reagents and asphalts using tests, as shown in Figure 1.

Next, statically indeterminate equations should be established to determine the surface energy parameters of asphalt using the least squares (LS) method [14]. The widespread prevalence of LS can be attributed to its straightforward computational framework and well-established theoretical underpinnings. As a cornerstone of linear regression analysis, LS finds extensive application in diverse fields, such as data modeling and trend forecasting. In the context of recent research on calculating asphalt surface energy, for instance, LS continues to be the method of choice, as evidenced by previous studies [3,4,15].

The surface energy parameters of known chemical reagents must be determined before solving the equations. In most studies, the surface energy parameters of various chemical reagents listed in the literature are often directly selected. Through the analysis of extensive data and long-term experimental research, it has been found that the surface energy parameters of the same asphalt measured with different reagent combinations are obviously different, and for some reagent combinations, the surface energy parameters cannot even be calculated [11,16,17]. The negative square root often appears in the solution process, which is contrary to the physical implication of asphalt surface energy. Some researchers believe that the negative square root is caused by ill-conditioned equations and try to solve it using singular value decomposition [18], but this problem persists; others directly set the surface energy parameter value of the negative square root to zero [19], which is inconsistent with the basic properties of asphalt.

At the same time, many studies show that the asphalt surface energy parameters calculated with different reagent combinations are considerably different. When selecting reagent combinations, some scholars treat the surface energy parameters of reagent combinations as a coefficient matrix and then calculate the condition number to compare the stability of asphalt surface energy parameters calculated with different reagent combinations [20,21]. However, when too many kinds of reagents are involved, a feasible reagent selection scheme cannot be determined by the condition number. Some scholars plot the curve of the relationship between the total surface energy $[eqn]$ of the chemical reagents and the value of $[eqn]$ and select reagent combinations according to the smoothness of the fitting line [22,23]. Nonetheless, even using this selection method, the calculated values of asphalt surface energy parameters still have large differences.

Therefore, based on the two perspectives of optimizing the equation calculation method and reagent combination, the main reasons for the differences between asphalt surface energy parameters measured using the Wilhelmy plate method can be summarized as follows:

• (1)When the LS method is used to solve asphalt surface energy parameters, the optimal fitting method is not used, which causes a large fitting error between the calculated and real values of asphalt surface energy parameters;
• (2)There is a lack of a reasonable and effective evaluation method for selecting different reagent combinations.

Therefore, although the surface energy is widely used to evaluate the adhesion of asphalt, problems remain in terms of reagent selection and the calculation method for determining asphalt surface energy, and there is still a lack of research in this area. In order to solve these problems, this study focused on the optimization of the calculation method and the selection of chemical reagents in surface energy tests of road asphalt materials. First, eight chemical reagents were selected to test their contact angles with asphalts using the Wilhelmy plate method. Second, by analyzing the basic principles and geometric interpretation of the LS method in the solution process, the main reasons for the errors were explored, and the total least squares (TLS) method was proposed to calculate asphalt surface energy parameters. Third, the best reagent combination was selected by using the leap degree method. Fourth and last, the main research results are summarized. The methodology flowchart is shown in Figure 2.

2. Contact Angle Test

2.1. Reagent Selection

In this study, eight chemical reagents were preliminarily selected: distilled water (W), formamide (F), ethylene glycol (E), glycerol (G), dimethyl sulfoxide (S), diiodomethane (D), benzyl alcohol (B), and n-octanol (N). The basic selection principles were as follows [6]:

• The chemical reagent is a single, homogeneous, and pure liquid reagent and does not dissolve or react with asphalt material.
• In order to solve asphalt surface energy parameters using simultaneous equations, the surface energy parameters of chemical reagents must be known;
• Chemical reagent droplets can form a stable contact angle on the surface of asphalt film slides; that is, the total surface energy of the chemical reagent is greater than that of asphalt material.

The above eight chemical reagents were selected because they were easily accessible in the laboratory and are consistent with the basic principles. The polarity parameters of the eight reagents cover a wide range, from highly polar to nonpolar, effectively constraining the surface energy parameters of asphalt. The abbreviations and surface energy parameters of all reagents are listed in Table 1, where W, F, E, G, S, and B are nonpolar reagents, while D and N are unipolar reagents.

The polarity parameters of the eight reagents cover a wide range, from highly polar to nonpolar (as shown in Table 1), effectively constraining the surface energy parameters of asphalt:

• Highly polar reagents (such as water and formamide): They dominate the polar interactions and sensitively reflect the polar components ( $[eqn]$ and $[eqn]$ ) of asphalt.
• Moderately polar reagents (such as ethylene glycol and glycerol): They balance the polar and nonpolar interactions and are used to adjust the condition number of the equation set.
• Weakly polar/nonpolar reagents (such as diiodomethane and n-octanol): They mainly contribute nonpolar dispersion forces and constrain the $[eqn]$ component of asphalt.

2.2. Asphalt Preparation

Two types of asphalt materials, namely, 70# matrix asphalt and SBS (I-D type)-modified asphalt, were selected for this study, and the basic characteristics of each are presented in Table 2. The test results all meet the technical requirements for asphalt pavement construction. Measurements were carried out with reference to the standards for asphalt testing, specifically JTG E20-2011 in China [24].

2.3. Measurement Method

Although the Wilhelmy plate method produces more stable results than the sessile drop method, contact angle measurements may still be affected by factors such as asphalt film flatness, reagent droplet acceleration, and instrument accuracy. Therefore, in order to reduce errors caused by the experiment, experimentation was performed in accordance with ASTM D1331-14 [25].

The equipment used in the test was an automatic surface tension meter. Five parallel tests were conducted for each chemical reagent as well as for asphalt, and the average value of three test measurements was calculated to ensure that the test error is within ±2°. The measured results are summarized in Table 3.

3. Calculation of Asphalt Surface Energy Parameters

3.1. Calculation Principle of Asphalt Surface Energy Parameters

After the contact angles between chemical reagents and asphalts are measured, the asphalt surface energy parameters can be calculated according to the Young–Dupre equation, as shown in Formula (1):

$$[eqn]$$

where $[eqn]$ , $[eqn]$ , $[eqn]$ , and $[eqn]$ are the three surface energy parameters and total surface energy of the selected chemical reagent, respectively (mJ/m^2^); $[eqn]$ , $[eqn]$ , and $[eqn]$ are the three surface energy parameters of asphalt, respectively (mJ/m^2^); and θ is the contact angle measured using the Wilhelmy plate method (°).

Next, the surface energy polar component $[eqn]$ and total surface energy $[eqn]$ of asphalt are calculated using the following formulas:

$$[eqn]$$ $$[eqn]$$

As mentioned above, it is imperative to select at least three chemical reagents with known surface energy parameters as reagent combinations to calculate the asphalt surface energy parameters according to the Young–Dupre equation. After contact angles are measured, statically indeterminate equations are established and solved, as shown in Formula (4):

$$[eqn]$$

where $[eqn]$ , $[eqn]$ , $[eqn]$ , and $[eqn]$ are the three surface energy parameters and total surface energy of the $[eqn]$ chemical reagent, respectively (mJ/m^2^), and $[eqn]$ ; $[eqn]$ is the value of the contact angle between the $[eqn]$ chemical reagent and asphalt (°); and $[eqn]$ , $[eqn]$ , and $[eqn]$ are the same as in Formula (1).

In the surface energy test of asphalt materials, Equation (4) is usually solved using the LS method, whose basic principle is to solve and optimize a linear regression equation according to the LS criterion [26,27]. In Formula (4), there are three unknown quantities: $[eqn]$ , $[eqn]$ , and $[eqn]$ . Consequently, asphalt surface energy parameters are calculated by solving ternary linear equations. The surface energy parameters of different chemical reagents are plotted in a three-dimensional coordinate system and expressed as points P_1_ ( $[eqn]$ , $[eqn]$ , $[eqn]$ ), P_2_ ( $[eqn]$ , $[eqn]$ , $[eqn]$ ), …, P_n_ ( $[eqn]$ , $[eqn]$ , $[eqn]$ ), and $[eqn]$ . The plane ABC is obtained by fitting according to LS, which minimizes the average value of the sum of distances P_n_D from all the above points, P_1_, P_2_, …, P_n_, to the plane ABC along the Z-axis, as shown in Figure 3.

Thus, the length of the fitting distance P_n_D from any point P_n_ to the plane ABC can be expressed as

$$[eqn]$$

However, in the actual process of measurement, the errors in the raw data are not limited to a single coordinate axis, so LS cannot obtain the shortest fitting distance. As a comparison, using the TLS method to solve Equation (4) can account for all the errors in the direction of all axes [28,29]. Based on the fitting principle of LS, the plane $[eqn]$ is obtained by fitting according to TLS, which minimizes the average value of the sum of the vertical distances $[eqn]$ from all the above points to the plane $[eqn]$ , as shown in Figure 4.

Likewise, the length of the fitting distance $[eqn]$ from any point P_n_ to the plane $[eqn]$ can be expressed as

$$[eqn]$$

By comparison, it can be seen that $[eqn]$ < $[eqn]$ ; that is, TLS can be used to determine the plane with smaller fitting errors to calculate more accurate asphalt surface energy parameters.

3.2. Calculation of Asphalt Surface Energy Parameters Using TLS

At least three chemical reagents are needed to calculate asphalt surface energy parameters using TLS. If three of the eight chemical reagents with known surface energy parameters are randomly selected to solve Equation (4), then there are $[eqn]$ reagent combinations; if four reagents are selected, then there are $[eqn]$ reagent combinations. Among combinations, WFE represents the reagent combination of distilled water + formamide + ethylene glycol, while WFEG represents the reagent combination of distilled water + formamide + ethylene glycol + glycerol, and so on. For each reagent combination, Microsoft Excel was used to calculate the asphalt surface energy parameters by using LS and TLS, respectively. All calculated values of the total surface energy of asphalt $[eqn]$ are summarized in Figure 5 and Figure 6.

According to the analysis of the data in the figures above for 70# matrix asphalt and SBS-modified asphalt, when the LS and TLS methods are used for calculations, the absolute values of the total average change rate of surface energy are 35.1% and 235.0%, indicating a significant distinction between the calculation results of the two methods. When the reagent combination is FSB, the maximum absolute values of the change rate are as high as 1316.2% and 23,060.0%. This shows that the data processing method has a significant impact on the results.

The surface energy of asphalt is a basic physical parameter that reflects its adhesion, and it will not change with different reagent combinations. However, in practical calculations, the total surface energy of asphalt calculated with different reagent combinations will still have large differences. Using the LS method, the coefficients of variation of total surface energy for the 70# matrix asphalt and SBS-modified asphalt are 20.32% and 26.35%, respectively; using the TLS method, the coefficients of variation of the total surface energy of the two asphalts are 15.92% and 24.65%, respectively. Compared to LS, the coefficient of variation of the calculated results of TLS is smaller, so the results calculated with the TLS method are more stable. The LS method only minimizes the single-axis (Z-axis) error, while the TLS method minimizes the full-dimensional vertical error, especially in ill-conditioned equations (for example, the change rate of the calculated values for the FSB combination is reduced from 23,060% to an acceptable range).

At the same time, if the total surface energy of the two kinds of asphalt is calculated using the TLS method, the coefficients of variation of the total surface energy of the 70# matrix asphalt with three- and four-reagent combinations are 14.22% and 12.57%, respectively, while the coefficients of variation of the total surface energy of SBS-modified asphalt are 24.77% and 14.15%, respectively. Undoubtedly, the total surface energy of asphalt calculated with four-reagent combinations is more stable. All conclusions above indicate that the selection of reagents will directly affect the accuracy and stability of results in surface energy tests of asphalt materials. Therefore, the following paragraphs concentrate on the analysis of the calculated results and the optimization of reagent combinations.

3.3. Analysis of Test Results

The analysis of surface energy tests for 70# asphalt and SBS-modified asphalt reveals significant differences between calculation methods. The total least squares (TLS) method outperforms the least squares (LS) method in reducing fitting errors and improving stability, as it considers full-dimensional errors rather than single-axis errors. The selected reagent combination also impacts the results: four-reagent combinations generally yield more stable surface energy parameter calculations than three-reagent ones, likely due to the redundancy of hyperstatic equations. This study highlights the importance of rational reagent selection (including nonpolar and polar reagents) and appropriate calculation methods to ensure accuracy in asphalt surface energy testing. Further optimization of reagent combinations and the validation of calculation methods are recommended for reliable adhesion performance evaluation.

4. Optimization of Reagent Combinations

In this study, the rationality and outliers of the calculated results were analyzed to determine the optimal reagent combination scheme.

4.1. Rationality Analysis of Calculated Results

Through observing the values of every component of asphalt surface energy calculated using the TLS method described above, it was revealed that calculations with some reagent combinations result in zero values for certain surface energy components of asphalt, as shown in Table 4.

As a complex mixture composed of various polar/nonpolar organic compounds, asphalt exhibits both dispersion forces (nonpolar interactions) and polar interactions (such as hydrogen bonding and dipole interactions) between molecules. Therefore, its surface energy must include nonzero polar and nonpolar components. It is impossible for asphalt to be a unipolar or nonpolar substance, and its three surface energy components should be greater than zero. By excluding reagent combinations with zero values in the calculation results (WFN, FED, etc.; see Table 4), the 21 effective combinations selected (Figure 7) all meet the condition of three nonzero components, and the proportion of γL in the total surface energy in the calculation results is over 80%, which conforms to the physical characteristics of asphalt with nonpolar interactions as the main component. Thus, the calculated results of the above reagent combinations violate the basic chemical properties of asphalt. Furthermore, in the process of optimization, these reagent combination schemes should be excluded first, and only the reagent combinations whose calculated surface energy components are not zero should be considered.

In addition, according to the calculation principle of the Young–Dupre equation, droplets of the selected chemical reagents should be able to form a stable contact angle on the surface of asphalt slides during the test [8]. If the calculated value of the total surface energy of asphalt is greater than that of the chemical reagent, then the reagent cannot form a stable contact angle on the surface of the asphalt slide, which is inconsistent with the calculation principle of the Young–Dupre equation. Consequently, it is necessary to exclude such reagent combination schemes. A total of 18 reagent combinations meet the above two criteria, namely, WFB, WFG, WGS, WGD, WGB, WGN, WSB, WDB, WDN, WEG, WEB, WBN, FGB, EGB, GSB, GDB, GDN, and GBN.

In addition to the above 18 reagent combinations, the remaining 21 reagent combinations and the calculated values of asphalt surface energy parameters are summarized in Figure 7.

As can be seen in Figure 7, the distribution ranges of calculated surface energy parameters of the two kinds of asphalt are as follows: 17 mJ/m^2^ ≤ $[eqn]$ ; $[eqn]$ ≤ 25 mJ/m^2^; 0 mJ/m^2^ ≤ $[eqn]$ ; and $[eqn]$ ≤ 3 mJ/m^2^. Moreover, most $[eqn]$ values are much greater than $[eqn]$ and $[eqn]$ and are close to $[eqn]$ , indicating that asphalt is a substance dominated by intermolecular nonpolar forces, which is consistent with the basic chemical properties of asphalt.

Meanwhile, comparing different reagent combinations that meet the requirements shows that some differences in the calculated values of asphalt surface energy parameters remain.

Therefore, in order to obtain more accurate and reasonable asphalt surface energy parameters, it is crucial to establish a quick and effective method for selecting reagent combinations.

4.2. Testing for Outliers in the Calculated Results

In the calculated results that meet the requirements for surface energy parameters of the same kind of asphalt, the differences between most of the data points are tiny; therefore, it can be assumed that the main factor leading to differences is the test error caused by different reagent combinations. However, a small number of data points are significantly different from other data and can be called outliers [30,31,32,33]. In order to further eliminate the interference of outliers in the raw data and ensure the accuracy and reliability of the calculated results, it was necessary to conduct an outlier test on the calculated results of the 21 reagent combinations in Figure 7.

The basic steps of the outlier test are as follows: the calculated results for surface energy parameters of a certain asphalt with 21 reagent combinations are arranged in order from small to large and set as $[eqn]$ . $[eqn]$ is a point estimate of the expectation μ that depends only on the sequence $[eqn]$ , and $[eqn]$ is called the leap degree of expectation μ at a point k, or the leap degree at point k for short. The calculation formula for the leap degree is as follows:

$$[eqn]$$ $$[eqn]$$ $$[eqn]$$

where $[eqn]$ and $[eqn]$ are both point estimates of the expectation μ; k is a sequence of calculated values for any of the asphalt surface energy parameters arranged in order, $[eqn]$ ; and $[eqn]$ is the leap degree at point k.

Thus, the existence of an outlier will certainly cause an intermittent leap in the point estimate of the expectation, and when more than one outlier exists, the maximum leap point of the point estimate of the expectation—namely, the point with the largest leap degree—is most likely to be the starting point of abnormal data [32]. After calculating the leap degree at every point using Equations (7)–(9), in order from left to right, when the first maximum value greater than one appears in a leap degree, the point corresponding to that leap degree is the end point of small outliers; that is, all the calculated values before that point (including that point) are small outliers. In order from right to left, when the first maximum value greater than one appears in a leap degree, the point corresponding to that leap degree is the starting point of large outliers; that is, all the calculated values after that point (excluding that point) are large outliers.

There are three types of cases in which a set of data has outliers—there are only small outliers; there are only large outliers; or there are both small and large outliers [29]. The reagent combinations and leap degrees corresponding to every point are provided in Appendix A Table A1 and were obtained according to the above test steps, and the numbers of outliers in the calculated results of the 21 reagent combinations are listed in Appendix A Table A2. In Table A1, the reagent combinations shaded gray are those with outliers in the calculated results.

As can be seen from Appendix A Table A1 and Table A2, for 70# matrix asphalt and SBS-modified asphalt, the number of outliers in the calculated results of the 21 reagent combinations varies from 0 to 4, with the lowest number (0) found for the reagent combination WFSD—that is, there are no outliers. Therefore, through the leap degree test, the best reagent combination scheme is finally determined to be WFSD.

4.3. Analysis of Results

In this study, the optimal reagent combination scheme was determined through the following three steps:

• (1)Based on three basic principles, eight chemical reagents with known surface energy parameters were preliminarily selected: W, F, E, G, S, D, B, and N. Then, 126 reagent combinations were formed randomly.
• (2)After the contact angles were measured using the Wilhelmy plate method and the asphalt surface energy parameters were determined using TLS, the rationality of the calculated results for all reagent combinations was analyzed, and those that were obviously unreasonable were eliminated. Finally, 21 reagent combinations that met the requirements were identified.
• (3)Through the leap degree test, the reagent combination scheme with the fewest outliers was selected from the 21 reagent combinations so as to achieve the research purpose of this study.

Through a series of laboratory tests, theoretical calculations, and data analyses combined, the combination scheme WFSD was selected, as there are no outliers in the calculated results of asphalt surface energy parameters with this combination, and the calculated results are closer to the real values. The calculated results of asphalt surface energy parameters with WFEG, the most commonly used reagent combination, are summarized in Table 5. According to the comparative analysis, for 70# matrix asphalt and SBS-modified asphalt, the error rate with the optimal reagent combination scheme, WFSD, in calculating the total surface energy of asphalt, $[eqn]$ , is reduced by 17.71% and 64.80%, respectively. Stability reflects the degree of dispersion, and it is necessary to use the homogeneity-of-variance test (F-test) to test the stability of the results of the two experimental methods. After calculation, the variance in the TLS group was significantly lower than that in the LS group (F(4,4) = 6.39, p = 0.006 < 0.05), indicating that TLS produces more stable results in calculating asphalt surface energy.

5. Conclusions

The main conclusions of this paper are as follows:

• (1)Compared to the LS method, the TLS method is more suitable for solving the linear equations established using the Young–Dupre equation. The TLS approach reduces the fitting error between the calculated and real values of asphalt surface energy parameters and improves the accuracy and stability of the calculated results. The coefficients of variation of the total surface energy of the 70# matrix asphalt and SBS-modified asphalt calculated using TLS are 15.92% and 24.65%, respectively, which are 21.65% and 6.45% reductions compared to those obtained with the LS method. The LS method only minimizes the single-axis (Z-axis) error, while the TLS method minimizes the full-dimensional vertical error, especially in ill-conditioned equations (for example, the change rate of the calculated values for the FSB combination is reduced from 23,060% to an acceptable range).
• (2)In this study, the number of outliers in the calculated results was determined by analyzing the rationality of the calculated results and using the leap degree test method. Finally, among 126 reagent combinations, the optimal reagent combination scheme was determined to be WFSD; there were no outliers in the calculated results of asphalt surface energy parameters, which were closer to the real values. For the 70# matrix asphalt and SBS-modified asphalt, compared to the most commonly used reagent combination, WFEG, the error rate of the reagent combination WFSD in calculating the total surface energy of asphalt was reduced by 17.71% and 64.80%, respectively.
• (3)The leap degree test method can quickly and accurately determine the reagent combination with the fewest outliers, and the selected resolution is more reasonable and reliable. Therefore, it can be used to select the optimal reagent combination. These research results can provide an experimental basis for the accurate calculation of asphalt surface energy parameters and assist in accurately analyzing the performance of adhesion between asphalt and aggregate. Thus, surface energy tests for more kinds of asphalt, asphalt binders, and aggregates will be conducted in the future, and the optimal reagent combination scheme suitable for surface energy tests of different road materials will be explored based on classification.

Through the optimization of calculation methods and reagent combinations, a more accurate framework for testing the surface energy of asphalt was established, which provides an experimental basis for the quantitative evaluation of the adhesion between asphalt and aggregate. In the future, this approach can be further extended to analyze more asphalt types (such as mortar and modified asphalt) and aggregate surface energy tests, explore the optimal reagent combination of different road materials based on classification, and improve the application of surface energy theory in road engineering.