Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear reaction”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Estimating viscoelastic compliance of desiccating cementitious materials using drying prism tests

Highlights: • Determining relative humidity and degree of reaction of desiccating cementitious materials via an integrated method • Extracting moisture diffusion coefficient of cement mortar using the mass loss data of the prism specimens • Estimating uniaxial viscoelastic compliance of cement mortar using the shrinkage strain data of the prism specimens An integrated modeling approach to determining internal relative humidity and degree of reaction of desiccating cementitious materials has been developed. By simultaneously accounting for both self-desiccation and external drying, the model is capable of predicting the evolution of material properties of desiccating cementitious materials. Using the mass loss data of the prism specimens, the moisture diffusion coefficient of a nonlinear model for cement mortar was determined. Additionally, a previous poroviscoelastic shrinkage model was modified and extended. Through comparison of predicted and measured axial strain of the prism specimens, the viscoelastic compliance of cement mortar was extracted. The present work revealed that it is possible to use a simple drying prism test to fit the moisture diffusion coefficient and creep compliance of cementitious materials, though the comparison of the extracted viscoelastic compliance with measured compliance was only reasonable if the outer layer of the drying mortar was presumed to have microcracking.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Density Functional Tight-Binding Enables Tractable Studies of Quantum Plasmonics

Routine investigations of plasmonic phenomena at the quantum level present a formidable computational challenge due to the large system sizes and ultrafast time scales involved. This Feature Article highlights the use of density functional tight-binding (DFTB), particularly its real-time time-dependent formulation (RT-TDDFTB), as a tractable approach to study plasmonic nanostructures from a quantum mechanical purview. We begin by outlining the theoretical framework and limitations of DFTB, emphasizing its efficiency in modeling systems with thousands of atoms over picosecond time scales. Applications of RT-TDDFTB are then explored in the context of optical absorption, nonlinear harmonic generation, and plasmon-mediated photocatalysis. We demonstrate how DFTB can reconcile classical and quantum descriptions of plasmonic behavior, capturing key phenomena such as size-dependent plasmon shifts and plasmon coupling in nanoparticle assemblies. Lastly, we showcase DFTB’s ability to model hot carrier generation and reaction dynamics in plasmon-driven H2 dissociation, underscoring its potential to model photocatalytic processes. Collectively, these studies establish DFTB as a powerful, yet computationally efficient tool to probe the emergent physics of materials at the limits of space and time.

Electrical energy↗
NASA NTRS1984

Direct numerical simulation of reacting flows

The objectives of this work are: (1) to extend the technique of direct numerical simulations to turbulent, chemically reacting flows, (2) to test the validity of the method by comparing computational results with laboratory data, and (3) to use the simulations to gain a better understanding of the effects of turbulence on chemical reactions. The effects of both the large scale structure and the smaller scale turbulence on the overall reaction rates are addressed. The relationship between infinite reaction rate and finite reaction rate chemistry is compared with some of the results of calculations with existing theories and laboratory data. The direct numerical simulation method involves the numerical solution of the detailed evolution of the complex turbulent velocity and concentration fields. Using very efficient numerical methods (e.g., pseudospectral methods), the fully nonlinear (possibly low pass filtered) equations of motion are solved and no closure assumptions or turbulence models are used. Statistical data are obtained by performing spatial, temporal, and/or ensemble averages over the computed flow fields.

Riley, J. J.↗
NASA NTRS1986

Do large rate coefficients for ion-polar neutral reactions have a serious effect on chemical models of dense clouds?

In order to incorporate large ion-polar neutral rate coefficients into existing gas phase reaction networks, it is necessary to utilize simplified theoretical treatments because of the significant number of rate coefficients needed. The authors have used two simple theoretical treatments: the locked dipole approach of Moran and Hamill for linear polar neutrals and the trajectory scaling approach of Su and Chesnavich for nonlinear polar neutrals. The former approach is suitable for linear species because in the interstellar medium these are rotationally relaxed to a large extent and the incoming charged reactants can lock their dipoles into the lowest energy configuration. The latter approach is a better approximation for nonlinear neutral species, in which rotational relaxation is normally less severe and the incoming charged reactants are not as effective at locking the dipoles. The treatments are in reasonable agreement with more detailed long range theories and predict an inverse square root dependence on kinetic temperature for the rate coefficient. Compared with the locked dipole method, the trajectory scaling approach results in rate coefficients smaller by a factor of approximately 2.5.

Herbst, E.↗
NASA NTRS1994

Kinetic Study of Methyl Acetate Oxidation in a Pt/Al2O3 Fixed-Bed Reactor

To support technology development for future long-term missions, a metabolic simulator will be used in a closed chamber to test the functions of a Controlled Ecological Life Support System (CELSS). Methyl acetate (MA) was selected as the fuel because its metabolic respiratory quotient is near that of humans. A kinetic study of the catalytic oxidation of MA over Pt/Al203 was then conducted to support the design and operation of the simulator. Kinetic data were obtained as a conversion percentage of MA versus retention time. The reaction was studied at one atmosphere and temperatures from 220 to 340 deg. C. The inlet MA concentration was varied from 100 to 2000 ppm with retention times from 0.01 to 10 sec. A first-order rate law and a Langmuir-Hinshelwood rate equation were tested by nonlinear regression of the kinetic data to estimate rate constants in the rate law. Regression results of the L-H equation explain the kinetic data better than the results of the first-order rate law. A Taguchi experimental design was used to study the effects of temperature, retention time, and concentrations of MA, CO2, and O2 on the conversion of MA. Results indicate that temperature has greatest effect, followed by retention time, and finally MA concentration. It was further determined that the effects of CO2 and O2 concentrations, and the cross effects, are negligible.

Hoy, Michael↗
NASA NTRS2006

Preliminary Study Using Forward Reaction Control System Jets During Space Shuttle Entry

Failure or degradation of the flight control system, or hull damage, can lead to loss of vehicle control during entry. Possible failure scenarios are debris impact and wing damage that could result in a large aerodynamic asymmetry which cannot be trimmed out without additional yaw control. Currently the space shuttle uses aerodynamic control surfaces and Reaction Control System jets to control attitude. The forward jets are used for orbital maneuvering only, while the aft jets are used for yaw control during entry. This paper develops a controller for using the forward reaction control system jets as an additional control during entry, and assesses its value and feasibility during failure situations. Forward-aft jet blending logic is created, and implemented on a simplified model of the space shuttle entry flight control system. The model is validated and verified on the nonlinear, six degree-of-freedom Shuttle Engineering Simulator. A rudimentary human factors study was undertaken using the forward cockpit simulator at Johnson Space Center, to assess flying qualities of the new system and pilot workload. Results presented in the paper show that the combination of forward and aft jets provides useful additional yaw control, in addition to potential fuel savings and the ability to balance the use of the fuel in the forward and aft tanks to meet availability constraints of both forward and aft fuel tanks. Piloted simulation studies indicated that using both sets of jets while flying a damaged space shuttle reduces pilot workload, and makes the vehicle more responsive.

Restrepo, Carolina↗
DOE OSTI2022

Fast-neutron-induced fission of Pu 240 and Pu 242

Herein we report the measurement of the total kinetic energy (TKE) release in the fast neutron induced fission of 240 Pu and 242 Pu. The results are compared to the predictions of the GEF model, the CGMF model, and the model of Denisov and Sedykh as well as previous exptl. work on these reactions. Our absolute measurements of the TKE release are in good agreement with the previous measurements of Nethaway et al. for the interaction of 14.8 MeV neutrons with 240 Pu [Phys. Rev. C16, 1907 (1977)] and of Winkelmann and Aumann for the interaction of 15 MeV neutrons with 242 Pu [Phys. Rev. C30, 934 (1984)]. The general trends of the measured TKE values agree with phenomenol. models but the variances of the TKE distributions are significantly less than predicted by various models. The mean postneutron emission TKE release decreases nonlinearly with increasing neutron energy and can be represented as TKE(MeV) = 175.8 ± 0.3 - (2.4 ± 0.8) lo g 10 E n - (1.4 ± 0.4) lo g 10 $E^{2}_{n}$ for 240 Pu and TKE(MeV) = 177.1 ± 0.3 - (1.2 ± 0.9) lo g 10 E n - (1.8 ± 0.5)lo g 10 $E_{n}^{2}$ for 242 Pu.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Jacobian-scaled K-means clustering for physics-informed segmentation of reacting flows

This work introduces Jacobian-scaled K-means (JSK-means) clustering, which is a physicsinformed clustering strategy centered on the K-means framework. The method allows for the injection of underlying physical knowledge into the clustering procedure through a distance function modification: instead of leveraging conventional Euclidean distance vectors, the JSKmeans procedure operates on distance vectors scaled by matrices obtained from dynamical system Jacobians evaluated at the cluster centroids. The goal of this work is to show how the JSKmeans algorithm - without modifying the input dataset - produces clusters that capture regions of dynamical similarity, in that the clusters are redistributed towards high-sensitivity regions in phase space and are described by similarity in the source terms of samples instead of the samples themselves. The algorithm is demonstrated on a complex reacting flow simulation dataset (a channel detonation configuration), where the dynamics in the thermochemical composition space are known through the highly nonlinear and stiff Arrhenius-based chemical source terms. Interpretations of cluster partitions in both physical space and composition space reveal how JSK-means shifts clusters produced by standard K-means towards regions of high chemical sensitivity (e.g., towards regions of peak heat release rate near the detonation reaction zone). Furthermore, the findings presented here illustrate the benefits of utilizing Jacobian-scaled distances in clustering techniques, and the JSK-means method in particular displays promising potential for improving former partition-based modeling strategies in reacting flow (and other multi-physics) applications.

Clustering↗
DOE OSTI2024

NMPC for Mode-Switching Operation of Reversible Solid Oxide Cell Systems

Solid oxide cells (SOCs) are a promising dual-mode technology that generates hydrogen through high-temperature water electrolysis and generates power through a fuel cell reaction that consumes hydrogen. Reversible operation of SOCs requires a transition between these two modes for hydrogen production setpoints as the demand and price of electricity fluctuate. Moreover, a well-functioning control system is important to avoid cell degradation during mode-switching operation. In this work, we apply nonlinear model predictive control (NMPC) to an SOC module and supporting equipment and compare NMPC performance to classical proportional integral (PI) control strategies, while ramping between the modes of hydrogen and power production. While both control methods provide similar performance in many metrics, NMPC significantly reduces cell thermal gradients and curvatures (mixed spatial-temporal partial derivatives) during mode switching. A dynamic process flowsheet of the reversible SOC system was developed in the open-source, equation-based IDAES modeling framework. Our IDAES dynamic simulation results show that NMPC can ramp the SOC system between hydrogen and power production targets within short mode-switching times. Moreover, NMPC can comply with operating limits in the SOC system more effectively than PI, and only NMPC can directly enforce user-specified limits for mixed spatial-temporal partial derivatives of temperature. This allows for management of the trade-off be-tween operating efficiency and cell degradation, which is dependent on these temperature curvatures.

Li, Mingrui↗
NASA NTRS1974

Water impact analysis of space shuttle solid rocket motor by the finite element method

Preliminary analysis showed that the doubly curved triangular shell elements were too stiff for these shell structures. The doubly curved quadrilateral shell elements were found to give much improved results. A total of six load cases were analyzed in this study. The load cases were either those resulting from a static test using reaction straps to simulate the drop conditions or under assumed hydrodynamic conditions resulting from a drop test. The latter hydrodynamic conditions were obtained through an emperical fit of available data. Results obtained from a linear analysis were found to be consistent with results obtained elsewhere with NASTRAN and BOSOR. The nonlinear analysis showed that the originally assumed loads would result in failure of the shell structures. The nonlinear analysis also showed that it was useful to apply internal pressure as a stabilizing influence on collapse. A final analysis with an updated estimate of load conditions resulted in linear behavior up to full load.

Buyukozturk, O.↗
DOE OSTI2022

Synthesis, Crystal and Electronic Structures, Nonlinear Optical Properties, and Magnetic Properties of Two Thiophosphates: KInP 2 S 7 and KCrP 2 S 7

Two thiophosphates, KInP 2 S 7 and KCrP 2 S 7 , were structurally characterized without investigating any optical properties. Herein in this work, KInP 2 S 7 and KCrP 2 S 7 were revisited to investigate their optical and magnetic properties, respectively. Pure polycrystalline samples and crystals of KInP 2 S 7 and KCrP 2 S 7 were grown by high temperature solid state reactions, where mm-sized crystals of KCrP 2 S 7 were collected. KCrP 2 S 7 is isostructural to KInP 2 S 7 , which features a layered structure. KInP 2 S 7 and KCrP 2 S 7 possess close relationship to the layered thiophosphate M 2 P 2 S 6 (M = Fe, Co, Zn, etc.). The bonding pictures of KInP 2 S 7 were studied using the electron localization function (ELF) coupled with crystal orbital Hamilton population (COHP) calculations. The intrinsically distorted [PS 4 ] tetrahedra and [InS 6 ] octahedra are made by strong covalent P-S interactions and ionic In-S interactions, respectively. Electronic structure analysis confirmed that the optical properties of KInP 2 S 7 are mainly contributed to by [PS 4 ] tetrahedra together with small amounts of the contributions coming from [InS 6 ] octahedra. Magnetic measurement on mm-sized crystals of KCrP 2 S 7 verified that there is an antiferromagnetic transition around 21 K, and the Cr atoms are trivalent. KInP 2 S 7 is predicated to be an indirect bandgap semiconductor of 2.38 eV, which is confirmed by the UV-Vis measurement of 2.4(1) eV. KInP 2 S 7 is not a type-I phase-matching material and exhibits moderate second harmonic generation (SHG) response (0.51 × AgGaS 2 , sample of particle size of 100 µm). The laser damage threshold (LDT) of KInP 2 S 7 is very high of 5.2 × AgGaS 2 . Bandgap engineering were undergone to enhance the SHG response of KInP 2 S 7 .

36 MATERIALS SCIENCE↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Kotov, Dmitry V.↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗
DOE OSTI2021

Dynamic Catalysis Fundamentals: I. Fast calculation of limit cycles in dynamic catalysis

Dynamic catalysis—the forced oscillation of catalytic reaction coordinate potential energy surfaces (PES)—has recently emerged as a promising method for the acceleration of heterogeneously-catalyzed reactions. Theoretical study of enhancement of rates and supra-equilibrium product yield via dynamic catalysis has, to-date, been severely limited by onerous computational demands of forward integration of stiff, coupled ordinary differential equations (ODEs) that are necessary to quantitatively describe periodic cycling between PESs. Here, we establish a new approach that reduces, by ≳108×, the computational cost of finding the time-averaged rate at dynamic steady state (i.e. the limit cycle for linear and nonlinear systems of kinetic equations). Our developments are motivated by and conceived from physical and mathematical insight drawn from examination of a simple, didactic case study for which closed-form solutions of rate enhancement are derived in explicit terms of periods of oscillation and elementary step rate constants. Generalization of such closed-form solutions to more complex catalytic systems is achieved by introducing a periodic boundary condition requiring the dynamic steady state solution to have the same periodicity as the kinetic oscillations and solving the corresponding differential equations by linear algebra or Newton-Raphson-based approaches. The methodology is well-suited to extension to non-linear systems for which we detail the potential for multiple solutions or solutions with different periodicities. For linear and non-linear systems alike, the acute decrement in computational expense enables rapid optimization of oscillation waveforms and, consequently, accelerates understanding of the key catalyst properties that enable maximization of reaction rates, conversions, and selectivities during dynamic catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1983

Nonlinear aerodynamic modeling of flap oscillations in transonic flow - A numerical validation

The regime of validity of a nonlinear aerodynamic force and moment formulation, based on concepts from nonlinear functional analysis and applicable to a transonic airfoil with a deflecting flap, is investigated. A time-dependent finite difference technique is used to evaluate the aerodynamic data of the formulation in terms of specified, characteristic motions. Flap-motion histories are generated from the flap inertial equations of motion, with aerodynamic reactions specified by the moment formulation. The motion histories depicting the cases of decaying and growing flap oscillations are compared with histories generated through simultaneous, coupled solution of the fluid-dynamic equations and flap inertial equations of motion. The range of applicability of the formulation is discussed.

Chyu, W. J.↗
DOE OSTI2023

Quantum interferometry and pathway selectivity in the nonlinear response of photosynthetic excitons

We propose a time–frequency resolved spectroscopic technique which employs nonlinear interferometers to study exciton–exciton scattering in molecular aggregates. A higher degree of control over the contributing Liouville pathways is obtained as compared to classical light. We show how the nonlinear response can be isolated from the orders-of-magnitude stronger linear background by either phase matching or polarization filtering. Both arise due to averaging the signal over a large number of noninteracting, randomly oriented molecules. We apply our technique to the Frenkel exciton model which excludes charge separation for the photosystem II reaction center. We show how the sum of the entangled photon frequencies can be used to select two-exciton resonances, while their delay times reveal the single-exciton levels involved in the optical process.

Science & Technology - Other Topics↗
NASA NTRS1993

Inertial-space disturbance rejection for space-based manipulators

The implementation of a disturbance rejection controller for a 6-DOF PUMA manipulator mounted on a 3-DOF platform was described. A control algorithm is designed to track the desired position and attitude of the end-effector in inertial space, subject to unknown disturbances in the platform axes. Experimental results are presented for step, sinusoidal, and random disturbances in the platform rotational axis and in the neighborhood of kinematic singularities. Robotic manipulators were proposed as a means of reducing the amount of extra vehicular activity time required for space station assembly and maintenance. The proposed scenario involves a robotic manipulator attached to some mobile platform, such as a spacecraft, satellite, or the space station itself. Disturbances in the platform position and attitude may prevent the manipulator from successfully completing the task. The possibility of using the manipulator to compensate for platform disturbances was explored. The problem of controlling a robotic manipulator on a mobile platform has received considerable attention in the past few years. Joshi and Desrochers designed a nonlinear feedback control law to carry out tasks (with respect to the robot base frame) in the presence of roll, pitch and yaw disturbances in the platform axes. Dubowsky, Vance, and Torres proposed a time-optimal planning algorithm for a robotic manipulator mounted on a spacecraft, subject to saturation limits in the attitude control reaction jets. Papadopoulos and Dubowsky developed a general framework for analyzing the control of free-floating space manipulator systems. Most recently, Torres and Dubowsky have presented a technique called the enhanced disturbance map to find manipulator trajectories that reduce the effect of disturbances in the spacecraft position and attitude. One common assumption in the literature is that the disturbance signal is exactly known. If this is the case, then the end-effector location can be calculated without relying on direct end-point sensing. However, this assumption is invalid if there is a significant delay in the platform position and attitude measurements, or if the kinematics of the platform are not well known, or if the platform is a non-rigid structure. In the more likely case that only the nominal platform location and upper bound on the disturbance signal are known, direct end-point sensing is needed to measure the end-effector location.

Holt, Kevin↗
DOE OSTI2021

Ion Pairing Mediates Molecular Organization Across Liquid/Liquid Interfaces

Liquid/liquid interfaces play a central role in scientific fields ranging from nanomaterial synthesis and soft matter electronics to nuclear waste remediation and chemical separations. This diversity of functions arises from an interface’s ability to respond to changing conditions in its neighboring bulk phases. Understanding what drives this interfacial flexibility can provide novel avenues for designing new functional interfaces. However, limiting this progress is an inadequate understanding of the subtle intermolecular and interphase interactions taking place at the molecular level. Here, we use surface-specific vibrational sum frequency generation spectroscopy combined with atomistic molecular dynamics simulations to investigate the self-assembly and structure of model ionic oligomers consisting of an oligodimethylsiloxane (ODMS) tail covalently attached to a positively charged methyl imidazolium (MIM + ) head group at buried oil/aqueous interfaces. We show how the presence of seemingly innocuous salts can impart dramatic changes to the ODMS tail conformations in the oil phase via specific ion effects and ion-pairing interactions taking place in the aqueous phase. These specific ion interactions are shown to drive enhanced amphiphile adsorption, induce morphological changes, and disrupt emergent hydrogen-bonding structures at the interface. Tuning these interactions allows for independent control over the oligomer structure in the oil phase versus interfacial population changes and represents key mechanistic insight that is needed to control chemical reactions at liquid/liquid interfaces.

36 MATERIALS SCIENCE↗