Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase space dynamics”

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

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

At least 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Approximate symmetries of guiding-centre motion

In a strong, inhomogeneous magnetic field, charged particle dynamics may be studied in the guiding-centre approximation, which is known to be Hamiltonian. When the magnetic field is quasisymmetric, the first-order guiding-centre (FGC) Hamiltonian structure admits a continuous symmetry, and therefore a conserved quantity in addition to the energy. Since the FGC system is only an approximation, it is also interesting to consider approximate symmetries of the guiding-centre Hamiltonian structure. We find that any approximate spatial symmetry coincides with quasisymmetry to leading order. For approximate phase-space symmetries, we derive weaker conditions than quasisymmetry. The latter include 'weak quasisymmetry' as a subcase, recently proposed by Rodríguez et al. Our results, however, show that weak quasisymmetry is necessarily non-spatial at first order. Finally, we demonstrate that if the magnetic field is constrained to satisfy magnetohydrostatic force balance then an approximate symmetry must agree with quasisymmetry to leading order.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Charged dust in planetary magnetospheres: Hamiltonian dynamics and numerical simulations for highly charged grains

We use a combination of analytical and numerical methods to investigate the dynamics of charged dust grains in planetary magnetospheres. Our emphasis is on obtaining results valid for particles that are not necessarily dominated either by gravitational or electromagnetic forces. A Hamiltonian formulation of the problem yields exact results, for all values of charge-to-mass ratio, when we introduce two constraints: particles remain in the equatorial plane and the magnetic field is taken as axially symmetric. In particular, we obtain locations of equilibrium points, the frequencies of stable periodic orbits, the topology of separatrices in phase space, and the rate of longitudinal drift. These results are significant for specific applications: motion in the nearly aligned dipolar field of Saturn, and the trajectories of arbitrarily charged particles in complex magnetic fields for limited periods of time after ejection from parent bodies. Since the model is restrictive, we also use numerical integrations of the full three-dimensional equations of motion and illustrate under what conditions the constrained problem yields reasonable results. We show that a large fraction of the intermediately charged and highly charged (gyrating) particles will always be lost to a planet's atmosphere within a few hundred hours, for motion through tilted-dipole magnetic fields. We find that grains must have a very high charge-to-mass ratio in order to be mirrored back to the ring plane. Thus, except perhaps at Saturn where the dipole tilt is very small, the likely inhabitants of the dusty ring systems are those particles that are either nearly Keplerian (weakly charged) grains or grains whose charges place them in the lower end of the intermediate charge zone. Fianlly, we demonstrate the effect of plasma drag on the orbits of gyrating particles to be a rapid decrease in gyroradius followed by a slow radial evolution of the guiding center.

Schaffer, L.↗
NASA NTRS2003

Defect, Kinetics and Heat Transfer of CDTE Bridgman Growth without Wall Contact

A detached growth mechanism has been proposed, which is similar to that proposed by Duffar et al. and used to study the current detached growth system. From numerical results, we can conclude that detached growth will more likely appear if the growth and wetting angles are large and meniscus is flat. Detached thickness is dependent on growth angle, wetting angle, and gap width and shape of the fins. The model can also explain why the detached growth will not happen for metals in which the growth angle is almost zero. Since the growth angle of CdZnTe cannot be changed, to promote detached growth, the number density of the fins should be low and the wetting angle should be high. Also, a much smaller gap width of the fins should be used in the ground experiment and the detached gap width is much smaller. The shape of the fins has minor influence on detached growth. An integrated numerical model for detached solidification has been developed combining a global heat transfer sub-model and a wall contact sub-model. The global heat transfer sub-model accounts for heat and mass transfer in the multiphase system, convection in the melt, macro-segregation, and interface dynamics. The location and dynamics of the solidification interface are accurately tracked by a multizone adaptive grid generation scheme. The wall contact sub-model accounts for the meniscus dynamics at the three-phase boundary. Simulations have been performed for crystal growth in a conventional ampoule and a designed ampoule to understand the benefits of detached solidification and its impacts on crystalline structural quality, e.g., stoichiometry, macro-segregation, and stress. From simulation results, both the Grashof and Marangoni numbers will have significant effects on the shape of growth front, Zn concentration distribution, and radial segregation. The integrated model can be used in designing apparatus and determining the optimal geometry for detached solidification in space and on the ground.

Larson, D. J., Jr.↗
DOE OSTI2023

Real-Space Infrared Spectroscopy of Ferroelectric Domain Walls in Multiferroic h -(Lu,Sc)FeO 3

Herein we employ synchrotron-based near-field infrared spectroscopy to image the phononic properties of ferroelectric domain walls in hexagonal ( h ) Lu 0.6 Sc 0.4 FeO 3 , and we compare our findings with a detailed symmetry analysis, lattice dynamics calculations, and prior models of domain-wall structure. Rather than metallic and atomically thin as observed in the rare-earth manganites, ferroelectric walls in h -Lu 0.6 Sc 0.4 FeO 3 are broad and semiconducting, a finding that we attribute to the presence of an A -site substitution-induced intermediate phase that reduces strain and renders the interior of the domain wall nonpolar. Mixed Lu/Sc occupation on the A site also provides compositional heterogeneity over micron-sized length scales, and we leverage the fact that Lu and Sc cluster in different ratios to demonstrate that the spectral characteristics at the wall are robust even in different compositional regimes. This work opens the door to broadband imaging of physical and chemical heterogeneity in ferroics and represents an important step toward revealing the rich properties of these flexible defect states.

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

Goal-oriented real-time Bayesian inference for linear autonomous dynamical systems with application to digital twins for tsunami early warning

We present a goal-oriented framework for constructing digital twins with the following properties: (1) they employ discretizations of high-fidelity partial differential equation (PDE) models governed by autonomous dynamical systems, leading to large-scale forward problems; (2) they solve a linear inverse problem to assimilate observational data to infer uncertain model components followed by a forward prediction of the evolving dynamics; and (3) the entire end-to-end, data-to-inference-to-prediction computation is carried out without approximation and in real time through a Bayesian framework that rigorously accounts for uncertainties. Several challenges must be overcome to realize this framework, including the large scale of the forward problem, the high dimensionality of the parameter space, and for a class of problems including those we target, the slow decay of the singular values of the parameter-to-observable map. Here we introduce a methodology to overcome these challenges by exploiting the autonomous structure of the forward model to decompose the solution of the inverse problem into a one-time-only offline phase in which the PDE model is solved a limited number of times (equal to the number of sensors), and an online phase that maps well onto GPUs and computes the parameter inference and prediction of quantities of interest in real time, given observational data. Our ultimate goal is to apply this framework to construct digital twins for subduction zones, including Cascadia, to provide early warning for tsunamis generated by megathrust earthquakes. To this end, we demonstrate how our methodology can be used to employ seafloor pressure observations, along with the coupled acoustic–gravity wave equations, to infer the earthquake-induced spatiotemporal seafloor motion (discretized with $\mathscr{O}$ (10 9 ) parameters) and forward predict the tsunami propagation. We present results of an end-to-end inference, prediction, and uncertainty quantification for a representative test problem with $\mathscr{O}$ (10 8 ) inversion parameters for which goal-oriented Bayesian inference is accomplished exactly and in real time, that is, in a matter of seconds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The Dynamics of Oxygen Ion Exchange in Epitaxial Strontium Cobaltite Bilayers

Abstract The exchange of ions across interfaces is key to the field of iontronics, where the properties of the device can be altered by the local ion concentration. This study investigates a complex oxide system where structural and electronic phase transitions can be driven by changes in the concentration of oxygen ions. In situ coherent X‐ray studies are conducted on epitaxial bilayers of insulating SrCoO 2.5 and metallic SrCoO 3 − δ . The diffusion of oxygen ions across the bilayer is studied with X‐ray photon correlation spectroscopy to capture the dynamical behavior of the interface in reducing and oxidizing environments. The behavior is strongly asymmetric, with much slower dynamics appearing in reducing versus oxidizing environments. According to the correlation functions determined from different points in reciprocal space, this study finds that the dynamics near the center of the SrCoO 2.5 crystal are generally similar to those near the heterointerfaces. The results suggest that the interface is stable and reversible, making SrCoO x a model system for the study of iontronic behavior.

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

Rate-dependent strength and deformation heterogeneity of B4C-reinforced Al composite: Time-resolved imaging with synchrotron X-rays

Dynamic (up to 5500 s -1 ) and quasi-static compression tests are conducted on a 15 wt% B4C particle-reinforced Al (B 4 C/Al) composite. In situ, high-speed synchrotron X-ray phase contrast imaging and digital image correlation are employed to map mesoscale deformation fields at μm and μs scales. The bulk stress–strain curves show significant strain and strain-rate hardening under dynamic compression. The strain-rate sensitivity exponent is an order of magnitude higher at high strain rates (> 10 3 s -1 ) than that at low strain rates (< 10 -2 s -1 ). Strain field mapping demonstrates distinct compressive strain localizations for both quasi-static and dynamic loading. Nevertheless, compressive strain localizations appear denser in spacing under dynamic loading, owing to spontaneous dislocation nucleation in both weak and strong zones. This results in a higher density of geometrically necessary dislocations, which contributes to the higher strain and strain-rate hardening of B 4 C/Al under dynamic loading. Additionally, the ratio of the maximum local strain to the bulk average is 1.5, and the local strain-rate enhancement cannot explain the increased rate sensitivity of B 4 C/Al under dynamic loading. Therefore, the rate-dependent deformation heterogeneity dominates the strain-rate hardening of B 4 C/Al. Postmortem analyses help correlate deformation features to particle- and grain-scale microstructures, yielding consistent results with mesoscale strain fields.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Influence of alloy solidification path on melt pool behavior in additive manufacturing

Numerical models used to study transport phenomena in laser-powder bed fusion processes often rely on assumptions and simplifications to reduce their computational expense. One common simplification is in the description of latent heat evolution during the solid-liquid phase change (i.e., the solidification pathway), justified by the fact that the mushy zone thickness is similar to the numerical grid spacing used for continuum transport models. The lack of resolution of transport phenomena in the mushy zone motivates the use of computationally convenient solidification paths such as linear or sigmoidal relationships over pathways derived from fundamental solidification theory such as equilibrium or Scheil models. In the present work, an uncertainty quantification (UQ) framework is used to analyze the influence of solidification pathway selection on the solidification dynamics and melt pool geometries in laser based additive manufacturing (AM) of IN625. Results show the solidification pathway has a quantifiable influence on the cooling rate at the liquidus isotherm, mushy zone thickness, and solidification time. Due to similarities in the latent heat evolution at the beginning of solidification, the equilibrium and Scheil models predict similar cooling rates near the liquidus isotherm, however the wider freezing range of Scheil leads to a wider mushy zone compared to equilibrium. The non-physical latent heat release profiles of sigmoidal and linear paths lead to significant overpredictions of cooling rates at the liquidus isotherm compared to equilibrium and Scheil. Finally, these results indicate that careful consideration should be given to the choice of solidification pathway to ensure reliable model predictions.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Vibration of a large space beam under gravity effect

The structural characteristics of a large simply supported beam subjected to gravity are described. The nonlinear governing equations for both the static and the dynamic response are derived and solved analytically. The results show the feasibility of verifying the on-orbit dynamic characteristics of a large space beam by utilizing ground test data of such a structure. It is noted that the gravity effect interacts mostly with the first vibration mode. It was also found that the system of a large space beam subjected to its own weight is a hardening system. The differential equation for the asymmetric mode is a Duffing type equation. However, the governing equation for the symmetric mode has an additional quadratic term. It is this term that causes the maximum vibration amplitudes at different phases to be non-identical.

Shih, C.-F.↗
NASA NTRS1992

A status of the Turbine Technology Team activities

The recent activities of the Turbine Technology Team of the Consortium for Computational Fluid Dynamics (CFD) Application in Propulsion Technology is presented. The team consists of members from the government, industry, and universities. The goal of this team is to demonstrate the benefits to the turbine design process attainable through the application of CFD. This goal is to be achieved by enhancing and validating turbine design tools for improved loading and flowfield definition and loss prediction, and transferring the advanced technology to the turbine design process. In order to demonstrate the advantages of using CFD early in the design phase, the Space Transportation Main Engine (STME) turbines for the National Launch System (NLS) were chosen on which to focus the team's efforts. The Turbine Team activities run parallel to the STME design work.

Griffin, Lisa W.↗
NASA NTRS2002

Dynamics of Sheared Granular Materials

This work focuses on the properties of sheared granular materials near the jamming transition. The project currently involves two aspects. The first of these is an experiment that is a prototype for a planned ISS (International Space Station) flight. The second is discrete element simulations (DES) that can give insight into the behavior one might expect in a reduced-g environment. The experimental arrangement consists of an annular channel that contains the granular material. One surface, say the upper surface, rotates so as to shear the material contained in the annulus. The lower surface controls the mean density/mean stress on the sample through an actuator or other control system. A novel feature under development is the ability to 'thermalize' the layer, i.e. create a larger amount of random motion in the material, by using the actuating system to provide vibrations as well control the mean volume of the annulus. The stress states of the system are determined by transducers on the non-rotating wall. These measure both shear and normal components of the stress on different size scales. Here, the idea is to characterize the system as the density varies through values spanning dense almost solid to relatively mobile granular states. This transition regime encompasses the regime usually thought of as the glass transition, and/or the jamming transition. Motivation for this experiment springs from ideas of a granular glass transition, a related jamming transition, and from recent experiments. In particular, we note recent experiments carried out by our group to characterize this type of transition and also to demonstrate/ characterize fluctuations in slowly sheared systems. These experiments give key insights into what one might expect in near-zero g. In particular, they show that the compressibility of granular systems diverges at a transition or critical point. It is this divergence, coupled to gravity, that makes it extremely difficult if not impossible to characterize the transition region in an earth-bound experiment. In the DE modeling, we analyze dynamics of a sheared granular system in Couette geometry in two (2D) and three (3D) space dimensions. Here, the idea is to both better understand what we might encounter in a reduced-g environment, and at a deeper level to deduce the physics of sheared systems in a density regime that has not been addressed by past experiments or simulations. One aspect of the simulations addresses sheared 2D system in zero-g environment. For low volume fractions, the expected dynamics of this type of system is relatively well understood. However, as the volume fraction is increased, the system undergoes a phase transition, as explained above. The DES concentrate on the evolution of the system as the solid volume fraction is slowly increased, and in particular on the behavior of very dense systems. For these configurations, the simulations show that polydispersity of the sheared particles is a crucial factor that determines the system response. Figures 1 and 2 below, that present the total force on each grain, show that even relatively small (10 %) nonuniformity of the size of the grains (expected in typical experiments) may lead to significant modifications of the system properties, such as velocity profiles, temperature, force propagation, and formation shear bands. The simulations are extended in a few other directions, in order to provide additional insight to the experimental system analyzed above. In one direction, both gravity, and driving due to vibrations are included. These simulations allow for predictions on the driving regime that is required in the experiments in order to analyze the jamming transition. Furthermore, direct comparison of experiments and DES will allow for verification of the modeling assumptions. We have also extended our modeling efforts to 3D. The (preliminary) results of these simulations of an annular system in zero-g environment will conclude the presentation.

Kondic, Lou↗
NASA NTRS2001

Twenty Years of Radiation Measurements in Low-Earth Orbit - What Have We Learned Space Radiation Environment?

The advent of the Space Shuttle program has made possible space radiation environment measurements spanning a wide range of altitudes and orbital inclinations over multiple solar cycles. These measurements range from routine integral dose measurements with thermoluminescent dosimeters to particle energy spectra measurements made with a charged particle telescope. This paper will review the new understanding about the space radiation environment gained from this diverse data set. Major findings from these measurements include: estimations of the westward drift rate of the South Atlantic Anomaly (SAA) of 0.28-0.49/y; evidence for a northward component to the SAA drift of 0.08-0.12/y; observation of the formation and decay of the pseudo-stable additional radiation belt following the Mar 1991 SPE and geomagnetic storm with an estimated decay e-folding time of 9-10 months; observation of a local geomagnetic east-west trapped proton exposure anisotropy with an estimated magnitude of 1.6-3.3; demonstration that the trapped proton exposure in low-Earth orbit (LEO) can be reasonably modeled as a power law function of atmospheric density in the SAA region, with best correlations obtained when the exospheric temperature saturates at 938-975 K; the actual solar cycle modulation of trapped proton exposure in LEO is less than predicted by the AP8 model; and the testing and validation of GCR flux models, radiation transport codes, and dynamic geomagnetic cutoff models. Long-term, time-resolved proportional counter measurements made aboard the Mir during the same period provides further demonstration of the solar cycle modulation of the trapped protons at low altitudes - the observed modulation is also well described as power law function of atmospheric density. These data and findings have helped to improve the overall accuracy of pre-mission crew exposure projections using various semi-empirical space environment models, radiation transport codes, and spacecraft radiation shielding models. During the rise phase of solar cycle 22 (1987-1991), the RMS error between preflight exposure projections and measured crew exposure was 73%. For the rise phase of cycle 23 (1997-2001), the preflight exposure projection RMS error has decreased to 23%. The launch and assembly of the Space Station has begun a new era of long-term LEO space environment monitoring. The radiation environment at the Space Station will be monitored with three external charged particle telescopes oriented in the velocity vector, anti-velocity vector, and zenith directions. Data from the telescopes will provide charge, mass, energy, and arrival direction for incident particles with energy to mass ratios of 13- 450 MeV/amu and Z of 1-24. The external environment data will be complimented by measurements from a portable charged particle telescope and proportional counter located inside the vehicle.

Golightly, Michael J.↗
NASA NTRS2013

Status of the Correlation Process of the V-HAB Simulation with Ground Tests and ISS Telemetry Data

The Virtual Habitat (V-HAB) is a dynamic Life Support System (LSS) simulation, created for investigation of future human spaceflight missions. It provides the capability to optimize LSS during early design phases. The focal point of the paper is the correlation and validation of V-HAB against ground test and flight data. In order to utilize V-HAB to design an Environmental Control and Life Support System (ECLSS) it is important to know the accuracy of simulations, strengths and weaknesses. Therefore, simulations of real systems are essential. The modeling of the International Space Station (ISS) ECLSS in terms of single technologies as well as an integrated system and correlation against ground and flight test data is described. The results of the simulations make it possible to prove the approach taken by V-HAB.

Ploetner, P.↗
NASA NTRS1976

Space telescope optical telescope assembly/scientific instruments. Phase B: Preliminary design and program definition study. Volume 2A(3): Astrometry

Wide field measurements, namely, measurements of relative angular separations between stars over a relatively wide field for parallax and proper motion determinations, were made with the third fine guidance sensor. Narrow field measurements, i.e., double star measurements, are accomplished primarily with the area photometer or faint object camera at f/96. The wavelength range required can be met by the fine guidance sensor which has a spectral coverage from 3000 to 7500 A. The field of view of the fine guidance sensor also exceeds that required for the wide field astrometric instrument. Requirements require a filter wheel for the wide field astrometer, and so one was incorporated into the design of the fine guidance sensor. The filter wheel probably would contain two neutral density filters to extend the dynamic range of the sensor and three spectral filters for narrowing effective double star magnitude difference.

Source record↗
NASA NTRS2014

Near-Surface Meteorology During the Arctic Summer Cloud Ocean Study (ASCOS): Evaluation of Reanalyses and Global Climate Models.

Atmospheric measurements from the Arctic Summer Cloud Ocean Study (ASCOS) are used to evaluate the performance of three atmospheric reanalyses (European Centre for Medium Range Weather Forecasting (ECMWF)- Interim reanalysis, National Center for Environmental Prediction (NCEP)-National Center for Atmospheric Research (NCAR) reanalysis, and NCEP-DOE (Department of Energy) reanalysis) and two global climate models (CAM5 (Community Atmosphere Model 5) and NASA GISS (Goddard Institute for Space Studies) ModelE2) in simulation of the high Arctic environment. Quantities analyzed include near surface meteorological variables such as temperature, pressure, humidity and winds, surface-based estimates of cloud and precipitation properties, the surface energy budget, and lower atmospheric temperature structure. In general, the models perform well in simulating large-scale dynamical quantities such as pressure and winds. Near-surface temperature and lower atmospheric stability, along with surface energy budget terms, are not as well represented due largely to errors in simulation of cloud occurrence, phase and altitude. Additionally, a development version of CAM5, which features improved handling of cloud macro physics, has demonstrated to improve simulation of cloud properties and liquid water amount. The ASCOS period additionally provides an excellent example of the benefits gained by evaluating individual budget terms, rather than simply evaluating the net end product, with large compensating errors between individual surface energy budget terms that result in the best net energy budget.

surface energy↗
DOE OSTI2020

Whose Gas is it anyway? Differentiating the Source of a Large Soil Vapor Plume beneath Two Adjacent Waste Sites - 20487

DOE contractor CH2M Hill Plateau Remediation Company is currently responsible for conducting groundwater contamination monitoring at several RCRA treatment, storage, and disposal units located on the Hanford Site in Richland, Washington State. The Nonradioactive Dangerous Waste Landfill treatment, storage, and disposal unit presents a distinct groundwater monitoring problem because of a large multi-contaminant soil vapor plume beneath it that is a likely source of low-level volatile organic compound groundwater contamination. Adjacent to Nonradioactive Dangerous Waste Landfill is the Solid Waste Landfill. Volatile organic compounds are inventory components of both the Nonradioactive Dangerous Waste Landfill and the Solid Waste Landfill. Therefore, it is possible that both sites could be contributing to the soil vapor plume. For regulatory purposes, it is important to differentiate which site is the primary contributor of volatile organic compounds to the plume. An approach was developed to identify the primary volatile organic compound source of the soil vapor plume beneath Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill. The site conceptual model hypothesis of vapor-phase volatile organic compound transport to the dissolved phase in groundwater was tested by a simple mathematical model of vapor/liquid equilibrium concentrations at the groundwater/air interface. Once it was shown that vapor-phase volatile organic compound transport to groundwater was a valid conceptual model for Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill, spatial and statistical methods were used to determine the primary site contributing to the majority of volatile organic compounds to the soil vapor plume. Average groundwater chloroform, tetrachloroethene, and trichloroethene concentrations from Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill monitoring network wells were plotted on maps of the facilities and immediate vicinities and compared to soil vapor sampling probe locations. Principal component analysis and mixing ratios were used to identify source contributions of each treatment, storage, and disposal unit to the plume. Results of the vapor/liquid equilibrium concentrations mathematical model showed that transport phenomena outweigh steady-state equilibria. Estimated vapor/liquid equilibrium concentrations were considerably lower than soil vapor measurements. The results indicate that dynamic vadose zone and groundwater factors such as decreased vapor concentrations with depth, vapor dilution from dispersion in the vadose zone, and advective and diffusional volatile organic compound dilution in groundwater result in groundwater volatile organic compound concentrations much less than would be measured under steady-state equilibrium conditions. Site source contribution differentiation by principal component analysis and mixing ratios was inconclusive using actual soil gas data because of the similarity in concentration values in both datasets for Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill. Similar data populations suggest mixing of the vapor contributions from both sites by dispersion through the soil matrix pore spaces. However, when groundwater volatile organic compound data were compared between the Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill monitoring networks, Solid Waste Landfill mean concentrations were higher, suggesting more vapor-phase volatile organic compound transport to groundwater at those locations. Simulated volatile organic compound soil vapor and groundwater datasets created to test the methods developed for this study show that the method can be successful in source differentiation when significantly different datasets are compared. This paper will describe a method of testing a conceptual model for vapor-phase contaminant transport to groundwater and for differentiating site sources of contaminants comprising a mixed-constituent soil vapor plume. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Null Raychaudhuri: canonical structure and the dressing time

Abstract We initiate a study of gravity focusing on generic null hypersurfaces, non-perturbatively in the Newton coupling. We present an off-shell account of the extended phase space of the theory, which includes the expected spin-2 data as well as spin-0, spin-1 and arbitrary matter degrees of freedom. We construct the charges and the corresponding kinematic Poisson brackets, employing a Beltrami parameterization of the spin-2 modes. We explicitly show that the constraint algebra closes, the details of which depend on the non-perturbative mixing between spin-0 and spin-2 modes. Finally we show that the spin zero sector encodes a notion of a clock, called dressing time, which is dynamical and conjugate to the constraint. It is well-known that the null Raychaudhuri equation describes how the geometric data of a null hypersurface evolve in null time in response to gravitational radiation and external matter. Our analysis leads to three complementary viewpoints on this equation. First, it can be understood as a Carrollian stress tensor conservation equation. Second, we construct spin-0, spin-2 and matter stress tensors that act as generators of null time reparametrizations for each sector. This leads to the perspective that the null Raychaudhuri equation can be understood as imposing that the sum of CFT-like stress tensors vanishes. Third, we solve the Raychaudhuri constraint non-perturbatively. The solution relates the dressing time to the spin-2 and matter boost charge operators. Finally we establish that the corner charge corresponding to the boost operator in the dressing time frame is monotonic. These results show that the notion of an observer can be thought of as emerging from the gravitational degrees of freedom themselves. We briefly mention that the construction offers new insights into focusing conjectures.

Physics↗