Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetically constrained models”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Validation of an integrated modeling framework for investigating 3D plasma responses in tokamak plasmas

As contemporary experimental tokamaks are pushed toward reactor-relevant operation, they provide essential testbeds for demonstrating and exploring ELM control strategies for deployment in future fusion pilot plants. Accurate predictions of full plasma responses are essential to guide and optimize these demonstrations. This paper introduces and validates an integrated modeling framework over a historical range of DIII-D operational space. The integrated modeling framework uses only scalar plasma parameters and optional reference boundary to self-consistently and flexibly scan through tokamak operational space and estimate the corresponding ELM-suppression relevant plasma response over 3D coil phase space. The framework generates tightly converged equilibria that satisfy a target set of plasma parameters (I p , β N , l i ), with kinetic profiles constrained by an EPED(NN)-computed pedestal and empirical core model. The plasma response of these modeled equilibria is calculated with GPEC. Validated against 55 distinct DIII-D equilibria, the framework consistently reproduces experimentally constrained equilibrium pressure, q, and other representative profiles using the time-varying information from only evolving scalar plasma parameters; plasma response validation is performed for the n=3 perturbation. In particular, the plasma response of a DIII-D discharge scanning q 95 in search of ELM suppression windows is quantitatively reproduced. By accurately reproducing experimental equilibria and plasma responses across wide parameter variation, the framework supports potential for synthetic parameter scans in key operational, stability, and plasma response dimensions to investigate RMP ELM-suppression experiments and inform predictive RMP scenario optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The role of the droplet interface in controlling the multiphase oxidation of thiosulfate by ozone

Predicting reaction kinetics in aqueous microdroplets, including aerosols and cloud droplets, is challenging due to the probability that the underlying reaction mechanism can occur both at the surface and in the interior of the droplet. Additionally, few studies directly measure the surface activities of doubly charged anions, despite their prevalence in the atmosphere. Here, deep-UV second harmonic generation spectroscopy is used to probe surface affinities of the doubly charged anions thiosulfate, sulfate, and sulfite, key species in the thiosulfate ozonation reaction mechanism. Thiosulfate has an appreciable surface affinity with a measured Gibbs free energy of adsorption of -7.3 ± 2.5 kJ mol -1 in neutral solution, while sulfate and sulfite exhibit negligible surface propensity. The Gibbs free energy is combined with data from liquid flat jet ambient pressure X-ray photoelectron spectroscopy to constrain the concentration of thiosulfate at the surface in our model. Stochastic kinetic simulations leveraging these novel measurements show that the primary reaction between thiosulfate and ozone occurs at the interface and in the bulk, with the contribution of the interface decreasing from ~65% at pH 5 to ~45% at pH 13. Additionally, sulfate, the major product of thiosulfate ozonation and an important species in atmospheric processes, can be produced by two different pathways at pH 5, one with a contribution from the interface of >70% and the other occurring predominantly in the bulk (>98%). The observations in this work have implications for mining wastewater remediation, atmospheric chemistry, and understanding other complex reaction mechanisms in multiphase environments. Future interfacial or microdroplet/aerosol chemistry studies should carefully consider the role of both surface and bulk chemistry.

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

Number of sites-based solver for determining coverages from steady-state mean-field micro-kinetic models

Kinetic models parameterized by ab-initio calculations have led to significant improvements in understanding chemical reactions in heterogeneous catalysis. These studies have been facilitated by implementations which determine steady-state coverages and rates of mean-field micro-kinetic models. As implemented in the open-source kinetic modeling program, CatMAP, the conventional solution strategy is to use a root-finding algorithm to determine the coverage of all intermediates through the steady-state expressions, constraining all coverages to be non-negative and to properly sum to unity. Though intuitive, this root-finding strategy causes issues with convergence to solution due to these imposed constraints. In this work, we avoid explicitly imposing these constraints, solving the mean-field steady-state micro-kinetic model in the space of number of sites instead of solving it in the space of coverages. Here, we transform the constrained root-finding problem to an unconstrained least-squares minimization problem, leading to significantly improved convergence in solving micro-kinetic models and thus enabling the efficient study of more complex catalytic reactions.

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

Design and application of a kinetic model of lipid metabolism in Saccharomyces cerevisiae

Lipid biosynthesis plays a vital role in living cells and has been increasingly engineered to overproduce various lipid-based chemicals. However, owing to the tightly constrained and interconnected nature of lipid biosynthesis, both understanding and engineering of lipid metabolism remain challenging, even with the help of mathematical models. Here we report the development of a kinetic metabolic model of lipid metabolism in Saccharomyces cerevisiae that integrates fatty acid biosynthesis, glycerophospholipid metabolism, sphingolipid metabolism, storage lipids, lumped sterol synthesis, and the synthesis and transport of relevant target-chemicals, such as fatty acids and fatty alcohols. The model was trained on lipidomic data of a reference S. cerevisiae strain, single knockout mutants, and lipid overproduction strains reported in literature. The model was used to design mutants for fatty alcohol overproduction and the lipidomic analysis of the resultant mutant strains coupled with model-guided hypothesis led to discovery of a futile cycle in the triacylglycerol biosynthesis pathway. In addition, the model was used to explain successful and unsuccessful mutant designs in metabolic engineering literature. Thus, this kinetic model of lipid metabolism can not only enable the discovery of new phenomenon in lipid metabolism but also the engineering of mutant strains for overproduction of lipids.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

A kinetic model for simulating non-equilibrium mass transport in oxides applied to hematite growth under irradiation

We propose an approach to simulating the dynamic evolution and transport of charged point defects within and through the oxide scales that form during corrosion. The method follows the cluster dynamics formalism widely adopted for radiation damage in solids, which can apply in both quasi-static and far from equilibrium conditions, such as in radiation environments. By treating each charged defect as a cluster of an atomic defect and associated unit charges, the proposed model flexibly allows charge state transitions and shifts in the Fermi level by absorption and emission of charge carriers from point defects in a reaction network with rates constrained by local equilibrium. Applying this model to hematite predicts changes in self-diffusion and oxidation kinetics in irradiation environments, surprisingly showing reduced oxidation rates in many conditions, despite enhanced self-diffusion. In conclusion, the origin of this effect lies in a change in Fermi level induced by excess vacancies formed under irradiation, which in turn suppresses the transport of cation interstitials which facilitate hematite growth.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Toward the performance assessment of advanced nuclear waste forms: temperature dependence of lanthanide borosilicate glass dissolution

Lanthanide borosilicate (LaBS) glasses are among the most promising waste forms for the immobilization of high-level radioactive waste generated from advanced nuclear fuel cycles. However, the temperature dependence of their dissolution kinetics remains poorly understood and constrained, limiting the integration of these materials into established performance assessment models. Here, we investigate the dissolution behavior of the legacy AmCm2-19 LaBS glass and the benchmark alkali aluminoborosilicate ISG-1 in deionized water between 50 °C and 250 °C using ASTM C1285 (Product Consistency Test-B) protocols. For AmCm2-19 LaBS glass, normalized elemental release rates for boron and silicon increase with temperature before plateauing near 150 °C, consistent with solubility-limited behavior. From data obtained at 50 °C and 100 °C, Arrhenius analysis yields activation energies of E a (B) = 24.8 ± 0.3 kJ mol⁻¹ and E a (Si) = 14.4 ± 0.2 kJ mol⁻¹, similar or slightly lower than those previously reported for two other compositions of LaBS glasses. No secondary phases or alteration layers were detected by SEM-EDX or pXRD. These results establish one of the first temperature-dependent kinetic datasets for LaBS glass dissolution, providing quantitative parameters to inform mechanistic corrosion models and predictive simulations of glass degradation in geological disposal environments.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Multiscale modeling of solute diffusion in triblock copolymer membranes

We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model’s predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Furthermore, our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Topological kinetic crossover in a nanomagnet array

Ergodic kinetics, which are critical to equilibrium thermodynamics, can be constrained by a system’s topology. Here, we studied a model nanomagnetic array in which such constraints visibly affect the behavior of the magnetic moments. In this system, magnetic excitations connect into thermally active one-dimensional strings whose motion can be imaged in real time. At high temperatures, our data showed the merging, breaking, and reconnecting of strings, resulting in the system transitioning between topologically distinct configurations. Below a crossover temperature, the string motion is dominated by simple changes in length and shape. In this low-temperature regime, the system is energetically stable because of its inability to explore all possible topological configurations. This kinetic crossover suggests a generalizable conception of topologically broken ergodicity and limited equilibration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Observationally Constrained Modeling of the Reactive Uptake of Isoprene-Derived Epoxydiols under Elevated Relative Humidity and Varying Acidity of Seed Aerosol Conditions

Isoprene is the non-methane volatile organic compound (VOC) emitted in largest amounts to the atmosphere, and it is a significant source of secondary organic aerosol (SOA) mass. The uptake of isoprene oxidation products followed by multiphase chemistry in fine particles is the key pathway to form isoprene epoxydiol-derived SOA (IEPOX-SOA). However, many parameters that relate to diffusion and reaction of IEPOX in the particle phase remain uncertain, since reaction kinetics previously measured in bulk aqueous phase solutions might be different from atmospheric aerosols. Here, we use simultaneous environmental chamber measurements of multiple parameters governing IEPOX-SOA formation at timescales of ~hours: particle size distribution, composition, and volatility of IEPOX-SOA to constrain the key parameters governing IE-POX-SOA formation under humid (i.e., 50% relative humidity, RH) and varying seed aerosol acidity conditions. Reducing the 2-methyltetrol (tetrol) reaction rate constants by a factor of 4 brings the model predictions in agreement with the IEPOX-SOA measurements with acidified ammonium bisulfate seed aerosols. For less acidic ammonium sulfate aerosols both the organo-sulfate (OS) and tetrol reaction rate constants need to be reduced to bring model predictions closer to chamber observations. Using the measured non-volatile content of IEPOX-SOA we constrain the oligomerization timescale of 2-methyltetrols. We find that the oligomerization timescale is 4 hours with acidified seed aerosols, but a much longer time scale of 24 hours is needed for non-acidified seed aerosols, indicating that the aerosol acidity greatly affects the oligomerization rate of tetrols. Here we show that the actual kinetics of IEPOX-SOA formation rate on aerosol seeds consisting of both ammonium bisulfate and ammonium sulfate are a factor of 4~5 slower under 50-60% RH conditions compared to their application in previous models, which were based on bulk aqueous solution measurements.

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

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗
DOE OSTI2025

Understanding the Competition between Alcohol Formation and Dimerization during Electrochemical Reduction of Aromatic Carbonyl Compounds

The electrochemical reductive dimerization of small aromatic carbonyl compounds derived from lignocellulosic biomass is a crucial C−C coupling reaction for upgrading small molecules to longchain hydrocarbons, particularly in the synthesis of drop-in sustainable aviation fuels. Although other electrochemical reduction reactions of these reactants (i.e., hydrogenation and hydrogenolysis) have undergone extensive mechanistic investigation, the understanding of dimerization remains relatively underdeveloped. Most importantly, there is a lack of understanding of the selectivity-determining step between dimerization and monomer reduction and critical factors that can affect this step. In this study, we provide a comprehensive mechanistic model to explain the competition between dimerization and monomer reduction of benzaldehyde under various conditions. Our model proposes that the selectivity between dimerization and monomer reduction depends on the competition between desorption of a ketyl radical from the electrode, necessary for dimerization, and further reduction of the ketyl radical to an alcohol on the electrode by proton-coupled electron transfer (PCET). Computationally comparing the adsorption/desorption energy and PCET activation barrier energy is challenging because conventional DFT calculations substantially underestimate the PCET kinetic barriers. In this study, we employed constrained DFTbased configuration interaction (CDFT-CI) to obtain a reliable comparison of these energies. Our mechanistic model was tested and supported by experimental results obtained with four electrodes (Cu, Pb, Bi, graphite), three pH conditions (acidic, neutral, basic), and three potentials. Our study offers a coherent mechanistic foundation that can explain how each of these conditions impacts the desorption and PCET processes and the selectivities for dimerization and alcohol production.

09 BIOMASS FUELS↗
DOE OSTI2024

Convergence in simulating global soil organic carbon by structurally different models after data assimilation

Abstract Current biogeochemical models produce carbon–climate feedback projections with large uncertainties, often attributed to their structural differences when simulating soil organic carbon (SOC) dynamics worldwide. However, choices of model parameter values that quantify the strength and represent properties of different soil carbon cycle processes could also contribute to model simulation uncertainties. Here, we demonstrate the critical role of using common observational data in reducing model uncertainty in estimates of global SOC storage. Two structurally different models featuring distinctive carbon pools, decomposition kinetics, and carbon transfer pathways simulate opposite global SOC distributions with their customary parameter values yet converge to similar results after being informed by the same global SOC database using a data assimilation approach. The converged spatial SOC simulations result from similar simulations in key model components such as carbon transfer efficiency, baseline decomposition rate, and environmental effects on carbon fluxes by these two models after data assimilation. Moreover, data assimilation results suggest equally effective simulations of SOC using models following either first‐order or Michaelis–Menten kinetics at the global scale. Nevertheless, a wider range of data with high‐quality control and assurance are needed to further constrain SOC dynamics simulations and reduce unconstrained parameters. New sets of data, such as microbial genomics‐function relationships, may also suggest novel structures to account for in future model development. Overall, our results highlight the importance of observational data in informing model development and constraining model predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Combined density functional theory/kinetic Monte Carlo investigation of surface morphology during cycling of Li-Cu electrodes

In this work, we employed a combined density functional theory and kinetic Monte Carlo approach in order to investigate the nanoscale behavior of lithium stripping and plating on lithium-copper nanoalloys. Density functional theory calculations were used to develop an atom-by-atom understanding of the nucleation behavior of lithium on copper surfaces with different surface morphologies. We found that nanometric pits resulted in significantly more favorable lithium deposition when compared to pristine or lightly-defective Cu surfaces, which is in agreement with previous experimental studies. We then used density functional theory and constrained ab initio molecular dynamics to develop interaction and reaction parameters for a larger-scale kinetic Monte Carlo model in order to provide insight into the stripping and plating behavior of nanostructured Li-Cu nanoalloy anodes. Substantial losses in capacity and strong indications of morphological change were observed for a wide range of Cu:Li ratios over 100 cycles. We observed a complex relationship between the nanoparticle composition and the properties of the formed solid-electrolyte interphase.

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

Theoretical and kinetic modeling study of H 2 S pyrolysis

Hydrogen sulfide pyrolysis was investigated theoretically and through chemical kinetic modeling. Reactions on the SHH potential energy surface, primarily S + H 2 (+Ar) ⇌ H 2 S (+ Ar) (R1) and S + H 2 ⇌ SH + H (R6b) were characterized by ab initio calculations. Results for k 1 were in good agreement with experiment, but deviated strongly below 2000 K from values previously used in modeling. Collider efficiencies for H 2 S, S 2 , and N 2 compared to Ar were calculated for R1. Hydrogen sulfide decomposition experiments reported in literature were re-examined in terms of an updated detailed chemical kinetic model. Concentration profiles for the atomic S at high temperature in shock tubes supported the present value of k 1 and served to constrain the rate constants for reaction of S with SH and H 2 S. To explain results from batch and flow reactors, conducted at high H 2 S concentrations in the 900–1400 K range, a very fast rate constant was required for HSS + H ⇌ SH + SH. Under dilute conditions, the gas-phase chemistry was too slow to compete and the decomposition of H 2 S was controlled by loss on the reactor surface.

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

Toward Higher Voltage Solid-State Batteries by Metastability and Kinetic Stability Design

The energy density of battery systems is limited largely by the electrochemical window of the electrolyte. Here, the combined thermodynamic and kinetic effects of mechanically induced metastability are shown to greatly widen the operational voltage window of solid-state batteries based on ceramic-sulfide electrolytes. Solid electrolyte voltage stability up to 10 V is achieved with minimal degradation, far beyond the capability of organic liquid electrolytes. Furthermore, combined experiment, ab initio computation, and theoretical modeling identify the nature of mechanically constrained Li 10 GeP 2 S 12 decomposition both within the bulk and at interfaces with cathode materials at very high voltages. Previously unclear kinetic processes are identified that, when properly implemented, can potentially allow solid-state full cells with remarkably high operational voltages.

25 ENERGY STORAGE↗
DOE OSTI2022

Structural magnetic glassiness in the spin ice Dy 2 Ti 2 O 7

The origin and nature of glassy dynamics presents one of the central enigmas of condensed-matter physics across a broad range of systems ranging from window glass to spin glasses. The spin-ice compound Dy 2 Ti 2 O 7 , which is perhaps best known as hosting a three-dimensional Coulomb spin liquid with magnetically charged monopole excitations, also falls out of equilibrium at low temperature. How and why it does so remains an open question. Based on an analysis of low-temperature diffuse neutron-scattering experiments employing different cooling protocols alongside recent magnetic noise studies, combined with extensive numerical modeling, we argue that upon cooling, the spins freeze into what may be termed a “structural magnetic glass,” without an a priori need for chemical or structural disorder. Specifically, our model indicates the presence of frustration on two levels, first producing a near-degenerate constrained manifold inside which phase ordering kinetics is in turn frustrated. A remarkable feature is that monopoles act as sole annealers of the spin network and their pathways and history encode the development of glass dynamics, allowing the glass formation to be visualized. Our results suggest that spin ice Dy 2 Ti 2 O 7 provides one prototype of magnetic glass formation specifically and a setting for the study of kinetically constrained systems more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

System-size dependence of the hadronic rescattering effect at energies available at the CERN Large Hadron Collider

The first measurements of K * ( 892 ) 0 resonance production as a function of charged-particle multiplicity in Xe-Xe collisions at s N N = 5.44 TeV and p p collisions at s = 5.02 TeV using the ALICE detector are presented. The resonance is reconstructed at midrapidity ( | y | < 0.5) using the hadronic decay channel K * 0 → K ± π ∓ . Measurements of transverse-momentum integrated yield, mean transverse-momentum, nuclear modification factor of K * 0 , and yield ratios of resonance to stable hadron ( K * 0 / K ) are compared across different collision systems ( p p , p -Pb, Xe-Xe, and Pb-Pb) at similar collision energies to investigate how the production of K * 0 resonances depends on the size of the system formed in these collisions. The hadronic rescattering effect is found to be independent of the size of colliding systems and mainly driven by the produced charged-particle multiplicity, which is a proxy of the volume of produced matter at the chemical freeze-out. In addition, the production yields of K * 0 in Xe-Xe collisions are utilized to constrain the dependence of the kinetic freeze-out temperature on the system size using the hadron resonance gas–partial chemical equilibrium model. ©2024 CERN, for the ALICE Collaboration 2024 CERN

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A hybrid-kinetic simulation tool for non-thermal warm x-ray z-pinch sources, with gas-puff and wire array exemplars

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗