Parameter Identification of the Fractional-Order Mathematical Model for Convective Mass Transfer in a Porous Medium

Fractional calculus is an essential tool in studying new phenomena in hydromechanics and heat and mass transfer, particularly anomalous hydromechanical advection–dispersion considering the fractal nature of the porous medium. They are valuable in solving the urgent problem of convective mass transfer in a porous medium (e.g., membranes, filters, nozzles, convective coolers, vibrational prillers, and so on). Its solution allows for improving chemical engineering and technology workflows, refining process models for obtaining porous granular materials, realizing the convective cooling of granular and grain materials, and ensuring the corresponding apparatuses’ environmental safety. The article aims to develop a reliable convective mass transfer model for a porous medium and proposes a practical approach for its parameter identification. As a result, a general scientific and methodological approach to parameter identification of the fractional convective mass transfer model in a porous medium was proposed based on available experimental data. It mainly used Riemann–Liouville fractional time and coordinate derivatives. The comprehensive application of the Laplace obtained the corresponding general solution transform with respect to time and a coordinate, the Mittag-Leffler function, and specialized functions. Different partial solutions in various application case studies proved this solution. Moreover, the algorithm for practically implementing the developed approach was proposed to evaluate parameters for the considered model by evaluation data. It was reduced to the two-parameter model and justified by the available experimental data.


Introduction
Research on the transport phenomena across porous membranes became a significant problem many years ago [1].This process was accompanied by the development of membrane science [2].Successful attempts to advance solving the problem of convective mass transfer in porous membranes have gained more importance recently [3].Differential equations that describe the convective mass transfer process present significant difficulties for their solution [4].An analytical solution to these equations can be practically expedient only for cases with significant simplifications [5].
Therefore, approaches based on the similarity theory are frequently used to solve this problem [6].However, using algebraic criterion equations significantly limits understanding of the physical meaning of mass transfer processes [7].
Moreover, using porous materials predetermines the need to consider the fractional derivatives in mass transfer equations [8].Standard partial differential equations cannot study relevant phenomena occurring in such media [9].
In this regard, the problem of considering fractional-order partial differential equations and their subsequent analytical solution for specific boundary and initial conditions is a fundamental scientific problem.Moreover, the practical significance of solving the parameter identification problem of the fractional convective mass transfer model in a porous medium is highlighted by the need to intensify heat and mass transfer processes in energy-efficient apparatuses of the chemical industry (e.g., convective coolers for granular and grain materials, vibrational prilling of fertilizers, and so on).
Numerous key publications in the related fields also highlight the significance of this problem.Li and Kurths [10] proposed a numerical algorithm for solving the fractional convective diffusion equation.As a result, the behavior of the fractional convection model was analyzed.Salomoni and De Marchi [11] described a nonlinear anomalous hydromechanical advection-dispersion model.As a result, fractional advection-diffusion equations using the Grünwald-Letnikow definition of a fractional derivative were solved numerically.
Allwright and Atangana [12] considered the fractional-order advection-dispersion equation for groundwater transport.As a result, the fractional integral of a system with self-similarities was calculated numerically.Anwar et al. [13] investigated heat transfer in a nanofluid flow through a porous medium.As a result, the friction coefficient and Nusselt number were calculated numerically for the fractional derivative model.
Madhura et al. [14] studied the impact of nanoparticle shapes on natural convection flow with heat and mass transfer rates of nanofluids with fractional derivatives.As a result, the governing fractional-order differential equations were solved using the Laplace transform with boundary conditions.However, this model does not consider fractional derivatives for porous media.
Borah et al. [15] numerically analyzed thermal transfer equations for non-Newtonian fluids.Based on Atangana-Baleanu and Caputo-Fabrizio fractional derivatives, the mixed convective unsteady flow through an infinite vertical porous plate was investigated for a uniform field.Li et al. [16] applied a fractional lattice Boltzmann method to solve multiphase flows with complex interfacial behavior and significant density contrast.As a result, a good agreement was obtained between the numerical results and the available numerical data.
Raza et al. [17] proposed a non-singular fractional approach for the thermal analysis of natural convection in nanofluids.As a result, an impact of the fractional parameter on temperature and velocity fields was obtained graphically.Khan et al. [18] considered free convection flow between two parallel plates.As a result, the coupled system of Fick's and Fourier's equations was solved numerically using Python 3.11.5 software.
Rehman et al. [19] proposed generalized Mittag-Leffler solutions of free convection heat and mass transfer flow of the Maxwell fluid with Newtonian heating.As a result of using the Prabhakar fractional derivative approach, a comparison between velocity and temperature fields for Maxwell and Newtonian fluids for both fractional and classical cases with and without slip conditions was provided.A similar problem was solved semi-analytically by Riaz et al. [20].
Ali et al. [21] solved a generalized magnetohydrodynamic two-phase free convection flow model of dusty Casson fluid between parallel plates.As a result, the temperature field was plotted using the MathCAD 14 software.However, these studies did not consider the fractional derivatives with respect to time coordinates in the governing equations.Meng et al. [22] numerically solved the fractional viscoelastic fluid problem.As a result of using the finite difference approach, the complexity of solving the coupled system of temperature and concentration equations was highlighted.However, fractional convection was not considered.
Convection in a porous medium is influenced by a significant number of factors (hydrodynamic, temperature, and so on).Their comprehensive impact on the features of the operating processes can be initially determined by the criteria of similarity (e.g., Reynolds, Nusselt, Prandtl) given in the work [23].However, the integer-order derivatives describing the acting forces in the process do not allow for considering either the phenomenon of anomalous convection at low Prandtl numbers described in the work [24] or the fractal nature of the medium [25].Therefore, the fractional-order mathematical model should be considered, and the parameter identification approach by available experimental results data should be developed.
Due to the abovementioned analysis, the following research gaps should be eliminated despite the numerous studies in this field.Firstly, preconditions for possible analytical solutions of the fractional convective mass transfer equation should be considered.The available partial solutions should justify these solutions for a continuous medium and integer time dimension.Secondly, a scientific approach for parameter identification of the stationary model should be applied.The corresponding methodology for evaluating the fractional order by experimental data of concentration changes should be proposed.It should also be extended for the case of the non-stationary model.
Therefore, the main aim of the research is to develop a reliable convective mass transfer model for a porous medium and propose an effective approach for its parameter identification based on available experimental data.

A Flowchart of the Research
A flowchart of the research methodology is presented in Figure 1.It consists of the following three main stages: analytical (both classical and a proposed new one), experimental (obtaining results in a hydrodynamic apparatus), and numerical (evaluating the parameters of the fractional order mathematical model).describing the acting forces in the process do not allow for considering either the phenomenon of anomalous convection at low Prandtl numbers described in the work [24] or the fractal nature of the medium [25].Therefore, the fractional-order mathematical model should be considered, and the parameter identification approach by available experimental results data should be developed.Due to the abovementioned analysis, the following research gaps should be eliminated despite the numerous studies in this field.Firstly, preconditions for possible analytical solutions of the fractional convective mass transfer equation should be considered.The available partial solutions should justify these solutions for a continuous medium and integer time dimension.Secondly, a scientific approach for parameter identification of the stationary model should be applied.The corresponding methodology for evaluating the fractional order by experimental data of concentration changes should be proposed.It should also be extended for the case of the non-stationary model.
Therefore, the main aim of the research is to develop a reliable convective mass transfer model for a porous medium and propose an effective approach for its parameter identification based on available experimental data.

A Flowchart of the Research
A flowchart of the research methodology is presented in Figure 1.It consists of the following three main stages: analytical (both classical and a proposed new one), experimental (obtaining results in a hydrodynamic apparatus), and numerical (evaluating the parameters of the fractional order mathematical model).The first (analytical) stage allows for identifying the problems with the traditional model (e.g., inability to describe certain hydrodynamic phenomena) and observing ways for its further extension.The first (analytical) stage allows for identifying the problems with the traditional model (e.g., inability to describe certain hydrodynamic phenomena) and observing ways for its further extension.
The second (experimental) stage allows for obtaining experimental results data for the practical case study (e.g., pneumatic classification of porous material in a rhomb-shaped hydrodynamic apparatus).
The last (numerical) stage allows for approximating the evaluated results data to make fundamentals of the parameter identification of the proposed fractional-order mathematical model.It should also justify the reliability of this model in terms of its compliance with available data.

The Traditional Model
The general convection-diffusion equation in a vector form is as follows [26]: where c( → r , t)-mass concentration field in a vector form, kg/kg; → r -radius-vector, m; t-time, s; Despite the complexity of this spatial-dimension model, the following one-dimensional model with unchanged velocity U 0 and diffusion coefficient D 0 is practically applied in studying problems of filtration through porous medium [27], the energy efficiency of liquidvapor ejectors [28], inertial-filtering separation of gas-liquid mixtures [29], and pneumatic classification of granular material [30]: where z-coordinate, m; c(z, t)-mass concentration, kg/kg; R(z, t)-specific driving force, kg/s.However, even this equation does not consider the fractional origin of the process in a porous medium.Therefore, it should be rewritten considering Riemann-Liouville fractional derivatives with respect to time and dimension for a porous medium despite standard integer dimensions.
The initial equations contain the integer-order partial differential equations.However, these models do not allow for considering either the anomalous convection phenomenon or the medium's fractal structure.This is because such phenomena can be studied by considering the fractional-order derivatives.Therefore, these equations should be extended to a general fractional-order mathematical model using fractional derivatives.Moreover, a scientific approach should be developed to identify the model's parameters by available experimental results data.

The Considered Fractional-Order Convection Model
For the reasons mentioned above, the following assumptions and hypotheses are proposed.Firstly, the process is predominantly convective rather than diffusive for significant values of the comprehensive criterion Re•Pr, where Re-Reynolds number; Pr-Prandtl number.In this case, the diffusion component can be neglected compared to the convective one.
Secondly, for a fractional convective process, the time derivative of the mass concentration should be rewritten as a fractional derivative ∂t α , where α = (0, 1]-the fractional order (for the standard integer order α = 1).Additionally, for a porous medium, the convective mass transfer component should be rewritten as −U 0 ∂z β , where β = (0, 1]-the fractional order (for a continuous medium β = 1); U 0 -convective ratio as a fractional analog of the carrier flow velocity.In this case, the following Riemann-Liouville fractional derivatives [31,32] are considered: where Γ(θ)-Euler's gamma function of the argument θ; τ, ζ-independent time and coordinate parameters, respectively.Fourthly, the arbitrary specific driving force of the mass transfer process can be presented in a linearized form as R(z, t) = K•[c s − c(z, t)], where K-mass transfer coefficient, kg/s; c s -limiting mass concentration of the dispersed phase under the saturation, kg/kg.In this case, a constant value of the mass transfer coefficient K is usually considered for the interfacial surface.
Moreover, since the convective mass transfer process occurs due to the change in limiting s s and current c concentrations and simplifying the below-mentioned transformations, zero boundary and initial conditions can be used: c(0, t) = c(z, 0) = 0. Furthermore, these conditions will, in the future, allow us not to complicate the problem with the traditional procedure of eliminating singularities associated with the application of Riemann-Liouville derivatives.
Finally, the stationary mode occurs for significant values of the comprehensive criterion Fo•Pr  Re , where Fo-diffusion Fourier number.This result entirely corresponds to the stationary convection equation describing the steady-state behavior of such systems.In this case, the time derivative can be neglected, and the concentration field leads to the limiting function c ∞ (z) with the analogous boundary condition c ∞ (0) = 0.
Under the assumptions mentioned above, the following differential equations as the fractional convective mass transfer models in a porous medium are considered: These equations should be solved analytically for the above boundary and initial conditions.

Analytical Solution of the Non-Stationary Equation
The consequent application of the Laplace transform is used to solve fractional convective mass transfer Equations ( 4) and (5).
Equation ( 4), considering the mentioned boundary and initial conditions, takes the following form: where s, τ-derivative operators with respect to longitudinal direction and time, respectively; C(s, τ) = L t {L s [c(z, t)]}-Laplace transform with respect to both coordinate and time.

Justification of the Obtained Solutions
Due to the complexity of the obtained general solution, its reliability should be proven for the following four reduced case studies.The first considers the traditional approach when α and β are integer numbers, and the second is for its stationary mode.The last two justifications should be based on reducing the general solution (8) to the stationary mode and the case without significant convection.
In justification of the reliability of the obtained solutions, Equation ( 8) should be a reduction to the traditional approach when α = β = 1.In this case, the following solution can be obtained: where H(t)-Heaviside function.This solution can be represented graphically in the following dimensionless form (Figure 2) for the case of K ≥ 0 and U 0 ≥ 0: where where Eα,1, Eβ,1-Mittag-Leffler functions: where Γ  1    -Euler's gamma function; i = [0, ∞ )-integer numb θ-independent variable.

Justification of the Obtained Solutions
Due to the complexity of the obtained general solution, its reliability should proven for the following four reduced case studies.The first considers the traditional a proach when α and β are integer numbers, and the second is for its stationary mode.T last two justifications should be based on reducing the general solution (8) to the statio ary mode and the case without significant convection.
In justification of the reliability of the obtained solutions, Equation (8) should be reduction to the traditional approach when α = β = 1.In this case, the following soluti can be obtained: where H(t)-Heaviside function.This solution can be represented graphically in the following dimensionless for (Figure 2) for the case of  0 and  0: Notably, the reduced solution (10) entirely corresponds with the available soluti [33] of Equation (2).Therefore, the obtained general solution of the fractional convecti mass transfer model in a porous medium is valid.Notably, the reduced solution (10) entirely corresponds with the available solution [33] of Equation (2).Therefore, the obtained general solution of the fractional convective mass transfer model in a porous medium is valid.
Moreover, for the stationary mode (t ≥ z U 0 ; c = c ∞ (z) = invar(t)), the last equation takes the following form: that corresponds to the saturation equation [34].This fact also proves the obtained general solution (8).
For the stationary mode, when t → ∞, two last inverse Laplace transforms in Equation (8) are reduced to 0, and concentration c(z, t) is reduced to its limiting function c ∞ (z) = lim t→∞ c(z, t).Therefore, the following solution for the stationary mode is obtained: This solution can also be obtained directly from Equation (5).Particularly, after the direct Laplace transform, that equation considering the boundary mentioned above condition takes the following form: where ]-the Laplace transform with respect to longitudinal coordinate.
The following result can be derived from the last equation: After the inverse Laplace transform with respect to coordinate, the stationary solution ] takes the same form as Equation ( 13).This fact finally proves the reliability of the general solution (8).
For the case of insufficient convection (z → ∞), Equation ( 4) can be rewritten as follows: where c 0 (t) = lim z→∞ c(z, t)-limiting concentration change in time.
After the direct Laplace transform, this equation considering the boundary condition mentioned above, takes the following form: where C 0 (τ) = L τ [c 0 (t)]-the Laplace transform with respect to time.
The following result can be derived from the last equation: After the inverse Laplace transform with respect to coordinate, the stationary solution c 0 (t) = L −1 τ [C 0 (τ)] takes the following form: This solution is a partial case of the general solution (7) reduced by the substitution of z → ∞, and c 0 (t) = lim z→∞ c(z, t).This fact finally proves the reliability of the general solution (7).

Parameter Identification of the Fractional-Order Model
For practical purposes, Equation (4) needs to evaluate five parameters: the limiting concentration c s , fractional orders α and β, as well as a convective ratio U 0 and mass transfer coefficient K.The following algorithm for their identification is proposed.Firstly, limiting concentration c s is evaluated under the saturation condition c s = lim z→∞ c ∞ (z).
Secondly, Equation ( 17) is rewritten in a dimensionless form (Figure 2): where c 0 (t) = c 0 (t) c s -dimensional time-varying concentration.Fractional order α and mass transfer coefficient K are evaluated based on the best fit of this solution with an available dataset of experimental parameters c <k> 0 for the k-th time moment t k (k = 1, N, where N is the total number of the time-varying dataset).For this purpose, the following least square error function should be minimized: Such minimization is based on the following necessary and sufficient conditions: After equal transformations and considering the first Expression ( 9), the following system of nonlinear equations can be obtained: where Ψ(θ) = d dθ {ln[Γ(θ)]}-the specialized function as a derivative of the natural loga- rithm of Euler's gamma function.
Equation ( 23) allows for identifying both the parameters α and K.
Next, the evaluated parameter K is substituted for Equation (13), which can be rewritten in the following dimensionless form: where c ∞ (z) = c ∞ (z) c s -dimensional concentration for the stationary mass transfer process.Fractional order β and convective ratio U 0 are evaluated based on the best fit of this stationary solution with an available dataset of experimental parameters c <j> ∞ for the j-th number of the coordinate z j (j = 1, n, where n is the total number of the dataset).For this purpose, the following least square error function should be minimized: Such minimization is based on the following necessary and sufficient conditions: Membranes 2023, 13, 819 9 of 16 After equal transformations and considering the second Expression (8), the following system of nonlinear equations can be obtained: This system allows for identifying both the parameters β and U 0 .Therefore, parameter identification of the main parameters (α, K, β, and U 0 ) of the fractional convective mass transfer model in a porous medium is based on solving Equations ( 22) and (26).However, due to the complexity of this approach, it should be generalized for both cases as below.

The Generalized Parameter Identification Approach
Despite the difference between these expressions, there is a way to present them similarly (Figure 3).For this purpose, the following single representation of solutions ( 20) and ( 24) is used: where κ-generalized coefficient; γ-index.
After equal transformations and considering the second Expression (8), the following system of nonlinear equations can be obtained: This system allows for identifying both the parameters β and  .Therefore, parameter identification of the main parameters (α, K, β, and  ) of the fractional convective mass transfer model in a porous medium is based on solving Equations (22) and (26).However, due to the complexity of this approach, it should be generalized for both cases as proposed below.

The Generalized Parameter Identification Approach
Despite the difference between these expressions, there is a way to present them similarly (Figure 3).For this purpose, the following single representation of solutions ( 20) and ( 24) is used: where κ-generalized coefficient; γ-index.If coefficient θ = K, index γ = α, and independent variable θ = t, this solution corresponds to (20).Otherwise, in the case of θ = K U 0 , γ = β, and θ = z, corresponding to (24).Therefore, generally, the parameter identification of the fractional convective mass transfer model in a porous medium is based on the minimization of the following least square error: where θ j , c j -experimental datasets for independent parameters and the dimensionless concentration, respectively.
Such minimization is based on the following necessary and sufficient conditions: or in the following form: After considering the following derivatives from the Mittag-Leffler function: the system of the nonlinear Equation ( 31) can be rewritten like ( 23) and ( 27): This system of equations allows for identifying both the parameters γ and κ.Therefore, the generalized scientific and methodological approach for evaluating the main parameters of the fractional convective mass transfer model in a porous medium is based on solving the generalized system of the nonlinear Equation (33).

The Evaluation Algorithm
The significance of such an approach is associated with the following algorithm of its realization.The initial value γ 0 = 1 of the parameter γ is chosen in the initial step.Then, the first approximation κ 1 of the parameter κ is evaluated numerically from the second Formula (33).Finally, the first approximation γ 1 of the parameter γ is evaluated numerically from the first Formula (33).This cycle should be realized during p ≥ 1 step till condition (29) is fulfilled.
A rhomb-shaped pneumatic classifier was considered to obtain the evaluation dataset.Such an apparatus allows for effectively separating samples of granular materials (porous ammonium nitrate NH 4 NO 3 , superphosphate Ca(H 2 PO 4 ) 2 •H 2 O, and carbamide CO(NH 2 ) 2 ), reaching relatively high values of the target fraction.
The experimental setup is presented in Figure 4.It consists of a pneumatic classifier, cyclone, gas blower for creating a suction flow, regulating valve, and separate tanks for collecting coarse and fine fractions.
A rhomb-shaped pneumatic classifier was considered to obtain the evaluation dataset.Such an apparatus allows for effectively separating samples of granular materials (porous ammonium nitrate NH4NO3, superphosphate Ca(H2PO4)2•H2O, and carbamide СО(NH2)2), reaching relatively high values of the target fraction.
The experimental setup is presented in Figure 4.It consists of a pneumatic classifier, cyclone, gas blower for creating a suction flow, regulating valve, and separate tanks for collecting coarse and fine fractions.Sieves of various sizes are used to obtain the different charges.The Digimizer 6.3.0 software can be additionally applied to decrease the root mean square deviation up to 0.014 mm.
The pressure and the airflow in the apparatus at different cross-sections were measured by the U-shaped liquid pressure gauges and the Pitot-Prandtl tube (Sumy, Ukraine), respectively.A sloping tube with a micromanometer MMN-240(5)-0.1 (Poznan, Poland) allowed for controlling the pressure difference at the cross-sections.
The experimental data (concentration for the considered case study) were determined by weighting batches of the porous granular material using the electronic scale Momert-6000 (Bratislava, Slovak Republic) with an accuracy of 0.1 g.
The dataset with the total number of the generated experimental n = 30 was considered for the case study of the pneumatic classification of the porous granular material.
The evaluation results are summarized in Table 1.Sieves of various sizes are used to obtain the different charges.The Digimizer 6.3.0 software can be additionally applied to decrease the root mean square deviation up to 0.014 mm.
The pressure and the airflow in the apparatus at different cross-sections were measured by the U-shaped liquid pressure gauges and the Pitot-Prandtl tube (Sumy, Ukraine), respectively.A sloping tube with a micromanometer MMN-240(5)-0.1 (Poznan, Poland) allowed for controlling the pressure difference at the cross-sections.
The experimental data (concentration for the considered case study) were determined by weighting batches of the porous granular material using the electronic scale Momert-6000 (Bratislava, Slovakia) with an accuracy of 0.1 g.
The dataset with the total number of the generated experimental n = 30 was considered for the case study of the pneumatic classification of the porous granular material.
The evaluation results are summarized in Table 1.The dataset and the corresponding parameter identification results are presented in Figure 5.The dataset and the corresponding parameter identification results are presented in Figure 5.The following parameters have been identified after p = 9 approximation cycles: γ = 0.587 and κ = 0.301.In this case, the minimized least square error (28) equals 0.038.
Overall, the reliability of the proposed procedure is proven by small values of the averaged relative error between the initial and evaluated datasets of concentrations.Particularly, the average relative error is equal to 2.1%.

Discussion
The obtained results significantly extend previous studies in the mathematical modeling of working processes in mass transfer equipment.Mainly, the mathematical model of the concentration changes [33] can been improved in considering the proposed model's fractional origin.This also makes it possible to improve practical ways for ensuring the reliability of the linear fractional advection-diffusion, advection-dispersion, and convection-dispersion models [35], i.e., through heteroporous membranes [36,37].It also allows for improving the pneumatic classification process for granular material in rhomb-shaped The following parameters have been identified after p = 9 approximation cycles: γ = 0.587 and κ = 0.301.In this case, the minimized least square error (28) equals 0.038.
Overall, the reliability of the proposed procedure is proven by small values of the averaged relative error between the initial and evaluated datasets of concentrations.Particularly, the average relative error is equal to 2.1%.

Discussion
The obtained results significantly extend previous studies in the mathematical modeling of working processes in mass transfer equipment.Mainly, the mathematical model of the concentration changes [33] can been improved in considering the proposed model's fractional origin.This also makes it possible to improve practical ways for ensuring the reliability of the linear fractional advection-diffusion, advection-dispersion, and convectiondispersion models [35], i.e., through heteroporous membranes [36,37].It also allows for improving the pneumatic classification process for granular material in rhomb-shaped apparatuses [38], intensifying the cooling process during obtaining granular materials [30], decreasing dust emission by monodisperse system technology [39], and reducing the technogenic impact on the environment from the emissions of heat power engineering by using highly efficient equipment [40].
The research findings are also helpful in solving numerous practical problems with mass transfer equipment.Particularly, the obtained solution allows for identifying parameters of the fractional-order mathematical model of particle sedimentation [41], convective heat transfer in porous insulation layers [9], and hydrodynamics of viscoelastic fluid flows [22].They are also applicable in practice to solve the coupled problem of heat and mass transfer with non-stationary boundaries for the consequent modeling of the oscillation's propagation in a working medium of a vibrational priller.
Moreover, the applied model can be practically applied to extend the application of a functionally graded fractional porous medium in biomechanical applications [42].
The practical significance of the obtained results corresponds to the significance of the convective mass transfer model in terms of theoretical substantiation of mathematical models of oil filtration through a porous medium [27].It also allows the implementation of the developed model to design inertial-filtering separators [29] with the consequent evaluation of their operating parameters.
The theoretical value of the research is in substantiating the scientific fundamentals of evaluating orders of fractional derivatives of the mathematical model for convective mass transfer.The developed parameter identification approach will allow for evaluating the porosity and fractional dimension of porous materials based on experimental data.
The limitations of the study are mainly predetermined by the case of a predominantly convective mode for significant values of the multiplication of the Reynolds and Prandtl numbers.However, further extending the proposed model to diffusion and Taylor dispersion in higher space dimensions using, e.g., Riesz-Feller fractional operators [43], will allow for eliminating this limitation.
More precise studies of the impact of operating parameters on the size of pores, as well as the physical and chemical properties of the resulting granules, will also be carried out further by using an electron microscope.

Conclusions
Thus, a fractional mathematical model of the convective mass transfer has been developed.It extends the available models and makes the fundamentals in considering new phenomena in hydromechanics and heat and mass transfer, considering the fractal nature of the porous medium.
During the study, the fractional convective mass transfer equation was solved analytically based on the Laplace transform with respect to time and coordinate and further representation using the Mittag-Leffler function and Duhamel integral.
The obtained generalized solution has been justified by four available partial solutions (i.e., integer parameters, stationary mode, saturation mode, and without significant convection).
The experimental results have been realized for the case of a hydrodynamic classifier for porous granular materials.They have been approximated by analytical curves.Such an approach has allowed for identifying parameters of the proposed fractional-order mathematical model.
As a result, the stationary and non-stationary convective mass transfer models in a porous medium have been solved.The corresponding parameter identification approach has been developed.Moreover, the impact of the fractional orders for a porous medium on the behavior of spatial surface and plan curves has been presented graphically for the dimensionless concentration.The reliability of the proposed approach has been proven by the average error of the approximation equal to 2.1%.
Overall, a scientific and methodological approach to parameter identification of the fractional convective mass transfer model has been developed.It has allowed for evaluating the parameters of this model by an experimental dataset of concentration changes during pneumoclassification of porous granular materials.The corresponding algorithm of such a parameter identification has been described.The evaluation example has proved the developed approach with the least square error of 0.038.
The scientific novelty of the research is in considering the fractional-order nature of the mathematical model for convective mass transfer, which has created the scientific foundation of further substantiating the anomalous hydromechanical advection-dispersion in porous structures of fractal nature.perspective plan of development of a scientific direction "Technical sciences" Sumy State University" (state registration number 0121U112684).The research was also supported by the International Association for Technological Development and Innovations.

Conflicts of Interest:
The authors declare no conflict of interest.

Figure 1 .
Figure 1.The flowchart of the research.

Figure 1 .
Figure 1.The flowchart of the research.

Table 1 .
The evaluation results data.

Table 1 .
The evaluation results data.