Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fractionalization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An optimization-based approach to parameter learning for fractional type nonlocal models

Nonlocal operators of fractional type are a popular modeling choice for applications that do not adhere to classical diffusive behavior; however, one major challenge in nonlocal simulations is the selection of model parameters. In this work we propose an optimization-based approach to parameter identification for fractional models with an optional truncation radius. We formulate the inference problem as an optimal control problem where the objective is to minimize the discrepancy between observed data and an approximate solution of the model, and the control variables are the fractional order and the truncation length. For the numerical solution of the minimization problem we propose a gradient-based approach, where we enhance the numerical performance by an approximation of the bilinear form of the state equation and its derivative with respect to the fractional order. Several numerical tests in one and two dimensions illustrate the theoretical results and show the robustness and applicability of our method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

An index to characterize gas‐solid and solid‐solid mixing from average volume fraction fields

Abstract A mixing index based on solid volume fraction fields is developed for gas‐solid flows. Conventional mixing indices are based on particle realizations of granular mixing and are applicable to experimental data or discrete element method simulations. However, these indices cannot be used as‐is for multifluid models, and an index for characterizing mixing in gas‐solid flows from continuous fields is needed. The performance of the new mixing index is tested in two applications. The first is a 3D simulation of the mixing of biomass and sand in a fluidized bed reactor, and the second is a 2D simulation of binary particle segregation in a fluidized bed. The simulations are performed using OpenFOAM®. The mixing index is used to quantify gas‐solid mixing using solid volume fractions and solid‐solid mixing using solid fractions. The formulation of conventional mixing indices is extended to be used with solid volume fractions fields, and methods for performance improvement are presented.

09 BIOMASS FUELS↗
DOE OSTI2025

Measurements of the branching fractions of ${\Xi }_{c}^{+}\to {\Sigma }^{+}{K}_{S}^{0}$, ${\Xi }_{c}^{+}\to {\Xi }^{0}{\pi }^{+}$, and ${\Xi }_{c}^{+}\to {\Xi }^{0}{K}+$ at Belle and Belle II

Using 983.0 fb −1 and 427.9 fb −1 data samples collected with the Belle and Belle II detectors at the KEKB and SuperKEKB asymmetric energy e + e − colliders, respectively, we present studies of the Cabibbo-favored ${\Xi }_{c}^{+}$ decays ${\Xi }_{c}^{+}\to {\Sigma }^{+}{K}_{S}^{0}$ and ${\Xi }_{c}^{+}\to {\Xi }^{0}{\pi }^{+}$, and the singly Cabibbo-suppressed decay ${\Xi }_{c}^{+}\to {\Xi }^{0}{K}^{+}$. The ratios of branching fractions of ${\Xi }_{c}^{+}\to {\Sigma }^{+}{K}_{S}^{0}$ and ${\Xi }_{c}^{+}\to {\Xi }^{0}{K}^{+}$ relative to that of ${\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+}$ are measured for the first time, while the ratio $\mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{0}{\pi }^{+})/\mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+})$ is also determined and improved by an order of magnitude in precision. The measured branching fraction ratios are $\begin{array}{c}\frac{\mathcal{B}\left({\Xi }_{c}^{+}\to {\Sigma }^{+}{K}_{S}^{0}\right)}{\mathcal{B}\left({\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+}\right)}=0.067\pm 0.007\pm 0.003,\\ \frac{\mathcal{B}\left({\Xi }_{c}^{+}\to {\Xi }^{0}{\pi }^{+}\right)}{\mathcal{B}\left({\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+}\right)}=0.251\pm 0.005\pm 0.010,\\ \frac{\mathcal{B}\left({\Xi }_{c}^{+}\to {\Xi }^{0}{K}^{+}\right)}{\mathcal{B}\left({\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+}\right)}=0.017\pm 0.003\pm 0.001.\end{array}$ Additionally, the ratio $\mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{0}{K}^{+})/\mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{0}{\pi }^{+})$ is measured to be 0.068 ± 0.010 ± 0.004. Here, the first and second uncertainties are statistical and systematic, respectively. Multiplying the ratios by the branching fraction of the normalization mode, $\mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+})=(2.9\pm 1.3)%$, we obtain the following absolute branching fractions $\begin{array}{c}\mathcal{B}({\Xi }_{c}^{+}\to {\Sigma }^{+}{K}_{S}^{0})=(0.194\pm 0.021\pm 0.009\pm 0.087)\text{%},\\ \mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{0}{\pi }^{+})=(0.728\pm 0.014\pm 0.027\pm 0.326)\text{%},\\ \mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{0}{K}^{+})=(0.049\pm 0.007\pm 0.003\pm 0.022)\text{%},\end{array}$ where the third uncertainties are from $\mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+})$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Solidification texture, variant selection, and phase fraction in a spot-melt electron-beam powder bed fusion processed Ti-6Al-4V

The effect of a spot melt strategy on the solidification texture, variant selection, phase fraction, and their variations along the build height was investigated in comparison to a conventional linear melt strategy. Ti-6Al-4V alloy samples were manufactured by an electron-beam powder bed fusion process using two different melt strategies and characterized using high-energy synchrotron x-ray diffraction. The linear and spot melt cases resulted in β<100> cube and fiber textures, respectively, with the cube/fiber axis nearly parallel to the build direction. The α phase exhibited Burgers orientation relationship (BOR) with the parent β phase for both melt strategies. Although all of the key texture components and the BOR between the α and β phases were consistently observed in both cases along the build height, there were several measurable differences between the two. At the top of the build height, the β<100> build texture intensity was stronger in the linear case. However, the texture intensity in the linear case decreased gradually from the top to near the substrate, whereas the spot melt case did not show this trend. Moreover, the BOR was examined along the build height and the results showed that there is a distinct variant selection in the spot melt case, unlike the linear case. Specifically, the planar variants sharing ($\bar101$) β || (0002) α with [$\bar2\bar2\bar2$]β and [$2\bar22$]β were the most prominent in the spot melt, but no preference was identified in the directional variants. The β phase fraction was slightly lower in the spot melt case compared to the linear case. In both the cases, the β phase fraction decreased moving towards the substrate but remained significantly higher than that of the powder feedstock. Overall, the novel spot melt strategy produced a more homogeneous microstructure in terms of both the phase fraction and texture across the build height.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Linking Manganese Fractions and Organic Carbon in Soils of Contrasting Land Use Systems

Manganese (Mn) is an essential micronutrient that influences carbon (C) cycling by binding or oxidizing soil organic matter. Mn fractions in soil and their plant availability depend largely on soil pH, which is commonly altered by agricultural practices. Fractions of Mn in soils range from readily available (e.g., bioavailable Mn, exchangeable Mn) to stabilized (e.g., Mn-oxide minerals, Mn contained in silicates). Furthermore, the distribution of soil Mn with depth was assessed in contrasting land use systems (organic agriculture, conventional agriculture, and unmanaged forest) using a sequential extraction method that targets Mn fractions ranging in bioavailability. Both agricultural sites had similar amounts of total Mn but had available Mn (1–7% of total Mn) lower than that of the unmanaged forested site (15% of total Mn). Manganese and organic C were generally positively correlated, but this relationship depended on soil depth, Mn fractions, and concentrations, while land management change had little influence.

Crops↗
DOE OSTI2021

Improved Li + Transport in Polyacetal Electrolytes: Conductivity and Current Fraction in a Series of Polymers

Polymer electrolytes mitigate safety concerns surrounding flammable liquid electrolytes in lithium-ion batteries. Poly(ethylene oxide) (PEO) electrolytes demonstrate viable conductivity values (~1 × 10 -3 S/cm) at elevated temperatures (>70 °C) but a relatively low Li + current fraction (≤0.2) because strong Li + coordination inhibits cation mobility. We have developed a series of polyacetal electrolytes by systematically varying methylene oxide (MO) and ethylene oxide (EO) units in the polymer backbone. These materials maintain high oxygen-to-carbon ratios like PEO but offer improved ion transport, revealing trends of decreasing conductivity and increasing current fraction with respect to polymer composition. In particular, the increasing current fraction measured via the Bruce-Vincent method suggests that MO units improve Li+ mobility relative to anion mobility. Here, we calculate an overall efficacy (product of conductivity and current fraction) for each polymer/salt composition and identify two polymers - P(EO-MO) and P(EO-2MO) - that outperform PEO at high and low salt concentrations, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Oxidative Catalytic Fractionation of Lignocellulosic Biomass under Non-alkaline Conditions

Biomass pretreatment methods are commonly used to isolate carbohydrates from biomass, but they often lead to modification, degradation, and/or low yields of lignin. Catalytic fractionation approaches provide a possible solution to these challenges by separating the polymeric sugar and lignin fractions in the presence of a catalyst that promotes cleavage of the lignin into aromatic monomers. In this work, we demonstrate an oxidative fractionation method conducted in the presence of a heterogeneous non-precious-metal Co-N-C catalyst and O 2 in acetone as the solvent. The process affords a 15 wt% yield of phenolic products bearing aldehydes (vanillin, syringaldehyde) and carboxylic acids (p-hydroxybenzoic acid, vanillic acid, syringic acid), complementing the alkylated phenols obtained from existing reductive catalytic fractionation methods. The oxygenated aromatics derived from this process have appealing features for use in polymer synthesis and/or biological funneling to value-added products, and the non-alkaline conditions associated with this process support preservation of the cellulose, which remains insoluble at reaction conditions and is recovered as a solid.

09 BIOMASS FUELS↗
DOE OSTI2022

Fractional correlated insulating states at one-third filled magic angle twisted bilayer graphene

Abstract The observation of superconductivity and correlated insulating states in twisted bilayer graphene has motivated much theoretical progress at integer fillings. However, little attention has been given to fractional fillings. Here we show that the three-peak structure of Wannier orbitals, dictated by the symmetry and topology of flat bands, facilitates the emergence of a state we name a “fractional correlated insulator” at commensurate fractional filling of ν = n ± 1/3. Specifically for the filling of 1/3 electrons per moiré unit cell, we show that short-range interactions lead to an extensive entropy due to the “breathing” degree of freedom of an irregular honeycomb lattice that emerges through defect lines. The leading further-range interaction lifts this degeneracy and selects a ferromagnetic nematic state that breaks AB/BA sublattice symmetry. The proposed fractional correlated insulating state might underlie the suppression of superconductivity at ν = 2 − 1/3 filling observed in ref. 1 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Generalized vs. fractional: a comparative analysis of Maxwell models applied to entangled polymer solutions

Fractional viscoelastic models provide an excellent description of rheological data for polymer systems with power-law behaviour. However, the physical interpretation of their model parameters, which carry fractional units of time, often remains elusive. We show that for poly(ethylene oxide) (PEO) solutions, the fractional Maxwell model (FMM) requires much fewer model parameters than the classical generalized Maxwell model for a good description of the data and that it can be applied consistently to solutions with varying degrees of viscoelasticity. Despite their fractional units, the parameters exhibit scaling laws similar to classical parameters as a function of polymer concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Integrating catalytic fractionation and microbial funneling to produce 2-pyrone-4,6-dicarboxylic acid and ethanol

Replacing biorefinery designs that use stepwise fractionation, depolymerization, and conversion processes with designed tandem process steps can greatly reduce a biorefinery's operating costs and environmental impact. Reductive Catalytic Fractionation (RCF) is a highly efficient lignin-first approach that combines biomass fractionation and lignin depolymerization to generate a hydrogenolysis oil and pulp. The oil is composed of a complex mixture of phenolic monomers, dimers, and oligomers. Intergrating this chemical deconstruction process with microbial funneling of phenolics can simplify the product mixture and make high-value products. We applied RCF to poplar biomass in a biomass-to-bioproduct processing chain in which the phenolics were funneled to 2-pyrone-4,6-dicarboxylic acid (PDC) by an engineered strain of Novosphingobium aromaticivorans DSM12444. The pulp was enzymatically digested and the glucose and xylose was funneled to ethanol by an engineered strain of Saccharomyces cerevisiae GLBRCY945. By combining biomass fractionation and lignin depolymerization we removed a costly processing step that directly translated into a 29% reduction in the minimum selling price of PDC. This work combines experimentation with process modeling of an integrated biorefinery design to show how utilizing tandem process steps can significantly reduce operating expenses and environmental impact of upgrading lignocellulosic biomass to a portfolio of high value products.

Sener, Canan [Great Lakes Bioenergy Research Cente↗
DOE OSTI2020

Shear-induced lift force on spheres in a viscous linear shear flow at finite volume fractions

Several studies have shown a significant increase in drag on a distribution of solid spherical particles within a fluid with increasing particle volume fraction. As a result, many empirical drag laws accounting for the dependence on the Reynolds number and volume fraction can be found in the literature. This study investigates the possibility of a similar effect of the particle volume fraction on the mean hydrodynamic lift force on randomly distributed spherical particles in a linear shear flow. Particle-resolved direct numerical simulations are performed to evaluate the mean lift force, and the results are compared with the case of an isolated particle in a linear shear flow for the same Reynolds number and shear rate. The mean lift force acting on the particles appears to remain nearly the same as that on an isolated particle. However, due to the influence of neighboring particles, there is a substantial force variation in transverse directions on each individual particle, whose magnitude is comparable to the mean drag force. The distribution of drag force in a linear shear flow is shown to be nearly the same as in a uniform flow at the same volume fraction and Reynolds number. A simple stochastic model based on a Gaussian distribution is presented for the lift force variation, and its performance is compared to the prediction of the deterministic pairwise interaction extended point-particle model.

42 ENGINEERING↗
DOE OSTI2023

Carbon isotope fractionation by an ancestral rubisco suggests that biological proxies for CO 2 through geologic time should be reevaluated

The history of Earth’s carbon cycle reflects trends in atmospheric composition convolved with the evolution of photosynthesis. Fortunately, key parts of the carbon cycle have been recorded in the carbon isotope ratios of sedimentary rocks. The dominant model used to interpret this record as a proxy for ancient atmospheric CO 2 is based on carbon isotope fractionations of modern photoautotrophs, and longstanding questions remain about how their evolution might have impacted the record. Therefore, we measured both biomass (ε p ) and enzymatic (ε Rubisco ) carbon isotope fractionations of a cyanobacterial strain ( Synechococcus elongatus PCC 7942) solely expressing a putative ancestral Form 1B rubisco dating to ≫1 Ga. This strain, nicknamed ANC, grows in ambient pCO 2 and displays larger ε p values than WT, despite having a much smaller ε Rubisco (17.23 ± 0.61‰ vs. 25.18 ± 0.31‰, respectively). Surprisingly, ANC ε p exceeded ANC ε Rubisco in all conditions tested, contradicting prevailing models of cyanobacterial carbon isotope fractionation. Such models can be rectified by introducing additional isotopic fractionation associated with powered inorganic carbon uptake mechanisms present in Cyanobacteria, but this amendment hinders the ability to accurately estimate historical pCO 2 from geological data. Understanding the evolution of rubisco and the CO 2 concentrating mechanism is therefore critical for interpreting the carbon isotope record, and fluctuations in the record may reflect the evolving efficiency of carbon fixing metabolisms in addition to changes in atmospheric CO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Constraining the local pulsar population with the cosmic-ray positron fraction

Observations of the TeV halos associated with nearby pulsars indicate that these objects inject significant fluxes of very high-energy electron-positrons pairs into the interstellar medium (ISM), thereby likely providing the dominant contribution to the cosmic-ray positron flux. In this paper, we use the cosmic-ray positron fraction as measured by the AMS-02 Collaboration to constrain the characteristics of the local pulsar population. For reasonable model parameters, we find that we can obtain good agreement with the measured positron fraction up to energies of E e ~ 300 GeV. At higher energies, the positron fraction is dominated by a small number of pulsars, making it difficult to reliably predict the shape of the expected positron fraction. The low-energy positron spectrum supports the conclusion that pulsars typically transfer approximately η ~ 5–20% of their total spindown power into the production of very high-energy electron-positron pairs, producing a spectrum of such particles with a hard spectral index, α ~ 1.5–1.7. Such pulsars typically spindown on a timescale on the order of τ ~ 10 4 years. Our best fits were obtained for models in which the radio and gamma-ray beams from pulsars are detectable to 28% and 62% of surrounding observers, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Polarization fraction of Planck Galactic cold clumps and forecasts for the Simons Observatory

ABSTRACT We measure the mean-squared polarization fraction of a sample of 6282 Galactic cold clumps at 353 GHz, consisting of Planck Galactic cold clump (PGCC) catalogue category 1 objects [flux densities measured with signal-to-noise ratio (S/N) > 4]. At 353 GHz, we find the mean-squared polarization fraction, which we define as the mean-squared polarization divided by the mean-squared intensity, to be (4.79 ± 0.44) × 10−4 equation to an $11\, \sigma$ detection of polarization. We test if the polarization fraction depends on the clumps’ physical properties, including flux density, luminosity, Galactic latitude, and physical distance. We see a trend towards increasing polarization fraction with increasing Galactic latitude, but find no evidence that polarization depends on the other tested properties. The Simons Observatory, with angular resolution of order 1 arcmin and noise levels between 22 and $54\, \mu$K−arcmin at high frequencies, will substantially enhance our ability to determine the magnetic field structure in Galactic cold clumps. At $\ge 5\, \sigma$ significance, we predict the Simons Observatory will detect at least ∼12 000 cold clumps in intensity and ∼430 cold clumps in polarization. This number of polarization detections would represent a two orders of magnitude increase over the current Planck results. We also release software that can be used to mask these Galactic cold clumps in other analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Superfluid condensate fraction and pairing wave function of the unitary Fermi gas

The unitary Fermi gas is a many-body system of two-component fermions with zero-range interactions tuned to infinite scattering length. Despite much activity and interest in unitary Fermi gas and its universal properties, there have been great difficulties in performing accurate calculations of the superfluid condensate fraction and pairing wave function. In this paper, we present auxiliary-field lattice Monte Carlo simulations using a lattice interaction which accelerates the approach to the continuum limit, thereby allowing for robust calculations of these difficult observables. As a benchmark test, we compute the ground-state energy of 33 spin-up and 33 spin-down particles. As a fraction of the free Fermi gas energy $E_{\text{FG}}$, we find $E_0/E_{\text{FG}}$ = 0.369(2), 0.372(2), using two different definitions of the finite-system energy ratio, in agreement with the latest theoretical and experimental results. We then determine the condensate fraction by measuring off-diagonal long-range order in the two-body density matrix. We find that the fraction of condensed pairs is α = 0.43(2). Further, we also extract the pairing wave function and find the pair correlation length to be $ζ_pk_F$ = 1.8(3)ℏ, where $k_F$ is the Fermi momentum. Provided that the simulations can be performed without severe sign oscillations, the methods we present here can be applied to superfluid neutron matter as well as more exotic $\textit{P}$-wave and $\textit{D}$-wave superfluids.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Edge-driven transition between extended quantum anomalous Hall crystal and fractional Chern insulator in rhombohedral graphene multilayers

Fractional Chern insulators (FCI) with fractionally quantized Hall conductance at fractional fillings and an extended quantum anomalous Hall (EQAH) crystal with an integer quantized Hall conductance over an extended region of doping were recently observed in pentalayer graphene. One particularly puzzling observation is the transition between the EQAH and FCI regimes, driven either by temperature or electrical current. Here, in this study, we propose a scenario to understand these transitions based on the topologically protected gapless edge modes that are present in both the FCI and EQAH phases and should be most relevant at temperature scales below the energy gap. Our consideration is based on the simple assumption that the edge velocity in FCI is smaller than that in EQAHE and thus contributes to a higher entropy. We further argue that domains with opposite fractionally quantized Hall conductance are ubiquitous in the devices due to disorder, which gives rise to a network of edge modes. The velocity of the edge modes between domains is further reduced due to edge reconstruction. The edge velocity can also be reduced by current when the occupation of the edge mode approaches the gap edge. The edge entropy therefore drives the transition from EQAH to FCI either by temperature or current at a nonzero temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Transverse momentum dependent feed-down fractions for bottomonium production

We extract transverse-momentum-dependent feed-down fractions for bottomonium production using a data-driven approach. We use data published by the ATLAS, CMS, and LHCb Collaborations for √𝑠 =7 TeV proton-proton collisions. Based on this collected data, we produce fits to the differential cross sections for the production of both 𝑆- and 𝑃-wave bottomonium states. Combining these fits with branching ratios for excited state decays from the Particle Data Group, we compute the feed-down fractions for both the ϒ⁡(1⁢𝑆) and ϒ⁡(2⁢𝑆) as a function of transverse momentum. Our results indicate a strong dependence on transverse momentum, which is consistent with prior extractions of the feed-down fractions. When evaluated at the average momentum of the states, we find that approximately 75% of ϒ⁡(1⁢𝑆) and ϒ⁡(2⁢𝑆) states are produced directly. Our results for the transverse-momentum-dependent feed-down fractions are provided in tabulated form so that they can be used by other research groups.

Astronomy & Astrophysics↗
DOE OSTI2024

Measurement of the absolute branching fraction of the three-body decay Λ c + → Ξ 0 K + π 0 and search for Λ c + → n K + π 0 , Σ 0 K + π 0 , and Λ K + π 0

The Cabbibo-favored decay Λ c + → Ξ 0 K + π 0 is studied for the first time using 6.1 fb − 1 of e + e − collision data at center-of-mass energies between 4.600 and 4.840 GeV, collected with the BESIII detector at the BEPCII collider. With a double-tag method, the branching fraction of the three-body decay Λ c + → Ξ 0 K + π 0 is measured to be ( 7.79 ± 1.46 ± 0.95 ) × 10 − 3 , where the first and second uncertainties are statistical and systematic, respectively. The branching fraction of the two-body decay Λ c + → Ξ ( 1530 ) 0 K + is ( 5.99 ± 1.04 ± 0.32 ) × 10 − 3 , which is consistent with the previous result of ( 5.02 ± 0.99 ± 0.31 ) × 10 − 3 . In addition, the upper limit on the branching fraction of the doubly Cabbibo-suppressed decay Λ c + → n K + π 0 is 7.1 × 10 − 4 at the 90% confidence level. The upper limits on the branching fractions of Λ c + → Σ 0 K + π 0 and Λ K + π 0 are also determined to be 1.8 × 10 − 3 and 2.0 × 10 − 3 , respectively. Published by the American Physical Society 2024

Ablikim, M.↗