Coal Refinery Process Absorbability Index Assessment against Foot Print of Air Pollution by Usage of Robust Optimization Algorithms: A Novel Green Environment Initiative

Coals are employed as fiery substance, and every day, millions of ton coal are consumed by coal users around the world. It is investigated that the millions of coal’s transportation/logistic till the coal user plants via road route and also inside the coal user plants (known as twice factors) not only enhance the air pollution but also cause the global warming. It is earlier known that coals emit the toxic pollutants and offensive gases such as sulfur dioxide, SO2; nitrogen oxides, NOX; hydrogen chloride arsenic; carbon monoxide, CO; methane; CH4; and CO2 on reacting with environmental O2 due to said twice factors, i.e., during the transportation from coal refinery spot to entry gate of coal user plants (another spot) and in process logistic/movement inside the coal user plants (loading to conveyor to coal fire tubes “attached with coal crushers”). Therefore, the coal refinery technique/process is found as the best practice to control air pollution under concerns of twice factors. The reliable and trustworthy coal refining technology improves the quality of coal by eradicating or eliminating the coating or layers of toxic particles from coal’s surface, which speedily crumble or decompose in reacting with environmental O2 under twice factors. As results, coal refining technology adds the green supply chain value into proposed twice factors and also save the world from breeding of ills and viruses. It is understood that the best coal refinery technique/process helps to overcome and reduce air pollution by responding discussed twice factors (accepted as research challenge and motivation of research). In the presented research work, the authors developed and proposed a dynamic multidimension Coal Refinery Process Absorbability Index (CRPAI) structure (consisted of coal refinery core dimension and subdimensions correspond to CRPA alternative techniques/processes) appended with Robust Optimization Algorithm (ROA) to be explored for opting the best CRP from available options. But due to inherent ambiguity, vagueness, and inconsistency involve in both dimensions of proposed structure, the assessment of expert’s panel is gathered in the terms of linguistic variable “appropriateness ratings” against the subdimensions of CRPAI structure corresponding to preferred CRP options. Next, assigned appropriateness ratings against the subdimensions are substituted by GIVFN. To arrive to core dimensions from subdimensions of CRPAI structure, a GITFN-OWGO (Ordered Weighted Geometric Operator) is investigated and modified as a Ordered Weighted Geometric Average Operator (OWGAO) to be applied for estimating the weights of subdimensions (core novelty of work). Finally, a ROA (consisted of MULTI-MOORA with dominance theory) is applied on the output of OWGAO for opting the viable and best CRP option. The positive effect of the dynamic multidimension CRPAI structure is that it helps the coal refinery companies to assess measure and evaluate the best and feasible coal refinery process under concern of twice factors using expert information. The research can be used to control the air pollution by responding aforesaid twice factors by single practice (the best coal refinery process/technique assessment and evaluation).


Introduction and State of Art
The CO 2 is ascertained as the foremost carbon intensive, which is accountable for more than 72% climate altering phenomenon and major source of air pollution and global warming vice versa. Today, coals are employed as burning material across the global coal user plants. The coal is found chemically intricate fuel as it librates flue gases, which might contain such as sulfur dioxide, SO 2 ; nitrogen oxides, NOX; hydrogen chloride arsenic; carbon monoxide, CO; methane; CH 4 ; and carbon dioxide, CO 2 [1][2][3][4][5][6][7][8][9][10]. Every day, millions of ton coal are fired by coal user plants. It is noticed that coal becomes extremely low contaminated/toxic gas emitter from coal refinery location/spot to entry gate of coal user plants and also throughout the logistic of coal inside the coal user plants if coal refinery is carried with best practices and technique at inception. Therefore, coal refinery is searched as the preeminent practice amongst many green practices to clean the environment and control air pollution. The coal refinery is respected as ecological balancer in the world, which function is to refine the coal to eliminate further air pollution ( [6,[11][12][13][14][15]; CURC 2004; [1-5, 7-10, 16, 17]).
Refined coals are the burnable goods, which are produced by application of a coal-upgrading practices-technology, which eliminate the moisture and certain pollutants from surface of coal stone and add the calorific values. Coal refining is precombustion treatments or processes, which changes the characteristics of coal prior to move is to be burned. In India nation, the coal refinery company use, such as dynamic process, is exhibited here: crushing and breaking-crushing and breaking devices are feeder breakers, rotary breakers, hammer mills, and roll crushers. Sizing them-it is sizing the coal under 20 to 45 mm diameters. Storage and stockpiling-coal is stored in silos or stockpiled before and after cleaning. Spherical, cylindrical, conical silos, or stockpiled is used. Density separation-coal is cleaned by separating the lower density organic material from the higher density by falling it from high 7-10 m height to earth. Washing-washing is done by quenching algorithm and passing through flowing water. Coal drying-thermal drying is lastly carried out by using container, drying in open overcast by sun to reduce the moisture content.
In investigating the CRPAI research documents (RDs), the authors found the many RDs in relation to analyzing the coal supplier performance and materializing the coal emission quantities around the coal user plants. So few RDs are found, dealt with mapping CRPAI under core dimension modeling by certain or quantitative data (except core dimension modeling by uncertain or qualitative data). In addition to that, few CRPA index documents are also grasped, dealt with CRPAI, but constrained to consider the core dimension (except both: core and subdimension simultaneously) [2-4, 7, 9, 16].
Furthermore, in investigating the RDs against ROAs, the authors determined that most of the Absorbability Index evaluation RDs dealt with the single optimization algorithm (could tackle certain data vs. core dimension only). The authors were not found still relevant research document, especially dealt with dynamic multidimension CRPAI struc-ture (consisted of coal refinery core dimension and subdimensions correspond to CRPA alternative processes) with ROAs (can tackle uncertain/qualitative data vs. both: core and subdimension only) to be used for opting the best CRP from available options. Therefore, the aforesaid research gaps motivated the authors to develop and propose a dynamic multidimension CRPAI structure (consisted of coal refinery core dimension and subdimensions) with ROA to opt the best CRP from available options under linguistic information of expert's panel.
To construct a dynamic multidimension CRPAI structure in addressing the qualitative modeling of core and subdimension simultaneously (CRPAI structure) and finding ROAs, the authors conduct a few literature surveys under aforesaid concerns: 1.1. A Literature over Framing the Dynamic Multidimension CRPAI Structure. [18][19][20][21], and [22] outlined the nature of the wastes; waste generating industries; waste characterization; health and environmental implications of waste management practices; steps towards planning; design and development of models for effective hazardous waste management; and treatment, approaches, and regulations for disposal of hazardous waste. Evaluation of the entire situation with reference to Indian scenario has been attempted in order that a better cost-effective strategy for waste management be evolved in the future. [23] proposed a decision model based on environmental performance indicators, which supported decision-making in supply chains in presence of environmental considerations. [24] implemented IVFN-TOPSIS (interval-valued fuzzy number-Technique for Order Preference Similar to Ideal Solution) approach in purpose to find out preference orders of manufacturing alternative firms in context of green supply chain management. [25] articulated an innovative, quantitative assessment of pollution avoidance attributable to environmental regulation enforced through integrated licensing, using Ireland's pharmaceutical-manufacturing sector as a case study. [16] argued the current coal consumption in the UK united with ecological impacts of coal mining, coal processing, and coal utilization for power generation. [7] shown impact of coal mining on land, water, health of workers, air, and social impact on surrounding villages and cities have evaluated. In addition to that, environment management plan of coal mining has also enlarged. [9] carried out a coal mining project in northwestern Bangladesh country in order to appraise the impacts of coal mining on the surrounding environment specifically on soil and water. The research of [26] was conducted to seek the miscellaneous paths to preserve the environmental and also put the environment out of fewer emission of carbon dioxide. The countless several conducted state of arts have been depicted as a citation in order that readers could seek the other way to diminish the universal warming issues in planet world ( [1,2,8,10,17,[27][28][29]; Sahu et al., 2012;Sahu et al., 2013;[7, 30]a, b; [9,22]).
Kalibatas and Turskis 2008 and [26] have also formed miscellaneous decision support system in context of manufacturing realm. [31] thermally treated the petroleum coke in the range of 1173-1673 K. The changes of the 2 Adsorption Science & Technology composition and crystalline structure of petroleum coke were investigated. ( [11,14,15]; ; [1-10, 13, 16, 17, 24, 30] a, b) worked to reduce air pollution by applying multivariable optimization methods. Venugopal et al. [32] stimulate the need for developing native solutions against the industrial challenges and revealed the importance of coal industry and coal washing process. Dhawan et al. [33] conducted an experiments for investigating the process parameters of catalytic coke pretreatment processes parameters. Sitorus et al. [34] applied a hybrid MCDM method to evaluate the coal mining and mineral processing units and discussed the pathways and future research directions of the proposed MCDM method. Pourabdollah et al. [35] studied the effect of thermal operation parameters against the graphitization of ethane-based catalytic coal. Dhawan et al. [36] investigated the usage of the catalyst K 2 CO 3 over raw coal and the use of KCO 3 over the raw coal as well as help to extract the CO 2 gas from raw coal. Shukla et al. [37] studied about the emission level of CO 2 in India from coal. Various advanced methods for capturing CO 2 are highlighted.
1.2. Literature over Framing the ROAs. [38] explored the application of multiobjective optimization on the basis of ratio analysis (MOORA) algorithm to solve different decision-making problems as frequently encountered in the real-time manufacturing environment. Six decision-making problems have been solved to prove the applicability, potentiality, and flexibility of this algorithm while solving various complex decision-making problems in present day manufacturing environment. [39] implemented the MOORA algorithm for solving multiple criteria (objective) optimization problem in milling process. Six decisionmaking problems which include selection of suitable milling process parameters in different milling processes are considered, and the obtained results almost match with those derived by the previous researchers which prove the applicability, potentiality, and flexibility of this algorithm in manufacturing environment subsequently. Sahu et al. 2012, Sahu et al. 2013, [30] a, b, [9,22] have conducted research in context of environmental subdimensions, which have been viewed as a case study. Sahu et al. [24] applied a weighted geometric aggregation operator on robot evaluation model to score the industrial value of robots. Bu et al. [40] proposed a GWmZd model consisted of green vendor evaluation parameters and solved it by knowledge-based grey-holistic technique with dominance theory for sustainable vendor evaluation. He et al. [41] proposed an advanced manufacturing machine tool evaluation model under objective information and solved it by grey relational analysis approach with dominance theory.
1.3. Summary of State of Art. As discussed, the authors only highlighted so few relevant literatures with brief description and cited more RDs to shorten research work. After literature survey, the authors consolidated the relevant CRPAI core and subdissensions and framed a dynamic multidimension CRPAI structure (consisted of coal refinery core dimension and subdimensions). In addition to that, the authors found a most significant ROA (MULTI-MOORA) and enable same to opt the best CRP from available options under linguistic information of expert's panel. But searched ROA was able for qualitative modeling of core dissensions. By using concept of fuzzy average rule, the authors made ROA able to tackle core and subdissensions simultaneously by introducing modification in the OWG operator (in previous researches, operator was able to act for qualitative modeling of only core level dimension). Moreover, on focusing over the most recent research of Bu et al. [40] and He et al. [41] for obtaining the reliable results, the authors applied the dominance theory for the comparison amongst the CRPAI preferences, provided by ROA.

Rationale to Conduct the Research
Coal pollutants mixed into the atmosphere, which disorder the ecosystem, are shown in Figure 1  Later, the authors preaudited the RDs of CRPAI and ascertained that there are still no RDs pertaining to CRPAI structure, which can deal with the qualitative modeling of both core and subdimensions. Up next, the authors found that there are a few RDs related to evaluation of CRPAI by executing single optimization algorithms (SOAs) with the OWG operator (operator was able to act for qualitative modeling of only core level dimension). The authors claimed that there is no RD yet, which can tackle subdimensions of CRPAI structure under qualitative modeling in application of the OWG operator appended with ROAs ( [30,[50][51][52]c, e). Therefore, the aforesaid rationale stimulated the authors to develop and propose a dynamic multidimension CRPAI structure (consisted of coal refinery core dimension and subdimensions) with constructing ROAs (MULTI-MOORA with modification of OWG operator by average fuzzy rule) with application of dominance theory (Bu et al., [40] and He et al., [41] for evaluation of the best CRP from available options under the linguistic information of expert's panel. This proposed research work would assist the managers to identify the best CRP practice and technique or to opt viable CRPs amongst preferred CRPs.

Rationale to Implement Interval-Valued Fuzzy Number Set (IVFNs) than Other Fuzzy Number Set
In the present context, authors fruitfully implemented IVFN except other fuzzy numbers, as IVFN deals with completely unknown information. In IVFNs, fuzzy membership function values are intended in between 0 and 1, aligned over

Fuzzy Preliminaries
Fuzzy sets and fuzzy logic are powerful mathematical tools, employed for modeling uncertain systems or qualitative information. A fuzzy set is an extension of a crisp set. A crisp set only allows full membership or nonmembership, while fuzzy sets allow the partial membership. The theoretical fundaments of fuzzy set theory are overviewed by [52,55,[57][58][59][60][61]. This section presents the generalized interval-valued trapezoidal fuzzy numbers. In addition, the arithmetic operations and aggregation of the generalized interval-valued trapezoidal fuzzy numbers are as follows.

The Generalized Trapezoidal Fuzzy Numbers. A fuzzy set
A in a universe of discourse X is characterized by a membership function μÃðxÞ which associates with each element x in X a real number in the interval [0, 1]. The function value μÃðxÞ is termed the grade of membership of x inÃ. A trapezoidal fuzzy number can be defined asÃ = ða 1 , a 2 , a 3 , a 4 ; wÃÞ as shown in Figure 2, and the membership function μÃðxÞ: R ⟶ ½0, 1 is defined as follows: Here, a 1 ≤ a 2 ≤ a 3 ≤ a 4 and wÃ ∈ ð0, 1Þ. Suppose thatã = ða 1 , a 2 , a 3 , a 4 ; wÃÞ andb = ðb 1 , b 2 , b 3 , b 4 ; wBÞ are two trapezoidal fuzzy numbers, then the operational rules of the trapezoidal fuzzy numbersã andb are shown as follows: [27] introduced the center of gravity (COG) measure for generalized trapez4oidal fuzzy numbers. Let there be a generalized trapezoidal fuzzy numberÃ = ða 1 , a 2 , a 3 , a 4 ; wÃÞ. Here, Assume that there are two IVFNs ẽ A and ẽ B, where From . The operation rules of interval-valued trapezoidal fuzzy numbers as given by [57] have been reproduced below.

Suppose that
and are the two interval-valued trapezoidal fuzzy numbers, where (i) The sum of two interval-valued trapezoidal fuzzy (ii) The difference of two interval-valued trapezoidal

Adsorption Science & Technology
(iii) The product of two interval-valued trapezoidal (iv) The product between an interval-valued trapezoidal fuzzy number and a constant λ ẽ A (v) The division between two interval-valued trapezoidal Here, and the operator "/" denotes exclusion of a certain term from sets U L andU U .
(vi) Rising to the power of a constant λ By considering (Equation (3)), we can define the COG point for an interval-valued trapezoidal fuzzy number ẽ (3) is employed to obtain the coordinates of the COG points for the lower and upper values of ẽ A, viz., ð , respectively. Secondly, the COG of ðx ẽ A , y ẽ A Þ is found as follows: Let there exist an interval-valued fuzzy number.
Adsorption Science & Technology respectively, are calculated by virtue of the Euclidean distance: The COG coordinates can also be employed when estimating the distance between two interval-valued trapezoidal fuzzy numbers say ẽ A and ẽ B [50,53].
Alternatively, one can employ the following technique (Liu and Jin, 2012):

The Generalized Interval-Valued Trapezoidal Fuzzy
Number Ordered Weighted Geometric Average Operator (OWGAO). [52] proposed an interesting and wellgrounded approach, named the ordered weighted average (OWA), which enabled to aggregate the variables in terms of their order in the set. Such an approach enables to avoid the subjectivity arising from group decision-making. [51] and later [24] introduced the generalized interval-valued trapezoidal fuzzy number ordered weighted geometric (GIFN-OWG) operator, which was able to tackle fuzzy ratings at 1 st level. But in the presented research work, the (GIFN-OWGA) is merged with average fuzzy rule and modified with renaming generalized interval-valued trapezoidal fuzzy number ordered weighted geometric average (GIFN-OWG) operator.

The Crisp MULTI-MOORA Algorithm
The Multiobjective Optimization by Ratio Analysis (MOORA) algorithm was introduced by [43], Sahu et al. [24], Bu et al. [40], and He et al. [41] and extended the algorithm to make it more robust as MULTI-MOORA (MOORA plus the full multiplicative form). The MOORA algorithm begins with matrix X where its elements x ij denote i th alternative of j th objective ði = 1, 2, ⋯, m ; j = 1, 2, ⋯, nÞ. The MOORA algorithm consists of two parts: the ratio system and the reference point approach. The MULTI-MOORA algorithm includes internal normalization and treats originally all the objectives equally important. In principle, all stakeholders interested in the issue only could give more importance to an objective. Therefore, they could either multiply the dimensionless number representing the response on an objective with a significance coefficient or they could decide beforehand to split an objective into different subobjectives.

The Ratio
System of MOORA. Ratio system defines data normalization by comparing alternative of an objective to all values of the objective: Here, x * ij denotes i th alternative of j th objective. Usually these numbers belong to the interval [0, 1]. These indicators are added (if desirable value of indicator is maximum) or subtracted (if desirable value is minimum); thus, the 7 Adsorption Science & Technology summarizing index of each alternative is derived in this way: Here, g = 1, ⋯, n denotes number of objectives to be maximized. Then, every ratio is given the rank: the higher the index, the higher the rank.
In some cases, it is often observed that some attributes are more important than the others. In order to give more importance to an attribute, it could be multiplied with its corresponding weight (significance coefficient) [38,42,62]. When these attribute weights are taken into consideration, Equation (28) becomes as follows: Here, w j is the weight of j th subdimensions.

The Reference Point of MOORA.
Reference point approach is based on the ratio system. The Maximal Objective Reference Point (vector) is found according to ratios found by employing Eq. The j th coordinate of the reference point can be described as (r j = max x * ij ) in the case of maximization. Every coordinate of this vector represents maximum or minimum of certain objective (indicator). Then, every element of normalized response matrix is recalculated, and the final rank is given according to deviation from the reference point and the Min-Max Metric of Tchebycheff: 5.3. The Full Multiplicative Form and MULTI-MOORA. [43] proposed MOORA to be updated by the full multiplicative form algorithm embodying maximization as well as minimization of purely multiplicative utility function. Overall utility of the i th alternative can be expressed as dimensionless number: Here, A i = Q g j=1 x ij ; i = 1, 2, ⋯, m denotes the product of objectives of the i th alternative to be maximized with g = 1, 2, ⋯, n being the number of objectives to be maximized and where B i = Q n j=g+1 x ij ; i = 1, 2, ⋯, m denotes the product of objectives of the i th alternative to be minimized with n − g being the number of objectives (indicators) to be minimized. Thus, MULTI-MOORA summarizes MOORA (i.e., ratio system and reference point) and the full multiplicative form.

MULTI-MOORA Algorithm Based upon Interval-Valued
Trapezoidal Fuzzy Numbers. Let E = 1, 2, ⋯, e denote the e th expert involved in a decision-making process. Suppose that the experts provide ratings for each i th alternative against each j th criterion with i = 1, 2, ⋯, m and j = 1, 2, ⋯, n. The set of criteria can be split into two subsets, namely, those of cost criteria, C, and benefit criteria, B. Cost criteria are to be minimized, whereas benefit criteria are to be maximized. Each criterion can be attributed with respective weight W j , such that W j ≥ 0, and ∑ j W j = 1.
Step 1. Each of decision-makers constructs his own decision ; w U ijk Þ being responses of alternatives on criteria.

Adsorption Science & Technology
Step 2. Individual decision matrices are aggregated by employing the GITFNOWGA operator.
Here, w k is the weight of the k th largest response obtained by Equation (26) and It must be noted that Equation (21) is employed to com-pare the values of the generalized interval-valued trapezoidal fuzzy numbers.
Step 3. In case some of criteria involve numeric data, the normalization has to be carried out. Here, , 3, 4g for ∀j = 1, 2, ⋯, n.
Step 4. The ratio system: Effective primary breaking with low gas emission, ξ 11 Effective final breaking with low gas emission, ξ 12 Absorbability index on sizing process, ξ 2 Effective machining for sizing with low gas emission, ξ 21 Absorbability index on storage and stockpiling process and density separation, ξ 3 Toxic gas emission on the stockpiling of coal, ξ 31 Emission of toxic gases on density separation, ξ 32 Absorbability index on density separation and washing process, ξ 4 Emission of toxic gases on quenching process, ξ 41 Water prequenching treatment level, ξ 42 Water final-quenching treatment level, ξ 43 Absorbability index on drying process, ξ 5 Moisture drying pretreatment level, ξ 42 , ξ 51 Moisture drying final treatment level, ξ 51 Table 2: Nine-member linguistic terms (rating scale) and their corresponding interval-valued fuzzy numbers. The normalized values are added up for the benefit criteria and subtracted for the cost criteria: Here, RS i denotes the overall utility of the i th alternative in terms of the ratio system. The alternatives are then ranked by measuring their distances from the origin point in the spirit of Equation (21). Specially, alternatives with higher distances receive higher ranks.
Step 5. The reference point approach: For the sake of convenience, one can employ the Maximal Utopian Reference Point (MURP), rather than the Maximal Objective Reference Point. In the case of the generalized interval-valued trapezoidal fuzzy numbers, MURP is defined as follows: Therefore, Equations (23) and (24) can be utilized to identify the maximal deviation from the MURP for each alternative: Then, the alternatives are ranked by minimizing the maximal deviances found in Equation (36).
Step 6. The full multiplicative form: The fuzzy utility of the i th alternative is obtained by employing Equations (12) and (14).

Empirical Research
In the presented research forum, the authors developed a standard dynamic multidimension Coal Refinery Process Absorbability Index (CRPAI) structure by conducting the real CRP audit and the prior state of arts and proposed the same to CRP companies for finding the best CRP technique/practice to control further air pollution-global warming under discussed twice factors, shown in Table 1.
To simulate the CRPAI, the authors used modified OWGA operator with application of the MULTI-MOORA algorithm embedded with dominance theory. To validate the CRPAI, the authors solved a problem of a CRP company, existing north part of India. The CRP company wanted to evaluate and finalize a reliable CRP alternative amongst available (θ 1 , θ 2 , θ 3 , and θ 4 ) for controlling further air pollution, occurs due to discussed twice factors. For data modeling of the subdimensions of CRPAI, a board of five highly experience decision-makers E 1 , E 2, E 3, E 4 , and E 5 are evaluated from the inventory, production, and environmental departments of CRP (case study) company. Next, the expert's panel is invited to sight, audit, and assess the   capability of alternative CRPs such as θ 1 , θ 2 , θ 3 , and θ 4 with high concern about the discussed twice factors. After that, decision-makers are assisted with linguistic variable scale, expressed 1-9 IVFN scale as shown in Table 2, so that decision-makers can assign appropriateness ratings against individual subdimensions.
The procedural steps are summarized below: Step 1. Gathering information from the expert group in relation to performance ratings of different evaluation dimensions and subdimensions using linguistic terms: In order to acquire appropriateness ratings against subdimensions, evaluated board expert's panel, e.g., E 1 , E 2 , E 3 , E 4 , and E 5 , instructed to assign appropriateness ratings (evaluation score) against subdimensions (associated in proposed structure), which is further transformed into IV-fuzzy number. The appropriateness ratings (in linguistic terms) against evaluated subdimensions assigned by the expert panel corresponding to preferred options, e.g., θ 1 , θ 2 , θ 3 , and θ 4 , have depicted in Table 3.
Step 2. Approximation of the linguistic evaluation information by IV trapezoidal fuzzy numbers: The assigned appropriateness ratings acquired in the form of IVFNs against subdimensions, shown in Table 3, are shifted into 1 st level (core dimensions) by using modified GITFN-OWGAO (generalized interval-valued trapezoidal fuzzy number ordered weighted geometric average operator); In Equations (25) and (26), consequently, the problem appeared as structure modeling, revealed in Table 4.
Step 3. Normalization: After forming structure matrix, the IV fuzzy set value against core dimensions is found out in between [0, 1] and core dimensions were beneficial in nature, Hence, structural matrix was not normalized with exploration of Equation (34).
Step 4. The ratio system: In the ratio system (RS), the values of IVFNs against core dimensions corresponding to CRP alternative or techniques are added up for the benefit criteria and subtracted for the cost criteria with exploration of Equation (35); the results are shown in Table 5. The higher scores are prioritized corresponds to higher preference orders.
Step 5. Reference point approach: We define the reference point (ascertained the core dimensions are benefit in nature) Thus, the preference order of preferred CRPs is decided in accordance to their distances from fuzzy reference points and is computed by using Equations (36) and (37). Smaller distances of subdimensions from fuzzy reference points are preferred first, shown in Table 6.
Step 6. The full multiplicative form: Eventually, Equation (38) is utilized to obtain preference orders of CRP options in accordance to full multiplicative form of MOORA, and the results are shown in Table 7.
Step 7. Final ranking orders by exploring of dominance theory: By passing through respective computational algorithms of MULTI-MOORA: ratio system analysis, reference point, and full multiplicative form, finally the dominance theory ( [45]; [40], [41]) is implemented to summarize the preference orders provided by three respective parts of MULTI-MOORA computational algorithms Bu et al. [40] and He et al. [41]. Hence, the final preference orders of preferred CRP options are presented in Table 8.
As results, the fourth CRP practice/technique/process ð θ 4 Þ is found more viable amongst preferred CRP practices/techniques/processes as per the subjective report of expert panel, whereas the first CRP practice/technique/process ðθ 1 Þ and subsequently ðθ 3 Þ are found as a second and third best  It was visibly found that ðθ 4 Þ CRP technique/process would be the best to diminish the air pollution, growing universal warming and preserving the ecological balance caused by CRPs said by ([7, 19, 20, 22, 28, 48, 63]; Langergraber).

Novelties of Conducted Research
(i) The authors served a research clue in the front of CRP researchers to overcome global warming issue by addressing the discussed twice factors (ii) The authors proposed and developed CRPAI structure (consisted of coal refinery core dimension and subdimensions corresponding to CRP alternative) to diminish the air pollution (iii) This research dealt with comprehensive qualitative modeling of dimension and subdimensions corresponding to CRP alternative using expert information (iv) These authors adapted and applied the three qualitative data optimization algorithms concurrently to robustly opt the viable CRP [40,41] (v) The authors modified the IVFN-OWGA (intervalvalued fuzzy number ordered weighted geometric average) operator by introducing the average fuzzy rule, depicted in Equations (25) and (26), and consequently, ROAs were made able to tackle core and subdissensions simultaneously by applying the modified OWGA operator (operator was able to act for qualitative modeling of only core level dimensions) (vi) The authors applied dominance theory to conduct robust comparison amongst the CRPAI preferences, received by the three parts of the MUTI-MOORA algorithm [40,41].
(vii) The modified OWGA operator can be explored in future to tackle the extended or advanced CRP core-subdimensions

Conclusions
The coal is determined as a fossil fuel, which releases the several offensive monoparticles and high magnitude of CO 2 to the atmosphere and also causes acid rain and high temperature ( [6,[11][12][13]; CURC 2004; [1,20]). It is well purported by the ecological scholars that CRP contributes in controlling the air pollution [2,7,18,19,22,26,28,48,63]. Therefore, to control such as air pollution and global warming vice versa, the presented research work enrolled a dynamic multidimension CRPAI structure accompanied with the ROA (modification of OWGAO appended with MULTI-MOORA with application of dominance theory) algorithm to scrutinize and benchmark the best CRP from available options under linguistic information of expert's panel. This proposed research work assisted the managers to identify the best CRP practice and technique or to opt viable CRPs amongst preferred CRPs that will reduce the air pollution, consequently weaken the virus and ills, and make living people happy. For making decision, manager can use the proposed dynamic multidimension CRPAI structure and invite the evaluated team of expert of own and other industries to view the preferred CRP alternatives one by one. Next, the expert would ask to provide the ratings against subdimensions of CRPAI structure by using proposed linguistic variable corresponding to IVFN. Next, assigned rating would be substituted by IVFN by using OWGOA. Later, the MUTI-MOORA algorithm would deliver results for further comparison by dominance theory and making final decision. This proposed research work assists the managers of CRP to acquire the many coal supplying orders from coal user plants (which strictly follow the green policies). The proposed research dossier helps to finalize the best CRPAI technique or CRPs amongst available for winning the discussed twice factors such as coal logistic from coal refinery location to entry gate of coal user power plants and in process logistic/movement inside the power plants, linked to diminish the air pollution. The wok cannot shape other problems such as liner programming and route algorithm.

Data Availability
The data used to support the findings of this study are available in Tables 5-7.