Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive transport 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Effect of reactive intermediates on the propagation of turbulent premixed flames

The most simple version of the Bray and Libby model of turbulent flame propagation has been extended to include a two-step reaction mechanism in which reactants decompose into reactive intermediates which then form products. A transport equation for the joint pdf of intermediate and product concentrations was derived, the pdf parameterized, and transport equations for the parameters obtained. The one-dimensional concentration equations were then solved for the parameters of the pdf and the turbulent kinetic energy. Flame speed and angle are recovered as eigenvalues of the problem. The results show the essential features observed in experiments. The product pdf is composed of peaks at zero and unity product concentration, with a continuous structure representing partially reacted gas. The intermediates pdf has a peak at zero concentration and a broad bell-shaped distribution at nonzero concentration, representing partially reacted gas. The results also show the effect of partitioning of energy release on the flame thickness and speed.

Azzazy, M.↗
DOE OSTI2023

Estimating source-sink distributions and fluxes of reactive nitrogen and sulfur within a mixed forest canopy

The vertical source-sink distribution of air pollutants within and above forested canopies is necessary for describing the biological, physical, and chemical processes influencing the soil-vegetation-atmosphere exchange. Here, this study implemented inverse modeling methods to estimate the source-sink and flux profiles of reactive nitrogen (N) and sulfur (S) compounds from measurements of the mean concentration profiles of ammonia (NH 3 ), nitric acid (HNO 3 ), sulfur dioxide (SO 2 ), and particulate ammonium (NH 4 + ), nitrate (NO 3 − ), and sulfate (SO 4 2− ) at a forest site in the southern Appalachian Mountains. Three inverse approaches utilizing different approximations to scalar transport within the canopy were developed and evaluated against sensible heat flux measurements. The Eulerian model (EUL), which incorporates vertical velocity skewness, performed well in reproducing the turbulent heat fluxes and was subsequently used to calculate the chemical source-sink and flux profiles. Above-canopy fluxes of NH 3 were downward, indicating that the forest was a net sink of NH 3 . The soil/litter layer was both a source and a sink for NH 3 but the exchange rate at the forest floor was small. Fluxes of HNO 3 , SO 2 , NO 3 − , NH 4 + , and SO 4 2- were uni-directional (deposition only) between the air and the canopy/ground and increased monotonically from the forest floor to the canopy top. Crown foliage dominated the uptake of reactive N and S during the growing season, accounting for 80–90% of the total canopy-scale flux. Fluxes and canopy-ground partitioning estimated using the resistance-based Surface Tiled Aerosol and Gas Exchange (STAGE) model were generally comparable to EUL. The comparison highlights the need for improved parameterizations of litter exchange and NH 3 compensation points in resistance models for forest ecosystems. The findings here benefit the application of critical loads in forest ecosystems and guide further development of resistance-based exchange models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Thermochemical reduction modeling in a high-temperature moving-bed reactor for energy storage: 1D model

The design of robust and efficient high-temperature thermochemical reactors and determination of operating conditions are critical steps toward enabling high-efficiency long-duration solar energy storage. This work presents a computational model for the thermal reduction of a metal oxide material (Mg-Mn-O) up to 1450 °C and the coupled complex transport phenomena in a novel tubular thermal reactor design that features the capability for a high extent-of-reduction (high energy storage density) and inherent heat recuperation. A one-dimensional model coupling counter-current gas–solid flow, two-phase heat transfer, thermochemical redox reactions, and species transport in a moving-bed reactor is developed. Simplified versions of the model are validated with published results in the literature for packed beds with both inert and reactive particles; the fully coupled model is also validated with experimental measurements of a moving-bed reactor in terms of local temperatures and oxygen release at the exit. Detailed comparisons on the effects of different boundary conditions in the reaction zone (prescribed wall temperature vs. heat flux conditions) and formulations based on a simple uniform flow assumption vs. plug flow using Ergun equation for gas flow are investigated. The results are compared with experimental measurements, and for all cases, the energy flow components in the reactor system and the thermal to chemical conversion efficiency and overall system efficiency are computed. Finally, the predicted high thermal-to-chemical efficiency ~95% and system efficiency ~30% agree with experimental measurements.

25 ENERGY STORAGE↗
DOE OSTI2023

UNF-ST&DARDS Enhancements for RCCA Data in As-loaded Dual Purpose Cask Models

This report summarizes the work performed to enable detailed modeling of rod cluster control assembly (RCCAs) in dual purpose canisters (DPCs) in as-loaded configurations using the Used Nuclear Fuel – Storage, Transportation & Disposal Analysis Resource and Data System (UNF-ST&DARDS). The goal of this project was to evaluate the reactivity impact that RCCAs have on k eff of DPCs with pressurized water reactor (PWR) fuel to potentially use the additional margin in future post-closure criticality safety analysis. This preliminary evaluation determined the number of DPCs that currently require compensatory actions prior to emplacement in a repository because of their high reactivity under post-closure criticality scenarios, which could be made acceptable by including the as-loaded RCCAs as specified in the Unified Database (UDB). This report briefly describes the modeling methods currently used within UNF-ST&DARDS and the modifications made to automate inclusion of the as-loaded RCCAs in the DPCs for post-closure criticality calculations for the loss of neutron absorber (NA) scenarios, or NA models. For the Zion site, the loss of basket scenarios, or degraded basket (DB) models, are also included. Discussion is also provided regarding the UDB data compared to available site-specific loading map data for the sites, which are specifically evaluated herein. This report compares the modifications to a similar evaluation for the Zion by site conducted by Walker as an initial validation. After the partial validation, all existing PWR sites with applicable NA models within the UDB were evaluated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Facilitated transport membranes for H 2 purification from coal-derived syngas: A techno-economic analysis

A single-stage membrane process has been designed for using facilitated transport membranes (FTMs) to decarbonize the coal-derived syngas from an integrated gasification combined cycle (IGCC) power plant. The necessary process model and costing method have also been developed to assess the technical feasibility and process economics. In order to account for the carrier saturation phenomenon associated with FTMs, a homogeneous reactive diffusion model is integrated into the process model. The techno-economic study reveals that the mitigated carrier saturation upon bulk CO 2 removal can lead to appreciable increases in the CO 2 permeance and CO 2 /H 2 selectivity, which can be utilized to achieve 95% CO 2 purity and 95% H 2 recovery with a CO 2 /H 2 selectivity of 50 at the complete carrier saturation. FTMs with different facilitated transport characteristics can also be arranged in a hybrid membrane configuration to render a H 2 recovery of 99% and a cost of electricity of $\$118.5$ /MWh, which is 12.5% lower than that of the benchmark Selexol process.

42 ENGINEERING↗
NASA NTRS2021

Global Impact of Lightning-Produced Oxidants

Lightning plays a major role in tropospheric oxidation, and its role on modulating tropospheric chemistry was thought to be emissions of nitrogen oxides (NOx). Recent field and laboratory measurements demonstrate that lightning generates extremely large amounts of oxidants, including hydrogen oxides (HOx) and O3. Here we implement these lightning-produced oxidants in a global chemical transport model to examine its global impact on tropospheric composition. We find that lightning-produced oxidants can increase global mass weighted OH by 0.3-10%, and affect CO, O3, and reactive nitrogen substantially, depending on the emission strength of oxidants from lightning. Our work highlights the importance and uncertainties of lightning-produced oxidants, as well as the need for rethinking the role of lightning in tropospheric oxidation chemistry.

Jingqiu Mao↗
DOE OSTI2021

Isotopic fractionation as in-situ sensor of subsurface reactive flow and precursor for rock failure

Greater utilization of subsurface reservoirs perturbs in-situ chemical-mechanical conditions with wide ranging consequences from decreased performance to project failure. Understanding the chemical precursors to rock deformation is critical to reducing the risks of these activities. To address this need, we investigated the coupled flow-dissolution- precipitation-adsorption reactions involving calcite and environmentally-relevant solid phases. Experimentally, we quantified (1) stable isotope fractionation processes for strontium during calcite nucleation and growth, and during reactive fluid flow; (2) consolidation behavior of calcite assemblages in the common brines. Numerically, we quantified water weakening of calcite using molecular dynamics simulations; and quantified the impact of calcite dissolution rate on macroscopic fracturing using finite element models. With microfluidic experiments and modeling, we show the effect of local flow fields on the dissolution kinetics of calcite. Taken together across a wide range of scales and methods, our studies allow us to separate the effects of reaction, flow, and transport, on calcite fracturing and the evolution of strontium isotopic signatures in the reactive fluids.

42 ENGINEERING↗
NASA NTRS2018

Decadal Changes in Summertime Reactive Oxidized Nitrogen and Surface Ozone over the Southeast United States

Widespread efforts to abate ozone (O3) smog have significantly reduced emissions of nitrogen oxides (NOx) over the past 2 decades in the Southeast US, a place heavily influenced by both anthropogenic and biogenic emissions. How reactive nitrogen speciation responds to the reduction in NOx emissions in this region remains to be elucidated. Here we exploit aircraft measurements from ICARTT (International Consortium for Atmospheric Research on Transport and Transformation - July-August 2004), SENEX (Southeast Nexus - June-July 2013), and SEAC4RS (Studies of Emissions and Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys - August-September 2013) and long-term ground measurement networks alongside a global chemistry-climate model to examine decadal changes in summertime reactive oxidized nitrogen (RON) and ozone over the Southeast US. We show that our model can reproduce the mean vertical profiles of major RON species and the total (NO (sub y)) in both 2004 and 2013. Among the major RON species, nitric acid (HNO3) is dominant (approximately 42-45 percent), followed by NOx (31 percent), total peroxy nitrates (Sigma PNs; 14 percent), and total alkyl nitrates (Sigma ANs; 9-12 percent) on a regional scale. We find that most RON species, including NOx, Sigma PNs, and HNO3, decline proportionally with decreasing NOx emissions in this region, leading to a similar decline in NO (sub y). This linear response might be in part due to the nearly constant summertime supply of biogenic VOC (Volatile Organic Compounds) emissions in this region. Our model captures the observed relative change in RON and surface ozone from 2004 to 2013. Model sensitivity tests indicate that further reductions of NOx emissions will lead to a continued decline in surface ozone and less frequent high-ozone events.

Li, Jingyi↗
DOE OSTI2025

Unraveling membrane electrode assembly design for electrochemical conversion of carbon dioxide to formate/formic acid

This work presents a one-dimensional continuum modeling approach to investigate various cell architectures used for electrochemical conversion of CO 2 to formate/formic acid. Ion transport is simulated by a system of generalized modified Poisson–Nernst–Planck (GMPNP) equations that reflect the reactive transport phenomena including steric effects as the electrolyte solutions become concentrated. In the cathode catalyst layer, ionic current contributions from both the supporting electrolyte and solid-state ionomer are considered. Voltage and CO 2 utilization breakdowns are utilized to deconvolute the impacts of the cell architecture. The origins of (bi)carbonate formation in the cathode are explored, as the subsequent decrease in CO 2 availability is a key reason for low faradaic efficiencies to formate/formic acid. In addition, the role of a supporting electrolyte (KOH) is investigated to understand its tradeoffs: while the K + ions can improve both conductivity and electrochemically active surface area in the cathode, the presence of OH − ions raises the pH and leads to deleterious formation of (bi)carbonates. To this end, we also present parametric studies on the concentration and flow rate of supplied KOH to the cell, to establish a path towards eliminating the need for a supporting electrolyte.

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

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2021

Chlamydomonas as a model for reactive oxygen species signaling and thiol redox regulation in the green lineage

While the maintenance of redox balance is important for all organisms, the energy balance in the chloroplast is particularly tightly linked to the overall cellular redox status in photosynthetic organisms. Excess light can cause overreduction of the photosynthetic electron transport chain, which, when unchecked, results in the generation of reactive oxygen species (ROS) that can oxidize cellular components including nucleic acids, proteins, and lipids (Li et al., 2009). A photosynthetic cell possesses multiple layers of defense against oxidative stress, such as antioxidant enzymes (Foyer and Noctor, 2005) and alternative electron transfer pathways (Alric and Johnson, 2017). At the same time, ROS can function as important signals that mediate responses to biotic and abiotic stresses (Apel and Hirt, 2004). Reversibly oxidizable small molecules and proteins can serve as ROS sinks and signals and are re-reduced by multiple mechanisms. For example, the ascorbate–glutathione (GSH) cycle is central in controlling ROS levels and regenerating/reducing oxidized cellular redox agents (Figure 1A; reviewed by Foyer and Noctor, 2011). The functions of many proteins can be regulated by the redox state of cysteine residues, such as the formation or reduction of disulfides by thioredoxins (TRX) (Buchanan, 1991) and glutaredoxins (Michelet et al., 2013b; Figure 1, B and C).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Pore-scale simulation of multiphase flow and reactive transport processes involved in geologic carbon sequestration

Multiphase flow and reactive transport are two essential physicochemical processes that govern the effectiveness of geological carbon sequestration (GCS). The interaction and feedback among different phases and components during intricate physicochemical processes hold great significance in understanding CO 2 sequestration. Pore-scale simulations can account for multiphase flow and reactive transport processes in porous media and obtain spatial distributions of parameters (density, velocity, concentration, etc.) in the pore space as well as their temporal evolutions. This proves especially valuable considering that experiments can be hindered by constraints in spatial and temporal resolution. The comprehensive insights garnered from pore-scale research can be leveraged for continuum modeling using the representative elementary volume (REV) concept. In this contribution, four sequential mechanisms of CO 2 -brine-rock interaction in three zones delineated by CO 2 saturation are elaborated to elucidate complicated physicochemical processes involved in GCS, which are followed by general descriptions of mathematical equations and pore-scale numerical methods. In addition, as interested and commonly encountered processes, leakage risks during GCS and CO 2 -enhanced oil recovery (CO 2 -EOR) processes are presented. The existing challenges and future directions are discussed for both the performance of the pore-scale models and the current gaps in the field of GCS. Importantly, we expect that this review will prove beneficial for researchers interested in pore-scale simulations, GCS, and related disciplines.

58 GEOSCIENCES↗
DOE OSTI2022

Design Considerations for Controlling Silicon Nanoparticle Nucleation and Growth in a Nonthermal Plasma

Controlling the nucleation and growth of nanoparticles in low temperature plasma systems is imperative for controlling nanoparticle size distributions; and for some applications such as deposition and etching in microelectronic processing, preventing particle contamination. In this work, silicon nanoparticle (NP) production from silane is used as a model system to investigate the nucleation process. Although the mechanisms responsible for silicon NP nucleation and growth have been studied, it is unclear how controllable system parameters (e.g., pressure, system geometry, and gas composition) can be used to inhibit or promote NP formation. For example, the transport of reactive silane species is expected to significantly affect the feed fraction of silane required to nucleate silicon NP (the nucleation onset fraction) due to losses at the reactor walls. In this work, NP mass density was determined as a function of system pressure, gas composition, and reactor diameter for a tubular flow-through radiofrequency (RF) plasma using Ar/H2/He/SiH4 gas mixtures. A quartz crystal microbalance (QCM) impactor was developed to measure the total aerosol mass density downstream of the plasma and thereby identify the nucleation onset and its dependence on process parameters. Finally, a reaction mechanism was developed and incorporated into a global plasma chemistry model to better understand the nucleation onset and NP growth.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Development of a three-dimensional APOLLO3 neutrons deterministic scheme for the CABRI reactor

CABRI is an experimental reactor to study the fuel behavior during reactivity injection transients. These transients being highly multiphysics, the development of suitable modeling and simulation tools to simulate them is important for the optimization of the tests and the control of the experimental conditions. This paper focuses on the development of an APOLLO3 deterministic core calculation tool dedicated to the CABRI transient analysis. It represents the first stage of the incremental process for the implementation of a multiphysics time-dependent modeling of the CABRI transient. The neutron calculation scheme is based on a classical two-step approach. The first step consists of a 281-energy group calculation flux with the TDT-MOC (Method Of Characteristics) solver for cross-section space and energy (23 groups) collapsing for the CABRI different assembly clusters. The bias on a 2D core neutron calculation due to the self-shielding calculation and collapsing on a restricted pattern are investigated thanks to a comparison with a direct full 2D calculation on a quarter of core. The second step relies on a pin-resolved transport 3D transport core calculation with the SN solver MINARET. A progressive numerical validation process is followed to quantify the calculation biases on reactivity and reaction rates at each step using reference calculations with the stochastic code TRIPOLI4. The next development stage toward a multiphysics scheme will be the implementation of the 3D-kinetics equation resolution and the coupling with a core thermal-hydraulics model. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Reactive chemical transport simulations of geologic carbon sequestration: Methods and applications

Chemical reaction simulations are considerably used to quantitatively assess the long-term geologic carbon sequestration (GCS), such as CO 2 sequestration capacity estimations, leakage pathway analyses, enhanced oil recovery (EOR) efficiency studies, and risk assessments of sealing formations (caprocks), wellbores, and overlying underground water resources. All these require a deep understanding of the CO 2 - associated chemical reactions. To ensure long-term, safe CO 2 sequestration in the intended formations, modeling is the only way to plausibly assess the CO 2 flow, reaction, and transport over thousands of years. This review summarizes the multiple methodologies for describing homogeneous and heterogeneous chemical reaction patterns and multiscale application examples, the recent progress and current status of chemical reaction simulations for GCS, and the impact of such simulations on geological CO2 sequestration performance. Technical gaps and future challenges are also discussed for further study. The trends and challenges of such studies include: (1) the combination of coupled chemical, mechanical, and transport processes with calibrated experiments and associated uncertainty/risk assessments; (2) enhancement of the ability to simulate detailed geophysical and geochemical equations to mimic in situ conditions; and (3) characterization of multiscale subsurface systems with detailed conceptual models and assignment of suitable boundary conditions for field-scale sequestration fields. One major gap remaining is the current lack of accurate (and scale-justified) kinetic and equilibrium chemical reaction parameters under reservoir conditions. Advanced models that couple chemical, mechanical, and transport processes with scale-justified parameters, from lab to field-scale experiments, are required for quantitative assessments of sequestration capacity and the long-term safety of GCS projects.

58 GEOSCIENCES↗
DOE OSTI2022

Effects of Chlorine Capture and a Proposed Density Law on the Reactivity of Plutonium Solution Systems

During fissionable material processing, all normal and credible abnormal conditions must remain safely subcritical. Nuclear Criticality Safety (NCS) uses a number of methods to determine subcriticality, one of which is the use of neutron transport codes such as MCNP6. In order to create models for use with MCNP6, both the geometry and materials in fissionable material processes must be known, or assumptions must be made and quantified for the impact to bias. One of the systems with a significant amount of bias due to material modeling assumptions is in the area of aqueous plutonium processing. These solutions are typically plutonium nitrate solutions or plutonium chloride solutions, which are modeled as fictitious plutonium metal-water mixtures because little is known about the actual density of the solution and there is no current predictive capability approved for use at Los Alamos National Laboratory (LANL) for modeling them. This research is currently underway to fill the gap and develop an algorithm for use with MCNP6 to model the density of plutonium chloride solutions. The method is to be validated with experimental data for density, and also validated with critical experiments using MCNP6. Note that the Chlorine Worth Study (CWS) was performed in December 2021 to help bridge the gap in chlorine data for critical experiments, and is currently awaiting International Criticality Safety Benchmark Evaluation Project (ICSBEP) review. This study was performed by LANL at the National Criticality Experiments Research Center (NCERC) at the Nevada National Security Site (NNSS). Additional information regarding this experiment may currently be found in LA-UR- 22-29180. Additionally, the Chemistry-Actinide Analytical Chemistry (C-AAC) at LANL has performed a number of solution density measurements for PuCl 3 -HC 1 -H 2 O, allowing for such data be used to create a semi-empirical density via the Pitzer method. The published dataset for the measurements is documented in LA-UR-22-25454. This method has already been tested successfully for aqueous plutonium nitrate solutions in SCALE. Current solution density measurements exist of plutonium concentrations of 0-~142g/L, all at 2M HC1, at temperatures 20-40°C. Additional data was taken for HC1-corrected density values, which essentially mimics the data for a pure PuCl x -water solution. The calculations in this report aim to support the current research by demonstrating the difference in system reactivity for the current modeling method when compared to the new proposed modeling with a density law implementation, which is being written as a Python tool to be used with MCNP6.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Sodium diffusion in heterogeneous porous media: Connecting laboratory experiments and simulations

Sodium has been suggested as a tracer for brine in reservoir formations where a significant amount of sodium ion is found mainly in the aqueous phase. Fortunately, the most abundant sodium isotope, 23 Na, is Nuclear Magnetic Resonance (NMR) active, making it possible to study the structure and dynamical properties of both bulk and pore confined fluid systems. We investigated the diffusion of key dissolved sodium species in bulk solution and porous media as a function of composition, ionic strength, and permeability of the carbonate samples using NMR relaxometry and pulsed-field gradient (PFG) NMR techniques. We use Indiana limestone as an example of natural porous media and water as a freely diffusible tracer and carrier. We demonstrate apparent diffusion measurement of sodium based on changes in spin–spin relaxation time (T 2 ) signal in pore confinement. The diffusion rate of sodium decreases with increasing counter-ion size. This effect is greater at higher ionic strengths and lower chemical potentials in porous media. The reactive transport code, CrunchFlow, was used to complement the NMR experiments to simulate diffusion behavior in porous media. Combining these two methods provides a powerful approach to estimating effective diffusion coefficients in heterogeneous matrices. The modeling considers the influence of physical properties (porosity and tortuosity) and chemical properties (geochemical composition and chemical gradients). The results presented in this work highlight the advantage of measuring apparent diffusivity using NMR T 2 relaxometry in conjunction with numerical simulation to derive effective diffusivity and the corresponding matrix properties (i.e., tortuosity) of the system.

58 GEOSCIENCES↗