AHybrid Neutrosophic-Grey Analytic Hierarchy Process Method: Decision-Making Modelling in Uncertain Environments

+e analytic hierarchy process (AHP) is recognised as one of the most commonly applied methods in the multiple attribute decision-making (MADM) literature. In the AHP, encompassing uncertainty feature necessitates using suitable uncertainty theories, since dealing efficiently with uncertainty in subjective judgements is of great importance in real-world decision-making problems. +e neutrosophic set (NS) theory and grey systems are two reliable uncertainty theories which can bring considerable benefits to uncertain decision-making.+e aim of this study is to improve uncertain decision-making by incorporating advantages of the NS and grey systems theories with the AHP in investigating sustainability through agility readiness evaluation in large manufacturing plants. +is study pioneers a combined neutrosophic-grey AHP (NG-AHP) method for uncertain decisionmaking modelling. +e applicability of the hybrid NG-AHP method is shown in an illustrative real-case study for agility evaluations in the Iranian steel industry. +e computational results indicate the effectiveness of the proposed method in adequately capturing uncertainty in the subjective judgements of decision makers. In addition, the results verify the significance of the research in group decision-making under uncertainty. +e practical outcome reveals that, to become a more sustainable agile steel producer in the case country, they should first focus on the “organisation management agility” as the most significant criterion in the assessment followed by “manufacturing process agility,” “product design agility,” “integration of information system,” and “partnership formation capability,” respectively.


Introduction
In recent years, corporations have moved to the centre of focus in the sustainability debate. e reason for this is that they are considered to be responsible for enormous negative impacts on the environment and society [1]. Sustainability aims to produce a dynamic balance between the three sustainability dimensions: environmental, economic, and social over time [2]. In the contexts of sustainable operations and agility, agile manufacturing (AM) should first be defined. AM has the property of robustness, which means that AM systems must be able to tolerate changes and interruptions within the given demand requirements. erefore, AM operations can be seen as inherently sustainable. In other words, agility and sustainability are interconnected in this sense, and an agile system can have the potential capability to work as a sustainable system. is link has not been adequately researched in the literature. In this study, we have considered the general definition of sustainability taking into account all three pillars of sustainability in the steel industry.
In addition to the predicted increase in the global population up to 9.7 billion people by 2050 [3], the steel requirement per capita is also expected to increase by 2050 to 11.8 tons. Steel production is responsible for around 7% of global greenhouse gas (GHG) emissions. Bearing in mind the increasing production volume per capita and the necessity to decrease our global greenhouse gas emissions to tackle climate change, it is evident that the steel industry has to shift to more sustainable processes [4]. Steel manufacturing needs more attention from the sustainable development perspective particularly in developing countries. Iskanius et al. [5] investigated the leading forces and abilities for agility in Finnish steel product manufacturing and concluded that the need for agility is clearly recognised in the traditional steel industry and it has to be considered in the long-term strategy planning of steel manufacturing enterprises. us, the research question to be addressed is as follows: RQ: how can agility readiness impact sustainable engineering decisions under uncertain decision-making environment in steel manufacturing?
Decision support tools such as multiple attribute decision-making (MADM) can now be identified as invaluable business analytic methods for helping large organisations to move forward towards providing sustainable operations by developing agility in their manufacturing processes. In the MADM literature, there are methods that have simple implementation and flexibility such as analytic hierarchy process (AHP), best-worst method (BWM) [6], level-based weight assessment (LBWA) [7], and full consistency method (FUCOM) [8].
e AHP is one of the most commonly practised MADM methods [9] mainly because of its ease of application and its flexibility for integration with various methods. An abundance of studies in the literature is focused specifically on applications of AHP such as traffic accessibility [10], advertising media selection [11], selecting e-purse smart card technology [12], and topic popularity selection [13], to name a few. AHP has been utilised to assess complicated multiattribute alternatives by collecting opinions of a group of decision makers (DMs). e feature of inclusion of subjective factors has been considered as one of the AHP's advancements compared to other MADM methods [14]. Many studies have focused on the fuzzy set (FS)-based extension of AHP, namely, fuzzy AHP (F-AHP), so as to capture uncertainty [15][16][17][18]. However, few studies have considered the extension of AHP simultaneously with other uncertainty theories such as grey systems and neutrosophic set (NS) theories, which are able to enhance decision-making process under uncertain environment. ere are only a few recent developments and applications of AHP and NS theory in the literature [19][20][21][22][23][24][25]. Furthermore, only limited research, which is directly related to grey AHP method, has been carried out [26][27][28][29][30][31][32][33].
is gap has motivated the current research to develop AHP under hybrid grey and NS decisionmaking environments to deal with uncertainty embedded in human subjective judgements which can incorporate the advantage of both in one decision-making model. e NS theory is able to independently quantify the indeterminacy membership function values. Unlike the FS theory, the NS theory has the capability to express the information about rejection. ere are growing applications of the NS theory in the decision-making literature [6,24,[34][35][36][37][38][39][40]. Smarandache [41] introduced the NS theory and in [42] thoroughly elaborated on the distinctions between NS and intuitionistic fuzzy set (IFS) theories by providing explanatory examples. It has improved the IFS theory which was initially introduced by Atanassov [43] as an extension of Zadeh's FS theory [44]. Besides the NS theory, Pythagorean fuzzy set (PFS) theory was introduced by Yager [45] and has been a recent extension of IFS theory which is drawing the attention of researchers in the realm of decision-making under uncertainty [46][47][48]. D-numbers are also introduced to deal with uncertainty in decision-making [49,50]. In addition to NS theory, grey system theory compared to many mainstream uncertainty theories, such as FS theory, has appreciable features, particularly when it is necessary to deal with uncertain data, and lack of information such as (1) generating satisfactory results utilising a relatively small data volume; (2) producing robust results regarding the noise, and lack of modelling information; and (3) yielding fairly flexible, nonparametric assumptions, and a general way to integrate fuzziness into a problem [16]. Smarandache [51] discussed the NS theory and grey systems all together. In most grey AHP studies, the utilised whitenisation functions cause information loss by converting grey information into crisp values. Moreover, calculating the consistency ratio (CR) to check the consistency among evaluations of DMs in pairwise comparison matrices is another cause of concern. Additionally, in several studies, the integration method is utilised as a combination of grey relational analysis (GRA), and AHP or grey incidence analysis (GIA) and AHP [52][53][54][55][56]. e GRA and GIA are characterised under the grey system theory concept as two distinct MADM methods. e main common feature of these studies is that AHP is applied for calculating criteria weights, and then either GRA or GIA is used for evaluation of alternatives. is category of studies should not be mixed up with the grey-based AHP method, because GRA/GIA-AHP methods do not apply AHP in combination with grey systems theory.
In this paper, we take advantage of operational research (OR) tools from the realm of MADM to evaluate agility readiness of Iranian steel manufacturing corporations with the aim of developing sustainable operations. Two methods, namely, G-AHP (i.e., grey AHP) and N-AHP (i.e., neutrosophic AHP), are combined, and the application of the model is demonstrated as a hybrid method NG-AHP (i.e., neutrosophic-grey AHP) in an agility evaluation case in the Iranian steel industry. It is believed that based on the provided method, the uncertainty of DMs can be best handled via hybrid neutrosophic-grey uncertainty theories. e N-AHP approach is an integration of the NS theory with the AHP, in which the single-valued trapezoidal neutrosophic numbers (SVTNNs) are utilised in the AHP calculations [19]. e proposed G-AHP method has furthered the existing grey AHP methods in two major ways. Firstly, it preserves the grey characteristics of grey numbers during calculation steps by reducing information loss, specifically by omitting the need for whitenisation function deployment. Secondly, it ensures assessment consistency in pairwise comparisons by introducing two importance rating scales and constructing the pairwise comparisons based on the suggested procedure. It also obtains the aggregated opinions of DMs efficiently, while also handling the inherent ambiguity in the subjective judgements of DMs through preserving grey values.
To the best of our knowledge, OR tools such as MADM methods have not been applied extensively for sustainable development and there is a gap in the literature in this matter [57]. Moreover, only a few studies have explored sustainability through the agility perspective in manufacturing settings. In other words, achieving sustainable operations in engineering through AM or sustainable agility has not been a well-researched topic in the literature. It is explained that agility and sustainability are closely connected, meaning that studying agility in manufacturing enterprises can thus lead to better understanding of sustainable development. To bridge this gap, we are exploring agility readiness in the steel manufacturing business by applying a novel combined MADM method. e current research features the following three specific contributions: (1) Investigating sustainable engineering by agility readiness evaluation in an Iranian steel manufacturing setting.
(2) Extending group AHP to a grey environment (G-AHP) while preserving the grey characteristics of the judgements with fully consistent evaluations in the pairwise comparison matrices. e first characteristic of the proposed G-AHP is that no whitenisation function is needed unlike most other grey AHP methods which would utilise whitenisation functions to get crisp values. e second feature is that no CR calculation is required in pairwise comparison matrices due to the way they are established, which also helps save time and cost. To the best of our knowledge, no grey AHP method in the literature has these two traits simultaneously and with a straightforward procedure.
(3) Integrating N-AHP [19] with G-AHP (i.e., NG-AHP) in a real-world agility evaluation case in the Iranian steel industry to illustrate its capability and versatility in one hybrid methodological application. e application of the hybrid methodology (i.e., NG-AHP) was shown to reveal the benefits of both methods in one single framework, while also emphasising the synergistic effects of the two in one single framework and also overcoming their drawbacks at the same time.
In Figure 1, a generic hierarchical structure is shown in which the applied levels of proposed methods are presented. e calculation steps of each method are shown in Figure 2. e proposed hybrid NG-AHP is comprised of two separate methods including integration of the N-AHP (Section 4), and G-AHP (Section 5). e criteria weights are calculated by the N-AHP, and the importance weights of alternatives are obtained by G-AHP. Ultimately, the weights are integrated, and alternatives are ranked to calculate the total weights of alternatives in the final decision matrix (Section 6). In Section 7, findings are discussed, and the paper is concluded in Section 8.

Sustainability and Agility
Agility is characterised as the ability to react to and handle unpredictable changes and encompasses cost reduction, quality improvement, delivery, and service improvement. Agility lies in the domain of AM which is the ability to meet volatile business requirements with adaptability and has been developed in response to lean manufacturing (LM) systems [58,59]. Leanness aims at maximising profit through cost reduction, while agility tries to maximise profit by providing precisely what a customer needs [60]. Agility is also considered as the interface between the company and the market [61].
Sustainability generally concentrates on protecting natural resources against exploitation via productivity and competitiveness by manufacturing and service organisations. However, the concept of sustainability includes two key aspects other than the environmental aspect, which are economic and social [62,63]. ereby, the three dimensions of sustainability (i.e., environmental, economic, and social) have to be considered and treated equally. Gunasekaran and Spalanzani [62] investigated sustainable business development (SBD) in manufacturing and services, which has been regarded as a critical issue due to many causes such as climate change and natural disasters. Sustainability efforts can be included in all stages of a supply chain from product design and manufacturing to the product end-of-life stage such as remanufacturing [64]. Rostamzadeh et al. [65] investigated sustainability issues in the supply chain risk management domain by applying an integrated fuzzy MCDM based on TOPSIS and criteria importance through intercriteria correlation (CRITIC). Ivory and Brooks [66] offered a conceptual framework illuminating the strategic agility metacapabilities (resource fluidity, collective commitment, and strategic sensitivity) and related practices/ processes that firms use to effectively deal with corporate sustainability with a paradoxical lens. ere would be an intuitive possible connection between agility and sustainability because more efficient and improved quality production by being quick and flexible in agile manufacturing potentially would lead to less production waste and carbon emissions and ultimately to more sustainable production. Carvalho et al. [60] recognised the trade-offs between lean, agile, resilient, and green (LARG) management systems as a probable pathway towards a more sustainable system. It is also indicated that agility and sustainability are regarded as performance measures for contemporary enterprises. In the current manufacturing scenario, agility needs to be matched with sustainability [67]. Pham and omas [68] suggested that for firms to be competitive, they should achieve an effective level of leanness, agility, and sustainability that associates with change and uncertainty in an operational system and the individual business environment. Flumerfelt et al. [59] investigated theories and practices of agile and lean manufacturing systems to gain an understanding of whether these employ Mathematical Problems in Engineering 3 sustainability or not. ey recognised AM operations are sustainable because sustainability means ability to endure, and AM systems must be robust which means they are capable to endure alterations under various demand circumstances.

Preliminaries
3.1. Neutrosophic Set eory. Some basic definitions of NS theory are provided in this section [69].
Definition 1 (see [41]). Let U be a finite set of objects and let x signify a generic element in U. e NS A in U is characterised by a truth-membership function T A (x), an indeterminacy-membership function I A (x), and a falsity- , and F A (x) are the elements of ]0 − , 1 + [. It can be shown as Definition 2 (see [70]). Let U be a finite set of elements, and let x signify a generic element in U. A single-valued neutrosophic set (SVNS) A in U is defined by a truth-membership function T A (x), an indeterminacy-membership function I A (x), and a falsity-membership function F A (x). T A (x), I A (x), and F A (x) are the elements of [0, 1]. It can be shown as Definition 3 (see [71]) and w a , u a , y a ∈ [0, 1] is a particular single-valued neutrosophic number (SVNN) whose T a (x), I a (x), and F a (x) are presented as the following equations, respectively: Definition 4 (see [72]). Given a � < (a 1 , b 1 , c 1 , d 1 ); w a , u a , y a > , and then equations (6) and (7) are true: When (8) and (9) are true: Definition 5 (see [72]). Given a � < (a, b, c, d); w a , u a , y a > and a, b, c, d > 0. en, the score function of a can be calculated in accordance with the following equation: Definition 6 (see [72]). In order to compare two SVTNNs ; w a j , u a j , y a j > (j � 1, 2, . . . , n) be a set of SVTNNs, then a trapezoidal neutrosophic weighted arithmetic averaging (TNWAA) operator is computed on the basis of [72] TNWAA a 1 , a 2 , . . . , a n � n j�1 p j a j where p j is the weight of a j (j � 1, 2, . . . , n) while p j > 0, and n j�1 p j � 1.

Subtraction, Division, and Inverse of SVTNNs.
e subtraction and division of simplified SVNNs (or singlevalued neutrosophic values) and SVNSs are introduced by Smarandache [73] and Ye [74], respectively. Rani and Garg [75] also studied subtraction and division operations on interval neutrosophic sets. In this section, subtraction, Mathematical Problems in Engineering division, and inverse of SVTNNs in general nonsimplified form are defined.

Subtraction of SVTNNs
Note: for a negative value, replace it with zero. For a value of over one, replace it with one.
Proof. Let us consider equation (13) where By adding neutrosophically, b to the sides of equation (13)-(16) results, en, and

then the division of the two SVTNNs is shown in
Note: for a negative value, replace it with zero. For a value of over one, replace it with one.
Proof. Let us consider equation (18) where By multiplying neutrosophically, b to the sides of equation (18)-(21) is obtained: en, and and w a , u a , y a , ∈ [0, 1], then the inverse of a is represented in Note: for a negative value, replace it with zero. For a value of over one, replace it with one.

Proof.
Let us consider equation (23), where a � < (a 1 , b 1 , c 1 , d 1 ); w a , u a , y a > , and a − 1 � < (x, r, z, s); w c , u c , y c > , en, based on the division rule of two SVTNNs referring to equation (17), the proof is provided. □ 3.3. Grey System eory. In this section, some basic definitions of grey systems theory are provided [69].
Definition 8. A grey number ⊗X is defined as an interval with known upper, and lower bounds which are shown by X and X, respectively, but there is no known distribution information for X [76,77], as represented in Definition 9. Given ⊗X 1 � [X 1 , X 1 ] and ⊗X 2 � [X 2 , X 2 ] are two grey numbers, then the basic operations of grey numbers can be defined as follows [78,79]: Definition 10. e length of a grey number ⊗X is defined as Definition 11. Comparison of grey numbers [80].
Given ⊗X 1 � [X 1 , X 1 ] and ⊗X 2 � [X 2 , X 2 ] are two grey numbers, the possibility degree of ⊗X 1 ≤ ⊗ X 2 can be defined as follows: where L * � L( ⊗ X 1 ) + L( ⊗ X 2 ). ere are four possible cases on the real number axis to determine the relationship between ⊗X 1 and ⊗X 2 : Definition 12 (see [81]). Whitenised (whitened or crisp value) of a grey number is a deterministic number with its value between the upper and lower bounds of a grey number ⊗X. e whitenised value x (λ) can be defined as equation (31) in which λ is whitening coefficient, and λ ∈ [0, 1]: For λ � 0.5, equation (32) will be resulted: Definition 13 (see [81,82]). Given ⊗X 1 � [X 1 , X 1 ] and ⊗X 2 � [X 2 , X 2 ] are two grey numbers, then the distance between ⊗X 1 and ⊗X 2 can be calculated as signed difference between their centres as shown in Definition 14 (see [79]). Given ⊗X � [X, X] is a grey number, and >0; then, equation (34) is resulted:

The N-AHP Method
e N-AHP method follows the steps below as introduced in [19].
Step 1 (hierarchical structure): it is an essential step to establish a hierarchy, representing the goal, criteria, and alternatives because it makes the problem more comprehensible.
Step 2 (pairwise comparison matrix): the DMs evaluate elements (i.e., alternatives or criteria), using the Saaty rating scale Table 1. In the experts' judgements questionnaire, DMs choose a linguistic phrase representing the importance degree of each element in comparison to others. Given C 1 , C 2 , . . . , C n signify the elements, and a ijk shows a quantified evaluation on a pair of elements, C i and C j by k th DM (k � 1, 2, . . . , p). is leads to a pairwise comparison matrix as represented in [84,85] Step 3 (calculating CR): referring to Saaty's suggestion [86], a consistency test has to be conducted to differentiate the consistent comparisons from the inconsistent comparisons. See equation (36) and Table 2. If the value CR ≥ 0.1, then the DMs have to do a revision in their evaluations [88]: Step 4 (replacing the linguistic information with the SVTNNs): the elements of the pairwise comparison matrices are replaced with the corresponding SVTNNs using the scale shown in Table 3 (see Section 3.2.3 for calculating inverse of an SVTNN).
Step 5 (aggregating the opinions of DMs in SVTNNs): to aggregate the opinions of DMs, the TNWAA operator is used, as shown in equation (11).
Step 6 (neutrosophic synthetic values): the neutrosophic synthetic value of each element (S i ) is computed based on where n is the number of elements and η ij is the (i, j) th element of the aggregated pairwise comparison matrix.
Step 7 (determining the final importance weights): this is calculated based on equation (38), and the final importance weights are shown by W i which are in SVTNNs. In order to compare weights, equation (10) is used:

The G-AHP Method
e proposed G-AHP is inspired by the fuzzy Delphi method in [84,85]. e main characteristics of the proposed G-AHP method compared to other similar grey AHP methods in the literature are as follows: (1) no whitenisation function is used; all the calculations from the beginning to the end are in grey numbers, and in accordance with basic grey operations rules (Section 3.3).
is preserves the grey characteristics of the values and judgements and helps reach a more valid outcome; (2) no consistency calculation is needed; the pairwise comparison matrices are constructed in a way that CR values of any pairwise comparison matrix are zero, and evaluations are fully consistent; (3) two judgement scales are introduced. Here, the 5-point judgement or importance scale (Table 4) is utilised by DMs to show the significance of each element individually. e 9-point relative importance scale (Table 5) is constructed for obtaining an importance comparison of each element compared to other elements in pairwise comparison matrices.
It is assumed that there are p DMs, and let ρ � (ρ 1 , ρ 2 , . . . , ρ p ) T be the importance weight vector of the DMs where p k�1 ρ k � 1, ρ k ≥ 0, k � 1, . . . , p. It is also given that the decision-making model includes two finite sets of alternatives and criteria which are shown by � x 1 , x 2 , . . . , x n and C � c 1 , c 2 , . . . , c m , respectively. e steps of the grey weights' calculation applied in the proposed G-AHP method are represented as follows: (i) Step 1 (constructing the hierarchical structure): at this initial step, the hierarchical structure of the decision-making problem including goal, criteria, alternatives, or subalternatives will be constructed. (ii) Step 2 (asking for experts' opinions): the DMs are asked to evaluate elements (criteria or alternatives) on the basis of their significance. e DMs determine the relative importance of each alternative x i over x j or each criterion c i over c j by using the importance scale (Table 4).  Strong plus importance 7 Very strong importance 8 Very very strong importance 9 Extreme importance  〈 (9,9,9,9); 1, 0, 0〉 9 8 Mathematical Problems in Engineering scale ( Table 4). As shown in equation (39), in the case of comparing criteria, n should be replaced with m. e a ijk represents the relative significance of element i over element j from the viewpoint of the k th DM: (iv) Step 4 (weighted pairwise comparison matrices): ρ k is the importance weight of the k th DM which belongs to the interval [0, 1], and the greater the weight value, the more significant the DM's opinion is. According to each DM's importance weight ρ k and elements of matrices of equation (39), the β ijk values can be calculated based on (v) In the case of equal importance weights of DMs, there is no need to calculate equations (40) and (41), and simply equation (39) can be used. (vi) Step 5 (grey number calculation): in order to calculate grey numbers ⊗a ij , all evaluations are taken into account, considering the importance weight of each DM as equations (42)-(44), where c ij ≥ α ij : (vii) According to the aforementioned explanations, the weighted grey pairwise comparison matrix for alternatives is defined in equations (45) and (46). In the case of criteria, n should be replaced with m in equations (45) and (46), where μ � min(ρ k ), ∀k � 1, . . . , p, and π � max(ρ k ), ∀ k � 1, . . . , p: (vii) Step 6 (grey weight calculation): grey weight of each alternative (i.e., ⊗ W i ) can be calculated using equations (47) and (48). For criteria, n should be replaced with m in the following equations: Fairly good (FG) 6 Good (G)

The Case Application
e proposed NG-AHP method was applied to agility evaluations in the Iranian steel industry. Five agility evaluation criteria were observed as evaluation criteria and were applied to four steel companies [89]. e criteria were organisation management agility (C1), product design agility (C2), manufacturing process agility (C3), partnership formation capability (C4), and integration of information system (C5). Based on the data collected from the chosen experts, the aim was to identify the most relevant agility criteria. Subsequently, the steel enterprises were ranked according to the agility readiness criteria. Four steel manufacturing companies were investigated in the present research, and their names were anonymized as SC1, SC2, SC3, and SC4.
Here, the expert selection process and their importance weight assignment task were carried out based on the experts' knowledge and expertise in the related steel industry. Six DMs who were steel industry experts and were available to provide insights on agility readiness criteria evaluation, as well as being independent from the four steel companies, were selected. Brief profiles of the experts are represented in Table 6. e importance weight of each DM is provided as ρ � (0.15, 0.30, 0.10, 0.25, 0.15, 0.05) T regarding their knowledge and experience. e experts were initially contacted to participate in the study by completing two types of questionnaires for N-AHP and G-AHP, based on the scales provided in Tables 3 and 4. e acquired data are presented in Tables 7 and 8. e hierarchical structure of this problem is depicted in Figure 3. e proposed N-AHP was applied so as to obtain weights for five criteria. ese weights were used later in the G-AHP method to acquire the final ranking of steel companies. In Table 7, the initial pairwise comparison matrices based on the opinions of six DMs using the NS rating scale (Table 3) are shown (A 1 , . . . , A 6 ). e calculated CRs for each pairwise comparison matrix were 2.23%, 7.66%, 2.36%, 3.99%, 6.56%, and 7.57%, respectively; they were all below 10% indicating cardinal output-based consistency. e aggregation neutrosophic matrix was calculated based on TNWAA operator, and then by applying equations (37) and (38), final weights were estimated (Table 9). rough G-AHP, opinions of DMs were obtained for the evaluation of each steel company (SC1, SC2, SC3, and SC4) against criteria based on the scale provided in Table 4. e numerical values in Table 4 then were substituted for linguistic phrases (Table 8).
Here, only the weight computations of four steel companies based on C1 (organisation management agility) are presented to show how the G-AHP method works. e resulted weights then make up the first column of the final decision matrix as shown in Table 10. e pairwise comparison matrices of four steel companies based on C1 (organisation management agility) according to opinions of six DMs are denoted as A 1 , A 2 , A 3 , A 4 , A 5 , and A 6 as presented in Table 11. All the CRs for comparative matrices will be equal to zero due to the applied method of acquiring opinions of DMs. e interpretation of the values in linguistic terms can be figured out based on the scale represented in Table 5. ese values range from 0.33 with the corresponding verbal term extremely less important to 3 with the corresponding verbal term extremely more important.
e ⊗Z i values were calculated according to equation (47) and grey weights of each steel company (i.e., ⊗ W i ) were calculated using equation (48) as shown in Table 13. e same steps as C1 for the other four criteria (C2, C3, C4, and C5) were carried out to obtain weights of the four alternatives, i.e., steel companies (SC1, SC2, SC3, and SC4).

Results and Discussion
e final grey decision matrix is shown in Table 10. Taking into account weights of each criterion obtained from N-AHP (Table 9) as shown in bold in Table 10, the final grey importance weights of each steel company can be calculated. e final grey importance weights of each steel company can be achieved by multiplying the weights and adding them up (Table 14). e four steel companies were ranked based on the obtained importance weights which are shown in Table 14.
Comparisons were carried out with reference to Definition 11. e final ranking of steel companies is obtained as In respect of the crisp weights obtained from N-AHP (Table 9), it was revealed that C1 (organisation management agility) with the weight of 0.3262 is the most significant criterion in the assessment followed by C3 (manufacturing process agility), C2 (product design agility), C5 (integration of information system), and C4 (partnership formation capability) with weights of 0.2643, 0.2123, 0.1337, and 0.0635, respectively. e obtained weights of five criteria were utilised in the G-AHP method to reach the final ranking of steel companies. After applying G-AHP, regarding the final obtained weights in grey numbers (Table 14), it was concluded that steel company four (SC4) has the highest competence for agile strategies based on the five evaluation criteria. e SC3, SC2, and SC1 lay second, third, and fourth in the final prioritisation order, respectively.
It should be noted that either G-AHP or N-AHP has the capability to be used for the calculation of all the AHP steps in any similar decision-making problem separately. Noticing the difference that the final weights in the N-AHP will be in crisp values, the G-AHP will obtain grey values of weights. It has been demonstrated that G-AHP and N-AHP can operate together in one decision-making framework, namely, NG-AHP. Apart from separate merits of each of the two methods, this integration can provide a synergic effect and can be more beneficial by incorporating advantages of the NS and grey systems theories, simultaneously.
Organisation management agility (C1) is recognised as the most significant criterion leading a manufacturing firm towards agility and subsequent potential sustainability. Sharifi and Zhang [90] indicated that uncertainty in the business environment has been considered as the root of most failures in manufacturing industry and agility has two main factors for coping with unexpected, uncertain changes: (1) responding to change (anticipated or unexpected) in suitable ways and in due time and (2) taking advantage of changes as opportunities. Organisation management agility refers to the capability that enables a company to rapidly adapt in response to market changes by improving the organisational management procedures. Lozano et al. [91] mentioned that corporate sustainability (CS) has been recognised as one way of incorporating a sustainability agenda in the activities of an organisation in order to address the negative impacts of its operations on the environment and society. It is indicated that CS changes need to be integrated in soft organisational issues including values, visions, policies, and change management practices which are related to organisational systems of a company [91].
Manufacturing process agility (C3) is the second highly important factor. In prior research, some scholars suggested that internal business processes could be significant factors linking information technology (IT) capability and organisational performance [92], and a notable aspect of internal business processes is business process agility [93,94]. Jayal et al. [95] indicated that achieving sustainability in manufacturing needs a holistic perspective spreading over both internal and external parts of an organisation. us, manufacturing process agility as an internal part can play a central role in attaining sustainability particularly in economic and environmental dimensions [96].
Findings revealed that product design agility (C2) is the third element that should be taken into consideration by practitioners to develop sustainable operations. It is in close connection with information systems (ISs) and computer technologies as they can play a considerable role in facilitating time reduction in product design and development [94,97]. Vinodh [67] indicated that integration of agility and sustainability results in various advantages and that one of these is product design and development. Economic and environmental sustainability are significantly influenced by product design. Concentrating simultaneously on the green design of products during the premanufacturing stage and also on design agility in order to respond to customer needs can have positive impacts on environmental sustainability [98].
Integration of the information system (C5) which is connected to IT is regarded as the fourth agility criterion leading to sustainability. IT has been widely recognised as being a crucial factor for the survival and growth of an organisation [93,99]. Aggoune et al. [100] in the evaluation of the relationship between the agility and sustainability of the information system indicated that there is no sustainability without agility. Frayret et al. [97] highlighted the significance of speed factor in agile manufacturing and declared that computer technologies are platforms for agility [101]. Information disclosure as one of the social sustainability indicators [52,102] can be more applicable through integrated information systems. e reason is that more related information regarding materials being used during the production process and information about carbon emissions would become available. Additionally, application of blockchain technology by providing decentralised and immutable storage of verified transaction data can be significant in this matter [19].
Finally, partnership formation capability (C4) is another component building sustainable processes with notable importance in agile manufacturing. Yusuf et al. [103] named partnerships as one of the attributes agile organisations have, including trust-based relationship with customers/suppliers, close relation with suppliers, strategic link with customers, and rapid partnership formation. is criterion can relate to social sustainability for instance, by extending partnership with external stakeholders through green outsourcing which would lead to higher satisfaction of the community that has an interest in the outcomes from the actions of an organisation [52]. It can also result in economic sustainability as an outcome of joint ventures or green procurement contracts [104,105].

Conclusions
In this article, we assess agility in the steel manufacturing industry with criteria which can help develop sustainable engineering operations in organisations. e trade-offs between lean, agile, resilient, and green (LARG) management systems as a probable pathway towards a more sustainable system is studied in the literature. It is also indicated that agility and sustainability are regarded as performance measures for contemporary enterprises. In this context, agility and sustainability are considered to be interconnected even though their link has not been adequately researched. In this study, a decision-modelling approach was introduced to show the application of the proposed hybrid method (i.e., NG-AHP) in an agility evaluation case in the Iranian steel industry. In the MADM field, a variety of uncertainty theories such as FS, IFS, PFS, and grey systems theories are applied to deal with the unsuitability of crisp values for the efficient modelling of reallife problems. Achieving sustainable operations via agile manufacturing is regarded as a vital management paradigm for making the production system more efficient and streamlined. Our proposed method is comprised of two MADM methods, namely, N-AHP and the G-AHP under uncertain decision environments. e importance weights of agility evaluation criteria in the Iranian steel industry were determined by applying the N-AHP, and then, the G-AHP was utilised to rank steel companies.
is study contributes to the literature in three main ways to answer the research question of how the agility readiness impacts sustainable engineering decisions under uncertain decision-making environment in steel manufacturing. First, it explored sustainable engineering by agility readiness evaluation in an Iranian steel manufacturing setting. Second, a new grey-based AHP method, namely, G-AHP, inspired by the fuzzy Delphi method was proposed. is provided a distinctive approach for the method compared to similar grey AHP methods in the literature in two main ways by introducing two judgement scales (Tables 4 and 5). e proposed G-AHP preserves the grey characteristics of the judgements in the AHP computing steps while preserving fully consistent evaluations in the pairwise comparison matrices. ird, a real-case example of agility evaluation in the steel industry was provided to demonstrate the joint applicability of the two methods as one decision-making methodology, namely, NG-AHP. Findings from the application revealed that in long-term strategy planning, steel manufacturing managers who are interested in agility, in the context of our study, should first deal with organisation management agility (C1) as the most significant criterion in the assessment followed by manufacturing process agility (C3), product design agility (C2), integration of information systems (C5), and partnership formation capability (C4), respectively. It was also concluded that steel company four (SC4) has the highest competence for agile strategies based on the five evaluation criteria followed by steel companies three (SC3), two (SC2), and one (SC1), respectively.
While this article offers contributions to the literature, it also has limitations which call for future research initiatives. Aspects of sustainability (i.e., economic, social, and environmental) to which our analysed criteria might have been more closely linked can be regarded as an interesting future research topic. In this study, however, we have considered the general definition of sustainability covering all three  pillars and discussed more closely connected sustainability dimensions based on the literature. e relation between lean manufacturing, agile manufacturing, and sustainability can be further explored in other manufacturing contexts to provide more insights on their relationships. In addition, recent mathematical developments in the realm of the NS theory can be applied in the MADM field to effectively capture uncertainty in future research such as α-discounting method for multicriteria decision-making [106]. For instance, comparing SVTNNs in the final weights of the N-AHP with no need to get the crisp values to make comparisons can help reduce information loss and reach a better evaluation. It would also be interesting to compare Pythagorean fuzzy AHP with the NG-AHP in numerical examples via simulation. As another future research direction, triangulation can be utilised by applying methods such as BWM, LBWA, or FUCOM to increase the validity of the proposed method.

Data Availability
e data used to support the findings of the study are made available by Amin Vafadarnikjoo on Mendeley Data under the licence CC BY 4.0 (https://doi.org/10.17632/ 8d6s67bgcw.1).

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.