Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial Temporal Model”

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

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

At least 235 records · Page 13

Variabilities in Pm2.5 and Black Carbon Surface Concentrations Reproduced by Aerosol Optical Properties Estimated by In-Situ Data, Ground Based Remote Sensing and Modeling

Because of the increased temporal and spatial resolutions of the sensors onboard recently launched satellites, satellite-based surface aerosol concentration, which is usually estimated from the aerosol optical depth (AOD), is expected to become a strategic tool for air quality studies in the future. By further exploring the relationships of aerosol concentrations and their optical properties using ground observations, the accuracies of these products can be improved. Here, we analyzed collocated observations of surface mass concentrations of fine particulate matter (PM2.5) and black carbon (BC), as well as columnar aerosol optical properties from a sky radiometer and aerosol extinction profiles obtained by multi-axis differential optical absorption spectroscopy (MAX-DOAS), during the 2019–2020 period. We focused the analyses on a daily scale, emphasizing the role of the ultraviolet (UV) spectral region. Generally, the correlation between the AOD of the fine fraction (i.e., fAOD) and the PM2.5 surface concentration was moderately strong, regardless of considerations of boundary layer humidity and altitude. In contrast, the fAOD of the partial column below 1 km, which was obtained by combining sky radiometer and MAX-DOAS retrievals, better reproduced the variability of the PM2.5 and resulted in a linear relationship. In the same manner, we demonstrated that the absorption AOD of the fine fraction (fAAOD) of the partial column was related to the variability of the BC concentration. Analogous analyses based on aerosol products from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) confirmed these findings and highlighted the importance of the shape of the aerosol profile. Overall, our results indicated a remarkable consistency among the retrieved datasets, and between the datasets and MERRA-2 products. These results confirmed the well-known sensitivity to aerosol absorption in the UV spectral region; they also highlighted the efficacy of combined MAX-DOAS and sky radiometer observations.

Alessandro Damiani↗

Wake Vortex Interaction of Urban Air Mobility Aircraft

This paper presents a study of wake vortex interaction modeling for urban air mobility. Wake vortex modeling approaches for fixed-wing aircraft and rotorcraft are investigated. A wake age model is developed that accounts for the temporal and spatial dissipation of the wake-induced downwash. A wake vortex model is developed for rotorcraft that accounts for blade flapping motion and cyclic pitch control to reduce the lift asymmetry. This in turn results in an equal lift circulation strength on the advancing side and retreating side of the rotor. Two wake vortex interaction simulations are performed to illustrate the modeling approaches.

Wake Vortex↗

Neural network for processing both spatial and temporal data with time based back-propagation

Neural networks are computing systems modeled after the paradigm of the biological brain. For years, researchers using various forms of neural networks have attempted to model the brain's information processing and decision-making capabilities. Neural network algorithms have impressively demonstrated the capability of modeling spatial information. On the other hand, the application of parallel distributed models to the processing of temporal data has been severely restricted. The invention introduces a novel technique which adds the dimension of time to the well known back-propagation neural network algorithm. In the space-time neural network disclosed herein, the synaptic weights between two artificial neurons (processing elements) are replaced with an adaptable-adjustable filter. Instead of a single synaptic weight, the invention provides a plurality of weights representing not only association, but also temporal dependencies. In this case, the synaptic weights are the coefficients to the adaptable digital filters. Novelty is believed to lie in the disclosure of a processing element and a network of the processing elements which are capable of processing temporal as well as spacial data.

Villarreal, James A.↗

Discrete Aurora on the Nightside of Mars: Occurrence Location and Probability

This paper represents the first attempt to predict the occurrence location and probability of discrete electron aurora on the nightside of Mars. We run a 3-D time-dependent magnetohydrodynamic model to characterize the spatial and temporal dynamics of magnetic field and plasma distributions over the course of one planetary rotation. We perform eight simulation cases under solar minimum quiet-solar-wind conditions (four equinox/solstice seasons, each with two interplanetary magnetic field polarities) and in an actual interplanetary coronal mass ejection (ICME) case to assess quiet and space weather situations, respectively. The occurrence of detectable discrete aurora is subject to the combination of the probabilities that (a) the ionosphere is magnetically connected with high altitudes through open field lines and (b) precipitating energy fluxes of >30 eV electrons exceed 0.1 erg/cm2/s. Our results show that during quiet solar activity, discrete aurora occurs likely on small-scale patches embedded inside strong crustal magnetic field regions (with a magnitude greater than 50 nT at 150 km), and the overall chance across the globe is ∼0.77%. The higher probability over strong crustal field regions is attributed to the stronger magnetic field convergence. Modeling shows the occurrence probability dramatically increases during the ICME event, particularly by more than an order of magnitude in weak crustal field regions. Our model results reasonably agree with NASA Mars Atmosphere and Volatile EvolutioN and Mars Express observations. Our study suggests that nightside discrete electron aurora is not caused by the direct entry of magnetosheath plasma in a cusp-like process but due to the recycling of nightside magnetospheric electrons.

Xiaohua Fang↗

Using a Support Vector Machine and a Land Surface Model to Estimate Large-Scale Passive Microwave Temperatures over Snow-Covered Land in North America

A support vector machine (SVM), a machine learning technique developed from statistical learning theory, is employed for the purpose of estimating passive microwave (PMW) brightness temperatures over snow-covered land in North America as observed by the Advanced Microwave Scanning Radiometer (AMSR-E) satellite sensor. The capability of the trained SVM is compared relative to the artificial neural network (ANN) estimates originally presented in [14]. The results suggest the SVM outperforms the ANN at 10.65 GHz, 18.7 GHz, and 36.5 GHz for both vertically and horizontally-polarized PMW radiation. When compared against daily AMSR-E measurements not used during the training procedure and subsequently averaged across the North American domain over the 9-year study period, the root mean squared error in the SVM output is 8 K or less while the anomaly correlation coefficient is 0.7 or greater. When compared relative to the results from the ANN at any of the six frequency and polarization combinations tested, the root mean squared error was reduced by more than 18 percent while the anomaly correlation coefficient was increased by more than 52 percent. Further, the temporal and spatial variability in the modeled brightness temperatures via the SVM more closely agrees with that found in the original AMSR-E measurements. These findings suggest the SVM is a superior alternative to the ANN for eventual use as a measurement operator within a data assimilation framework.

Modeling↗

Optimization as a Tool for Consistency Maintenance in Multi-Resolution Simulation

The need for new approaches to the consistent simulation of related phenomena at multiple levels of resolution is great. While many fields of application would benefit from a complete and approachable solution to this problem, such solutions have proven extremely difficult. We present a multi-resolution simulation methodology that uses numerical optimization as a tool for maintaining external consistency between models of the same phenomena operating at different levels of temporal and/or spatial resolution. Our approach follows from previous work in the disparate fields of inverse modeling and spacetime constraint-based animation. As a case study, our methodology is applied to two environmental models of forest canopy processes that make overlapping predictions under unique sets of operating assumptions, and which execute at different temporal resolutions. Experimental results are presented and future directions are addressed.

TEMPORAL RESOLUTION↗

Stochastic simulations of temporally and spatially variable Fe cycling within a floodplain aquifer

This data package contains model input and output files for local and global sensitivity analysis (SA) of a reactive transport model simulating redox cycling within a floodplain aquifer. The reaction network is implemented in CrunchFlow and focuses on spatially and temporally heterogeneous Fe cycling, though it contains several other species and redox pathways (37 aqueous species and 7 minerals in total). The package contains data from 600 Monte Carlo simulations (for global SA) and 104 one-at-a-time simulations (for local SA). It includes all input files required to run the CrunchFlow model, as well as limited model output – .csv files containing time series of primary species concentrations at 21 grid cells within the domain. These time series files contain the following species: H+, CO32-, SO42-, Cl-, Ca2+, Mg2+, Na+, K+, Fe(II), Fe(III), HS-, CH3COO-, O2(aq), S(aq), and Br-.These data were used to demonstrate the application of a novel form of global SA, distance-based generalized sensitivity analysis (DGSA) to a reactive transport model. The 2D model simulates the export of reduced species from pockets of fine-grained, organic-rich sediments embedded in a coarse sand aquifer. Stochastic simulations jointly varied 17 key model input parameters across several orders of magnitude, including the spatial variability of each parameter. DGSA on the model results reveal that the amplitude and variability of exported Fe(II) are most sensitive to the interaction between sand permeability and individual reaction rates. By contrast, the propagation of reducing conditions downgradient from the fine-grained lenses depends most heavily on the rate of dissolved organic carbon and sulfate release from the fines.

54 ENVIRONMENTAL SCIENCES↗

Modeling and analysis of heat emissions from buildings to ambient air

Heat emissions from buildings is a significant source of anthropogenic heat influencing the urban microclimate; however, they are usually oversimplified in urban climate and microclimate modeling. This study developed a bottom-up physics-based approach to calculate heat emissions from buildings to the ambient air and implemented the approach in EnergyPlus. In this work, a simple result verification was conducted by comparing the EnergyPlus simulated results against the spreadsheet calculations. Simulations covering 16 commercial building types, four climates, and two energy efficiency levels were conducted to understand and evaluate the building heat emissions and their temporal patterns as well as three major components: (1) building envelope (convective heat transfer to ambient air), (2) zones (air exfiltration and exhaust air), and (3) HVAC systems (relief air and heat rejection from condensers or cooling towers). The main findings are: (1) heat emissions are usually higher than the site energy use (about 2.5 times), and their dynamics should be considered; (2) building characteristics and their energy systems lead to differences in heat emission contributions from the three components, and their dynamics, for example, in the warehouse models, the envelope component accounts for 90.4%, while it is 12.7% for the large office models; (3) for most building typologies, the climate has a strong impact on heat emissions, for example, buildings with dominant heat emissions from the zone exhaust air and/or the HVAC reject heat, a general decrease in heat emissions in hotter climates is observed, while envelope-dominated buildings show the opposite; and (4) building technologies that reduce energy use in buildings may perform differently in reducing heat emissions. The developed heat emissions calculation method can be adopted in EnergyPlus and most other building energy modeling programs. It can provide dynamic building heat emissions as an input to urban climate computational fluid dynamics (CFD) models at a higher spatial and temporal resolution than is currently available, to improve the simulation accuracy of the urban microclimate and capture the urban heat island effect and urban overheating.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Flow Characterization of the NASA Langley Unitary Plan Wind Tunnel, Test Section 2: Computational Results

This is an abstract for an invited paper at the AIAA Aviation Conference, June 2021. The work described here is part of an effort of coordinated experiments in the Unitary Plan Wind Tunnel (UPWT) facility at the NASA Langley Research Center (LaRC) and matching CFD simulations. The primary goal of the work is to assess the productivity and true predictive accuracy of CFD, absent any guidance from experiment, in the high supersonic speed range as compared to experiments performed at the NASA LaRC’s UPWT facility. This report concerns CFD simulation of the primary flow-path in the empty wind tunnel, including the settling chamber, nozzle, test section, and some of the tunnel downstream of the test section. The Mach number in the test section ranges from M~2.4 to M~4.6, and the required area ratio variation is achieved by translation of a nozzle block which constricts the area of a loosely S-shaped throat. Flow past protuberances in the settling chamber and into this S-bend throat are predicted by CFD to generate streamwise vorticity that may, or may not, persist through the throat and into the test section as coherent vortices. Some flow conditions are notably unsteady at frequencies well below those of turbulence, due to unsteady separated flow ahead of the nozzle block. The bulk flow moves at velocities ranging from 'walking speed' in the settling chamber to M~4.6 in the test section. Heat transfer to the settling chamber walls and buoyancy are significant at high Mach number. Subtle variations in surface curvature in the nozzle generate Mach waves that propagate into the test section. CFD of the empty tunnel serves two purposes. Firstly, the full-tunnel simulations are used to provide upstream boundary conditions for CFD of vehicle aerodynamics which are generally performed in a domain confined to the wind tunnel test section; these companion studies are addressed in other papers. Secondly, it is a challenging test for CFD to resolve all of the empty tunnel flow phenomena relevant to flow in the test section. It requires a more complete definition of geometry than was originally anticipated. In addition, it requires good spatial and temporal accuracy, and turbulence modeling that performs well on specific phenomena such as corner flows. The boundary conditions and solution algorithms must perform well from incompressible to almost hypersonic speeds. The final state of the pre-test CFD was a product of an iterative self-improvement process. The initial simulations of the empty tunnel were deficient in many respects, but hints to those deficiencies were recognized in the solutions, and remedies were implemented. Possible further improvements will be studied in the post-test phase when comparisons with experimental data are possible. The CFD was performed by five separate collaborative teams using four different flow solvers: FUN3D, Overflow, Star-CCM+ and USM3D. The level of effort of these teams varied significantly, but each made important contributions to the goals of the work. A concerted effort to use uncertainty quantification methods (UQ) in CFD is also a goal of this work. To this end, variations in CFD results due to grid refinement, turbulence modeling, and boundary conditions have been characterized. Code-to-code variation is another means of assessing CFD uncertainty. All CFD solvers predict similarity among the primary flow features; these include the variation of Mach number due to changes in Reynolds number, and the bulk flow angularity due to tunnel-wall curvature. All CFD solvers also predict similar trends in secondary flows, such as the downwash in the side-wall boundary layers. Three of the high-spatial resolution simulations give similar predictions of a complex secondary flow phenomena, streamwise vortices generated in the S-bend throat that persist into the side-wall boundary layers of the test section. Two of the highest-resolution simulations, run with the same turbulence model in different CFD solvers, gave encouragingly similar predictions of a complex tertiary flow phenomenon, small transient "sprites" of upwash flows, resulting from vortices that presumably originate in the separated flow near the leading edge of the nozzle block. While the CFD was done in a "blind pre-test" mode, requests from the experimental team for CFD results pertaining to unsteadiness and total temperature variations in the test section resulted in CFD runs that included a cooled wall in settling chamber. This then led to a change in the standard practice for running the Overflow results, which would not have occurred without this "release" of this experimental information. CFD was also used to guide some measurements. The paper will focus on establishing the consensus among CFD results and understanding differences among those results. Initial findings from the efforts to characterize CFD uncertainty have been done and will be included in the paper.

Robert Edward Childs↗

NASA Tech Briefs, December 2008

Topics covered include: Crew Activity Analyzer; Distributing Data to Hand-Held Devices in a Wireless Network; Reducing Surface Clutter in Cloud Profiling Radar Data; MODIS Atmospheric Data Handler; Multibeam Altimeter Navigation Update Using Faceted Shape Model; Spaceborne Hybrid-FPGA System for Processing FTIR Data; FPGA Coprocessor for Accelerated Classification of Images; SiC JFET Transistor Circuit Model for Extreme Temperature Range; TDR Using Autocorrelation and Varying-Duration Pulses; Update on Development of SiC Multi-Chip Power Modules; Radio Ranging System for Guidance of Approaching Spacecraft; Electromagnetically Clean Solar Arrays; Improved Short-Circuit Protection for Power Cells in Series; Electromagnetically Clean Solar Arrays; Logic Gates Made of N-Channel JFETs and Epitaxial Resistors; Improved Short-Circuit Protection for Power Cells in Series; Communication Limits Due to Photon-Detector Jitter; System for Removing Pollutants from Incinerator Exhaust; Sealing and External Sterilization of a Sample Container; Converting EOS Data from HDF-EOS to netCDF; HDF-EOS 2 and HDF-EOS 5 Compatibility Library; HDF-EOS Web Server; HDF-EOS 5 Validator; XML DTD and Schemas for HDF-EOS; Converting from XML to HDF-EOS; Simulating Attitudes and Trajectories of Multiple Spacecraft; Specialized Color Function for Display of Signed Data; Delivering Alert Messages to Members of a Work Force; Delivering Images for Mars Rover Science Planning; Oxide Fiber Cathode Materials for Rechargeable Lithium Cells; Electrocatalytic Reduction of Carbon Dioxide to Methane; Heterogeneous Superconducting Low-Noise Sensing Coils; Progress toward Making Epoxy/Carbon-Nanotube Composites; Predicting Properties of Unidirectional-Nanofiber Composites; Deployable Crew Quarters; Nonventing, Regenerable, Lightweight Heat Absorber; Miniature High-Force, Long-Stroke SMA Linear Actuators; "Bootstrap" Configuration for Multistage Pulse-Tube Coolers; Reducing Liquid Loss during Ullage Venting in Microgravity; Ka-Band Transponder for Deep-Space Radio Science; Replication of Space-Shuttle Computers in FPGAs and ASICs; Demisable Reaction-Wheel Assembly; Spatial and Temporal Low-Dimensional Models for Fluid Flow; Advanced Land Imager Assessment System; Range Imaging without Moving Parts.

Source record↗

Expanding the Time Dimension of Hyperspectral Infrared Sounding Observations: Designing the NUCAPS-Forecast System

Low Earth orbit (LEO) hyper-spectral infrared (IR) sounders have significant yet untapped potential for characterizing thermodynamic environments of convective initiation and ongoing convection. While LEO soundings are of enormous value to weather forecasters, they cannot provide the temporal resolution needed to resolve the rapidly evolving thermodynamics of the convective environment. A novel methodology was developed to extend snapshots of LEO soundings forward in time up to six hours into the future and has been made available within the NOAA National Weather Service Advanced Weather Interactive Processing System (AWIPS) for testing in a quasi-operational setting. This methodology is based on parcel forward-trajectory calculations from the satellite observing time to recreate future soundings of temperature and moisture at regularly gridded intervals in space and time. The soundings are based on NOAA-Unique Combined Atmospheric Processing System (NUCAPS) retrievals from the Suomi NPP and NOAA-20 satellite platforms to allow for rapid near-real time product processing. Science code was transitioned to operational production in 2019 and the product, NUCAPS-Forecast, was tested at the NOAA Hazardous Weather Testbed Spring Experiment in 2019 and 2021. Forecaster feedback has led to improving the product and processing system to address product limitations and improve its applicability for pre-convective forecasting. This presentation highlights the design of the NUCAPS-Forecast system. Solutions were tested and implemented to remove excess moisture and conjoin forward-trajected parcels into vertical profiles for calculating grids of stability indices. Additionally, use cases that explore the tendencies of derived convective available potential energy (CAPE) and convective inhibition (CIN) were investigated. When evaluated against gridded, hourly accumulated rainfall obtained from the Multi-Radar Multi-Sensor (MRMS) observational data set areas with elevated CAPE (reduced CIN) are shown to be associated with areas of precipitation and increase confidence the LEO soundings advected forward in time physically represent atmospheric processes. These results suggest that adiabatic parcel advection of LEO satellite sounding snapshots forward in time exhibits skill in discriminating locations and times that are more likely to lead to convective initiation. Leveraging trajectory modeling to increase the spatial and temporal resolution of LEO IR sounding observations today demonstrates the value of these observations from the future constellation of JPSS satellites and GeoXO mission.

weather↗

Modeling Forest Biomass and Growth: Coupling Long-Term Inventory and Lidar Data

Combining spatially-explicit long-term forest inventory and remotely sensed information from Light Detection and Ranging (LiDAR) datasets through statistical models can be a powerful tool for predicting and mapping above-ground biomass (AGB) at a range of geographic scales. We present and examine a novel modeling approach to improve prediction of AGB and estimate AGB growth using LiDAR data. The proposed model accommodates temporal misalignment between field measurements and remotely sensed data-a problem pervasive in such settings-by including multiple time-indexed measurements at plot locations to estimate AGB growth. We pursue a Bayesian modeling framework that allows for appropriately complex parameter associations and uncertainty propagation through to prediction. Specifically, we identify a space-varying coefficients model to predict and map AGB and its associated growth simultaneously. The proposed model is assessed using LiDAR data acquired from NASA Goddard's LiDAR, Hyper-spectral & Thermal imager and field inventory data from the Penobscot Experimental Forest in Bradley, Maine. The proposed model outperformed the time-invariant counterpart models in predictive performance as indicated by a substantial reduction in root mean squared error. The proposed model adequately accounts for temporal misalignment through the estimation of forest AGB growth and accommodates residual spatial dependence. Results from this analysis suggest that future AGB models informed using remotely sensed data, such as LiDAR, may be improved by adapting traditional modeling frameworks to account for temporal misalignment and spatial dependence using random effects.

Babcock, Chad↗

Evaluating three evapotranspiration estimates from model of different complexity over China using the ILAMB benchmarking system

The land surface models range in complexity of terrestrial evapotranspiration, yet it is unknown how model complexity translates to accuracy of modeled evapotranspiration estimates. Here, we use the International Land Model Benchmarking system to assess ET estimates from three models of varying complexity driven by the same forcing datasets: an earth system model, a terrestrial biosphere model, and a stand-alone ET model. The performance assessment includes both temporal and spatial evaluation, and different plant functional types across China. Our results indicate that the most complex model, an earth system model, performed best against the benchmarking datasets and metrics. Terrestrial biosphere model performed best in simulating inter-annual variability of ET, while earth system model performed best in simulating the seasonal cycle. The more complex models (earth system model and terrestrial biosphere model) perform better in forest, shrub and crop ecosystems, while the simpler model (stand-alone ET model) perform better in grass ecosystems. Our study demonstrates the impact of model complexity on ET estimates and highlights directions for future ET model improvements.

54 ENVIRONMENTAL SCIENCES↗

A New Capability of E4D For 3D Parallel Joint Inversion of DC Resistivity And Traveltime Data on Unstructured Mesh

A major challenge in interpreting geophysical data is how to derive consistent three-dimensional (3D) earth models of different physical properties from spatially and temporally limited measurements. Joint inversion with cross-gradient constraints is an approach to find such models by imposing structural similarities between different physical parameters. We have developed a parallel distributed-memory joint inversion code for direct-current (DC) resistivity and traveltime data using the cross-gradient constraint on unstructured mesh. The code utilizes existing E4D framework for parallel forward simulation, distributed storage and computation of the Jacobian matrix of forward operator, and parallel execution of matrix-vector multiplication during inversion. Besides, the joint inversion is solved by nonlinear conjugate gradient algorithm parallelized for DC resistivity and traveltime data. The joint inversion capability of E4D was tested using synthetic data from cross-borehole DC resistivity and traveltime data. The results indicate that the shape and size of the anomalies from the joint inversion are more reliable than those from separate inversions.

58 GEOSCIENCES↗

Modelling Airborne Transmission and Ventilation Impacts of a COVID-19 Outbreak in a Restaurant in Guangzhou, China

Computational fluid dynamics (CFD) modelling was performed to simulate spatial and temporal airborne pathogen concentrations during an observed COVID-19 outbreak in a restaurant in Guangzhou, China. The reported seating configuration, overlap durations, room ventilation, layout, and dimensions were modelled in the CFD simulations to determine relative exposures and probabilities of infection. Results showed that the trends in the simulated probabilities of infection were consistent with the observed rates of infection at each of the tables surrounding the index patient. Alternative configurations that investigated different boundary conditions and ventilation conditions were also simulated. Increasing the fresh-air percentage to 10%, 50%, and 100% of the supply air reduced the accumulated pathogen mass in the room by an average of ~30%, ~70%, and ~80%, respectively, over 73 min. Overall, the probability of infection was reduced by ~10%, 40%, and 50%, respectively.

60 APPLIED LIFE SCIENCES↗

Coincident Occurrences of Tropical Individual Cirrus Clouds and Deep Convective Systems Derived from TRMM Observations

Measurements of cloud properties and atmospheric radiation taken between January and August 1998 by the Tropical Rainfall Measuring Mission (TRMM) satellite were used to investigate the effect of spatial and temporal scales on the coincident occurrences of tropical individual cirrus clouds (ICCs) and deep convective systems (DCSs). It is found that there is little or even negative correlation between instantaneous occurrences of ICC and DCS in small areas, in which both types of clouds cannot grow and expand simultaneously. When spatial and temporal domains are increased, ICCs become more dependent on DCSs due to the origination of many ICCs from DCSs and moisture supply from the DCS in the upper troposphere for the ICCs to grow, resulting in significant positive correlation between the two types of tropical high clouds in large spatial and long temporal scales. This result may suggest that the decrease of tropical high clouds with SST from model simulations is likely caused by restricted spatial domains and limited temporal periods. Finally, the radiative feedback due to the change in tropical high cloud area coverage with sea surface temperature appears small and about -0.14 W/sq m per degree Kelvin.

Lin, Bing↗

High fidelity simulations of contaminant dispersion in an urban environment with comparison to magnetic resonance imaging measurements

The dispersion of a contaminant in an urban environment has the potential to impact a large population of people. In this work, a complex urban canopy flow based on the Oklahoma City downtown business district circa 2003 is studied using Magnetic Resonance Imaging (MRI) and high-fidelity Large Eddy Simulations (LES). MRI is a novel experimental technique that can provide high-resolution measurements in four dimensions (three spatial and temporal) for lab scale models. The experiments and simulations use the same geometry and boundary conditions providing a one-to-one comparison of the two methods. Results are presented on the time-averaged velocity and concentration fields, the temporal dynamics of the concentration plumes for a transient release, and a novel Cloud Identification Algorithm that can separate plumes produced by periodic contaminant releases used for ensemble averaging over many releases. The MRI and LES datasets both include millions of measurement voxels and the comparisons highlight the complex 3D nature of the flow including strong vertical velocities in spanwise street canyons and flow acceleration in streamwise street canyons. The concentration fields are qualitatively similar albeit the LES shows larger dispersion. A quantitative analysis with performance measures compares the datasets pointwise and demonstrates that the two 3D datasets are similar with respect to many measures including a fractional bias of 0.02 (ideal=0.0), correlation coefficient of 0.87 (ideal = 1.0), and the fraction points within a factor of 2 is 0.98 (ideal = 1.0). Plume analysis compares the arrival and residence time of contaminant and is found to vary significantly with location within the urban environment with arrival times between 0 and 1.25 and differences within the contaminant cloud less than 10% at most locations.

54 ENVIRONMENTAL SCIENCES↗