Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactive model”

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 613 records · Page 34

Evaluation of Joint Cyber/Safety Risk in Nuclear Power Systems

This report presents an analysis of the Emergency Core Cooling System (ECCS) for a generic Boiling Water Reactor (BWR)-4 NPP. The Electric Power Research Institute (EPRI) developed Hazards and Consequences Analysis for Digital Systems (HAZCADS) process is applied to the ECCS and its subsystems to identify unsafe control actions (UCAs) which act as possible cyber events of concern. The analysis is performed for two design basis events: Small-break Loss of Coolant Accident (SLOCA) and general transients (TRANS), such as unintended reactor trip. In previous work, HAZCADS UCAs were combined with other cyber-attack analysis to develop a risk-informed approach; however, this was for a single system. This report explores advanced systems engineering modeling approaches to model the interactions between digital assets across multiple systems which may be targeted by cyber adversaries. The complex and interdependent design of digital systems has the potential to introduce emergent cyber properties that are generally not covered by hazard analyses nor formal nuclear Probabilistic Risk Assessment (PRA). The R&D and supporting analysis presented here explores approaches to predict and manage how interdependent system properties effect risk. To show the potential impact of a successful cyber-attack to formal PRA event tree probabilities, HAZCADS analysis was also used. HAZCADS was also used to model the automatic depressurization system (ADS) automatic actuation. This analysis extended to an integrated system analysis for common-cause failure (CCF). In this aspect, the HAZCADS analysis continued by analyzing plant design details for system connectivity in support of critical plant functions. A dependency matrix was developed to depict the integrated functionality of the interconnected systems. Areas of potential CCF are indicated. Future work could include adversary attack development to show how CCF could be caused, resulting in PRA events. Across the multiple systems that comprise the ECCS, the analysis shows that the change in such probabilities was very different between systems. This indicates that some systems have a larger potential risk impact from successful cyber-attack or digital failure, which indicates a need for these systems to have a higher priority for design and defensive measures. Furthermore, we were able to establish that a risk analysis using any arbitrary threat model establishes an ordering of components with regard to cyber-risk. This ordering can be used to influence the overall system design with an eye to lowering risk, or as a way to understand real-time risk to operational systems based on a current threat landscape. Expert knowledge of both the analysis process and the system being analyzed is required to perform a HAZCADS analysis. The need for a tiered risk analysis is demonstrated by the results of this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Rotational transition, domain formation, dislocations, and defects in vortex systems with combined sixfold and twelvefold anisotropic interactions

We introduce a phenomenological model for a pairwise repulsive interaction potential of vortices in a type-II superconductor, consisting of superimposed sixfold and twelvefold anisotropies. Using numerical simulations we study how the vortex lattice configuration varies as the magnitudes of the two anisotropic interaction terms change. A triangular lattice appears for all values, and rotates through 30 ° as the ratio of the sixfold and twlevefold anisotropy amplitudes is varied, in agreement with experimental results. The transition causes the vortex lattice to split into domains that have rotated clockwise or counterclockwise, with grain boundaries that are “decorated” by dislocations consisting of fivefold and sevenfold coordinated vortices. We also find intradomain dislocations and defects, and characterize them in terms of their energy cost. We discuss how this model could be generalized to other particle-based systems with anisotropic interactions, such as colloids, and consider the limit of very large anisotropy where it is possible to create cluster crystal states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Non-Hermitian approach for quantum plasmonics

Herein, we examine the limits of applicability of a simple non-Hermitian model for exciton/plasmon interactions in the presence of dissipation and dephasing. The model can be used as an alternative to the more complete Lindblad density matrix approach and is computationally and conceptually simpler. We find that optical spectra in the linear regime can be adequately described by this approach. The model can fail, however, under continuous optical driving in some circumstances. In the case of two quantum dots or excitons interacting with a plasmon, the model can also describe coherences and entanglement qualitatively when both dissipation and dephasing are present and quantitatively in the limit with no dephasing. The approach, within a single excitation manifold, is also applied to assess the role of disorder for 50 quantum dots interacting with a plasmon, where we find that, on average, large enough disorder can help stabilize the ensemble average of the open quantum system toward a dark quasi-steady-state much faster than without disorder. Furthermore, while such single excitation manifold calculations in this size limit can readily be done with either the non-Hermitian or Lindblad forms, as one goes to larger Hilbert space sizes, the computational and storage advantages of the non-Hermitian approach can become more useful. Published under license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bacterial exometabolites influence Chlamydomonas cell cycle and double algal productivity

Algal-bacterial interactions provide clues to algal physiology, but mutualistic interactions are complicated by dynamic exchange. In this study, we characterized the response of Chlamydomonas reinhardtii to the presence of a putative alga-benefitting commensal bacterium ( Arthrobacter strain ‘P2b’). Co-cultivation promoted chlorophyll content, biomass, average cell size, and number of dividing cells, relative to axenic cultures. Addition of bacterial spent medium (whole, size-fractionated and heat-treated) had similar effects, indicating P2b does not require algal interaction to promote growth. Nutrients and pH were excluded as putative effectors, collectively indicating a commensal interaction mediated by Arthrobacter -released small exometabolite(s). Proteogenomic comparison revealed similar response to co-cultivation and spent media, including differential cell cycle regulation, extensive downregulation of flagellar genes and histones, carbonic anhydrase and RubisCO downregulation, upregulation of some chlorophyll, amino acid and carbohydrate biosynthesis genes, and changes to redox and Fe homeostasis. Further, Arthrobacter protein expression indicated some highly expressed putative secondary metabolites. Together, these results revealed that low molecular weight bacterial metabolites can elicit major physiological changes in algal cell cycle regulation, perhaps through a more productive G1 phase, that lead to substantial increases in photosynthetically-produced biomass. This work illustrates that model commensal interactions can be used to shed light on algal response to stimulating bacteria.

59 BASIC BIOLOGICAL SCIENCES↗

β -delayed neutron emission of Mn 64 , Cr 62 , and Fe 65

The β-decay properties of nuclei near the second nuclear “island of inversion” around neutron rich nuclei with neutron number 40 are important tests of nuclear structure models and interactions. In particular, the β-delayed neutron emission branch (Ρ n ), is useful for investigating β-strength and neutron-γ competition above the neutron separation energies of the daughter nuclei. We report new constraints for Ρ n values for three nuclei in the region: 62 Cr (Ρ n <1%), 64 Mn (Ρ n = 1.5⁢(6)%), and 65 Fe (Ρ n < 1%), measured with the Neutron Emission Ratio Observer (NERO) neutron long counter system and the Beta Counting Station (BCS) at the National Superconducting Cyclotron Laboratory (NSCL). Our results resolve the large discrepancy between previous direct and indirect measurements for 64 Mn and confirm the predictions of global theoretical models when a statistical treatment of the γ and neutron decays of the daughter states is included. Here we also obtain improved half-lives for 62 Cr [206(5) ms] and the short-lived isomer in the 62 Fe daughter [112(7) ms] from β-delayed γ emission data obtained in the same experiment with the Summing NaI (SuN) total absorption spectrometer. Finally, we use γ emission data to obtain a new upper limit for the 62 Cr β-decay population of the long-lived isomeric state in 62 Mn.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Towards a predictive, physics-based friction model for the dynamics of jointed structures

Bolted connections are ubiquitous in mechanical designs and pose a significant challenge to understanding and predicting the vibration response of assembled structures. The present paper develops a physics-based rough contact model of the frictional interactions within a joint. This model sums over the probable interactions of asperities – defined as locally maximum surface features – to determine the contact forces. Here, the tangential contact forces vary smoothly between sticking and slipping and allow the model to better capture the qualitative trends of experimental amplitude dependent frequency and damping than previous studies. Furthermore, the novel model is generalized to allow for arbitrarily varying normal pressure including potential separation to better represent the interfacial dynamics. This includes developing a new, computationally tractable approximation to the analytical Mindlin partial slip solution for tangential loading of contacting spheres. Further, the results highlight the importance of accurately characterizing the as-built topology of the interface, the plastic behavior of the contacting asperities, the relevant length scale of asperities, and the eccentricity of asperities. A predictive friction coefficient based on plasticity provides a poor match to experiments, so fitting the friction coefficient is also considered. Numerical results are compared to experiments on the Brake-Reuß Beam to assess the predictive potential of the models. While blind predictions over-predict the slip limit, the current model presents a significant improvement in physics-based modeling and highlights areas for ongoing research.

42 ENGINEERING↗

Comparison of temperature adaptive calibration methods for laser induced fluorescence based fuel-in-oil instrument

Fuel dilution of engine oil (or fuel-in-oil, FiO) is an important issue as multiple and late-cycle fuel injection, integral to many combustion efficiency and emissions improvements (e.g., downsized boosted gasoline engines and catalyst thermal management) increases FiO rate. In addition to causing general wear and corrosion in engine due to decreased oil viscosity and pH buffering, FiO is also believed to cause destructive low speed pre-ignition (or super knock) in boosted SI engines. To understand the effects of engine operating conditions on the FiO rate, an optical diagnostic capable of measuring transient FiO on minute timescales was recently developed and demonstrated on a modified GM Ecotech engine system (Neupane et al., Applied Spectroscopy 2021). The measurement is based on adding a dye to the fuel and monitoring for its presence in oil via laser-induced fluorescence (LIF). The measured LIF signal is related to FiO concentration via pre-determined calibration factors using a multivariate classical least square (CLS) method.Since fluorescence quantum yield is a function of temperature, measured LIF intensity not only depends on FiO concentration but also oil temperature. To expand the applicability of the FiO diagnostic to transient oil-temperature conditions (e.g., cold start in practical engines), this study develops a method to account for oil-temperature variations. The effect of oil temperature (20°C - 95°C) on the LIF spectra of eight FiO samples ranging from ~0.8-15% was investigated. LIF intensities of the FiO samples decreased linearly with increasing temperature; the reductions being more significant at dye peaks. We develop a new calibration model (T-adaptive CLS) incorporating the temperature (T) effects on LIF intensity that enables simultaneous calculation of FiO and oil temperature. The improved FiO diagnostic with T-adaptive calibration is more robust, and applicable to varying oil-temperature conditions. For example, when strategies such as multiple/late fuel injections are applied to overcome cold-start instability due to use of low vaporization bio-based fuels such as ethanol, FiO rate is expected to be very high; the improved T-adaptive FiO diagnostic is hence relevant for engine- and fuel-system calibration and optimization. The diagnostic could also provide validation data for flow-field and spray interaction CFD models, further broadening the diagnostic’s utility for advancing engine technology and efficiency.

Neupane, Sneha↗

Majorana Edge Modes in Isolated Wires

Topological superconductors are believed to host exotic quasiparticle excitations known as Majorana zero modes (MZMs), with much of the evidence based on BCS mean-field theory. The direct application of mean-field arguments is tenuous in finite, isolated systems relevant in some experiments. Here, we develop a new correlation-based method for identifying MZMs in interacting, number-conserving systems. Using the density matrix renormalization group, we study fermion number-conserving models with long-range interactions, which under periodic boundary conditions exhibit robust topological and nontopological superconductivity, tuned by the strength of interaction [Ortiz et al., Phys. Rev. Lett. 113, 267002 (2014)]. We find evidence that, on the topological side, Majorana edge modes appear in open chains, manifesting as the vanishing of the energy splitting between odd- and even-parity ground states with increasing system size. Additionally, off-diagonal two-point correlation functions show nonlocal, parity-dependent edge effects. These correlations reveal the spatial structure of Majorana modes in the many-body wave function. We show that the correlation diagnostic applies broadly, including to short-range interacting models, where topological superconductivity is more fragile due to the absence of a bulk excitation gap.

Thomas-Markarian, Jaden↗

Crowdsourcing the Frontier: Advancing Hybrid Physics‐ML Climate Simulation via a $\$$50,000 Kaggle Competition

Subgrid machine-learning (machine learning [ML]) parameterizations have the potential to introduce a new generation of climate models that incorporate the effects of higher-resolution physics without incurring the prohibitive computational cost associated with more explicit physics-based simulations. However, important issues, ranging from online instability to inconsistent online performance, have limited their operational use for long-term climate projections. To more rapidly drive progress in solving these issues, domain scientists and ML researchers opened up the offline aspect of this problem to the broader ML and data science community with the release of ClimSim, a NeurIPS Data sets and Benchmarks publication, and an associated Kaggle competition. This paper reports on the downstream results of the Kaggle competition by coupling emulators inspired by the winning teams' architectures to an interactive climate model (including full cloud microphysics, a regime historically prone to online instability) and systematically evaluating their online performance. Our results demonstrate that online stability in the low-resolution real-geography setting is reproducible across multiple diverse architectures, which we consider a key milestone. All tested architectures exhibit strikingly similar offline and online biases, though their responses to architecture-agnostic design choices (e.g., expanding the list of input variables) can differ significantly. Multiple Kaggle-inspired architectures achieve state-of-the-art results on certain metrics such as zonal mean bias patterns and global Root Mean Squared Error, indicating that crowdsourcing the essence of the offline problem is one path to improving online performance in hybrid physics-AI climate simulation.

Environmental sciences↗

Statistical analysis for the neutrinoless double-β-decay matrix element of 48 Ca

Neutrinoless double-β-decay (0⁢vββ) nuclear matrix elements (NME) are the object of many theoretical calculation methods, and are very important for analysis and guidance of a large number of experimental efforts. However, there are large discrepancies between the NME values provided by different methods. Here, in this paper, we propose a statistical analysis of the 48 Ca 0v⁢ββ NME using the interacting shell model, emphasizing the range of the NME probable values and their correlations with observables that can be obtained from the existing nuclear data. Based on this statistical analysis with three independent effective Hamiltonians, we propose a common probability distribution function for the 0⁢vββ NME, which has a range of (0.45–0.95) at 90% confidence level, and a mean value of 0.68.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Acceleration mechanisms of energetic ion debris in laser-driven tin plasma EUV sources

Laser-driven tin plasmas are driving new-generation nanolithography as sources of extreme ultraviolet (EUV) radiation centered at 13.5 nm. A major challenge facing industrial EUV source development is predicting energetic ion debris produced during the plasma expansion that may damage the sensitive EUV channeling multilayer optics. Gaining a detailed understanding of the plasma dynamics and ion acceleration mechanisms in these sources could provide critical insights for designing debris mitigation strategies in future high-power EUV sources. Here, we develop a fully kinetic model of tin-EUV sources using one-dimensional particle-in-cell simulations to study ion debris acceleration, which will be valuable for cross-validation of radiation-hydrodynamic simulations. An inverse-bremsstrahlung heating operator is used to model the interaction of a tin target with an Nd:YAG laser, and thermal conduction is included through a Monte Carlo Coulomb collision operator. While the large-scale evolution is in reasonable agreement with analogous hydrodynamic simulations, the significant timescale for collisional equilibration between electrons and ions allows for the development of prominent two-temperature features. A collimated flow of energetic ions is produced with a spectrum that is significantly enhanced at high energies compared to fluid simulations. The dominant acceleration mechanism is found to be a large-scale electric field supported mainly by the electron pressure gradient, which is enhanced in the kinetic simulations due to the increased electron temperature. We discuss the implications of these results for future modeling of tin-EUV sources and the development of debris mitigation schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development and transferability of neural-network models for plasma-surface interactions

Plasma-surface interactions are increasingly critical to modern technologies; yet, accurate molecular dynamics simulations remain limited by the capabilities of interatomic potentials. Deep Potentials (DPs) promise to revolutionize the field by providing a systematic method for producing accurate interatomic potentials. The primary challenge of DP development is selecting a dataset, which efficiently spans the set of atomic environments one expects to encounter in the subsequent molecular dynamics simulations. The computational cost of density functional theory calculations, which are the typical basis for DP development, makes it impossible to directly verify the quality of a given DP. To address this challenge, we explore the development of a deep-learned interatomic potential, “DeepREBO,” trained to reproduce the behavior of the REBO2 empirical potential, enabling direct validation of training methodology and transferability. Using an active learning framework, we begin with a minimal dataset and iteratively expand it to train a Deep Potential-Smooth Edition model that faithfully reproduces REBO2 results for 25 eV hydrogen bombardment of diamond (001), a particularly challenging case. We show that small, carefully curated datasets can outperform large, unguided ones, with effective models requiring fewer than 15 000 snapshots. Subsequent transferability tests demonstrate that while DeepREBO generalizes well to diamond (111) surfaces, performance degrades for amorphous carbon or higher-energy impacts, highlighting the need for use-case-specific training data. We also evaluate methods to improve short-range repulsion. This study outlines best practices for training robust deep potentials and underscores the importance of dataset design for predictive plasma simulations.

Ab-initio molecular dynamics↗

Chaotic Dynamics Driven by Particle-Core Interactions

High-intensity beams in modern linacs are frequently encircled by diffuse halos, which drive sustained particle losses and result in gradual degradation of accelerating structures. In large part, the growth of halos is facilitated by internal space-charge forces within the beams, and detailed characterization of this process constitutes an active area of ongoing research. A partial understanding of dynamics that ensue within space-charge dominated beams is presented by the particle-core interaction paradigm – a mathematical model wherein single particle dynamics, subject to the collective potential of the core, are treated as a proxy for the broader behavior of the beam. In this work, we investigate the conditions for the onset of large-scale chaos within the framework of this model, and demonstrate that the propensity towards stochastic evolution is strongly dependent upon the charge distribution of the beam. In particular, we show that while particle motion within a uniformly charged beam is dominantly regular, rapid deterministic chaos readily arises within space-charge dominated Gaussian beams. Importantly, we find that for sufficiently high values of the beam’s space charge and beam pulsation amplitude, enhanced chaotic mixing between the core and the halo can lead to an enhanced radial diffusion of charged particles. We explain our results from analytic grounds by demonstrating that chaotic motion is driven by the intersection of two principal resonances of the system, and derive the relevant overlap conditions. Additionally, our analysis illuminates a close connection between the mathematical formulation of the particle-core interaction model and the Andoyer family of integrable Hamiltonians

43 PARTICLE ACCELERATORS↗

Superfluid effective field theory for dark matter direct detection

We develop an effective field theory (EFT) framework for superfluid 4 He to model the interactions among quasiparticles, helium atoms and probe particles. Our effective field theory approach brings together symmetry arguments and power-counting and matches to classical fluid dynamics. We then present the decay and scattering rates for the relevant processes involving quasiparticles and helium atoms. The presented EFT framework and results can be used to understand the dynamics of thermalization in the superfluid, and can be further applied to sub-GeV dark matter direct detection with superfluid 4 He.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optimization of the BDX experiment for Light Dark Matter searches at JLab

The Light Dark Matter (LDM) hypothesis postulates the existence of a new class of sub GeV particles, neutral under Standard Model (SM) interactions. In its simplest form, LDM consists of particles ¿ with masses below 1 GeV/c2, interacting with SM particles via a new force mediated by a light, spin-1 boson A0, commonly referred to as “Dark Photon”. This framework envisions a distinct “Dark Sector” with its own particles and interactions, offering a theoretically well motivated explanation for Dark Matter, consistent with astrophysical observations and a thermal production mechanism. The Beam Dump eXperiment (BDX) is an approved experiment at Jefferson Lab designed to search for Light Dark Matter. BDX will utilize an 11 GeV electron beam impinging on a thick target to produce a forward-boosted secondary beam of Light Dark Matter particles, which will then be detected by a dedicated downstream detector. Approved in 2018, BDX is expected to be commissioned in 2026 and run in 2027-2029. My thesis focuses on the preparatory work to deploy the BDX experiment. The detector design has been optimized to balance practicality with enhanced Light Dark Matter detection capabilities. Extensive characterization of detector components has been performed to ensure a precise understanding of detector response. A custom Monte Carlo framework has been developed to simulate Light Dark Matter signal and explore various theoretical models of interest. Additionally, a comprehensive data analysis framework has been developed to maximize the experiment sensitivity to Light Dark Matter, with optimizations based on expected detector performance. The ultimate goal of this thesis is to optimize BDX as a flagship experiment in Light Dark Matter searches, enabling it to probe different Dark Matter models.

Spreafico, M. [Univ. of Genova (Italy)]↗

Development of a kinetic model to describe six types of symbiotic interactions in a formate utilizing microalgae-bacteria cultivation system

This study investigated an algae/formate-utilizing-bacteria system that has been developed for carbon capture. This photomixotrophic consortium consumed formate to support bacterial growth so that the resulting respiration CO 2 could be simultaneously used by algae to avoid CO 2 gas-lipid mass transfer limitation. To understand biomass growth and population interactions in this unique system, a kinetic model has been developed to describe algae and bacteria multiplication, formate and nitrogen utilizations, CO 2 mass transfer, O 2 generation and consumption, lighting condition and shading effect. The simulation indicated that this ecosystem could form six types of interactions (mutualism, commensalism, parasitism, neutralism, amensalism, and competition) depending on the light intensity and nutrient availability. Furthermore, the simulation of both batch culture and chemostat, along with parameter sensitivity tests and experimental observations, offered insights into optimal applications of algae-bacteria consortium based photobiorefinery.

59 BASIC BIOLOGICAL SCIENCES↗

Indirect climate impacts of the Hunga eruption

Injection of sulfur and water vapour by the Hunga volcanic eruption significantly altered chemical composition and radiative budget of the stratosphere. Yet, whether the eruption could also affect surface climate, especially via indirect pathways, remains poorly understood. Here we investigate these effects using large ensembles of simulations with the CESM2(WACCM6) Earth system model, incorporating interactive chemistry and aerosols in both coupled ocean and atmosphere-only configurations. We find some statistically significant extratropical regional climate responses to the eruption driven by circulation changes; these are partially linked to the modulation of El Niño–Southern Oscillation, and its associated teleconnections, and to perturbations of the stratospheric polar vortex in both hemispheres. The stratospheric anomalies affect surface climate through modulating the North Atlantic Oscillation in the Northern Hemisphere (up to three boreal winters following the eruption) and the Southern Annular Mode in the Southern Hemisphere in late 2023. The latter is partly related to a concurrent reduction in Antarctic ozone, as increased stratospheric aerosols and water vapor reach the polar vortex. Our study suggests that the eruption could have had a non-negligible influence on regional surface climate, and discusses the mechanisms via which such an influence could occur. However, the results also highlight that this forcing is relatively weak compared to interannual variability, and is subject to model uncertainties in the representation of key processes. More research is thus needed before definitive statements on the role of the eruption in contributing to surface climate and weather events in the following years are made.

Geosciences↗

Identifying mitigation strategies for $\mathrm{COVID}$-$\mathrm{19}$ superspreading on flights using models that account for passenger movement

Despite commercial airlines mandating masks, there have been multiple documented events of COVID-19 superspreading on flights. Conventional models do not adequately explain superspreading patterns on flights, with infection spread wider than expected from proximity based on passenger seating. An important reason for this is that models typically do not consider the movement of passengers during the flight, boarding, or deplaning. Understanding the risks for each of these aspects could provide insight into effective mitigation measures. We modeled infection risk from seating and fine-grained movement patterns – boarding, deplaning, and inflight movement. We estimated infection model parameters from a prior superspreading event. We validated the model and the impact of interventions using available data from three flights, including cabin layout and seat locations of infected and uninfected passengers, to suggest interventions to mitigate COVID-19 superspreading events during air travel. Specifically, we studied: 1) London to Hanoi with 201 passengers, including 13 secondary infections among passengers; 2) Singapore to Hangzhou with 321 passengers, including 12 to 14 secondary infections; 3) a non-superspreading event on a private jet in Japan with 9 passengers and no secondary infections. Here our results show that the inclusion of passenger movement better explains the infection spread patterns than conventional models do. We also found that FFP2/N95 mask usage would have reduced infection by 95–100%, while cloth masks would have reduced it by only 40–80%. Results indicate that leaving the middle seat vacant is effective in reducing infection, and the effectiveness increases when combined with good quality masks. However, with a good mask, the risk is quite low even without the middle seats being empty. Our results suggest the need for more stringent guidelines to reduce aviation-related superspreading events of COVID-19.

54 ENVIRONMENTAL SCIENCES↗