Design of Experiments for Evaluating the Relevance of Change in Test Method for Kinematic Viscosity of Opaque Oils

: Viscosity is a physicochemical property that evaluates the resistance that fuel offers to flow, influencing the engine’s operation and combustion process. Its control is aimed at good fuel atomization and the preservation of lubricating characteristics. Changes in viscosity can lead to wear on various parts of the engine. Viscometers typically measure the viscosity of fuels in the oil and gas industry. These instruments can measure the time it takes for a fluid to move a given distance through a pipe or the time it takes for an object of a given size and density to pass through the liquid. The traditional test method, ASTM D445, differentiates the procedure for opaque liquids from transparent ones; that is, it requires a warm-up of the sample between 60 ◦ C and 65 ◦ C for 1 h. This additional step can overload laboratory routines, although it is not guaranteed to have a metrologically significant effect on the final result. Thus, this study evaluated the relevance of complying with this step in the test method for the kinematic viscosity of opaque liquids using a 3 2 factorial experimental design. Based on the F test, p -value, confidence intervals, and percentage contribution of the sum of squares approaches concerning the regression analysis, one concluded that the warm-up time was not a relevant factor in the kinematic viscosity, specifically of very low sulphur fuel oil, Brazilian fuel oil, and atmospheric residue diluted with diesel oil, which are fluids at room temperature.


Introduction
Many petroleum derivatives are used as lubricants, essential for the correct functioning of some equipment that depends on the appropriate viscosity of the fluid used.Therefore, the viscosity of these petroleum fuels is important for estimating ideal storage, handling, and operating conditions.Therefore, reliable viscosity determination is essential for many product specifications [1].
Kinematic viscosity is the relationship between momentum transport and momentum storage.Such relationships are called diffusivities, with dimensions of length squared divided by time, and the SI unit is a square metre divided per second (m 2 /s).However, in the oil and gas industry, it is normally converted to cSt (1 mm 2 /s = 10 −6 m 2 /s = 1 cSt).Among the transport properties for heat, mass, and momentum transfer, kinematic viscosity is momentum diffusivity.
The literature presents some studies concerning the kinematic viscosity determined with ASTM D445 [1], such as the investigation of the difference in vehicle engine performance due to degradation of engine oil properties [2]; the kinematic viscosity of engine oils is reduced by the influence of the accumulation of gasoline-bioethanol mixtures [3]; the characterisation of used oil distillate [4]; and the investigation of kinematic viscosity and other physicochemical properties arising from storage in biodiesel mixtures derived from used cooking oils [5].Few studies shed light on the kinematic viscosity of opaque liquids determined using ASTM D445 [6][7][8].Recently and in a still timid manner, some studies using the design of experiments have been applied to fuel and biofuels, such as the prediction of the physical properties of biofuels [9], the evaluation of the correlation between the ratio of light biofuel and refined palm oil blends on their density and kinematic viscosity properties [10], and the minimisation of the kinematic viscosity and maximisation of the biodiesel yield during a transesterification reaction [11].Furthermore, kinematic viscosity has been used in several studies regarding the quality of biodiesel [12][13][14][15].
The analytical procedure for testing residual fuel oils and opaque liquids recommends heating the sample in the original container at 60 • C and 65 • C for 1 h, unlike the procedure for transparent liquids [1].
Uniform heating time for opaque oils may be important; however, from a physical point of view, if there are no chemical transformations in the liquid or no dissolution process, the heating rate does not affect the physicochemical properties.On the other hand, the purpose of warming up the product before testing is to facilitate its pouring into the viscometer tube and the dissolution of paraffins, if present.
Therefore, this study aimed to evaluate metrologically, through the design of experiments, whether there are significant differences in the kinematic viscosity results, specifically of very low sulphur fuel oil (VLSFO), Brazilian fuel oil, and atmospheric residue diluted with diesel oil (RAT), which are fluids at room temperature, when they receive the same treatment as transparent liquids; that is, they are not heated.

Design of Experiments
Many experiments aim to evaluate the effects and possible correlations of two or more factors.Whenever it is necessary to reduce the number of experiments and increase the maximum amount of information about the system, factorial designs are recommended, as they are more efficient and effective.When using a factorial design, all possible combinations of factor levels are examined in each complete trial or replication of the experiment.For instance, for levels a of factor A and b of factor B, in addition to information regarding their own levels, it is possible to have information about combinations ab, which are normally called interactions.
Variations in the response (output magnitude) can be observed when each factor (input quantity) passes from one level to another, which can be defined as the effect of the factor.When this effect refers to the main factors of interest in the experiment, it is called the main effect.The simplest case is a two-factor factorial experiment with both design factors at two levels.The low level is also called "−1", while the high level is named "+1".The main effect of factor A in a two-level design can be understood as the difference perceived in the average output quantity when going from the low level to the high level of A [16].
The simplest case is a factorial design with two factors.To move on to the general case, let yijk be the response observed when factor A is at the i-th level (i = 1, 2, ..., a) and factor B is at the j-th level (j = 1, 2, ..., b) for the k-th replica (k = 1, 2, ..., n).In general, a two-factorial design will appear as in Table 1.To guarantee randomness to the process, abn observations are performed in a completely random sequence., and  .define the total of all observations related to the ith level of factor A, the jth level of the factor B, and its grand total, respectively.And their respective averages can be calculated as follows: Finally, the corrected total sum of squares is written as The total sum of squares has been spread over a sum of squares due to the "rows", or factor A (SSA); a sum of squares due to "columns", or factor B (SSB); a sum of squares related to the interaction between A and B (SSAB); and a sum of squares concerning the residual or error (SSE).This information is available in algorithms for two-way ANOVA with interaction.
Presuming that the data are normally homoscedastic and independently distributed, the relationship between variances, mean squares, and their respective degrees of freedom can be used to evaluate the relevance of each parameter (Table 2).

Material and Methods
In this section, the products and test method used are described.
A y a11 , y a12 , . . ., y a1n y a21 , y a22 , . . ., y a2n y ab1 , y ab2 , . . ., y abn y i.. , y .j. , and y ij.define the total of all observations related to the ith level of factor A, the jth level of the factor B, and its grand total, respectively.And their respective averages can be calculated as follows:

(1)
The total sum of squares has been spread over a sum of squares due to the "rows", or factor A (SS A ); a sum of squares due to "columns", or factor B (SS B ); a sum of squares related to the interaction between A and B (SS AB ); and a sum of squares concerning the residual or error (SS E ).This information is available in algorithms for two-way ANOVA with interaction.
Presuming that the data are normally homoscedastic and independently distributed, the relationship between variances, mean squares, and their respective degrees of freedom can be used to evaluate the relevance of each parameter (Table 2).

Squares (MS) F
Sample (A)

Material and Methods
In this section, the products and test method used are described.

Experimental
The study, relating to samples using VLSFO (very low sulphur fuel oil), fuel oil (OCA1), and atmospheric residue diluted with diesel oil (RAT) (Table 3), was carried out in August 2022 in a Brazilian fuel storage terminal without heating, with heating between 60 • C and 65 • C for 30 min, and heating between 60 • C and 65 • C for 60 min.The viscometers used were the routine Cannon-Fenske reverse-flow type.The glass capillary viscometers, the temperature measuring, and the timing devices, Figure 1, were calibrated within their expiration dates.

Experimental
The study, relating to samples using VLSFO (very low sulphur fuel oil), fuel oil (OCA1), and atmospheric residue diluted with diesel oil (RAT) (Table 3), was carried out in August 2022 in a Brazilian fuel storage terminal without heating, with heating between 60 °C and 65 °C for 30 min, and heating between 60 °C and 65 °C for 60 min.The viscometers used were the routine Cannon-Fenske reverse-flow type.The glass capillary viscometers, the temperature measuring, and the timing devices, Figure 1, were calibrated within their expiration dates.
In a simplified way, the procedure used for determining the kinematic viscosity of these opaque oils was as follows: (i) we charged the reverse-flow-type viscometer with the opaque oils; (ii) we placed the viscometer in the viscometer bath at the desired test temperature; (iii) we let the charged viscometer remain in the bath for enough time to reach the test temperature; and (iv) after ensuring that the opaque oil flowed freely, we measured to the nearest 0.1 s the time required for the sample to flow between the calibrated marks on the viscosimetric tube.
At least two kinematic viscosity determinations are necessary to consider a valid result.The test results for the kinematic viscosity were reported in four significant figures.

Test Method
The test method details an experimental procedure for evaluating the kinematic viscosity of opaque and transparent liquid petroleum products.The analytical principle In a simplified way, the procedure used for determining the kinematic viscosity of these opaque oils was as follows: (i) we charged the reverse-flow-type viscometer with the opaque oils; (ii) we placed the viscometer in the viscometer bath at the desired test temperature; (iii) we let the charged viscometer remain in the bath for enough time to reach the test temperature; and (iv) after ensuring that the opaque oil flowed freely, we measured to the nearest 0.1 s the time required for the sample to flow between the calibrated marks on the viscosimetric tube.
At least two kinematic viscosity determinations are necessary to consider a valid result.The test results for the kinematic viscosity were reported in four significant figures.

Test Method
The test method details an experimental procedure for evaluating the kinematic viscosity of opaque and transparent liquid petroleum products.The analytical principle is based on measuring the time it takes a portion of the sample to flow under gravity using a calibrated glass capillary viscometer [1].
The time is measured for a constant volume of liquid to flow under gravity through the capillary of a calibrated viscometer under a reproducible drive head and at a known, tightly controlled temperature.The kinematic viscosity (determined value) is the product of the measured flow time and the viscometer calibration constant.Two such determinations are required to calculate a kinematic viscosity result, i.e., the average of two acceptable determined values.

Results and Discussion
The matrix, 3 2 , with the experiment and design-considering two factors and three levels-is detailed in Table 4.A is the product and B is the warm-up time, with two replicates of the output quantity, kinematic viscosity, R1, and R2.The graphical representation of the matrix 3 2 is shown in Figure 2.
Metrology 2024, 4, FOR PEER REVIEW 5 based on measuring the time it takes a portion of the sample to flow under gravity using a calibrated glass capillary viscometer [1].The time is measured for a constant volume of liquid to flow under gravity through the capillary of a calibrated viscometer under a reproducible drive head and at a known, tightly controlled temperature.The kinematic viscosity (determined value) is the product of the measured flow time and the viscometer calibration constant.Two such determinations are required to calculate a kinematic viscosity result, i.e., the average of two acceptable determined values.

Results and Discussion
The matrix, 3 , with the experiment and design-considering two factors and three levels-is detailed in Table 4.A is the product and B is the warm-up time, with two replicates of the output quantity, kinematic viscosity, R1, and R2.The graphical representation of the matrix 3 is shown in Figure 2. To help interpret the results of this experiment, it is useful to plot the average responses for each treatment combination in Figure 3  To help interpret the results of this experiment, it is useful to plot the average responses for each treatment combination in Figure 3.The absence of significant interaction is indicated by the parallelism of the lines.Visually, the kinematic viscosity is not significantly affected by the warm-up time.
The contrasts-linear combination of parameters; effects, change suffered by the response variable when moving from the low level of the factor to the highest level; coefficients (and their standard errors); and the sum of squares (SS)-were calculated considering an  of 1418.615(Table 2).For factor A (product),  The absence of significant interaction is indicated by the parallelism of the lines.Visually, the kinematic viscosity is not significantly affected by the warm-up time.
The contrasts-linear combination of parameters; effects, change suffered by the response variable when moving from the low level of the factor to the highest level; coefficients (and their standard errors); and the sum of squares (SS)-were calculated considering an MS error of 1418.615(Table 2).
These data were subjected to a two-way ANOVA with interaction (Table 5).Based on the p-value, factor A, the product has significantly impacted the kinematic viscosity output quantity since the p-value was less than 5% [18,19].However, factor B, warm-up time, did not significantly impact the kinematic viscosity output quantity since the p-value is greater than 5%.A lack of interaction between the product and warm-up time variables can also be observed.
A third approach used the confidence intervals (CI) of the regression coefficients, which were calculated as follows: CI = coe f f icient ± t n−2 × standard error [20].
CI A = (361.9748, 408.0735);CI B = (−35.0535, 11.0452); and CI AB = (−40.8608, 15.5983).One only considers statistically significant effects whose estimates (obtained in the experiment) are greater in absolute value than the product of the standard error and critical values of the Student's t-distribution because, this way, the confidence interval does not include the zero value.Thus, only factor A proved relevant, confirming the two previous approaches.
Lastly, using the percentage contribution of the sum of squares approach [17], the relevance of each factor was evaluated as follows: SS Regression = 1729.201,781,928.9= 0.096%; and SS AB SS Regression = 1276.3881,781,928.9= 0.071%.The warm-up time and its interaction with the product did not have a significant impact on the kinematic viscosity output quantity since their percentages were less than 5%.

Table 1 .
Experimental matrix for a factorial design with two factors.

Table 2 .
Algorithms for two-way ANOVA with interaction.

Table 2 .
Algorithms for two-way ANOVA with interaction.