Scientific Research

An Academic Publisher

**The Movement of a Mixture of Cotton with an Air Stream during Pneumatic Transport by Pipeline of Variable Cross Section** ()

Keywords

Share and Cite:

*Engineering*,

**11**, 531-540. doi: 10.4236/eng.2019.118037.

1. Introduction

Depending on the number of components in the pipeline transport, one-component and multi-component media are distinguished. A single-component medium is a medium consisting of a single material or substance, which is called a homogeneous mixture, and a multi-component—of several materials or substances that differ in physical, mechanical, chemical and other properties, is called a heterogeneous mixture [1] [2] . The mathematical description of the motion of a multicomponent medium is complicated by the difference in the reaction of each component to the movement of another. Therefore, to simplify the process of theoretical consideration, multicomponent media are taken as a one-component or two-component medium consisting of a continuous phase—air and a discrete phase—material, which ultimately gives more general and less accurate results.

Especially, when considering the processes of pneumatic transportation of materials, including raw cotton, the theory of motion of a discrete medium in a stationary or moving continuous medium has been adopted to date, which will give a more or less distorted idea of the process of moving air and material inside the pipeline [3] . More suitable, in our opinion, to study the process of pneumatic transportation theory, is the theory of multi-speed systems, according to which, all components are moved through the pipeline separately, with interdependent parameters, including speed, which is a new approach to the study of the processes of pneumatic transportation of materials. To ensure the universality of the laws, it should consider the motion of the mixture inside the pipeline with a variable cross section [4] [5] .

The mixture of air with fiber particles is assumed to be a heterogeneous mixture, to describe the motion of which you can use the theory of multi-speed systems proposed in the work of Kh. A. Rakhmatullin [6] .

Heterogeneous mixtures, as a rule, are described by a multi-speed model taking into account the dynamic effects arising due to the mismatch of the speeds of the individual phases. At the same time, we consider the air to be an ideal fluid and the internal force of interaction is determined through the normal pressure, which is common to the whole mixture, opposite to the direction of movement of the particles of the fibrous medium [7] [8] .

The following designations are used in the research:

${u}_{0}\left(x\right)$ —air velocity (index-0) in an arbitrary section of the pipeline;

${u}_{1}\left(x\right)$ —the speed of the particles of raw cotton (index-1) in an arbitrary section of the pipeline;

${u}_{00}$ —air flow velocity in the pipeline section $x=0$ ;

${\rho}_{0}$ —reduced air density in an arbitrary section of the pipeline, kg/m^{3};

${\rho}_{0}^{\left(0\right)}$ —the true density of air, kg/m^{3};

${\rho}_{00}$ —values of air density in section $x=0$ ;

${\rho}_{1}$ —the reduced density of the fibrous mass in an arbitrary section of the pipeline, kg/m^{3};

${\rho}_{1}^{\left(0\right)}$ —their true density of fibrous mass, kg/m^{3};

${\rho}_{10}$ —values of the density of the fibrous mass in the cross section;

${\rho}_{0}\left(x\right)$ , ${\rho}_{1}\left(x\right)$ and ${u}_{0}\left(x\right)$ , ${u}_{1}\left(x\right)$ —current density and velocity values of the components;

m—the porosity of cotton in an arbitrary section of the pipe;

${m}_{0}$ —porosity of cotton in section $x=0$ ;

k—coefficient of resistance of cotton to air movement, kg/s;

${k}_{\u043f}$ —the boundary value of the coefficient, kg/s;

${p}_{0}$ —the pressure of the mixture in the cross section $x=0$ ;

$q={k}_{0}p$ —lateral flow pressure acting on the walls of the pipeline, PA;

${k}_{0}$ —lateral pressure coefficient;

${F}_{tr}=fqL$ —the friction force between the raw and the surface of the pipe, N;

f—the coefficient of friction between the outside surface of the pipe;

L—length of the pipe cross-section contour), m.

2. Study of the Motion of a Multicomponent Medium

Set the origin of the coordinates in the initial section of the pipeline. We direct the axis along the axis of the pipeline, the cross-sectional area of which varies according to the law, considering the process stationary, denoted by, respectively, the air velocity (index-) and cotton particles of raw mass (index-1) in an arbitrary section of the pipeline. We believe that in the cross section of the pipeline the air flow with speed acts on the moving mass of the mixture. The equations of one-dimensional motion of the components of the mixture and the laws of mass conservation, according to, are written in the form

${\rho}_{0}{u}_{0}s\frac{\text{d}{u}_{0}}{\text{d}x}=-\frac{{\rho}_{0}}{{\rho}_{0}^{\left(0\right)}}\frac{\text{d}\left(sp\right)}{\text{d}x}+sk\left({u}_{1}-{u}_{0}\right)$ (1)

${\rho}_{1}{u}_{1}s\frac{\text{d}{u}_{1}}{\text{d}x}=-\frac{{\rho}_{1}}{{\rho}_{1}^{\left(0\right)}}\frac{\text{d}\left(sp\right)}{\text{d}x}-sk\left({u}_{1}-{u}_{0}\right)$ (2)

${\rho}_{0}{u}_{0}s={u}_{00}{\rho}_{00}{s}_{0}=const$ , ${\rho}_{1}{u}_{1}s={u}_{10}{\rho}_{10}{s}_{0}=const$ ,

${\rho}_{0}=m{\rho}_{0}^{\left(0\right)}$ , ${\rho}_{1}=\left(1-m\right){\rho}_{1}^{(0)}$

${\rho}_{0}=\frac{m}{{m}_{0}}{\rho}_{00}$ , ${\rho}_{1}=\frac{1-m}{1-{m}_{0}}{\rho}_{10}$ (3)

From (3) the following relations follow

${u}_{0}=\frac{{m}_{0}}{m}\frac{{s}_{0}}{s}{u}_{00}$ , ${u}_{1}=\frac{1-{m}_{0}}{1-m}\frac{{s}_{0}}{s}{u}_{10}$ , $\frac{{\rho}_{0}^{\left(0\right)}}{{\rho}_{0}}+\frac{{\rho}_{1}^{\left(0\right)}}{{\rho}_{1}}=1$ (4)

Equations (1) and (2) we bring to the form

${\rho}_{0}^{\left(0\right)}{u}_{0}s\frac{\text{d}{u}_{0}}{\text{d}x}=-\frac{\text{d}\left(sp\right)}{\text{d}x}+\frac{{\rho}_{0}^{\left(0\right)}}{{\rho}_{0}}sk\left({u}_{1}-{u}_{0}\right)$ (5)

${\rho}_{1}^{\left(0\right)}{u}_{1}s\frac{\text{d}{u}_{1}}{\text{d}x}=-\frac{\text{d}\left(sp\right)}{\text{d}x}-\frac{{\rho}_{1}^{\left(0\right)}}{{\rho}_{1}}sk\left({u}_{1}-{u}_{0}\right)$ (6)

Excluding the derivative from the system $\frac{\text{d}\left(sp\right)}{\text{d}x}$ , we get

${\rho}_{0}^{\left(0\right)}{u}_{0}\frac{\text{d}{u}_{0}}{\text{d}x}-{\rho}_{1}^{\left(0\right)}{u}_{1}\frac{\text{d}{u}_{1}}{\text{d}x}=\left(\frac{{\rho}_{0}^{\left(0\right)}}{{\rho}_{0}}+\frac{{\rho}_{1}^{\left(0\right)}}{{\rho}_{1}}\right)k\left({u}_{1}-{u}_{0}\right)$ (7)

Equations (3) and (4) imply that the densities and velocities of the components in an arbitrary section of the pipeline are expressed in terms of porosity $m\left(x\right)$ , which after excluding speeds ${u}_{0}\left(x\right)$ , ${u}_{1}\left(x\right)$ from Equation (7) with the help of relations (4), satisfies the equation:

$\frac{\text{d}m}{\text{d}x}=\frac{{F}_{1}\left(m\right)}{{F}_{0}\left(m\right)}\frac{{s}^{\prime}\left(x\right)}{s\left(x\right)}+\frac{{F}_{2}\left(m\right)}{{F}_{0}\left(m\right)}\frac{{s}_{0}}{s\left(x\right)}$ (8)

where:

${F}_{0}={\rho}_{0}^{\left(0\right)}{u}_{00}^{2}{m}_{0}^{3}{\left(1-m\right)}^{3}+{\rho}_{1}^{\left(0\right)}{u}_{10}^{2}{m}^{3}{\left(1-{m}_{0}\right)}^{3}$ ,

${F}_{1}=-m\left(1-m\right)\left[{\rho}_{0}^{\left(0\right)}{u}_{00}^{2}{m}_{0}^{2}{\left(1-m\right)}^{2}-{\rho}_{1}^{\left(0\right)}{u}_{10}^{2}{m}^{2}{\left(1-{m}_{0}\right)}^{2}\right]$ ,

${F}_{1}={m}^{2}{\left(1-m\right)}^{2}k\left[{u}_{00}{m}_{0}\left(1-m\right)-{u}_{10}m\left(1-{m}_{0}\right)\right]$

The pressure $p\left(x\right)$ is determined by the formula

$p=\frac{{s}_{0}}{s}\left[{p}_{0}+{u}_{00}^{2}{\rho}_{00}\left(1-\frac{{m}_{0}}{m}\frac{{s}_{0}}{s}\right)+{u}_{10}^{2}{\rho}_{10}\left(1-\frac{1-{m}_{0}}{1-m}\frac{{s}_{0}}{s}\right)\right]$ (9)

The densities ${\rho}_{0}\left(x\right)$ , ${\rho}_{1}\left(x\right)$ and velocities of the components ${u}_{0}\left(x\right)$ , ${u}_{1}\left(x\right)$ are expressed in terms of porosity by the Formulas (3) and (4).

In the process of transporting the mixture on the inner wall from the side of the particles of the solid component (cotton particles of raw) acts lateral pressure $q={k}_{0}p$ ( ${k}_{0}$ —coefficient of lateral pressure), which in the direction of motion creates a frictional force $fqL$ (f—coefficient of friction between the raw and the surface of the pipe, $L=l\left(x\right)$ —the length of the contour of the cross-section of the pipe), determined, for example, by the law of dry friction Kulon. In this case, Equations (1) and (2) are written in the form:

${\rho}_{0}{u}_{0}s\frac{\text{d}{u}_{0}}{\text{d}x}=-\frac{{\rho}_{0}}{{\rho}_{0}^{\left(0\right)}}\frac{\text{d}\left(sp\right)}{\text{d}x}+sk\left({u}_{1}-{u}_{0}\right)$ (10)

${\rho}_{1}{u}_{1}s\frac{\text{d}{u}_{1}}{\text{d}x}=-\frac{{\rho}_{1}}{{\rho}_{1}^{\left(0\right)}}\frac{\text{d}\left(sp\right)}{\text{d}x}-sk\left({u}_{1}-{u}_{0}\right)-{k}_{0}lfp\left(1-m\right)$ (11)

After excluding the derivative $\frac{\text{d}\left(sp\right)}{\text{d}x}$ from system (1) and (11), we obtain

${\rho}_{0}^{\left(0\right)}{u}_{0}\frac{\text{d}{u}_{0}}{\text{d}x}-{\rho}_{1}^{\left(0\right)}{u}_{1}\frac{\text{d}{u}_{1}}{\text{d}x}=k\left({u}_{1}-{u}_{0}\right)+\frac{{\rho}_{1}^{\left(0\right)}}{s{\rho}_{1}}\left(1-m\right)l{k}_{0}fp$

After the statement of the expressions of derivatives $\frac{\text{d}{u}_{0}}{\text{d}x}$ and $\frac{\text{d}{u}_{1}}{\text{d}x}$ taking into account dependencies (4), we obtain

$\frac{\text{d}m}{\text{d}x}=Z\left(m,x,p\right)$

$Z\left(m,x,\u0440\right)=\frac{{F}_{1}\left(m\right)}{{F}_{0}\left(m\right)}\frac{{s}^{\prime}\left(x\right)}{s\left(x\right)}+\frac{{F}_{2}\left(m\right)}{{F}_{0}\left(m\right)}\frac{{s}_{0}}{s\left(x\right)}-\frac{{k}_{0}fl\left(x\right)s\left(x\right){m}^{3}{\left(1-m\right)}^{3}}{{s}_{0}^{2}}p$

To determine the derivative $\frac{\text{d}\left(sp\right)}{\text{d}x}$ we add Equations (10) and (11) and, taking into account (3), we have

${\rho}_{00}{u}_{00}\frac{\text{d}{u}_{0}}{\text{d}x}+{\rho}_{10}{u}_{10}\frac{\text{d}{u}_{1}}{\text{d}x}=-\frac{\text{d}p}{\text{d}x}-l{k}_{0}fp\left(1-m\right)/s.$

Adding the values of the derivatives $\frac{\text{d}{u}_{0}}{\text{d}x}$ and $\frac{\text{d}{u}_{1}}{\text{d}x}$ we get:

$\frac{\text{d}p}{\text{d}x}+{Z}_{0}\left(m,x\right)p={Z}_{1}\left(x,m,p\right)$ (12)

Here

${Z}_{0}=\frac{l\left(x\right){k}_{0}f\left(1-m\right)}{s\left(x\right)}$ , ${Z}_{1}\left(x,m,p\right)={F}_{3}\left(x,m,p\right)+{F}_{4}\left(x,m\right)$ ,

${F}_{3}=\frac{{s}_{0}^{2}}{{s}^{2}\left(x\right)}\frac{Z\left(x,m,p\right)}{{m}^{2}{\left(1-m\right)}^{2}}\left[{\rho}_{00}{u}_{00}^{2}{m}_{0}^{2}\left(1-m\right)-{\rho}_{10}{u}_{10}^{2}{\left(1-{m}_{0}\right)}^{2}m\right]$ ,

${F}_{4}=\frac{{s}_{0}^{2}}{{s}^{2}\left(x\right)}\frac{{s}^{\prime}\left(x\right)}{s\left(x\right)}\frac{1}{m\left(1-m\right)}\left[{\rho}_{00}{u}_{0}^{2}{m}_{0}\left(1-m\right)+{\rho}_{10}{u}_{10}^{2}m\left(1-{m}_{0}\right)\right]$

Equations (11) and (12) together describe a system that satisfies the following initial conditions for determining cotton porosity and pressure:

At $x=0$ , $m={m}_{0}$ , $p={p}_{0}={\rho}_{00}{u}_{00}^{2}/2$ .

3. Computer Processing and Analysis of the Obtained Regularities

The equations are processed by the Maple 2015 program with the following parameters for a constant-section pipeline and the corresponding graphical dependencies are obtained: ${m}_{0}=0.4$ , ${u}_{00}=20\text{\hspace{0.17em}}\text{m}/\text{s}$ , ${\rho}_{00}=1.2\text{\hspace{0.17em}}\text{\u043a\u0433}/{\text{\u043c}}^{\text{3}}$ , ${Q}_{0}=10000\text{\hspace{0.17em}}\text{\u043a\u0433}/\text{hour}$ , ${\rho}_{10}=38\text{\hspace{0.17em}}\text{\u043a\u0433}/{\text{\u043c}}^{\text{3}}$ , $s={s}_{0}=\text{\pi}{D}^{2}/4$ , $l=\text{\pi}D$ , ${k}_{0}=0.6$ . The current coefficient of porosity of cotton is equal $m=0.6$ , and the boundary values of the coefficient of interaction are respectively equal ${k}_{\u043f}=50\text{\hspace{0.17em}}\text{H\u0441}/{\text{\u043c}}^{\text{4}}$ , ${k}_{\u043f}=100\text{\hspace{0.17em}}\text{H\u0441}/{\text{\u043c}}^{\text{4}}$ , ${k}_{\u043f}=200\text{\hspace{0.17em}}\text{H\u0441}/{\text{\u043c}}^{\text{4}}$ , ${k}_{\u043f}=300\text{\hspace{0.17em}}\text{H\u0441}/{\text{\u043c}}^{\text{4}}$ [9] [10] .

Figures 1-4 present the results of an analysis of the equations for the diameter of the pipeline $D=0.4\text{\hspace{0.17em}}\text{\u043c}$ : the change in the porosity of the cotton $m\left(x\right)$ , the velocity of the components ${u}_{0}\left(x\right)$ , ${u}_{1}\left(x\right)$ , and the pressure $p\left(x\right)$ along the length of the pipeline.

To establish the influence of the friction coefficient f on the transportation process, the calculations were made for its various values. In the absence of friction, the movement of the transported material occurs at all values of the interaction coefficient large from zero (
$k>0$ ). And, in the presence of friction, the movement occurs at
$k\ge {k}_{\u043f}$ (k_{n}—boundary value of the coefficient). At smaller values of this coefficient is k. At smaller values of this coefficient (
$k<{k}_{\u043f}$ ) in some sections of the pipeline, the pressure breaks down and movement will not occur. For such cases, the adopted model does not work. For all other cases, the equations obtained correctly describe the process of moving a two-phase flow.

For example, with $f=0.05$ , $f=0.1$ , $f=0.2$ and $f=0.3$ , the interaction coefficient will have the following values: ${k}_{\u043f}=8.7\text{\hspace{0.17em}}\text{H\u0441}/{\text{\u043c}}^{\text{4}}$ , ${k}_{\u043f}=16.8\text{\hspace{0.17em}}\text{H\u0441}/{\text{\u043c}}^{\text{4}}$ , ${k}_{\u043f}=35.4\text{\hspace{0.17em}}\text{H\u0441}/{\text{\u043c}}^{\text{4}}$ , ${k}_{\u043f}=50.5\text{\hspace{0.17em}}\text{H\u0441}/{\text{\u043c}}^{\text{4}}$ .

The analysis shows that the value of porosity (Figure 1), retains the property of a regular total increase for large values k_{n} (the graph on the right) is relatively smaller (the graph on the left) and is about 0.85. At the same time, when
$k=400\text{\hspace{0.17em}}\text{H\u0441}/{\text{\u043c}}^{\text{4}}$ the increase is much more intense at the beginning of the process, and after passing through 24 - 25 meters, it remains relatively constant.

On the pressure graph (Figure 2) there is a difference with the initial parameters—the pressure at the beginning of the pipeline is about 75 Pa for smaller (graph on the left), 95 Pa (graph on the right) for larger values k_{n}. However, the pressure value at lower values k_{n} in the initial 1 - 2 meters of the pipeline decreases sharply. When reaching 20 meters the intensity drops, then the pressure is relatively stabilized. At high values k_{n}, the pressure decreases with almost the same intensity.

The graph of air velocity (Figure 3) shows the property of decline with the same intensity for all values k_{n}. When it reaches 50 m, the air velocity drops to 14 - 12 m/s, then almost does not change. And the speed of cotton (Figure 4), starting from the 0th value, gradually increases. And when it reaches 50 m, it is relatively stabilized with a value of 7 - 8 m/s.

Then, it remains almost unchanged, and the speed of cotton (Figure 4), on the contrary, is of increasing nature. If you look at the charts, the speed of the cotton in the initial 20 meters of the pipeline length changes dramatically, then more stably. At large values of the interaction coefficient, although at the beginning of the pipeline the speed changes greatly, after passing about 30 m distance

Figure 1. Changing the porosity of cotton along of the pipeline.

Figure 2. The change in air pressure along the length of the pipeline.

Figure 3. The change in air velocity along the length of the pipeline.

Figure 4. Change the speed of cotton along the length of the pipeline.

there is no significant change in speed. At large values of the interaction coefficient, the nature of the change of parameters is preserved. The difference is only the intensity of the change. Figure 5 represents the change in the porosity of cotton in both cases at a distance of 50 m reaches the same value, 0.75 - 0.80.

Only with a large k, the growth of porosity in the initial 10 - 20 meters is much stronger and when it reaches 30 m, the change is significantly stabilized. A, the change in pressure, on the contrary (Figure 6), with smaller k decreases more strongly than relatively large k. This shows that k negatively affects the resistance force, i.e., with its increase, the resistance force decreases, and with decreasing, on the contrary, it increases.

In both cases (Figure 7) and values, the air velocity shows a character of decline with almost the same intensity. Only if in the first graphs the stabilization of air velocity occurs at a value of 9 - 10 m/s, then in the following graphs, i.e., at large values, the air velocity remains without changes after the value of 15 m/s.

The speed of cotton (Figure 8) also changes according to the previous pattern: at smaller values k_{n} the intensity of the increase in speed is relatively low, but this increase is maintained for a long time, and at large values k_{n} the growth rate is much higher, but after 30 meters the speed of cotton is almost stabilized. At the same time, the speed is about 9 - 10 m/s. If you pay attention to the air speed, you can see its value—at a distance of 40 - 50 meters it is about 14 - 15 m/s, that is, the air is ahead of cotton with a relative speed equal to the difference of their speeds. The ratio of cotton velocity to air velocity will give the latency coefficient, which in our case is equal to the tabular values of k = 0.7 - 0.75 [8] . This

Figure 5. Changing the porosity of cotton along of the pipeline.

Figure 6. The change in air pressure along the length of the pipeline.

Figure 7. The change in air velocity along the length of the pipeline.

Figure 8. Change the speed of cotton along the length of the pipeline.

ratio shows how late cotton is from air. The fact of the equality of theoretical values with the table values, which are established experimentally, proves the compliance of the established theoretical pattern with the natural pattern of movement of the cotton-air mixture in the transport pipelines of pneumatic installations and will provide an opportunity to offer the established formulas for the calculation of pneumatic transport and its design [11] .

The resulting model is universal and can be used for calculations of pipelines with variable cross-section.

4. Findings

1) Research has established that during pneumatic transportation, in the initial 20 - 25 meters of the pipeline, the porosity of transported cotton greatly increases, which indicates its exposure to tensile deformation and proves that cotton is loosened, fluffed, and trash is separated from it.

2) The entry of cotton into the pneumatic conveying pipeline leads to a sharp decrease in pressure and air velocity, and the speed of cotton in the initial 15 to 20 meters of the pipeline increases strongly, then with more moderate intensity.

3) The multi-speed model of the movement of a cotton-air mixture as a heterogeneous medium in pneumotransport pipelines more correctly describes the process of moving air and material during pneumatic transportation, which makes it possible to suggest its use for calculating pneumatic transport and in its design.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

[1] | Mardanov, B. and Sarimsakov, O. (2017) Theoretical Study of the Process of Mechanical Breakdown of Cotton Seeds. Journale of Fergana Politechniall I, No. 1, 125-127. |

[2] | Abbazov, I., Sarimsakov, O., Khodjiev, M. and Mardonov, B. (2018) Waste Produced at Cotton Waste Factories. American Journal of ASCIT Communications, 5, 22-28. |

[3] | Sarimsakov, O. and Ruzmetov, M. (2016) Dynamics of the Interaction of the Worker Organs with Cotton in the Bulletproof Machine. Journal of Mechanical Problems, No. 4, 31-34. |

[4] |
Abdukarimovich, M.O., Ibragimovich, A.K. and Sharipjanovich, S.O. (2018) Designing a New Design of a Loading Cylinder for Pneumomechanical Spinning Machines. Engineering, 10, 345-356. https://doi.org/10.4236/eng.2018.106025 |

[5] | Kholmirzaev, F., Azimov, S., Abdurahimov, K. and Sarimsakov, O. (2019) Investigation of the Loss of Air Pressure in the Pipeline of the Cotton Pneumatic Conveying. Saudi Journal of Engineering and Technology Dubai, United Arab Emirates, 4, 23-27. |

[6] | Muradov, R., Burnashev, R.Z. and Sarimsakov, O.Sh. (2001) The Dynamic Task of the Interaction of Raw Cotton with the Working Bodies of Processing Machines. Journal of Mechanic Problems, No. 34, 55-58. |

[7] |
Obidov, A., Akhmedkhodjaev, Kh., Sarimsakov, O. and Holikov, Q. (2018) Investigation of the Properties of Fibrous Cotton Seeds, for Sorting on a Mesh Surface. The Journal of Engineering, 10, 572-578. https://doi.org/10.4236/eng.2018.109041 |

[8] |
Sarimsakov, O. (2016) The Possibility of Reducing Cotton Consumption in Cotton. American Journal of Science and Technology, 4, 68-72.
http://www.aascit.org/journal/et |

[9] | Sarimsakov, O. (2019) Improving the Process of Filing Cotton in Pneumatic Transport and Its Transportation by Air. “Navruz” Publishing House, Namangan, 252. |

[10] | Sarimsakov, O. and Gaybnazarov, E. (2016) About Energy Consumption in Pneumatic Conveying of Raw Cotton. American Journal of Energy and Power Engineering, 3, 26-29. |

[11] |
Sarimsakov, O., Xusanov, C. and Muradov, R. (2016) The Change in Air Pressure Along the Length of the Pipeline Installation for Pneumatic Conveying of Raw Cotton. Journal of Engineering and Technology, 3, 89-92.
http://www.aascit.org/journal/et |

Copyright © 2020 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.