Understanding the Mechanism of Abrasive-Based Finishing Processes Using Mathematical Modeling and Numerical Simulation

Recent advances in technology and refinement of available computational resources paved the way for the extensive use of computers to model and simulate complex real-world problems difficult to solve analytically. The appeal of simulations lies in the ability to predict the significance of a change to the system under study. The simulated results can be of great benefit in predicting various behaviors, such as the wind pattern in a particular region, the ability of a material to withstand a dynamic load, or even the behavior of a workpiece under a particular type of machining. This paper deals with the mathematical modeling and simulation techniques used in abrasive-based machining processes such as abrasive flow machining (AFM), magnetic-based finishing processes, i.e., magnetic abrasive finishing (MAF) process, magnetorheological finishing (MRF) process, and ball-end type magnetorheological finishing process (BEMRF). The paper also aims to highlight the advances and obstacles associated with these techniques and their applications in flow machining. This study contributes the better understanding by examining the available modeling and simulation techniques such as Molecular Dynamic Simulation (MDS), Computational Fluid Dynamics (CFD), Finite Element Method (FEM), Discrete Element Method (DEM), Multivariable Regression Analysis (MVRA), Artificial Neural Network (ANN), Response Surface Analysis (RSA), Stochastic Modeling and Simulation by Data Dependent System (DDS). Among these methods, CFD and FEM can be performed with the available commercial software, while DEM and MDS performed using the computer programming-based platform, i.e., “LAMMPS Molecular Dynamics Simulator,” or C, C++, or Python programming, and these methods seem more promising techniques for modeling and simulation of loose abrasive-based machining processes. The other four methods (MVRA, ANN, RSA, and DDS) are experimental and based on statistical approaches that can be used for mathematical modeling of loose abrasive-based machining processes. Additionally, it suggests areas for further investigation and offers a priceless bibliography of earlier studies on the modeling and simulation techniques for abrasive-based machining processes. Researchers studying mathematical modeling of various microand nanofinishing techniques for different applications may find this review article to


Introduction
Metal cutting is considered the most significant procedure where the material removal process is performed by more than one edge to fulfill some of the specific criteria needed to shape the final workpiece in the form of particular geometrical dimensions and surface finish. Even though there has occurred a lot of advancement in the development of machining tools and process controlling and monitoring, there still exists room for improvement in metal cutting processes and procedures for which significant research can be conducted by various metal-forming industries [1]. This may occur due to the optimizing procedures' intricacy, particularly when the materials with higher "mechanical characteristics" are developed and the machinability has not yet been inspected fully [2].
In addition to the material removal and metal cutting process, there is a fine requirement that product surface quality is considerably better since they play a major role in aerospace, automotive, and biomedical industries [3]. Even a small scratch or an uneven burr can cause huge damage to the engine. It may fail the aerospace gadgets or malfunction any of its components, etc. These industries are trying their best to spend a huge amount on the finishing of such components in order to make components that are cutting/burrs marks-free [4]. For the last 40 years, various manufacturing industries have been using conventional procedures for finishing the components, such as grinding, honing, lapping, etc., to get the machining components' desired finishing. However, these conventional procedures of finishing are restricted to very few geometries and cannot work on complex and intricate geometries as well as complicated profiles for finishing of high level, which is required while the operation of the component is in process. These limitations and restrictions in the finishing process have led the industries to develop advanced finishing procedures, known as "Abrasive flow machining (AFM)" [5]. AFM process is used to perform the finishing process of the machining components and its internal features in several engineering materials such as alloys, ceramics, non-ferrous, super-alloys, refractory materials, semiconductors, carbides, quartz, and various other composites, etc., that are considered to be difficult to be done using the traditional processes economically and proficiently [6]. The significance of this AFM process is that it generates the nano-level finishing of the machining components that are essentially desirable at this time. The concept of the AFM finishing procedure was first given by the "Extrude Hone corporation of the USA" in 1960 for the finishing of aerospace components to achieve the desired accuracy. These days, the AFM process is the best technique for finishing intricate geometries, which were never achievable by any traditional finishing tools. Even a lot of study has been conducted by researchers in order to improve the finishing process performed by the AFM [7]. The AFM process efficiently improves the surface finish of the complicated geometry like gears, trim-dies, turbine blades, bio-medical implants, etc. [8,9]. Many researchers have highlighted the potential of the AFM process to finish any complex and freeform surfaces [10,11]. Researchers have also developed a hybrid form of AFM processes, i.e., rotational type abrasive flow machining (R-AFM), to finish the asymmetrical complicated workpiece [12]. Magnetorheological abrasive flow finishing (MRAFF) has also been developed under the umbrella of the hybrid AFM process to finish optical surfaces [13]. Many researchers have employed the AFM and its hybrid variants for surface finishing of metal matrix composite materials, i.e., Al/SiC MMCs [14,15]. Hybrid manufacturing processes, i.e., ultrasonic-assisted machining processes, are also a good solution due to less cutting force required during the machining operation [16]. Ultrasonic-assisted magnetic abrasive finishing (UAMAF) significantly improves the surface roughness as well as the hardness of the workpiece [17,18]. Centrifugal force-assisted abrasive flow machining (CFAFM) is an efficient hybrid machining process to finish cylindrical surfaces with better surface roughness [19].

Loose Abrasive-Based Machining Processes
The loose AFM process is considered the most significant approach to finishing the intricate and complex geometries and advanced materials used within industrial practices [20,21]. It is considered to be the recently developed "surface finishing process" that has been used widely in various industrial applications such as defense, aerospace, biomedical, tool and die, etc. [22]. In loose AFM, no structure exists that connects or links the grains together in the matrix in bonded "abrasive process." The various standard processes in Loose abrasive processes include AFM, polishing, lapping honing, and abrasive blasting. AFM is also one of the finishing techniques that entail the oil for the particle of abrasive and carrier polymer and assists in eliminating the barriers to acquiring the desired size and shape [14]. Thus, any complicated shape can be easily finished by Abrasive Flow Machining abbreviated AFM. Abrasive flow machining helped us to tackle the shortcomings of conventional processes like "grinding," "lapping," and "honing" [20]. The problem associated with AFM is its high cost and "low material removal rate MRR," which makes it laborious and strenuous. AFM processes of finishing techniques have also emerged recently to tackle the challenges and issues such as processes requiring high machining time [23][24][25][26]. Figure 1 below illustrates the "loose abrasive-based machining" (LAM) process. abrasive flow machining (CFAFM) is an efficient hybrid machining process to finish cylindrical surfaces with better surface roughness [19].

Loose Abrasive-Based Machining Processes
The loose AFM process is considered the most significant approach to finishing the intricate and complex geometries and advanced materials used within industrial practices [20,21]. It is considered to be the recently developed "surface finishing process" that has been used widely in various industrial applications such as defense, aerospace, bio-medical, tool and die, etc. [22]. In loose AFM, no structure exists that connects or links the grains together in the matrix in bonded "abrasive process." The various standard processes in Loose abrasive processes include AFM, polishing, lapping honing, and abrasive blasting. AFM is also one of the finishing techniques that entail the oil for the particle of abrasive and carrier polymer and assists in eliminating the barriers to acquiring the desired size and shape [14]. Thus, any complicated shape can be easily finished by Abrasive Flow Machining abbreviated AFM. Abrasive flow machining helped us to tackle the shortcomings of conventional processes like "grinding," "lapping," and "honing" [20]. The problem associated with AFM is its high cost and "low material removal rate MRR," which makes it laborious and strenuous. AFM processes of finishing techniques have also emerged recently to tackle the challenges and issues such as processes requiring high machining time [23][24][25][26]. Figure 1 below illustrates the "loose abrasive-based machining" (LAM) process.    The machining accuracy of different traditional and advanced machining processes is described in Figure 3.
A summary of loose abrasive-based machining processes about their basic principle and process parameters is shown in Table 1.
The comparison of surface roughness of traditional finishing techniques s shown in Figure 4 below. The machining accuracy of different traditional and advanced machining processes is described in Figure 3.  [27].
A summary of loose abrasive-based machining processes about their basic principle and process parameters is shown in Table 1. The machining accuracy of different traditional and advanced machining processes is described in Figure 3.  [27].
A summary of loose abrasive-based machining processes about their basic principle and process parameters is shown in Table 1.  The comparison of surface roughness of advanced finishing techniques, i.e., nan ishing techniques, is shown in Figure 5 below.  Table 1. Literature summary of loose abrasive-based Machining process with basic principle and process parameters.

S. No.
Loose Abrasive-Based Method Typical Materials Principle Average Surface Roughness (Ra µm) Application Ref.

Traditional Finishing Processes
Lapping Metals Harden the trapped partials between the surface of the workpiece and a soft counter formal surface 0.05 µm to 2.5 µm Improving the surface finish in loose abrasive-based machining [28] Buffing and Polishing Metal and wood Finishing the surface is performed through a wheel or an abrasive 0.1 µm to 0.41 µm Smoothing the material [29] Super abrasive machining Titanium, nickel-based alloys and metal matrix composites, etc.

Abrasive Flow Machining (AFM)
Hardened steel Removing material from the surface of the material using sliding of abrasive particles of another material 0.184 µm Can deal with complex components [26,28,30] Elastic Emission Machining (EEM) silicon removing the material atom by atom with no cracks, deformation, or deep indentations Less than 0.0005 µm Removing material in loose abrasive-based machining [28,30,31] Chemical Mechanical Polishing (CMP) silicon A layer is formed between the workpiece and the slurry by chemical reaction. This layer is softer than the material of the original workpiece, which can be easily removed.
0.005 µm to 0.01 µm used in the finishing of the wafers made from silicon [28,30,32] Magnetic Abrasive Finishing (MAF) silicon nitride In this process, ferromagnetic particles are mixed with the particles of the workpiece. This mixture is brought close to the surface of the workpiece to be finished.
Less than 0.01 µm Finishing the surface of nano-chips [28,30] Magnetorheological Finishing (MRF) Flat BG7 glass To finish the surface in this process, magnetorheological fluid (MR fluid) is used. This fluid consists of iron particles, carbonyl iron particles (CIPs), carrier fluid, and abrasive particles. The viscosity of this fluid is increased when it comes under the effect of the magnetic field. The CIPS arrange along the magnetic force line, and then the abrasive particles are entangled within the chains, and their motion causes the material removal Less than 0.001 µm Finishing micro/Nano-chips [28,30] Magnetorheological Abrasive Flow Finishing (MRAFF) This method is an extended version of MRF. the extrusion of the medium of the magnetorheological polishing is performed through the surface of the workpiece Less than 0.001 µm Finishing softer material than MRF [30,33,34] Ball end Magnetorheological Finishing (BEMRF) Metal mirror and glass for the finishing, this process uses the rotating spot of the stiffened MR fluid Less than 0.001 µm Used in the finishing of the complex 3d geometry [35] Metals 2022, 12, 1328 6 of 65 The comparison of surface roughness of advanced finishing techniques, i.e., nano finishing techniques, is shown in Figure 5 below. The comparison of surface roughness of advanced finishing techniques, i.e., nano finishing techniques, is shown in Figure 5 below.

Abrasive Flow Machining (AFM)
AFM is considered to be the non-conventional process of finishing that deburrs, polishes, alleviates the recast layers, and has the potential to generate the "compressive residual stresses" at the part where it is difficult to reach with the conventional finishing process. Few of the major applications of AFM include the finishing of various aerospace and medical components, automotive parts, and high-volume generation of electronic parts. Several materials such as soft aluminum, ceramics, tough nickel, and carbides have been micro-machined successfully through this process [36]. The process of AFM gives the highlevel surface finish as well as close geometric tolerance at reasonable and economical rates for a significant range of industrial parts and components [37]. One of the basic reasons behind implementing the AFM process is the ability of its media to finish the intricate points in the components where it is difficult to reach using conventional methods. Also, it can follow the "complex contours" and work simultaneously on several edges and surfaces, thus making it more desirable in comparison to the other finishing processes [38]. In the manufacturing process, several mechanical parts undergo one machining operation at least once in order to fulfill a certain function within the environment [39]. Some criteria for the quality of such pieces are desired, such as dimensional, roughness, as well as geometric specifications, etc.
AFM alleviates small quantities of the material through "semi-solid abrasive laden media" across the workpiece. In this process, two vertical cylinders opposing each other extrude the media in both directions (back and forth) across the passage in the workpiece as well as tooling, as shown in Figure 6. Here, the media comprises the abrasive grains and that of the semi-solid carrier, where the media acts as a "self-deformable stone" possessing protruding abrasive particles that act as the cutting tools.
ness, as well as geometric specifications, etc.
AFM alleviates small quantities of the material through "semi-solid abrasive laden media" across the workpiece. In this process, two vertical cylinders opposing each other extrude the media in both directions (back and forth) across the passage in the workpiece as well as tooling, as shown in Figure 6. Here, the media comprises the abrasive grains and that of the semi-solid carrier, where the media acts as a "self-deformable stone" possessing protruding abrasive particles that act as the cutting tools.

Developments in AFM Process
The abrasive flow machining process can provide good results for the machining components that need the removal of the defects created by the mechanical or

Developments in AFM Process
The abrasive flow machining process can provide good results for the machining components that need the removal of the defects created by the mechanical or manufacturing processes. The AFM process makes ideal for relief of surface stress, radiusing, polishing, deburring and geometric optimization [41]. Being an essential requirement of industry these days, a large amount of research work has been reported in this domain in literature.
AFM is considered to generate the "self-deforming tool" that accurately removes the extra workpiece material and then finishes the surface from the areas where it is difficult to do so with the conventional machining process [42].
The material removal process comprises three modes of deformation, which are as follows: • "Elastic deformation": it correlates with rubbing • "Plastic deformation (ploughing)": it correlates with the material being displaced without being alleviated. • "Micro-cutting" [43][44][45] Initially, in [46], the authors developed the one-way process of AFM where the medium traveled in a single direction and termed it to be the simplest type of AFM that consumes the least time during the finishing process [47][48][49][50]. It was observed that surface roughness could be minimized to 90% on the machined, cast, or EDM surface having the "dimensional tolerance" up to or ±0.005 mm [51,52].
In order to get a better radiusing and finishing action not only in the inner but outer surface as well as the components, two way AFM process of finishing was developed [53,54]. Its basic operation involves two horizontal or vertical hydraulic cylinders that lie opposite each other and are placed in between the work piece with the help of suitable fixtures, as illustrated in Figure 7. The pistons in the cylinders are used to make the "abrasive laden medium" move back and forth on the surface that needs to be finished. In this way, the finishing process is carried out in a two-way AFM operation.
In order to get a better radiusing and finishing action not only in the inner but outer surface as well as the components, two way AFM process of finishing was developed [53,54]. Its basic operation involves two horizontal or vertical hydraulic cylinders that lie opposite each other and are placed in between the work piece with the help of suitable fixtures, as illustrated in Figure 7. The pistons in the cylinders are used to make the "abrasive laden medium" move back and forth on the surface that needs to be finished. In this way, the finishing process is carried out in a two-way AFM operation. In [55], another development was made regarding the oscillatory motion of the workpiece. With oscillations to keep going, the work piece hits the "abrasive medium" with the "eccentric path" that causes the intricate and complex structure to make an interaction fully and completely with that of the abrasive medium, which results in the same (equal) abrasion on each side of the components. In [56], an AFM setup has been developed on the lathe to carry the experiment on brass as well as aluminum. It was observed that the dominant factors are the "abrasive concentration" in the medium, followed by the mesh size, speed of the medium flow, and the number of cycles. A lot of improvement could be In [55], another development was made regarding the oscillatory motion of the workpiece. With oscillations to keep going, the work piece hits the "abrasive medium" with the "eccentric path" that causes the intricate and complex structure to make an interaction fully and completely with that of the abrasive medium, which results in the same (equal) abrasion on each side of the components. In [56], an AFM setup has been developed on the lathe to carry the experiment on brass as well as aluminum. It was observed that the dominant factors are the "abrasive concentration" in the medium, followed by the mesh size, speed of the medium flow, and the number of cycles. A lot of improvement could be seen in the case of materials that were soft, considering the surface roughness enhancement. In [57], the temperature has been considered to be the most influential parameter in the AFM process regarding work efficiency; with the escalation in no. of cycles, the medium temperature increases, which means the decrease in the medium viscosity occurs with the increase in no. of cycles. In work, it has been presented that in AFM testing, an increase in no. of cycles significantly decreases the material removal rate (MRR) as well as surface roughness, thus reducing the efficiency [58].
In order to enhance the process performance and control the rheological properties, scientists and researchers have designed and developed several AFM process variants [59,60] by combining and incorporating the basic AFM operation with that of the traditional/nontraditional flow machining operation. The process of AFM has been classified as per the energy and forces used within the process [61].
Other than the above-mentioned four types of AFM process, there have also been developed some other AFM variants in order to compensate for the MRR and also for a longer finishing time. The variants of the AFM process have few external assistants, such as ultrasonic vibrations [62][63][64], the rotational effect [65][66][67], magnetic field assistance [39,68,69], and also the hybrid variant of the AFM process [54,[70][71][72]. The external assistance is found to show the increase in the finishing force [73][74][75], abrasive velocity, enhancement in the contact length among the work surface as well as abrasive, active density of abrasive, or utilization of the synergic effect of two processes in order to eradicate the material [76]. These effects result in enhancing the finishing time, MRR [62,63], and surface roughness. The further classification of the developed AFM variants has been shown in Figure 8 below.
A summary of the differences and similarities of the different AFM process variants is shown in Table 2.
tance is found to show the increase in the finishing force [73][74][75], abrasive velocity, enhancement in the contact length among the work surface as well as abrasive, active density of abrasive, or utilization of the synergic effect of two processes in order to eradicate the material [76]. These effects result in enhancing the finishing time, MRR [62,63], and surface roughness. The further classification of the developed AFM variants has been shown in Figure 8 below. A summary of the differences and similarities of the different AFM process variants is shown in Table 2.  In order to give the best control on the AFM processes, the selection of the parameters of the process is necessarily required. This review paper has segmented the AFM process parameters into three main categories. The classification of these three main AFM parameters, including machine, medium, and workpiece [77][78][79]. Here, the machine decides the level to which the abrasion can be processed through "viable extrusion pressure," flow rate, flow volume, and also the number of cycles [80,81]. "Rheological properties" of "abrasive laden medium" are the ones that signify the amount of abrasion. The grit size, viscosity, temperature, and polymer to dilute are considered to be the parameters for comprehending the rheological properties of the MR fluid. Next comes the workpiece that is finished by the AFM process. The properties of the workpiece, such as material type, roughness value, and also geometry, tell the machining time and abrasive type to be used [82].
Summary of the performance analysis of the AFM process regarding their process parameters is shown in Table 3. Table 3. Literature summary of performance analysis of the AFM process regarding their process parameters.

S. No.
Ref The typical defects during different manufacturing processes for electronic components, pipes, welding parts, and textile materials can be improved using the abrasive-based machining process, i.e., the AFM process and its variant processes, as shown in Figure 9. The final surface roughness was measured to be 35-78 nm.
The typical defects during different manufacturing processes for electronic components, pipes, welding parts, and textile materials can be improved using the abrasivebased machining process, i.e., the AFM process and its variant processes, as shown in Figure 9. Hashmi et al. also have discussed the various surface defects occurs during additive manufacturing processes i.e., staircasing/stair-stepping effects. Authors have discussed the various post-processing techniques including abrasive-based machining methods for removing the surface defects of these additively manufactured parts. Authors have also discussed the machine vision technique or artificial intelligence-based technique for meas- Hashmi et al. also have discussed the various surface defects occurs during additive manufacturing processes i.e., staircasing/stair-stepping effects. Authors have discussed the various post-processing techniques including abrasive-based machining methods for removing the surface defects of these additively manufactured parts. Authors have also discussed the machine vision technique or artificial intelligence-based technique for measuring the parameters of machining processes [88][89][90][91][92][93]. Authors have also discussed the experimental investigation of the AFM process to improve the surface finish of polymer-based AM parts [94][95][96][97].
A summary of the process performance analysis of abrasive-based machining processes in different applications is shown in Table 4. In the present research work, different research related to mathematical modeling and simulation methods for abrasive-based machining processes is categorized to provide a valuable study for the researchers in the exciting field of computational analysis. As a result, new ideas in the mathematical modeling and simulation methods and research gaps in the existing literature are obtained. Future research works are also suggested to push forward this interesting research field. Section 1 presents a review of research works related to the general overview of abrasive-based machining processes and their types. In Section 2, research works are classified according to the different topics in research to the mathematical modeling and simulation approaches for abrasive-based machining processes. In Section 3, research works are classified according to the various topics in research to the simulation techniques for abrasive-based machining processes. In Section 4, research works are classified according to the various topics in research to the mathematical modeling for abrasive-based machining processes. In Section 5, the conclusions of this literature review have been discussed, and future research works in the field of mathematical modeling and simulation of abrasive-based machining processes are also suggested. This review paper provides deep insight into the current developments in AFM processes, abrasive media, their modeling as well as optimization, and various promising application domains.
Since the paper mainly focuses on the modeling and simulation techniques of abrasivebased machining, the next sections describe various other modeling types and simulation processes of AFM.

Mathematical Modeling and Simulation Approaches for Abrasive Machining Processes
Recently, computer simulation applications have proved their efficiency in solving complex problems in Engineering. Different techniques of modeling and simulations are considered widely in the Industries, such as "computational fluid dynamics (CFD)" [83,104], "discrete element method (DEM)," "finite element method (FEM)," "artificial neural network (ANN)" [14], "multi-variable regression analysis (MVRA)", "response surface analysis (RSA)" and "Stochastic modeling (SM)" [105]. These techniques have several applications in the industrial process, such as predicting the behavior of the material removal in loose abrasive-based machining [106]. The results obtained from the modeling and simulations pave how to optimize the process parameters that affect the finish of the workpiece. Several researchers have prominently worked on the "mathematical modeling and simulation" of "loose abrasive-based finishing processes." Several experimental as well as computational techniques and approaches have been performed to analyze the ideal and optimal procedural parameters [103,[107][108][109].

Classification of Modeling and Simulation Techniques for Abrasive-Based Machining Processes
The modeling and simulation-based techniques for LAM may be classified based on the used techniques, i.e., "computational techniques and statistical techniques Petare et al. have also classified several simulation techniques [110], which are shown in Figure 10. processes, abrasive media, their modeling as well as optimization, and various promising application domains.
Since the paper mainly focuses on the modeling and simulation techniques of abrasive-based machining, the next sections describe various other modeling types and simulation processes of AFM.

Mathematical Modeling and Simulation Approaches for Abrasive Machining Processes
Recently, computer simulation applications have proved their efficiency in solving complex problems in Engineering. Different techniques of modeling and simulations are considered widely in the Industries, such as "computational fluid dynamics (CFD)" [83,104], "discrete element method (DEM)," "finite element method (FEM)," "artificial neural network (ANN)" [14], "multi-variable regression analysis (MVRA)", "response surface analysis (RSA)" and "Stochastic modeling (SM)" [105]. These techniques have several applications in the industrial process, such as predicting the behavior of the material removal in loose abrasive-based machining [106]. The results obtained from the modeling and simulations pave how to optimize the process parameters that affect the finish of the workpiece. Several researchers have prominently worked on the "mathematical modeling and simulation" of "loose abrasive-based finishing processes." Several experimental as well as computational techniques and approaches have been performed to analyze the ideal and optimal procedural parameters [103,[107][108][109].

Classification of Modeling and Simulation Techniques for Abrasive-Based Machining Processes
The modeling and simulation-based techniques for LAM may be classified based on the used techniques, i.e., "computational techniques and statistical techniques Petare et al. have also classified several simulation techniques [110], which are shown in Figure 10. Based on the classification presented in Figure 9 above, the detailed description of each classification has been elaborated in the next sections.
A summary of the shortcomings of simulation techniques for Abrasive-based machining processes is shown in Table 5. Based on the classification presented in Figure 9 above, the detailed description of each classification has been elaborated in the next sections.
A summary of the shortcomings of simulation techniques for Abrasive-based machining processes is shown in Table 5.

Computational Fluid Dynamics (CFD)
This approach is applied to simulate the fluid flow problems in any system aided by numerical methods by integrating the problem's equations. The characteristics of shearthinning, Newtonian, non-Newtonian, or the Maxwell fluid are considered while choosing the equations of Navier-stokes to model the flow problem. This method is also applied to take various flow channel geometry into account, locate the optimum parameters for the process, and reduce cost. The simulation performed by Computational fluid dynamics is performed using the available commercial software such as ANSYS ACE, ANSYS CFX, and ANSYS FLUENT using Ansys 2021 R1 software version developed by "Ansys", was founded in 1970 in western Pennsylvania, United States of America. The general working steps that happen inside the software to apply the CFD approach are shown in Figure 11 below.

Computational Fluid Dynamics (CFD)
This approach is applied to simulate the fluid flow problems in any system aided by numerical methods by integrating the problem's equations. The characteristics of shearthinning, Newtonian, non-Newtonian, or the Maxwell fluid are considered while choosing the equations of Navier-stokes to model the flow problem. This method is also applied to take various flow channel geometry into account, locate the optimum parameters for the process, and reduce cost. The simulation performed by Computational fluid dynamics is performed using the available commercial software such as ANSYS ACE, ANSYS CFX, and ANSYS FLUENT using Ansys 2021 R1 software version developed by "Ansys", was founded in 1970 in western Pennsylvania, United States of America. The general working steps that happen inside the software to apply the CFD approach are shown in Figure 11 below. Figure 11. Steps for applying the CFD in software.

Pre-Processing
In this step, the model of CAD that defines the fluid domain in our problem, including any solid boundaries, is prepared first. This step is followed by giving the CFD solver's required initial and boundary conditions to run the simulation. The other conditions, including the physics needed for the problem along with the meshing conditions, are also fed to the commercial code used for the CFD simulation of the problem. This initial step of creating the CAD model can be performed under the same roof within the software

Pre-Processing
Meshing Solution Post-Processing Figure 11. Steps for applying the CFD in software.

Pre-Processing
In this step, the model of CAD that defines the fluid domain in our problem, including any solid boundaries, is prepared first. This step is followed by giving the CFD solver's required initial and boundary conditions to run the simulation. The other conditions, including the physics needed for the problem along with the meshing conditions, are also fed to the commercial code used for the CFD simulation of the problem. This initial step of creating the CAD model can be performed under the same roof within the software used for the flow simulation itself [111].

Meshing
Meshing is a process of discretization where the fluid volume is converted into many more minor elements or volumes. However, many discretization methods are available, and most CFD software uses the Finite Volume Method (FVM), while some exceptional CFD tool uses the Finite Element Method (FEM). The Meshing mechanism involves subdividing the extensive computational domain into many discrete cells called finite volumes (or finite elements if FEM occurs) depending on the CFD toolbox used. The mesh can range from structured or unstructured, uniform or non-uniform, and even with multiple elements like hexahedral or polyhedral in the mesh domain. Meshing remains the main paramount parameter for a CFD simulation to work as desired. Any irregularity in the mesh will crash the simulation when run by the solver [111].
Due to the complexity of Partial differential equations, it is difficult for them to be solved using the computer; therefore, they are converted to another form to be solved in the computer. This process is called discretization, which provides getting continuous equations. By solving a matrix, the algebraic equations can be obtained from the Partial differential equations, and to do that, it is required to identify points. The discretization method is also known to be the "Finite volume method."

Solving
This step is where the actual application of the CFD technique happens. During this step, the CFD software solves the unknowns present in the Navier Stokes equations. The linear system of equations is solved iteratively, coming up with a better value after each iteration. When explicating the integration scheme, the uncoupled equations can result from the linear system, whereas the coupled equations can result from the implicit equation.
With turbulent flow problems, k-epsilon and k-w equations are solved used, and with a laminar flow, the normal Navier-stokes equation is solved.

Post-Processing
The results are obtained in this final step, where the pictorial depiction of the flow parameters is exported and used to study the problem further. Many plots from various sections depicting the distribution or the flow parameters are conducted in post-processing using inbuilt tools or open-source tools like ParaView.
With the advancements coming up in the industrial process, there is an essential role for the process safety in operation and the design of the process facilities. A significant effort was paid to ensuring that "chemical process industries" are at the safest and most secure workplaces among several other industries. However, there is a requirement to analyze the safety of the process mechanism to predict the scenarios. Using experiments to quantify the process safety mechanism is considered to be insufficient as it is highly "resource-intensive." Fortunately, with the increase in the capabilities of the computation and the enhancement in the "Computational Fluid Dynamics (CFD)" tools, the study of hazardous scenarios and the fundamental mechanics have been improved [112] and has become sophisticated [113]. This approach can model the fluid by solving numerical equations besides capturing the effect of viscosity. CFD can also predict chemical reactions, mass and heat transfer, and fluid flow [114,115]. CFD also has several applications in modeling and simulation in the engineering field [116].

Applications
The thermodynamic and hydrodynamic performance of crude oil before heating in refineries is impaired by fouling. There are inappreciable fouling levels at the early stage of fouling which cannot be determined due to the lack of data from the experiments. Therefore, computational fluid dynamics (CFD) is applied for a better understanding. The CFD technology can account for the chemical reaction, the turbulent flow, and the aging in the tubes present in the heat exchanger. The study [117] simulated three-dimensional CFD under different operating conditions to predict the deposition rates and determine the importance of the fouling mechanism. The interaction between the processes of precipitation fouling and the chemical processes can be characterized by the "interference factor" used to evaluate the suppression extent of the chemical mechanism [117].
The firing range depends critically on the ventilation systems for exposure to ammunition refuse personnel. CFD is used to simulate the ventilation system and predict airflow patterns. The purpose of this model is to identify back circulation, which emits fire. The validation of the solution of the model and the input parameters of the model are the values measured at several locations. It has been evaluated using the analysis of the massless particles fired at various heights corresponding to kneeling, standing, and prone firing positions. Using the recommended technical guidance for the ventilation of the firing ranges of the small arms compared to the dynamic flow of the system, different predictions are finished. The results showed that the particles return behind the shooter [118]. The air velocities of lower supply correspond to the particle's successful movement and the uniform flow under investigation in a small firing range. It has been demonstrated that the models of CFD are advantageous as it provides optimization in the air velocity of the supply and enhances the design standards [119]. CFD is also used to model the 3D Concrete Printing to predict the shape of cross-sectional 3D printed segments, which can be achieved through the simulations of "virtual printing." It has been found that the numerical results of CFD are appropriate to the results of the performed experiments [120].
CFD has the edge over the other modeling methods because of the sophistication available at its disposal. The angle of the impact made by the abrasive particles upon hitting the target surface in Ultrasonic Assisted Abrasive Flow Machining (UAAFM) could be visualized from the simulation results, and its outcome is therefore evaluated. With the results obtained from the CFD model of the UAAFM, the effect of the temperature rise on the media stability can be predicted, and major decisions could be undertaken [121]. The abrasive wear on a turbine blade's surface can also be evaluated, helped by CFD simulation. The study helps to identify the effects of the size and shapes of the quartz particles on the abrasive and corrosive attrition of the blade [122]. Kim et al. have investigated the deburring of AL6061 material using the AFM process. The authors have conducted the experimental study and three-dimensional CFD simulation of the deburring process and compared both experimental as well as numerical results [123]. Fu et al. have also investigated the rheological and finishing behavior of AFM media using the experimental study as well as CFD simulation study [124]. Pradhan has compared the experimental and CFD simulation results to predict erosion behavior during hot abrasive jet machining (HAJMing) [125]. Similarly, authors have modeled the HAJMing process using CFD analysis [126]. Amar et al. developed the analytical and numerical models for abrasive water slurry jet machining (AWSJM) to fabricate the micro-channel using experimental study and CFD simulation, respectively [127]. Zou et al. simulated the precision machining of straight internal gear using the large eddy simulation method on the CFD platform [128]. Chen et al. [129] have characterized the motion of abrasive particles using CFD-DEM modeling for the abrasive air jet machining process. The authors have analyzed the effect of nozzle pressure ratio (NPR) on the flow field of abrasive particles and stress distribution on the workpiece for the optimum parameters of abrasive air jet machining. The mesh model of the physical model used for CFD-DEM modeling and a few important simulation study results are shown in  machining of straight internal gear using the large eddy simulation method on the CFD platform [128]. Chen et al. [129] have characterized the motion of abrasive particles using CFD-DEM modeling for the abrasive air jet machining process. The authors have analyzed the effect of nozzle pressure ratio (NPR) on the flow field of abrasive particles and stress distribution on the workpiece for the optimum parameters of abrasive air jet machining. The mesh model of the physical model used for CFD-DEM modeling and a few important simulation study results are shown in Figures 12-14.   [130]. Similarly, the authors have also simulated the flow of abrasive media using the CFD simulation technique by adopting the carreauyasuda and wall slipping models. The rheological behavior of AFM media is affected using a material removal mechanism [131]. The mesh model and a few results of these CFD simulations are shown in Figures 15 and 16 below.
Zhang et al. [132] numerical investigated the machining (AFM) of micro-porous structures using the CFD simulation technique by Brid-Carreau model coupled with mixture and discrete phase models. The authors have demonstrated the simulation results and surface morphology for the finishing of the micro-hole, as shown in .     structures using the CFD simulation technique. The authors have concluded that the flow time (21 s), which was calculated at a velocity of 1.2 mm/s, further demonstrated the polymer chains' remaining stretched state [130]. Similarly, the authors have also simulated the flow of abrasive media using the CFD simulation technique by adopting the carreauyasuda and wall slipping models. The rheological behavior of AFM media is affected using a material removal mechanism [131]. The mesh model and a few results of these CFD simulations are shown in Figures 15 and 16 below.   Zhang et al. have simulated the behavior of visco-elastic AFM media on micro-slit structures using the CFD simulation technique. The authors have concluded that the flow time (21 s), which was calculated at a velocity of 1.2 mm/s, further demonstrated the polymer chains' remaining stretched state [130]. Similarly, the authors have also simulated the flow of abrasive media using the CFD simulation technique by adopting the carreauyasuda and wall slipping models. The rheological behavior of AFM media is affected using a material removal mechanism [131]. The mesh model and a few results of these CFD simulations are shown in Figures 15 and 16 below.   Zhang et al. [132] numerical investigated the machining (AFM) of micro-porous structures using the CFD simulation technique by Brid-Carreau model coupled with mixture and discrete phase models. The authors have demonstrated the simulation results and surface morphology for the finishing of the micro-hole, as shown in Figures 17-19.   Zhang et al. [132] numerical investigated the machining (AFM) of micro-porous structures using the CFD simulation technique by Brid-Carreau model coupled with mixture and discrete phase models. The authors have demonstrated the simulation results and surface morphology for the finishing of the micro-hole, as shown in Figures 17-19.   Zhu et al. [133]. used the double cosine kernel function (SCKF) to develop semi-resolved CFD-DEM for investigating abrasive air jet (AAJ) and studied the particle-scale distribution information in the nozzle of the AAJ machine. Some important results of this simulation study are shown in Figures 20 and 21 below.  Zhu et al. [133]. used the double cosine kernel function (SCKF) to develop semi-resolved CFD-DEM for investigating abrasive air jet (AAJ) and studied the particle-scale distribution information in the nozzle of the AAJ machine. Some important results of this simulation study are shown in Figures 20 and 21 below.  The flow pattern inside the AAJ nozzle and (b) the particle distribution after the particles leave the nozzle exit [133].
The summary of CFD-based modeling and simulation of the loose abrasive-based method is shown in Table 6. [115] -Abrasive flow machining (AFM) "Density, the extrusion pressure, piston velocity, abrasive hardness,

ANSYS FLUENT
The process was a precision finishing operation. Therefore, the removed material is low. This amount can be increased The summary of CFD-based modeling and simulation of the loose abrasive-based method is shown in Table 6.
Even though the CFD simulation model has the edge over other processing models, a lot of advancement has occurred in other simulation models. The finite element methodbased simulation model is one of them. The next section describes the significance of FEM simulation-based processes.  "Density, the extrusion pressure, piston velocity, abrasive hardness, particle hardness and media viscosity, workpiece hardness."

ANSYS FLUENT
The process was a precision finishing operation. Therefore, the removed material is low. This amount can be increased with the increase in the number of cycles. The modeling and the simulation are very important due to the high number of parameters [121] Visco-elastic materials "Ultrasonic assisted abrasive flow machining (UAAFM)" "Fluid pressure, the velocity profile of the fluid, temperature distribution in the working fluid, wall shear, angle of impact, and finishing rate" Commercial simulation tool for CFD The impact angle 'θ' plays a major role in the machining process and improves effectiveness. Wall shear helps predict if the process has better finishing rates.
[134] Turbine blades Abrasive wear Size and shape of quartz particles, leakage flow through clearance gaps, blade life, and efficiency.
When compared with the data from the actual erosion in turbines, the CFD results give more information. The leakage flow through clearance gaps of guide vanes causes erosion in the runner blade inlet.

Finite Element Method (FEM)
The techniques from the Finite Element Method are adopted to predict the surface roughness in a more efficient method of magnetic abrasive finishing (MAF), named "ultrasonic-assisted magnetic abrasive finishing (UAMAF)" [135]. The effect of strain hardening is not considered in different modeling techniques [136]. However, there is a significant role in this effect to play. Also, the effect of heat transfer is neglected in other techniques, and only the mechanical problem is considered. The finite element technique can overcome this problem besides solving the nonlinear problems. Also, it is possible to analyze thermomechanical through the advancement in the software. This approach utilizes the techniques of "numerical practices in metal machining". Replacing the range by "finite elements" forming a mesh is considered to be the basic principle of this approach called discretization. The accuracy of the simulation of FEM does not depend on the digit numbers; however, with the increase in the number of the equation, the accuracy is also increased. For the FEM modeling, the available commercial softwares are more trusted than the manual programs. The steps of applying the FEM approach are represented in Figure 22.
particle hardness and media viscosity, workpiece hardness." with the increase in the number of cycles. The modeling and the simulation are very important due to the high number of parameters [121] Visco-elastic materials "Ultrasonic assisted abrasive flow machining (UAAFM)" "Fluid pressure, the velocity profile of the fluid, temperature distribution in the working fluid, wall shear, angle of impact, and finishing rate" Commercial simulation tool for CFD The impact angle 'θ' plays a major role in the machining process and improves effectiveness. Wall shear helps predict if the process has better finishing rates. [134] Turbine blades Abrasive wear Size and shape of quartz particles, leakage flow through clearance gaps, blade life, and efficiency.
When compared with the data from the actual erosion in turbines, the CFD results give more information. The leakage flow through clearance gaps of guide vanes causes erosion in the runner blade inlet.
Even though the CFD simulation model has the edge over other processing models, a lot of advancement has occurred in other simulation models. The finite element methodbased simulation model is one of them. The next section describes the significance of FEM simulation-based processes.

Finite Element Method (FEM)
The techniques from the Finite Element Method are adopted to predict the surface roughness in a more efficient method of magnetic abrasive finishing (MAF), named "ultrasonic-assisted magnetic abrasive finishing (UAMAF)" [135]. The effect of strain hardening is not considered in different modeling techniques [136]. However, there is a significant role in this effect to play. Also, the effect of heat transfer is neglected in other techniques, and only the mechanical problem is considered. The finite element technique can overcome this problem besides solving the nonlinear problems. Also, it is possible to analyze thermomechanical through the advancement in the software. This approach utilizes the techniques of "numerical practices in metal machining". Replacing the range by "finite elements" forming a mesh is considered to be the basic principle of this approach called discretization. The accuracy of the simulation of FEM does not depend on the digit numbers; however, with the increase in the number of the equation, the accuracy is also increased. For the FEM modeling, the available commercial softwares are more trusted than the manual programs. The steps of applying the FEM approach are represented in Figure  22.

Pre-Processing
Processing Post-Processing The Finite Element Method (FEM) was developed in the 1960s by Clough in his book "The finite element method in plane stress analysis". This approach is used in spatial modeling and has several applications, such as modeling aircraft structures [137,138]. This method can be defined as a computational technique used to obtain approximate solutions in engineering for boundary value problems. This approach allows permitting the continuum to be discretized into finite elements. By assembling the discretized elements as per node properties, it is easy to estimate the continuous domain. Therefore, the FEM has been applied to solve several problems in engineering and applied science [139][140][141].
Several researchers have used FEM analysis for the simulation of the process of solid particle erosion that occurs in abrasive machining [142][143][144][145][146]. In [147], the author has simulated the impact of several spherical particles using the 3D FEM model and estimated erosion rates by calculating the average possible damage as well as extrapolating to level 1.0. However, the obtained results were in good agreement with the experimental results that have already been published in the literature.
In [148], the author used the FEM to simulate the "particle erosion process" by making use of the model of "elastoplastic material" with "failure criteria." The changes in erosion rate with that of the particle velocity and impact angle coincided with the outcome obtained from the analytical and experimental procedures. Residual stress has also been analyzed, which is generated during the period of the occurring erosion process, and it was observed that the profile was similar to the one seen in the "shot peening process." In [149], the author has analyzed the influence of the "abrasive water jet (AWJ)" parameters on "geometric characteristics" of the erosion caters by making use of the explicit FEM and also verified the results by performing the experiments.

Applications
Boundary layer flow is one of the most important mechanics of the applications of fluid and nanofluid. These applications recently are limited to one-term approximations of similarity. The non-similar fluid boundary layer problems are used widely in the industrial field. The finite element method promises to solve ordinary differential equations with a boundary layer that is not similar. The study aims to complement and enhance the numerical heat transfer by using the technique of FEM as it can solve more than one system of boundary layer-derived equations. This approach's obtained result is a robust technique when solving the boundary layer equations [150].
This approach has been commonly used in the test of beam bending to analyze the behavior of the cracking of asphalt mixtures. The internal structure of the Asphalt mixture has a great effect on the propagation of the crack. Besides the computational resources' limitation, the asphalt mixture's complex mesostructure leads to using the 2D simulation of fem in contrast to the 3D simulation. However, 2D modeling reflects the characteristics only of the microstructure of the section of one specimen, which precludes comparing the laboratory test results with the simulation results. To obtain several numbers of images of internal mesostructured, the technology of x-ray CT scanning was used. To construct the 2D fem model, image-based modeling has been used to simulate the tests of beam bending. It has been found that FEM modeling and simulations can predict the crack position [151][152][153].
Recently, adhesive joints have been used more extensively than mechanical joints because of the uniform distributions of stresses and lighter weight compared to mechanical joints. It has been found that adhesive joints are more appropriate with composite materials. To study the mechanical behavior of these joints, several simulations were performed. In these joints between the substrate in geometry and the adhesive, there is a scale difference. Therefore, the analysis using mesh generation is very difficult, and there is a need for the techniques of the manual mesh, which is sometimes not efficient, and the quality of the elements becomes worse. To overcome this problem, fem modeling and simulation are used. This method allows us to produce fine mesh by overlaying the local mesh patch on the existing mesh. Therefore, there is no need for re-meshing [154,155].
The Finite Element Analysis is conducted on the magnetic field to estimate the magnetic force needed for the hollow on the workpiece. The volume of material being removed by a particular indented abrasive is calculated to gauge the Surface Roughness. The study focuses on the relation between the Material Removal Rate (MRR) and the other parameters like the finishing speed, the nature of abrasive material, its size, and the concentration [156]. The finite Element Method is also the answer when it comes to modeling complex magnetic abrasive machining processes. It is being deployed to analyze the Cylindrical Electro-Chemical Magnetic Abrasive Machining (C-EMAM) performed on stainless steel (AISI-420) that is magnetic. C-EMAM helps achieve higher efficiency finishing when working on advanced material cylindrical jobs by abrasion and electrochemical dissolution. FEM aids in the calculation of the magnetic field distribution where the cylindrical workpiece is kept. The simulation results help in comparison with the experimental results of the material removed and the surface roughness and its validation [157].
The Finite Element Analysis also finds its use in the free abrasive polishing in an anhydrous environment to achieve stability during the polishing process. In the FEM Analysis, helped by an Infra-red camera, the effect of the pressure applied and the speed of rotation on the material removal is studied, along with the surface roughness. Other parameters such as heat generated are also evaluated and studied in this article [158]. Since its inception, the aerospace industry has seen many technological advancements, and current research focuses on advanced material alloys and composites. Due to the stringent requirements for the components to be used in the aerospace industry, exceptional accuracy, fatigue strength, and surface integrity are needed. Aerospace materials like Stainless Steel AISI 4130, Ultra high strength steel AF1410, and the AerMet 100 are employed due to their high fracture toughness and corrosion stress resistance. To retain the surface integrity of these substances, abrasive machining is used. The Finite Element simulations of these machining processes are used to predict and compare the results with their existing data from grinding models available [159]. The flow chart of FEM analysis is shown in Figure 23 below.
The Finite Element Analysis is conducted on the magnetic field to estimate the magnetic force needed for the hollow on the workpiece. The volume of material being removed by a particular indented abrasive is calculated to gauge the Surface Roughness. The study focuses on the relation between the Material Removal Rate (MRR) and the other parameters like the finishing speed, the nature of abrasive material, its size, and the concentration [156]. The finite Element Method is also the answer when it comes to modeling complex magnetic abrasive machining processes. It is being deployed to analyze the Cylindrical Electro-Chemical Magnetic Abrasive Machining (C-EMAM) performed on stainless steel (AISI-420) that is magnetic. C-EMAM helps achieve higher efficiency finishing when working on advanced material cylindrical jobs by abrasion and electrochemical dissolution. FEM aids in the calculation of the magnetic field distribution where the cylindrical workpiece is kept. The simulation results help in comparison with the experimental results of the material removed and the surface roughness and its validation [157].
The Finite Element Analysis also finds its use in the free abrasive polishing in an anhydrous environment to achieve stability during the polishing process. In the FEM Analysis, helped by an Infra-red camera, the effect of the pressure applied and the speed of rotation on the material removal is studied, along with the surface roughness. Other parameters such as heat generated are also evaluated and studied in this article [158]. Since its inception, the aerospace industry has seen many technological advancements, and current research focuses on advanced material alloys and composites. Due to the stringent requirements for the components to be used in the aerospace industry, exceptional accuracy, fatigue strength, and surface integrity are needed. Aerospace materials like Stainless Steel AISI 4130, Ultra high strength steel AF1410, and the AerMet 100 are employed due to their high fracture toughness and corrosion stress resistance. To retain the surface integrity of these substances, abrasive machining is used. The Finite Element simulations of these machining processes are used to predict and compare the results with their existing data from grinding models available [159]. The flow chart of FEM analysis is shown in Figure 23 below.  have simulated the impact of high-speed abrasive particles on the tensile block of Al6061-T6 and Ti-6Al-4V materials using a hybrid approach of SPH-FEM. The authors have simulated the results of the material properties changes from the single impact of abrasive particles during erosion [162]. Similarly, Du et al. also have a hybrid simulation technique, i.e., the SPH-DEM-FEM method, to simulate the abrasive particle movement and wear of the nozzle during the AWJ machining process. The results showed that the mixing chamber and focusing tube intersection experienced the most significant wear. [163]. Vasudevan et al. [164] have simulated the AWJ drilling process using a hybrid method of FEA and SPH models. The authors' proposed methodology to simulate the process is shown in Figures 24 and 25 below: Zhang et al. have modeled the process of material removal on thermal barrier coatings (TBCs) using the abrasive water jet (AWJ) by new smoothed particle hydrodynamics (SPH)-based model. The results showed that There is good agreement between the analytical model and the erosion rate as the incidence angle increases [160]. Ozcan et al. have Modeled and simulated the controlled depth of abrasive water jet machining (AWJM) of free-form surfaces using the FEM method by signed distance field (SDF) approach. The authors have validated the simulated results from literature and experimental results [161]. Liu et al. have simulated the impact of high-speed abrasive particles on the tensile block of Al6061-T6 and Ti-6Al-4V materials using a hybrid approach of SPH-FEM. The authors have simulated the results of the material properties changes from the single impact of abrasive particles during erosion [162]. Similarly, Du et al. also have a hybrid simulation technique, i.e., the SPH-DEM-FEM method, to simulate the abrasive particle movement and wear of the nozzle during the AWJ machining process. The results showed that the mixing chamber and focusing tube intersection experienced the most significant wear. [163]. Vasudevan et al. [164] have simulated the AWJ drilling process using a hybrid method of FEA and SPH models. The authors' proposed methodology to simulate the process is shown in Figures 24 and 25 below:  Beigmoradi et al. [165] have investigated the polishing of 2024 aluminum alloy using acoustics energy by a hybrid approach of the FEM-BEM-DEM process. The proposed methodology for acoustic abrasive polishing using the hybrid FEM-BEM-DEM method is shown in Figures 26 and 27 below. Beigmoradi et al. [165] have investigated the polishing of 2024 aluminum alloy using acoustics energy by a hybrid approach of the FEM-BEM-DEM process. The proposed methodology for acoustic abrasive polishing using the hybrid FEM-BEM-DEM method is shown in Figures 26 and 27 below. Beigmoradi et al. [165] have investigated the polishing of 2024 aluminum alloy using acoustics energy by a hybrid approach of the FEM-BEM-DEM process. The proposed methodology for acoustic abrasive polishing using the hybrid FEM-BEM-DEM method is shown in Figures 26 and 27 below.         The summary of FEM-based modeling and simulation of the loose abrasive-based method is shown in Table 7.
FEM is considered to be the numerical method for the continuum materials; however, in DEM, there is no need for such an assumption. Furthermore, the FEM is considered to be the element-based method, while the other one is considered to be the particle-based approach. In this way, DEM is considered the most accurate method that provides results efficiently compared to the FEM because of the compatibility equation, linear response in FEM, and the constitutive relation. The detailed literature analysis of AFM methods using DEM has been discussed in the section below.

Discrete Element Method (DEM)
This method can be defined as a numerical particle size method applied to model the granular material bulk behavior. In this method, a numerical representation for each particle can be identified through specific properties. There are many physical phenomena in abrasive flow machining. By using numerical simulation, a piece of valuable information can be obtained to understand these phenomena. This method was developed in 1977 by Cundell to solve the problems of rock mechanisms.
There are three stages of the DEM mechanism [169].

1.
Pre-treatment: to identify the arrangement of particles, geometry, and field and specify the physical field characteristics.

2.
Dynamic calculation: the calculation of position, velocity, acceleration, and forces of the particles.

3.
Post-treatment: numerical images and cutting force of the unidirectional composites machining during orthogonal cutting. Figure 29 shows the steps involved in DEM. 1. Pre-treatment: to identify the arrangement of particles, geometry, and field and specify the physical field characteristics. 2. Dynamic calculation: the calculation of position, velocity, acceleration, and forces of the particles. 3. Post-treatment: numerical images and cutting force of the unidirectional composites machining during orthogonal cutting. Figure 29 shows the steps involved in DEM. For modeling the microscale of the granular material, DEM can be used. The practical is evaluated as a separate entity through this approach, and the Newtonian motion equations are solved for each practical. The major challenge facing their usage of them is dealing with several objects. However, due to using the graphics processing unit (GPU) and to perform modern computer systems, this approach is used to solve many problems. Various types of spatial-temporal data can be resulted by using the simulation of DEM. This data has information about the particle velocities and positions. The calculation of DEM is performed in two-stage. These stages are the process simulation and the analysis of the data [170][171][172].
For the industrialized gross national product, the cost of abrasive wear is about 1-4%. Therefore, several studies were investigated to classify the various wear. These types are fretting wear, impact wear, fluid erosion, sliding wear, and fretting wear. To estimate the wear in various industries, especially in mining machinery, different types of numerical models were developed, consisting of two parts, an assessment of the wear and a simulation of the particle motion. The discrete element method has proved its efficiency

Pre-treatment
Dynamic Calculation Post-treatment For modeling the microscale of the granular material, DEM can be used. The practical is evaluated as a separate entity through this approach, and the Newtonian motion equations are solved for each practical. The major challenge facing their usage of them is dealing with several objects. However, due to using the graphics processing unit (GPU) and to perform modern computer systems, this approach is used to solve many problems. Various types of spatial-temporal data can be resulted by using the simulation of DEM. This data has information about the particle velocities and positions. The calculation of DEM is performed in two-stage. These stages are the process simulation and the analysis of the data [170][171][172].
For the industrialized gross national product, the cost of abrasive wear is about 1-4%. Therefore, several studies were investigated to classify the various wear. These types are fretting wear, impact wear, fluid erosion, sliding wear, and fretting wear. To estimate the wear in various industries, especially in mining machinery, different types of numerical models were developed, consisting of two parts, an assessment of the wear and a simulation of the particle motion. The discrete element method has proved its efficiency in protecting the wear in different industries [173]. The Discrete Element Model allows directly modeling micro-mechanical parameters, such as the block's behavior and the building geometry. This method also allows simulation of progressive masonry failure, the masonry structure dynamic behavior, and softening. The DEM simulations are performed using 3DEC software [174].
In [175], the authors have worked on the development of DEM simulations as well as coupled SPH for the scrutiny of the abrasive processes regarding the technical fluids in industrial applications. Here, these hold-ups are intermingled with that of the small solid particles that are held responsible for abrasive wear. Furthermore, in various other industrial applications, the suspensions represent viscoelastic behavior. In the presented research, numerical evaluation of the fluids based on the Oldroyd-B fluid has also been executed using the "SimPARTIXR SPH code". Nevertheless, from the simulations of the "oscillator rheometer," it has been observed that such simulation models are only for the crude approximation to the analysis of the intricate polymer fluids. In addition to these, a technique has also been presented enabling to numerically build the whole 3-D geometry of the single (abrasive) grains by making use of the DEM. The DEM-based abrasive grains have been used in coupled SPH simulation processes of AFM. Furthermore, the rheological simulations being performed for the suspension comprising the Newtonian fluid have also presented that executed simulations are in accordance with the experimentally collected data. This has made use of the suspension for the simulation of AFM.
In [176], inherent laws of the friction mechanism have been analyzed while the process of abrading is going on. The research presents that in the flow code modeling, abrading procedure is easy to simplify to the movement of such grains and particles in a "parallelplate shear friction system". For the said research, PFC2D software has been used to build the "particle flow friction system" along with the set of "parallel plate" as well as model parameters as per the equipment of the abrading procedure as well as the processing materials, and also through control simulation of single variation. Furthermore, a comparison has been made of the output data in order to approximate the impact of variation in parameters on the force chain. From the simulations, it has been analyzed that shear dilatancy is possible to be segmented into three stages: "plastic strain", "macroscopic failure", and "granular recombination". Here, the load rates and the distribution of weak force chains rely on velocity, load, friction coefficient among the granular diameter, granules, and the number of layers of the granular. It has been observed that the load and granular layers upsurge the cause of the "direction of force chain" to be oriented with that of the vertical direction, and thus makes the force chain move in the horizontal direction with the increase in the velocity. With an increase in load, there occurs no variation in shear dilatancy stage, friction coefficient, and velocity. The flow chart of the DEM simulation-based model is presented in Figure 30 below. weak force chains rely on velocity, load, friction coefficient among the granular diameter, granules, and the number of layers of the granular. It has been observed that the load and granular layers upsurge the cause of the "direction of force chain" to be oriented with that of the vertical direction, and thus makes the force chain move in the horizontal direction with the increase in the velocity. With an increase in load, there occurs no variation in shear dilatancy stage, friction coefficient, and velocity. The flow chart of the DEM simulation-based model is presented in Figure 30 below.

Applications
The propagation modeling of cross-contamination of bacteria within bulk material is very important to validate the safety of the microbial of the tree nuts complex post-harvest processing. It has been found that the discrete element method is very efficient in modeling and simulating the propagation of complex bacteria during the material processing of bulk food [177]. Recently, the DEM approach became a promising alternative to the

Applications
The propagation modeling of cross-contamination of bacteria within bulk material is very important to validate the safety of the microbial of the tree nuts complex postharvest processing. It has been found that the discrete element method is very efficient in modeling and simulating the propagation of complex bacteria during the material processing of bulk food [177]. Recently, the DEM approach became a promising alternative to the measurements of the physical properties of the systems of the thermal particulate. DEM has been widely used in several industrial applications. However, the simulation of DEM has not been used in a theory review that underpins the thermal simulation of DEM. Therefore, the study presents a comprehensive review of all thermal model majors and the heat transfer mechanisms using the DEM simulations [178].
Permeability of the bed is considered an important parameter in green iron ore pellet's induration in the grate furnaces. For achieving higher quality, lower fuel consumption, and higher productivity, optimization and control must be pursued. The discrete element method is used to improve the permeability of the bed to control the green pellets distribution size to feed the furnace of the induration [56]. Handling the wet particles is very difficult in chemical engineering because of the cohesive acting force in these particles, which causes the aggregation in these particles, unlike the dry particles. Therefore, it has been found that the discrete element method is very efficient in the investigation of the wet particles' behavior [179]. The Discrete Element Method (DEM) is used alongside Finite Element Analysis or alone in drilling and abrasive machining processes. Initially, the interaction between the cutter and the surface material was analytically modeled. But, with the advancement in DEM processes, they are used to model and explore the cutting action of the PDC drill bit and its drilling mechanism [180].
Another application of the Discrete Element Method (DEM) is when fused Silica Glass undergoes Abrasive grinding. The subsurface (SSD) damage is measured aided by the wet etch technique. This article helps to recognize the impact of process parameters on the SSD depth. The DEM overthrows the other simulation and modeling methods due to its ability to calculate a larger number of interfaces and contacts. The results obtained from the Discrete Element method based on its microscopic model go well with the measurements derived from the macroscopic model, and hence optimization of the SSC depths is now made possible [181].
Beigmoradi et al. studied the abrasive drag finishing using the DEM method [182]. Similarly, Ge et al. used the DEM method to simulate the cavitation effect of traditional soft abrasive flow (SAF) polishing for large workpieces [183]. Huang et al. investigated the effect of wear on the sand blasting process by hybrid CFD-DEM coupling method. The results showed that the increases in particle diameter, density, shear modulus, stacking height, and inlet pressure are associated with an increase in the sandblasting machine's erosion volume [184]. Zhu et al. [185] have developed the computational model by hybridization of the Monte Carlo (MC) method and discrete element method (DEM), to study the erosion of glass during micro-hole fabrication using the abrasive air jet (AAJ) process. Results show that the particle interference effects, which are crucial for regulating material removal and microstructure shape in AAJ machining, may be captured by the provided DEM model, and the proposed methodology and computation are shown in Figures 31 and 32 below.
The model setup and computational results of the above MC-DEM simulation are shown in Figures 33 and 34 below.
The summary of DEM-based modeling and simulation of the loose abrasive-based method is shown in Table 8.  Figure 31. The flowchart of the modeling framework of the MC method coupled with the DEM for multi-particle impact erosion [185]. Figure 31. The flowchart of the modeling framework of the MC method coupled with the DEM for multi-particle impact erosion [185]. Figure 31. The flowchart of the modeling framework of the MC method coupled with the DEM for multi-particle impact erosion [185].   Model setup for the abrasive particle jet against the solid specimen in the DEM simulation [185]. Figure 33. Model setup for the abrasive particle jet against the solid specimen in the DEM simulation [185]. DEM methods are considered to be relatively intensive in computational terms, which limits and restricts either the particle numbers or the run length of the simulation. Therefore, in order to tackle the issue of a longer time requirement, another simulation named Molecular Dynamic simulation has been introduced. DEM is closely in accordance with the MDS. However, the method is usually distinguished through the rotational DOF and state contact and is usually considered useful for complicated geometries. The next section describes the detailed analysis and literature on Molecular dynamic simulation.

Molecular Dynamic Simulation (MDS)
This method can be defined as a computational tool that is used to refer to how the orientations, velocity, and position of molecules are developed. The force and energy are related to the configuration by the simulation depending on the interaction of molecular scale. Using newton's second law, the acceleration of the particle can be calculated. The "molecular dynamic simulation" gives an advanced "spatial and temporal resolution of the cutting process" than the "continuum mechanics approach". The atom position can be obtained with "high time resolution (smaller than the atom vibration motion", which is of the order 1 × 10 −15 s or 1 fs) by solving the "Newton equation of motion". The molecular dynamic simulation is computationally intense. "Molecular statics analysis" follows the point where only the resultant force (Fr) is equal to zero. Here, the atom follows the "minimum potential energy position". It is thus a "quasi-static method". The general steps of the MDS technique are shown in Figure 35.
DEM methods are considered to be relatively intensive in computational terms, which limits and restricts either the particle numbers or the run length of the simulation. Therefore, in order to tackle the issue of a longer time requirement, another simulation named Molecular Dynamic simulation has been introduced. DEM is closely in accordance with the MDS. However, the method is usually distinguished through the rotational DOF and state contact and is usually considered useful for complicated geometries. The next section describes the detailed analysis and literature on Molecular dynamic simulation.

Molecular Dynamic Simulation (MDS)
This method can be defined as a computational tool that is used to refer to how the orientations, velocity, and position of molecules are developed. The force and energy are related to the configuration by the simulation depending on the interaction of molecular scale. Using newton's second law, the acceleration of the particle can be calculated. The "molecular dynamic simulation" gives an advanced "spatial and temporal resolution of the cutting process" than the "continuum mechanics approach". The atom position can be obtained with "high time resolution (smaller than the atom vibration motion", which is of the order 1 × 10 −15 s or 1 fs) by solving the "Newton equation of motion". The molecular dynamic simulation is computationally intense. "Molecular statics analysis" follows the point where only the resultant force (Fr) is equal to zero. Here, the atom follows the "minimum potential energy position". It is thus a "quasi-static method". The general steps of the MDS technique are shown in Figure 35. This method is used to see the transformation of the crystal talking in the workpiece and tool during machining and provide knowing the science behind the process. The software used for the simulation is "large atomic and molecular massively parallel simulator (lamps)" machining, and for the visualizing and results, an available "visualization tool (Ovito)" is used. Ovito is a 3D "visualization software" that has been established for "postprocessing" of the data acquired from the "molecular dynamics" or "Monte-Carlo This method is used to see the transformation of the crystal talking in the workpiece and tool during machining and provide knowing the science behind the process. The software used for the simulation is "large atomic and molecular massively parallel simulator (lamps)" machining, and for the visualizing and results, an available "visualization tool (Ovito)" is used. Ovito is a 3D "visualization software" that has been established for "post-processing" of the data acquired from the "molecular dynamics" or "Monte-Carlo simulation." This software uses "object-oriented C++" and can be regulated with the help of "python scripts". It can be easily extended utilizing a "plug-in" interface. This approach is defined as a theoretical method that provides spatial resolution besides high temporal. This simulation technique is an important method used in the atomic scale process due to its high efficiency in this field [187]. This approach is used to model and simulate scratching the polycrystalline and monocrystalline silicon carbide with diamond grit to predict the behavior of the material removal [188][189][190]. This method is also used in simulating the process of copper milling as it has proved its efficiency in investigating the deformation in copper machining [191].
"Vibration Assisted Nano Impact-machining by Loose Abrasives (VANILA)" is a nanomachining process that includes the principles of tip-based nanomachining and the machining of vibration-assisted abrasives. MDS is used in various fields to investigate the parameter of this process, such as the particle size and impact velocity [192,193]. With the further advancements in computational mechanics, Molecular Dynamics Simulations (MDS) are recognized and accepted in nanoengineering and scientific fields to study the changes and movements that happen at such small scales. The MDS is used to evaluate the monocrystalline gallium nitride (GaN) that undergoes abrasive machining at a nanoscale. This gives an advantage over the other modeling techniques as it gives insights into the monocrystalline material's deformation behavior. The MDS enriches the nano enthusiasts to have insights into the effects due to changes in temperature and stress on the velocity of the cutting, Depth of Cut (DOC), or the usage of cub-octahedral abrasive for the machining process. The nano-scale deformation behaviour of the GaN material when such high-precision machining processes are studied [194].
The FEM applied in the abrasive polishing of materials is Molecular Dynamic Simulation (MDS), which finds its way useful for further studies. Chemical Mechanical Polishing (CMP) is one of the noteworthy methods available to achieve the highest precision required to polish diamonds. But, the internal frictional behavior during the machining process needed some explanation, and MDS sheds some light on the frictional behaviour during different pressure conditions. Molecular Dynamics supports the theory behind the damage reduction when machining the diamond and its precision [195].
Asadollahi et al. investigated the abrasive water jet machining of Al 2 O 3 , Graphite, and Silicon Carbide reinforced on aluminum alloy 7075 metal matrix (MMC) using an experimental method and MDS simulation. According to the experimental findings, the kerf taper angle decreased as the flow rate increased (from 305 to 470 g min1), traverse speed decreased (from 167.64 to 55.08 mm min1), weight percentage increased (7.5 to 2.5), water pressure decreased (270 to 200 bar), and SiC was used as a matrix reinforcement particle. This indicates that the machining was of high quality [196]. Meng et al. have also investigated the abrasive-based polishing of single crystal silicon using MDS simulation [197]. Zhou et al. investigated the mechanism of ultrasonic-assisted nano-cutting of sapphire using the MDS simulation technique. According to the results of the MDS simulation, cutting with tangential and radial vibrations has lower tangential, normal, and overall forces than cutting with conventional methods [198]. Ranjan et al. [199] have investigated the mechanical polishing of polycrystalline Cu using MDS simulation. Figures 36 and 37 below illustrate the material modeling and the impact of abrasive size and cutting velocity on material re-removal in nano-cutting.  Su et al. [200] simulated the cutting characteristics of monocrystalline silicon in the elliptical vibration nano cutting (EVNC) process using MDS simulation. The authors examined the amplitude ratio mechanism's effect on hydrostatic pressure's thermal properties of the elliptical vibration nano-cutting force. The MDS model and comparative analysis of cutting forces during EVNC are shown in Figures 38 and 39 below.
Liu et al. [201] simulated the ductile deformation for nano cutting of Silicon Carbide using MDS simulation. The MDS model of nanometric cutting on monocrystalline 3C-SiC and displacement vectors of the workpiece atoms during nanometric cutting is shown in Figures 40 and 41 below.  Figure 36. (a) Lattice optimization of the copper lattice with respect to cohesive energy, (b) dimensions of the atomistic model of p-Cu work material (a rectangular block of size 150 Å × 80 Å × 80 Å), and abrasive particle (spheres of sizes Ø43 Å and Ø57 Å) in perspective view, (c) front sectional view of the atomistic model of grains and their boundaries with initial height of abrasive and depth of nanoindentation and (d) equilibration graphs for cohesive energy and temperature for initial 25 ps [199].            Similarly, Liu et al. [202] simulated the multi-angle precision micro-cutting of a single-crystal copper surface using the MDS simulation technique. Numerical analysis model of boron nitride particles cutting and their respective cutting forces in different directions under various micro cutting angles are shown in Figures 42 and 43 below.  The summary of MDS-based modeling and simulation of the loose abrasive-based method is shown in Table 9.  The summary of MDS-based modeling and simulation of the loose abrasive-based method is shown in Table 9. From the above literature, it has been observed that the MDS method can possibly be distinguished through the rotational DOF and state contact and can be considered to be useful for complicated geometries. However, it does not assist in understanding the relationships among different variables. For this purpose, a Multivariate regression model has been presented that tells the relationship or a link among the variables within the data. It assists in comprehending the correlation among independent as well as dependent variables. The next section tells a detailed study of the literature review regarding the multi-variable regression analysis.

Multi-Variable Regression Analysis
Several researchers in "loose abrasive-based machining" use experimental data. The researchers developed a mathematical model using various techniques called "multivariable regression analysis (MVRA)". MVRA is considered the "statistical modeling process" used to estimate the association between the output machining parameters and input machining parameters. This model is stated through the following equation [204].
where y represents the dependent output and x 1 , x 2 , . . . . x k represents the input. The nonlinear relationship between variables that may be dependent or independent can be modeled using this approach. Regression analysis is considered a necessary means of this approach is regression analysis [205]. The discrete data can be approximated by curve fitting of the analytic expression, which obtains a mathematical correlation using the given discrete point. By evaluating the tiny linear segments, these interpolation points can be connected to curves. The smooth curve can be formed when selecting the interpolation points properly. The curve straightening is a mean of the curve fitting. The linear equations can be obtained from the nonlinear functions by simple variable transformations such as the S-curve function, power function, and exponential function. There are quadratic and cubic fitting methods for quadratic functions and cubic functions [206][207][208].

Applications
Within machine learning and the statistical literature, the techniques of regularization are used to construct sparse models. Most strategies of the regularization only work for the data when treating all the predictors identically. However, several problems include various types of predictors and need a term of tailored regularization. It has been found this approach is efficient in overcoming this problem [209]. The surface irregularities are one of the major players that determine the quality of processed engineering materials. The surface roughness of the Magnesium alloy (AZ 91) is modeled helped by regression analysis and Fuzzy logic as the alloy undergoes an Abrasive Flow Machining (AFM) process [210]. The pressure of water, speed of traverse, and the stand-off distance are noted as the input variables for Regression analysis. Using the regression equation, the Surface Roughness values are projected, and it is compared with the values obtained from the Fuzzy Logic prediction. Investigation resulted in the fact that the maximization of pressure and minimization of stand-off distance caused the surface roughness to reduce during the Abrasive jet machining process [128].
The precision machining of ceramic materials cannot be left unnoticed due to its application in various fields, including biomedical and aerospace. Despite the difficulty faced during the machining of ceramic materials, unconventional machining, such as Hot Abrasive Jet machining, is employed to machine advanced ceramic materials like Silicon Nitride (Si 3 N 4 ). The difficulty in the machining process requires an optimization of the process not only to ease them but also to increase safety. The process parameters like air pressure, the distance of nozzle tip, and the abrasive temperature are analyzed using variance, and Regression models are obtained. The Regression analysis and CFD analysis are carried out to discover the optimum machining parameters for HAJM [126]. The summary of MVRA-based modeling and simulation of the loose abrasive-based method is shown in Table 10.
Silicon Nitride (Si 3 N 4 ) Even though MVRA is a process that speeds up everything quickly considering the mathematical model, there is still a need to develop the model and technique that can process the abrasive mechanism in a more advanced way. For this purpose, the ANN technique has been developed. The details of this technique in terms of AFM have been described in detail in the next section.

Artificial Neural Network (ANN)
The Artificial Neural Network (ANN) can be defined as an artificial intelligence-based technique like the biological neural network of the human brain used in problem-solving. Several researchers use this technique to develop a mathematical model to analyze the parameters of loose abrasive-based machining. This technique works as a parallel massively distributed processor with a strong neural network for performance calculations for the mathematical model or specified task. This method is developed to provide computational techniques that can store and acquire information using trained algorithms such as unsupervised learning, reinforcement learning, and supervised learning. For loose abrasive-based machining, the modeling of the backpropagation neural network is used [212,213]. It is called a multi-layered network that has some hidden layers. The structure of this approach is shown in Figure 44. Even though MVRA is a process that speeds up everything quickly considering the mathematical model, there is still a need to develop the model and technique that can process the abrasive mechanism in a more advanced way. For this purpose, the ANN technique has been developed. The details of this technique in terms of AFM have been described in detail in the next section.

Artificial Neural Network (ANN)
The Artificial Neural Network (ANN) can be defined as an artificial intelligencebased technique like the biological neural network of the human brain used in problemsolving. Several researchers use this technique to develop a mathematical model to analyze the parameters of loose abrasive-based machining. This technique works as a parallel massively distributed processor with a strong neural network for performance calculations for the mathematical model or specified task. This method is developed to provide computational techniques that can store and acquire information using trained algorithms such as unsupervised learning, reinforcement learning, and supervised learning. For loose abrasive-based machining, the modeling of the backpropagation neural network is used [212,213]. It is called a multi-layered network that has some hidden layers. The structure of this approach is shown in Figure 44. This approach is based on the mathematical model construction that considers various kinds of uncertainties [214,215]. Artificial Neural Networks (ANN) allow modeling of cutting in an Artificial Neural Networks [216], which applies the known empirical dependencies and accumulated research experience. This approach can be used in creating adaptive models to control and direct the cutting process, predict the parameter of the process, and select the optimal mode for cutting [217]. The ANN approach seems efficient by considering the processors of computations for data compression, various operations, data integration using multiple sensors, noise filtering, and modeling and forecasting. There has already been experience gained in using the ANN method in grinding processes to predict the wear of the abrasive tool, the appearance of cauterization, cutting force, and This approach is based on the mathematical model construction that considers various kinds of uncertainties [214,215]. Artificial Neural Networks (ANN) allow modeling of cutting in an Artificial Neural Networks [216], which applies the known empirical dependencies and accumulated research experience. This approach can be used in creating adaptive models to control and direct the cutting process, predict the parameter of the process, and select the optimal mode for cutting [217]. The ANN approach seems efficient by considering the processors of computations for data compression, various operations, data integration using multiple sensors, noise filtering, and modeling and forecasting. There has already been experience gained in using the ANN method in grinding processes to predict the wear of the abrasive tool, the appearance of cauterization, cutting force, and roughness [218,219]. This approach is an effective method to deal with the characteristics of the data-driven. Between the output and the input, ANN is considered a universal approximation. The nonlinear behaviour of the input can be approximated using ANN after determining the parameter [220].
The memory requirement of the ANN method is dependent only on the scale of the dataset, despite the dimension and structure. The characteristics embodied can be represented by the ANN in real-time simulation on a large scale of the dataset. One advantage of ANN is that the parallel structure of ANN makes it much more compatible with the FPGA, which can provide reducing the overall computing latency [221][222][223].

Applications
This approach is used in predicting the components and the temperature of cutting force when using the minimum quantity lubrication technology to grind the nickel. The simulation results prove that the ANN is very efficient in predicting the optimal output [218]. Understanding the reservoir parameters optimizes the designated reservoir usage. There are four main parameters: injection rate, injection temperature, fracture permeability, and well spacing. It has been found that the ANN technology can be used in analyzing the effect of the different factors on the granodiorite reservoir thermal performance. The ANN allowed the enhancement of an efficient system in lesser time. It has been proved that applying the ANN has great accuracy in predicting the production temperature [224].
Sensible Heat (H) and Latent Heat (LE) are two major components of the balance of the energy on the earth's surface, which have a great role in global warming besides water cycling. Several methods are used for the measurement or estimation of these quantities, but they have several limitations. It has been found that the technology of ANN can investigate the relationship between Sensible Heat (H) and Latent Heat (LE) besides inspecting the empirical data [225]. Besides the above-said methods to model the Abrasive Machining processes, Artificial Intelligence in Neural Networks (ANN) makes its contribution [226]. Ti-alloys like Ti-6Al-4V are widely used in the aerospace and medical industries because of their rare quality of low density but higher strength and resistance to corrosion. Due to the lack of sophistication in Conventional machining, Abrasive Waterjet Machining (AWJM) is performed to machine the alloy with lesser heated zones. To come up with the right process parameters to attain the best quality, the parameters are used to develop mathematical equations aided by the Regression Investigation Method (RIM) and Artificial Neural Network (ANN). The models made from the input and output data help us to explore deeper into the control methods and optimization possibilities [227].
With the vision to improve Abrasive Waterjet Machining (AWJM) performance, Taguchi's Design of Experiments method is utilized to choose the input process parameters. The vast set of experimental data obtained is fed into the Artificial Neural Network (ANN) model, and the MATLAB programming language is used to train the data. The ANN resulted in forecasting the data to predict the effect of the input parameters on the system's stability [228].
Deshpande et al. optimized the surface quality and kerf taper angle of AWJ machining of AISI 1018 steel using ANN-based modeling. Surface roughness (Ra) = 2.46 m, taper angle () = 1.25 • , flow rate (Af) = 450 g/min, stand-off distance (Sd) = 3 mm, and traverse speed (Tv) = 85 mm/min were determined to be the optimal response solution [229]. Saini et al. have also modeled the response parameters of abrasive jet machining of composite using ANN-based modeling [230]. Similarly, Chawla et al. have investigated the magnetic abrasive flow machining of MMCs of hybrid reinforced Al/SiC/B4C using ANN-based modeling [48].
Salinas et al. investigated the drag finishing abrasive effect for cutting edge preparation in broaching tools using ANN modeling. The results showed that the reproducibility accuracy of re was found in low training and prediction errors of 1.22% and 0.77% [231]. González et al. have investigated the super abrasive machining of integral rotary parts by grinding flank tools. Authors have modeled the cutting forces using ANN modeling and experimental study [232]. Similarly, authors have compared the flank super abrasive machining vs. flank milling on inconel ® 718 surfaces [233].
The summary of ANN-based modeling and simulation of the loose abrasive-based method is shown in Table 11.

MATLAB
Feeding the data to the ANN is followed by a back propagation algorithm to train the model.
The literature has analyzed that ANN is useful in many ways; however, it does not consider the relationships among several variables. For this purpose, Response surface analysis has been developed, and the details of this analysis model are presented in the next section.

Response Surface Analysis (RSA)
Response surface analysis is also known to be the "statistical technique" for modeling numerical practices. "Response surface analysis (RSA)" is an empirical modeling method used for the "design of experiments (DoE)". The output components of the AFM process for optimization rely on various independent input variables of AFM techniques after the careful DOE. A DOE has tests, known as runs, where transformations are made during the "input process parameters" in order to find out the impact of these variations and transformations in the "output process parameters". Applying RSA is thought to be beneficial for minimizing the expenses of costly analysis approaches (e.g., "finite element method or CFD analysis") and also their "associated numerical noise". This approach was developed by Box and Wilson in the early 1950s, which can be defined as a collection of statistical and mathematical techniques used in the optimization and approximation of stochastic models [141,235,236]. With such a model, the objective function is subjected to random noise, and it is called a stochastic or noisy objective function. This approach has several applications in the development, formulation, and design of new products besides enhancing existing design [237]. The RSA method is also useful as a tool to analyze probabilistically. This approach can be defined as a method used to approximate an unknown function with few computed values [238]. It is used in different processes in experimental optimization as it optimizes the mixture design for consideration of various attributes like workability, durability, cost, and environmental impact [239]. This method is commonly used to analyze slope reliability because of its high efficiency [240].

Applications
In capturing CO 2 , several studies were investigated to maximize the capacity of capture to reach the amine temperature to save energy and reduce cost. The RSM is a tool for the optimization of this process. It has been proved that because of its efficiency, the simulation results showed that the absorber's energy consumption is reduced by 16% in optimum conditions [241]. Researchers study the different parameters affecting the machining process, making improvements for better machining a priority for the Abrasive Waterjet Machining (AWJM). The Abrasive Waterjet Machining of AA6082 is evaluated to improve its efficiency to attain minimal surface roughness with the maximum material being removed. The Response Surface Methodology is employed to optimize the machining process by considering the parameters, including the abrasive feed, stand-off distance, and the nozzle speed [242]. The summary of RSA-based modeling and simulation of the loose abrasive-based method is shown in Table 12.

Stochastic Modeling and Simulation (Data Dependent System)
There is a relation between the word "Stochastic "and a random variable. This approach can be defined as a mathematical tool developed to estimate the probability distributions of the potential outcomes by allowing random variation in the input. The random variation depends on the fluctuations seen in historical data by making use of standard time-series techniques for a selected period. To research the "surface finish" after applying AFM, Williams et al. utilized the modeling and simulation of stochastic approaches named the "data-dependent system (DDS)" [79,243]. It has been concluded that the media topography should be understood to control the procedure and enhance productivity and efficiency with the least cost. Jain et al. have established the "stochastic mathematical model" for simulation purposes in order to forecast the density of the "active grain on the media surface," making a connection with the surface of the work [244]. The summary of DDS-based modeling and simulation of the loose abrasive-based method is shown in Table 13.
From the literature presented above, various models and simulation techniques have been analyzed and discussed in detail. Further, in the next sections, the mathematical modeling of these techniques will be elaborated on and discussed, taking into consideration the abrasive-based machining process. AFM is defined as the process of non-traditional-finishing, which is used in surface finishing, generating compressive residual stresses, and edge contouring [168,208,245]. Some factors affect the finishing of the surface besides the material removal in the abrasive flow machining process, such as extrusion pressure, piston velocity, media viscosity, and particle size. It is not easy to perform experiments to analyze each factor's effect in a short time. To find the optimum parameter, the simulations of CFD were created [7,26,84].
The definition of abrasion in AFM is that the solid material gets removed from the material's surface by sliding abrasive particles from another material [166]. Penetration on the surface of the workpiece will be caused by force acting on the spherical grain. When the grain goes in the horizontal direction, it removes the material. The stock removal amount is the same as the produced grooves' volume on the workpiece. The volume of material removed (MRR), surface roughness and stock removal can be obtained if the depth of the groove and the abrasive grains numbers are known. The roughness profile (assumed as a triangular shape) and the abrasive grain (assumed as a spherical shape) are presented in Figure 45 [28]. The 2D model simulation is performed using the CFD tool using a multiphase m ture model, and the validation of the CFD tool is presented in [28]. The validation of t present CFD simulation with respect to their boundary conditions is shown in Table 14   Table 14. The validation of the CFD tool [28].  The 2D model simulation is performed using the CFD tool using a multiphase mixture model, and the validation of the CFD tool is presented in [247]. The validation of the present CFD simulation with respect to their boundary conditions is shown in Table 14. The simulation results show the effects of the abrasive particles' volume fraction on parameters such as strain rate, horizontal velocity, and dynamic pressure, which are represented in Figures 46-48 [247]. Figure 46 interpreted the simulation results, and it showed that the change of strain rate is inversely proportional to the distance to the surface. It means that the strain rate decreases with the increased distance of abrasive grain from the surface. Strain rate is directly proportional to the material removal rate (MRR), so MRR increases with reducing the distance to the surface.  Figure 47 interpreted the simulation results, and it showed that the velocity of abrasive grains is directly proportional to the radial distance to the surface. It means that the velocity of abrasive particles increases with the increased radial distance of abrasive grain from the surface. The velocity of abrasive grain is directly proportional to the material removal rate (MRR), so MRR increases with increasing the abrasive radial distance to the surface. Radial distance to the surface depends on the viscosity of AFM media. Figure 46. Impact of volume fraction in strain rate [247]. . Impact of volume fraction in horizontal velocity [247]. Figure 48 interpreted the simulation results, and it showed that the dynamic pressure of AFM media is directly proportional to the radial distance to the surface. It means that when the dynamic pressure of AFM media increases, the radial distance of abrasive grain from the surface increases. The dynamic pressure of AFM media is directly proportional to the material removal rate (MRR), so MRR increases with increasing the dynamic pressure of AFM media, and radial distance to the surface depends on the viscosity of AFM Figure 47. Impact of volume fraction in horizontal velocity [247]. Figure 47. Impact of volume fraction in horizontal velocity [247]. Figure 48 interpreted the simulation results, and it showed that the dynamic pressure of AFM media is directly proportional to the radial distance to the surface. It means that when the dynamic pressure of AFM media increases, the radial distance of abrasive grain from the surface increases. The dynamic pressure of AFM media is directly proportional to the material removal rate (MRR), so MRR increases with increasing the dynamic pressure of AFM media, and radial distance to the surface depends on the viscosity of AFM media. Figure 48. Impact of volume fraction in dynamic pressure [247].

Boundary Conditions Applied Parameters
As shown in the graphs, it has been found that when the fraction volume increases, the dynamic pressure is increased towards the wall from the center, and the strain rate increases with the increase in the volume fraction, which increases the shear stress and the deformation per unit time in the near-wall region. For the velocity parameter, the maximum velocity is obtained at the center of the flow domain. The magnitude of the velocity is also increased with the increase in volume fraction [247].
A study was proposed to investigate the different machining parameters and their effect on a novel nanomachining process by Loose Abrasives (VANILA). In this process, an "atomic force microscope (AFM)" is made to be the platform. Between the silicon workpiece and the vibrating AFM probe, the Nano abrasives slurry is injected. The acceleration Figure 48. Impact of volume fraction in dynamic pressure [247]. Figure 47 interpreted the simulation results, and it showed that the velocity of abrasive grains is directly proportional to the radial distance to the surface. It means that the velocity of abrasive particles increases with the increased radial distance of abrasive grain from the surface. The velocity of abrasive grain is directly proportional to the material removal rate (MRR), so MRR increases with increasing the abrasive radial distance to the surface. Radial distance to the surface depends on the viscosity of AFM media. Figure 48 interpreted the simulation results, and it showed that the dynamic pressure of AFM media is directly proportional to the radial distance to the surface. It means that when the dynamic pressure of AFM media increases, the radial distance of abrasive grain from the surface increases. The dynamic pressure of AFM media is directly proportional to the material removal rate (MRR), so MRR increases with increasing the dynamic pressure of AFM media, and radial distance to the surface depends on the viscosity of AFM media.
As shown in the graphs, it has been found that when the fraction volume increases, the dynamic pressure is increased towards the wall from the center, and the strain rate increases with the increase in the volume fraction, which increases the shear stress and the deformation per unit time in the near-wall region. For the velocity parameter, the maximum velocity is obtained at the center of the flow domain. The magnitude of the velocity is also increased with the increase in volume fraction [247].
A study was proposed to investigate the different machining parameters and their effect on a novel nanomachining process by Loose Abrasives (VANILA). In this process, an "atomic force microscope (AFM)" is made to be the platform. Between the silicon workpiece and the vibrating AFM probe, the Nano abrasives slurry is injected. The acceleration of the Nano abrasives is performed by generating Nano cavities on the workpiece surface, hitting the workpiece, and the vibration of the AFM. As loose abrasives, diamond particles are used. To create a numerical model, the (FEM) is used. The validation of the FEM model is obtained by comparing the analytical model results with the simulation results. The FEM simulation, the spherical diamond abrasive, the silicon workpiece, and the boundary conditions for the symmetric surface are represented in Figure 49 [248][249][250].
The results obtained from the simulation showed that for different fraction coefficients, 0.05, 0.3, and 0.5 with several hits (0-25). The operating temperature, besides impact speed, is constant. However, with the increase in the number of hits, the damage volume also increases with the increase in the friction coefficients. It happened because of the different stress triaxiality ratios resulting from the different coefficients. These results are presented below in Figure 50 [248].
It has been found there is a substantial influence of the operating temperature on the volume of the material removal. For a temperature of 20 • C and 100 • C, the critical value 1 could not be reached by the damage variable. However, when the temperature increases to 200 • C, after 16 hits, the damaged volume occurs. When the temperature is increased to 800 • C, the removal of undesired material occurs [141]. piece and the vibrating AFM probe, the Nano abrasives slurry is injected. The acceleration of the Nano abrasives is performed by generating Nano cavities on the workpiece surface, hitting the workpiece, and the vibration of the AFM. As loose abrasives, diamond particles are used. To create a numerical model, the (FEM) is used. The validation of the FEM model is obtained by comparing the analytical model results with the simulation results. The FEM simulation, the spherical diamond abrasive, the silicon workpiece, and the boundary conditions for the symmetric surface are represented in Figure 49 [247][248][249]. The results obtained from the simulation showed that for different fraction coefficients, 0.05, 0.3, and 0.5 with several hits (0-25). The operating temperature, besides impact speed, is constant. However, with the increase in the number of hits, the damage volume also increases with the increase in the friction coefficients. It happened because of the different stress triaxiality ratios resulting from the different coefficients. These results are presented below in Figure 50 [247]. It has been found there is a substantial influence of the operating temperature on the volume of the material removal. For a temperature of 20 °C and 100 °C, the critical value 1 could not be reached by the damage variable. However, when the temperature increases to 200 °C, after 16 hits, the damaged volume occurs. When the temperature is increased to 800 °C, the removal of undesired material occurs [141].
AFM can be implemented to the manufactured additive components with high roughness [236]. However, the obtained geometry may not meet the dimension required because of the inhomogeneous material removal (MR). A solution has been developed by the proposed study [250] to predict the material removal distribution of a component. During the design, the materials are compensated as they are added to the component. This method's advantage and the feasibility were shown on a laser sintered test similar to the nozzle guide vane (NGV). The CFD model was used in this study to predict compensation for the material removal [6,[250][251][252][253][254].
The materials of the rock are widely used in several industrial applications. The mechanical and physical properties of natural rock differ from one another. To explore the removal mechanism besides rock breakage during the machining process, several experiments were performed, including linear cutting, disc cutting, and drilling. The simulation and modeling technologies have proved their efficiencies in exploring the rock material fracture mechanism except for the experiments. These technologies include the discrete element method (DEM), finite element method (FEM), and hybrid finite-discrete element method (FDEM). The FDEM principles are shown in Figure 51. First, the object was divided into solid elements, and then between the adjacent solid elements, the cohesive elements with zero thickness were inserted. The solid elements determine the combined model deformation, and the cohesive elements control the failure process. Traction separation law (TSL) describes the relationship between the relative displacement and the AFM can be implemented to the manufactured additive components with high roughness [236]. However, the obtained geometry may not meet the dimension required because of the inhomogeneous material removal (MR). A solution has been developed by the proposed study [251] to predict the material removal distribution of a component. During the design, the materials are compensated as they are added to the component. This method's advantage and the feasibility were shown on a laser sintered test similar to the nozzle guide vane (NGV). The CFD model was used in this study to predict compensation for the material removal [6,[251][252][253][254][255].
The materials of the rock are widely used in several industrial applications. The mechanical and physical properties of natural rock differ from one another. To explore the removal mechanism besides rock breakage during the machining process, several experiments were performed, including linear cutting, disc cutting, and drilling. The simulation and modeling technologies have proved their efficiencies in exploring the rock material fracture mechanism except for the experiments. These technologies include the discrete element method (DEM), finite element method (FEM), and hybrid finitediscrete element method (FDEM). The FDEM principles are shown in Figure 51. First, the object was divided into solid elements, and then between the adjacent solid elements, the cohesive elements with zero thickness were inserted. The solid elements determine the combined model deformation, and the cohesive elements control the failure process. Traction separation law (TSL) describes the relationship between the relative displacement and the stress of the cohesive elements during the simulation [256]. element method (DEM), finite element method (FEM), and hybrid finite-discrete element method (FDEM). The FDEM principles are shown in Figure 51. First, the object was divided into solid elements, and then between the adjacent solid elements, the cohesive elements with zero thickness were inserted. The solid elements determine the combined model deformation, and the cohesive elements control the failure process. Traction separation law (TSL) describes the relationship between the relative displacement and the stress of the cohesive elements during the simulation [255]. Figure 51. The principles of FDEM [255]. Figure 51. The principles of FDEM [256].
The schematic diagram of the formation of the crack under tensile load is represented in Figure 52 When the lower and upper surfaces of these elements are stretched. Stress occurs to prevent separation. The damage occurs in the cohesive elements when the stress exceeds its maximum value. At a certain value of the damage, this element will be damaged completely, and the displacement between the lower and the upper surface of these elements is named final failure displacement. The schematic diagram of the formation of the crack under tensile load is represented in Figure 52 When the lower and upper surfaces of these elements are stretched. Stress occurs to prevent separation. The damage occurs in the cohesive elements when the stress exceeds its maximum value. At a certain value of the damage, this element will be damaged completely, and the displacement between the lower and the upper surface of these elements is named final failure displacement.

Mathematical Modeling of MAF Process
The size effect of MAP on the material removal besides surface finished was studied, and it has been found that the forces acting in X and Y directions are given by [28].
and = ℎ( ℎ/ ) V represents the ferromagnetic particle's volume, xm represents the susceptibility of MAPs, and h represents the magnetic field strength at point A, as shown in Figure 53 [28].

Mathematical Modeling of MAF Process
The size effect of MAP on the material removal besides surface finished was studied, and it has been found that the forces acting in X and Y directions are given by [28].
and Fy = xm Vh(dh/dy) V represents the ferromagnetic particle's volume, x m represents the susceptibility of MAPs, and h represents the magnetic field strength at point A, as shown in Figure 53 and = ℎ( ℎ/ ) V represents the ferromagnetic particle's volume, xm represents the susceptibility of MAPs, and h represents the magnetic field strength at point A, as shown in Figure 53 [28]. Figure 53. The magnetic forces that act in the X and Y directions on abrasive particles [256]. Figure 53. The magnetic forces that act in the X and Y directions on abrasive particles [257].
The magnetic pressure can be calculated by P = 3B 2 Wπ 4u0 ((ur − 1) (3(ur + 2) + π(ur − 1)W)) (4) where B represents the density of magnetic flux ur represents the reliability, and W represents the volumetric ratio. The force acting on the abrasive can be calculated through the following equation [1]. ∆F = KD/n (5) D is the MAP diameter, K is constant, and the total is represented by n. the model representing the output response is represented in MR and R and can be calculated through these equations [33].
and R a = R 0 − ΠDPtNnNs 2ll lw where h d represents the indentation diameter, dm represents the magnetic particle diameter, li represents the finishing area length and the width is l w , Ns represents the rotational speed and R 0 is the initial surface roughness [28].

Mathematical Modeling of MRF Process
The parameters that affect the process are shown in Figure 54  The magnetic interaction force is measured using the following equation.
where M represents the intensity, r represents the radius of CIP, P represents the distance between the CIPs, and K is a constant [9]. The shear stress can be calculated from the following equation where H represents the strength of the magnetic field, h is the gab, and is the plastic viscosity. The Centrifugal Force is given by [9].
where rb is the mean radius, Pp is the density of the CIPs, w is the rotational speed and rh is the radius of the MR fluid. The material removal can be expressed as [9].
where Bnd and BCL are coefficients, is CIP size, Cnd represents nondiamond concentration, C is CIP concentration, D represents the penetration depth of the MR fluid and represents the size of non-diamond particles [9]. The magnetic interaction force is measured using the following equation. F = u0 π 9 r 2 KM P 2 (8) where M represents the intensity, r represents the radius of CIP, P represents the distance between the CIPs, and K is a constant [9]. The shear stress can be calculated from the following equation where H represents the strength of the magnetic field, h is the gab, and η is the plastic viscosity. The Centrifugal Force is given by [9]. f (r) = P p 4r 3 rh − r 2rb rh + r 2rb ω 2 (10) where rb is the mean radius, P p is the density of the CIPs, w is the rotational speed and rh is the radius of the MR fluid. The material removal can be expressed as [9]. VRR = Bnd· ϕnd −1/3 · Cnd 1/3 + BCL· ϕCL −4/3 · CCL E/KHv 2 tDV\ where B nd and B CL are coefficients, ϕCL is CIP size, C nd represents nondiamond concentration, C is CIP concentration, D represents the penetration depth of the MR fluid and ϕCL represents the size of non-diamond particles [9].

Mathematical Modeling of EEM Process
Isoviscous-elastic (IE) regime: the solid's elastic deformation is important relative to the fluid film thickness; however, the viscosity cannot be increased due to the low contact pressure. Viscous-elastic regime: the solid's elastic deformation is important relative to the fluid film thickness that separates them in electrohydrodynamic lubrication; besides that, the contact pressure is high, and thus the viscosity is increased. The minimum thickness of the fluid can be described by [28]. hIR = CIR W 2 µ ω 2 (11) and hIE = CIE ω 0.66 µ 0.66 W −0.211 E where the normal load is represented by w, E represents the effective elastic modulus, C IE and C IR are geometric parameters. The shear stress in different regions can be represented by tIR = r ω −1 w 2 CIR −1 µ −1 (13) tIE = r ω 0.34 w 0.211 CIR −1 µ 0.34 E 0.4422 (14)

Mathematical Modeling of MRAFF Process
From Figure 15, the cross-sectional area of penetrated abrasive particle can be obtained through the following equation where D is the diameter of the abrasive particle and t is the depth of indentation. The volume of the material removed can be calculated through the following equation where l w represents the workpiece total length and R i a represents the surface roughness after i cycles [34].

Mathematical Modeling of Ball-End Type MRF (BEMRF) Process
The force on the CIPs can be calculated from the following equation where B represents the density of the applied magnetic field, m represents the CIPs mass, and ∇B represents the gradient of the magnetic field. The "depth of indentation" is calculated through the following equation [33]: where D represents the abrasive particle diameter and D i is the indentation diameter. The final roughness after performing the process is then can be calculated as follows [33]:

Modeling and Simulation of AFM Process by FEM Method
Understanding the process is essential to successfully applying the process as it provides the optimum results with a minimum effort. There are three main components in abrasive flow finishing (AFF): setup, tooling, and medium. The process modeling can be performed mainly in four ways: theoretical, numerical simulation, soft computing, and empirical. The success of modeling the process lies in matching its results with the results obtained from the experiments. To reduce the complexity of AFF, several assumptions are made by the authors in three governing equations, namely momentum, continuity, and constitutive equations of steady-state and axisymmetry, that are used to build the model of FEM. After applying the weighted residual method of Galerkin and weighing functions, the governing equations take a weak formulation. After applying the boundary conditions, the final global FEM equation is [28].
Here, B represents the matrix with different coefficients of the shape functions for the velocities. Np represents the column vector for the approximation of pressure, Nb represents the one-dimensional matrix of the functions of biquadratic shape, and tb represents the traction vectors. The generated stresses in the medium can be obtained by the finite element analysis, which is used to calculate the MR rate besides the roughness of the surface. A simplified abrasive particle model is shown in Figure 15a. It is considered that the surface MR is taking place from the action of micro-cutting. The roughness is considered uniform and triangular as shown in Figure 15b. The total volumetric MR in ith stroke can be calculated through the following equation [28]. Vi = 2π Na Ls Rm 2 Rw dp 2 4 sin −1 2 √ t(dp−t) dp − t(dp − t) dp (25) where Rw represents the workpiece radius, Rm represents the medium cylinder radius, Na represents the active abrasive particle number, Ls is the length of the stock, do is the diameter of the abrasive particle, and t represents the depth. This model is validated by comparing its results with the experimental results. It has been found that the finite element method can obtain the stresses at the surface of the workpiece with high accuracy.

Conclusions
The main objective of this paper is to highlight the developments and challenges related to these techniques and their use in flow machining. By examining the various modeling and simulation methods, such as Molecular Dynamic Simulation (MDS), Computational Fluid Dynamics (CFD), Finite Element Method (FEM), Discrete Element Method (DEM), Multivariable Regression Analysis (MVRA), Artificial Neural Network (ANN), Response Surface Analysis (RSA), Stochastic Modeling and Simulation by Data Dependent System, this study contributes to a better understanding (DDS). CFD and FEM can be carried out using commercially available software. In contrast, DEM and MDS are carried out using a platform based on computer programming, such as "LAMMPS Molecular Dynamics Simulator" or C, C++, or Python programming. The modeling and simulation of loose abrasive-based machining operations seem to be more successful with these techniques. The other four methods-MVRA, ANN, RSA, and DDS-are experimental and based on statistical techniques that can be used to represent loose abrasive-based machining operations mathematically. Additionally, it includes a priceless bibliography of past works on the modelling and simulation methods for abrasive-based machining processes as well as areas for additional research. Researchers studying mathematical modeling of various micro-and nano-finishing techniques for different applications may find this review article to be of great help. The process of loose abrasive-based machining is an emerging field in machining. The simulation and modeling of this process are still in a budding phase. Several studies investigated on modeling and simulation of this process using different techniques like MDS, CFD, FEM, DEM, MVRA, ANN, RSA, and DDS. The following conclusions could draw based on the present study of various simulation techniques for abrasive-based machining processes: 1.
MDS simulation is used to predict how disruption will affect a molecular system. So, it simulated the results in depth observations. 2.
The accuracy of MDS simulation is much higher than other computational techniques i.e., FEM and CFD simulation, because it predicts the simulation results at a molecular level.

3.
CFD simulation is used to study intricate fluid-fluid, fluid-solid, or fluid-gas interaction problems. So, the CFD simulation can simulate any kind of abrasive-based machining process involving the interaction fluid and solid domains. 4.
CFD simulation predicted the simulation results to observe the flow field in depth and its entirety without any interference. So, the accuracy of CFD simulation results is much more accurate than the FEA simulation results. 5.
FEM simulation is used to study any physical problem by the system's disintegration, i.e., meshing, applying boundary conditions, and solving the equations using commercial FEM simulation software, i.e., Ansys, Abaqus, etc. 6.
DEM method is similar to the MDS simulation technique, and it allows a more detailed study of the micro-dynamics of powder flows than is often possible using physical experiments. 7.
Many researchers have adopted hybrid methods, i.e., DEM-FEM or DEM_CFD coupled method, to solve the complex modeling problem with high accuracy. 8.
MDS, DEM, and CFD need high configured and advanced computation platforms to perform the fast simulation. 9.
The accuracy of simulation results obtained from MDS and DEM is much higher than CFD or FEM simulation results. 10. Out of the different techniques explored in the purpose study, MDS and DEM seem to be the most promising technology for this process due to their higher accuracy in the results compared to the other techniques. 11. Like applying MDS in other domains, abrasive grain particle tracking can also be performed using this approach, which predicts the machining behaviors. 12. However, higher computational power is required for the DEM and MDS than the other approaches limiting their use and application. 13. MVRA, ANN, RSA, DDS are the experimental based statistical techniques to predict the results of the investigation of abrasive-based machining processes 14. The availability of the statistical techniques helps determine the parameters that affect the procedure as well as the extent of their influence. 15. The research in the purpose study provides knowledge to explore and predict the optimum process parameters for abrasive-based machining processes such as AFM, MRF, MAF, MRAFF, and BEMRF.