卷 21, 编号 1 (2020)
- 年: 2020
- ##issue.datePublished##: 25.03.2020
- 文章: 17
- URL: https://journals.eco-vector.com/2712-8970/issue/view/7906
完整期次
Section 1. Computer Science, Computer Engineering and Management
Semi-analytical method for calculating elastic-hydrodynamic contact
摘要
A semi-analytical method for calculating elastic-hydrodynamic contact based on the partial use of Computer Aided Design / Computer Aided Engineering (CAD / CAE) packages and solutions of the integral equation of functional relationship between pressure and deformation have been described. The pressure in the lubricating layer is described by solving the modernized Reynolds equation taking into account the factors such as elastic deformation of surfaces in the contact zone, cavitation effect in the low-pressure region, and variable viscosity of the lubricant layer, which depends on thermodynamic parameters. Based on the stationary solution, a tensor damping coefficient has been obtained, with the help of which calculations of transient non-stationary modes that occur in cases of a sharp change in the external load have been further performed. A comparison of the results of modeling a plain bearing obtained by using the proposed semi-analytical method has been made and the full calculation performed using CAD / CAE programs such as ANSYS and COMSOL Multiphysics. The comparison showed good convergence of all numerical methods. At the same time, the “hybrid” method showed a number of advantages over direct calculations in CAD / CAE packages, such as: faster calculation speed, low requirements for computing resources and accounting for the cavitation effect. The described semi-analytical method allows to create digital twins of bearing units, centrifugal pumps and hydraulic supports used in satellite cooling systems and in rotary mechanisms of ground-based satellite dishes.



Development of method for increasing sensitivity in wireless optical data transmission channels in visible wavelength range
摘要
The original method for encoding binary data streams based on QPSK quadrature phase shift keying in a wireless optical communication channel in the visible range is suggested. The algorithm for analyzing signals in the receiving tract is presented. It allows to analyze the presence of two or three pulses of different colors at the input, which will signal the presence of interference or the occurrence of "illumination". In addition, the algorithm provides a possibility of dynamic compensation of external "illumination" by changing the gain of the photodetectors and adjusting the brightness of emitting LEDs. The functional scheme of the device for realization of the offered coding method in the wireless channel on the basis of optical radiation has been developed. Given that most photodiodes are sufficiently wide-band in the visible range of light waves, to increase sensitivity of each color channel and selectivity of the receiving tract it is necessary to apply optical filters for each color channel. The most effective are interference filters made of optically transparent materials with different physical characteristics. The approach for calculating optical filters has been presented.



Peculiarities of design software architecture of adaptive information processing, modeling and control systems
摘要
The article proposes an approach to developing the architecture of a service-oriented information processing system, modeling and process control. The system, which is being developed, is a tool for identifying, predicting and controlling discrete-continuous processes. Its mathematical apparatus is based on nonparametric algorithms of identification and control. The software architecture includes the following main modules: the module for processing data, modeling and forecasting output process variables and the process control module. The first module includes data preprocessing algorithms: normalization, centering and analysis of outliers and omissions. The modeling module is an algorithm for research and recovery dependencies between process variables, process identification using nonparametric estimation of the regression function from observations. The last module is an implementation of nonparametric dual control algorithms. Control devices built on the basis of these algorithms perform functions of both object control and its study.
The article discusses the application of architectural solutions based on two proven approaches in the field of software development: the composite approach and the service- oriented approach.. The main principles of composite architecture as a set of software systems with many characteristics that perform a specific task and service-oriented architecture as a modular approach to software development are described. The advantages of the applied composite service-oriented architecture over other variants of software architecture for control systems are shown, in particular, monolithic software architecture is compared with composite service-oriented architecture. This means that a researcher can use a single operation, which is a logically isolated, repeated task related to the production process of the enterprise. At the same time, it is necessary to ensure positive results when integrating with existing software products of enterprises which greatly complicates and requires the development of new components, as well as support for the "inherited" parts of the system.



Research of options of transition from unlimited to limited parallelism on the example of matrix multiplication
摘要
Today, there are many approaches to developing parallel programs. It is considered that it is more efficient to write such programs for a particular computing system. The article proposes to ignore the features of a particular computing system and outline plans for the development of a certain automated system that allows trying to improve code efficiency by developing programs with unlimited parallelism, as well as explore the possibility of developing more efficient programs using the restrictions imposed on maximum parallelism. This approach was demonstrated on the example of the analysis of various matrix multiplication algorithms. As a mathematical apparatus, the study considered various approaches to the description of algorithms to increase their implementation, including an approach based on unlimited parallelism and, also, an approach based on various restrictions on parallelism is proposed. In the course of the work, sequential and parallel methods of matrix multiplication were studied in detail, including tape and block algorithms. As a result of the study, various matrix multiplication methods (sequential, with left and right recursion, parallel methods) were studied and more effective ones were found in terms of the resources used and the restrictions imposed on parallelism. A sequential method and a cascade summation scheme were analyzed and proposed as possible ways of convolving the results of solving the problem obtained after the decomposition stage. Also, a number of programs with different levels of parallelism were developed and implemented in the functional-stream parallel programming language. In the future, such transformations can be carried out formally, relying on a knowledge base and a language that allows equivalent transformations of the original program in accordance with the axioms and algebra of transformations laid down in it, as well as replacing functions that are equivalent in results and have different levels of parallelization. These studies can be used to increase the efficiency of developed programs in terms of resource use in many branches of science, including in the field of software development for the needs of astronomy and rocket science.



The use of the inverse transformation method for time series analysis
摘要
In modern conditions of technology development, signs of systemacity are manifested to one degree or another in all areas, so the use of system analysis is an urgent task. In this case, the main factors in this situation are data processing and prediction of the state of a system. Mathematical modeling is used as a prediction method for a given subject area. A mathematical model is a universal tool for describing complex systems representing the approximate description of the class of phenomena of the external world expressed by mathematical concepts and language. The mathematical model can be represented as a set of systematic components and a random component. In this article, the object of prediction is the irregular random component of a model, which reflects the impact of numerous random factors. The origin, nature and laws of variation of the random variable are known, therefore, to simulate its behavior or predict its future value, one needs high degree of certainty to establish the form of continuous distribution function of the random variable. The empirical distribution function is calculated using the sample of random variable values. This empirical function is close to the values of the desired unknown function of distribution. The resulting empirical function is discrete, therefore it is necessary to apply piecewise linear interpolation to obtain a continuous distribution function.
The predicted random component of time series has been included in the initial regression model. In order to compare augmented and initial regression models, several values were excluded from the time series and new prediction was built. The value of the average approximation error for assessing the quality of the model is calculated. The augmented regression model proved to be more effective than the original one.



About non-parametric identification of partial-parametred discrete-continuous process
摘要
The paper considers a new class of models under conditions of incomplete information. We are talking about multidimensional discrete-continuous processes for the case where the components of the vector of output variables are stochastically dependent. The nature of this dependence is a priori unknown, but for some channels the a priori information corresponds to both nonparametric and parametric type of the initial data in the process under study. Such a situation leads to a system of nonlinear equations, some of which will be unknown, while others are known accurate to the vector of parameters.
The main purpose of the model is to determine the forecast of output variables with known input, and for implicit nonlinear equations it is only known that one or another component of the output depends on other variables that determine the state of the object.
Thus, a rather nontrivial situation arises when solving a system of implicit nonlinear equations under conditions where in one channel of a multidimensional system equations themselves are not in the usual sense, while in others they are known up to parameters. Therefore, an object model cannot be constructed using the methods of the existing identification theory as a result of a lack of a priori information. If it was possible to parameterize the system of nonlinear equations, then with a known input this system should be solved, since it is known and the parameterization stage is over. However, in this case it is still necessary to evaluate parameters. The main content of this article is the solution of the identification problem in the presence of a partially-parameterized discrete-continuous process, despite the fact that the parameterization stage cannot be overcome without additional a priori information on the process under study.
In this regard, the scheme for solving the system of nonlinear equations can be represented as a certain sequential algorithmic chain. First, on the basis of the available training sample, including all components of the input and output variables observation, a residual vector is formed. After that, an estimate of the object output with known values of the input variables is constructed based on the estimates of Nadarai-Watson. Thus, for given values of the input variables of such a process, it is proposed to carry out a procedure for evaluating the forecast of output variables in accordance with the developed algorithmic chain.
Numerous computational experiments, studying the proposed models of partially-parameterized discrete-continuous processes have shown their rather high efficiency. The article presents the results of computational experiments illustrating the effectiveness of the proposed technology for predicting values of output variables from known input variables.



On the function of time distribution of a complex computing system uptime
摘要
Any space computing complex is a complicated system. A complicated system is understood as a set of functionally related heterogeneous devices designed to perform certain functions and solve problems facing the system. One of the important characteristics of a system is its uptime. This characteristic is often considered to be a random variable. However, such a mathematical model is quite limited, since the uptime depends on many characteristics (parameters) that describe a system. Therefore, the uptime can be assumed to be a continuous random field (that is, a random function of many variables). It is this approach that is used in this work. If there are certain restrictions on the uptime of a computing system, upper estimates are found for the distributions of a random number of system failures. Therefore, the problem of estimating Gaussian field distribution in Hilbert space arises.
Two theorems that allow calculating the probability of a Gaussian vector falling into a sphere of a given radius are proved in the paper.
The paper is devoted to the reliability of a computing system. The random number of a computing system failures v(r) is a characteristic of its reliability. The v(r) distribution is the distribution of the sum of a computing system random uptime. It is impossible to write down the distribution v (r) explicitly. Therefore, one has to look for an estimate of these distributions from above. Assuming that the uptime of a computing system is the sum of many variables, the authors of the paper obtained the following results: it is shown that the problem of estimating the distributions of a random number of system failures can be considered as the problem of estimating the convergence rate in the central limit theorem in Banach spaces; if there are certain restrictions on the uptime of a computing system, upper estimates are found for the distributions of a random number of system failures. The estimates obtained can be used for further research in the theory of computing systems reliability. Knowing these upper estimates, it is possible to predict the level of average costs for computer systems restoration, as well as for the development of special mathematical and algorithmic support for analysis systems, for management, decision-making and information processing tasks.



Section 2. Aviation and Space Technology
Analysis of the ADS-B airspace monitoring system
摘要
One of the most important aspects of flight safety is awareness of AC air position (AC is the short for aircraft). The leading method of stating AC airspace location is the use of radar systems – primary, secondary, combined primary – secondary surveillance radars-though radar systems have significant drawbacks. However, at present, more advanced technologies are also in use, for example, ADS-B and multilateration. This article is focused on ADS-B broadcasting. Global coverage, low cost, great amount of obtainable information makes Automatic Dependent Surveillance – Broadcast a highly efficient system. Application of the method for AC air positioning is equally effective for helicopters, especially for those operated by special emergency services. As for the infrastructure of air navigation, the research in this sphere is focused on surveillance systems necessary for reliable control of increasing air traffic. The problem of better awareness of AC air position is still acute and has always been the object of extensive research. At present, home-manufactured civil aviation helicopters are practically never equipped with ADS-B transponders, and hardly ever use the available resources of transceiver-based surveillance systems. The objective of the analysis presented is to demonstrate the applicability of Flightradar system options, as well as implementation of ADS – B transponders for helicopter fleet. Operating surveillance systems like Flightradar may considerably increase flight safety by improving the awareness of helicopters current air position.



Features and modernization methods of thrust measurement devices for liquid rocket engine test stands
摘要
During the liquid rocket engines (LRE) testing, direct thrust measurement is carried out using thrust measurement devices. The aim of the work was, on the basis of existing data from the theory of tests and test stands devices, to highlight the design features of the thrust measurement devices and propose an option to improve the performance of this stand system. The work considers the basic circuit power schemes of thrust measurement devices by the example of power measuring systems of existing fire test stands and the features of work on preparing systems for testing. The types of calibration systems worked out in practice, their advantages and disadvantages, which constitute calibration errors, are considered. An option is proposed to modernize thrust measurement devices, in particular, through implementing of an electromechanical drive based on a planetary roller-screw mechanism as a force setting element into the calibration system. A possible general conceptual diagram of the power drive operation as a part of the calibration system of the thrust measurement devices is given. The advantages and disadvantages, the predicted effects of implementation are considered.
A more detailed analysis of this proposal may serve as an occasion for the modernization of the specific operational thrust measurement devices design at the fire test stand for LREs or may be a working option when designing a new thrust measurement device.



Elaboration and testing of the algorithm which ensures an achievement of minimal deviation angle of flying model’s main centroidal axis of inertia during her counterbalancing in a sole correction flatness
摘要
High complexity and cost of developing flying models necessitate the use of such design and production techniques that would ensure the best flight technical and technological characteristics of the model also would raise of it operation effectiveness. These techniques include the experimental control method of flying model’s mass-inertia asymmetry parameters during final assembly of the model. Solution of the problem of optimization the process of bringing parameters of mass-inertia asymmetry of the conical flying model to specified standards is considered in the article. The only correction plane is designed to be positioned close to cone face, away from the center mass of the flying model. The flying model as a component of prefabricated rotor is being balanced in dynamic mode on a low-frequency dynamic vertical stand, which based on gas bearings. Before balancing experiment the weigh, longitudinal center of mass and inertia moments of the flying model have to be controlled with use of another measurement equipment. As a criterion of optimization is sorted the reaching of minimum of the angle of deviation of principal longitudinal centroidal axis of inertia from geometrical axis of the flying model. But simultaneously the pre-set standard of center-mass shift from the geometrical axis must be ensured. Balancing algorithm, easy-to-realized by modern computers, is presented. Numerical illustration of balancing is given. The algorithm enables omitting intermediate steps of balancing, reducing them to one step (as a rule), and shortening the balancing time, as well. In one step of balancing the engineering model permits either bringing parameters of mass-inertia asymmetry of the flying model to specified standards, or diagnosing impossibility of attaining the specified standards with available design of flying model. The algorithm and balancing method are experimentally tested at newly-designed vertical dynamic stand on conical gas bearings. It’s high precision and efficiency are corroborated.



To the question of forecasting the technical condition of low-thrust liquid rocket engines
摘要
In the rapidly developing space and rocket industry, spacecrafts are being equipped with low-thrust liquid rocket engines. High requirements are imposed on the reliability, efficiency and economy of fuel use for this type of rocket engine. To ensure monitoring of the characteristics of spacecrafts, a functional diagnostic system is used, which includes telemetry and analytical data processing. Telemetry performs the functions of receiving and transmitting information. Information processing is carried out in computer centers located on the spacecraft and the Earth. The most promising computing tool capable of predicting time series and classifying a large amount of interconnected data is considered an artificial neural network. In this regard, the subject of research in the work is data processing methods based on an artificial neural network. The purpose of the work is to develop a method for forecasting the technical condition of low-thrust liquid rocket engines using an artificial neural network.
The relevance of research on the use of a neural network in the system of functional diagnostics of low-thrust liquid rocket engines for spacecraft is explained in the introduction. In the main part, an analysis of many telemetric data of the rocket engine is carried out and their strength in the forecast of the main diagnostic parameters is determined. It is proposed to use traction, specific impulse, and temperature of the structure as diagnostic parameters. The prognostic capabilities of the neural network were investigated and a schematic diagram of a method for predicting the technical condition of a low-thrust liquid rocket engine was developed. In the developed method, at the first stage, the neural network performs the approximation of the function and extrapolates the time series of telemetric data; the second stage determines the probable class of the technical condition of the engine.
The conclusion outlines a plan for further experimental research in the study area and provides recommendations on the development and improvement of algorithms for functioning of artificial neural networks as part of the functional diagnostics system of the spacecraft. Due to the generalized nature of the methodological schemes, the results of the work can be applied to any type of rocket engines and used at all enterprises of the rocket and space industry of the corresponding profile.



Ensuring extreme regulation of power of primary energy sources at their joint operation for total load
摘要
Heterogeneous energy sources and homogeneous energy sources with different characteristics are frequently used in autonomous power supply systems. Solar batteries are widely used as primary energy sources for on-board power supply systems of spacecrafts, unmanned and manned aircrafts. Renewable energy sources such as solar, wind, geothermal and hydro energy, serve as primary energy sources of terrestrial autonomous power supply systems.
Matching primary energy sources with different characteristics and operating conditions within a unified power supply system leads to problems connected with the power control of energy sources, which determines the relevance of the considered problems.
The main aim of the study is to develop a combination of primary energy sources and control techniques which allow using primary energy sources with different characteristics and operating conditions in unified autonomous power supply system.
The objectives of the study are to create the simulation model of a power supply system using MATLAB/Simulink software; to develop and test control algorithms for primary energy source controllers that would allow to maintain the needed battery charging current; to develop and test control algorithms for primary energy source controllers that would allow the primary energy sources to operate in the maximum power point tracking mode and to minimize the maximum power point search time.
Methods used in the study: the simulation of a power supply system using MATLAB 7.9 Simulink software.
Results: the simulation model of a power supply system including two primary energy sources with different characteristics is designed. In the case of excess power generating by the primary energy source, its controller operates in the battery charging mode. When the primary source power shortage occurs, its controller operates in the maximum power point tracking mode. The proposed power supply system structure allows controlling two energy sources independently, thus the primary energy source controllers can operate in different modes. This provides flexibility of the power supply system. The use of fuzzy logic control algorithm increases the accuracy and search speed of the maximum power point tracking algorithm. Simulation results confirmed the efficiency of the proposed solar controller operation algorithms in all modes stated above. The efficiency of controller operation modes selection algorithm was confirmed in different operating conditions. The proposed algorithms allow implementing the effective control of primary power sources depending on power supply system operating conditions.



Elaboration and testing of the algorithm which ensures an achievement of minimal deviation angle of flying model’s main centroidal axis of inertia during her counterbalancing in a sole correction flatness
摘要
High complexity and cost of developing flying models necessitate the use of such design and production techniques that would ensure the best flight technical and technological characteristics of the model also would raise of it operation effectiveness. These techniques include the experimental control method of flying model’s mass-inertia asymmetry parameters during final assembly of the model. Solution of the problem of optimization the process of bringing parameters of mass-inertia asymmetry of the conical flying model to specified standards is considered in the article. The only correction plane is designed to be positioned close to cone face, away from the center mass of the flying model. The flying model as a component of prefabricated rotor is being balanced in dynamic mode on a low-frequency dynamic vertical stand, which based on gas bearings. Before balancing experiment the weigh, longitudinal center of mass and inertia moments of the flying model have to be controlled with use of another measurement equipment. As a criterion of optimization is sorted the reaching of minimum of the angle of deviation of principal longitudinal centroidal axis of inertia from geometrical axis of the flying model. But simultaneously the pre-set standard of center-mass shift from the geometrical axis must be ensured. Balancing algorithm, easy-to-realized by modern computers, is presented. Numerical illustration of balancing is given. The algorithm enables omitting intermediate steps of balancing, reducing them to one step (as a rule), and shortening the balancing time, as well. In one step of balancing the engineering model permits either bringing parameters of mass-inertia asymmetry of the flying model to specified standards, or diagnosing impossibility of attaining the specified standards with available design of flying model. The algorithm and balancing method are experimentally tested at newly-designed vertical dynamic stand on conical gas bearings. It’s high precision and efficiency are corroborated.



Section 3. Technological Processes and Materials
Research on electrical properties of manganese sulphides doped by thulium and ytterbium ions
摘要
Materials exhibiting connection between electrical and magnetic properties are attractive for possible use as an element base in microelectronics, spintronics, and sensor devices. Compounds with mixed valence exhibit a number of metal-insulator phase transitions, magnetic phase transitions, including changes in magnetic properties without changing magnetic symmetry.
Promising materials for studying these effects are cation-substituted Mn1−xRexS compounds (Re = 4f elements) synthesized on the basis of the antiferromagnetic semiconductor of manganese monosulfide. The latter is of practical importance in the development of new materials for temperature sensors, widely used in the metallurgical industry.
The structural and electrical properties of compounds with mixed valences TmXMn1-ХS (0 ≤ X ≤ 0.15) and TmXMn1-ХS (0 ≤ X ≤ 0.25) were studied in the temperature range 80–1100K. The regions of existence of solid solutions of TmXMn1-XS sulfides with an fcc (face-centered cubic) lattice of the NaCl type were determined. It was found that conductivity decreases upon the substitution of manganese cations with thulium ions and the lattice constant increases more sharply in comparison with Vegard’s law. When ytterbium ions are substituted, the conductivity increases with increasing concentration and the temperature dependence has the form typical of semiconductors.



The influence of prefinishing operations at titanium alloys on the characteristics of mao coatings
摘要
Improving the reliability, service life and operational safety of titanium alloy structures exposed to thermal, chemical and mechanical stresses can be achieved by applying various protective coatings.
One of the effective methods of protecting such alloys is the formation on their surface of oxide coatings that are resistant to external factors. Of great interest from this point of view is the method of micro-arc oxidation (MAO), which allows one to obtain multifunctional ceramic-like oxide coatings with unique properties. Such coatings can be used to create a durable heat and electrical insulating layer on parts, protect surfaces from erosion in high-speed gas flows, corrosion in aggressive environments and wear by friction, to increase the surface emissivity, etc.
This method is well established for the oxidation of aluminum alloys. Despite the fact that the mechanism of coating formation during MAO is the same for aluminum and titanium alloys, there are certain differences in the structure and characteristics of the resulting coating. For example, it is believed that during the MAO treatment of aluminum alloys, preliminary surface preparation is not required and the adhesive strength is comparable with the strength of the substrate material. However, when processing titanium alloys, we noted cases of a significant decrease in adhesive strength. One of the reasons may be the lack of preliminary surface preparation before coating.
Therefore, studies aimed at studying the influence of the method of surface preparation and the resulting roughness on the characteristics of the applied coatings are relevant.



Peculiar properties of technological improvement and optimization of production costs of 3D-configuration pipes
摘要
The analyzes of the requirements to 3D-configuration pipelines production at the rocket and space industry enterprises is done. A review of different approaches to pipe bending technology (with heat treatment and without heat treatment) is carried out. The object of the study is the bending process and a universal bending machine for pipelines’ production of complex configuration. The article is divided into four sections, which consider the key factors, causing directly the effectiveness of the technological operation of pipeline bending of a complex 3D-trajectory. An overview of no-temperature shaping of the pipeline is given in the first section. The requirements to the technology, excluding: corrugation, flattening, stretching and thinning of pipeline walls during their bending, are considered. The actual regulatory documents and industry aerospace standards, regulating production of pneumatic and hydraulic pipelines are given. An example of calculating the minimal allowable bend radius of the pipe, depending on the diameter and thickness of the pipe wall, is given. The requirements to unification of the pipe size production and gaps are listed. The dependence of the maximal allowable internal pressure in the pipeline is shown. The requirements to equipment, used in pipeline bending and to the design of the pipe bending machine are considered. In the second section, the possibilities of temperature influence on the pipe bending process are viewed. The analysis of patent and technical literature and six possible methods of effective thermal effects are presented: heating of the whole pipeline length, narrow zone heating of the bend pipe place, water cooling with nitrogen in the pipe, laser-cooling of atoms of the pipes, application of the petroleum products on the place of heating of the pipe and using of modern fillers inside the pipe to change its temperature. In the third section the tasks of the development of a universal bending machine are set; the system of the algorithm of the universal bending machine operation is considered; the system of algorithm of the bending machine operating with CNC is shown. The General functional scheme of the bending machine and the sequence diagram of the equipment operation is given.



Influence of a constant electric field on the adsorption purification of water from iron ions
摘要
Using electrochemical action (ECA) to treat water was first proposed in UK in 1889. At present, many methods of ECA are known (electro flotation, electro coagulation, electro osmosis, electrophoresis, etc.).
In the production of rocket and space technology, galvanic technologies are used, as a result of which waste water is contaminated with metal ions. Known methods of wastewater treatment do not allow to ensure the maximum permissible concentration of metal ions in treated water, or are expensive or difficult to operate in industry. Iron ions are among the most polluting components of wastewater of most industries. So increased control and the development of effective methods of wastewater treatment are necessary. Iron affects the intensity of phytoplankton development and the qualitative composition of microflora in reservoirs. The toxicity of iron compounds in water depends on the hydrogen index of water. The alkaline environment dramatically increases the risk of fish poisoning, as in such conditions, iron hydroxides are formed, which are deposited on the gills, clog and corrode them. In addition, iron compounds bind oxygen dissolved in water, which leads to the mass death of fish and other hydrobionts.
The article presents the method of conducting experiments, the methods of sorption, electrochemical and combined water treatment, including electrochemical action and adsorption. The results of studies of these methods of water purification from iron ions are presented. The dependence of the degree of purification on the electric field strength, interelectrode distance and water treatment time is revealed. With an electric field strength of 5.16 V/mm, a temperature of 20–22 °C using quartz sand as an adsorbent and a processing time of 1 minute, the concentration of iron ions decreased from 2.5 to 0.25 mg/l (at MPC = 0.3 mg/l). The proposed combined cleaning method requires inexpensive and affordable materials and is easy to operate.


