The aim of the current study is to reduce marijuana use among the general population. Because marijuana is an illegal narcotic with numerous negative health effects, it continues to pose a severe threat to public health in emerging nations. In this article, a modified mathematical model of the non-users, experimental users, recreational users, and addict's (NERA) model for marijuana consumption is established by incorporating a new compartment that represents the individuals who are being moved to jail by police intervention. The overall population of humans is divided into five main components: the non-smoker's compartment, experimental smoker's compartment, recreational smoker's compartment, addicted smoker's compartment, and prisoner's compartment. The novelty of this work is to modify the NERA model for marijuana consumption and validate the modified model. Furthermore, with the help of sensitivity analysis, control strategies for marijuana consumption in the population are addressed. The invariant region and the basic reproductive number (R0) are those parts that are needed for the validation of the proposed model. For the numerical simulation of the given model, the 4th-order Runge Kutta method will be used with the help of MATLAB to examine how the control strategies will play a role in marijuana consumption.
Citation: Atta Ullah, Hamzah Sakidin, Shehza Gul, Kamal Shah, Yaman Hamed, Maggie Aphane, Thabet Abdeljawad. Sensitivity analysis-based control strategies of a mathematical model for reducing marijuana smoking[J]. AIMS Bioengineering, 2023, 10(4): 491-510. doi: 10.3934/bioeng.2023028
| [1] | M. Botros, E. A. A. Ziada, I. L. EL-Kalla . Semi-analytic solutions of nonlinear multidimensional fractional differential equations. Mathematical Biosciences and Engineering, 2022, 19(12): 13306-13320. doi: 10.3934/mbe.2022623 |
| [2] | Jian Huang, Zhongdi Cen, Aimin Xu . An efficient numerical method for a time-fractional telegraph equation. Mathematical Biosciences and Engineering, 2022, 19(5): 4672-4689. doi: 10.3934/mbe.2022217 |
| [3] | Yang Liu, Jie Liu, Tao Yu . Sharp conditions for a class of nonlinear Schrödinger equations. Mathematical Biosciences and Engineering, 2023, 20(2): 3721-3730. doi: 10.3934/mbe.2023174 |
| [4] | Muhammad Nadeem, Ji-Huan He, Hamid. M. Sedighi . Numerical analysis of multi-dimensional time-fractional diffusion problems under the Atangana-Baleanu Caputo derivative. Mathematical Biosciences and Engineering, 2023, 20(5): 8190-8207. doi: 10.3934/mbe.2023356 |
| [5] | Xiaomeng Ma, Zhanbing Bai, Sujing Sun . Stability and bifurcation control for a fractional-order chemostat model with time delays and incommensurate orders. Mathematical Biosciences and Engineering, 2023, 20(1): 437-455. doi: 10.3934/mbe.2023020 |
| [6] | Noura Laksaci, Ahmed Boudaoui, Seham Mahyoub Al-Mekhlafi, Abdon Atangana . Mathematical analysis and numerical simulation for fractal-fractional cancer model. Mathematical Biosciences and Engineering, 2023, 20(10): 18083-18103. doi: 10.3934/mbe.2023803 |
| [7] | Ilse Domínguez-Alemán, Itzel Domínguez-Alemán, Juan Carlos Hernández-Gómez, Francisco J. Ariza-Hernández . A predator-prey fractional model with disease in the prey species. Mathematical Biosciences and Engineering, 2024, 21(3): 3713-3741. doi: 10.3934/mbe.2024164 |
| [8] | Swadesh Pal, Malay Banerjee, Vitaly Volpert . Spatio-temporal Bazykin’s model with space-time nonlocality. Mathematical Biosciences and Engineering, 2020, 17(5): 4801-4824. doi: 10.3934/mbe.2020262 |
| [9] | Xian Zhang, Chen Huang . Existence, multiplicity and non-existence of solutions for modified Schrödinger-Poisson systems. Mathematical Biosciences and Engineering, 2023, 20(2): 3482-3503. doi: 10.3934/mbe.2023163 |
| [10] | Debao Yan . Existence results of fractional differential equations with nonlocal double-integral boundary conditions. Mathematical Biosciences and Engineering, 2023, 20(3): 4437-4454. doi: 10.3934/mbe.2023206 |
The aim of the current study is to reduce marijuana use among the general population. Because marijuana is an illegal narcotic with numerous negative health effects, it continues to pose a severe threat to public health in emerging nations. In this article, a modified mathematical model of the non-users, experimental users, recreational users, and addict's (NERA) model for marijuana consumption is established by incorporating a new compartment that represents the individuals who are being moved to jail by police intervention. The overall population of humans is divided into five main components: the non-smoker's compartment, experimental smoker's compartment, recreational smoker's compartment, addicted smoker's compartment, and prisoner's compartment. The novelty of this work is to modify the NERA model for marijuana consumption and validate the modified model. Furthermore, with the help of sensitivity analysis, control strategies for marijuana consumption in the population are addressed. The invariant region and the basic reproductive number (R0) are those parts that are needed for the validation of the proposed model. For the numerical simulation of the given model, the 4th-order Runge Kutta method will be used with the help of MATLAB to examine how the control strategies will play a role in marijuana consumption.
Atmospheric corrosion of steel embedded in reinforced concrete is the phenomenon that most influences the deterioration of structures, mainly those built in coastal areas due to the effect of marine spray. This is a worldwide problem that is very difficult to solve, and many researchers have raised the need to carry out studies that lead to prevent this phenomenon and its consequences. Real case studies are extremely important to indicate the factors of influence to be considered as variables in these models [1].
It must be taken into account that rebars embedded in concrete could be exposed not only to chlorides present in the concrete itself, but also to chlorides found in the environment that penetrate the concrete, being-in both cases-a durability problem for the building' service life. Additionally, carbonation of concrete or penetration of acidic gases into the concrete is other causes of reinforcement corrosion to consider [2].
The cost of global steel and concrete infrastructures is expected to exceed ＄300 billion by 2025 [3] and the National Association of Corrosion Engineers estimated the global cost of corrosion around ＄3 trillion per year, or 3.4% of global GDP [4], clearly demonstrating the need to address the problem scientifically.
As a guide, it should be noted that although the ambient temperature conditions are not particularly aggressive, the humidity and salinity conditions are deemed to be aggressive. The records obtained from Cádiz city observatory, provided by AEMET (State Meteorology Agency of Spain), show that annual temperatures range between 7.5 and 32.10 ℃, averaging 18.52 ℃, and relative humidity varies between 32 and 99%, with a relatively high average of 73.08%. In spite of these environmental conditions, the correct preservation of the reinforced concrete structures is guaranteed as long as the architectural project fulfils with quality and execution requirements indicated in the standards in force.
According to the limited data available (data from the original project could not be obtained), the concrete used is 25 N/mm2, manufactured with type-Ⅱ white cement with limestone addition whose compressive strength was 35 N/mm2. According to the legislation in force when the building was built [5], for concrete exposed to aerial marine environments, the cement content should have been 300 kg/m3, the water/cement ratio 0.55 and the reinforcements should have been arranged at a minimum depth of 40 mm. The knowledge of the thickness of the reinforcement coating is essential to determine the level of protection of the steel and identify potential causes of deterioration, since the barrier effect of the concrete protecting the reinforcement has a decisive impact in the deterioration processes over time.
The lack of coating would require an intervention on the façade that would cause the regeneration or strengthening of the protective barrier layer on the building as a whole. It is important to note that the aesthetic quality of the white concrete exposed to the open air is seriously affected by contamination due to the development of biological stains. In this case, the material has behaved as a low bioreceptivity concrete [6], since it does not present alterations due to colonization by biotic organisms.
In general, the specific conditions that allow chlorides to reach the limit state to initiate corrosion are not clearly established, since this depends on different circumstances, so it was not possible to introduce a single representative value in the models. Although numerous studies have been carried out to determine the critical level of chlorides in the reinforcement capable of initiating corrosion, the suggested range of variation is extremely broad, which makes it difficult to define a single value [7,8,9]. The chloride content in the concrete that reaches the reinforcement depends on the ion diffusivity-linked to the nature of the cement and the compactness of the concrete [10], and the solid geometry determined by the thickness of the reinforced concrete covering and the micro-cracking of the concrete [11].
The corrosion of the reinforcements usually involves two stages. During the first phase, two reactions take place simultaneously: the metal (Fe0) dissolves in certain areas, while at other points a cathodic reaction occurs, e.g., from hydroxyl ions (OH−) if oxygen content is high enough in interstitial water. In the second phase, the dissolved metal ion (Fe++) reacts with the anions in the water to form a precipitate that covers the steel. Observations show that if the concrete is sufficiently sound, i.e. if its pH remains high (pH > 11) and is not contaminated by substances such as chlorides, the coating product on the steel will be sufficiently watertight to prevent water from reaching the metallic surface, since the high alkalinity of the concrete protects the steel [12]. In the presence of chlorides, the passive layer of the steel can no longer be watertight; in this way, the steel corrodes due to dissolution, giving rise to corrosion if chloride content exceeds a certain threshold determined by the pH of the concrete [9].
Taking into account the studied characteristics of the concrete, an increase of w/c ratio contributes to lower chloride thresholds [13], due to the different conditions for the passivation of the reinforcement, higher moisture content and oxygen availability, lower binding capacity and lower hydroxide content in a porous network solution [8].
Taking as a reference the penetration models indicated in Annex 9 of the Code on Structural Concrete (EHE-2008), "Additional considerations on durability" [14], in case of corrosion the total time needed for a significant attack or degradation can be expressed as Eq 1:
| tL=ti+tp | (1) |
tL: Total time required for the attack or degradation to become significant.
ti: Time taken by the penetration front of the aggressive agent to reach the reinforcement thereby causing the onset of corrosion.
tp: Corrosion propagation period until the structural elements suffers significant degradation.
Chloride concentration is expressed in % of concrete weight, using the cement content in kg/m3 and the actual density of the concrete. The critical chloride concentration, in accordance with current legislation, is a value of 0.6% of the cement weight for checking the limit state in relation to the corrosion of the passive reinforcements. The limit state concept implies that, if it is exceeded, a structural failure may occur because the project does not reach the end of its useful life as aggressive agents generate a degree of deterioration that prevents the structure from behaving in accordance with the requirements according to which it has been planned [14].
It is well known that carbon dioxide (CO2) contained in air penetrates over time through the porous system of the concrete, combining with calcium hydroxide to form calcium carbonate. This circumstance causes concrete pH to decrease from the outside inward, eventually causing the disappearance of the passivating oxide layer that naturally protects the steel reinforcement. The speed of the process depends on these variables: relative humidity, open porosity, cement dosage and age of the concrete [15]. Therefore, under the action of mainly ambient humidity and oxygen, the effect of corrosion can develop more quickly; even more in the presence of chlorides, as is the case of the marine spray that affects the building studied [16,17].
The object of this study was a building constructed in 1996 on the seafront of Cádiz, a city located in southern Spain, between 50 and 100 m from the sea. The material used in the building envelope was white concrete reinforced with conventional carbon steel reinforcement. This characteristic makes the building unique and represents an added maintenance difficulty since the degradation processes are more evident in white walls. The proximity of the sea and the subsequent action of marine spray and saline condensation on the surface has had an important impact on the evolution and conservation of the façade materials (Figure 1). Thus, undocumented repairs have been carried out in some areas consisting in the application of external mortar coating in highly localized areas.
Some evidence of concrete fragments becoming detached due to reinforcement corrosion (Figure 2), prompting the need for a general study of the building's façades to evaluate the condition of the concrete and the reinforcement making up the exterior construction elements in order to diagnose their current state, forecast the behaviour of the materials and report on the necessary rehabilitation actions.
The purpose of this work was to determine the state of the materials comprising the exposed concrete façade of the building studied by extracting cores from different orientations of the building, determining their density and open porosity, the concrete coatings of the reinforcements, the depth of carbonation and chloride penetration, in order to evaluate the development of corrosion processes, as well as their origin and intensity. This study serves as a useful tool for the technicians entrusted with the maintenance and rehabilitation of the building, supporting decision-making in areas already damaged and in other areas where there are not yet external signs of deterioration.
A general prior analysis of the reinforcement was carried out in the studied areas using a GPR (ground penetrating radar) model GP8800 manufactured by PROCEQ with a defined frequency range between 400–6000 MHz and a penetration capacity for measuring up to 65 cm. The data were processed using Proceq GPR Live software enabling results to be read in situ. GPR is a non-destructive method that uses UHF (Ultra High Frequency) electromagnetic waves designed to detect the internal characteristics of a specific medium and basically responds to the principle of wave reflection in dielectric discontinuities [18]. In this case, it was used to detect the rebar location in the reinforced concrete, areas of discontinuity in reinforcement, lack of reinforcing bars or overlapping reinforcements. Horizontal and vertical façade readings were taken and 3D analyses were performed to obtain a simulated image of the layout and state of the reinforcement [19].
The measurements taken with the GPR were complemented with measurements taken on the surface with the Proetisa Elcometer 331 BH detector to accurately confirm the position and depth of the steel bars and determine average coating thicknesses in the sampled areas and surrounding areas.
To carry out the investigation, in agreement with the technicians responsible for the conservation and maintenance of the building, a total of six cylindrical cores measuring 10 cm in diameter by 15–18 cm in depth were extracted, analysed and tested (Figure 3). They were obtained and stored properly in accordance with the criteria and measures required by the UNE-EN 12504-1 standard [20]. Three two-centimetre-thick slices were cut from each core, from the outside to the inside, in order to characterize the chloride penetration front. The cuts were carried out in a dry state to avoid altering the concentration of soluble salts in the samples [21].
Real density, apparent density and open porosity of the concrete in vacuum conditions were determined in accordance with the UNE-EN 83980 standard [22]. From the extracted cores, samples were taken from the slices cut at different depths (2, 4 and 6 cm), analysing the penetration profile in each sample and chloride concentration by cement weight (Figure 4).
Chloride analysis was performed by ion chromatography on a Metrohm 930 Compact IC Flex (Herisau, Switzerland). For the calculation of the percentage of chlorides in cement weight, "Chlorides in Concrete" CTI Technical Note C2 (September 2002) [23,24] was used as a reference, which considered real density calculated in the porosity test. In accordance with this methodology, the samples were prepared by grinding 50 g of the respective samples to a size of less than 125 µm. These were then quartered to extract an amount of 20 g, which was stirred for 60 min in deionized water at room temperature and subsequently filtered until get a volume of 250 mL.
Mineralogical analysis was carried out by x-ray diffraction (XRD) in a Bruker-AXS model D8I-A25 diffractometer, equipped with a Cu Kα copper filament (λ = 1.5405 Ǻ), with Bragg-Brentano θ-θ configuration, filter nickel and Lynxeye linear detector, using the powder technique; and thermogravimetry (TGA) on a TA Instruments Q600 SDT Differential Scanning Calorimeter.
Finally, the carbonation front of the cores was measured using the phenolphthalein test, in accordance with the UNE 112011 standard [25], and carried out inside the specimens-cores, one for each extraction, previously broken by the generatrix using the indirect traction test, in accordance with the UNE-EN 12390-6 standard [26]. The pH assessment of the cores produces a magenta colour on the concrete if its pH is greater than 8.5 and is colourless if the pH is less than 8. The test was carried out in the laboratory on samples that were packed immediately after extraction and isolated to avoid their superficial carbonation (Figure 5).
Table 1 shows the average coating thicknesses measured in the sampling areas.
| Sample | Location | Average thickness (mm) |
| 1. | Inner courtyard façade (east) | ≥45 |
| 2. | Main façade (west) | 15–20 |
| 3. | Main façade (south) | 30–40 |
| 4. | Exterior ramp façade (west) | 15–20 |
| 5. | Roof façade (east) | ≥50 |
| 6. | Roof façade (west) | 15–20 |
DownLoad:
CSV
The minimum coating value recorded on the construction date of the building was 40 mm [1] Therefore, it is evident that, generally speaking, this parameter has not been fulfilled. These measurements may be one of the potential causes of the alterations present.
The readings taken with the GPR corroborated that the external reinforcement generally consisted of horizontal rebars arranged every 20 cm and vertical rebars every 20–25 cm. However, areas of discontinuity in the horizontal reinforcement were detected due to the lack of reinforcing bars between the different sections of wall executed during construction works. The GPR also allow to identify areas with defective bonding, thus voids were observed in the radargrams. These voids were mainly found on the south-west façade -exposed to greater action of the marine spray- and some of these areas presented an imminent risk of detachment. In addition, approximate measurements of the existing coatings were obtained in each reading, being < 3 cm in most of the areas studied which is less than recommended by national standards at the time of construction (EHE-2008) (Table 1) (Figures 6 and 7).
Figure 6 presents the readings of the vertical radargrams in four areas, showing the following items: (ⅰ) the horizontal rebars from the top reinforcement (hyperbolas with higher contrast) and the bottom reinforcement, (ⅱ) their poor coating, (ⅲ) the wall thickness and (ⅳ) significant alterations like voids (Figure 6B, C) and lack of reinforcement in specific areas (Figure 6D). Additionally, it could be seeing the location of sample 4 (Figure 6B) and sample 2 (Figure 6D), due to the different concrete used to replace the horizontal coring previously extracted.
In the case of Figure 7, it is shown the vertical reinforcement and clearly identify differences in the coatings thickness-this fact is obvious, in the distance between rebars from the same wall section, voids and even double reinforcements, probably due to previous reinforcement interventions.
All the extracted samples could be classified as porous and permeable concretes, unsuitable for exposure to a marine environment as all exceeded 15% porosity. This high porosity may have its origin in a high w/c ratio and poor execution due to insufficient compaction [27,28]. The results obtained for the physical properties of the elements tested are shown in Table 2.
| Sample | Real density (g/cm3) | Apparent density (g/cm3) | Open porosity (%) |
| 1. | 2.65 | 2.15 | 18.8 |
| 2. | 2.65 | 2.15 | 18.9 |
| 3. | 2.66 | 2.16 | 19.0 |
| 4. | 2.77 | 2.22 | 19.7 |
| 5. | 2.64 | 2.20 | 16.5 |
| 6. | 2.63 | 2.12 | 19.3 |
DownLoad:
CSV
It should be noted that the density results obtained are fully consistent values, since the actual density of calcite is 2.71 g/cm3 and concrete is mainly composed of crushed limestone aggregate.
The available data makes impossible to estimate the times indicated in Eq 1 in this building. Although the EHE-08 Code indicates that the critical concentration of chlorides in the reinforcement area is 0.6% by cement weight for passive reinforcements [14]. If this threshold value is taken as a reference, the data obtained in the analysis carried out are clearly well below that value. The results for chloride concentration obtained are shown in Table 3.
| Cores | Sample location | Chloride concentration (%) |
| 1 | 1.1 External sample | 0.039 |
| 1.2 Intermediate sample | 0.023 | |
| 1.3 Interior sample | 0.023 | |
| 2 | 2.1 External sample | 0.144 |
| 2.2 Intermediate sample | 0.061 | |
| 2.3 Interior sample | 0.019 | |
| 3 | 3.1 External sample | 0.250 |
| 3.2 Intermediate sample | 0.032 | |
| 3.3 Interior sample | 0.019 | |
| 4 | 4.1 External sample | 0.181 |
| 4.2 Intermediate sample | 0.095 | |
| 4.3 Interior sample | 0.079 | |
| 5 | 5.1 External sample | 0.058 |
| 5.2 Intermediate sample | 0.023 | |
| 5.3 Interior sample | 0.020 | |
| 6 | 6.1 External sample | 0.081 |
| 6.2 Intermediate sample | 0.095 | |
| 6.3 Interior sample | 0.043 |
DownLoad:
CSV
It is important to highlight that the samples from the south-west and south-east façades had higher chloride contents than the rest. The south-east façade is in a cul-de-sac that favours condensation. While the south-west façade is directly facing the sea, and both are subject to the prevailing winds at this location (Figure 1).
The diffractograms obtained for the samples taken from the exterior and interior areas of the cores did not allow the identification of any mineral phase that can be associated with chlorides (Figure 8), thus making it possible to consider that the soluble salts which have penetrated through the porous system of the material may have reacted with tricalcium aluminate, but not enough chloroaluminate has been produced to develop visible peaks in the XRD spectra.
In this case, it was very important to determine the penetration depth of the carbonation front to confirm if the area where the metal reinforcements were located had been reached. This fact occurs in all samples except in 1. and 5. leading to their depassivation and their subsequent oxidation and expansion due to the formation of iron hydroxides. It was observed that this process is causing cracking in these elements and detachment of fragments of this material, due to the expansion of the products resulting from the corrosion process. The results of the carbonation depth and tests are shown in Table 4.
| Sample | Carbonation depth (mm) | Reinforcement coating (mm) |
| 1. | 20–25 | 40–50 |
| 2. | 10–15 | 15–20 |
| 3. | 30 | 30–40 |
| 4. | 10–15 | 15–20 |
| 5. | 15–20 | ≥50 |
| 6. | 20 | 15–20 |
DownLoad:
CSV
As described by Sanjuan et al. (2015) [29] and considering the characteristic of the concrete used in this case study (CEM Ⅱ cement and a w/c ratio of 0.5), the estimated chloride penetration coefficient (KCl) was 8.05 mm2/year. The most widely used model to evaluate the degree of penetration of chlorine salts in concrete consists in the application of Fick's second law of diffusion of aggressive chlorides through the coating. In this regard, the Eq 2 retrieved from the Annex 9 of the EHE-2008 expresses in a simplified way the relationship between the depth reached by a given concentration of aggressive agent over time and permits calculating the initiation period:
| d=KCl⋅√t | (2) |
Equation 2 indicates the time required for critical chloride concentration to be reached at a given distance from the surface of the element by taking into account the chloride diffusion coefficient. According to this value, based on the current age of the building (23 years).The depth reached by the penetration front of Cl− would be 38.6 mm.
Considering in accordance with Annex 9 of the EHE-2008, a limit value of 0.6% in cement weight to check the limit state in relation to the corrosion of passive reinforcements, it can be observed that this value has not been reached in the walls analysed in this building. However, an additional circumstance has seriously conditioned the behaviour of the reinforcements and concrete: the resulting carbonation and the consequent depassivation of reinforcements placed at an insufficient distance from the surface.
If the values shown in Table 4 are analysed, it can be observed that in samples 2, 3, 4 and 6 the carbonation front has gone beyond the position of the reinforcements or limits with said position. Given this circumstance, it may be concluded that. Although chloride concentration values never exceed the service limit state value [30], chloride concentration is sufficient for significant corrosion of the reinforcements in those areas where carbonation has caused depassivation and the consequent lack of protection of the steel elements. Reinforcements were even observed where chloride action has triggered pitting corrosion processes in non-carbonated areas (Figure 9).
Previous studies on the structural concrete cores of the building's façade have revealed the start principle of carbonation, depassivation on the reinforcements and their subsequent corrosion as the main causes of deterioration. These processes derive from, among other factors, the environmental conditions to which the building is exposed, aggravated by the physical characteristics of the concrete and the insufficient thickness of the reinforcement coating in certain areas. It should be noted that the problem detected in this building stems not only from the marine environment, but also from poor execution and poor emplacement of the steel reinforcement bars. So, it is essential in these cases a careful execution following the technical specifications and standards to guarantee the maximum durability of the constructive elements, particularly in the façades because they are more exposed.
The comparative test of properties such as density and porosity allowed us to evaluate their suitability to provide adequate protection, based on our estimates of the factors responsible for the damage. It is important to highlight that it was possible to demonstrate that non-compliance with the requirements stipulated in current regulations has negatively influenced the behaviour and durability of the reinforced concrete elements.
In reinforced concrete structures close to the coast, the action of marine spray is often assumed to be the cause of corrosion of rebars, and carbonation processes are not adequately considered. It has been shown that carbonation depassivation phenomena in concrete greatly accelerate the aggressive action of chlorides. It was confirmed that the action of chlorides -by itself- has not been a determining factor in the degradation of the building's façade. The depassivation of the steel reinforcements, due to the carbonation of the concrete, was confirmed as a necessary factor for the start of the alteration of the reinforcements due to chloride activity.
The study of the carbonation front and the thickness of the coating is an essential tool for technicians that make it easier to establish intervention criteria for the repair of the façade. Based on the scientific and technical results obtained, a proposal for the protection of the façade surfaces can be defined and a waterproofing agent for its application on concrete elements can be selected. These measurements would prevent the humidity penetration and condensation processes, that favour the migration of a greater quantity of chlorides towards the interior. Additionally, a palliative treatment to inhibit corrosion should be considered. Finally, attention must be drawn to the importance of considering a periodic protective maintenance programme, considering the complex environmental conditions to which the building is exposed. In accordance with the products selected for protecting of the façade and the durability indicated by the supplier, the stakeholders and technicians responsible for the retrofitting must establish a maintenance programme that guarantees the effectiveness of the executed work over time.
This research has been carried out thanks to the financing of project PGC2018-093470-B-I00 of the Ministry of Science, Innovation and Universities of the Government of Spain. Authors would like to thank the Ⅵ PPIT-2021-I.3 from the University of Seville and the Research, Technology and Innovation Centre (CITIUS) of the University of Seville, in whose laboratories the characterization analyses were carried out.
All authors declare no conflicts of interest in this paper.
| [1] | Algahtani OJ, Zeb A, Zaman G, et al. (2015) Mathematical study of smoking model by incorporating campaign class. Wulfenia 22: 205-216. |
| [2] | Leong DP, Teo KK, Rangarajan S, et al. (2018) World Population Prospects 2019. Department of Economic and Social Affairs Population Dynamics. New York (NY): United Nations; 2019 (https://population.un.org/wpp/Download/, accessed 20 September 2020). The decade of healthy ageing. Geneva: World Health Organization. World 73: 362k2469. |
| [3] |
Fligiel SEG, Roth MD, Kleerup EC, et al. (1997) Tracheobronchial histopathology in habitual smokers of cocaine, marijuana, and/or tobacco. Chest 112: 319-326. https://doi.org/10.1378/chest.112.2.319 
|
| [4] |
BALDWIN GC, Tashkin DP, Buckley DM, et al. (1997) Marijuana and cocaine impair alveolar macrophage function and cytokine production. Am J Respir Crit Care Med 156: 1606-1613. https://doi.org/10.1164/ajrccm.156.5.9704146
|
| [5] |
Shay AH, Choi R, Whittaker K, et al. (2003) Impairment of antimicrobial activity and nitric oxide production in alveolar macrophages from smokers of marijuana and cocaine. J Infect Dis 187: 700-704. https://doi.org/10.1086/368370
|
| [6] |
Naz H, Dumrongpokaphan T, Sitthiwirattham T, et al. (2023) A numerical scheme for fractional order mortgage model of economics. Results Appl Math 18: 100367. https://doi.org/10.1016/j.rinam.2023.100367
|
| [7] | Alrabaiah H, Ullah Z, Ahmad I, et al. Iterative investigation of korteweg–de vries equation using ab derivative in caputo sense (2023). https://doi.org/10.1142/s0218348x23400959 |
| [8] |
Moore THM, Zammit S, Lingford-Hughes A, et al. (2007) Cannabis use and risk of psychotic or affective mental health outcomes: a systematic review. Lancet 370: 319-328. https://doi.org/10.1016/s0140-6736(07)61162-3
|
| [9] |
Vozoris NT, Zhu J, Ryan CM, et al. (2022) Cannabis use and risks of respiratory and all-cause morbidity and mortality: a population-based, data-linkage, cohort study. BMJ Open Respir Res 9: e001216. https://doi.org/10.1136/bmjresp-2022-001216
|
| [10] | Tashkin DP, Coulson AH, Clark VA, et al. (1987) Respiratory symptoms and lung function in habitual heavy smokers of marijuana alone, smokers of marijuana and tobacco, smokers of tobacco alone, and nonsmokers. Am Rev Respir Dis 135: 209-216. https://doi.org/10.1378/chest.78.5.699 |
| [11] | Polen MR, Sidney S, Tekawa IS, et al. (1993) Health care use by frequent marijuana smokers who do not smoke tobacco. West J Med 158: 596. https://doi.org/10.1007/978-1-4615-1907-2_99 |
| [12] | Merz F (2018) United nations office on drugs and crime: world drug report 2017. 2017. SIRIUS-Zeitschr Strateg Anal 2: 85-86. https://doi.org/10.1515/sirius-2018-0016 |
| [13] |
Degenhardt L, Charlson F, Ferrari A, et al. (2018) The global burden of disease attributable to alcohol and drug use in 195 countries and territories, 1990–2016: a systematic analysis for the global burden of disease study 2016. Lancet Psychiatry 5: 987-1012. https://doi.org/10.1016/s2215-0366(18)30337-7
|
| [14] | Al Juboori R, Subramaniam DS, Hinyard L (2022) Understanding the role of adult mental health and substance abuse in perpetrating violent acts: In the presence of unmet needs for mental health services. Int J Ment Health Addiction 21: 1-9. https://doi.org/10.1007/s11469-022-00778-1 |
| [15] |
Lee J, Thrul J (2021) Trends in opioid misuse by cigarette smoking status among US adolescents: Results from national survey on drug use and health 2015–2018. Prev Med 153: 106829. https://doi.org/10.1016/j.ypmed.2021.106829
|
| [16] |
Testai FD, Gorelick PB, Aparicio HJ, et al. (2022) Use of marijuana: effect on brain health: a scientific statement from the American Heart Association. Stroke 53: e176-e187. https://doi.org/10.1161/str.0000000000000396
|
| [17] |
Beutler JA, Marderosian AH (1978) Chemotaxonomy of Cannabis I. Crossbreeding between Cannabis sativa and C. ruderalis, with analysis of cannabinoid content. Econ Bot 32: 387-394. https://doi.org/10.1007/bf02907934
|
| [18] |
Di Forti M, Morgan C, Dazzan P, et al. (2009) High-potency cannabis and the risk of psychosis. Br J Psychiatry 195: 488-491. https://doi.org/10.1192/bjp.bp.109.064220
|
| [19] | Mousavi SE, Bozorgian A (2020) Investigation the kinetics of CO2 hydrate formation in the water system+ CTAB+ TBAF+ ZnO. Int J New Chem 7: 195-219. https://doi.org/10.33945/sami/ecc.2020.4.10 |
| [20] |
Ellis RJ (2009) Smoked medicinal cannabis for neuropathic pain in HIV: a randomized, crossover clinical trial. Neuropsychopharmacology 34: 672-680. https://doi.org/10.1038/npp.2008.120
|
| [21] |
Russo E (2001) Cannabinoids in pain management: Study was bound to conclude that cannabinoids had limited efficacy. BMJ 323: 1249-1250. https://doi.org/10.1136/bmj.323.7323.1249
|
| [22] |
Russo EB (2008) Cannabinoids in the management of difficult to treat pain. Ther Clin Risk Manag 4: 245-259. https://doi.org/10.2147/tcrm.s1928
|
| [23] |
Cascini F, Aiello C, Di Tanna GL (2012) Increasing delta-9-tetrahydrocannabinol (Delta-9-THC) content in herbal cannabis over time: systematic review and meta-analysis. Curr Drug Abuse Rev 5: 32-40. https://doi.org/10.2174/1874473711205010032
|
| [24] | Savaki HE, Cunha J, Carlini EA, et al. (1976) Pharmacological activity of three fractions obtained by smoking cannabis through a water pipe. Bull Narc 28: 49-56. https://doi.org/10.1007/978-3-642-51624-5_39 |
| [25] | Samimi A, Zarinabadi S, Bozorgian A (2021) Optimization of corrosion information in oil and gas wells using electrochemical experiments. Int J New Chem 8: 149-163. https://doi.org/10.33945/sami/jcr.2020.2.5 |
| [26] |
Izadi M, Yüzbaşı Ş, Ansari KJ (2021) Application of Vieta–Lucas series to solve a class of multi-pantograph delay differential equations with singularity. Symmetry 13: 2370. https://doi.org/10.3390/sym13122370
|
| [27] |
Tashkin DP, Gliederer F, Rose J, et al. (1991) Tar, CO and Δ9-THC delivery from the 1st and 2nd halves of a marijuana cigarette. Pharmacol Biochem Behav 40: 657-661. https://doi.org/10.1016/0091-3057(91)90378-f
|
| [28] |
Aldington S, Williams M, Nowitz M, et al. (2007) Effects of cannabis on pulmonary structure, function and symptoms. Thorax 62: 1058-1063. https://doi.org/10.1136/thx.2006.077081
|
| [29] |
Bloom JW, Kaltenborn WT, Paoletti P, et al. (1987) Respiratory effects of non-tobacco cigarettes. Br Med J (Clin Res Ed) 295: 1516-1518. https://doi.org/10.1136/bmj.295.6612.1516
|
| [30] | Sadr MB, Bozorgian A (2021) An overview of gas overflow in gaseous hydrates. J Chem Rev 3: 66-82. https://doi.org/10.22034/jcr.2021.118870 |
| [31] |
Cruickshank EK (1976) Physical assessment of 30 chronic cannabis users and 30 matched controls. Ann N Y Acad Sci 282: 162-167. https://doi.org/10.1111/j.1749-6632.1976.tb49895.x
|
| [32] |
Tashkin DP, Simmons MS, Sherrill DL, et al. (1997) Heavy habitual marijuana smoking does not cause an accelerated decline in FEV1 with age. Am J Respir Crit Care Med 155: 141-148. https://doi.org/10.1164/ajrccm.155.1.9001303
|
| [33] | Tashkin DP, Reiss S, Shapiro BJ, et al. (1977) Bronchial effects of aerosolized Δ9-tetrahydrocannabinol in healthy and asthmatic subjects. Am Rev Respir Dis 115: 57-65. https://doi.org/10.1007/978-1-4613-4286-1_8 |
| [34] |
Din A, Li Y, Liu Q (2020) Viral dynamics and control of hepatitis B virus (HBV) using an epidemic model. Alexandria Eng J 59: 667-679. https://doi.org/10.1016/j.aej.2020.01.034
|
| [35] |
McNulty W, Usmani OS (2014) Techniques of assessing small airways dysfunction. Eur Clin Respir J 1: 25898. https://doi.org/10.3402/ecrj.v1.25898
|
| [36] |
Roth MD, Arora A, Barsky SH, et al. (1998) Airway inflammation in young marijuana and tobacco smokers. Am J Respir Crit Care Med 157: 928-937. https://doi.org/10.1164/ajrccm.157.3.9701026
|
| [37] | Bozorgian A (2021) Exergy analysis for evaluation of energy consumptions in hydrocarbon plants. Int J New Chem 8: 329-344. https://doi.org/10.22034/ijnc.2020.123715.1105 |
| [38] | US Department of Health and Human Services.The benefits of smoking cessation. A report of the surgeon general (DHHS publication no. CDC 90-8416) (1990) . |
| [39] |
Khan H, Alzabut J, Shah A, et al. (2023) On fractal-fractional waterborne disease model: A study on theoretical and numerical aspects of solutions via simulations. Fractals 2023: 2340055. https://doi.org/10.1142/s0218348x23400558
|
| [40] |
Khan H, Ahmed S, Alzabut J, et al. (2023) A generalized coupled system of fractional differential equations with application to finite time sliding mode control for Leukemia therapy. Chaos Solitons Fractals 174: 113901. https://doi.org/10.1016/j.chaos.2023.113901
|
| [41] | Bozorgian A (2020) Investigation of hydrate formation phenomenon and hydrate inhibitors. J Eng Ind Res 1: 99-110. https://doi.org/10.22034/ijnc.2020.123715.1105 |
| [42] |
Din A, Li Y, Khan T, et al. (2020) Mathematical analysis of spread and control of the novel corona virus (COVID-19) in China. Chaos Solitons Fractals 141: 110286. https://doi.org/10.1016/j.chaos.2020.110286
|
| [43] |
Tashkin DP (2015) The respiratory health benefits of quitting cannabis use. Eur Respir Soc 46: 1-4. https://doi.org/10.1183/09031936.00034515
|
| [44] | Bozorgian A, Ghanavati B (2020) Removal of copper II from industrial effluent with beta zeolite nanocryst. https://doi.org/10.22034/pcbr.2022.328704.1213 |
| [45] |
Ahmed S, Azar AT, Tounsi M (2022) Design of adaptive fractional-order fixed-time sliding mode control for robotic manipulators. Entropy 24: 1838. https://doi.org/https://doi.org/10.3390/e24121838
|
| [46] | Rabipour S, Abdulkareem Mahmood E, Afsharkhas M (2022) Medicinal use of marijuana and its impacts on respiratory system. J Chem Lett 3: 86-94. https://doi.org/10.22034/jchemlett.2022.347365.1072 |
| [47] |
Biehl JR, Burnham EL (2015) Cannabis smoking in 2015: a concern for lung health?. Chest 148: 596-606. https://doi.org/10.1378/chest.15-0447
|
| [48] |
Ryan SA, Ammerman SD, O'Connor ME, et al. (2018) Marijuana use during pregnancy and breastfeeding: implications for neonatal and childhood outcomes. Pediatrics 142: e20181889. https://doi.org/10.1542/peds.2018-1889
|
| [49] |
Dauhoo M, Korimboccus B, Issack S (2013) On the dynamics of illicit drug consumption in a given population. IMA J Appl Math 78: 432-448. https://doi.org/10.1093/imamat/hxr058
|
| [50] | Ginoux JM, Naeck R, Ruhomally YB, et al. (2019) Chaos in a predator–prey-based mathematical model for illicit drug consumption. Appl Math Comput 347: 502-513. https://doi.org/10.1016/j.amc.2018.10.089 |
| [51] | Ullah A, Sakidin H, Gul S, et al. (2023) Sensitivity analysis-based validation of the modified NERA model for improved performance. J Adv Res Appl Sci Eng Technol 32: 1-11. https://doi.org/10.37934/araset.32.3.111 |
| [52] |
Zamir M, Zaman G, Alshomrani AS (2016) Sensitivity analysis and optimal control of anthroponotic cutaneous leishmania. PloS one 11: e0160513. https://doi.org/10.1371/journal.pone.0160513
|
| [53] |
Zamir M, Abdeljawad T, Nadeem F, et al. (2021) An optimal control analysis of a COVID-19 model. Alexandria Eng J 60: 2875-2884. https://doi.org/10.1016/j.aej.2021.01.022
|
| [54] |
Liu P, Din A, Zarin R (2022) Numerical dynamics and fractional modeling of hepatitis B virus model with non-singular and non-local kernels. Results Phys 39: 105757. https://doi.org/10.1016/j.rinp.2022.105757
|
| [55] |
Zamir M, Zaman G, Alshomrani AS (2017) Control strategies and sensitivity analysis of anthroponotic visceral leishmaniasis model. J Biol Dyn 11: 323-338. https://doi.org/10.1080/17513758.2017.1339835
|
| [56] | Zamir M, Sultana R, Ali R, et al. (2015) Study on the threshold condition for infection of visceral leishmaniasis. Sindh Univ Res J 47: 619-622. |
| [57] |
Liu P, Huang X, Zarin R, et al. (2023) Modeling and numerical analysis of a fractional order model for dual variants of SARS-CoV-2. Alexandria Eng J 65: 427-442. https://doi.org/10.1016/j.aej.2022.10.025
|
| [58] |
Zamir M, Nadeem F, Alqudah MA, et al. (2022) Future implications of covid-19 through mathematical modeling. Results Phys 33: 105097. https://doi.org/10.1016/j.rinp.2021.105097
|
| [59] |
Madhusudanan V, Srinivas MN, Murthy BSN, et al. (2023) The influence of time delay and Gaussian white noise on the dynamics of tobacco smoking model. Chaos Solitons Fractals 173: 113616. https://doi.org/10.1016/j.chaos.2023.113616
|
| [60] |
Bushnaq S, Shah K, Tahir S, et al. (2022) Computation of numerical solutions to variable order fractional differential equations by using non-orthogonal basis. AIMS Math 7: 10917-10938. https://doi.org/10.3934/math.2022610
|
| [61] | Guo X, Guo Y, Zhao Z, et al. (2022) Computing R0 of dynamic models by a definition-based method. Infect Dis Model 7: 196-210. https://doi.org/10.1016/j.idm.2022.05.004 |
| [62] | Van den Driessche P (2017) Reproduction numbers of infectious disease models. Infect Dis Model 2: 288-303. https://doi.org/10.1016/j.idm.2017.06.002 |
| [63] |
Rodrigues HS, Monteiro MTT, Torres DFM (2013) Sensitivity analysis in a dengue epidemiological model. Conference Papers in Mathematics 2013: 1-7. https://doi.org/10.1155/2013/721406
|
| [64] | Ahmad N, Charan S, Singh VP (2015) Study of numerical accuracy of Runge-Kutta second, third and fourth order method. Int J Comput Math Sci 4: 111-118. |
| [65] | Islam MA (2015) Accurate solutions of initial value problems for ordinary differential equations with the fourth order Runge Kutta method. J Math Res 7: 41. https://doi.org/10.5539/jmr.v7n3p41 |
| [66] | Aida M (2022) Fourth-order runge-kutta method for solving applications of system of first-order ordinary differential equations. Enhanced Knowl Sci Technol 2: 517-526. https://doi.org/10.1007/s40995-017-0258-1 |
| 1. | Yasin Asghari, Mostafa Eslami, Mashallah Matinfar, Hadi Rezazadeh, Novel soliton solution of discrete nonlinear Schrödinger system in nonlinear optical fiber, 2024, 90, 11100168, 7, 10.1016/j.aej.2024.01.024 | |
| 2. | Bo Liu, Zhou-Bo Duan, Li-Fang Niu, Soliton solutions of the TWPA-SNAIL transmission line circuit equation under continuum approximation via the Jacobi elliptic function expansion method, 2024, 98, 0973-7111, 10.1007/s12043-024-02791-6 |
Figures(6) / Tables(3)
Atta Ullah, Hamzah Sakidin, Shehza Gul, Kamal Shah, Yaman Hamed, Maggie Aphane, Thabet Abdeljawad. Sensitivity analysis-based control strategies of a mathematical model for reducing marijuana smoking[J]. AIMS Bioengineering, 2023, 10(4): 491-510. doi: 10.3934/bioeng.2023028
| Sample | Location | Average thickness (mm) |
| 1. | Inner courtyard façade (east) | ≥45 |
| 2. | Main façade (west) | 15–20 |
| 3. | Main façade (south) | 30–40 |
| 4. | Exterior ramp façade (west) | 15–20 |
| 5. | Roof façade (east) | ≥50 |
| 6. | Roof façade (west) | 15–20 |
DownLoad:
CSV
| Sample | Real density (g/cm3) | Apparent density (g/cm3) | Open porosity (%) |
| 1. | 2.65 | 2.15 | 18.8 |
| 2. | 2.65 | 2.15 | 18.9 |
| 3. | 2.66 | 2.16 | 19.0 |
| 4. | 2.77 | 2.22 | 19.7 |
| 5. | 2.64 | 2.20 | 16.5 |
| 6. | 2.63 | 2.12 | 19.3 |
DownLoad:
CSV
| Cores | Sample location | Chloride concentration (%) |
| 1 | 1.1 External sample | 0.039 |
| 1.2 Intermediate sample | 0.023 | |
| 1.3 Interior sample | 0.023 | |
| 2 | 2.1 External sample | 0.144 |
| 2.2 Intermediate sample | 0.061 | |
| 2.3 Interior sample | 0.019 | |
| 3 | 3.1 External sample | 0.250 |
| 3.2 Intermediate sample | 0.032 | |
| 3.3 Interior sample | 0.019 | |
| 4 | 4.1 External sample | 0.181 |
| 4.2 Intermediate sample | 0.095 | |
| 4.3 Interior sample | 0.079 | |
| 5 | 5.1 External sample | 0.058 |
| 5.2 Intermediate sample | 0.023 | |
| 5.3 Interior sample | 0.020 | |
| 6 | 6.1 External sample | 0.081 |
| 6.2 Intermediate sample | 0.095 | |
| 6.3 Interior sample | 0.043 |
DownLoad:
CSV
| Sample | Carbonation depth (mm) | Reinforcement coating (mm) |
| 1. | 20–25 | 40–50 |
| 2. | 10–15 | 15–20 |
| 3. | 30 | 30–40 |
| 4. | 10–15 | 15–20 |
| 5. | 15–20 | ≥50 |
| 6. | 20 | 15–20 |
DownLoad:
CSV
| Sample | Location | Average thickness (mm) |
| 1. | Inner courtyard façade (east) | ≥45 |
| 2. | Main façade (west) | 15–20 |
| 3. | Main façade (south) | 30–40 |
| 4. | Exterior ramp façade (west) | 15–20 |
| 5. | Roof façade (east) | ≥50 |
| 6. | Roof façade (west) | 15–20 |
| Sample | Real density (g/cm3) | Apparent density (g/cm3) | Open porosity (%) |
| 1. | 2.65 | 2.15 | 18.8 |
| 2. | 2.65 | 2.15 | 18.9 |
| 3. | 2.66 | 2.16 | 19.0 |
| 4. | 2.77 | 2.22 | 19.7 |
| 5. | 2.64 | 2.20 | 16.5 |
| 6. | 2.63 | 2.12 | 19.3 |
| Cores | Sample location | Chloride concentration (%) |
| 1 | 1.1 External sample | 0.039 |
| 1.2 Intermediate sample | 0.023 | |
| 1.3 Interior sample | 0.023 | |
| 2 | 2.1 External sample | 0.144 |
| 2.2 Intermediate sample | 0.061 | |
| 2.3 Interior sample | 0.019 | |
| 3 | 3.1 External sample | 0.250 |
| 3.2 Intermediate sample | 0.032 | |
| 3.3 Interior sample | 0.019 | |
| 4 | 4.1 External sample | 0.181 |
| 4.2 Intermediate sample | 0.095 | |
| 4.3 Interior sample | 0.079 | |
| 5 | 5.1 External sample | 0.058 |
| 5.2 Intermediate sample | 0.023 | |
| 5.3 Interior sample | 0.020 | |
| 6 | 6.1 External sample | 0.081 |
| 6.2 Intermediate sample | 0.095 | |
| 6.3 Interior sample | 0.043 |
| Sample | Carbonation depth (mm) | Reinforcement coating (mm) |
| 1. | 20–25 | 40–50 |
| 2. | 10–15 | 15–20 |
| 3. | 30 | 30–40 |
| 4. | 10–15 | 15–20 |
| 5. | 15–20 | ≥50 |
| 6. | 20 | 15–20 |
