Research article Special Issues

A computational study of time-fractional gas dynamics models by means of conformable finite difference method

  • Received: 21 February 2024 Revised: 19 May 2024 Accepted: 27 May 2024 Published: 19 June 2024
  • MSC : 26A33, 35R11, 76M20, 76N10

  • This paper introduces a novel numerical scheme, the conformable finite difference method (CFDM), for solving time-fractional gas dynamics equations. The method was developed by integrating the finite difference method with conformable derivatives, offering a unique approach to tackle the challenges posed by time-fractional gas dynamics models. The study explores the significance of such equations in capturing physical phenomena like explosions, detonation, condensation in a moving flow, and combustion. The numerical stability of the proposed scheme is rigorously investigated, revealing its conditional stability under certain constraints. A comparative analysis is conducted by benchmarking the CFDM against existing methodologies, including the quadratic B-spline Galerkin and the trigonometric B-spline functions methods. The comparisons are performed using L2 and L norms to assess the accuracy and efficiency of the proposed method. To demonstrate the effectiveness of the CFDM, several illustrative examples are solved, and the results are presented graphically. Through these examples, the paper showcases the capability of the proposed methodology to accurately capture the behavior of time-fractional gas dynamics equations. The findings underscore the versatility and computational efficiency of the CFDM in addressing complex phenomena. In conclusion, the study affirms that the conformable finite difference method is well-suited for solving differential equations with time-fractional derivatives arising in the physical model.

    Citation: Majeed A. Yousif, Juan L. G. Guirao, Pshtiwan Othman Mohammed, Nejmeddine Chorfi, Dumitru Baleanu. A computational study of time-fractional gas dynamics models by means of conformable finite difference method[J]. AIMS Mathematics, 2024, 9(7): 19843-19858. doi: 10.3934/math.2024969

    Related Papers:

    [1] Gang Chen, Jinbo Ni, Xinyu Fu . Existence, and Ulam's types stability of higher-order fractional Langevin equations on a star graph. AIMS Mathematics, 2024, 9(5): 11877-11909. doi: 10.3934/math.2024581
    [2] Jiali Yu, Yadong Shang, Huafei Di . Existence and nonexistence of global solutions to the Cauchy problem of thenonlinear hyperbolic equation with damping term. AIMS Mathematics, 2018, 3(2): 322-342. doi: 10.3934/Math.2018.2.322
    [3] Abdissalam Sarsenbi, Abdizhahan Sarsenbi . Boundary value problems for a second-order differential equation with involution in the second derivative and their solvability. AIMS Mathematics, 2023, 8(11): 26275-26289. doi: 10.3934/math.20231340
    [4] Hawsar Ali Hama Rashid, Mudhafar Fattah Hama . Approximate solutions for a class of nonlinear Volterra-Fredholm integro-differential equations under Dirichlet boundary conditions. AIMS Mathematics, 2023, 8(1): 463-483. doi: 10.3934/math.2023022
    [5] Teeranush Suebcharoen, Watchareepan Atiponrat, Khuanchanok Chaichana . Fixed point theorems via auxiliary functions with applications to two-term fractional differential equations with nonlocal boundary conditions. AIMS Mathematics, 2023, 8(3): 7394-7418. doi: 10.3934/math.2023372
    [6] Mian Bahadur Zada, Muhammad Sarwar, Reny George, Zoran D. Mitrović . Darbo-Type Zm and Lm contractions and its applications to Caputo fractional integro-differential equations. AIMS Mathematics, 2021, 6(6): 6340-6355. doi: 10.3934/math.2021372
    [7] Hanadi Zahed, Ahmed Al-Rawashdeh, Jamshaid Ahmad . Common fixed point results in F-metric spaces with application to nonlinear neutral differential equation. AIMS Mathematics, 2023, 8(2): 4786-4805. doi: 10.3934/math.2023237
    [8] Muath Awadalla, Mohammed S. Abdo, Hanan A. Wahash, Kinda Abuasbeh . Qualitative study of linear and nonlinear relaxation equations with ψ-Riemann-Liouville fractional derivatives. AIMS Mathematics, 2022, 7(11): 20275-20291. doi: 10.3934/math.20221110
    [9] Apassara Suechoei, Parinya Sa Ngiamsunthorn . Extremal solutions of φCaputo fractional evolution equations involving integral kernels. AIMS Mathematics, 2021, 6(5): 4734-4757. doi: 10.3934/math.2021278
    [10] Hongzhuan Wang, Xianhao Shi, Ber-Lin Yu . On the eccentric connectivity coindex in graphs. AIMS Mathematics, 2022, 7(1): 651-666. doi: 10.3934/math.2022041
  • This paper introduces a novel numerical scheme, the conformable finite difference method (CFDM), for solving time-fractional gas dynamics equations. The method was developed by integrating the finite difference method with conformable derivatives, offering a unique approach to tackle the challenges posed by time-fractional gas dynamics models. The study explores the significance of such equations in capturing physical phenomena like explosions, detonation, condensation in a moving flow, and combustion. The numerical stability of the proposed scheme is rigorously investigated, revealing its conditional stability under certain constraints. A comparative analysis is conducted by benchmarking the CFDM against existing methodologies, including the quadratic B-spline Galerkin and the trigonometric B-spline functions methods. The comparisons are performed using L2 and L norms to assess the accuracy and efficiency of the proposed method. To demonstrate the effectiveness of the CFDM, several illustrative examples are solved, and the results are presented graphically. Through these examples, the paper showcases the capability of the proposed methodology to accurately capture the behavior of time-fractional gas dynamics equations. The findings underscore the versatility and computational efficiency of the CFDM in addressing complex phenomena. In conclusion, the study affirms that the conformable finite difference method is well-suited for solving differential equations with time-fractional derivatives arising in the physical model.



    In recent decades, the rapid development of industry leads to high activities in commercial and agriculture causing an increase in the urban population. Consequently, the high utilization of natural resources leads to the emission of pollutants grew rapidly affect the air quality in the surrounding areas [1,2,3,4]. Almost all air pollutants such as nitrogen oxides (NOx), sulfur oxides (SOx), polycyclic aromatic hydrocarbons (PAHs), carbon monoxide (CO), volatile organic compounds (VOCs), and heavy metals come from road transport, which is generated by the combustion engines of the vehicles and other components such as brake and clutch linings and pads, tires, and fuel tanks. Moreover, particulate matter (PM) has been detected which is also formed in the atmosphere and is a major risk to human health [5]. Particulate matter has been classified by particle size as coarse (2.5–10 μm), fine (0.1–2.5 μm), and ultrafine (<0.1 μm), also called nanoparticles because of their size restrict to the 100 nm or smaller particles. Their toxicity depends on the size of the particles because they are retained longer in the human body [6,7,8].

    In recent years, many publications of the study on air pollution control have been published and various technologies have addressed the removal or reduction of toxic pollutants. There are different methods to treat the air pollutants such as the use of electrokinetic treatment [9], bioelectrochemical technology in the form of air-cathode single chamber microbial fuel cells (SC-MFCs) [10], nano-bioremediation technologies [11,12], bio-nanomaterials (nanofibrous air filters, nano-adsorptive materials) [13,14]. Nanotechnology's air pollution control system has been proposed in three categories: remediation and treatment, detection and sensing, and pollution prevention [15]. Even though the nanomaterials presented a high performance for air pollutant treatment, the synthesis of NMs is quite expensive and involves hazardous chemical substances. Therefore, phytoremediation technology has been established, the specific plants are employed to decontaminate the pollutants. The distinct plants accumulate pollutants in their parts and mobilize them into the tissues [16]. Phytoremediation is driven by solar energy, a natural method, and an inexpensive process; however, it is a long-term process.

    Traditionally, air pollution assessment is based on monitoring of pollutants in the surrounding area and the analysis measurement related to chemical substances' effects on both the environment and organisms. The direct measurements of particulate matter have several techniques as measuring concentration (gravimetric, optical, and microbalance), size distribution Scanning Mobility Particle Sizer (SMPS), and Electrical Low-Pressure Impactor (ELPI) [17,18,19]. Therefore, plants are also used to be biomonitoring of air pollutants and are suitable tools for evaluating the levels of air pollutants, having several advantages compared with the conventional methods [20,21]. Tree bark (e.g., black poplar tree, pine tree) is also used to examine the levels of toxic metals and trap the particulate matter in the atmosphere because their porous surface can enable efficient accumulation and retention of heavy metals [22,23,24]. The leaves of an ornamental tree (Tibouchina granulosa (Desr.) Cong.) with various types of trichomes located on the adaxial and abaxial surfaces have been investigated to adsorb particulate matter (PM) revealed that the adsorbed particles size range 2.5–100 μm and did not return to the environment under normal weather conditions [25].

    Furthermore, lower plants like mosses and lichens can be used as atmospheric pollutant monitors. Lichens are used to be biomonitors of air pollutants such as particulate matter, Persistent Organic Pollutants (POPs), SO2, fluorides, ozone, nitrogen oxides, peroxyacetyl nitrate (PAN), and heavy metals due to their extraordinary capability to accumulate mineral elements and their occurrence at a large geographical range [26,27,28,29]. Mosses are effective in adsorbing and accumulating particulate matter (e.g., benxohexachloride isomers and polyaromatic hydrocarbons) in the air [30]. Tillandsia usneoides L. (Spanish moss) has been reported as a biomonitor assessing the particulate matter and examining particulate matter deposited on moss samples by SEM-EDS, ICP-MS, atomic absorption spectrometry [31], neutron activation analysis (NAA), etc. The Spanish moss grows in a zigzagging pattern, linear, twisted, and densely with scaled leaves, in addition, their trichomes include a stalk and flexible wings that can absorb water and nutrients directly from the air instead of their root, it does only adhere to a substratum. Consequently, Spanish moss is more effective to inhabit extreme conditions and is suitable for monitoring air pollution [32,33,34]. This method is inexpensive, easy to use, and is an alternative way to be developed on-site where there is no infrastructure or instruments for air quality monitoring. However, the adsorption and accumulation of the particulate matter probably depend on the morphology of the species due to the spaces between leaves and stems, so the focus should be on selecting plants for biomonitoring technology.

    Poor air quality in Thailand is attributed to PM2.5 which occurs from burned crop fields to harvest sugarcane or remove remnant biomass in the burning season, moreover, industries and vehicles are generated air pollutants [35]. Nevertheless, the research on mosses as biomonitoring is still low in numbers. Native thalli of the epiphytic lichen Parmotrema tinctorum were used as an effective tool for assessing metal pollutants in a Khao Yai Natural Park [36]. Thai native mosses have been reported as a bio-adsorbent to assess the heavy metal in the atmosphere [37]. However, the moss-containing building model used in particulate matter treatment had never been studied, just only the design of the air treatment plant provides only architectural aesthetics.

    This study focuses on the use of Riccia, a genus of liverworts in the family Ricciaceae of division Marchantiophyta as biofilters in the adsorption chamber to treat the particulate matter. The alignment of plants used to treat air pollutants using flow dynamic mathematical models was performed using three different arrangements including blocking panels, zigzags panels, and parallel panels. The zigzag panels model revealed maximum absorption efficiency and were suitable for plant filters due to the high surface area to adsorb the air pollutants [38]. In this work, we built the absorption chamber to investigate the efficiency of Riccia for treating the particulate matter generated by a diesel smoke generator. Moreover, the second aim of this work was also to optimize the alignment of Riccia in the air chamber. The sample was identified based on morphological descriptions and taxonomic keys. The accumulated PM deposited on the thalli of Riccia was evaluated through SEM observations. We assumed that the development of bio-innovative research could help to solve and reduce urban air pollution.

    Glutaraldehyde, sodium cacodylate, and osmium tetroxide were used as received. Distilled water was used for plant sample preparation and washing purposes. The Riccia was obtained from the Moss Land shop, which is the online shop of the greenhouse garden in Nonthaburi Province in central Thailand. The species was identified by Assoc. Prof. Sahut Chantanaorrapint, Prince of Songkla University herbarium (PSU), Division of Biological Science, Faculty of Science, Prince of Songkla University, specializes in Bryophytes in Thailand. The ash content and moisture content were evaluated. The Air pollution tolerance index (APTI) which included four parameters; total chlorophyll, relative water content (RWC), Ascorbic acid content, and pH of extract plant was also measured. All samples were cut into similar pieces and had a dimension of 5 × 5 cm2. The study specimens were washed with deionized water and dried in the open air at room temperature. The absorption efficiency was performed in the treatment chamber with simulated diesel smoke. The concentration of particulate matter was measured with various airflow speeds of 0.5, 1.0, 1.5, and 2.0 m/s using the Laser dust sensor PMS3003 Arduino. Data were presented as mean values with standard deviation (SD). All graphs were produced with Origin 2021. In addition, the volatile organic compounds (VOCs) and polyaromatic hydrocarbons (PAHs) of the treated plant were analysed using Headspace gas chromatography/ mass spectroscopy (HS-GC-MS). The metals were analysed using Flame Atomic absorption spectroscopy (FAAS). A scanning electron microscopy (SEM) image was collected by the Hitachi S2500 scanning electron microscope to analyze the morphology of liverwort samples.

    In this research, the selected plant act as a biofilter is Riccia due to its morphology and ecology. Riccia is a genus of liverworts belonging to the Ricciaceae family of the Marchantiophyta division. It is complex thalloid liverworts that its thalli are internally differentiated. Its dorsal side tissue is green, while the ventral tissue side is colorless. The rhizoid walls are mostly densely papillose, with numerous papilla-like projections on the inner walls, and the ventral surface of the thallus is usually covered by multicellular scales, in two or more rows, which are often deeply purple. The reproduction can be either sexual or asexual. The genus Riccia is terrestrial, mostly growing on moist damp soil, rock, the bark of trees, and rotten wood, and adsorbs nutrients and water from the air through its outer surfaces instead of using roots, as many plants do [39]. It prefers moist habitats such as tropical montane cloud forests, moist damp soil, or shady places, and is used to be a good indicator of humidity because it can absorb large quantities of rainwater caused by thick layers of the thallus. In tropical cloud forests, the rainfall is captured by the liverworts up to 20%–40% [40]. It plays an important role in the water balance and nutrient cycles of the forest. Therefore, it has the possibility of the use Riccia as a biofilter. Recently, particular attention is focus on absorbing particulate matter from air pollution by living biomass such as mosses [41,42], and lichens [43].

    The ash content is used to determine the total composition of biomass samples, which consists of the mineral content and other inorganic matter in biomass. The sample was weighed 0.5 g, to the nearest 0.1 mg of a test specimen into the tared crucible which was prepared by placing them in the furnace at 600 ℃ for a minimum of three hours and allowed to cool at room temperature in a desiccator. Place the sample into the furnace at 600 ℃ for a minimum of three hours and dry to constant weight. Constant weight is defined as less than ± 0.3 mg change in the weight upon one hour of re-heating the crucible. The ash content was calculated by using Eq (1) as follows:

    Ash content (%) = WCA-WCWix100,  (1)

    where WCA, WC, and Wi are the weight of the crucible plus ash, the weight of the crucible, and the dry weight of the sample, respectively.

    The measurement of moisture content by direct measurement is determined by removing moisture in an oven and then by measuring weight loss. The sample was weighed 0.5 g, to the nearest 0.1 mg of a test specimen into the beaker which was prepared by placing them in the oven at 105 ℃ for a minimum of six hours and allowed to cool in a desiccator. Place the sample into the oven at 105 ℃ for a minimum of six hours and dry to constant weight. Constant weight is defined as less than ± 0.3 mg change in the weight upon one hour of re-heating the beaker. The moisture content was calculated by using Eq (2) as follows:

    Moisture content (%) = Wi- WtWix 100, (2)

    where Wi and Wt are the initial weight of the plant, and the terminated weight of the plant, respectively.

    1.0 g of the study plant (fresh weight; FW) was floated in distilled water in darkness at 4 ℃ for 24 h. Then, dried with filter paper and its turgidity weight (TW) was determined. Place it into the oven at 70 ℃ for 24 h to determine its dry weight [44]. RWC was calculated by using Eq (3) as follows:

    RWC = [FW-DW][TW-DW]x 100. (3)

    1.0 g of the study plant was crushed by porcelain mortar and homogenized in 10 mL of deionized distilled water. The sample was centrifuged at 2,400 rpm for 10 min and the pH was measured using a pH meter.

    1.0 g of the study plant was crushed and 10 mL of 80% acetone was added. Allow to room temperature for 15 min for thorough extraction. The extract was centrifuged at 2,400 rpm for 10 min. Then, the absorbance was measured at 645 nm and 663 nm using a spectrophotometer after calibrating with the reagent blank. The chlorophyll A, chlorophyll B, and total chlorophyll content were calculated using Eqs (4)–(6) as follows [45]:

    Chlorophyll A = [12.72 (A663) - 2.69 (A645)] x[V1,000W], (4)
    Chlorophyll B = [22.9 (A645) - 4.68 (A663)] x[V1,000W], (5)
    Total chlorophyll = [20.2 (A645) + 8.02 (A663)] x[V1,000W], (6)

    where A645, A663, V, and W are the absorbance at 645 nm, absorbance at 663 nm, the total volume of the extract (mL), and the weight of the sample (g), respectively.

    The ascorbic acid content was measured by the method of titration with 2, 6-dichlorophenol indophenol (DCIP). 1.0 g of the study plant was homogenized by grinding in a pestle with a mortar and 5% oxalic acid was added to avoid such a loss of ascorbic acid. The extract was filtered, and the volume is made up to 100 mL with 5% oxalic acid. Then, 10 mL of the unknown solution was titrated with the DCIP dye, and the appearance of pink color indicated the endpoint. 10 mL of ascorbic acid solution and 5% oxalic acid were titrated with DCIP dye for standard and blank solutions. The ascorbic acid content was calculated by using Eq (7) as follows:

    Ascorbic acid = [U-BS-B]xA total volume of unknownThe volume has takenfor each samplex Amount of AA in 10 mL of standard, (7)

    where U, S, and B are the titre value of unknown, titre value of standard, and titre value of blank, respectively.

    APTI of the plant was measured by four parameters of RWC, pH, TCH, and AA using Eq (8) as follows:

    APTI = [A(T+P) + (R)]10, (8)

    where A, T, P, and R are ascorbic acid content (mg/g of dry weight), total chlorophyll content (mg/g), plant extract pH, and relative water content (%), respectively.

    APTI values were categorized into four tolerance levels, there are level (1) tolerant (T) species, level (2) moderately tolerant (MT) species, level (3) intermediate (I) species, and level (4) sensitive (S) species [45].

    Three replicates of plant specimens were cut into similar pieces with the dimension of 1 × 1 cm and transferred to the zip bag which was weight before use. The dust was generated from talcum powder and the particle size was measured using the Laser dust sensor PMS3003 Arduino. The sizes ranged from 1.0–30.0 μm. The dust was added to the zip bag and the weight of the plant with dust and zip lock bag was measured. Then, the dust was removed and the weight of the plant after removing of dust was determined. The dust-capturing potential of the plant was calculated using Eq (9) as follows:

    W = (W1-W2-W3)A, (9)

    where W, W1, W2, W3, and A are dust-capturing potential (g/cm2), the weight of plant with dust and zip lock bag (g), the weight of plant after removing dust (g), weight of zip lock bag (g), and total area of the plant (cm2).

    The particulate matter used in the experiment was generated by simulated fuel combustion smoke using a diesel smoke generator and loaded into the tank with a capacity of 200 liters by a fan. The generated particulate matter was distributed by the size ranges (in micrometers) as follows: PM1, PM2.5, and PM10. The concentration was set at 5,000 ± 250 μg/m3 using automatic adjustment of the smoke pipe valve from the controlled program. The smoke was released into the adsorption chamber through a 1-meter-long pipe via the PM2.5 concentration sensor (Laser dust sensor PMS3003 Arduino) with the flow speed of air pollutants at the inlet of the box of 0.5, 1.0, 1.5, and 2.0 m/s.

    All plant samples were prepared in the open-air room and cut into a similar pieces with an area of 5 × 5 cm2. The adsorption efficiency of the PM adsorbed on the biofilter with the various wind speed conditions was measured. The inflow air velocity was set to four different velocities of 0.5, 1.0, 1.5, and 2.0 m/s for one hour each. The initial concentration of PM2.5 and the remaining PM2.5 were measured before and after entering the treatment box, respectively. Four panels set up the chamber with the different alignments of sample plants, vertical, zigzag, and horizontal. Each panel has three points in a vertical line and four points in a horizontal line and the distance from each other is 10 centimeters, therefore, there are twelve points of each panel in the treatment box as described by Meesang et al. [38]. The control system measured the concentration of PM2.5 before entering and after entering through the adsorbed plants in the treatment box. All data were selected using IBM SPSS statistics 25 followed by Microsoft excel for calculating the adsorption efficiency value. The scattering graphs were plotted using Origin 2021. The adsorption efficiency was expressed as a percentage of the residue matter compared to its initial concentration by using Eq (9) as follows [46]:

    η(%) = [PM2.5IN-PM2.5OUTPM2.5IN]x 100, (9)

    where PM2.5IN and PM2.5OUT are the initial and remaining concentrations of a particulate matter 2.5, respectively.

    The indirect determination of volatile organic compounds (VOCs) and polyaromatic hydrocarbons (PAHs) constituents in plant samples by analysing the vapor phase that is in thermodynamic equilibrium with the sample in a closed system using headspace gas chromatography/ mass spectroscopy (HS-GC-MS). The extraction was performed in an HS vial filled with plant samples and sealed with a silicone septum cap, thermostated for 15 min at 105 ℃. Once equilibrium between the matrix and the gaseous phase was reached, the HS needle descended, pierced the septum cap, and pressurized with nitrogen at 45 psi for 3 min. Then, the over-pressured gas sample in the HS vial was applied to the GC column at 120 ℃ with a flow rate of 1 mL/min. Subsequently, the GC column was used to separate the VOCs and PAHs. The mass spectrometer was used as the detector. The data was carried out for 14 min at a scan rate of 1 scan s−1. Peak identification in the HS-GC-MS chromatograms was done by comparing their mass spectra with those contained in the NIST (National Institute of Standards and Technology) library.

    For analysis of the samples with atomic absorption spectrophotometry, plant samples of around 0.5 g were digested in a solution of HNO3:H2O2 (7:2) in Teflon digestion vessels and exposed to microwaves (METASH, MWD-610). The digestion was carried out in 2 steps, in the first step, the solution was heated at 180 ℃ for 20 min ramp time, and the next step was carried out at 180 ℃ for 20 min hold time with a power of 800 W. Finally, the vessels were allowed to cool and carefully opened, the digest solution transferred into 50 mL of volumetric flask and volume adjusted by added the distilled water. Metal concentrations were measured by atomic absorption spectrophotometry (Perkin-Elmer, PinAAcle 900F). Samples were analyzed in triplicates and the results were given in mg/g dry weight.

    All plant samples were cut into pieces, approximately 4 mm long, and immediately fixed in glutaraldehyde (2.5%) in sodium cacodylate buffer (0.1 M, pH 7.2–7.4) at 4 ℃ for 4 hours, rinsed three times in sodium cacodylate / Phosphate buffer (pH 7.2–7.4) at 4 ℃ for 15 minutes each. Post-fixed sample with 1.0% Osmium tetroxide in buffer at 4 ℃ for 1 hour. Dehydrated tissue with a series of ethanol from 30% v/v-100% v/v ethanol at 4 ℃, 2 times for 15 minutes each. Dried by the CO2 critical point (Hitachi HCP-2), the samples were mounted on stubs with double-coated conductive carbon tape. Then samples were coated with Pt-Pd for 4 minutes in Hitachi E-102 ion sputter. The plants were observed and photographed with Hitachi S2500 Scanning Electron Microscope. The size distribution curve of particles adsorbed on the sample plant was measured and analyzed by ImageJ software [47]. The scale bar of each image was used to calibrate the software before particle measurements. All particles size were obtained from a mean value of two values equivalent to the orthogonal axes in micrometers unit.

    The ash content and moisture content were 14.44 ± 1.54 and 6.32 ± 0.34, respectively. The results revealed that the lower the dry matter the higher the ash due to the moisture content affecting the initial weight of the sample as shown in Table 1. Owing to the ash content is the amount of inorganic noncombustible material, which is residues after a completely burnt, such as oxides of the inorganic elements including metal salts which are important for photosynthesis processes requiring ions such as Na+ (Sodium), K+ (Potassium), and Ca2+ (Calcium) [48], and trace minerals which are required for unique molecules, such as chlorophyll and hemoglobin [49]. The obtained value (14.44%) shows a higher as compared to Spanish mosses (3.7%–7.0%) [50]. It means Riccia can be absorbed some minerals or trace elements and are essential for many processes. Moreover, the plants may have a short or long shelf life depending on the moisture content. The plant saps reflect the ability of a plant to exist under widely different environmental conditions, and the ability to adjust its tissue fluids to the new environment [51]. Riccia has a higher moisture content as compared to some plants such as Rauvolfia serpentina (0.73%) [52] and Sphagnum spp. (0.70%) [53]. From the above mentioned, Riccia has the potential to absorb the particulate matter and has a high tolerance to the stress condition.

    Table 1.  The ash content and moisture content.
    Parameter Value
    Ash content 14.44 ± 1.54
    Moisture content 6.32 ± 0.34

     | Show Table
    DownLoad: CSV

    The RWC showed a high value of 87.32 ± 4.15%, which RWC referred to the current water content of the leaf tissue relative to the full turgid water content, normally RWC in leaves ranges between 98%, in turgid and transpiring leaves to around 20%–50%, and RWC at initial wilting is about 60% to 70% [54]. The high level of relative water content leads to a high tolerance ability to air pollutants owing to high RWC to retain its physiological balance and prevent the loss of water and nutrients [55]. The pH of an extracted plant is neutral with a value of 7.2 ± 0.1. Plants with a low pH showed high sensitivity to air pollutants, especially SO2 and NO2. Their stomata were faster closed when exposed to air pollutants [31]. pH adjustment plays an important role in the movement of vesicles and intracellular compounds. When plants are exposed to air pollutants (especially SO2), large amounts of H+ are produced in the intercellular fluid to react with SO2, therefore this H+ may produce H2SO4, leading to decreased pH of the leaf extract. Riccia exhibited a low chlorophyll concentration of 0.1542 ± 0.0124 mg/g fr wt. Air pollution has a great influence on the total chlorophyll of plants because the photosynthetic pigment has been degraded when exposed to air pollutants. Therefore, the total chlorophyll of plants decreases as air pollution increases. Ascorbic acid content was found less with a value of 0.74 ± 0.11 mg/g fr wt. Ascorbic acid revealed a strong reducer and plays a role in photosynthetic carbon fixation, cell wall synthesis, defence system, and cell division. The antioxidant capacity of ascorbic acid affects its resistance to air pollution. Plant with a low ascorbic acid content has a high sensitivity to air pollution. In the study, all four parameters can be calculated the APTI value which indicated that Riccia showed a sensitive level to air pollutants with a value of 9.27. It can be used as a bioindicator of air pollution. Physiological and biochemical properties are presented in Table 2.

    Table 2.  Physiological and biochemical properties.
    Parameter Value
    Relative water content (RWC) 87.32 ± 4.15%
    Plant extract pH 7.2 ± 0.1
    Total chlorophyll content (TCH) 0.1542 ± 0.0124 mg/g fr wt.
    Chlorophyll A 0.1142 ± 0.0090 mg/g fr wt.
    Chlorophyll B 0.0400 ± 0.0034 mg/g fr wt.
    Ascorbic acid (AA) 0.74 ± 0.11 mg/g fr wt.
    Air pollution tolerance index (APTI) 9.27

     | Show Table
    DownLoad: CSV

    Dust collecting capacity of plants depends on their morphological traits like phyllotaxy, leaf shape, size and texture, presence or absence of hairs, stomatal frequency, surface roughness, stomatal density, and trichome length. The median dust-capturing potential of Riccia has a higher value of 0.0667 ± 0.0172 g/cm2 (Table 3) as compared to greening plant species (0.387–0.451 g/cm2) [56]. The Riccia revealed dust-capturing potential to reduce particulate matter, therefore it acts as a biological air filter.

    Table 3.  Dust capturing potential.
    Dust capturing potential Value
    Median 0.0667 ± 0.0172 g/cm2
    Interquartile range, IQR 0.0121

     | Show Table
    DownLoad: CSV

    The adsorption efficiency of Riccia with four different velocities of 0.5, 1.0, 1.5, and 2.0 m/s was in the range of 2.3%–49.6% as shown in Table 4. The highest efficiency values for PM1, PM2.5, and PM10 were 31.4, 40.1, and 49.6, respectively, which belonged to the horizontal panel with a velocity of 2.0 m/s. As compared to the previous results from Epipremnum aureum with the removal efficiencies of 75.2% for PM2.5, and 71.9% for PM10 [57] and Tillandsia usneoides exhibited the removal efficiency of 16.5% for PM2.5 from incense and 9.2% for A1 rigid, respectively [58]. From our previous work [38], the CFD calculation shows the optimum condition for trapped particulate matter was alternating zigzags panels due to the air pollution in all air inlets that can be contacted with the surface area of the study plant at all twelve positions. On the other hand, the results show that the horizontal panel has the highest values might be from the air circulation has no blocking from the panel, while the zigzags have partially blocked. Owing a higher airflow rate creates turbulent conditions in the treatment box, in which the motion particles display transverse motion, all particles oscillate along paths at right angles to the direction of the wave's advance, and as a result, all 12 points of the sample plant in the horizontal panel were contacted with the particulate matter leading to high adsorption capacity [59,60]. It can be concluded that the alignment of the panel affects the efficiency. Moreover, the air velocities affect the adsorption efficiency. Figure 1 demonstrates the effect of different air velocities on the adsorption efficiency which showed that the increased velocity increased the efficiency. Generally, higher velocities in the gas phase refer to the diffusion rates being higher leading to a shorter residence time, however, the results showed that the increased air velocities increased adsorption efficiency because high airflow velocity causes airflow turbulence and generates transverse motion which enhances the rate of energy and momentum exchange increasing heat transfer and friction coefficient leading to gas-solid mass transfer [61]. As a result, air circulation around the treatment box contributes to mass transfer between air pollutants and the study plant at all positions affecting the higher dust absorption.

    Table 4.  The adsorption efficiency with four different velocities of 0.5, 1.0, 1.5, and 2.0 m/s.
    Air velocity (m/s) Panel alignment Adsorption efficiency (%)
    PM1 PM2.5 PM10
    0.5 Vertical 2.3 ± 4.8 24.4 ± 1.7 29.6 ± 3.0
    Zigzag 3.3 ± 4.4 3.7 ± 2.0 7.1 ± 5.9
    Horizontal 16.6 ± 3.7 31.5 ± 12.6 15.5 ± 4.7
    1.0 Vertical 7.3 ± 3.6 14.9 ± 7.2 30.2 ± 4.3
    Zigzag 15.6 ± 3.3 26.5 ± 6.9 36.2 ± 4.6
    Horizontal 12.5 ± 5.6 22.3 ± 10.3 39.9 ± 6.2
    1.5 Vertical 19.6 ± 4.9 21.3 ± 6.0 23.6 ± 5.2
    Zigzag 3.3 ± 3.5 6.1 ± 3.1 18.6 ± 2.6
    Horizontal 11.0 ± 10.5 29.0 ± 12.7 40.3 ± 14.9
    2.0 Vertical 12.0 ± 4.7 21.8 ± 9.0 25.8 ± 1.7
    Zigzag 12.6 ± 4.0 15.9 ± 1.6 36.2 ± 4.5
    Horizontal 31.4 ± 5.2 40.1 ± 10.8 49.6 ± 9.6

     | Show Table
    DownLoad: CSV
    Figure 1.  The absorption efficiency of Riccia for capturing PM1, PM2.5, and PM10 using vertical, zigzag, and horizontal airflow patterns at a velocity of (a) 0.5, (b) 1.0, (c) 1.5, and (d) 2.0 m/s.

    The chromatogram of the particulate matter treated and untreated Riccia was displayed in Figure 2. The chromatogram of the particulate matter treated shows several peaks of an interesting compound with the retention time as shown in Figure 2(a). The mass spectra of interested compounds adsorbed on Riccia were identified by comparison with the NIST database. The interesting chemical constituent of substances adsorbed on Riccia with retention time 2.964–9.892 min and predominant fragment at m/z 156–m/z 278 were shown in Table 5. According to the result of the HS-GC-MS chromatogram and mass spectra, Riccia can be adsorbed by the particulate matter of several compounds related to air pollutants.

    Figure 2.  HS-GC-MS chromatogram (a) particulate matter treated (b) particulate matter untreated.
    Table 5.  The chemical constituent of substances adsorbed on Riccia compared to the NIST database.
    RT (min) Compound name m/z CAS number Structure
    2.964 Undecane 156 1120-21-4
    3.670 Butylbenzene 134 104-51-8
    4.099 Dodecane 170 112-40-3
    4.244 2, 6-dimethyl-Undecane 184 17301-23-4
    4.710 2-methyl-1-octanol 143 818-81-5
    4.865 1, 2, 3, 4-tetrahydro-6-methyl- Naphthalene 146 1680-51-9
    5.310 Tridecane 184 629-50-5
    6.210 2, 6, 10-trimethyl-Dodecane 212 3891-98-3
    6.515 Tetradecane 198 629-59-4
    7.686 Pentadecane 212 629-62-9
    8.817 Sulfurous acid, hexyl octyl ester 278 Not available
    9.892 Sulfurous acid, 2-ethyl hexyl isohexyl ester 278 Not available

     | Show Table
    DownLoad: CSV

    Cd, Pb, and Na were determined using a Flame atomic absorption spectrophotometer (Perkin-Elmer, PinAAcle 900F). The hollow cathode lamps for Cd, Pb, and Na were used as radiation sources and the fuel was air acetylene. All the standards and samples were run in triplicates. Metals concentrations of Riccia and its blank were determined using FAAS and the results obtained were shown in Table 6. Riccia shows the quantity of Cd, Pb, and Na 0.0044 ± 0.0069 mg/gDW, 0.0208 ± 0.0278 mg/gDW, and 0.9395 ± 0.1009 mg/gDW, respectively. The results suggest that Riccia can be adsorbed with some heavy metal in the environment. However, the quantity of analyte metal from particulate matter treated and untreated Riccia shows the result is not different.

    Table 6.  The concentration of metal analyte in plant samples.
    Metal Plant samples
    Riccia Blank of Riccia
    Cd (mg/gDW) 0.0044 ± 0.0069 0.0048 ± 0.0019
    Pb (mg/gDW) 0.0208 ± 0.0278 0.0320 ± 0.0464
    Na (mg/gDW) 0.9395 ± 0.1009 0.9800 ± 0.2700

     | Show Table
    DownLoad: CSV

    The feature of Riccia, a non-vascular land plant, classified into the family Ricciaceae of Marchantiophyta, is distinctive in having oil bodies and non-sporogenous cells. It can be easily identified by a flat, dichotomously branched, ribbon-like, and fleshy thallus in the gametophytic stage as shown in Figure 3 (the upper left and the middle). Scanning electron micrographs of liverworts exhibited cells occupied by arbuscules in the main colonization zone within the thallus tissue which is a flattened mass of tissue. After air pollutants were deposited on plant samples in the gas form via wet deposition, the SEM image revealed that the particles covered the thallus surface. The particle size was obtained from a mean value of two values equivalent to the orthogonal axes. The final values of all particles were distributed by the size ranging from less than 1–4 μm.

    Figure 3.  The photographs of Riccia show a portion of fleshy thalli (the upper left), a close-up of thallus showing oil bodies (the upper middle), SEM of the dorsal view of thalli before being treated showing globose-subglobose epidermal cells (the upper right), SEM photographs of thalli in cross-section after being treated (the lower left and middle) and the distribution curve showing the particle size covered on the thallus surface (the lower right). Scanning electron micrograph before deposited air pollutants at 500× magnification, SEM photograph of the particles deposited over the surface at 250× magnification and 2000× magnification.

    In this study, the effectiveness of Riccia as a biofilter for adsorbed particulate matter in the urban ecosystem was investigated. Due to its morphology which is a complex thalloid, oil leaf, and non-sporogenous cells, and its ecology, which adsorbs the nutrients and water from the air through its outer surfaces [39,40], therefore Riccia was selected as a study plant as a biofilter of pollutants. The physical properties, physiological, and biological were measured to evaluate the particulate matter capture potential. The results revealed Riccia can be absorbed some minerals or trace elements through its surface via the stomata. Furthermore, Riccia can exist in unfavourable environmental conditions and adjust its tissue fluids to the new environment. Besides, Riccia was classified in a sensitive group with a lower air pollution tolerance index (APTI) value of 9.27 which is indicated that Riccia can be used as a biomonitoring plant for air pollution. According to the previous study, the ability to capture the elements of plants depending on their leaf properties [62] such as phyllotaxy, leaf shape, size and texture, presence or absence of hairs, stomatal frequency, stomatal density, surface roughness, and trichome length [63]. The biofiltration process consists of two major steps; are sorption and biodegradation in the removal of the contaminants from air or water. In the sorption step, there are two major mechanisms including adsorption by cuticle and uptake by opened stomata [64]. Epipremnum aureum was via stomata of plants and elected to use as a biofilter with the removal efficiencies of 85% for total suspended particles, 75.2% for PM2.5, and 71.9% for PM10 [57]. The use of Tillandsia usneoides as biofilter exhibited higher removal efficiency of 16.5% and 9.2% in PM2.5 from incense and A1 rigid, respectively. Moreover, T. usneoides with trichome structure exhibited larger removal efficiencies than without trichome, which indicates that the increase in the total effective surface was effective for the deposition of PM [58].

    Herein, the removal efficiency of the particulate matter was determined and expressed to the adsorption efficiency. The alignment of the plant in the treatment box is three types, horizontal panels, zigzags panels, and vertical panels with four different velocities of 0.5, 1.0, 1.5, and 2.0 m/s. The highest efficiency values for PM1, PM2.5, and PM10 were 31.4, 40.1, and 49.6, respectively, which belonged to the horizontal panel with a velocity of 2.0 m/s due to the high airflow velocity creating turbulence and generating transverse motion which enhances the rate of energy and momentum exchange increasing heat transfer and friction coefficient leading to gas-solid mass transfer between air pollutants and sample plants via stomata. The VOCs and PAHs were analyzed from the plants after being treated with air pollutants generates by diesel B7. The result revealed twelve compounds from treated Riccia as compared to untreated Riccia, which indicate Riccia can be adsorbed with VOCs and PAHs. Consequently, metals such as Cd, Pb, and Na were determined and exhibited the analyte metal in treated Riccia, however, the quantity of analyte metal from treated and untreated Riccia is not significantly different. The morphology of Riccia shows a portion of fleshy thalli, oil bodies, globose-subglobose epidermal cells, and roughness surface which is characteristic of Riccia. Additionally, the treated Riccia shows the particulate matter cover on the thallus surface, and the particle sizes were in the range of 1–4 μm. These results indicate that the removal efficiency of particulate matter using Riccia is effective owing to this method being sustainable, easy-to-manage, and eco-friendly practices.

    Due to peculiar morphological and physiological characteristics, mosses are very suitable adsorbents for a wide variety of pollutants (i.e., metals and metalloids, PAHs, radionuclides) and, when used as transplants in nylon bags, they allow easy monitoring potentially of any site, with a highly dense sampling network. Nevertheless, the technique uses moss-bag methods for the biomonitoring of air pollution but there are still limitations and less metal uptake occurs. The other variables (climate, bag size, shape, exposure time, and height) had low or no influence at all [65]. Innovation of this study used Computational Fluid Dynamics (CFD) as a modeling technique that can predict in terms of area and time. The research on airflow and indoor air circulation has been studied and applied [66,67]. CFDs can provide accurate information about flow distribution and concentrations in modeling domains, rather than only target areas for data collection [68,69]. Model boxes were constructed to study airflow data through 12 locations in the treatment box model, such as air velocity, flow characteristics, and air temperature distribution. Air in treatment box model 2 was circulated through every location. The airflow average that passed through 12 locations was similar throughout the treatment mode at all velocities of inlet air, airflow passed through all locations and contacted the surfaces of every panel. As a result, there was a high being absorbed by air pollution treatment by plants such as mosses and lichen to assess the concentration of heavy metal in the atmosphere [38]. In the future, possibility to apply Riccia as a biofilter of a large-scale urban site to detect emission sources and quantitative analysis of airborne persistent pollutants.

    The authors would like to express their appreciation for the grant from the National Research Council of Thailand (NRCT) and Udon Thani Rajabhat University which all supported this paper. We would like to thank Prof. Vinich Promarak, the Smart Materials and Semiconductor Devices Laboratory, Vidyasirimedhi Institute of Science and Technology (VISTEC) for SEM analysis assistance. We also thank Assoc. Prof. Sahut Chantanaorrapint, PSU Herbarium, Division of Biological Science, Faculty of Science Prince of Songkla University for plant identification.

    The authors declare no conflict of interest.



    [1] L. Glass, J. D. Murray, Interdisciplinary Applied Mathematics: Mathematical Biology I, New York: Springer, 2001.
    [2] M. Constantin, D. Gheorghe, J. Tenreiro, Introduction to Fractional Differential Equations, New York: Springer, 2019.
    [3] P. O. Mohammed, R. P. Agarwal, I. Brevik, M. Abdelwahed, A. Kashuri, M. A. Yousif, On Multiple-Type Wave Solutions for the Nonlinear Coupled Time-Fractional Schrödinger Model, Symmetry, 16 (2024), 553.
    [4] L. Sadek, A Cotangent Fractional Derivative with the Application, Fractal Fractional, 7 (2023), 444.
    [5] L. Sadek, B. Abouzaid, E. M. Sadek, H. T. Alaoui, Controllability, observability and fractional linear-quadratic problem for fractional linear systems with conformable fractional derivatives and some applications, Int. J. Dyn. Control, 11 (2023), 214–228.
    [6] L. Sadek, T. A. Lazǎr, On Hilfer cotangent fractional derivative and a particular class of fractional problems, AIMS Mathematics, 8 (2023), 28334–28352. https://doi.org/10.3934/math.20231450 doi: 10.3934/math.20231450
    [7] A. M. S. Mahdy, Numerical solutions for solving model time-fractional Fokker–Planck equation, Numer. Methods Partial Differ. Eq., 37 (2021), 1120–1135.
    [8] S. Noor, B. M. Alotaibi, R. Shah, S. M. Ismaeel, On the Solitary Waves and Nonlinear Oscillations to the Fractional Schrödinger–KdV Equation in the Framework of the Caputo Operator, Symmetry, 15 (2023), 1616.
    [9] S. Noor, A. S. Alshehry, N. H. Aljahdaly, H. M. Dutt, I. Khan, R. Shah, Investigating the Impact of Fractional Non-Linearity in the Klein–Fock–Gordon Equation on Quantum Dynamics, Symmetry, 15 (2023), 881. https://doi.org/10.3390/sym15040881 doi: 10.3390/sym15040881
    [10] S. Noor, M. A. Hammad, R. Shah, A. W. Alrowaily, S. A. El-Tantawy, Numerical Investigation of Fractional-Order Fornberg–Whitham Equations in the Framework of Aboodh Transformation, Symmetry, 15 (2023), 1353.
    [11] N. Attia, A. Akgül, D. Seba, A. Nour, On solutions of time-fractional advection–diffusion equation, Numer. Methods Partial Differ. Eq., 39 (2023), 4489–4516.
    [12] M. Mulimani, S. Kumbinarasaiah, Numerical solution of time-fractional telegraph equations using wavelet transform, Int. J. Dynam. Control, 2023. https://doi.org/10.1007/s40435-023-01318-y doi: 10.1007/s40435-023-01318-y
    [13] S. O. Edeki, G. O. Akinlabi, S. A. Adeosun, Analytic and Numerical Solutions of Time-Fractional Linear Schrödinger Equation, Commun. Math. Appl., 7 (2016), 1–10.
    [14] D. Li, W. Sun, C. Wu, A novel numerical approach to time-fractional parabolic equations with nonsmooth solutions, Numer. Math., 14 (2021), 355–376.
    [15] A. A. Alderremy, R. Shah, N. A. Shah, S. Aly, K. Nonlaopon, Comparison of two modified analytical approaches for the systems of time fractional partial differential equations, AIMS Mathematics, 8 (2023), 7142–7162. http://doi.org/%2010.3934/math.2023360
    [16] Z. Guang-an, Numerical solutions to time-fractional stochastic partial differential equations, Nume. Algorithms, 82 (2019), 553–571.
    [17] M. Alaroud, O. Ababneh, N. Tahat, S. Al-Omari, Analytic technique for solving temporal time-fractional gas dynamics equations with Caputo fractional derivative, AIMS Mathematics, 7 (2022), 17647–17669. http://doi.org/10.3934/math.2022972 doi: 10.3934/math.2022972
    [18] S. Das, R. Kumar, Approximate analytical solutions of fractional gas dynamic equations, Appl. Math. Comput., 217 (2011), 9905–9915.
    [19] A. Esen, O. Tasbozan, Cubic B-spline collocation method for solving time fractional gas dynamics equation, Tbilisi Math. J., 8 (2015), 221–231.
    [20] A. Esen, O. Tasbozan, An approach to time fractional gas dynamics equation: Quadratic B-spline Galerkin method, Appl. Math. Comput., 261 (2015), 330–336.
    [21] R. Noureen, M. N. Naeem, D. Baleanu, P. O. Mohammed, M. Y. Almusawa, Application of trigonometric B-spline functions for solving Caputo time fractional gas dynamics equation, AIMS Mathematics, 8 (2023), 25343–25370. http://dx.doi.org/%2010.3934/math.20231293
    [22] K. Shah, T. Singh, A. Kılıçman, Combination of integral and projected differential transform methods for time-fractional gas dynamics equations, Ain Shams Eng. J., 9 (2018), 1683–1688.
    [23] A. Prakash, M. Kumar, Numerical Method for Time-Fractional Gas Dynamic Equations, Proc. Nat. Acad. Sci. Ind. Sec. A – Phys. Sci., 89 (2019), 559–570.
    [24] K. M. Saad, E. H. AL-Shareef, M. S. Mohamed, X. Yang, Optimal q-homotopy analysis method for time-space fractional gas dynamics equation, Eur. Phys. J. Plus, 132 (2017), 23.
    [25] S. S. Zhou, N. A. Shah, I. Dassios, S. Saleem, K. Nonlaopon, A comparative analysis of fractional-order gas dynamics equations via analytical techniques, Mathematics, 9 (2021), 1735. https://doi.org/10.3390/math9151735 doi: 10.3390/math9151735
    [26] R. Khalil, M. Al Horani, A. Yousef, M. Sababheh, A new definition of fractional derivative, J. Comput. Appl. Math., 264 (2014), 65–70.
    [27] T. Abdeljawad, On conformable fractional calculus, J. Comput. Appl. Math., 279 (2015), 57–66.
    [28] M. A. Yousif, F. K. Hamasalh, Conformable non-polynomial spline method: A robust and accurate numerical technique, Ain Shams Eng. J., 15 (2024), 102415.
    [29] A. Jhangeer, M. Muddassar, M. Kousar, B. Infal, Multistability and Dynamics of Fractional Regularized Long Wave equation with Conformable Fractional Derivatives, Ain Shams Eng. J., 12 (2021), 2153–2169.
    [30] L. Pedram, D. Rostamy, Numerical solutions of the initial boundary value problem for the perturbed conformable time Korteweg-de Vries equation by using the finite element method, Numer. Methods Partial Differ. Eq., 37 (2021), 1449–1463.
  • This article has been cited by:

    1. Doaa Filali, Mohammad Dilshad, Mohammad Akram, Nonlinear Contractions Employing Digraphs and Comparison Functions with an Application to Singular Fractional Differential Equations, 2024, 13, 2075-1680, 477, 10.3390/axioms13070477
    2. Doaa Filali, Mohammad Dilshad, Mohammad Akram, Weak ψ-Contractions on Directed Graphs with Applications to Integral Equations, 2024, 12, 2227-7390, 2675, 10.3390/math12172675
    3. Doaa Filali, Mohammad Akram, Mohammad Dilshad, A Class of φ-Contractions in Orthogonal Metric Spaces with an Application, 2024, 16, 2073-8994, 1462, 10.3390/sym16111462
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1152) PDF downloads(44) Cited by(7)

Figures and Tables

Figures(15)  /  Tables(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog