Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spreading models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,405 records · Page 78

The influence of the directional energy distribution on the nonlinear dispersion relation in a random gravity wave field

The influence of the directional distribution of wave energy on the dispersion relation is calculated numerically using various directional wave spectrum models. The results indicate that the dispersion relation varies both as a function of the directional energy distribution and the direction of propagation of the wave component under consideration. Furthermore, both the mean deviation and the random scatter from the linear approximation increase as the energy spreading decreases. Limited observational data are compared with the theoretical results. The agreement is favorable.

Huang, N. E.↗

Modeling and technical use of gas evolving electrodes. Part 2: Modeling of gas-evolving electrolyzers with free electrolyte circulation

In an electrochemical reactor with gas-evolving electrodes, the transporting action of the gas bubbles can be used to move the electrolyte in a cycle flow, when the structure of the flow channels is suitable. For an electrolysis cell with such a circulation system a mathematic model was set up and evaluated. It is shown that in this manner, a rapid flow through the electrode gap can be achieved without additional energy consumption, in addition to a low gas fraction and a low cell voltage. The cell voltage and the attainable cycle spread are investigated as a function of the geometric parameters for their optimum values.

Schleiff, M.↗

Axial expansion of the electron pulse in a free-electron laser

The necessity of including longitudinal electrostatic fields in free electron laser theory is investigated. A simple model is proposed which examines the axial expansion of a free electron beam of finite length, and criteria are obtained for including macroscopic fields. It is shown that longitudinal Coulomb forces do not pose a problem in current free electron laser experiments. Transverse space-charge spreading is neglected, and negligible emittance is assumed.

Speziale, T.↗

Ridge-trench collision in Archean and Post-Archean crustal growth: Evidence from southern Chile

The growth of continental crust at convergent plate margins involves both continuous and episodic processes. Ridge-trench collision is one episodic process that can cause significant magmatic and tectonic effects on convergent plate margins. Because the sites of ridge collision (ridge-trench triple junctions) generally migrate along convergent plate boundaries, the effects of ridge collision will be highly diachronous in Andean-type orogenic belts and may not be adequately recognized in the geologic record. The Chile margin triple junction (CMTJ, 46 deg S), where the actively spreading Chile rise is colliding with the sediment-filled Peru-Chile trench, is geometrically and kinematically the simplest modern example of ridge collision. The south Chile margin illustrates the importance of the ridge-collision tectonic setting in crustal evolution at convergent margins. Similarities between ridge-collision features in southern Chile and features of Archean greenstone belts raise the question of the importance of ridge collision in Archean crustal growth. Archean plate tectonic processes were probably different than today; these differences may have affected the nature and importance of ridge collision during Archean crustal growth. In conclusion, it is suggested that smaller plates, greater ridge length, and/or faster spreading all point to the likelihood that ridge collision played a greater role in crustal growth and development of the greenstone-granite terranes during the Archean. However, the effects of modern ridge collision, and the processes involved, are not well enough known to develop specific models for the Archean ridge collison.

Nelson, E. P.↗

Effects of core turbulence on jet excitability

The effects of varying freestream core turbulence on the evolution of a circular jet with and without tonal excitation are examined. Measurements are made on an 8.8 cm diameter jet at a Mach number of 0.3. The jet is excitated by plane waves at Strouhal number 0.5. For the excited and unexcited cases the turbulence level is varied by screens and grids placed upstream of the nozzle exit. The experiment results are compared with a theoretical model which incorporates a variable core turbulence and considers the energy interactions between the mean flow, the turbulence and the forced component. Both data and theory indicate that increasing the freestream turbulence diminishes the excitability of the jet and reduces the effect of excitation on the spreading rate of the jet.

Mankbadi, Reda R.↗

Effects of core turbulence on jet excitability

The effects of varying freestream core turbulence on the evolution of a circular jet with and without tonal excitation are examined. Measurements are made on an 8.8 cm diameter jet at a Mach number of 0.3. The jet is excited by plane waves at Strouhal number 0.5. For the excited and unexcited cases the turbulence level is varied by screens and grids placed upstream of the nozzle exit. The experiment results are compared with a theoretical model which incorporates a variable core turbulence and considers the energy interactions between the mean flow, the turbulence and the forced component. Both data and theory indicate that increasing the freestream turbulence diminishes the excitability of the jet and reduces the effect of excitation on the spreading rate of the jet.

Mankbadi, Reda R.↗

Quickbird Satellite in-orbit Modulation Transfer Function (MTF) Measurement Using Edge, Pulse and Impulse Methods for Summer 2003

The spatial characteristics of an imaging system cannot be expressed by a single number or simple statement. However, the Modulation Transfer Function (MTF) is one approach to measure the spatial quality of an imaging system. Basically, MTF is the normalized spatial frequency response of an imaging system. The frequency response of the system can be evaluated by applying an impulse input. The resulting impulse response is termed the Point Spread function (PSF). This function is a measure of the amount of blurring present in the imaging system and is itself a useful measure of spatial quality. An underlying assumption is that the imaging system is linear and shift-independent. The Fourier transform of the PSF is called the Optical Transfer Function (OTF) and the normalized magnitude of the OTF is the MTF. In addition to using an impulse input, a knife-edge in technique has also been used in this project. The sharp edge exercises an imaging system at all spatial frequencies. The profile of an edge response from an imaging system is called an Edge Spread Function (ESF). Differentiation of the ESF results in a one-dimensional version of the Point Spread Function (PSF). Finally, MTF can be calculated through use of Fourier transform of the PSF as stated previously. Every image includes noise in some degree which makes MTF of PSF estimation more difficult. To avoid the noise effects, many MTF estimation approaches use smooth numerical models. Historically, Gaussian models and Fermi functions were applied to reduce the random noise in the output profiles. The pulse-input method was used to measure the MTF of the Landsat Thematic Mapper (TM) using 8th order even functions over the San Mateo Bridge in San Francisco, California. Because the bridge width was smaller than the 30-meter ground sample distance (GSD) of the TM, the Nyquist frequency was located before the first zero-crossing point of the sinc function from the Fourier transformation of the bridge pulse. To avoid the zero-crossing points in the frequency domain from a pulse, the pulse width should be less than the width of two pixels (or 2 GSD's), but the short extent of the pulse results in a poor signal-to-noise ratio. Similarly, for a high-resolution satellite imaging system such as Quickbird, the input pulse width was critical because of the zero crossing points and noise present in the background area. It is important, therefore, that the width of the input pulse be appropriately sized. Finally, the MTF was calculated by taking ratio between Fourier transform of output and Fourier transform of input. Regardless of whether the edge, pulse and impulse target method is used, the orientation of the targets is critical in order to obtain uniformly spaced sub-pixel data points. When the orientation is incorrect, sample data points tend to be located in clusters that result in poor reconstruction of the edge or pulse profiles. Thus, a compromise orientation must be selected so that all spectral bands can be accommodated. This report continues by outlining the objectives in Section 2, procedures followed in Section 3, descriptions of the field campaigns in Section 4, results in Section 5, and a brief summary in Section 6.

Helder, Dennis↗

Single-column Model Intercomparison for a Stably Stratified Atmospheric Boundary Layer

The parameterization of the stably stratified atmospheric boundary layer is a difficult issue, which has a large impact on the medium-range weather forecasts and on climate integrations. A non-strongly stratified arctic case is simulated by nineteen single-column turbulence schemes. The statistics from the Large-eddy simulation (LES) intercomparison made for the same case by eight different models are used as a guiding reference. The single-column parameterizations include research schemes and operational schemes from major forecast and climate research centres. First order schemes, a large number of turbulence kinetic energy closures, and other proposals have submitted results. There is a large spread in the results; in general, the operational schemes mix more efficiently than the research ones, and the TKE and other higher order closures give results closer to the LES statistics. The sensitivities of the schemes to the parameters of their turbulence closures are partially explored.

Cuxart, J.↗

Impact Assessment of Biomass Burning on Air Quality in Southeast and East Asia During BASE-ASIA

A synergy of numerical simulation, ground-based measurement and satellite observation was applied to evaluate the impact of biomass burning originating from Southeast Asia (SE Asia) within the framework of NASA's 2006 Biomass burning Aerosols in Southeast Asia: Smoke Impact Assessment (BASE-ASIA). Biomass burning emissions in the spring of 2006 peaked in MarcheApril when most intense biomass burning occurred in Myanmar, northern Thailand, Laos, and parts of Vietnam and Cambodia. Model performances were reasonably validated by comparing to both satellite and ground-based observations despite overestimation or underestimation occurring in specific regions due to high uncertainties of biomass burning emission. Chemical tracers of particulate K(+), OC concentrations, and OC/EC ratios showed distinct regional characteristics, suggesting biomass burning and local emission dominated the aerosol chemistry. CMAQ modeled aerosol chemical components were underestimated at most circumstances and the converted AOD values from CMAQ were biased low at about a factor of 2, probably due to the underestimation of biomass emissions. Scenario simulation indicated that the impact of biomass burning to the downwind regions spread over a large area via the Asian spring monsoon, which included Southern China, South China Sea, and Taiwan Strait. Comparison of AERONET aerosol optical properties with simulation at multi-sites clearly demonstrated the biomass burning impact via longrange transport. In the source region, the contribution from biomass burning to AOD was estimated to be over 56%. While in the downwind regions, the contribution was still significant within the range of 26%-62%.

biomass burning↗

Diffusional growth of cloud droplets in homogeneous isotropic turbulence: DNS, scaled-up DNS, and stochastic model

This paper presents a novel methodology to use direct numerical simulation (DNS) to study the impact of isotropic homogeneous turbulence on the condensational growth of cloud droplets. As shown by previous DNS studies, the impact of turbulence increases with the computational domain size, that is, with the Reynolds number, because larger eddies generate higher and longer-lasting supersaturation fluctuations that affect growth of individual cloud droplets. The traditional DNS can only simulate a limited range of scales because of the excessive computational cost that comes from resolving all scales involved, that is, from large scales at which the turbulent kinetic energy (TKE) is introduced down to the Kolmogorov microscale, and from following every single droplet. The novel approach is referred to as the “scaled-up DNS”. The scaling up is done in two parts, first by increasing both the computational domain and the Kolmogorov microscale and second by using super-droplets instead of real droplets. To ensure proper dissipation of TKE and scalar variance at small scales, molecular transport coefficients are appropriately scaled up with the grid length. For the scaled-up domains, say, meters and tens of meters, one needs to follow billions of real droplets. This is not computationally feasible, and so-called super-droplets are applied in scaled-up DNS simulations. Each super-droplet represents an ensemble of identical real droplets, and the number of real droplets represented by a super-droplet is referred to as the multiplicity attribute. After simple tests showing the validity of the methodology, scaled-up DNS simulations are conducted for five domains, the largest of 64 3 m 3 volume using a DNS of 256 3 grid points and various multiplicities. All simulations are carried out with vanishing mean vertical velocity and with no mean supersaturation, similarly to past DNS studies. As expected, the supersaturation fluctuations as well as the spread in droplet size distribution increase with the domain size, with the droplet radius variance increasing in time t as t 1∕2 as identified in previous DNS studies. Scaled-up simulations with different multiplicities document numerical convergence of the scaled-up solutions. Finally, we compare the scaled-up DNS results with a simple stochastic model that calculates supersaturation fluctuations based on the vertical velocity fluctuations updated using the Langevin equation. Overall, the results document similar scaling to previous small-domain DNS simulations and support the notion that the stochastic subgrid-scale model is a valuable tool for the multi-scale simulation of droplet spectral evolution applying a large-eddy simulation model.

54 ENVIRONMENTAL SCIENCES↗

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

Boundary layer transition↗

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition↗

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition↗

Hubble Space Telescope Images of Nearby Luminous Quasars: Results for Eight Quasars and Tests of the Detection Sensitivity - 2

Observations with the Wide-Field Camera of the Hubble Space Telescope (HST) are presented for eight intrinsically luminous quasars with redshifts between 0.16 and 0.29. These observations, when combined with a similar HST study of the quasar PKS 2349-014, show that luminous nearby quasars exist in a variety of environments. Seven companion galaxies brighter than M(V) = 16.5 (H(sub 0) = 100 km s(sup -1) Mpc(sup -1), Omega(sub 0) = 1.0) lie within a projected distance of 25 kpc of the quasars; three of the companions are located closer than 3'' (6 kpc projected distance) from the quasars, well within the volume that would be enclosed by a typical L* host galaxy. The observed association of quasars and companion galaxies is statistically significant and may he an important element in the luminous-quasar phenomenon. Apparent host galaxies are detected for three of the quasars: PG 1116+215, 3C 273, and PG 1444+407; the hosts have an average absolute magnitude of about 0.6 mag brighter than L*. The agreement between the previously published major-axis directions in ground-based images and in the present HST images of 3C 273 and PG 1444+407 constitutes important evidence supporting the reality of these candidate host galaxies. Upper limits are placed on the visual-band brightnesses of representative galactic hosts for all the quasars. These limits are established by placing galaxy images obtained with HST underneath the quasars and measuring at what faintness level the known galaxies are detected. On average, the HST spirals would have been detected if they were as faint as 1 mag below L*, and the early-type galaxies could have been detected down to a brightness level of about L*, where L* is the Schechter characteristic luminosity of field galaxies. Smooth, featureless galaxy models (exponential disks or de Vaucouleurs profiles) are fitted to the residual light after a best-fitting point source is subtracted from the quasar images. The results show that smooth spiral galaxies brighter than, on average, about L*, would have been detected. These upper limits, or possible detections, are consistent with, for example, the eight luminous quasars studied in this paper, occurring in host galaxies that have a Shechter luminosity function with a lower cutoff in the range 0.01-0.1 L*. Tests are performed to determine if our failure to detect, in some cases, luminous host galaxies could be an artifact caused by our analysis procedures. These tests include comparing the measured point-spread function (PSF) for our HST observations with the PSFs used in previous ground-based studies of host galaxies, measuring the fluctuations in the sky signals that were subtracted from the quasar images, evaluating empirically the effects of using different stellar PSFs in the analysis, carrying out the subtraction of the stellar (nuclear) source in different ways, creating and analyzing artificial active galactic nuclei (AGNs) with known surface brightnesses, and fitting the observed quasar light to an analytic model that includes a host galaxy.

Bahcall, John N.↗

A scale model wind tunnel study of dispersion in the Cleveland area. Laboratory simulation of lake breeze effects on diffusion from ground level emissions

A wind tunnel simulation of the diffusion patterns in a sea breeze was attempted. The results indicate that the low level onshore flow was well simulated for neutral, stable, unstable, and elevated inversion conditions. Velocity, turbulence, shear stress, and temperature data were taken, and the spread of emissions from ground level sources was investigated. Comparison is made with theoretical predictions by E. Inoue and with the open, homogeneous plane field results of Pasquill. Agreement with the predictions by Inoue is good, and the comparison with Pasquill's results shows that the wind tunnel flows are shifted two categories towards more stable. The discrepancy may be explained as a matter of averaging time.

Hoydysh, W. G.↗

A model for Io's atmosphere and sodium cloud

Observations of the sodium D-line emissions associated with the Galilean satellite Io suggest that the emitting region has two components. The first is an atmosphere which extends 1000 km above Io's surface. This is a relatively dense region which is optically thick in the sodium D lines. The second component is a surrounding optically thin cloud spread out along Io's orbit to form part of a torus. No emission was detected during any of three eclipses observed, although Io's sodium emission was consistently present when Io was near orbital elongation. This indicates that resonant scattering is the principal emission mechanism in both components. For the atmosphere, measurements of the D-line intensity, D-2/D-1 ratio, and line width indicate a temperature of 500 to 1000 K and a sodium column abundance of 3 by 10 to the 13th power per sq cm. Upper limits to the abundance of Mg, K, Ca, and Li indicate a high sodium abundance ratio in the atmosphere. This may also be true for Io's surface if sputtering by energetic particles is the source of the atmosphere. The sodium cloud has been detected as far as 60 arcsec from Io in the orbital plane and at least 6 arcsec above and below this plane. The cloud can be sustained by atoms escaping from the critical level at the top of the atmosphere.

Macy, W., Jr.↗

Future Atmospheric Rivers and Impacts on Precipitation: Overview of the ARTMIP Tier 2 High–Resolution Global Warming Experiment

Atmospheric rivers (ARs) are long, narrow synoptic scale weather features important for Earth’s hydrological cycle typically transporting water vapor poleward, delivering precipitation important for local climates. Understanding ARs in a warming climate is problematic because the AR response to climate change is tied to how the feature is defined. The Atmospheric River Tracking Method Intercomparison Project (ARTMIP) provides insights into this problem by comparing 16 atmospheric river detection tools (ARDTs) to a common data set consisting of high resolution climate change simulations from a global atmospheric general circulation model. ARDTs mostly show increases in frequency and intensity, but the scale of the response is largely dependent on algorithmic criteria. Across ARDTs, bulk characteristics suggest intensity and spatial footprint are inversely correlated, and most focus regions experience increases in precipitation volume coming from extreme ARs. The spread of the AR precipitation response under climate change is large and dependent on ARDT selection.

54 ENVIRONMENTAL SCIENCES↗

Investigating the Diurnal Radiative, Turbulent, and Biophysical Processes in the Amazonian Canopy‐Atmosphere Interface by Combining LES Simulations and Observations

Abstract We investigate the atmospheric diurnal variability inside and above the Amazonian rainforest for a representative day during the dry season. To this end, we combine high‐resolution large‐eddy simulations that are constrained and evaluated against a comprehensive observation set, including CO 2 concentrations, gathered during GoAmazon2014/15. We design systematic numerical experiments to quantify whether a multilayer approach in solving the explicit canopy improves our canopy‐atmosphere representation. We particularly focus on the relationship between photosynthesis and plant transpiration, and their distribution at leaf and canopy scales. We found the variability of photosynthesis drivers like vapor pressure deficit and leaf temperature to be about 3 times larger for sunlit leaves compared to shaded leaves. This leads to a large spread on leaf stomatal conductance values with minimum and maximum values varying more than 100%. Regarding the turbulent structure, we find wind‐driven stripe‐like shapes at the canopy top and structures resembling convective cells at the canopy. Wind‐related variables provide the best spatiotemporal agreement between model and observations. The potential temperature and heat flux profiles agree with an observed decoupling near the canopy top interface, although with less variability and cold biases of up to 3 K. The increasing complexity on the biophysical processes leads to the largest disagreements for evaporation, CO 2 plant assimilation and soil efflux. The model is able to capture the correct dependences and trends with the magnitudes still differing. We finally discuss the need to revise leaf and soil models and to complete the observations at leaf and canopy levels.

54 ENVIRONMENTAL SCIENCES↗