Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface profiling”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Overdriven-detonation states produced by spherically diverging waves

A series of experiments are currently underway at the Detonation Science and Technology group, within Los Alamos National Laboratory, to study the overdriven-detonation states achievable in energetic materials as a result of detonation-wave interactions. A multi-component, energetic-material array was designed to study the amplification of velocity and pressure states produced by spherically diverging detonation waves in pentaerythritol tetranitrate (PETN) output charges. The unique geometry provides a low-jitter, highly controlled series of interactions between three independent-detonation inputs. Streak-camera imaging was performed on the output face of PETN pellets ranging in thickness from 2.5 – 10 mm to characterize the resulting breakout profile. Additionally, photonic Doppler velocimetry (PDV) measurements were collected at the output-pellet surface to determine simultaneity within the system. Detonation-wave velocities upwards of 16 mm/µs were measured, as compared to a steady-state detonation velocity of 7.9 mm/µs for the PETN pressing density investigated. Finally, additional experiments are being conducted to measure the pressure amplification generated at key areas of interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Inverse metal-assisted chemical etching of germanium with gold and hydrogen peroxide

Abstract Metal-assisted chemical etching (MACE) is a flexible technique for texturing the surface of semiconductors. In this work, we study the spatial variation of the etch profile, the effect of angular orientation relative to the crystallographic planes, and the effect of doping type. We employ gold in direct contact with germanium as the metal catalyst, and dilute hydrogen peroxide solution as the chemical etchant. With this catalyst-etchant combination, we observe inverse-MACE, where the area directly under gold is not etched, but the neighboring, exposed germanium experiences enhanced etching. This enhancement in etching decays exponentially with the lateral distance from the gold structure. An empirical formula for the gold-enhanced etching depth as a function of lateral distance from the edge of the gold film is extracted from the experimentally measured etch profiles. The lateral range of enhanced etching is approximately 10–20 μ m and is independent of etchant concentration. At length scales beyond a few microns, the etching enhancement is independent of the orientation with respect to the germanium crystallographic planes. The etch rate as a function of etchant concentration follows a power law with exponent smaller than 1. The observed etch rates and profiles are independent of whether the germanium substrate is n-type, p-type, or nearly intrinsic.

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

Resistive wall tearing mode disruptions

Abstract This paper deals with resistive wall tearing mode (RWTM) disruptions. RWTMs are closely related to resistive wall modes. RWTMs are tearing modes whose linear and nonlinear behavior is strongly dependent on the resistive wall outside the plasma. The consequence for ITER, is that the thermal quench timescale could be much longer than previously conjectured. Active feedback stabilization is another possible way to mitigate or prevent RWTM disruptions. Simulations of RWTM disruptions are reviewed for DIII-D and MST. MST has a longer resistive wall time than ITER, and disruptions are not observed experimentally when MST is operated as a standard tokamak. Simulations indicate that the RWTM disruption time scale is longer than the experimental shot time. Edge cooling causes contraction of the current profile, which can destabilize RWTMs. The equilibria studied here have the q = 2 rational surface close to the edge of the plasma, and low current density between the q = 2 surface and the wall. A sequence of low edge current model equilibria has major disruptions only for a resistive, not ideal, wall, and edge q ⩽ 3.4 . This is consistent with regimes of tokamak disruptivity, suggesting that tokamak disruptions caused by edge cooling at low edge q could be RWTMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Electroluminescence Characterization of Recombination in Back Junction Silicon Heterojunction Test Structures: Role of the Inversion Layer

Electroluminescence (EL) characterization technique is used to quantify recombination properties in back junction silicon (Si) heterojunction (HJ) solar cell test structures with interdigitated back contact (IBC) emitters and metal contacts. These test structures are three terminal devices that allow the selective biasing and injection of minority carriers into spatially isolated emitter strips having different widths. This allows the study of minority carriers' current transport, flow, and recombination in these regions, which can, then, be applied to full IBC cells. EL is used to verify that emission is proportional to the minority carriers' concentration in the corresponding regions of the device, thus providing information relevant to minority carriers' collection when compared with laser beam induced current technique in considerably less time and with a higher spatial resolution. Current transport through the inversion layer that occurs due to the band bending at the n crystalline Si and p amorphous Si interface is discussed relying on the distance that minority carriers travel in different configurations and on their diffusion length. Laser defect spots (LDS) were, then, introduced in the emitter regions to analyze the spatial profile of minority carriers as they recombine at the LDS due to their high surface recombination velocity. Ideality factors in different configurations are spatially mapped from EL before and after introducing LDS and are used to show the effect of the inversion layer on the different recombination mechanisms. Current transport through the inversion layer is studied at different biasing conditions and correlated with current transport after LDS.

42 ENGINEERING↗
DOE OSTI2023

A Methodology for Estimating the Energy and Moisture Budget of the Convective Boundary Layer Using Continuous Ground-Based Infrared Spectrometer Observations

Abstract Land–atmosphere interactions play a critical role in both the atmospheric water and energy cycles. Changes in soil moisture and vegetation alter the partitioning of surface water and energy fluxes, influencing diurnal evolution of the planetary boundary layer (PBL). The mixing-diagram framework has proven useful in understanding the evolution of the heat and moisture budget within the convective boundary layer (CBL). We demonstrate that observations from the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains site provide all of the needed inputs needed for the mixing-diagram framework, allowing us to quantify the impact from the surface fluxes, advection, radiative heating, encroachment, and entrainment on the evolution of the CBL. Profiles of temperature and humidity retrieved from the ground-based infrared spectrometer [Atmospheric Emitted Radiance Interferometer (AERI)] are a critical component in this analysis. Large-eddy simulation results demonstrate that mean mixed-layer values derived are shown to be critical to close the energy and moisture budgets. A novel approach demonstrated here is the use of network of AERIs and Doppler lidars to quantify the advective fluxes of heat and moisture. The framework enables the estimation of the entrainment fluxes as a residual, providing a way to observe the entrainment fluxes without using multiple lidar systems. The high temporal resolution of the AERI observations enables the morning, midday, and afternoon evolution of the CBL to be quantified. This work provides a new way to use observations in this framework to evaluate weather and climate models. Significance Statement The energy and moisture budget of the planetary boundary layer (PBL) is influenced by multiple sources, and accurately representing this evolution in numerical models is critical for weather forecasts and climate predictions. The mixing-diagram approach, driven by profiling observations as illustrated here, provides a powerful way to quantify the contributions from each of these sources. In particular, the energy and moisture mixed into the PBL from above the PBL can be determined accurately from ground-based remote sensors using this approach.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

RF surface resistance tuning of superconducting niobium via thermal diffusion of native oxide

Recently, Nb superconducting radio frequency cavities vacuum heat treated between 300 and 400°C for a few hours have exhibited very high quality factors (~5 × 10 10 at 2.0 K). Secondary ion mass spectrometry measurements of O, N, and C show that this enhancement in RF surface conductivity is primarily associated with interstitial O alloying via dissolution and diffusion of the native oxide. We use a theory of oxide decomposition and O diffusion to quantify previously unknown parameters crucial in modeling this process. RF measurements of a vacuum heat-treated Nb superconducting radio frequency cavity confirm the minimized surface resistance (higher Q 0 ) previously expected only from 800 °C diffusive alloying with N.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Sodium tracer measurements of an expanded dense aluminum plasma from e-beam isochoric heating

Spatially and temporally resolved visible absorption spectroscopy is performed on sodium D-lines present as surface contaminants on an expanded dense aluminum plasma plume. An 80-ns FWHM, intense, relativistic electron beam deposits 5.4 J into a 100- μ m-thick Al foil, which isochorically heats and subsequently hydrodynamically expands the material through the warm dense matter state and into a classical-like plasma state, with a coupling parameter of approximately 0.2 and a degeneracy parameter of approximately 270. The Na contamination, carried along with the expanding plume, shows saturated absorption features in the dense Al continuum for λ > 450 nm. X-ray photoelectron spectroscopy and laser-induced breakdown spectroscopy confirm Na is a surface contaminant with an atomic concentration of ~ 0.1 % when interrogating identical foil samples. Additioanlly, a spectroscopic-quality radiation transport model is used to post-process 2D hydrodynamic simulations to interpret the plasma conditions based on the measured Na 3p-3s doublet line profiles. A sodium number density of 3 × 10 15 cm -3 best matches the experimental spectra, which originate from a dense surface plasma with n e = 3.0 ± 0.8 × 10 18 cm -3 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Continental-scale patterns of extracellular enzyme activity in the subsoil: an overlooked reservoir of microbial activity

Abstract Chemical stabilization of microbial-derived products such as extracellular enzymes (EE) onto mineral surfaces has gained attention as a possibly important mechanism leading to the persistence of soil organic carbon (SOC). While the controls on EE activities and their stabilization in the surface soil are reasonably well-understood, how these activities change with soil depth and possibly diverge from those at the soil surface due to distinct physical, chemical, and biotic conditions remains unclear. We assessed EE activity to a depth of 1 m (10 cm increments) in 19 soil profiles across the Critical Zone Observatory Network, which represents a wide range of climates, soil orders, and vegetation types. For all EEs, activities per mass of soil correlated positively with microbial biomass (MB) and SOC, and all three of these variables decreased logarithmically with depth ( p < 0.05). Across all sites, over half of the potential EE activities per mass soil consistently occurred below 20 cm for all measured EEs. Activities per unit MB or SOC were substantially higher at depth (soils below 20 cm accounted for 80% of whole-profile EE activity), suggesting an accumulation of stabilized (i.e. mineral sorbed) EEs in subsoil horizons. The pronounced enzyme stabilization in subsurface horizons was corroborated by mixed-effects models that showed a significant, positive relationship between clay concentration and MB-normalized EE activities in the subsoil. Furthermore, the negative relationships between soil C, N, and P and C-, N-, and P-acquiring EEs found in the surface soil decoupled below 20 cm, which could have also been caused by EE stabilization. This finding suggests that EEs may not reflect soil nutrient availabilities deeper in the soil profile. Taken together, our results suggest that deeper soil horizons hold a significant reservoir of EEs, and that the controls of subsoil EEs differ from their surface soil counterparts.

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

Edge ICRF simulations in 3D geometry: From MHD equilibrium to coupling determination

We present in this work a consistent numerical scheme that allows the computation of 3D magnetic fields a nd 3D density profiles and their usage in ion cyclotron range of frequencies (ICRF) coupling simulations. We first utilize the PARVMEC code to compute the 3D free-boundary plasma equilibrium in the ideal magnetohydrodynamic (MHD) approximation. Since the PARVMEC solution is only defined within the last closed flux surface (LCFS), the magnetic field domain is extended to the scrape-off layer (SOL) via the BMW code, which computes a divergence-free magnetic field solution arising from the external conductors’ vacuum field a nd t he PARVMEC flux surface currents. This magnetic reconstruction is then used in the EMC3-EIRENE transport code in order to compute 3D density profiles. In the last step, the RAPLICASOL code is utilized to compute the ICRF antenna S-matrices resulting from the 3D density profiles. We exemplify this scheme for the ASDEX Upgrade tokamak. A new implementation of a curved model for the ASDEX Upgrade ICRF 2-strap antenna in RAPLICASOL allows simulations in realistic geometry, without any coordinate transformations.

Lopez, Graham↗
DOE OSTI2020

Apparent dust size discrepancy in aerosol reanalysis in north African dust after long-range transport

North African dust reaches the southeastern United States every summer. Size-resolved dust mass measurements taken in Miami, Florida, indicate that more than one-half of the surface dust mass concentrations reside in particles with geometric diameters less than 2.1 µm, while vertical profiles of micropulse lidar depolarization ratios show dust reaching above 4 km during pronounced events. These observations are compared to the representation of dust in the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) aerosol reanalysis and closely related Goddard Earth Observing System model version 5 (GEOS-5) Forward Processing (FP) aerosol product, both of which assimilate satellite-derived aerosol optical depths using a similar protocol and inputs. These capture the day-to-day variability in aerosol optical depth well, in a comparison to an independent sun-photometer-derived aerosol optical depth dataset. Most of the modeled dust mass resides in diameters between 2 and 6 µm, in contrast to the measurements. Model-specified mass extinction efficiencies equate light extinction with approximately 3 times as much aerosol mass, in this size range, compared to the measured dust sizes. GEOS-5 FP surface-layer sea salt mass concentrations greatly exceed observed values, despite realistic winds and relative humidities. In combination, these observations help explain why, despite realistic total aerosol optical depths, (1) free-tropospheric model volume extinction coefficients are lower than those retrieved from the micro-pulse lidar, suggesting too-low model dust loadings in the free troposphere, and (2) model dust mass concentrations near the surface can be higher than those measured. The modeled vertical distribution of dust, when captured, is reasonable. Large, aspherical particles exceeding the modeled dust sizes are also occasionally present, but dust particles with diameters exceeding 10 µm contribute little to the measured total dust mass concentrations after such long-range transport. Remaining uncertainties warrant a further integrated assessment to confirm this study's interpretations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Heating effects on pyrogenic organic matter mass loss and temperature in a sand matrix in a lab setting from 2022

This dataset contains data associated with the preprint “A laboratory method for simulating the effects of subsequent fires on pyrogenic organic matter at different exposure depths in a sand matrix” (https://doi.org/10.1101/2024.07.30.605814), which has been accepted for publication in the International Journal of Wildland Fire. The dataset was generated to investigate the effects of subsequent fires on preexisting pyrogenic organic matter (PyOM) at different soil depths. Specifically, we designed a laboratory method to simulate soil heating from above, applying realistic heat flux profiles (High, Low, Control) derived from jack pine (Pinus banksiana Lamb.) log burns. Using a cone calorimeter, these profiles were applied to buried jack pine PyOM (produced at 350°C) in a sand matrix pyrocosm to examine PyOM mass loss at different depths (surface, 1 cm, 5 cm). Each treatment (e.g., High heat flux at 1 cm depth) was applied once to five replicates in a pyrocosm.The dataset includes mass loss fraction of PyOM and time-series temperature profiles (measured with type K thermocouples). Additionally, calibration datasets for high and low heat flux profiles are provided, including programmed heat flux profile data and three sets of calibration results. All datasets are in .csv format and were analyzed in R using code found at https://github.com/MengmengLuo/Reburning-pyrogenic-organic-matter-A-laboratory-method-for-dosing-dynamic-heat-fluxes-from-above. This dataset supports research on PyOM transformations during subsequent fires, and methodological advancements in simulating soil heating. It can be used to explore fire-induced PyOM alterations across different fire intensities and soil depths.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dark matter cores and cusps in spiral galaxies and their explanations

We compare proposed solutions to the core vs cusp issue of spiral galaxies, which has also been framed as a diversity problem, and demonstrate that the cuspiness of dark matter halos is correlated with the stellar surface brightness. We compare the rotation curve fits to the SPARC sample from a self-interacting dark matter (SIDM) model, which self-consistently includes the impact of baryons on the halo profile, and hydrodynamical N-body simulations with cold dark matter (CDM). The SIDM model predicts a strong correlation between the core size and the stellar surface density, and it provides the best global fit to the data. The CDM simulations without strong baryonic feedback effects fail to explain the large dark matter cores seen in low surface brightness galaxies. On the other hand, with strong feedback, CDM simulations do not produce galaxy analogs with high stellar and dark matter densities, and therefore they have trouble in explaining the rotation curves of high surface brightness galaxies. This implies that current feedback implementations need to be modified. We also explicitly show how the concentration-mass and stellar-to-halo mass relations together lead to a radial acceleration relation (RAR) in an averaged sense, and reiterate the point that the RAR does not capture the diversity of galaxy rotation curves in the inner regions. These results make a strong case for SIDM as the explanation for the cores and cusps of field galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Influence of plateau, slope, and valley on soil hydrology during the dry season in a Central Amazon old‐growth forest

Soil moisture regulates plant water supply and drought sensitivity in tropical forests, yet its vertical and topographic variation remains poorly characterized. We combined high-frequency time-domain reflectometry measurements from 5 to 100 cm across plateau, slope, and valley landforms at the Zona Florestal 2 research site north of Manaus, Central Amazonia, to quantify how soil moisture memory, timing of responses to rainfall, dry-down rates (τ), and soil–water depletion vary across these contrasting landforms. Landform-specific soil moisture calibration curves ensured accurate volumetric water content estimates in these highly weathered soils. During the 2023 dry-to-wet transition (August–November), soil moisture memory showed strong topographic contrasts, with valley profiles increasing from ∼47 h at 5 cm to ∼154 h at 100 cm, while plateaus exhibited higher near-surface persistence (∼124 h at 5 cm) but weaker memory at depth. Dry-down behavior reinforced these differences as valley soils exhibited τ values exceeding ∼200 h, more than double the characteristic τ of plateau soils (∼90 h). Rainfall–soil moisture correlations indicated immediate responses at shallow depths in valleys and progressively longer lags with depth on plateaus and slopes. These hydrologic patterns were mirrored in depletion profiles, which declined sharply below 30 cm on plateaus but remained high and sustained throughout the upper meter in slopes and valleys. Together, these findings provide the first depth-resolved field measurements of soil moisture memory, rainfall coupling, dry-down constants, and depletion dynamics across major upland landforms in Central Amazonia and offer clear observational benchmarks for improving land-surface and ecosystem model representations of soil–water processes.

Hillslope↗
DOE OSTI2025

Tailoring the Physicochemical Properties of Nb Thin Films via Surface Engineering Methods

The modification of surface oxide layers formed on niobium (Nb) thin films via chemical mechanical planarization (CMP) and accelerated neutral atom beam (ANAB) processing provides a promising route toward tailoring their emergent properties and performance when used as superconducting qubits. Here, in this study, we show that CMP- and ANAB-formed Nb oxides are significantly thinner and smoother than the native oxide, as revealed by transmission electron microscopy (TEM) and atomic force microscopy. Scanning TEM and energy-dispersive X-ray spectroscopy along with X-ray photoelectron spectroscopy identified an oxidation gradient within the native and surface-engineered oxides. The topside layer is dominated by Nb 5+ (Nb 2 O 5 ), with various Nb suboxides present closer to the oxide/metal interface. Time-of-flight secondary ion mass spectrometry (ToF-SIMS) depth profiling confirmed the presence of an oxygen content gradient and demonstrated the enhanced resistance of the CMP- and ANAB-formed oxides to oxygen surface exchange and subsequent diffusion via 18 O 2 isotopic labeling experiments. ToF-SIMS also identified an interfacial layer containing trapped hydrogen (H)-containing species at the Nb oxide/metal interface. In situ ToF-SIMS and TEM revealed migration of the H/OH interfacial layer coinciding with decomposition of the surface oxide. Furthermore, our density functional theory calculations indicated that both H from moisture present in ambient air and bulk H in Nb films tend to segregate at the interface. These findings underscore the importance of understanding surface oxidation mechanisms, hydrogen incorporation, and their impact on the designed functionalities of Nb-based devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluating wind profiles in a numerical weather prediction model with Doppler lidar

We use Doppler lidar wind profiles from six locations around the globe to evaluate the wind profile forecasts in the boundary layer generated by the operational global Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts (ECMWF). The six locations selected cover a variety of surfaces with different characteristics (rural, marine, mountainous urban, and coastal urban). We first validated the Doppler lidar observations at four locations by comparison with co-located radiosonde profiles to ensure that the Doppler lidar observations were of sufficient quality. The two observation types agree well, with the mean absolute error (MAE) in wind speed almost always less than 1 m s –1 . Large deviations in the wind direction were usually only seen for low wind speeds and are due to the wind direction uncertainty increasing rapidly as the wind speed tends to zero. Time–height composites of the wind evaluation with 1 h resolution were generated, and evaluation of the model winds showed that the IFS model performs best over marine and coastal locations, where the mean absolute wind vector error was usually less than 3 m s –1 at all heights within the boundary layer. Larger errors were seen in locations where the surface was more complex, especially in the wind direction. For example, in Granada, which is near a high mountain range, the IFS model failed to capture a commonly occurring mountain breeze, which is highly dependent on the sub-grid-size terrain features that are not resolved by the model. The uncertainty in the wind forecasts increased with forecast lead time, but no increase in the bias was seen. At one location, we conditionally performed the wind evaluation based on the presence or absence of a low-level jet diagnosed from the Doppler lidar observations. The model was able to reproduce the presence of the low-level jet, but the wind speed maximum was about 2 m s –1 lower than observed. This is attributed to the effective vertical resolution of the model being too coarse to create the strong gradients in wind speed observed. Our results show that Doppler lidar is a suitable instrument for evaluating the boundary layer wind profiles in atmospheric models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The ability of Salmonella enterica subsp. enterica strains to form biofilms on abiotic surfaces and their susceptibility to selected essential oil components

Abstract The ability of Salmonella enterica subsp. enterica to persist and form biofilms on different surfaces can constitute a source of food contamination, being an issue of global concern. The objective of this study was to understand the biofilm formation profile of 14 S. enterica strains among different serovars and sources and to evaluate the ability of essential oil (EO) components (carveol, citronellol, and citronellal) to disinfect the biofilms formed on stainless steel and polypropylene surfaces. All the strains were able to form biofilms with counts between 5.34 to 6.78 log CFU cm−2. Then, the minimum inhibitory concentration (MIC) and minimum bactericidal concentration (MBC) of EO components were evaluated on two selected strains. All compounds inhibited the growth of Salmonella Typhimurium (strain 1; MIC = 800–1000 µg ml−1) and Salmonella Enteritidis (strain 5; MIC = 400–1000 µg ml−1) and only carveol showed bactericidal activity against strains 1 and 5 (MBC = 1200 µg ml−1). Biofilms were exposed to the EO components at 10 × MIC for 30 min and polypropylene surfaces were more difficult to disinfect showing reductions between 0.9 and <1.2 log CFU cm−2. In general, the S. enterica biofilms demonstrated a significant tolerance to disinfection, demonstrating their high degree of recalcitrance on food processing surfaces.

Pérez-Lavalle, Liliana (ORCID:0000000344825031)↗
DOE OSTI2020

Retention of Long-period Gas Giant Planets: Type II Migration Revisited

During their formation, emerging protoplanets tidally interact with their natal disks. Proto–gas giant planets, with Hill radii larger than the disk thickness, open gaps and quench gas flow in the vicinity of their orbits. It is usually assumed that their type II migration is coupled to the viscous evolution of the disk. Although this hypothesis provides an explanation for the origin of close-in planets, it also encounters a predicament on the retention of long-period orbits for most gas giant planets. Moreover, numerical simulations indicate that the planets' migrations are not solely determined by the viscous diffusion of their natal disk. In this work, we carry out a series of hydrodynamic simulations combined with analytic studies to examine the transition between different paradigms of type II migration. We find a range of planetary mass for which gas continues to flow through a severely depleted gap so that the surface density distribution in the disk region beyond the gap is maintained in a quasi-steady state. The associated gap profile modifies the location of corotation and Lindblad resonances. In the proximity of the planet's orbit, high-order Lindblad and corotation torque are weakened by the gas depletion in the gap, while low-order Lindblad torques near the gap walls preserve their magnitude. Consequently, the intrinsic surface density distribution of the disk delicately determines both the pace and direction of the planets' type II migration. We show that this effect might stall the inward migration of giant planets and preserve them in disk regions where the surface density is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Algebraic expressions for estimating the impact depths of a surface barrier over a homogeneous soil

Engineered surface barriers are used to isolate subsurface contaminants for effective long-term containment of municipal solid waste, other nonhazardous solid and liquid waste, hazardous and toxic wastes, and radioactive waste. The impact depths of a surface barrier are affected by the pre-barrier recharge rate and the properties of the soil beneath the barrier. In this paper, the pore-size-specific (PSS) water velocity is defined and an algebraic expression of PSS velocity is derived based on the stream tube concept and the Brooks and Corey hydraulic retention model. Algebraic expressions are developed to estimate drainage velocities and barrier impact depths after the emplacement of a surface barrier. Four impact depth terms are used to convey the protective effect: drainage front, average drainage, the location with 50% impact, and drainage tail. The drainage front depth is the deepest point at which the barrier has a detectable impact at a specific time (also called the near zero-impact depth). The average-impact depth is the depth at which average drainage occurs. At the 50% impact depth, the water flux rate is reduced by half because of the surface barrier. Lastly, the drainage tail depth (also called the full-impact depth) is the deepest depth at which the water conditions above it are in equilibrium with the barrier. The algebraic expressions show that the average-impact depth is no more than 1/3 of near zero-impact depth, while the 50% impact depth is slightly larger than 1/2 of the near zero-impact depth. The full-impact depth, depending on the final recharge rate from the surface barrier, is usually much smaller than the other impact depths. These differences lead to a very large transition zone beneath a surface barrier. Numerical simulations were conducted to replicate the same conditions. The numerical results corroborated the analytical models by predicting very similar water content profiles and near zero-, average-, 50%, and full-impact depths. The algebraic expressions provided in this paper are useful for quickly identifying sites where the depth of the existing contaminants could be beyond the protection of a surface barrier.

54 ENVIRONMENTAL SCIENCES↗