Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “point clouds”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

ARM shortwave spectrometers to study the clear-cloud transition zone and mixing processes

The proposed research is a collaborative effort between NASA/Goddard Space Flight Center and the Hebrew University of Jerusalem. While the NASA team focused on analyzing ground-based hyper-spectral radiance observations to understand cloud edge properties and their connection to mixing processes, the Hebrew University team tackled the problem through cloud modeling activities. By approximating the shortwave spectra in the cloud-clear transition zone as a linear combination of purely clear and purely cloudy spectra we can characterize the variations of cloud optical thickness and cloud droplet effective radius in the transition zone. When applying this method to the measurements of a ground-based shortwave spectroradiometer at the ARM’s SGP site, representing continental conditions, and MAGIC field campaign between Log Angeles, California and Honolulu, Hawaii, representing maritime scenarios, we found that cloud optical depth consistently decreases in both cases, but droplet size decreases much more substantially for the continental regime, suggesting different mixing processes for the continental and maritime conditions. The investigation and measurements of radiation clouds were coupled with a unique cloud modeling. A novel spectral bin microphysics was developed and implemented to the System of Atmospheric Modeling (SAM). In order to resolve cloud transition zones with high spatial gradients of microphysical variables a unique high resolution (10 m) was used in simulations. The model calculates droplet size distributions in each grid point. The model output was transferred to the NASA/GSFC team for utilization in radiative calculations and testing of both radiative algorithm and model representation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Dissipation of Mixed-Phase Arctic Clouds and Its Relationship to Aerosol Properties - Final Technical Report

The overarching objective of this project was to further our understanding of low-level Arctic cloud behavior in ultra-low aerosol concentration environments. To this end, we used a combination of modeling and analysis of DOE ARM datasets to understand the frequency of aerosol-limited cloud dissipation events, to understand the microphysical processes that lead to aerosol-limited dissipation, and to understand how above-cloud aerosol particles can work to maintain low-level clouds that exist in boundary layers with ultra-low aerosol concentrations. We found that aerosol-limited dissipation is unlikely to occur at the DOE long-term observation facilities at North Slope Alaska and Oliktok Point but possibly does occur at higher latitudes in the Arctic. At these higher latitudes, aerosol concentrations are frequently observed to be higher, often substantially higher, than in the boundary layer. These reservoirs can be sufficient to maintain low-level clouds. When they are insufficient, clouds appear to dissipate due to a depletion of liquid water through precipitation that is reinforced by a shutdown of the dynamical processes that support these clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Kinematical Analysis of Substructure in the Southern Periphery of the Large Magellanic Cloud

We report the first 3D kinematical measurements of 88 stars in the direction of several recently discovered substructures in the southern periphery of the Large Magellanic Cloud (LMC) using a combination of Gaia proper motions and radial velocities from the APOGEE-2 survey. More specifically, we explore stars in assorted APOGEE-2 pointings in a region of the LMC periphery where various overdensities of stars have previously been identified in maps of stars from Gaia and DECam. By using a model of the LMC disk rotation, we find that a sizable fraction of the APOGEE-2 stars have extreme space velocities that are distinct from, and not a simple extension of, the LMC disk. Using N-body hydrodynamical simulations of the past dynamical evolution and interaction of the LMC and Small Magellanic Cloud (SMC), we explore whether the extreme-velocity stars may be accounted for as tidal debris created in the course of that interaction. We conclude that the combination of LMC and SMC debris produced from their interaction is a promising explanation, although we cannot rule out other possible origins, and that these new data should be used to constrain future simulations of the LMC–SMC interaction. We also conclude that many of the stars in the southern periphery of the LMC lie outside of the LMC plane by several kiloparsecs. Given that the metallicity of these stars suggests that they are likely of Magellanic origin, our results suggest that a wider exploration of the past interaction history of the Magellanic Clouds is needed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A framework to measure the properties of intergalactic metal systems with two-point flux statistics

ABSTRACT The abundance, temperature, and clustering of metals in the intergalactic medium are important parameters for understanding their cosmic evolution and quantifying their impact on cosmological analysis with the Ly α forest. The properties of these systems are typically measured from individual quasar spectra redward of the quasar’s Ly α emission line, yet that approach may provide biased results due to selection effects. We present an alternative approach to measure these properties in an unbiased manner with the two-point statistics commonly employed to quantify large-scale structure. Our model treats the observed flux of a large sample of quasar spectra as a continuous field and describes the one-dimensional, two-point statistics of this field with three parameters per ion: the abundance (column density distribution), temperature (Doppler parameter), and clustering (cloud–cloud correlation function). We demonstrate this approach on multiple ions (e.g. ${\rm C\, \small {\rm IV}}$ , ${\rm Si\, \small {\rm IV}}$ , and ${\rm Mg\, \small {\rm II}}$ ) with early data from the Dark Energy Spectroscopic Instrument (DESI) and high-resolution spectra from the literature. Our initial results show some evidence that the ${\rm C\, \small {\rm IV}}$ abundance is higher than previous measurements and evidence for abundance evolution over time. The first full year of DESI observations will have over an order of magnitude more quasar spectra than this study. In a future paper, we will use those data to measure the growth of clustering and its impact on the Ly α forest, as well as test other DESI analysis infrastructure such as the pipeline noise estimates and the resolution matrix.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Occurrence of seeding multi-layer clouds in the Arctic from ground-based observations

Studies of Arctic clouds often focus on low-level single-layer clouds (SLCs). Here, we use combined observations of soundings and cloud radar during the MOSAiC, ACSE, and AO2018 research cruises as well as from long-term observations at Ny-Ålesund, Svalbard and Utqiagvik, Alaska to investigate the occurrence of SLCs and multi-layer clouds (MLCs) in the Arctic and to assess the rate of ice-crystal seeding in cold MLCs. MOSAiC observations show cloudy conditions in between 70 % and 90 % of sounding-radar cases. SLCs show occurrence rates of 30 % to 40 % with the highest value of 45 % during October. Cold MLCs are most abundant from November to June (40 % to 55 % of cases). Seeding occurs in about half to two thirds of the identified cold MLCs during MOSAiC for which the sub-saturated layer extends between 100 and 1000 m. The seeding rate increases by as much as 20 percentage points as the assumed size of the falling ice crystals is increased from 100 to 400 µm. The observations reveal a stable rate of cloud-free conditions of around 20 % over the covered latitude range. Cloud occurrence during MOSAiC and at Ny-Ålesund in July, when the geographical distance between observations was minimal, shows reasonable agreement. Comparisons of MOSAiC and other research cruises to the central Arctic also indicate consistent occurrence rates of different cloud types despite the likely effect of year-to-year variability. The comparison of data from ship campaigns and land sites suggests that the latter are not necessarily a good indicator of cloud occurrence in the high Arctic.

Achtert, Peggy [Leipzig Univ. (Germany); Meteorolo↗
DOE OSTI2026

Grid Operator Analytics and Assessment Tools for Inverter- Based Resources Dominated Grid (GOAAT-IBR) Project Update

This presentation provides an update on the OPTIMA GOAAT project, with emphasis on the cloud-native data platform developed in-house to ingest, manage, and operationalize high-resolution power system data. Since our last NASPI presentation, accessible via OSTI ID #2671437, the project team advanced the design and deployment of a scalable architecture capable of handling both synchronized and non-synchronized streams, including PMU, point-on-wave (POW), COMTRADE, and SCADA data. These materials review the project status, recent progress, and key lessons learned. The core of the presentation examines the architecture and engineering of our cloud-native ingestion and data management platform. We then explain how pipelines were designed to collect, normalize, time-align, store, and serve heterogeneous data at scale. We will discuss design choices such as data models, streaming versus batch ingestion, storage tiers, and interoperability with analytics applications. Practical experiences with cloud-native technologies were shared during the event, including benefits, limitations, and integration challenges in a utility environment, along with methods used to improve performance, reduce latency, and optimize resource usage. The presentation also showcases user interface designs and visualization tools that convert raw measurements and analytics results into intuitive, actionable insights for operators and engineers. During the presentation examples were provided demonstrating how visualization, event views, and summarized analytics enhance situational awareness and support operational decision-making. These use cases illustrate how a well-designed data infrastructure can bridge the gap between high-volume measurements and practical grid operations.

Aminifar, Farrokh↗
DOE OSTI2021

Total X-Ray Emission from the LMC Observed with HaloSat

We constrain the global parameters of the Large Magellanic Cloud (LMC) by modeling the first soft X-ray (0.4–6.8 keV) observations of the entire LMC produced in a single pointing with moderate spectral resolution by the HaloSat CubeSat. These data are best fit with the sum of two thermal plasma components in collisional ionization equilibrium and a power law. We find cool (0.210 ± 0.014 keV) and hot (0.89 ± 0.14 keV) components. The total X-ray luminosity of the LMC is (1.08 ± 0.14) × 10{sup 39} erg s{sup −1}. X-ray binaries make up a large fraction of the emission with a luminosity of (6.0 ± 0.8) × 10{sup 38} erg s{sup −1}, followed by cool gas from superbubbles, supernovae, and diffuse emission with a luminosity of (3.0 ± 0.3) × 10{sup 38} erg s{sup −1}. The hot gas from star formation contributes the smallest fraction, with a luminosity of (1.9 ± 0.5) × 10{sup 38} erg s{sup −1}. We estimate the total volume of the cool gas to be between (0.2–1.2) × 10{sup 10} pc{sup 3} and the hot gas to be between (1.0–5) × 10{sup 7} pc{sup 3} for filling factors of f = 1 and 0.2. These volumes result in a total thermal energy for the cool gas between (1.4–3) × 10{sup 54} erg for electron densities of 0.017–0.04 cm{sup −3}, and a thermal energy for the hot gas between (1.7–4) × 10{sup 53} erg for electron densities of 0.12–0.3 cm{sup −3}. This yields cooling timescales for the cool and hot gas of (1.5–3) × 10{sup 8} yr and (1.8–6) × 10{sup 7} yr, respectively.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Crossroads Able "Gilda"

At 5:55 am local time on June 30, 1946, a B-29 bomber took to the air carrying a Fat Man bomb, christened Gilda. Named after a Rita Haworth movie character, Gilda was released over a target array of ninety-three ships, including the Japanese battleship Nagato, the German cruiser Prinz Eugen, and the United States aircraft carrier Independence, all of which carried varying amounts of munitions, fuel, and animals. Although Gilda missed its aiming point, the explosion sank five ships and severely damaged eight others. A large number of pigs, goats, and rats were exposed to lethal doses of ionizing radiation. Visually, Gilda was not that impressive, particularly to those expecting an apocalyptic event. As the New York Times reporter William Laurence recalled, “To some of the newspaper men aboard, keyed up to the point of expecting the observer ship to be blown out of the water, the spectacle, obscured somewhat by an intervening white cloud, was a disappointment. To me, who could distinguish between the natural cloud and the atomic cloud, the sight was awesome and spine chilling.”

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Characteristics and sources of fluorescent aerosols in the central Arctic Ocean

The Arctic is sensitive to cloud radiative forcing. Due to the limited number of aerosols present throughout much of the year, cloud formation is susceptible to the presence of cloud condensation nuclei and ice nucleating particles (INPs). Primary biological aerosol particles (PBAP) contribute to INPs and can impact cloud phase, lifetime, and radiative properties. We present yearlong observations of hyperfluorescent aerosols (HFA), tracers for PBAP, conducted with a Wideband Integrated Bioaerosol Sensor, New Electronics Option during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition (October 2019–September 2020) in the central Arctic. We investigate the influence of potential anthropogenic and natural sources on the characteristics of the HFA and relate our measurements to INP observations during MOSAiC. Anthropogenic sources influenced HFA during the Arctic haze period. But surprisingly, we also found sporadic “bursts” of HFA with the characteristics of PBAP during this time, albeit with unclear origin. The characteristics of HFA between May and August 2020 and in October 2019 indicate a strong contribution of PBAP to HFA. Notably from May to August, PBAP coincided with the presence of INPs nucleating at elevated temperatures, that is, >–9°C, suggesting that HFA contributed to the “warm INP” concentration. The air mass residence time and area between May and August and in October were dominated by the open ocean and sea ice, pointing toward PBAP sources from within the Arctic Ocean. As the central Arctic changes drastically due to climate warming with expected implications on aerosol–cloud interactions, we recommend targeted observations of PBAP that reveal their nature (e.g., bacteria, diatoms, fungal spores) in the atmosphere and in relevant surface sources, such as the sea ice, snow on sea ice, melt ponds, leads, and open water, to gain further insights into the relevant source processes and how they might change in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Analysis of the EMCCD point-source response using x-rays

Electron Multiplying Charge Coupled Devices, EMCCD are used as x-ray detectors. The NSLS-II Soft Inelastic x-ray Scattering (SIX) beam line has two EMCCDs for x-ray detection in the spectrometer arm. The spectrometer with high resolving power disperses x-rays vertically. The x-ray vertical position on the sensor plane is related to its energy. This allows for very accurate x-ray energy measurements through x-ray coordinates. X-rays interact with silicon and create a number of electron-hole pairs proportional to the x-ray energy. Electrons drift and diffuse toward pixel gates and are collected there. The diffused electrons form a charge cloud distributed over several neighboring pixels. This charge sharing enables coordinate measurements with accuracy better than the pixel pitch. The charge distribution shape has to be taken into account to achieve ultimate accuracy in coordinate measurements. In this paper, we present a method of the charge distribution shape analysis and demonstrate its applications. The drift and diffusion of electrons ftom the point of generation to pixel gates results in the bell-shaped electron cloud usually approximated by Gaussian shape. The number of electrons collected under a pixel is proportional to the shape function integral. These electron packets get transferred to the sense node of the output amplifier. The transfer process could introduce distortions to the original charge distribution. For example, during transfers, electrons in the packet could be exposed to traps if they are present in the sensor. The trapping and later the release processes distort the apparent shape of the charge distribution. Therefore, deviations of the charge distribution shape ftom the originally symmetrical form can indicate the presence of trap centers in the sensor and can be used for sensor diagnostics.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Responses of the AC/DC Global Electric Circuit to Volcanic Electrical Activity in the Hunga Tonga-Hunga Ha'apai Eruption on 15 January 2022

Responses of the AC and DC global electric circuits (GECs) to the large eruption of the Hunga Tonga-Hunga Ha'apai (HT-HH) volcano on 15 January 2022 are discussed. The AC-related investigation is based on Schumann resonance (SR) measurements from six stations on four continents. The DC-related investigation utilizes atmospheric electric field (potential gradient, PG) measurements from six recording stations in Europe and the USA. According to data from the GLD360 and WWLLN lightning detection networks, the peak lightning stroke rate, 83/s, was dominated by negative polarity lightning, but the distributions of positive and negative lightning discharges in latitude and longitude around the volcano differed. A global intensification of SR is apparent in connection with the enhanced lightning activity caused by the eruption. SR data-based results confirm that the lightning activity in the eruption dominated the naturally occurring global activity for a period of about 1 hr. The highly localized increase in lightning activity over HT-HH was a unique point source of SR excitation. PG measurements suggest that impulse-like charging of the DC GEC, by ~15%, via negative cloud-to-ground lightning strokes took place twice during the eruption. A time constant of 7 or 8 min has been inferred for near-surface PG changes due to these enhancements. This could be the first direct measurement of the time constant of the GEC near the Earth's surface, as well as the first observation of the direct charging of the DC GEC by a unique atmospheric electrified source.

58 GEOSCIENCES↗
DOE OSTI2021

Polarimetric Radar Variables in Snowfall at Ka- and W-Band Frequency Bands: A Comparative Analysis

Dual-frequency millimeter-wavelength radar observations in snowfall are analyzed in order to evaluate differences in conventional polarimetric radar variables such as differential reflectivity (ZDR) specific differential phase shift (K DP ) and linear depolarization ratio (LDR) at traditional cloud radar frequencies at Ka and W bands (~35 and ~94 GHz, correspondingly). Low radar beam elevation (~5°) measurements were performed at Oliktok Point, Alaska, with a scanning fully polarimetric radar operating in the horizontal–vertical polarization basis. This radar has the same gate spacing and very close beam widths at both frequencies, which largely alleviates uncertainties associated with spatial and temporal data matching. It is shown that observed Ka- and W-band Z DR differences are, on average, less than about 0.5 dB and do not have a pronounced trend as a function of snowfall reflectivity. The observed Z DR differences agree well with modeling results obtained using integration over nonspherical ice particle size distributions. For higher signal-to-noise ratios, KDP data derived from differential phase measurements are approximately scaled as reciprocals of corresponding radar frequencies indicating that the influence of non-Rayleigh scattering effects on this variable is rather limited. This result is also in satisfactory agreement with data obtained by modeling using realistic particle size distributions. Observed Ka- and W-band LDR differences are strongly affected by the radar hardware system polarization “leak” and are generally less than 4 dB. Smaller differences are observed for higher depolarizations, where the polarization “leak” is less pronounced. Realistic assumptions about particle canting and the system polarization isolation lead to modeling results that satisfactorily agree with observational dual-frequency LDR data.

42 ENGINEERING↗
DOE OSTI2026

A point-particle-based hydride shell-shedding model for ejecta particle transport in reactive environments

A shock wave passing over a rough or perturbed metal surface will induce a limiting case of Richtmyer–Meshkov instability and will cause small particles to eject from the surface and transport into the surrounding medium. These particles are known as ejecta and can be either solid or liquid in nature. Recent experiments have shown that liquid cerium ejecta clouds exhibit unexpected non-monotonic acceleration behaviors as well as temperature plateaus after a brief temperature rise if they are transporting in a chemically reactive, hydrogen-based medium while they act as expected in an inert medium. This work details a point-particle model developed for reactive cerium ejecta transport, which attempts to account for these new physics through the behavior of a developing solid hydride shell, which is believed to form as a product of the reaction. The overall model incorporates the effects of the reaction on the particle properties as well as the effects of potential shedding of the shell into sub-micrometer scale flakes and potential phase change of the hydride if the ejecta particles reach the melt point of the hydride layer. The model is tested by performing simulations of the original motivating experiments and comparing quantities, such as ejected mass, velocimetry, and temperature profiles, against the experimental data. While the model is able to capture many general features of the observed anomalies, some inaccuracies still exist. These point to both missing physics in the model (such as a deuterium adsorption mechanism on the hydride layer) as well as a lack of knowledge of certain material properties (such as the strength of cerium hydride to determine dynamic fracture thicknesses) needed to fully reduce the uncertainties in the model by up to an order of magnitude and perform a true attempt at model validation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Influences of the South American Low-Level Jet on the Convective Environment in Central Argentina Using a Convection-Permitting Simulation

This study documents the spatial and temporal distribution of the South American low-level jet (SALLJ) and quantifies its impact on the convective environment using a 6.5-month convection-permitting simulation during the Remote Sensing of Electrification, Lightning, And Mesoscale/Microscale Processes with Adaptive Ground Observations and Clouds, Aerosols, and Complex Terrain Interactions (RELAMPAGO-CACTI) campaigns. Overall, the simulation reproduces the observed SALLJ characteristics in central Argentina near the Sierras de Córdoba (SDC), a focal point for terrain-focused upscale growth. SALLJs most frequently occur in the summer with maxima to the northwest and east of the SDC and minima over the higher terrain. The shallower SALLJs (<1750 m) have a strong overnight skew, while the elevated jets are more equally spread throughout the day. SALLJ periods often have higher amounts of low-level moisture and instability compared to non-SALLJ periods, with these impacts increasing over time when the SALLJ is present and decreasing afterward. The SALLJ may enhance low-level wind shear magnitudes (particularly when accounting for the jet height); however, enhancement is somewhat limited due to the presence of speed shear in most situations. SALLJ periods are associated with low-level directional shear favorable for organized convection and an orientation of cloud-layer wind shear parallel to the terrain, which could favor upscale growth. A case study is shown in which the SALLJ influenced both the magnitude and direction of wind shear concurrent with convective upscale growth near the SDC. This study highlights the complex relationship between the SALLJ and its impacts during periods of widespread convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Sizing ice hydrometeor populations using the dual-wavelength radar ratio

Dual-wavelength (3.2 and 0.32 cm, i.e., X- and W-radar bands) radar ratio (DWR) measurements in ice clouds and precipitation using Canada's National Research Council Institute for Aerospace Research airborne radar are compared to closely collocated particle microphysical in situ sampling data in order to develop relations between DWR and characteristic hydrometeor size. This study uses the radar and in situ data sets collected during the In-Cloud ICing and Large-drop Experiment (ICICLE) campaign in midlatitude frontal clouds. Since atmospheric particle scattering at X band is predominantly in the Rayleigh regime and the W-band frequency is the highest frequency usually used for hydrometeor remote sensing, the X–W-band combination provides a relatively strong dual-wavelength reflectivity difference. This study considers radar and in situ measurements conducted in relatively homogeneous cloud and precipitation conditions. Measurements show that under these conditions, the difference between the X-band radar reflectivities observed with vertical and horizontal pointing of the radar beam are generally small and often negligible. However, W-band reflectivities at vertical beam pointing are, on average, larger than those for horizontal beam pointing by about 4 dB, which is a non-Rayleigh scattering effect from preferentially oriented non-spherical particles. A horizontal radar beam DWR–mean volume particle size relation, D v , provides robust estimates of this characteristic size for populations of particles with different habits. Uncertainties of D v retrievals using DWR are around 0.6 mm when D v is greater than approximately 1 mm. Size estimates using vertical radar beam DWRs have larger uncertainties due to smaller dual-wavelength signals and stronger influences of hydrometeor habits and orientations at this geometry of beam pointing. Mean relations among different characteristic sizes, which describe the entire particle size distribution (PSD), such as D v , and other sizes used in various applications (e.g., the mean, effective, and median sizes) are derived, so the results of this study can be used for estimating different PSD characteristic sizes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Oscilloscope Data Push Program

This paper details the development of a Python program designed to automate the data acquisition and conversion for an oscilloscope for the purposes of a one-off/temporary data acquisition system for users that readily need data, and do not have the option of obtaining a Data Acquisition (DAQ) solution. Creating DAQ systems for analyzing a system requires expensive electronics and a dedicated team of engineers for support. Traditionally, manual data collection and processing are time consuming and prone to error. By automating these processes, the cost, efficiency and accuracy of data handling are improved upon. This project involves the creation of a program that interacts with the oscilloscope. During this interaction, there are various functions being performed such as the acquisition of waveform data via floating points, generating plots with the acquired wave points, and storing of floating points in a CSV file format for future reference and plotting purposes. While the initial aim of the project included continuous logging to a cloud database, this was deferred due to time constraints. The results portrayed an almost-instant rate of data collection with a buffer time, showcasing the potential for further integration and real-time data processing.

Osei-Tutu, Jason↗
DOE OSTI2020

Comparison of Eulerian Bin and Lagrangian Particle-Based Microphysics in Simulations of Nonprecipitating Cumulus

A single nonprecipitating cumulus congestus setup is applied to compare droplet spectra grown by the diffusion of water vapor in Eulerian bin and particle-based Lagrangian microphysics schemes. Bin microphysics represent droplet spectral evolution applying the spectral density function. In the Lagrangian microphysics, computational particles referred to as superdroplets are followed in time and space with each superdroplet representing a multiplicity of natural cloud droplets. The same cloud condensation nuclei (CCN) activation and identical representation of the droplet diffusional growth allow the comparison. The piggybacking method is used with the two schemes operating in a single simulation, one scheme driving the dynamics and the other one piggybacking the simulated flow. Piggybacking allows point-by-point comparison of droplet spectra predicted by the two schemes. Overall, the results show the impact of inherent limitations of the two microphysics simulation methods, numerical diffusion in the Eulerian scheme and a limited number of superdroplets in the Lagrangian scheme. Numerical diffusion in the Eulerian scheme results in a more dilution of the cloud upper half and thus smaller cloud droplet mean radius. The Lagrangian scheme typically has larger spatial fluctuations of droplet spectral properties. A significantly larger mean spectral width in the bin microphysics across the entire cloud depth is the largest difference between the two schemes. A fourfold increase of the number of superdroplets per grid volume and a twofold increase of the spectral resolution and thus the number of bins have small impact on the results and provide only minor changes to the comparison between simulated cloud properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

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↗