Next Article in Journal
Study of Diels–Alder Reactions of Purpurogallin Tetraacetate with Various Dienophiles
Previous Article in Journal
Photovoltaic Properties and Microstructures of Polysilane-Added Perovskite Solar Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Targeting Desired Particle Size for Improved Dissolution and Manufacturability of Mefenamic Acid †

Continuous Manufacturing and Advanced Crystallization (CMAC) Future Research Hub, Loughborough University, Loughborough LE11 3TU, UK
*
Author to whom correspondence should be addressed.
Presented at the 3rd International Online Conference on Crystals, 15–30 January 2022; Available online: https://iocc_2022.sciforum.net/.
Chem. Proc. 2022, 9(1), 21; https://doi.org/10.3390/IOCC_2022-12165
Published: 15 January 2022
(This article belongs to the Proceedings of The 3rd International Online Conference on Crystals)

Abstract

:
This work aims at the production of uniformly sized, equant shaped crystals using a direct nucleation control (DNC) strategy to help achieve enhanced manufacturability and bioavailability of the active pharmaceutical ingredients (API). A crystallization system equipped with in situ and external wet-milling devices was used. The experiments were conducted using mefenamic acid (MFA) mixed into solvents of 2-butanol and heptane with a mass ratio of 90:10. Firstly, a linear cooling crystallization was conducted as a reference case and also to identify suitable DNC set-points, which were then implemented using closed-loop feedback control based on the measurement of the FBRM counts. The in situ and external wet-milling based DNC approaches were compared in terms of temperature cycling efficiency and final crystal size and shape. This integrated wet-milling crystallization system demonstrated enhanced control over the crystal size distribution (CSD) and morphology compared to conventional cooling crystallization processes. In addition, the expansion of the design space was achieved by integrating wet-mill to DNC crystallization compared to the cooling crystallization-only process. In addition, enhanced control of the rate of nucleation and attrition is accomplished through combining wet-milling and DNC.

1. Introduction

Crystallization is a key purification process which commonly follows the synthetic process where the active pharmaceutical ingredients (APIs) are produced [1,2,3]. However, obtaining the crystalline product with consistent/reproducible purity, size, and shape, and yet better-quality pharmaceutical manufacturing, is always challenging. Traditional batch crystallization is subject to inconsistent product quality, which translates into inconstant crystal size and shape distribution, which may affect the downstream processing such as filtration and drying. Mefenamic acid (MFA) exhibits high hydrophobicity and a propensity to stick to surfaces [4], which may cause significant fouling and blockage issues. Furthermore, MFA shows poor water solubility and slow dissolution rates, leading to significant challenges in achieving adequate bioavailability. Many efforts were made to produce MFA nanocrystals using antisolvent precipitation by accelerating the nucleation rate in the presence of additives [5]. However, these additives can further cause a problem for impurity removal. This paper describes the application of the Direct Nucleation Control (DNC) of a combined crystallization and wet-milling system to fine-tune the crystal size and shape distribution of MFA. The rotor–stator type of wet-milling has been widely used in the manufacturing process of API. This method is preferred as it can modify the crystal morphology more effectively compared to changing the solvents or seeking effective additives that may cause further purification or bioavailability issues. For years, wet-milling has been used for size reduction and continuous seed generators. However, an in-depth understanding of the mechanisms of wet milling is still very limited.
Firstly, the linear cooling crystallization of MFA in 2-butanol/heptane was studied with a batch crystallization platform which is subject to crystal agglomeration. Secondly, the DNC of a crystallization was studied under different targeted counts to investigate the crystal size and CSD. Finally, the DNC of a combined crystallization and wet-milling process was then experimentally investigated in both external and in situ wet-milling. The DNC with integrated crystallization and wet-milling aimed at reducing particle size and improving the size distribution. One of the key objectives was to investigate ways to use fluid shear imparted by the wet-mill to influence nucleation and breakage kinetics, thereby controlling the size and shape distribution of the crystal product. The wet-mill was either incorporated in a recycle loop with a one-stage MSMPR or used as an in situ devise that promotes nucleation and attrition operation. Process analytical technologies (PAT) were used to monitor the solute concentration, particle count, morphology, and size.

2. Material and Experimental Method

2.1. Materials

The API used in this study was MFA from Sigma-Aldrich (Gillingham, UK) with a purity of 98+%. The crystallization was intended to operate with an initial solvent mass ratio of 2-butanol (purity = 99%, ACROS Organics, Fisher Scientific) and heptane (purity > 99%, ACROS Organics, Fisher Scientific, Loughborough, UK) of 9:1. The batch linear cooling crystallization indicated that no other polymorphs crystallized in the 2-butanol/heptane system except the stable α form.

2.2. Experimental Setup

The experimental rig was built using a batch crystallizer (400 mL) combined with an in situ or external wet-mill device, as presented in Figure 1.
The process temperature in the jacketed glass crystallizer (Labtex Limited, Huddersfield, UK) was measured by a Pt100 thermocouple and controlled by a water bath (Ministat 230, Pilot ONE Huber, Offenburg, Germany). The PAT tools were applied to monitor the process variables. A Mettler Toledo (Redmond, WA, USA)’s G400 focused beam reflectance measurement (FBRM) probe (icFBRM software, version 4.3.377) was used to measure the crystal chord length distribution (CLD) and the number of crystals; the measurement interval was set as 15 seconds. A V918 particle vision and measurement (PVM) probe (Redmond, WA, USA) (Mettler Toledo PVM image acquisition software, version 8.3) with a high-resolution camera was used to record crystal images.
The particle size was analyzed by postprocessing the captured images in real-time using a built-in Blob image analysis method to detect particles and extract the particle size, aspect ratio (AR), and particle number (counts/s). The analyzed results were obtained from Lasentec PVM Stat Acquisition 6.0 software for frame edge lengths, which provides information about the size and AR from the particle projection. Solute concentration was monitored using the intensity of a specific peak in the spectra measured by an ATR–UV/vis probe (MCS651 UV, Carl Zeiss Microcopy, Gena, Germany). The UV/vis measurements’ interval was 5.3 seconds, and data were collected through LabVIEW based software. The real-time data were integrated and controlled through the LabVIEW based software, Crystallization Monitoring and Control (CryMOCO) [6]. The in situ and external devices consisted in a wet-mill of IKA ULTRA-TURRAX T25 basic (Werke, Germany) and wet-mill of IKA Magic Lab Module ULTRA-TURRAX (Werke, Germany). The external wet-mill used a medium type of generator that has two rows with a gap thickness of 0.5 mm. The working rotor diameters were 25 mm and 30 mm for the external wet-mill and 13.4 mm for the in situ mill. The gaps between rotor and stator were 0.5 mm and 0.25 mm for the external and in situ mill, individually. The slurry was circulated from the crystallizer to the external mill unit using a peristaltic pump (MasterFlex L/S, Thermo Fisher Scientific, Madison, WI, USA).

2.3. Experimental Methods

The initial solute concentration for experiments LC-1, LC-2, LC-3, DNC-1, and DNC-2 in Table 1 was prepared at a saturation temperature of 60 °C, based on a 0.35 kg 2-butanol/heptane solvent ratio at 90/10 for all experiments. The starting solution was heated up to 10 °C above the saturation temperature to ensure the dissolution of all MFA particles and cooled back to the saturation temperature without nucleation taking place. The effects of different cooling rates, ranging from 0.3 to 0.8 K/min (LC-1 to LC-3), as well as DNC crystallization with different final target counts of 5000 and 3000 #/s on crystal habit, were investigated in DNC-1 and DNC-2.
A lower initial concentration was chosen for the mill-aided DNC experiments, as milling significantly increases the number of crystals in the system, which may degrade the particle image capabilities due to overlapped particles, leading to inaccurate size and shape measurements. In this case, the experiments IMDNC-1 to IMDNC-2 and EMDNC were prepared at a saturation temperature of 40 °C. The wet-mill was conducted throughout the entire DNC heating and cooling cycles. The initial solution temperature was set at 40 °C for IMDNC-1, while it was forced to cool down to 10 °C at static conditions for IMDNC-2. The in situ milling speed was set to 2# gear of 9500 rpm and the overhead stirring speed at 400 rpm. Both IMDNC1 and IMDNC2 were carried out by setting the target counts at 25,000 #/s and the temperature range from 40 to 10 °C with the heating and cooling rates of 1 °C/min. EMDNC was carried out with a milling speed of 9600 rpm and 200 mL/min flow rate of the recirculation loop. EMDNC also used the same target counts of 25,000 #/s, but with a temperature ranging from 45 to 10 °C at heating and cooling rates of 1 °C/min. The list of the experiments conducted in this work is shown in Table 1 along with the crystallization methods, as well as a brief description of the target counts settings for DNC runs.

3. Results and Discussion

3.1. Batch Crystallization

To investigate the impact of the different cooling rates on batch cooling crystallization, linear cooling rates of 0.3, 0.5, and 0.8 K/min were used. Figure 2a shows that the nucleation temperatures are 50.1, 49.5, and 49 °C, respectively. This indicates a very narrow metastable zone with (MSDW) for MFA in the 2-butanol system. The temperature at which nucleation occurs decreases with the increased cooling rates, and, as a result, the system nucleates at a higher supersaturation at faster cooling rates. Consequently, higher nucleation rates and larger FBRM counts were obtained at faster cooling rates of 0.8 K/min. The PVM images of the crystals obtained with the different cooling rates depicted in Figure 2b–d show that elongated and slightly agglomerated crystals are formed in all cases. However, the crystal shape that evolves is plate-like, to form elongated crystals with the increased cooling rates. This indicated that the crystal growth of MFA is supersaturation dependent, relying on the depletion rate of the solute concentration [7]. These observations are consistent with the FBRM results, where higher FBRM counts were observed and more elongated crystal habits were obtained for faster cooling rates, such as in LC-3 (Figure 2d), while lower FBRM counts with more a plate-like crystal habit corresponded to slower cooling rates of LC-1 (Figure 2b).

3.2. DNC

The standard DNC approach aims at the control of the number of counts of particles to achieve a target value (or set-point) during the crystallization through the implementation of heating and cooling cycles. Thus, the number of counts is controlled by manipulating the operating conditions on either side of the solubility line to effectively control the dissolution of fine crystals during the heating cycles while promoting growth and preventing nucleation during the cooling cycles. The temperature, solute concentration, and the FBRM particle count profiles associated with the DNC experiment are shown in Figure 3. The DNC controller made the system converge to the target counts of 5000 #/s after two dissolution cycles, corresponding to a three-hour batch time (Figure 3a). A different set point of 3000 #/s was set to reduce the fine particles further.
Moreover, the fine particles were removed from the product, and there is evident growth between the heating and subsequent cooling stages after three dissolution cycles for approximately four hours (Figure 3c). Hence, more uniform CSD can be achieved by lowering the targeted counts from 5000 to 3000 #/s, by comparing the PVM images shown in Figure 3b and 3d. As mentioned in Section 3.1, the supersaturation dependence growth mechanism facilitates larger MFA crystals’ growth, hence the shape. Thus, the DNC approach enabled the production of less elongated crystals with fewer fines compared to the standard linear cooling of LC-1 to LC-3.

3.3. Mill-Aided DNC

The DNC operating cycles, crystal size, and habit were compared in both in situ and external mill-aided DNC crystallization. In situ mill-aided DNC seems to reach a pseudo-steady-state (Figure 4) where the FBRM counts fluctuate between 3000 and 4000 #/s, which is more than the preset target margin of 10%, 25,000 ± 2500 #/s.
Both in situ wet-mill aided DNC crystallizations can achieve a very uniform distribution with relatively small crystals. IMDNC-1 and IMDNC-2 with different starting temperatures show slight differences in their crystal habits. IMDNC-1 started with a lower initial temperature of 10 °C involved a higher-level of supersaturation, leading to a lower AR than that of IMDNC-2, which started with a saturation temperature of 40 °C.
The external wet-mill aided DNC run targeting 25,000 #/s is shown in Figure 5. The upper temperature of 45 °C, chosen here, is slightly higher than the saturated temperature due to the elevated nucleation rate triggered by the high dissipation rate of the external wet-mill, where the rate of dissolution cannot match the rate of nucleation and attrition. The preliminary experimental trials proved that 40 °C is not enough to carry on the DNC runs to regulate the final counts within the set of target counts at 25,000 ± 2500 #/s during the first heating cycle.
The external mill-aided DNC can further reduce particle sizes and narrow the particle size distribution, as shown in Figure 5b. The final product of EMDNC has fewer fines than those in IMDNC-1 or IMDNC-2 (Figure 4). Moreover, the AR of the crystals is smaller with a more equal particle shape. By comparing the two types of wet-mills, as mentioned in Section 2.2, the working rotor–stator diameter for the external mill is doubled, hence the mill tip speed, as it is proportional to the rotor diameter. Therefore, the mill-aided DNC can accomplish the size reduction and improve the shape of the MFA crystals by setting a higher target count.
Figure 6 shows the particle size distribution for all experiments. It can be seen that experiment DNC-2 displays the maximized particle size and EMDNC shows the minimized particle size. More fines have been observed showing bimodal size distribution for the experiments IMDNC-1 and IMDNC-2. In these two experiments, the pseudo-steady-state crystal products consisted of the first peak from crystal breakage and secondary nucleation due to milling and the second peak of larger crystals grown during the DNC cooling cycles. The experiment EMDNC presents relatively uniform CSD due to the larger dissipation rate, which enabled milling of most of the large-sized particles. Hence, large crystal size and enhanced AR can be obtained by optimizing unseeded DNC operating parameters, setting a different target count, and adjusting temperature ramp rates, while small particle size can be obtained by tuning the milling speed and solution flow-rate through the mill device.
The quantile particle sizes of D10, D50, and D90 and ARs were also estimated from image analysis of the PVM pictures, as summarized in Table 2. These data show that original cooling crystallization from 2-butanol/heptane solvents favors the crystal growth at major facets of (0 0 1), which is related to the degree of elongation [8,9]. MFA crystals grown from nonpolar solvents result in needle-shaped crystals, polar aprotic solvents result in elongated plate-shaped crystals, and polar protic solvents result in hexagonal cuboid-shaped crystals. In this work, 2-butanol has the nonpolar side of -C4H9 and a polar side of -OH. The linear cooling MFA product crystals LC-3 show mixtures of elongated and plate-like crystal shapes. For the standard DNC, the growth of facet (1 1 2) is favored and, hence, low AR crystals under a lower level of supersaturation for the rest of the cooling cycles after the first dissolution of DNC crystallization. Product crystals of DNC-1 and DNC-2 from standard DNC crystallization show larger particle sizes, D90, than the linear cooling crystalline products of LC-1 to LC-3, which is consistent with the final FRBM counts. All three mill-aided crystallization experiments have a similar medium size, but crystals from in situ mill-aided crystallization showed a broader distribution, with bimodal distribution as shown in Figure 6.

4. Conclusions

Much greater control over aspect ratio and crystal size distribution was demonstrated through DNC crystallization alone and in the presence of wet-milling. The wet-milling helped produce initial crystals with a reduced aspect ratio and minimal particle size by breaking the elongated crystals along their major axis. The supersaturation dependence growth rates at different facets of MFA crystals provide the potential to use DNC cycles to dissolve the fines and foster the directional growth to the minor axis of the remaining large crystals. In the last set of experiments, an in situ mill-aided DNC and external mill-aided DNC were implemented and compared against normal DNC (without wet-milling). The DNC significantly outperforms linear cooling crystallization, and the control of both aspect ratio and particle size distribution proved effective using DNC to maximize size and minimize fines. Furthermore, combined DNC with wet-milling proved effective to control particle size and particle size distribution. Hence, the dissolution rate of MFA can be greatly improved due to minimizing the particle size, which can be achieved by wet-milling and DNC, which reduces crystal fines.

Author Contributions

Conceptualization, W.L., C.D.R., B.B.; methodology, W.L., C.D.R., B.B.; software; validation, W.L., formal analysis, W.L.; investigation, W.L.; resources, C.D.R., B.B.; data curation, W.L.; writing—original draft preparation, W.L.; writing—review and editing, B.B. and C.D.R.; visualization, W.L.; supervision, B.B. and C.D.R.; project administration, B.B. and C.D.R.; funding acquisition, B.B. and C.D.R. All authors have read and agreed to the published version of the manuscript.

Funding

UK Engineering and Physical Research Council for the Future Continuous Manufacturing and Advanced Crystallisation Research Hub Grant EP/P006965/1.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors acknowledge funding from the UK Engineering and Physical Research Council for the Future Continuous Manufacturing and Advanced Crystallisation Research Hub Grant EP/P006965/1 (https://www.cmac.ac.uk/).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Su, Q.; Benyahia, B.; Nagy, Z.K.; Rielly, C.D. Mathematical Modeling, Design, and Optimization of a Multisegment Multiaddition Plug-Flow Crystallizer for Antisolvent Crystallizations. Org. Process Res. Dev. 2015, 19, 1859–1870. [Google Scholar] [CrossRef]
  2. Fysikopoulos, D.; Benyahia, B.; Borsos, A.; Nagy, Z.K.; Rielly, C.D. A framework for model reliability and estimability analysis of crystallization processes with multi-impurity multi-dimensional population balance models. Comput. Chem. Eng. 2019, 122, 275–292. [Google Scholar] [CrossRef]
  3. Zhou, L.; Su, M.; Benyahia, B.; Singh, A.; Barton, P.I.; Trout, B.L.; Myerson, A.S.; Braatz, R.D. Mathematical modeling and design of layer crystallization in a concentric annulus with and without recirculation. AIChE J. 2013, 59, 1308–1321. [Google Scholar] [CrossRef]
  4. Cesur, S.; Gokbel, S. Crystallization of mefenamic acid and polymorphs. Cryst. Res. Technol. 2008, 43, 720–728. [Google Scholar] [CrossRef]
  5. Bodnár, K.; Hudson, S.P.; Rasmuson, Å.C. Promotion of Mefenamic Acid Nucleation by a Surfactant Additive, Docusate Sodium. Cryst. Growth Des. 2019, 19, 519–603. [Google Scholar] [CrossRef]
  6. Borsos, Á.; Szilágyi, B.; Agachi, P.Ş.; Nagy, Z.K. Real-Time Image Processing Based Online Feedback Control System for Cooling Batch Crystallization. Org. Process Res. Dev. 2017, 21, 511–519. [Google Scholar] [CrossRef]
  7. Lovette, M.A.; Doherty, M.F. Predictive modeling of supersaturation-dependent crystal shapes. Cryst. Growth Des. 2012, 12, 656–669. [Google Scholar] [CrossRef]
  8. Shah, U.V.; Olusanmi, D.; Narang, A.S.; Hussain, M.A.; Gamble, J.F.; Tobyn, M.J.; Heng, J.Y.Y. Effect of crystal habits on the surface energy and cohesion of crystalline powders. Int. J. Pharm. 2014, 472, 140–147. [Google Scholar] [CrossRef] [PubMed]
  9. Hatcher, L.E.; Li, W.; Payne, P.; Benyahia, B.; Rielly, C.D.; Wilson, C.C. Tuning Morphology in Active Pharmaceutical Ingredients: Controlling the Crystal Habit of Lovastatin through Solvent Choice and Non-Size-Matched Polymer Additives. Cryst. Growth Des. 2020, 20, 5854–5862. [Google Scholar] [CrossRef]
Figure 1. Schematic of the experimental set-up for (a) the in situ wet-mill crystallization unit and (b) external wet-mill crystallization with recirculation.
Figure 1. Schematic of the experimental set-up for (a) the in situ wet-mill crystallization unit and (b) external wet-mill crystallization with recirculation.
Chemproc 09 00021 g001
Figure 2. (a) Temperature , FBRM counts , and concentration profiles for experiments LC-1 to LC-3 at different cooling rates of □ 0.3 K/min, ◯ 0.5 K/min, and △ 0.8 K/min; PVM image of final product crystals of (b) LC-1, (c) LC-2, and (d) LC-3.
Figure 2. (a) Temperature , FBRM counts , and concentration profiles for experiments LC-1 to LC-3 at different cooling rates of □ 0.3 K/min, ◯ 0.5 K/min, and △ 0.8 K/min; PVM image of final product crystals of (b) LC-1, (c) LC-2, and (d) LC-3.
Chemproc 09 00021 g002
Figure 3. (a) Temperature , FBRM counts , and concentration profiles for experiment DNC-1 with target counts of 5000 #/s and (b) PVM image of final product crystals of DNC-1. (c) Temperature , FBRM counts , and concentration profiles for experiment DNC-2 with target counts of 3000 #/s, and (d) PVM image of final product crystals for DNC-2.
Figure 3. (a) Temperature , FBRM counts , and concentration profiles for experiment DNC-1 with target counts of 5000 #/s and (b) PVM image of final product crystals of DNC-1. (c) Temperature , FBRM counts , and concentration profiles for experiment DNC-2 with target counts of 3000 #/s, and (d) PVM image of final product crystals for DNC-2.
Chemproc 09 00021 g003
Figure 4. (a) Temperature , FBRM counts , and concentration profiles for experiment IMDNC-1, and (b) PVM image of final product crystals of IMDNC-1. (c) Temperature , FBRM counts , and concentration profiles for experiment IMDNC-2, and (d) PVM image of final product crystals for IMDNC-2.
Figure 4. (a) Temperature , FBRM counts , and concentration profiles for experiment IMDNC-1, and (b) PVM image of final product crystals of IMDNC-1. (c) Temperature , FBRM counts , and concentration profiles for experiment IMDNC-2, and (d) PVM image of final product crystals for IMDNC-2.
Chemproc 09 00021 g004
Figure 5. (a) Temperature , FBRM counts , and concentration profiles for experiment EMDNC-1 with final target counts of 25,000 #/s, and (b) PVM image of final product crystals of EMDNC.
Figure 5. (a) Temperature , FBRM counts , and concentration profiles for experiment EMDNC-1 with final target counts of 25,000 #/s, and (b) PVM image of final product crystals of EMDNC.
Chemproc 09 00021 g005
Figure 6. Particle size distribution for experiments LC-1, LC-2, LC3, DNC-1, DNC-2, IMDNC-1, IMDNC-2, EMDNC.
Figure 6. Particle size distribution for experiments LC-1, LC-2, LC3, DNC-1, DNC-2, IMDNC-1, IMDNC-2, EMDNC.
Chemproc 09 00021 g006
Table 1. Process conditions for the crystallization experiments.
Table 1. Process conditions for the crystallization experiments.
ExperimentCrystallization MethodOperating T (°C)Ramp Rate (K/min)Target Counts (#/s)
LC-1Linear cooling60–200.3-
LC-2Linear cooling60–200.5-
LC-3Linear cooling60–200.8-
DNC-1DNC60–2015000
DNC-2DNC60–2013000
IMDNC-1In situ Milling DNC40–10; initial T = 10125,000
IMDNC-2In situ Milling DNC40–10; initial T = 40125,000
EMDNCExternal Milling DNC45–10; initial T = 45125,000
Table 2. Process conditions of crystallization experiments.
Table 2. Process conditions of crystallization experiments.
ExperimentLC-1LC-2LC-3DNC-1DNC-2IMDNC-1IMDNC-2EMDNC
D10 (μm)7234307187232422
D50 (μm)1778478213245505441
D90 (μm)32016214843543616117165
AR2.32.83.02.72.51.92.01.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, W.; Rielly, C.D.; Benyahia, B. Targeting Desired Particle Size for Improved Dissolution and Manufacturability of Mefenamic Acid. Chem. Proc. 2022, 9, 21. https://doi.org/10.3390/IOCC_2022-12165

AMA Style

Li W, Rielly CD, Benyahia B. Targeting Desired Particle Size for Improved Dissolution and Manufacturability of Mefenamic Acid. Chemistry Proceedings. 2022; 9(1):21. https://doi.org/10.3390/IOCC_2022-12165

Chicago/Turabian Style

Li, Wei, Chris D. Rielly, and Brahim Benyahia. 2022. "Targeting Desired Particle Size for Improved Dissolution and Manufacturability of Mefenamic Acid" Chemistry Proceedings 9, no. 1: 21. https://doi.org/10.3390/IOCC_2022-12165

APA Style

Li, W., Rielly, C. D., & Benyahia, B. (2022). Targeting Desired Particle Size for Improved Dissolution and Manufacturability of Mefenamic Acid. Chemistry Proceedings, 9(1), 21. https://doi.org/10.3390/IOCC_2022-12165

Article Metrics

Back to TopTop