Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian particle tracking”

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.

89 records · Page 5

Ice Cloud Formation and Dehydration in the Tropical Tropopause Layer

Stratospheric water vapor is important not only for its greenhouse forcing, but also because it plays a significant role in stratospheric chemistry. Several recent studies have focused on the potential for dehydration due to ice cloud formation in air rising slowly through the tropical tropopause layer (TTL). Holton and Gettelman showed that temperature variations associated with horizontal transport of air in the TTL can drive ice cloud formation and dehydration, and Gettelman et al. recently examined the cloud formation and dehydration along kinematic trajectories using simple assumptions about the cloud properties. In this study, a Lagrangian, one-dimensional cloud model has been used to further investigate cloud formation and dehydration as air is transported horizontally and vertically through the TTL. Time-height curtains of temperature are extracted from meteorological analyses. The model tracks the growth, advection, and sedimentation of individual cloud particles. The regional distribution of clouds simulated in the model is comparable to the subvisible cirrus distribution indicated by SAGE II. The simulated cloud properties and cloud frequencies depend strongly on the assumed supersaturation threshold for ice nucleation. The clouds typically do not dehydrate the air along trajectories down to the temperature minimum saturation mixing ratio. Rather the water vapor mixing ratio crossing the tropopause along trajectories is 10-50% larger than the saturation mixing ratio. I will also discuss the impacts of Kelvin waves and gravity waves on cloud properties and dehydration efficiency. These simulations can be used to determine whether observed lower stratospheric water vapor mixing ratios can be explained by dehydration associated with in situ TTL cloud formation alone.

Jensen, Eric↗

Numerical investigation of inertial sphere impaction on a cylinder in subcritical flow

We investigated the accuracy of a Lagrangian scheme for simulating both the motion and wall impaction of neutrally buoyant inertial spheres (8.85 ≤ St ≤ 18.04) moving in the subcritical regime (Re = 10,972 and 22,366) of flows past a circular cylinder. The accuracy of the eddy-resolving simulations of the flow field was verified based on available data at nearly comparable flow conditions. The accuracy of the eddy-resolving simulations of the flow field was verified based on available data at nearly comparable flow conditions. The accuracy of the particle tracking scheme was verified with respect to corresponding experimental data in a laboratory flume. Particle dispersion en route to the circular cylinder is mostly driven by the flow disturbance from the particle injector tube. The trends in collision rates are largely explained by the rate of particle entertainment into the boundary layer around the cylinder and the inertial condition of the spheres in this proximity region—both of which can be better assessed by post-processing of the numerical outcomes. Here, the primary incentive to undertake this work was the need for validating a particle-based method for quantifying the hydraulic stressors leading to bodily damage of fish passing flows influenced by hydraulic infrastructure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Assessment of Oceanographic Conditions during the North Atlantic EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) Field Campaign

This manuscript presents an overview of NASA’s EXport Processes in the Ocean from Remote Sensing 2021 Field Campaign in the North Atlantic (EXPORTS NA) and provides quantitative and dynamical descriptions of the physical processes modulating water transformations during the study. A major programmatic goal was to conduct the sampling in a Lagrangian mode so that ocean ecological and biogeochemical changes can be observed independent from physical advective processes. To accomplish this goal, EXPORTS NA conducted a multi-ship, multi-asset field sampling program within a retentive, anticyclonic mode water eddy. Beneath depths of ∼100 m, Lagrangian sampling assets remained within the eddy core waters (ECWs) throughout the experiment, demonstrating that the ECWs within the mode water eddy were retentive. However, strong westerly winds from four storm events deepened the mixed layer (ML) of the surface core waters (SCWs) above the eddy’s mode water core by 25-40 m and exchanged some of the SCWs with surface waters outside of the eddy via Ekman transport. Estimates of flushing times ranged from 5-8 days, with surface exchange fractions ranging from 20-75% and were consistent with particle tracking advected by combined geostrophic and Ekman velocities. The relative contributions of horizontal and vertical advection on changes in ECW tracers depended on the horizontal and vertical gradients of that tracer. For example, horizontal advection played a large role in ECW salinity fluxes, while vertical entrainment played a larger role in the fluxes of nutrients into SCW ML. Each storm injected nutrients and low oxygen waters into the ML, after which the surface ocean ecosystem responded by reducing nutrient concentrations and increasing %O2 saturation levels. Overall, ECW values of chlorophyll and POC were the largest at the onset of the field program and decreased throughout the campaign. The analysis presented provides a physical oceanographic context for the many measurements made during the EXPORTS NA field campaign while illustrating the many challenges of conducting a production-flux experiment, even in a Lagrangian frame, and the inherent uncertainties of interpreting biological carbon pump observations that were collected in a Eulerian frame of reference.

oceanography↗

Pore-Scale Modeling of Mineral Growth and Nucleation in Reactive Flow

A fundamental understanding of mineral precipitation kinetics relies largely on microscopic observations of the dynamics of mineral surfaces exposed to supersaturated solutions. Deconvolution of tightly bound transport, surface reaction, and crystal nucleation phenomena still remains one of the main challenges. Particularly, the influence of these processes on texture and morphology of mineral precipitate remains unclear. This study presents a coupling of pore-scale reactive transport modeling with the Arbitrary Lagrangian-Eulerian approach for tracking evolution of explicit solid interface during mineral precipitation. It incorporates a heterogeneous nucleation mechanism according to Classical Nucleation Theory which can be turned “on” or “off.” This approach allows us to demonstrate the role of nucleation on precipitate texture with a focus at micrometer scale. In this work precipitate formation is modeled on a 10 micrometer radius particle in reactive flow. The evolution of explicit interface accounts for the surface curvature which is crucial at this scale in the regime of emerging instabilities. The results illustrate how the surface reaction and reactive fluid flow affect the shape of precipitate on a solid particle. It is shown that nucleation promotes the formation of irregularly shaped precipitate and diminishes the effect of the flow on the asymmetry of precipitation around the particle. The observed differences in precipitate structure are expected to be an important benchmark for reaction-driven precipitation in natural environments.

58 GEOSCIENCES↗

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗

Markovian transport processes in a heterogeneous, variably saturated watershed: A multi-domain spatial Markov model

Transport processes in watersheds remain challenging to describe accurately with upscaled models because of their intricate complexities at multiple scales that can lead to multi-modal, non-Fickian breakthrough curves. One of the recent advances in solute transport modeling has been the Lagrangian spatial Markov model (SMM), which describes transport using velocity increments and models the correlated transitions between them. The method has been applied previously to saturated porous media and fractured media, but these successful applications suggest it may be useful for describing watershed scale processes as well. An existing, 3-D, heterogeneous variably saturated model of a headwaters catchment and classical random walk particle tracking were used to generate a forward simulation and the resulting transport was analyzed with a spatial Markov framework. Existing SMM methods were used to analyze the transport paths through the complex, unconfined domain and we found that the vadose zone plays an important role in the correlation structure. As such, a new kind of SMM is proposed, termed a multi-domain SMM (MD-SMM), that uses separate correlation models for auto-transitions in the vadose zone and the saturated zone, and for the cross-transitions between them. The MD-SMM improves representation of transport relative to the single domain model and can naturally delineate times when a particle is in the vadose zone or the saturated zone, adding an additional degree of realism to the upscaled model. Lastly, we show that the transition probability matrix for transport in the saturated portion of this unconfined system is approximately Markovian, which is an important validation of the proposed framework.

54 ENVIRONMENTAL SCIENCES↗

LeWRON: Agentic Analysis of Electroweak Phase Transitions

The electroweak phase transition (EWPT) is a central topic in particle physics and cosmology, connecting collider phenomenology, baryogenesis, and gravitational-wave observatories. Its analysis requires a technically demanding, convention-sensitive, and model-dependent pipeline, from constructing the finite-temperature effective potential to tracking thermal histories, computing bubble nucleation rates, and predicting gravitational-wave spectra. We present LeWRON (Learning ElectroWeak phase tRansitiON), an agentic framework that orchestrates this pipeline starting from an input Lagrangian. LeWRON combines audited toolbox construction with an Explorer module that uses the generated model-specific code for further analysis, including scans and plots. Intermediate analytic outputs are checked by auditor agents and stored as structured artifacts, enabling reproducible human inspection and downstream use through both a command-line interface and a public Python API. The framework supports a reproduction mode, which infers conventions from the literature and reproduces published results, and a discovery mode, which guides users through structured checkpoints for new models. We demonstrate LeWRON across representative beyond-the-Standard-Model scenarios and release the code on GitHub.

Wang, Isaac R. [Fermilab] (ORCID:000000030789218X)↗

Long range frequency chirping of Alfvén eigenmodes

In this work, a theoretical framework has been developed for an NBI scenario to model the hard non-linear evolution of global Alfvén eigenmodes (GAEs) where the adiabatic motion of phase-space structures (holes and clumps), associated with frequency chirping, occurs in generalized phase-space of slowing down energetic particles. The radial profile of the GAE is expanded using finite elements, which allows an update of the mode structure as the mode frequency chirps. Constants of motion are introduced to track the dynamics of energetic particles during frequency chirping by implementing proper action-angle variables and canonical transformations which reduce the dynamics essentially to 1D. Consequently, we specify whether the particles are drifting inward/outward as the frequency deviates from the initial MHD eigenfrequency. Using the principle of least action, we have derived the non-linear equation describing the evolution of the radial profile by varying the total Lagrangian of the system with respect to the weights of the finite elements. For the choice of parameters in this work, it is shown that the peak of the radial profile is shifted and also broadens due to frequency chirping. The time rate of frequency change is also calculated using the energy balance and we show that the adiabatic condition remains valid once it is satisfied. This model clearly illustrates the theoretical treatment to study the long-range adiabatic frequency sweeping events observed for Alfvén gap modes in real experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Face that Launched a Thousand Slips

Helene, (approximately 17.6 kilometers mean radius) is an L4 Trojan co-orbital of Saturn's moon Dione. Its hemisphere features an unusual morphology consisting of broad depressions and a generally smooth surface patterned with streaks and grooves. The streaks appear to be oriented down-gradient, as are the grooves. This pattern suggests intensive mass-wasting as a dominant process on the leading hemisphere. Kilometer-scale impact craters are very sparse on the leading hemisphere other than the degraded kilometer-scale basins defining the overall satellite shape, and many small craters have a diffuse appearance suggesting ongoing mass wasting. Thus mass wasting must dominate surface-modifying processes at present. In fact, the mass wasting appears to have been sufficient in magnitude to narrow the divides between adjacent basins to narrow septa, similar, but in lower relief, to the honeycomb pattern of Hyperion. The prominent groves occur primarily near topographic divides and appear have cut into a broad, slightly lower albedo surface largely conforming to the present topography but elevated a few meters above the smooth surfaces undergoing mass wasting flow. Low ridges and albedo markings on the surface suggest surface flow of materials traveling up to several kilometers. Diffusive mass wasting produces smooth surfaces - such a pattern characterizes most of the low-lying surfaces. The grooves, however, imply that the transport process is advective at those locations where they occur, that is, erosion tends to concentrate along linear pathways separated by divides. In fact, in many places grooves have a fairly regular spacing of 125-160 meters, defining a characteristic erosional scale. Several questions are prompted by the unusual morphology of Helene: 1) What is the nature of the surface materials? 2) Are the transport processes gradual or catastrophic motion from one or a few events? 3) What mechanisms drive mass wasting and groove development? 4) Have the formative processes been active in the recent past? 5) Finally, is the surface accreting or eroding? The smooth character of the leading edge hemisphere of Helene and the dominance of mass wasting suggest that the surface is composed of fine-grained debris, probably dominated by dust-size to small gravel particles. The Lagrangian points of saturnian satellites are locations where planetesimals might have accreted to form co-orbital satellites such as and they may capture ejecta from their master moon. Published models suggest that Helene is a site of net accretion, but we find no extant explanation for the dominance of fine grain sizes for the surface (and probable subsurface) composition of Helene and the other Lagrangian satellites. Observation of the mass wasting tracks on Helene suggests the presence of well-defined streams of debris with low bordering levees that may be depositional features or remnants of the dissected higher surface. Some flows in grazing illumination appear to have a convex cross-section. This mass-flow morphology might be consistent with the occurrence of large-scale flow events, but which might have occurred through rapid emplacement or slow glacier-like creep. On the other hand, small craters appear to have been "softened" by creep-like processes, indicating ongoing mass wasting.

honeycomb pattern of Hyperion↗

Comparison of Integrated Radiation Transport Models with TEPC Measurements for the Average Quality Factors in Spaceflights

The purpose of this work is to test our theoretical model for the interpretation of radiation data measured in space. During the space missions astronauts are exposed to the complex field of radiation type and kinetic energies from galactic cosmic rays (GCR), trapped protons, and sometimes solar particle events (SPEs). The tissue equivalent proportional counter (TEPC) is a simple time-dependent approach for radiation monitoring for astronauts on board the International Space Station. Another and a newer approach to Microdosimetry is the use of silicon-on-insulator (SOI) technology launched on the MidSTAR-1 mission in low Earth orbit (LEO). In the radiation protection practice, the average quality factor of a radiation field is defined as a function of linear energy transfer (LET), Q(sub ave)(LET). However, TEPC measures the average quality factor as a function of the lineal energy y, Q(sub ave)(y), defined as the average energy deposition in a volume divided by the average chord length of the volume. Lineal energy, y, deviates from LET due to energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector volume. Monte Carlo track structure simulation was employed to obtain the response of a TEPC irradiated with charged particle for an equivalent site diameter of 1 micron of wall-less counter. The calculated data of the energy absorption in the wall-less counter were compiled for various y values for several ion types at various discrete projectile energy levels. For the simulation of TEPC response from the mixed radiation environments inside a spacecraft, such as, Space Shuttle and International Space Station, the complete microdosimetric TEPC response, f( y, E, Z), were calculated with the Monte Carlo theoretical results by using the first order Lagrangian interpolation for a monovariate function at a given y value (y = 0.1 keV/micron 5000 keV/micron) at any projectile energy level (E = 0.01 MeV/u to 50,000 MeV/u) of each specific radiation type (Z = 1 to 28). Because the anomalous response has been observed at large event sizes in the experiment due to the escape of energy out of sensitive volume by delta-rays and the entry of delta-rays from the high-density wall into the low-density gas-volume cavity, Monte Carlo simulation was also made for the response of a walled-TEPC with wall thickness 2 mm and density 1 g/cm(exp 3). The radius of cavity was set to 6.35 mm and a gas density 7.874 x 10(exp -5) g/cm(exp 3). The response of the walled- and the wall-less counters were compared. The average quality factor Q(sub ave)(y) for trapped protons on STS-89 demonstrated the good agreement between the model calculations and flight TEPC data as shown. Using an integrated space radiation model (this includes the transport codes HZETRN and BRYNTRN, the quantum nuclear interaction model QMSFRG) and the resultant response distribution functions of walled-TEPC from Monte-Carlo track simulations, we compared model calculations with walled-TEPC measurements from NASA missions in LEO and made predictions for the lunar and the Mars missions. The Q(sub ave)(y) values for the trapped or the solar protons ranged from 1.9-2.5. This over-estimates the Qave(LET) values which ranged from 1.4-1.6. Both quantities increase with shield thickness due to nuclear fragmentation. The Q(sub ave)(LET) for the complete GCR spectra was found to be 3.5-4.5, while flight TEPCs measured 2.9-3.4 for Q(sub ave)(y). The GCR values are decreasing with the shield thickness. Our analysis for a proper interpretation of data supports the use of TEPCs for monitoring space radiation environment.

Kim, Myung-Hee Y.↗

A method for generating moving, orthogonal, area preserving polygonal meshes

A new method for generating locally orthogonal polygonal meshes from a set of generator points is presented in which polygon areas are a constraint. The area constraint property is particularly useful for particle methods where moving polygons track a discrete portion of material. Because Voronoi polygon meshes have some very attractive mathematical and numerical properties for numerical computation, a generalization of Voronoi polygon meshes was formulated that enforces a polygon area constraint. Area constrained moving polygonal meshes allow one to develop hybrid particle-mesh numerical methods that display some of the most attractive features of each approach. It is shown that this mesh construction method can continuously reconnect a moving, unstructured polygonal mesh in a pseudo-Lagrangian fashion without change in cell area/volume, and the method's ability to simulate various physical scenarios is shown. Overall, the advantages are identified for incompressible fluid flow calculations, with demonstration cases that include material discontinuities of all three phases of matter and large density jumps.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Oklahoma Health and Air Quality: Mapping Air Quality Using NASA Earth Observations to Investigate Recent Increases in Ozone Concentrations

Tropospheric ozone (O3) is formed by anthropogenic pollutants interacting with sunlight and is considered harmful to human health in high concentrations. In the summer of 2018, the Oklahoma Department of Environmental Quality (DEQ) measured unexpected spikes in O3 in Seiling, Oklahoma, with concentrations exceeding those measured in bustling Oklahoma City and Tulsa. The DEQ tracks air quality using ground monitors and does not utilize Earth observation data in its monitoring or analysis. This project used remotely sensed data to investigate these2018 air quality anomalies, identifying possible causes. We analyzed atmospheric data from Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS), and Sentinel-5P Tropospheric Ozone Monitoring Instrument (TROPOMI) in conjunction with ground-based measurements of tropospheric ozone (O3), nitrogen dioxide (NO2), methane (CH4), carbon monoxide (CO), formaldehyde (HCHO) and aerosol optical depth (AOD). We compared Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model simulations and Earth observation visualizations to pinpoint ozone spike causes. We also generated models to identify contributing factors to variations in ground ozone concentrations in our study area. The results point to a variety of ozone spike causes, primarily from outside of the state, and support the placement of additionalNO2, O3, and CO monitors to the southeast of Seiling. These analyses can help guide the placement of future monitors in the ground monitoring network and inform air quality regulations in Oklahoma.

Carolina Rosales↗

Tomographic X-ray particle tracking velocimetry

We investigate the feasibility of in-laboratory tomographic X-ray particle tracking velocimetry (TXPTV) and consider creeping flows with nearly density matched flow tracers. Specifically, in these proof-of-concept experiments we examined a Poiseuille flow, flow through porous media and a multiphase flow with a Taylor bubble. For a full 360° computed tomography (CT) scan we show that the specially selected 60 micron tracer particles could be imaged in less than 3 seconds with a signal-to-noise ratio between the tracers and the fluid of 2.5, sufficient to achieve proper volumetric segmentation at each time step. In the pipe flow, continuous Lagrangian particle trajectories were obtained, after which all the standard techniques used for PTV or PIV (taken at visible wave lengths) could also be employed for TXPTV data. And, with TXPTV we can examine flows inaccessible with visible wave lengths due to opaque media or numerous refractive interfaces. In the case of opaque porous media we were able to observe material accumulation and pore clogging, and for flow with Taylor bubble we can trace the particles and hence obtain velocities in the liquid film between the wall and bubble, with thickness of liquid film itself also simultaneously obtained from the volumetric reconstruction after segmentation. While improvements in scan speed are anticipated due to continuing improvements in CT system components, we show that for the flows examined even the presently available CT systems could yield quantitative flow data with the primary limitation being the quality of available flow tracers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Small Mission Concept to the Sun-Earth Lagrangian L5 Point for Innovative Solar, Heliospheric and Space Weather Science

We present a concept for a small mission to the Sun-Earth Lagrangian L5 point for innovative solar, heliospheric and space weather science. The proposed INvestigation of Solar-Terrestrial Activity aNd Transients (INSTANT) mission is designed to identify how solar coronal magnetic fields drive eruptions, mass transport and particle acceleration that impact the Earth and the heliosphere. INSTANT is the first mission designed to (1) obtain measurements of coronal magnetic fields from space and (2) determine coronal mass ejection (CME) kinematics with unparalleled accuracy. Thanks to innovative instrumentation at a vantage point that provides the most suitable perspective view of the Sun-Earth system, INSTANT would uniquely track the whole chain of fundamental processes driving space weather at Earth. We present the science requirements, payload and mission profile that fulfill ambitious science objectives within small mission programmatic boundary conditions.

Lavraud, B.↗

A Computational Information Criterion for Particle-Tracking with Sparse or Noisy Data

Traditional probabilistic methods for the simulation of advection-diffusion equations (ADEs) often overlook the entropic contribution of the discretization, e.g., the number of particles, within associated numerical methods. Many times, the gain in accuracy of a highly discretized numerical model is outweighed by its associated computational costs or the noise within the data. Herein, we address the question of how many particles are needed in a simulation to best approximate and estimate parameters in one-dimensional advective-diffusive transport. To do so, we use the well-known Akaike Information Criterion (AIC) and a recently-developed correction called the Computational Information Criterion (COMIC) to guide the model selection process. Random-walk and mass-transfer particle tracking methods are employed to solve the model equations at various levels of discretization. Numerical results demonstrate that the COMIC provides an optimal number of particles that can describe a more efficient model in terms of parameter estimation and model prediction compared to the model selected by the AIC even when the data is sparse or noisy, the sampling volume is not uniform throughout the physical domain, or the error distribution of the data is non-IID Gaussian.

97 MATHEMATICS AND COMPUTING↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

Reactive particle-tracking solutions to a benchmark problem on heavy metal cycling in lake sediments

Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.

54 ENVIRONMENTAL SCIENCES↗