Search NASA⌕ Search

SEARCH · Search NASA

Results for “dominant mode”

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

The MODIS Aerosol Algorithm, Products, Validation and Applications

The MODerate resolution Imaging Spectroradiometer (MODIS) currently aboard both the Terra and Aqua satellites produces a suite of products designed to characterize global aerosol distribution, optical thickness and particle size. Never before has a space-borne instrument been able to provide such detailed information, complementing field and modeling efforts to produce a comprehensive picture of aerosol characteristics. The three years of Terra-MODIS data have been validated by comparing with co-located AERONET observations of aerosol optical thickness and derivations of aerosol size parameters. Some 8000 comparison points located at 133 AERONET sites around the globe show that the MODIS aerosol optical thickness retrievals are accurate to within the pre-launch expectations. MODIS-derived size parameters are also compared with AERONET retrievals and found to agree well for fine-mode dominated aerosol regimes. Aerosol regimes dominated by dust aerosol are less accurate, attributed to what is thought to be nonsphericity. Errors due to nonsphericity will be reduced by introducing a new set of empirical phase functions, derived without any assumptions of particle shape. The major innovation that MODIS bring to the field of remote sensing of aerosol is the measure of particle size and the separation of finemode and coarsemode dominated aerosol regimes. Particle size can separate finemode man-made aerosols created during combustion, from larger natural aerosols originating from salt spray or wind erosion. This separation allows for the calculation of aerosol radiative effect and the estimation of the man-made aerosol radiative forcing. MODIS can also be used in regional studies of aerosol-cloud interaction that affect the global radiative and hydrological cycles.

Remer, L. A.↗

Crossflow-Vortex Breakdown on Swept Wings: Correlation of Nonlinear Physics

The spatial evolution of cross flow-vortex packets in a laminar boundary layer on a swept wing are computed by the direct numerical simulation of the incompressible Navier- Stokes equations. A wall-normal velocity distribution of steady suction and blowing at the wing surface is used to generate a strip of equally spaced and periodic disturbances along the span. Three simulations are conducted to study the effect of initial amplitude on the disturbance evolution, to determine the role of traveling cross ow modes in transition, and to devise a correlation function to guide theories of transition prediction. In each simulation, the vortex packets first enter a chordwise region of linear independent growth, then, the individual packets coalesce downstream and interact with adjacent packets, and, finally, the vortex packets nonlinearly interact to generate inflectional velocity profiles. As the initial amplitude of the disturbance is increased, the length of the evolution to breakdown decreases. For this pressure gradient, stationary modes dominate the disturbance evolution. A two-coeffcient function was devised to correlate the simulation results. The coefficients, combined with a single simulation result, provide sufficient information to generate the evolution pattern for disturbances of any initial amplitude.

Joslin, R. D.↗

Unsteady Analyses of Valve Systems in Rocket Engine Testing Environments

This paper discusses simulation technology used to support the testing of rocket propulsion systems by performing high fidelity analyses of feed system components. A generalized multi-element framework has been used to perform simulations of control valve systems. This framework provides the flexibility to resolve the structural and functional complexities typically associated with valve-based high pressure feed systems that are difficult to deal with using traditional Computational Fluid Dynamics (CFD) methods. In order to validate this framework for control valve systems, results are presented for simulations of a cryogenic control valve at various plug settings and compared to both experimental data and simulation results obtained at NASA Stennis Space Center. A detailed unsteady analysis has also been performed for a pressure regulator type control valve used to support rocket engine and component testing at Stennis Space Center. The transient simulation captures the onset of a modal instability that has been observed in the operation of the valve. A discussion of the flow physics responsible for the instability and a prediction of the dominant modes associated with the fluctuations is presented.

Shipman, Jeremy↗

Millimeter/Submillimeter Wave Communications Via Ceramic Ribbon

Finding very low loss waveguides in the millimeter/submillimeter wave range has been a problem of considerabel interest for many years. Starting from fundamentals, we have found a new way to design a waveguide structure which is capable of providing an attenuation coefficient of less than 10dB/km for the guided dominant mode.

Ceramic Ribbon millimeter wave waveguide wave rang↗

Curved Duct Noise Prediction Using the Fast Scattering Code

Results of a study to validate the Fast Scattering Code (FSC) as a duct noise predictor, including the effects of curvature, finite impedance on the walls, and uniform background flow, are presented in this paper. Infinite duct theory was used to generate the modal content of the sound propagating within the duct. Liner effects were incorporated via a sound absorbing boundary condition on the scattering surfaces. Simulations for a rectangular duct of constant cross-sectional area have been compared to analytical solutions and experimental data. Comparisons with analytical results indicate that the code can properly calculate a given dominant mode for hardwall surfaces. Simulated acoustic behavior in the presence of lined walls (using hardwall duct modes as incident sound) is consistent with expected trends. Duct curvature was found to enhance weaker modes and reduce pressure amplitude. Agreement between simulated and experimental results for a straight duct with hard walls (no flow) was excellent.

Dunn, M. H.↗

Multiscale Modeling of Intergranular Fracture in Aluminum: Constitutive Relation For Interface Debonding

Intergranular fracture is a dominant mode of failure in ultrafine grained materials. In the present study, the atomistic mechanisms of grain-boundary debonding during intergranular fracture in aluminum are modeled using a coupled molecular dynamics finite element simulation. Using a statistical mechanics approach, a cohesive-zone law in the form of a traction-displacement constitutive relationship, characterizing the load transfer across the plane of a growing edge crack, is extracted from atomistic simulations and then recast in a form suitable for inclusion within a continuum finite element model. The cohesive-zone law derived by the presented technique is free of finite size effects and is statistically representative for describing the interfacial debonding of a grain boundary (GB) interface examined at atomic length scales. By incorporating the cohesive-zone law in cohesive-zone finite elements, the debonding of a GB interface can be simulated in a coupled continuum-atomistic model, in which a crack starts in the continuum environment, smoothly penetrates the continuum-atomistic interface, and continues its propagation in the atomistic environment. This study is a step towards relating atomistically derived decohesion laws to macroscopic predictions of fracture and constructing multiscale models for nanocrystalline and ultrafine grained materials.

Yamakov, V.↗

Anatomy of the Binary Black Hole Recoil: A Multipolar Analysis

We present a multipolar analysis of the recoil velocity computed in recent numerical simulations of binary black hole coalescence, for both unequal masses and non-zero, non-precessing spins. We show that multipole moments up to and including 1 = 4 are sufficient to accurately reproduce the final recoil velocity (= 98%) and that only a few dominant modes contribute significantly to it (2 95%). We describe how the relative amplitude, and more importantly, the relative phase, of these few modes control the way in which the recoil builds up throughout the inspiral, merger, and ring-down phases. We also find that the numerical results can be reproduced, to a high level of accuracy, by an effective Newtonian formula for the multipole moments obtained by replacing in the Newtonian formula the radial separation with an effective radius computed from the numerical data. Beyond the merger, the numerical results are reproduced by a superposition of three Kerr quasi-normal modes. Analytic formulae, obtained by expressing the multipole moments in terms of the fundamental QNMs of a Kerr BH, are able to explain the onset and amount of '.anti-kick" for each of the simulations. Lastly, we apply this multipolar analysis to understand the remarkable difference between the amplitudes of planar and non-planar kicks for equal-mass spinning black holes.

Schnittman, Jeremy↗

Stability Analysis of Roughness Array Wake in a High-Speed Boundary Layer

Computations are performed to examine the effects of both an isolated and spanwise periodic array of trip elements on a high-speed laminar boundary layer, so as to identify the potential physical mechanisms underlying an earlier transition to turbulence as a result of the trip(s). In the context of a 0.333 scale model of the Hyper-X forebody configuration, the time accurate solution for an array of ramp shaped trips asymptotes to a stationary field at large times, indicating the likely absence of a strong absolute instability in the mildly separated flow due to the trips. A prominent feature of the wake flow behind the trip array corresponds to streamwise streaks that are further amplified in passing through the compression corner. Stability analysis of the streaks using a spatial, 2D eigenvalue approach reveals the potential for a strong convective instability that might explain the earlier onset of turbulence within the array wake. The dominant modes of streak instability are primarily sustained by the spanwise gradients associated with the streaks and lead to integrated logarithmic amplification factors (N factors) approaching 7 over the first ramp of the scaled Hyper-X forebody, and substantially higher over the second ramp. Additional computations are presented to shed further light on the effects of both trip geometry and the presence of a compression corner on the evolution of the streaks.

Choudhari, Meelan↗

Low-dimensional Representation of Error Covariance

Ensemble and reduced-rank approaches to prediction and assimilation rely on low-dimensional approximations of the estimation error covariances. Here stability properties of the forecast/analysis cycle for linear, time-independent systems are used to identify factors that cause the steady-state analysis error covariance to admit a low-dimensional representation. A useful measure of forecast/analysis cycle stability is the bound matrix, a function of the dynamics, observation operator and assimilation method. Upper and lower estimates for the steady-state analysis error covariance matrix eigenvalues are derived from the bound matrix. The estimates generalize to time-dependent systems. If much of the steady-state analysis error variance is due to a few dominant modes, the leading eigenvectors of the bound matrix approximate those of the steady-state analysis error covariance matrix. The analytical results are illustrated in two numerical examples where the Kalman filter is carried to steady state. The first example uses the dynamics of a generalized advection equation exhibiting nonmodal transient growth. Failure to observe growing modes leads to increased steady-state analysis error variances. Leading eigenvectors of the steady-state analysis error covariance matrix are well approximated by leading eigenvectors of the bound matrix. The second example uses the dynamics of a damped baroclinic wave model. The leading eigenvectors of a lowest-order approximation of the bound matrix are shown to approximate well the leading eigenvectors of the steady-state analysis error covariance matrix.

Tippett, Michael K.↗

Heat Entrapment Effects Within Liquid Acquisition Devices

We introduce a model problem to address heat entrapment effects or the local accumulation of thermal energy within liquid acquisition devices. We show that the parametric space consists of six parameters, namely the Rayleigh and Prandtl numbers, the aspect ratio, and heat flux ratios for the bottom, side, and top boundaries of the enclosure. For the range of Ra considered 1 to 10(sup 9), beyond Ra on the order of 10(sup 5), convective instability is the dominant mode of convection in comparison to natural convection. The flow field transitions to asymmetric modes at Ra on the order of 10(sup 7). Direct numerical simulation of a large geometric length scale prototype for Ra on the order of 10(sup 9) shows that the flow field evolves from small wavelength instability which gives rise to nonlinear growth of thermals, propagation of the instability occurs via growth of secondary and tertiary modes, and a travelling wave mode occurs prior to asymmetry. The effect of a large aspect ratio is to increase the number of modes in the vertical direction. Due to the slow diffusion of heat in the prototype, asymptotic states are not readily attained, we show that dynamical similarity can be used for a model which allows the attainment of asymptotic states and that transition to a chaotic state occurs for Ra on the order of 10(sup 9) via a broadband power spectrum. These dynamical events show that for the baseline condition in which heat is absorbed from background laboratory environment, higher heat flux is absorbed at the top and bottom boundaries of the enclosure than a nominal value of 34.9 ergs per square centimeter -second.

Duval, W. M. B.↗

Assessing Uncertainties in Boundary Layer Transition Predictions for HIFiRE-1 at Non-zero Angles of Attack

Boundary layer stability was analyzed for the HIFiRE-1 flight vehicle geometry for ground tests conducted at the CUBRC LENS I hypersonic shock test facility and the Langley Research Center (LaRC) 20- inch Mach 6 Tunnel. Boundary layer stability results were compared to transition onset location obtained from discrete heat transfer measurements from thin film gauges during the CUBRC test and spatially continuous heat transfer measurements from thermal phosphor paint data during the LaRC test. The focus of this analysis was on conditions at non-zero angles of attack as stability analysis has already been performed at zero degrees angle of attack. Also, the transition onset data obtained during flight testing was at nonzero angles of attack, so this analysis could be expanded in the future to include the results of the flight test data. Stability analysis was performed using the 2D parabolized stability software suite STABL (Stability and Transition Analysis for Hypersonic Boundary Layers) developed at the University of Minnesota and the mean flow solutions were computed using the DPLR finite volume Navier-Stokes computational fluid dynamics (CFD) solver. A center line slice of the 3D mean flow solution was used for the stability analysis to incorporate the angle of attack effects while still taking advantage of the 2D STABL software suite. The N-factors at transition onset and the value of Re(sub theta)/M(sub e), commonly used to predict boundary layer transition onset, were compared for all conditions analyzed. Ground test data was analyzed at Mach 7.2 and Mach 6.0 and angles of attack of 1deg, 3deg and 5deg. At these conditions, the flow was found to be second mode dominant for the HIFiRE-1 slender cone geometry. On the leeward side of the vehicle, a strong trend of transition onset location with angle of attack was observed as the boundary layer on the leeward side of the vehicle developed inflection points at streamwise positions on the vehicle that correlated to angle of attack. Inflection points are a strong instability mechanism that lead to rapid breakdown and transition to turbulence. The transition onset location on the windward side of the vehicle displayed no trend with angle of attack or freestream Reynolds number and transition was observed farther down the vehicle than observed on the leeward side of the vehicle. In analysis of both windward and leeward sides of the vehicle, use of the N factor methodology to develop trends to predict boundary layer transition onset showed improvements over the Re(sub theta)/M(sub e) empirical correlation methodology. Stronger correlations and less scatter in the data were observed when using the N factor method for these cases.

Marek, Lindsay C.↗

Modeling and Observations of the Response of Tropical Tropospheric Ozone to ENSO

The El Nino-Southern Oscillation (ENSO) is the dominant mode of tropical variability on interannual time scales. ENSO appears to extend its influence into the chemical composition of the tropical troposphere, Recent results have revealed an ENSO induced wave-1 anomaly in observed tropical tropospheric column ozone, This results in a dipole over the western and eastern tropical Pacific, whereby differencing the two regions produces an ozone anomaly with an extremely high correlation to the Nino 3.4 Index. We have successfully reproduced this result using the Goddard Earth Observing System Version 5 (GEOS-5) general circulation model coupled to a comprehensive stratospheric and tropospheric chemical mechanism forced with observed sea surface temperatures over the past 25 years, An examination of the modeled ozone field reveals the vertical contributions of tropospheric ozone to the column over the western and eastern Pacific region, We will show targeted comparisons with SHADOZ ozonesondes over these regions to provide insight into the vertical structure. Also, comparisons with NASA's Aura satellite Microwave Limb Sounder (MLS) and Tropospheric Emissions Spectrometer (TES) instruments and other appropriate data sets will be shown. In addition, the water vapor response to ENSO will be compared to help illuminate its role relative to dynamics in impacting ozone concentrations. These results indicate that the tropospheric ozone response to ENSO is potentially a very useful chemistry-climate diagnostic and should be considered in future modeling assessments.

Oman, L. D.↗

The Response of Tropospheric Ozone to ENSO in Observations and a Chemistry-Climate Simulation

The El Nino-Southern Oscillation (ENSO) is the dominant mode of tropical variability on interannual time scales. ENSO appears to extend its influence into the chemical composition of the tropical troposphere. Recent results have revealed an ENSO induced wave-l anomaly in observed tropical tropospheric column ozone. This results in a dipole over the western and eastern tropical Pacific, whereby differencing the two regions produces an ozone anomaly with an extremely high correlation to the Nino 3.4 Index. We have successfully reproduced this result using the Goddard Earth Observing System Version 5 (GEOS-5) general circulation model coupled to a comprehensive stratospheric and tropospheric chemical mechanism forced with observed sea surface temperatures over the past 25 years. An examination of the modeled ozone field reveals the vertical contributions of tropospheric ozone to the column over the western and eastern Pacific region. We will show targeted comparisons with observations from NASA's Aura satellite Microwave Limb Sounder (MLS), and the Tropospheric Emissions Spectrometer (TES) to provide insight into the vertical structure of ozone changes. The tropospheric ozone response to ENSO could be a useful chemistry-climate model evaluation tool and should be considered in future modeling assessments.

Oman, L. D.↗

El Nino-Southern Oscillation Correlated Aerosol Angstrom Exponent Anomaly Over the Tropical Pacific Discovered in Satellite Measurements

El Nino.Southern Oscillation (ENSO) is the dominant mode of interannual variability in the tropical atmosphere. ENSO could potentially impact local and global aerosol properties through atmospheric circulation anomalies and teleconnections. By analyzing aerosol properties, including aerosol optical depth (AOD) and Angstrom exponent (AE; often used as a qualitative indicator of aerosol particle size) from the Moderate Resolution Imaging Spectrometer, the Multiangle Imaging Spectroradiometer and the Sea ]viewing Wide Field ]of ]view Sensor for the period 2000.2011, we find a strong correlation between the AE data and the multivariate ENSO index (MEI) over the tropical Pacific. Over the western tropical Pacific (WTP), AE increases during El Nino events and decreases during La Nina events, while the opposite is true over the eastern tropical Pacific (ETP). The difference between AE anomalies in the WTP and ETP has a higher correlation coefficient (>0.7) with the MEI than the individual time series and could be considered another type of ENSO index. As no significant ENSO correlation is found in AOD over the same region, the change in AE (and hence aerosol size) is likely to be associated with aerosol composition changes due to anomalous meteorological conditions induced by the ENSO. Several physical parameters or mechanisms that might be responsible for the correlation are discussed. Preliminary analysis indicates surface wind anomaly might be the major contributor, as it reduces sea ]salt production and aerosol transport during El Nino events. Precipitation and cloud fraction are also found to be correlated with tropical Pacific AE. Possible mechanisms, including wet removal and cloud shielding effects, are considered. Variations in relative humidity, tropospheric ozone concentration, and ocean color during El Nino have been ruled out. Further investigation is needed to fully understand this AE ]ENSO covariability and the underlying physical processes responsible for it.

Li, Jing↗

Characterization of Deficiencies in the Frequency Domain Forced Response Analysis Technique for Supersonic Turbine Bladed Disks

Turbine blades in rocket and jet engine turbomachinery experience enormous harmonic loading conditions. These loads result from the integer number of upstream and downstream stator vanes as well as the other turbine stages. Assessing the blade structural integrity is a complex task requiring an initial characterization of whether resonance is possible and then performing a forced response analysis if that condition is met. The standard technique for forced response analysis in rocket engine turbines is to decompose a computational fluid dynamics (CFD).generated flow field into its harmonic components, and to then perform a frequency response analysis at the problematic natural frequencies using cyclically symmetric structural dynamic models. Recent CFD analysis and water-flow testing at NASA/MSFC, though, indicates that this technique may miss substantial harmonic and non ]harmonic excitation sources that become present in complex flows. This complex content can only be captured by a CFD flow field encompassing at least an entire revolution. A substantial development effort to create a series of software programs to enable application of the 360 degree forcing function in a frequency response analysis on cyclic symmetric models has been completed (to be described in a future paper), but the question still remains whether the frequency response analysis itself is capable of capturing the excitation content sufficiently. Two studies comparing frequency response analysis with transient response analysis, therefore, of bladed-disks undergoing this complex flow environment have been performed. The first is of a bladed disk with each blade modeled by simple beam elements and the disk modeled with plates (using the finite element code MSC/NASTRAN). The focus of this model is to be representative of response of realistic bladed disks, and so the dimensions are roughly equivalent to the new J2X rocket engine 1st stage fuel pump turbine. The simplicity of the model allows the CFD load to be able to be readily applied, along with analytical and experimental variations in both the temporal and spatial fourier components of the excitation. In addition, this model is a first step in identifying response differences between transient and frequency forced response analysis techniques. The second phase assesses this difference for a much more realistic solid model of a bladed-disk in order to evaluate the effect of the spatial variation in loading on blade dominated modes. Neither research on the accuracy of the frequency response method when used in this context or a comprehensive study of the effect of test-observed variation on blade forced response have been found in the literature, so this research is a new contribution to practical structural dynamic analysis of gas turbines. The primary excitation of the upstream nozzles interacts with the blades on fuel pump of the J2X causes the 5th Nodal diameter modes to be excited, as explained by Tyler and Sofrin1, so a modal analysis was first performed on the beam/plate model and the 5ND bladed-disk mode at 40167 hz was identified and chosen to be the one excited at resonance (see figure 1). The first forced response analysis with this model focuses on identifying differences between frequency and transient response analyses. A hypothesis going into the analysis was that perhaps the frequency response was enforcing a temporal periodicity that did not really exist, and so therefore it would overestimate the response. As high dynamic response was a considerable source of stress in the J2X, examining this concept could potentially be beneficial for the program.

Brown, Andrew M.↗

Analysis of Fine-Mode Aerosol Retrieval Capabilities by Different Passive Remote Sensing Instrument Designs

Remote sensing of aerosol optical properties is difficult, but multi-angle, multi-spectral, polarimetric instruments have the potential to retrieve sufficient information about aerosols that they can be used to improve global climate models. However, the complexity of these instruments means that it is difficult to intuitively understand the relationship between instrument design and retrieval success. We apply a Bayesian statistical technique that relates instrument characteristics to the information contained in an observation. Using realistic simulations of fine size mode dominated spherical aerosols, we investigate three instrument designs. Two of these represent instruments currently in orbit: the Multiangle Imaging SpectroRadiometer (MISR) and the POLarization and Directionality of the Earths Reflectances (POLDER). The third is the Aerosol Polarimetry Sensor (APS), which failed to reach orbit during recent launch, but represents a viable design for future instruments. The results show fundamental differences between the three, and offer suggestions for future instrument design and the optimal retrieval strategy for current instruments. Generally, our results agree with previous validation efforts of POLDER and airborne prototypes of APS, but show that the MISR aerosol optical thickness uncertainty characterization is possibly underestimated.

Knobelspiesse, Kirk↗

The Response of Lower Atmospheric Ozone to ENSO in Aura Measurements and a Chemistry-Climate Simulation

The El Nino-Southern Oscillation (ENSO) is the dominant mode of tropical variability on interannual time scales. ENSO appears to extend its influence into the chemical composition of the tropical troposphere. Recent work has revealed an ENSO-induced wave-1 anomaly in observed tropical tropospheric column ozone. This results in a dipole over the western and eastern tropical Pacific, whereby differencing the two regions produces an ozone anomaly with an extremely high correlation to the Nino 3.4 Index. We have successfully reproduced this feature using the Goddard Earth Observing System Version 5 (GEOS-5) general circulation model coupled to a comprehensive stratospheric and tropospheric chemical mechanism forced with observed sea surface temperatures over the past 25 years. An examination of the modeled ozone field reveals the vertical contributions of tropospheric ozone to the column over the western and eastern Pacific region. We will show composition sensitivity in observations from NASA s Aura satellite Microwave Limb Sounder (MLS) and the Tropospheric Emissions Spectrometer (TES) and a simulation to provide insight into the vertical structure of these ENSO-induced ozone changes. The ozone changes due to the Quasi-Biennial Oscillation (QBO) in the extra-polar upper troposphere and lower stratosphere in MLS measurements will also be discussed.

Oman, L. D.↗

Simultaneous Retrieval of Aerosol Properties and Clear-Sky Direct Radiative Effect Over the Global Ocean from MODIS

A unified satellite algorithm is presented to simultaneously retrieve aerosol properties (aerosol optical depth; AOD and aerosol type) and clear-sky shortwave direct radiative effect (hereafter, DREA) over ocean. The algorithm is applied to Moderate Resolution Imaging spectroradiometer (MODIS) observations for a period from 2003 to 2010 to assess the DREA over the global ocean. The simultaneous retrieval utilizes lookup table (LUT) containing both spectral reflectances and solar irradiances calculated using a single radiative transfer model with the same aerosol input data. This study finds that aerosols cool the top-of-atmosphere (TOA) and bottom-of-atmosphere (BOA) by 5.2 +/- 0.5 W/sq m and 8.3 W/sq m, respectively, and correspondingly warm the atmosphere (hereafter, ATM) by 3.1 W/sq m. These quantities, solely based on the MODIS observations, are consistent with those of previous studies incorporating chemical transport model simulations and satellite observations. However, the DREAs at BOA and ATM are expected to be less accurate compared to that of TOA due to low sensitivity in retrieving aerosol type information, which is related with the atmospheric heating by aerosols, particularly in low AOD conditions; consequently, the uncertainties could not be quantified. Despite the issue in the aerosol type information, the present method allows us to confine the DREA attributed only to fine-mode dominant aerosols, which are expected to be mostly anthropogenic origin, in the range from 1.1 W/sq m to 1.3 W/sq m at TOA. Improvements in size-resolved AOD and SSA retrievals from current and upcoming satellite instruments are suggested to better assess the DREA, particularly at BOA and ATM, where aerosol absorptivity induces substantial uncertainty.

MODIS↗