Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RANS”

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,117 records · Page 62

Slow Solar Wind Connection Science during Solar Orbiter's First Close Perihelion Passage

The Slow Solar Wind Connection Solar Orbiter Observing Plan (Slow Wind SOOP) was developed to utilize the extensive suite of remote-sensing and in situ instruments on board the ESA/NASA Solar Orbiter mission to answer significant outstanding questions regarding the origin and formation of the slow solar wind. The Slow Wind SOOP was designed to link remote-sensing and in situ measurements of slow wind originating at open–closed magnetic field boundaries. The SOOP ran just prior to Solar Orbiter's first close perihelion passage during two remote-sensing windows (RSW1 and RSW2) between 2022 March 3–6 and 2022 March 17–22, while Solar Orbiter was at respective heliocentric distances of 0.55–0.51 and 0.38–0.34 au from the Sun. Coordinated observation campaigns were also conducted by Hinode and IRIS. The magnetic connectivity tool was used, along with low-latency in situ data and full-disk remote-sensing observations, to guide the target pointing of Solar Orbiter. Solar Orbiter targeted an active region complex during RSW1, the boundary of a coronal hole, and the periphery of a decayed active region during RSW2. Postobservation analysis using the magnetic connectivity tool, along with in situ measurements from MAG and SWA/PAS, showed that slow solar wind originating from two out of three of the target regions arrived at the spacecraft with velocities between ∼210 and 600 km s −1 . The Slow Wind SOOP, despite presenting many challenges, was very successful, providing a blueprint for planning future observation campaigns that rely on the magnetic connectivity of Solar Orbiter.

Slow solar wind↗

Cryogenic Extension of NASA Species Polynomials Using Hydrogen and Oxygen at Stoichiometry

NASA has been conducting research in Rotating Detonation Rocket Engines (RDREs) for several years. A recent test at the Marshall Spaceflight Center has successfully shown operation of a RDRE that ran for 251 seconds which had combustion chamber cryogenic inlet conditions [1]. However, an overall pressure gain was not observed in the performance of the engine. Such an indication could mean an optimal chamber design is yet to be discovered. Such a design could significantly reduce parasitic pressure loses and improve overall engine efficiency beyond that of conventional rockets. To investigate different chamber configurations, new CFD capability is needed to accurately capture cryogenic thermophysical and transport property data. Such an effort is currently part of NASA’s Early Career Initiative (ECI) program. This paper focuses on the thermophysical property mixture model approach that makes use of NASA’s polynomial fits of species and another cryogenic model chosen from NIST’s Refprop program. Presented in this paper are results of a 1D denotation CFD simulation of stoichiometric Hydrogen and Oxygen mixture at a cryogenic upstream condition.

Combustion↗

Quantifying the Ability of Magnetohydrodynamic Models to Reproduce Observed Birkeland Current and Auroral Electrojet Magnitudes

Although global magnetohydrodynamic (MHD) models have increased in sophistication and are now at the forefront of modeling Space Weather, there is still no clear understanding of how well these models replicate the observed ionospheric current systems. Without a full understanding and treatment of the ionospheric current systems, global models will have significant shortcomings that will limit their use. In this study we focus on reproducing observed seasonal interhemispheric asymmetry in ionospheric currents using the Space Weather Modeling Framework (SWMF). We find that SWMF does reproduce the linear relationship between the electrojets and the FACs, despite the underestimation of the currents’ magnitudes. Quantitatively, we find that at best SWMF is only capturing approximately 60% of the observed current. We also investigate how varying F10.7 effects the ionospheric potential and currents during the summer and winter. We find that simulations ran with higher F10.7 result in lower ionospheric potentials. Additionally, we find that the models do not always replicate the expected behavior of the currents with varying F10.7. This work points to a needed improvement in ionospheric conductance models.

Tre'Shunda James↗

Heated Supersonic Axisymmetric Jet Cases for the NASA Turbulence Modeling Resource

This paper will describe three new supersonic jet cases for the NASA Turbulence Modeling Resource (TMR). These cases were taken from the 6th AIAA Propulsion Aerodynamics Workshop (PAW) nozzle test problem, which utilized data taken at NASA GRC consisting of non-intrusive measurements of jet plume velocities and temperatures, including both mean values and turbulence statistics. The three jet cases all used a Mach 1.63 axisymmetric nozzle: A heated jet for each at on-design conditions, a temperature-matched on-design condition at Mach 1.63, and an off-design heated condition for the Mach 1.63 nozzle. Computational fluid dynamics (CFD) solutions are obtained using three established Reynolds-averaged Navier-Stokes (RANS) codes: Wind-US, FUN3D, and Vulcan, all using the Menter Shear Stress Transport (SST) k-ω turbulence model. The NASA TMR requires at least two CFD codes yielding essentially identical results to certify that the approaches are verified, that is - providing solutions that may not necessarily provide close comparisons with experimental data, but solve the posed CFD equations as intended.

nozzle↗

NASA HECC Geometry and Performance Review Part 1: Validation of a Computational Model for the Vaneless Diffuser Configuration with As-Manufactured Impeller Geometry

An investigation of the NASA High Efficiency Centrifugal Compressor (HECC) vaneless diffuser configuration was performed. This multipart investigation focused on validating the computational model against experimental data. The validated model was then used to explore the effects of observed geometric differences between the as-manufactured and design-intent impellers. Lastly, the validated model was used for comparison against an experimental tip clearance study. Details regarding each of the parts are detailed below: A multipart investigation of the NASA High Efficiency Centrifugal Compressor (HECC) vaneless diffuser configuration was performed. Part I of the investigation focused on describing the HECC vaneless computational model and its validation against the experimental data published in GT2023-103128. To generate the model, a verification of the manufactured impeller was performed and was confirmed to match the impeller CAD. Section data from the impeller CAD, referred to as As-Manufactured, were used in the model. Steady RANS simulations using a multiblock structured mesh were used to conduct a mesh convergence study. The resultant mesh was then used for simulations at both design-speed and off-design speed conditions. The simulations compared well to data in both 1D and 2D terms. Based on the observed results the model was deemed validated and was used in the latter two parts of the multipart investigation.

centrifugal compressor↗

Multidisciplinary Optimization of a Transonic Truss Braced Wing Aircraft with Hybrid-Electric Propulsion

Next generation aircraft concepts have subsystems that are increasingly inter-connected. This necessitates advanced design tools using gradient based optimization to properly account for strong subsystem coupling in these aircraft. Additionally, these design tools require accurate modeling capability to support the increased interest in electrified propulsion. Design studies for these aircraft must consider the trade-offs between electric propulsion and turbojet engines at all points in the flight envelope to determine the optimal balance between propulsion options. This paper describes the development of an electric propulsion subsystem model designed to work within Aviary - an open source aircraft design tool. We will demonstrate the electric propulsion model operating with Aviary by showing results from an optimized Transonic Truss-Braced Wing (TTBW) concept using assisted electric propulsion during climb. First, we present results with and without electrification to show the overall system impacts of hybrid electric propulsion during climb. We also present two distinct battery models compatible with this optimization framework, and compare results using both. Next, we ran the same problem using a slightly different cell type, and demonstrate a considerable change in the result. Finally, we vary the cell energy density of the batteries, and provide an illustrative trend for the system-level impact for improving cell technology.

optimization↗

Classical Flutter Analysis of X-57 Aircraft Mod II, III, and IV Configurations

To support airworthiness assessment of the X-57 electric propulsion demonstrator aircraft, the National Aeronautics and Space Administration (NASA) completed classical flutter analysis of the Mod II, Mod III, and Mod IV aircraft configurations. Aeroelastic analysis was primarily performed using the ZAERO analysis code developed by Zona Technology. Supplemental computational fluid dynamics (CFD) flutter analysis for the Mod III and Mod IV configurations was performed using NASA’s FUN3D software to verify the aeroelastic behavior at deployed flaps conditions. The NASTRAN finite element models (FEMs) used for aeroelastic analysis were developed by NASA. The final Mod II FEM utilized shell elements for the wing and fuselage, and beam elements for the tail and wing control surfaces. For Mod III and Mod IV aircraft models, shell element wing FEMs developed throughout the Mod III/IV wing design cycle were utilized. The fuselage and empennage FEM components were shared between the Mod II, III, and IV aircraft models. ZAERO analyses predicted flutter boundaries were >60% beyond the X-57 flight-test envelope never exceed airspeed (V NE ) for all configurations. Aeroelastic Reynolds-averaged Navier–Stokes (RANS) CFD simulations using FUN3D of the Mod III/IV configurations predicted flap mode instability did not occur until at least speeds of 60% higher than the Mod III/IV maximum flap deployment airspeed (V FE ).

Keerti K. Bhamidipati↗

Novel Turbulence Model Tailored for Complex Rocket Motor Exhaust Jets

This paper describes a new RANS turbulence model intended to give improved accuracy in turbulence that is affected by compressibility, rotation and curvature, and variable density, making it well-suited to complex rocket motor exhaust jets. Four key elements of the model are: the use of simple-averaged (not mass-averaged) turbulence decomposition; a second-order Taylor series for the equilibrated turbulent fluxes; empirical modeling that accounts for the combined effects of aforementioned turbulence modifiers; and model tuning on large-eddy simulation data for mixing layers for a wide range of these turbulence modifiers. The model is implemented as an extension of the k-w SST model. Results include comparisons to detached eddy simulation results for hot and cold co-flowing exhaust jets, and to experimental data for supersonic retropropulsion on the CobraMRV Mars-lander concept vehicle. The model is named TMRC for its constituents elements, the Taylor series with Mach number, rotation and curvature effects.

turbulence modeling↗

Wall-Modeled Large Eddy Simulations of Transonic Buffet Over a Supercritical Airfoil

A series of scale-resolving simulations of flow over the ONERA OAT15A airfoil have been performed at an angle of attack of 3.5◦, just past the onset of buffet. The focus of this study is to document the sensitivity of the wall-modeled large eddy simulation (WMLES) methodology for curvilinear structured overset grids within the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework to mesh spacing, mesh distribution, and domain size. A secondary purpose of the study is to compare the results from WMLES to unsteady Reynolds-averaged Navier Stokes (URANS) simulations and hybrid RANS-LES (HRLES) within the same LAVA solver framework. The study provides a unique perspective regarding comparisons between different turbulence modeling approaches, time-integration methods, and computational performance since many of the same numerical routines are used for all three types of simulations. The results are compared with experiments and previous numerical studies of the same geometry and flow conditions.

TTT↗

Computationally Efficient Frequency Domain Method for Dynamic Control Derivative Estimation with Application to Transonic Truss-Braced Wing

This paper presents a computationally efficient method for dynamic control derivative estimation technique via control surface oscillation numerical experiments. Unsteady RANS CFD simulations in FUN3D are performed to simulate the control surface oscillations with a prescribed truncated square wave containing sufficient frequencies of interest. The truncated square wave oscillation offers the computational efficiency which reduces the computational cost by almost an order of magnitude compared to a sine wave oscillation. The nonlinear effect of large control surface oscillation amplitudes creates a spillover effect whereby the frequency response at the same input frequencies contains not only the linear aerodynamic response but also nonlinear aerodynamic response. A correction procedure is developed to remove the spillover effect from the linear aerodynamic response. A frequency-domain regression is performed to estimate the dynamic control derivatives after the correction. The results generally agree with the previous results obtained from the sine wave oscillation.

Stability and Control↗

Cali Urban Development II: Investigating the Impacts of Land Use Change on Urban Heat and Social Vulnerability in Cali, Colombia

The surface urban heat island (SUHI) effect is an environmental phenomenon resulting in cities with higher temperatures than rural areas due to increased pavement and decreased cooling from vegetation. The city of Santiago de Cali in Colombia faces SUHI challenges exacerbated by land use change. The Cali municipal government agency, Departamento Administrativo de Gestión del Medio Ambiente, and the community organization Fundacion Dinamizadores Ambientales partnered with NASA DEVELOP to evaluate communities in Cali most vulnerable to urban heat. This project illustrated the utility of using NASA Earth observations to evaluate the relationship between land use, temperature, and social factors in Cali, Colombia between 2013 and 2023. The team used Landsat 7 Enhanced Thematic Mapper Plus (ETM+), Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS), and Landsat 9 OLI-2/TIRS-2 to generate land surface temperature (LST) maps in Google Earth Engine through NASA DEVELOP’s Urban Heat Exposure Assessment Tempe 2.0 tool. Cloud cover limited the project feasibility, but it improved with Landsat 9 data. In ArcGIS Pro, the team found that LST was significantly higher in urban areas than in wetlands or forests. Using R Studio, the team ran a principal component analysis and found that health care and green space access were negatively correlated, and Afro-Colombian ethnicity was positively correlated with LST. With knowledge of the most impacted and vulnerable regions, the partner organizations can prioritize establishing healthcare facilities and green spaces in those areas to reduce the impacts of urban heat.

vegetation loss↗

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

In-Depth Modeling and Simulation Analysis of Artemis Missions Using the Impact Probabilistic Risk Assessment Tool

BACKGROUND The Artemis campaign is a Moon exploration program with a series of six planned missions, five of which will be crewed. These five crewed missions will contain a single mission segment (space flight), or multiple mission segments involving space flight (Orion), lunar landing (LTV) and/or space habitat (Gateway). Each crewed segment faces the risk of unique medical conditions, necessitating medical sets/kits tailored to those specificities. To support and enable a data-driven and evidence-based decision-making process through out a mission’s life cycle, a software tool called IMPACT was developed. Using probabilistic risk assessment (PRA) methodologies, IMPACT (Informing Mission Planning via Analysis of Complex Tradespaces) is a novel tool built for analyzing the possibility of encountering complex medical risks during space flight, and for identifying the medical resources and capabilities needed to treat those potential at-risk medical conditions. IMPACT achieves this by performing hundreds of thousands of Monte Carlo simulations of missions to build aggregate pictures of medical risk. During an extended simulation modeling phase, IMPACT generated analytical results for medical risks, and the medical resources and capabilities to address those risks, for every segment of every crewed Artemis mission. This presentation will highlight the reliability, consistency and validity of IMPACT’s computational modeling techniques and will showcase the library of analytical outcomes generated for the Artemis missions. OVERVIEW During the early stages of IMPACT’s design, architecture and technical requirements collection, “scenarios” (use cases) - achievement goals required for acceptance testing, were identified by stakeholders. IMPACT successfully completed the scenario testing requirements and undertook an extensive operational run phase utilizing a wide range of input combinations with a goal of delivering a cohesive, trustworthy, reliable, vast, and diverse body of evidence. The intent of these modeling runs was to validate consistency in output, ensure solidity of executable operations and to streamline processes by identifying areas requiring efficiency improvements. Using the many missions of Artemis, IMPACT ran variations of operational runs to assess the output for acceptable, as well as unusual characteristics. This rigorous long-term “shakedown” analysis was implemented to help build a collective body of evidence to aid in securing a high level of confidence, reliability, and validity in the output, whether from the applicational components of IMPACT, or the entirety of the operational process. ANTICIPATED ANALYSIS AND CONCLUSION This presentation will discuss the various categories of input criteria; the comparisons in the application of these input criteria to various Artemis missions; the preparation and collection of the body of evidence, and reliability of the computational modeling techniques. This paper serves as an initial analytical overview of IMPACT’s probabilistic risk assessment (PRA) medical risk outputs covering Artemis missions and is not intended to be deemed the official medical response for the Artemis campaign.

Crew Composition↗

System Modeling of a Lunar Molten Regolith Electrolysis Plant

Introduction: In-Situ Resource Utilization (ISRU) is the process of extracting local resources to produce commodities for propulsion, life support systems, and off-planet construction rather than transporting consumables from Earth. Molten Regolith Electrolysis (MRE) is a novel ISRU method of extracting oxygen gas and metal alloy from lunar regolith. The MRE process involves placing lunar regolith between two electrodes, through which current is passed, to melt the regolith and reduce the metal oxide constituents by direct electrolysis (e.g. FeO, SiO2, MgO, Al2O3) into oxygen gas and metal alloys. The oxygen is liquefied and used as propellant for landers, while the metals (e.g. Ferro-alloys) are further processed and used in structural building materials and parts manufacturing. A system model was developed that accounted for the major processes of an MRE plant (from excavation of raw materials to storage of products) to assess the feasibility of a lunar MRE plant. The System Engineering and Integration (SE&I) ISRU Modeling and Analysis (SIMA) team utilized its previously documented system sizing model, the Mission Analysis and Integration Tool (MAIT) [1] as framework of the system model. MAIT uses MATLAB/Simulink to integrate subsystem models into a complete system model of the MRE plant. Total mass, volume, and power requirements were computed for numerous iterations of a MRE plant. System Model: Figure 1: MRE Plant Block Diagram The regolith excavation model determines the mass and power needed to excavate sufficient regolith. The preheating auger initiates the regolith heating process before regolith enters the MRE re-actor to reduce the energy required to turn the solid into a molten liquid. The MRE reactor is modeled in COMSOL Multiphysics and based on the research by Dominguez, Sibille, and Schreiner [2, 3, 4]. This preliminary reactor model provides an accurate calculation of thermal equilibrium during electrochemical operation of the reactor system to assess the optimal mass and power required to process the inlet flow of regolith. The model also computes the outlet flowrates of oxygen and molten products. For this analysis, the primary components of the metal alloy considered were iron and silicon. The oxygen is then purified using an Yttrium Stabilized Zirconia (YSZ) electrode, followed by liquefaction using a 90K cryocooler to be stored as liquid oxygen in insulated cylindrical tanks. In future iterations of the system model, the molten metal tapped from the MRE reactor will undergo additional processing or refinement. However, downstream handling of metals is currently a technology gap that is missing a high TRL subsystem model. Therefore, for this analysis, the accumulated metal alloy stream terminates after leaving the MRE reactor. Study Goals: This analysis investigates multiple input variables to the system to determine the sensitivity of a (near) complete plant at full-scale. This preliminary investigation ran parametric sweeps on the MRE reactor geometry, electrical current supply, layers of multi-layer insulation (MLI) on the reactor, size of the electrodes in the oxygen purification model, and regolith composition (based on landing site location). Three production targets of oxygen (1,000, 10,000, and 50,000 kg/yr) were investigated for this analysis. The parametric sweeps conducted in this analysis provide valuable insight into the expected impact of the various model inputs on plant size. This information can be used to identify the most critical components of the plant and guide future decisions on allocating funding for research and development, providing subsystem developers with appropriate interfaces with downstream and upstream processes, and assessing the overall feasibility of MRE when compared to other ISRU plants. References: [1] Carlson, A. et al. (2024) ICES, ICES-2024-53. [2] Dominguez, D.A., and Sibille, L. (2011) AIAA, AIAA-2011-700. [3] Schreiner, S.S. (2015) MIT, Dissertation. [4] Schreiner, S.S. et al. (2016) ASR, 57(7), pp.1585-1603.

ISRU↗

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

High-Order Methods↗

Heated Supersonic Axisymmetric Jet Cases for the NASA Turbulence Modeling Resource

Three new supersonic jet cases were built for the NASA Turbulence Modeling Resource (TMR). These cases were taken from the 6th AIAA Propulsion Aerodynamics Workshop (PAW) nozzle test problem, which utilized data taken at NASA GRC consisting of non-intrusive measurements of jet plume velocities and temperatures, including both mean values and turbulence statistics. The three jet cases all used a Mach 1.63 axisymmetric nozzle: a heated jet for each at on-design conditions, a temperature-matched on-design condition at Mach 1.63, and an off-design heated condition for the Mach 1.63 nozzle. Computational fluid dynamics (CFD) solutions are obtained using three established Reynolds-averaged Navier-Stokes (RANS) codes: Wind-US, FUN3D, and VULCAN-CFD, all using the Menter Shear Stress Transport k- turbulence model with vorticity source term (SST-V). The NASA TMR requires at least two CFD codes yielding essentially identical results to certify that the approaches are verified, meaning that they solve the posed CFD equations as intended. The three codes generally provided very close agreement with each other for jet plume quantities, with the only exception being static temperature. This discrepancy was determined to be the result of whether or not turbulent kinetic energy was considered in the definition of total energy. Validation of the SST-V turbulence model for this case (its level of agreement/disagreement with experimental data) is also addressed.

turbulance↗

Investigating the Impacts of Land Use Change on Urban Heat and Vulnerability in Cali, Columbia

The urban heat island effect (UHI) is an environmental phenomenon where cities experience higher temperatures than rural areas due to increased pavement and decreased cooling from vegetation. Approximately 76% of people in Colombia live in urban areas, and the city of Santiago de Cali is facing UHI challenges exacerbated by land use change. Wetlands and forests formerly surrounded the city but were replaced by development and agriculture. The Colombian municipal government agency Departamento Administrativo de Gestión del Medio Ambiente and the community organization Fundacion Dinamizadores Ambientales partnered with NASA DEVELOP to evaluate communities in Cali most vulnerable to urban heat. This project illustrated the utility of using NASA Earth observations to evaluate the relationship between land use, temperature, and social factors in Cali, Colombia between 2013 and 2023. The team used Landsat 7 Enhanced Thematic Mapper Plus (ETM+), Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS), and Landsat 9 OLI-2/TIRS-2 to generate land surface temperature (LST), normalized difference vegetation index (NDVI), and albedo maps in Google Earth Engine. Heavy cloud cover limited the accuracy of the LST but incorporating up to three satellites for a median image reduced potential errors. Through further analysis in ArcGIS Pro, the team classified land use change using a deep learning model and found that LST was significantly higher in urban areas than in wetlands or forests. Using R studio, the team ran a principal component analysis to determine which social factors had the strongest correlation with LST. The team found that health care and green space access were negatively correlated, and Afro-Colombian ethnicity was positively correlated with LST. With awareness of the most impacted and vulnerable regions, the partner organizations can work to prioritize green space establishment in those areas to reduce the impacts of urban heat. Addressing the urban heat island effect will reduce environmental justice concerns within the city and improve overall health, air, and water quality for those who live there.

Brenna Bruffey↗