Search NASA⌕ Search

SEARCH · Search NASA

Results for “Temporal modeling”

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 649 records · Page 36

Machine Learning Surrogates of a Fuel Matrix Degradation Process Model for Performance Assessment of a Nuclear Waste Repository

Spent nuclear fuel repository simulations are currently not able to incorporate detailed fuel matrix degradation (FMD) process models due to their computational cost, especially when large numbers of waste packages breach. The current paper uses machine learning to develop artificial neural network and k-nearest neighbor regression surrogate models that approximate the detailed FMD process model while being computationally much faster to evaluate. Further, using fuel cask temperature, dose rate, and the environmental concentrations of CO 3 2- , O 2 , Fe 2+ , and H 2 as inputs, these surrogates show good agreement with the FMD process model predictions of the UO 2 degradation rate for conditions within the range of the training data. A demonstration in a full-scale shale repository reference case simulation shows that the incorporation of the surrogate models captures local and temporal environmental effects on fuel degradation rates while retaining good computational efficiency.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Deeper snow increases the net soil organic carbon accrual rate in moist acidic tussock tundra: 210Pb evidence from Arctic Alaska

The net change in the carbon inventory of arctic tundra remains uncertain as global warming leads to shifts in arctic water and carbon cycles. To better understand the response of arctic tundra carbon to changes in winter precipitation amount, we investigated soil depth profiles of carbon concentration and radionuclide activities ( 7 Be, 137 Cs, 210 Pb, and 241 Am) in the active layer of a twenty-two-year winter snow depth manipulation experiment in moist acidic tussock tundra at Toolik Lake, Alaska. Depth correlations of cumulative carbon dry mass (g cm -2 ) vs. unsupported 210 Pb activity (mBq g -1 ) were examined using a modified constant rate of supply (CRS) model. Results were best fit by two-slope CRS models indicating an apparent step temporal increase in the accumulation rate of soil organic carbon. Most of the best-fit model chronologies indicated that the increase in carbon accumulation rate apparently began and persisted after snow fence construction in 1994. The inhomogeneous nature of permafrost soils and their relatively low net carbon accumulation rates make it challenging to establish robust chronologic records. Nonetheless, the data obtained in this study support a decadal-scale increase in net soil organic carbon accumulation rate in the active layer of arctic moist acidic tussock tundra under conditions of increased winter precipitation.

54 ENVIRONMENTAL SCIENCES↗

The use of LANDSAT-4 MSS digital data in temporal data sets and the evaluation of scene-to-scene registration accuracy

The MSS sensor on LANDSAT 4 is, in certain performance aspects, different from those on LANDSATS 1 through 3. These differences created some concern in the NASA research community as to whether individual data sets can be registered accurately enough to produce acceptable data sets for multitemporal data analysis. The use of LANDSAT 4 MSS digital data in temporal data sets is examined and a method is presented for estimating temporal registration accuracy based on the use of an X-Y digitizer and grey tone electrostatic plots. Results indicate that the RMS temporal registration errors are not significantly different from the temporal data sets generated using LANDSAT 4 and LANDSAT 2 data (33.35 meters) and the temporal data set constructed from two LANDSAT 2 data sets (33.61 meters). A derivation of the model used to evaluate the temporal registration is included.

Anderson, J. E.↗

The use of Landsat-4 MSS digital data in temporal data sets and the evaluation of scene-to-scene registration accuracy

The MSS sensor on Landsat 4 is, in certain performance aspects, diferent from those of Landsats 1 through 3. These differences created some concern in the NASA research community as to whether individual data sets can be registered accurately enough to produce acceptable data sets for multitemporal data analysis. The use of Landsat 4 MSS digital data in temporal data sets is examined and a method is presented for estimating temporal registration accuracy based on the use of an X-Y digitizer and grey tone electrostatic plots. Results indicate that the RMS temporal registration errors are not significantly different from the temporal data sets generated using Landsat 4 and Landsat 2 data (33.35 meters) and the temporal data set constructed from two Landsat 2 data sets (33.61 meters). A derivation of the model used to evaluate the temporal registration is included.

Anderson, J. E.↗

The Radiation, Interplanetary Shocks, and Coronal Sources (RISCS) Toolset

We outline a plan to develop a physics based predictive toolset RISCS to describe the interplanetary energetic particle and radiation environment throughout the inner heliosphere, including at the Earth. To forecast and "nowcast" the radiation environment requires the fusing of three components: 1) the ability to provide probabilities for incipient solar activity; 2) the use of these probabilities and daily coronal and solar wind observations to model the 3D spatial and temporal heliosphere, including magnetic field structure and transients, within 10 AU; and 3) the ability to model the acceleration and transport of energetic particles based on current and anticipated coronal and heliospheric conditions. We describe how to address 1) - 3) based on our existing, well developed, and validated codes and models. The goal of RISCS toolset is to provide an operational forecast and "nowcast" capability that will a) predict solar energetic particle (SEP) intensities; b) spectra for protons and heavy ions; c) predict maximum energies and their duration; d) SEP composition; e) cosmic ray intensities, and f) plasma parameters, including shock arrival times, strength and obliquity at any given heliospheric location and time. The toolset would have a 72 hour predicative capability, with associated probabilistic bounds, that would be updated hourly thereafter to improve the predicted event(s) and reduce the associated probability bounds. The RISCS toolset would be highly adaptable and portable, capable of running on a variety of platforms to accommodate various operational needs and requirements.

Zank, G. P.↗

Optimal A-Train Data Utilization: A Use Case of Aura OMI L2G and MERRA-2 Aerosol Products

Ozone Monitoring Instrument (OMI) aboard NASA's Aura mission measures ozone column and profile, aerosols, clouds, surface UV irradiance, and the trace gases including NO2, SO2, HCHO, BrO, and OClO using UltraViolet electromagnetic spectrum (280 - 400 nm) with a daily global coverage and a pixel spatial resolution of 13 km × 24 km at nadir, and it's been one of the key instruments to study the Earth's atmospheric composition and chemistry. The second Modern-Era Retrospective analysis for Research and Applications (MERRA-2) is NASA's atmospheric reanalysis using an upgraded version of Goddard Earth Observing System Model, version 5 (GEOS-5) data assimilation system. Compared to its predecessor MERRA, MERRA-2 is enhanced with more aspects of the Earth system among which is aerosol assimilation. When comparing between satellite pixel measurements and modeled grid data, how to properly handle counterpart pairing is critical considering their spatial and temporal variations. The comparison between satellite and model data by simply using Level 3 (L3) products may result biases due to lack of detailed temporal information. It has been preferred to inter-compare or implement satellite derived physical quantity (i.e., Level 2 (L2) Swath type) directly with/to model measurements with higher temporal and spatial resolution as possible. However, this has posed a challenge in the community to handle. Rather than directly handling the L2 or L3 data, there is a Level 2G (L2G) product conserving L2 pixel scientific data quality but in Grid type with the global coverage. In this presentation, we would like to demonstrate the optimal utilization of OMI L2G daily aerosol products by comparing with MERRA-2 hourly aerosol simulations matched well in both space and time.

MERRA-2 reanalysis↗

Qualitative and temporal reasoning in engine behavior analysis

Numerical simulation models, engine experts, and experimental data are used to generate qualitative and temporal representations of abnormal engine behavior. Engine parameters monitored during operation are used to generate qualitative and temporal representations of actual engine behavior. Similarities between the representations of failure scenarios and the actual engine behavior are used to diagnose fault conditions which have already occurred, or are about to occur; to increase the surveillance by the monitoring system of relevant engine parameters; and to predict likely future engine behavior.

Dietz, W. E.↗

Trend Estimates of AERONET-Observed and Model-Simulated AOTs Between 1993 and 2013

Recently, temporal changes in Aerosol Optical Thickness (AOT) have been investigated based on model simulations, satellite and ground-based observations. Most AOT trend studies used monthly or annual arithmetic means that discard details of the generally right-skewed AOT distributions. Potentially, such results can be biased by extreme values (including outliers). This study additionally uses percentiles (i.e., the lowest 5%, 25%, 50%, 75% and 95% of the monthly cumulative distributions fitted to Aerosol Robotic Network (AERONET)-observed and ECHAM/MESSy Atmospheric Chemistry (EMAC)-model simulated AOTs) that are less affected by outliers caused by measurement error, cloud contamination and occasional extreme aerosol events. Since the limited statistical representativeness of monthly percentiles and means can lead to bias, this study adopts the number of observations as a weighting factor, which improves the statistical robustness of trend estimates. By analyzing the aerosol composition of AERONET-observed and EMAC-simulated AOTs in selected regions of interest, we distinguish the dominant aerosol types and investigate the causes of regional AOT trends. The simulated and observed trends are generally consistent with a high correlation coefficient (R = 0.89) and small bias (slope+/-2(sigma) = 0.75 +/- 0.19). A significant decrease in EMAC-decomposed AOTs by water-soluble compounds and black carbon is found over the USA and the EU due to environmental regulation. In particular, a clear reversal in the AERONET AOT trend percentiles is found over the USA, probably related to the AOT diurnal cycle and the frequency of wildfires. In most of the selected regions of interest, EMAC-simulated trends are mainly attributed to the significant changes of the dominant aerosols; e.g., significant decrease in sea salt and water soluble compounds over Central America, increase in dust over Northern Africa and Middle East, and decrease in black carbon and organic carbon over Australia.

AERONET↗

Modeling fluvial erosion on regional to continental scales

The fluvial system is a major concern in modeling landform evolution in response to tectonic deformation. Three stream bed types (bedrock, coarse-bed alluvial, and fine-bed alluvial) differ in factors controlling their occurrence and evolution and in appropriate modeling approaches. Spatial and temporal transitions among bed types occur in response to changes in sediment characteristics and tectonic deformation. Erosion in bedrock channels depends upon the ability to scour or pluck bed material; this detachment capacity is often a power function of drainage area and gradient. Exposure of bedrock in channel beds, due to rapid downcutting or resistant rock, slows the response of headwater catchments to downstream baselevel changes. Sediment routing through alluvial channels must account for supply from slope erosion, transport rates, abrasion, and sorting. In regional landform modeling, implicit rate laws must be developed for sediment production from erosion of sub-grid-scale slopes and small channels.

Howard, Alan D.↗

TBX3, a downstream target of TGF-β1, inhibits mesangial cell apoptosis

Chronic kidney disease (CKD) is an increasingly common condition characterized by progressive loss of functional nephrons leading to renal failure. TGF-β1-induced mesangial cell (MC) phenotype alterations have been linked to the genesis of CKD. Here we show that TGF-β1 regulates TBX3 gene expression in MC. This gene encodes for two main isoforms, TBX3.1 and TBX3+2α. TBX3.1 has been implicated in cell immortalization, proliferation and apoptosis by inhibiting p14{sup ARF}-Mdm2-p53 pathway, while TBX3+2α role has not been defined. We demonstrated that TBX3 overexpression abrogated MC apoptosis induced by serum deprivation. Moreover, we observed an enhancement in TBX3 protein expression both in glomerular and tubular regions in the model of 5/6 nephrectomy, temporally related to increased expression of TGF-β1, type IV collagen and fibronectin. Our results indicate that TBX3 acts as an anti-apoptotic factor in MC in vitro and may be involved in the mechanism by which TGF-β1 induces glomerulosclerosis and tubular fibrosis during the progression of nephropathies. - Highlights: • TBX3 isoforms are upregulated by TGF-b1 in mesangial cells. • TBX3 isoforms have different subcellular distribution profile in mesangial cells. • TBX3 isoforms exhibit antiapoptotic action in mesangial cells. • TBX3 protein is overexpressed in a model of nephropathy (5/6 nephrectomy)

60 APPLIED LIFE SCIENCES↗

One-dimensional hybrid-direct kinetic simulation of the discharge plasma in a Hall thruster

In order to model the non-equilibrium plasma within the discharge region of a Hall thruster, the velocity distribution functions (VDFs) must be obtained accurately. A direct kinetic (DK) simulation method that directly solves the plasma Boltzmann equation can achieve better resolution of VDFs in comparison to particle simulations, such as the particle-in-cell (PIC) method that inherently include statistical noise. In this paper, a one-dimensional hybrid-DK simulation, which uses a DK simulation for heavy species and a fluid model for electrons, is developed and compared to a hybrid-PIC simulation. Time-averaged results obtained from the hybrid-DK simulation are in good agreement with hybrid-PIC results and experimental data. It is shown from a comparison of using a kinetic simulation and solving the continuity equation that modeling of the neutral atoms plays an important role for simulations of the Hall thruster discharge plasma. In addition, low and high frequency plasma oscillations are observed. Although the kinetic nature of electrons is not resolved due to the use of a fluid model, the hybrid-DK model provides spatially and temporally well-resolved plasma properties and an improved resolution of VDFs for heavy species with less statistical noise in comparison to the hybrid-PIC method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Building a Prospective LCA Framework to Analyze Emerging Technologies in a Dynamic System Context

NREL's Lifecycle Analysis Integration into Opensource Numerical models (LiAISON) framework computes temporally explicit life cycle impacts and resource uses for specific technologies (foreground) in a dynamic system context (background). LiAISON computes results for a series of environmental mid-points enabling an analysis of prospective tradeoffs of emerging technologies toward 2100. This prospective feature is of critical importance when analyzing present-day emerging technologies whose large-scale impacts during deployment phases will occur in different, future system contexts. LiAISON systematically accounts for dynamic system changes by applying an integrated background of the future energy-economy-land-climate system, generated by exogenous integrated assessment models (IAMs). Using IAM scenarios, LiAISON generates a time-series of life cycle inventory (LCI) databases, which are then used to calculate the impacts per functional unit per time step. This expands current practice of using static, future system assumptions, e.g., a specific grid-mix each year. Further, IAM scenarios are provided in a standardized format of shared-socioeconomic pathways (SSP) and representative concentration pathways (RCP) combinations. These are coherent, regularly, published, and peer-reviewed scenario combinations that establish a reproducible and standardized societal and climate mitigation futures context. They are comparable across IAMs and expand the system boundary of the traditional LCA by including dimensions such as societal and behavioral changes. We apply the framework to assess two emerging Power-to-Hydrogen processes, high temperature electrolysis using solid oxide fuel cell (HT-SOE) and polymer electrolyte membrane electrolysis (PEME). We compare the technologies to a baseline Hydrogen production process via steam methane reforming. Despite the decarbonized electricity systems' beneficial effects on the PtH2 processes' carbon intensities, we find environmental tradeoffs, which require technology improvements via learning-by-doing to be alleviated. Future work via ongoing collaborations will focus on linking the framework to other energy-economy-land-climate models and open-source life cycle inventory databases.

ENERGY PLANNING, POLICY, AND ECONOMY↗

Evaluation of preindustrial to present-day black carbon and its albedo forcing from Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP)

As part of the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), we evaluate the historical black carbon (BC) aerosols simulated by 8 ACCMIP models against observations including 12 ice core records, long-term surface mass concentrations, and recent Arctic BC snowpack measurements. We also estimate BC albedo forcing by performing additional simulations using offline models with prescribed meteorology from 1996-2000. We evaluate the vertical profile of BC snow concentrations from these offline simulations using the recent BC snowpack measurements. Despite using the same BC emissions, the global BC burden differs by approximately a factor of 3 among models due to differences in aerosol removal parameterizations and simulated meteorology: 34 Gg to 103 Gg in 1850 and 82 Gg to 315 Gg in 2000. However, the global BC burden from preindustrial to present-day increases by 2.5-3 times with little variation among models, roughly matching the 2.5-fold increase in total BC emissions during the same period.We find a large divergence among models at both Northern Hemisphere (NH) and Southern Hemisphere (SH) high latitude regions for BC burden and at SH high latitude regions for deposition fluxes. The ACCMIP simulations match the observed BC surface mass concentrations well in Europe and North America except at Ispra. However, the models fail to predict the Arctic BC seasonality due to severe underestimations during winter and spring. The simulated vertically resolved BC snow concentrations are, on average, within a factor of 2-3 of the BC snowpack measurements except for Greenland and the Arctic Ocean. For the ice core evaluation, models tend to adequately capture both the observed temporal trends and the magnitudes at Greenland sites. However, models fail to predict the decreasing trend of BC depositions/ice core concentrations from the 1950s to the 1970s in most Tibetan Plateau ice cores. The distinct temporal trend at the Tibetan Plateau ice cores indicates a strong influence from Western Europe, but the modeled BC increases in that period are consistent with the emission changes in Eastern Europe, the Middle East, South and East Asia. At the Alps site, the simulated BC suggests a strong influence from Europe, which agrees with the Alps ice core observations. At Zuoqiupu on the Tibetan Plateau, models successfully simulate the higher BC concentrations observed during the non-monsoon season compared to the monsoon season but overpredict BC in both seasons. Despite a large divergence in BC deposition at two Antarctic ice core sites, some models with a BC lifetime of less than 7 days are able to capture the observed concentrations. In 2000 relative to 1850, globally and annually averaged BC surface albedo forcing from the offline simulations ranges from 0.014 to 0.019Wm−2 among the ACCMIP models. Comparing offline and online BC albedo forcings computed by some of the same models, we find that the global annual mean can vary by up to a factor of two because of different aerosol models or different BC-snow parameterizations and snow cover. The spatial distributions of the offline BC albedo forcing in 2000 show especially high BC forcing (i.e., over 0.1W/sq. m) over Manchuria, Karakoram, and most of the Former USSR. Models predict the highest global annual mean BC forcing in 1980 rather than 2000, mostly driven by the high fossil fuel and biofuel emissions in the Former USSR in 1980.

snow cover↗

Ionospheric response to the sustained high geomagnetic activity during the March '89 great storm

A simulation was conducted to model the high-latitude ionospheric to the sustaied level of high geomagnetic activity for the great magnetic storm period of March 13-14, 1989. The geomagnetic and solar activity indices and the Defense Meterological Satellite Program (DMSP) F8 and F9 satellite data for particle precipitation and high-latitude convection were used as inputs to a time-dependent ionospheric model (TDIM). The results of the TDIM were compared to both DMSP plasma density data and ground-based total electron content (TEC) measurements for the great storm period as well as with earlier storm observations. The comparisons show that the overall structure of the high-latitude ionosphere was dominated by an increased convection speed within the polar cap that led to increased ion temperatures. In turn, this enhanced the NO(+) density, raised the atomic-to-molecular ion transition height to over 300 km, decreased N(sub m)F(sub 2), increased h(sub m)F(sub 2), and in places either increased n(sub e) at 800 km or slightly decreased it. The morphology of the ionosphere under these extreme conditions was considerably different than that modeled for less distributed intervals. These differences included the character of the dayside tongue of ionization that no longer extended deep into the polar cap. Instead, as a result of the ion heating and consequent reduction in N(sub m)F(sub 2), a large polar hole occupied much of the polar region. This polar hole extended beyond the auroral oval and merged with the night sector midatitude trough. The limitaions associated with the applicability of the TDIM to the geomagnetic conditions present on March 13 and 14 are discussed. The primary limitations of the TDIM derive from the limited temporal resolution of the model input parameters and the lack of suitably dynamic thermospheric specification for the great storm conditions. These limitations leads to midlatitude ionospheric storm phases that do no follow those observed.

Sojka, J. J.↗

System and Method for Finite Element Simulation of Helicopter Turbulence

The present invention provides a turbulence model that has been developed for blade-element helicopter simulation. This model uses an innovative temporal and geometrical distribution algorithm that preserves the statistical characteristics of the turbulence spectra over the rotor disc, while providing velocity components in real time to each of five blade-element stations along each of four blades. for a total of twenty blade-element stations. The simulator system includes a software implementation of flight dynamics that adheres to the guidelines for turbulence set forth in military specifications. One of the features of the present simulator system is that it applies simulated turbulence to the rotor blades of the helicopter, rather than to its center of gravity. The simulator system accurately models the rotor penetration into a gust field. It includes time correlation between the front and rear of the main rotor, as well as between the side forces felt at the center of gravity and at the tail rotor. It also includes features for added realism, such as patchy turbulence and vertical gusts in to which the rotor disc penetrates. These features are realized by a unique real time implementation of the turbulence filters. The new simulator system uses two arrays one on either side of the main rotor to record the turbulence field and to produce time-correlation from the front to the rear of the rotor disc. The use of Gaussian Interpolation between the two arrays maintains the statistical properties of the turbulence across the rotor disc. The present simulator system and method may be used in future and existing real-time helicopter simulations with minimal increase in computational workload.

McFarland, R. E.↗

Pressure Sensitive Paint Applied to Flexible Models Project

One gap in current pressure-measurement technology is a high-spatial-resolution method for accurately measuring pressures on spatially and temporally varying wind-tunnel models such as Inflatable Aerodynamic Decelerators (IADs), parachutes, and sails. Conventional pressure taps only provide sparse measurements at discrete points and are difficult to integrate with the model structure without altering structural properties. Pressure Sensitive Paint (PSP) provides pressure measurements with high spatial resolution, but its use has been limited to rigid or semi-rigid models. Extending the use of PSP from rigid surfaces to flexible surfaces would allow direct, high-spatial-resolution measurements of the unsteady surface pressure distribution. Once developed, this new capability will be combined with existing stereo photogrammetry methods to simultaneously measure the shape of a dynamically deforming model in a wind tunnel. Presented here are the results and methodology for using PSP on flexible surfaces.

model deformation↗

Ensemble cure kinetics network (ECK-Net): A method to derive cure kinetics of thermosetting resin

This paper introduces an Ensemble Cure Kinetics Network (ECK-Net), a neural network (NN)–based framework for modeling the cure kinetics of thermosetting resins within a phenomenological context. ECK-Net replaces traditional analytic models, which require extensive chemical insight and multiple isothermal/non-isothermal experiments, with a data-driven surrogate that maps nonlinear relationships between temperature, degree of cure, and reaction rate from differential scanning calorimetry data. The proposed approach predicts input-dependent kinetic coefficients of a generalized nth-order reaction equation rather than reaction rates directly, enabling a single unified model to represent various epoxy systems without relying on iso-conversional analysis or predefined functional forms. To ensure robustness, multiple independently trained networks under different random initializations are blended through an ensemble strategy, effectively mitigating the stochastic variability inherent to neural networks. The framework is validated using experimental datasets from multiple resin systems, including aerospace-grade materials (Toray 3900-2, Cycom 5320-1, and Hexcel 8552) and a windmill-grade resin (RIMR 035c). The model accurately reproduces the temporal evolution of the degree of cure under manufacturers’ recommended cure cycles across all tested resins systems, yielding Pearson’s correlation coefficients of 0.992, 0.994, 0.993, 0.997, respectively. To demonstrate process-level applicability, the trained network was implemented within the Abaqus environment to simulate out-of-autoclave (OOA) curing process of the CFRP panel composed of Toray T830H-6K/3900-2D prepreg. The simulation results showed excellent agreement with experimental temperature response (maximum peak temperature, simulation: 189.6 °C, experiment: 188.5 °C) and the final degree of cure (simulation: 0.948, experiment: 0.960 ± 0.013), confirming ECK-Net’s capability as a reliable alternative to conventional cure kinetics modeling methods.

Composite curing↗

Machine Learning Design of Perovskite Catalytic Properties

Abstract Discovering new materials that efficiently catalyze the oxygen reduction and evolution reactions is critical for facilitating the widespread adoption of solid oxide fuel cell and electrolyzer (SOFC/SOEC) technologies. Here, machine learning (ML) models are developed to predict perovskite catalytic properties critical for SOFC/SOEC applications, including oxygen surface exchange, oxygen diffusivity, and area specific resistance (ASR). The models are based on trivial‐to‐calculate elemental features and are more accurate and dramatically faster than the best models based on ab initio‐derived features, potentially eliminating the need for ab initio calculations in descriptor‐based screening. The model of ASR enables temperature‐dependent predictions, has well calibrated uncertainty estimates and online accessibility. Use of temporal cross‐validation reveals the model to be effective at discovering new promising materials prior to their initial discovery, demonstrating the model can make meaningful predictions. Using the SHapley Additive ExPlanations (SHAP) approach, detailed discussion of different approaches of model featurization is provided for ML property prediction. Finally, the model is used to screen more than 19 million perovskites to develop a list of promising cheap, earth‐abundant, stable, and high performing materials, and find some top materials contain mixtures of less‐explored elements (e.g., K, Bi, Y, Ni, Cu) worth exploring in more detail.

25 ENERGY STORAGE↗