Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Debye temperature”

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
NASA NTRS1992

Collision integrals for the interaction of the ions of nitrogen and oxygen in a plasma at high temperatures and pressures

The corrections to the transport cross-sections and collision integrals for Coulomb interactions arising from the application of realistic interaction energies of the ions of nitrogen and oxygen are investigated. Accurate potential-energy curves from an ab initio electronic-structure calculation and a semiclassical description of the scattering are used to determine the difference between the cross-sections for the real interaction forces and a Coulomb force for large values of the Debye shielding parameter. Graphs of the correction to the diffusion and viscosity-collision integrals are presented for temperatures from about 10,000 K to 150,000 K. This correction can be combined with tabulations of the collision integrals for shielded Coulomb potentials to determine the contribution of N(+)-N(+), N(+)-O(+), and O(+)-O(+) interactions to the transport properties of high-temperature air. Analytical forms are fitted to the calculated results to assist this application.

Stallcop, James R.↗
DOE OSTI2021

\( \mathcal{N} \) = 4 supersymmetric Yang-Mills thermodynamics to order λ 2

We calculate the resummed perturbative free energy of \( \mathcal{N} \) = 4 supersymmetric Yang-Mills in four spacetime dimensions (SYM 4 , 4 ) through second order in the ’t Hooft coupling λ at finite temperature and zero chemical potential. Our final result is ultraviolet finite and all infrared divergences generated at three-loop level are canceled by summing over SYM 4 , 4 ring diagrams. Non-analytic terms at \( \mathcal{O} \) ( λ 3 / 2 ) and \( \mathcal{O} \) ( λ 2 log λ ) are generated by dressing the A 0 and scalar propagators. The gauge-field Debye mass m D and the scalar thermal mass M D are determined from their corresponding finite-temperature self-energies. Based on this, we obtain the three-loop thermodynamic functions of SYM 4 , 4 to \( \mathcal{O} \) ( λ 2 ). We compare our final result with prior results obtained in the weak- and strong-coupling limits and construct a generalized Padé approximant that interpolates between the weak-coupling result and the large- N c strong-coupling result. Our results suggest that the \( \mathcal{O} \) ( λ 2 ) weak-coupling result for the scaled entropy density is a quantitatively reliable approximation to the scaled entropy density for 0 ≤ λ ≲ 2.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A gyrokinetic simulation model for 2D equilibrium potential in the scrape-off layer of a field-reversed configuration

The equilibrium potential structure in the scrape-off layer (SOL) of the field-reversed configuration (FRC) can be affected by the penetration of edge biasing applied at the divertor ends. The primary focus of the paper is to establish a formulation that accurately captures both parallel and radial variations of the two-dimensional (2D) potential in SOL. The formulation mainly describes a quasi-neutral plasma with a logical sheath boundary. A full-f gyrokinetic ion model and a massless electron model are implemented in the GTC-X code to solve for the self-consistent equilibrium potential, given fixed radial potential profiles at the boundaries. The first essential point of this 2D model lies in its ability to couple radial and parallel dynamics stemming from resistive currents and drag force on ions. The model successfully recovers the fluid force balance and continuity equations. These collisional effects on 2D potential mainly appear through the density profile changes, modifying the potential through electron pressure gradient. This means an accurate prescription of electron density and temperature profiles is important in predicting the potential structure in the FRC SOL. The Debye sheath potential and the potential profiles applied at the boundaries can be additional factors contributing to the 2D variations in SOL. This comprehensive full-f scheme holds promise for future investigations into turbulent transport in the presence of the self-consistent 2D potential together with the non-Maxwellian distributions and open boundary conditions in the FRC SOL.

Physics↗
NASA NTRS2021

Seawater Debye Model Function at L-Band and Its Impact on Salinity Retrieval From Aquarius Satellite Data

A model function of seawater, which specifies the dielectric constant of seawater as a function of salinity, temperature, and frequency, is important for the retrieval of sea surface salinity using satellite data. In 2017, a model function has been developed based on measurement data at 1.4134 GHz using a third-order polynomial expression in salinity ( S ) and temperature ( T ). Although the model showed improvements in salinity retrieval, it had an inconsistent behavior between partitioned salinities. To improve the stability of the model, new dielectric measurements of seawater have been made recently over a broad range of salinities and temperatures to expand the data set used for developing the model function. The structure of the model function has been changed from a polynomial expansion in S and T to a physics-based model consisting of a Debye molecular resonance term plus a conductivity term. Each unknown parameter is expressed in S and T based on the expanded measurement data set. Physical arguments have been used to limit the number of unknown coefficients in these expressions to improve the stability of the model function. The new model function has been employed in the retrieval algorithm of the Aquarius satellite mission to obtain a global salinity map. The retrieved salinity using a different model function is compared with in situ data collected by Argo floats to evaluate the impact and the performance of model functions. The results indicate that the new model function has significant improvements in salinity retrieval compared with other existing models.

Seawater↗
DOE OSTI2024

Phase Transitions and Thermal Equation of State of Fe‐9wt.%Si Applied to the Moon and Mercury

Abstract Accurate knowledge of the phase transitions and thermoelastic properties of candidate iron alloys, such as Fe‐Si alloys, is essential for understanding the nature and dynamics of planetary cores. The phase diagrams of some Fe‐Si alloys between 1 atm and 16 GPa have been back‐extrapolated from higher pressures, but the resulting phase diagram of Fe 83.6 Si 16.4 (9 wt.% Si) is inconsistent with temperature‐induced changes in its electrical resistivity between 6 and 8 GPa. This study reports in situ synchrotron X‐ray diffraction (XRD) measurements on pre‐melted and powder Fe 83.6 Si 16.4 samples from ambient conditions to 60 GPa and 900 K using an externally heated diamond‐anvil cell. Upon compression at 300 K, the bcc phase persisted up to ∼38 GPa. The hcp phase appeared near 8 GPa in the pre‐melted sample, and near 17 GPa in the powder sample. The appearance of the hcp phase in the pre‐melted sample reconciles the reported changes in electrical resistivity of a similar sample, thus resolving the low‐pressure region of the phase diagram. The resulting high‐temperature Birch‐Murnaghan equation of state (EoS) and thermal EoS based on the Mie‐Gruneisen‐Debye model of the bcc and hcp structures are consistent with, and complement the literature data at higher pressures. The calculated densities based on the thermal EoS of Fe‐9wt.%Si indicate that both bcc and hcp phases agree with the reported core density estimates for the Moon and Mercury.

Berrada, Meryem↗
NASA NTRS1977

Direct measurement of the plasma screening length and surface potential near the lunar terminator

Direct measurement of the lunar dayside surface potential and screening length has been made by the suprathermal ion detector experiment (Side) near the terminator. In a region 20-30 deg from the terminator at the Apollo 14 and 15 sites the surface potential is found to be approximately 50 V negative, and the screening length to be about 1 km. This value of the screening length is more than 2 orders of magnitude greater than the solar wind 'Debye' length. The strong negative surface potential in this region may be due to enhanced temperature and density of the solar wind plasma.

Benson, J.↗
DOE OSTI2021

Large gyro-orbit model of ion velocity distribution in plasma near a wall in a grazing-angle magnetic field

A model is presented for the ion distribution function in a plasma at a solid target with a magnetic field $\boldsymbol {B}$ inclined at a small angle, $\alpha \ll 1$ (in radians), to the target. Adiabatic electrons are assumed, requiring $\alpha \gg \sqrt {Zm_{e}/m_{i}}$ , where $m_{e}$ and $m_{i}$ are the electron and ion mass, respectively, and $Z$ is the charge state of the ion. An electric field $\boldsymbol {E}$ is present to repel electrons, and so the characteristic size of the electrostatic potential $\phi$ is set by the electron temperature $T_{e}$ , $e\phi \sim T_{e}$ , where $e$ is the proton charge. An asymptotic scale separation between the Debye length $\lambda _{D} = \sqrt {\epsilon _0 T_{{e}} / e^{2} n_{{e}} }$ , the ion sound gyro-radius $\rho _{s} = \sqrt { m_{i} ( ZT_{e} + T_{i} ) } / (ZeB)$ and the size of the collisional region $d_{c} = \alpha \lambda _{\textrm {mfp}}$ is assumed, $\lambda _{D} \ll \rho _{s} \ll d_{c}$ . Here $\epsilon _0$ is the permittivity of free space, $n_{e}$ is the electron density, $T_{i}$ is the ion temperature, $B= |\boldsymbol {B}|$ and $\lambda _{\textrm {mfp}}$ is the collisional mean free path of an ion. The form of the ion distribution function is assumed at distances $x$ from the wall such that $\rho _{s} \ll x \ll d_{c}$ , that is, collisions are not treated. A self-consistent solution of the electrostatic potential for $x \sim \rho _{s}$ is required to solve for the quasi-periodic ion trajectories and for the ion distribution function at the target. The large gyro-orbit model presented here allows to bypass the numerical solution of $\phi (x)$ and results in an analytical expression for the ion distribution function at the target. It assumes that $\tau =T_{i}/(ZT_{e})\gg 1$ , and ignores the electric force on the quasi-periodic ion trajectory until close to the target. For $\tau \gtrsim 1$ , the model provides an extremely fast approximation to energy–angle distributions of ions at the target. These can be used to make sputtering predictions.

Physics↗
NASA NTRS1974

Low-temperature thermostatics of face-centered-cubic metallic hydrogen

The thermostatic properties of a high-symmetry phase of metallic hydrogen with atomic sphere radius between 0.1 and 1.5 bohr are studied, with special emphasis accorded to electronic screening and quantum proton motion. The electron-proton and proton-proton interactions receive a perturbation treatment based on the Singwi dielectric function, while the proton motion is handled by self-consistent harmonic approximation. Quantum behavior is found to be less pronounced than expected, and nuclear magnetism is absent. The phonon spectrum is, however, affected by screening and large proton motion. The zero-point vibrational energy and the superconducting critical temperature are below previous estimates. The crystalline-defect formation energies are a few times the Debye energy, which implies that defects contribute significantly to melting at the lower particle densities.

Caron, L. G.↗
DOE OSTI2023

A prediction of the thermodynamic, thermophysical, and mechanical properties of CrTaO 4 from first principles

It is reported that the self-forming CrTaO 4 oxide scale can protect refractory high-entropy alloys from oxidation, superior to Cr 2 O 3 . In this paper, the phase stability, mechanical, and thermal properties of three polymorphous phases of CrTaO 4 are systematically investigated from first-principles density functional theory calculations. The mechanical properties predicted using the strain–energy methods indicated that all three phases are mechanically stable. The temperature dependence of elastic constants and polycrystalline moduli of three phases demonstrated the thermal softening as temperature increase. The Helmholtz-free energies as a function of volume and temperature are derived from phonon dispersions within the quasi-harmonic approximation at six strained volumes. The calculated apparent bulk coefficients of thermal expansion of these three phases are evaluated, the highest value approximately 13.4× 10 –6 K –1 within a temperature range of 500–2000 K for the rutile I4 1 md phase. The lattice thermal conductivity calculated by the Debye–Callaway model suggested that the rutile type I4 1 md phase has the lowest value of approximately 2.1 W/m/K at 1800 K. The other two phases, C2/m and P2/c, exhibit higher values due to relatively lower Grüneisen parameters and larger phonon velocities. The melting point of CrTaO 4 is predicted to be between 1975 and 2449 K using ab initio molecular dynamics simulations. Finally, this work provides a comprehensive theoretical understanding of the thermodynamic, mechanical, and thermal properties for the new material CrTaO 4 and serves as an example of a viable computational design strategy for improved oxidation resistance of refractory alloys at high temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Short-range ferromagnetic order due to Ir substitutions in single-crystalline Ba(Co 1-x Ir x ) 2 As 2 (0 ≤ x ≤ 0.25)

The ternary-arsenide compound BaCo 2 As 2 was previously proposed to be in proximity to a quantum-critical point where long-range ferromagnetic (FM) order is suppressed by quantum fluctuations. In this paper we report the effect of Ir substitutions for Co on the magnetic and thermal properties of Ba(Co 1-x Ir x ) 2 As 2 (0 ≤ x ≤ 0.25) single crystals. These compositions all crystallize in an uncollapsed body-centered-tetragonal ThCr 2 Si 2 structure with space group I4/mmm. Magnetic susceptibility measurements reveal clear signatures of short-range FM ordering for x ≥ 0.11 below a nearly composition-independent characteristic temperature T cl ≈ 13 K. The small variation of T cl with x, thermomagnetic irreversibility between zero-field-cooled and field-cooled magnetic susceptibility versus T, the occurrence of hysteresis in magnetization versus field isotherms at low field and temperature, and very small spontaneous and remanent magnetizations <0.01 μ B /f.u. together indicate that the FM response arises from short-range FM ordering of FM spin clusters as previously inferred to occur in Ca(Co 1–x Ir x ) 2–y As 2 . Heat-capacity Cp(T) data do not exhibit any clear feature around Tcl, consistent with the very small moments of the FM clusters. The Cp(T) in the paramagnetic temperature regime 25–300 K is well described by the sum of a Sommerfeld electronic contribution and Debye and Einstein lattice contributions where the latter lattice contribution suggests the presence of low-frequency optic modes associated with the heavy Ba atoms in the crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

In situ high-temperature X-ray diffraction study of Sc-doped titanium oxide nanocrystallites

Titanium dioxide is an inexpensive wide-gap highly ionic semiconductor with striking photocatalytic capabilities in several heterogeneous photoredox reactions. A small crystal size is desirable to maximize the surface area, since photocatalytic reactions occur at the surface of a photocatalyst. Presented in this work are the synthesis and microstructural characterization of 4 at.% Sc-doped TiO 2 (4SDT) prepared by water-based co-precipitation. The crystal structure of 4SDT was examined via in situ high-temperature powder X-ray diffraction experiments from 25 to 1200°C. Rietveld analysis revealed single-phase anatase up to 875°C, while at 900°C the anatase-to-rutile phase transformation occurred and at higher temperatures additional reflections of Sc-rich phases (Sc 2 TiO 5 from 975°C and Ti 3 Sc 4 O 12 or Sc 2 O 3 at 1200°C) were observed. Debye function analysis (DFA) was applied to model the total scattering pattern directly in reciprocal space, allowing the reconstruction of Ti vacancies. Both Rietveld and DFA methods were applied to estimate the nanocrystallite size and shape with consistent growth in crystallite size with temperature: an ellipsoid shape with equatorial ~4.7 nm / axial (001) ~6.9 nm at 25°C to equatorial ~27.9 nm / axial (001) ~39.6 nm at 900°C refined by Rietveld analysis, versus a cylinder shape with D a , b = 4.3 nm and size dispersion σ ab = 1.5 nm, L c = 4.9 nm and σ c = 2.3 nm at 25°C to D a , b = 21.4 nm, σ ab = 8.3 nm, L c = 23.9 and σ c = 10.9 nm at 900°C estimated by DFA. The microstructural changes obtained by Rietveld and DFA methods were supported by high-resolution transmission electron microscopy image analysis, as well as by the less direct nitrogen sorption techniques that provide information on the size of non-agglomerated and dense particles. The Ti site-occupancy factor showed a linear increase from 0.6–0.8 at 25°C to unity at 900°C for anatase, and from ~0.7 at 900°C to unity at 1200°C for rutile, via Rietveld analysis and DFA.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Local structure, bonding, and asymmetry of $\mathrm{((NH_2)_2CH)PbBr_3}$, $\mathrm{CsPbBr_3}$, and $\mathrm{(CH_3MH_3)PbBr_3}$

Here we report local structure measurements for CsPbBr 3 and ((NH 2 ) 2 CH)PbBr 3 (FAPbBr 3 ) and compare them with recent results for (CH 3 NH 3 )PbBr 3 (MAPbBr 3 ). The Pb-Br bonding is similar for all three systems; the effective spring constants, κ, are comparable (ranging from 1.20 to 1.95 eV/Å 2 ), but small in magnitude indicating very soft materials. However, there are also important differences between the three systems. Static disorder is very small for CsPbBr3 but increases somewhat with the size of the organic molecular ions MA + and FA + . At room temperature, dynamic disorder dominates in all compounds. The thermal disorder of the Pb-Br pair distribution function (PDF), i.e., the Debye-Waller factor σ 2 follows a correlated Debye or Einstein model up to 300 K in CsPbBr 3 (orthorhombic phase), but for FAPbBr 3 and MAPbBr 3 , there is a break in the σ 2 (T) curve at the orthorhombic-tetragonal transition (o-t) near 150 K, indicating a small change in the spring constant κ. κ increases for MAPbBr 3 but decreases in FAPbBr 3 at this transition. These changes are attributed to changes in the H-bonding between Br - and MA + or FA + at this transition, as a result of librations or rotations of these molecular cations. In addition, the Pb-Br PDF becomes asymmetric at a relatively low temperature for FAPbBr 3 and MAPbBr 3 , while this effect is significantly smaller for CsPbBr 3 . Finally, we address the question of a model to explain the asymmetric PDF. Two main models are discussed in the literature, an anharmonic pair potential and a split-pair distribution, possibly driven by the presence of a lone pair on the Pb ion. We show that the fourth cumulant C 4 can differentiate between these two models and other possible models. Experimentally C 4 is positive at 250 K and above, for all three systems and that is inconsistent with a split-peak model, for which C 4 is negative for splittings larger than 0.12 Å.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Structural and thermal properties of Eu 2 Ga 11 Sn 35

Clathrates have been reported to form in a variety of different structure types; however, inorganic clathrate-I materials with a low-cation concentration have yet to be investigated. Furthermore, tin-based compositions have been much less investigated as compared to silicon or germanium analogs. We report the temperature-dependent structural and thermal properties of single-crystal Eu 2 Ga 11 Sn 35 revealing the effect of structure and composition on the thermal properties of this low-cation clathrate-I material. Specifically, low-temperature heat capacity, thermal conductivity, and synchrotron single-crystal x-ray diffraction reveal a departure from Debye-like behavior, a glass-like phonon mean-free path for this crystalline material, and a relatively large Grüneisen parameter due to the dominance of low-frequency Einstein modes. In conclusion, our analyses indicate thermal properties that are a direct result of the structure and composition of this clathrate-I material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analysis of coherent Thomson scattering from a low temperature plasma

In this work, the spectrum of coherent Thomson scattering (CTS) induced by a periodic ponderomotive perturbation in a low-density low temperature plasma is considered. The analysis is performed for the case when the period of the resulting optical lattice is less than the Debye screening length in the plasma by solving an electron Boltzmann equation, where the total force is the sum of the periodic force due to the optical lattice and the electrostatic force due to self-consistent electric field in the plasma. An analogy between the CTS spectra calculated here and coherent Rayleigh scattering spectra in a neutral gas is established. For relatively low intensity for the optical lattice, the calculated CTS spectra are nearly Gaussian with widths slightly wider than the incoherent Thomson widths. We demonstrate that at higher intensities the line shape narrows and saturates to a width approximately half of that found at low lattice intensities. The proportionality of the spectral width to the square root of the electron temperature allows one to extract the electron temperature from the saturated spectra. Possible application of CTS for remote measuring the electron temperature in plasma is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Thermoelastic properties of zircon: Implications for geothermobarometry

Abstract A thermal-pressure equation of state has been determined for zircon (ZrSiO4) that characterizes its thermoelastic behavior at metamorphic conditions. New pressure-volume (P-V) data from a “Mud Tank” zircon have been collected from 1 bar to 8.47(1) GPa using X-ray diffraction, and elastic moduli were measured from room temperature up to 1172 K by resonance ultrasound spectroscopy. These data were fitted simultaneously with temperature-volume (T-V) data from the literature in EosFit7c using a new scaling technique. The parameters of a third-order Birch-Murnaghan EoS with a Mie-Grüneisen-Debye model for thermal pressure have compressional EoS parameters K0T = 224.5(1.2) GPa, K0T' = 4.90(31) with a fixed initial molar volume V0 = 39.26 cm3/mol and thermal parameters γ0 = 0.868(15), q = 2.37(80), and ΘD = 848(38) K. EoS parameters that describe the variation of unit-cell parameters with pressure and temperature were determined using an isothermal-type EoS. This new EoS confirms that zircons are stiffer than garnets and exhibit a much lower thermal expansion. This results in steep isomekes between zircon and garnets, which makes zircon trapped as inclusions in garnets at metamorphic conditions a good piezothermometer.

Geochemistry & Geophysics↗
DOE OSTI2022

Accurate quantification of lattice temperature dynamics from ultrafast electron diffraction of single-crystal films using dynamical scattering simulations

In ultrafast electron diffraction (UED) experiments, accurate retrieval of time-resolved structural parameters, such as atomic coordinates and thermal displacement parameters, requires an accurate scattering model. Unfortunately, kinematical models are often inaccurate even for relativistic electron probes, especially for dense, oriented single crystals where strong channeling and multiple scattering effects are present. This article introduces and demonstrates dynamical scattering models tailored for quantitative analysis of UED experiments performed on single-crystal films. As a case study, we examine ultrafast laser heating of single-crystal gold films. Comparison of kinematical and dynamical models reveals the strong effects of dynamical scattering within nm-scale films and their dependence on sample topography and probe kinetic energy. Applying to UED experiments on an 11 nm thick film using 750 keV electron probe pulses, the dynamical models provide a tenfold improvement over a comparable kinematical model in matching the measured UED patterns. Also, the retrieved lattice temperature rise is in very good agreement with predictions based on previously measured optical constants of gold, whereas fitting the Debye–Waller factor retrieves values that are more than three times lower. Altogether, these results show the importance of a dynamical scattering theory for quantitative analysis of UED and demonstrate models that can be practically applied to single-crystal materials and heterostructures.

36 MATERIALS SCIENCE↗