SEARCH · Engineering Papers

Results for “ADVECTION”

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.

SourceResearch connectionFrom yearThrough year

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

At least 811 records · Page 45

All sources NASA NTRS DOE OSTI

NASA NTRS2016

On the Synergy Between Numerics and Subgrid Scale Modeling in LES of Stratified Flows: Grid Convergence of a Stratocumulus-Topped Boundary Layer

The effectiveness of a linear upwinding scalar advection scheme to suppress numerical dispersion errors near sharp inversions in large-eddy simulations of a nocturnal stratocumulus-topped boundary layer is assessed. Linear upwinding is a trade-off between non-dissipative and non-linear positive definite advection schemes. It is shown that linear upwinding does not negatively impact the model's grid convergence properties and a sharp inversion free of numerical artifacts is maintained. Even though mean profiles and turbulence fluxes show good grid convergence characteristics the liquid water amount varies significantly with grid resolution. The entrainment rate is identical for all resolutions and independent of the liquid water amount. For the present stratocumulus case, the impact of cloud-top radiative cooling is negligible and turbulence is largely driven by convection emanating from the surface.

Matheou, Georgios↗

NASA NTRS2019

Possible Impacts of Snow Darkening Effects on the Hydrological Cycle over Western Eurasia and East Asia

In this paper, we investigated the possible impact of snow darkening effect (SDE) by light-absorbing aerosols on the regional changes of the hydrological cycle over Eurasia using the NASA GEOS-5 Model with aerosol tracers and a state-of-the-art snow darkening module, the Goddard SnoW Impurity Module (GOSWIM) for the land surface. Two sets of ten-member ensemble experiments for 10 years were carried out forced by prescribed sea surface temperature (2002–2011) with different atmospheric initial conditions, with and without SDE, respectively. Results show that SDE can exert a significant regional influence in partitioning the contributions of evaporative and advective processes on the hydrological cycle, during spring and summer season. Over western Eurasia, SDE-induced rainfall increase during early spring can be largely explained by the increased evaporation from snowmelt. Rainfall, however, decreases in early summer due to the reduced evaporation as well as moisture divergence and atmospheric subsidence associated with the development of an anomalous mid- to upper-tropospheric anticyclonic circulation. On the other hand, in the East Asian monsoon region, moisture advection from the adjacent ocean is a main contributor to rainfall increase in the melting season. A warmer land-surface caused by earlier snowmelt and subsequent drying further increases moisture transport and convergence significantly enhancing rainfall over the region. Our findings suggest that the SDE may play an important role in leading to hotter and drier summers over western Eurasia, through coupled land-atmosphere interaction, while enhancing East Asian summer monsoonal precipitation via enhanced land-ocean thermal contrast and moisture transport due to the SDE-induced warmer Eurasian continent.

snow darkening effect↗

NASA NTRS2021

Large-Eddy Simulations of Marine Boundary-Layer Clouds Associated With Cold Air Outbreaks During the Activate Campaign – Part 1: Case Setup and Sensitivities to Large-Scale Forcings

Large-eddy simulation (LES) is able to capture key boundary-layer (BL) turbulence and cloud processes. Yet, large-scale forcing and surface turbulent fluxes of sensible and latent heat are often poorly prescribed for LES simulations. We derive these quantities from measurements and reanalysis obtained for two cold air outbreak (CAO) events during Phase I of the Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) in February-March 2020. We study the two contrasting CAO cases by performing LES and test the sensitivity of BL structure and clouds to large-scale forcings and turbulent heat fluxes. Profiles of atmospheric state and large-scale divergence and surface turbulent heat fluxes obtained from the reanalysis data ERA5 agree reasonablywell with those derived fromACTIVATE field measurements for both cases at the sampling time and location. Therefore, we adopt the time evolving heat fluxes, wind and advective tendencies profiles from ERA5 reanalysis data to drive the LES.We find that large-scale thermodynamic advective tendencies and wind relaxations are important for the LES to capture the evolving observed BL meteorological states characterized by the hourly ERA5 reanalysis data and validated by the observations. We show that the divergence (or vertical velocity) is important in regulating the BL growth driven by surface heat fluxes in LES simulations. The evolution of liquid water path is largely affected by the evolution of surface heat fluxes. The liquid water path imulated in LES agrees reasonably well with the ACTIVATE measurements. This study paves the path to investigate aerosol-cloud-meteorology interactions using LES informed and evaluated by ACTIVATE field measurements.

Large-eddy simulation↗

NASA NTRS2021

Quasi-two-day wave modulation of carbon dioxide in the mesosphere and lower thermosphere

In this paper, we report the quasi-two-day wave (QTDW) in daytime only SABER CO2 in January 2003. The QTDW in CO2 of ∼7 ppmv or 3–4% shows double peak structures at mid and low latitudes of both hemispheres. We compare the SABER measured QTDW disturbance in CO2 and temperature to those in the eCMAM (the extended Canadian Middle Atmosphere Model). Good agreement of QTDW CO2 spatial and temporal variations is seen, especially at southern mid-latitudes. eCMAM outputs are further analyzed to show that the QTDW in CO2 mixing ratio is mainly driven by both vertical transport and meridional advection. Because waves in vertical transport can be derived by temperature disturbance via a so-called adiabatic displacement approach, we can estimate the QTDW in CO2 mixing ratio from the SABER QTDW in temperature assuming meridional advection is weak at low and mid-latitudes.

planetary wave; carbon dioxide↗

NASA NTRSDate not supplied

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

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

weather↗

NASA NTRSDate not supplied

A Generalized, Compactly-Supported Correlation Function for Data Assimilation Applications

Correlation functions play an essential role in modern data assimilation, where they are used to model covariances given a set of tunable parameters or applied as tapering functions to localize covariances in ensemble-based schemes. One of the most widely-used correlation functions in data assimilation is the Gaspari and Cohn (1999) piecewise-rational, compactly-supported parametric correlation function (hereafter referred to as GC99). The GC99 correlation function is useful due to its tunable cut-off parameter c and Gaussian-like shape achieved when the parameter a is set to one-half. These properties are attractive for tapering functions in data assimilation applications. However, the GC99 correlation function is homogeneous over Euclidean 3-space and isotropic when restricted to the sphere, properties that may be less than ideal for some geophysical applications. GC99 is also compactly-supported on a sphere of fixed radius, which requires tuning of the cut-off parameter c that can depend on the specific application. This work presents a generalization of the GC99 correlation function that allows the cut-off parameter c and shape parameter a to vary over space to gain more flexibility in shape while maintaining its compact support property. The function, which we call the Generalized Gaspari Cohn (GenGC) correlation function, introduces inhomogeneity in Euclidean 3-space and anisotropy when restricted to the sphere by allowing both parameters c and a to vary, as functions, over the spatial domain. The GC99 correlation function is a special case of GenGC where the functions c and a are held constant, as fixed parameters rather than functions. The GenGC correlation function also generalizes the follow-on to the work of Gaspari and Cohn (1999) presented in Gaspari et al. (2006), which allowed a to vary while keeping c fixed. We illustrate through simple one- and two-dimensional examples the variety of inhomogeneous and anisotropic correlation functions GenGC can produce by varying c and a over space, and suggest applications where they may be useful in data assimilation, such as covariance modeling or localization. In particular, we describe how the GenGC correlation function can be used to construct covariances using correlation length and variance fields derived from dynamics. For example, the correlation length field for advective dynamics is governed by a partial differential equation (PDE) in N spatial dimensions, where N is the number of space dimensions of the state. Correlation length fields can be determined from this PDE and used with GenGC to construct the corresponding correlations. We can then approximate the full covariance by rescaling by the variance, which also satisfies a PDE in N spatial dimensions for advective dynamics. Thus we can approximate the full covariance without solving the covariance PDE, which is in 2N spatial dimensions, by solving just two PDEs each in only N spatial dimensions. This approach to evolving the correlation length and variance fields, then reconstructing the correlations using GenGC, is suggested as an alternative to current methods of covariance modeling in data assimilation algorithms.

GC99↗

NASA NTRS2022

Winds at the Mars 2020 Landing Site. 2. Wind Variability and Turbulence

Wind speeds measured by the Mars 2020 Perseverance rover in Jezero crater were fitted as a Weibull distribution. InSight wind data acquired in Elysium Planitia were also used to contextualize observations. Jezero winds were found to be much calmer on average than in previous landing sites, despite the intense aeolian activity observed. However, a great influence of turbulence and wave activity was observed in the wind speed variations, thus driving the probability of reaching the highest wind speeds at Jezero, instead of sustained winds driven by local, regional, or large-scale circulation. The power spectral density of wind speed fluctuations follows a power-law, whose slope deviates depending on the time of day from that predicted considering homogeneous and isotropic turbulence. Daytime wave activity is related to convection cells and smaller eddies in the boundary layer, advected over the crater. The signature of convection cells was also found during dust storm conditions, when prevailing winds were consistent with a tidal drive. Nighttime fluctuations were also intense, suggesting strong mechanical turbulence. Convective vortices were usually involved in rapid wind fluctuations and extreme winds, with variations peaking at 9.2 times the background winds. Transient high wind events by vortex-passages, turbulence, and wave activity could be driving aeolian activity at Jezero. We report the detection of a strong dust cloud of 0.75–1.5 km in length passing over the rover. The observed aeolian activity had major implications for instrumentation, with the wind sensor suffering damage throughout the mission, probably due to flying debris advected by winds.

D. Viúdez-moreiras↗

NASA NTRS2022

Estimating the Migrating Diurnal Tide Component of Mesospheric Water Vapor

This work presents a method for estimating the migrating diurnal tide (DW1) component of mesospheric H2O from observations of the temperature tide and zonal-mean H2O made by the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument. This work first shows that a 2D least-squares fit on SABER H2O yields an erroneous DW1 due to retrieval algorithm problems. Consequently, the adiabatic displacement method is used. Applying the method to SABER observations and comparing it with Microwave Limb Sounder (MLS) H2O observations reveals that the method reproduces the MLS-observed H2O DW1 component best in March and June over low-latitudes in the altitude range between 65 and 75 km where errors range between +5% and +20%. Applying the method for simulations using the Specified Dynamics—Whole Atmosphere Community Climate Model with Ionosphere/Thermosphere eXtension (SD-WACCM-X) indicates that in the model, it is best only in March and not in June. Model simulations further showed that in March, theoretical errors due to the vertical advection assumption and aliasing are best at less than +13% for the low-latitudes and altitudes between 65 and 75 km. These results, therefore, lead to the conclusion that despite the errors in the vertical advection assumption and aliasing, the adiabatic displacement method performs best in estimating mesospheric H2O DW1 with SABER observations during both March equinox and June solstice over the low latitudes between 65 and 75 km. Results also show that SD-WACCM-X satisfactorily simulates the H2O DW1 magnitude and mechanisms only in the March equinox period.

h20↗

NASA NTRS2023

Spatial and Temporal Characterization of Cyanobacteria Blooms in the Mississippi Sound and Their Relationship to the Bonnet Carré Spillway Openings

During the spring and summer of 2019, an unprecedented cyanobacterial harmful algal bloom (cyanoHAB) was responsible for beach advisories on 25 beaches along the Mississippi Sound for over 3 months. Due to the preceding heavy rainfall and flooding within the Mississippi River watershed, for the first time in history, the Bonnet Carré Spillway (BCS) opened twice in one year during 2019. The coastal cyanoHAB coincided with the second BCS opening. The main objectives of this study were: (1) to investigate the potential for using the National Aeronautics and Space Administration (NASA) ocean color standard Cyanobacteria Index (CIcyano) algorithm to characterize the spatial and temporal extent of the 2019 cyanoHAB; (2) to couple the CIcyano data with river discharge, salinity, and modeled-wind data to study the conditions leading to the cyanoHAB and factors aiding the advection and persistence of the bloom within the Mississippi Sound, including a possible relationship to the BCS; (3) to further investigate the relationship with the BCS by repeating the methods using data from 2018, which was a year when the BCS was opened but no evidence of cyanoHABs was reported along the Mississippi coast. Weekly means and monthly frequency CIcyano images, river discharge, salinity, and modeled-wind data from February to September of 2018 and 2019 were analyzed, which coincide with three BCS openings. In March 2018, a cyanobacteria bloom was observed within Lake Pontchartrain coinciding with the BCS opening; however, the month-long bloom was contained to the lake. Two distinct cyanoHABs were observed in 2019 and both blooms were advected into the Mississippi Sound, and likely contributed to the 3-month-long beach water advisories of 2019 along the Mississippi coastline. From March to mid-July 2019, salinity at stations within the Mississippi Sound was consistently near zero indicating high levels of freshwater. During that time, winds were predominantly northwestward, preventing the BCS waters from flushing into the Mississippi Shelf and resulting in BCS waters remaining longer within the estuarine lakes and Mississippi Sound. Although the BCS had an undeniable impact on the presence of the coastal cyanoHAB of 2019, other variables including wind direction, water flow, mixing, and persistence of freshwater within the Sound can determine the intensity and extent of the cyanoHABs. Coupling in situ phytoplankton information from freshwater water bodies to the marine continuum along with water flow, wind data, and satellite imagery could help identify cyanoHABs at early stages and forecast their trajectory and potential impacts on coastal areas.

Bonnet Carré Spillway↗

NASA NTRS2023

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

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 (SGP) 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 (called the Atmospheric Emitted Radiance Interferometer, or 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. Finally, the high temporal resolution of the AERI observations enable 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.

Boundary layer↗

NASA NTRS2023

Distributions of Mesoscale Periodic Structures in the Elemental and Ionic Composition of the Solar Wind

Multiple statistical and event studies based on in situ observations have shown that the solar wind contains mesoscale (∼ 100 – 10,000 Mm) periodic structures in the proton number density. Remote observations of such structures and event studies of concurrent variations in composition have demonstrated that they can form in the solar atmosphere and be preserved while advecting outwards through the heliosphere. Viall, Kepko, and Spence (J. Geophys. Res. (Space Phys.) 113, A07101, 2008; J. Geophys. Res. (Space Phys.) 114, A01201, 2009) and Kepko, Viall, and Wolfinger (J. Geophys. Res. (Space Phys.) 125(8), e28037, 2020) have reported that periodic proton density structures preferentially occur at specific radial length scales and have published their distributions from Wind measurements near L1. Here, we conduct a statistical study of 14 years (1998 – 2011) of 12-minute composition data measured by the Solar Wind Ion Composition Spectrometer instrument aboard the Advanced Composition Explorer spacecraft. We found that the elemental and ionic composition also contain statistically significant mesoscale periodic structures and, for the first time, present occurrence distributions for elemental abundance ratios with low, intermediate, and high first ionization potentials as well as for key solar wind charge states. These distributions set important constraints on solar wind formation in general and the formation of periodic mesoscale solar wind structures specifically, as the elemental and ionic composition are known to be determined at the Sun and to not evolve during advection.

Solar wind↗

DOE OSTI2015

Studying methane migration mechanisms at Walker Ridge, Gulf of Mexico, via 3D methane hydrate reservoir modeling

We have developed a 3D methane hydrate reservoir simulator to model marine methane hydrate systems. Our simulator couples highly nonlinear heat and mass transport equations and includes heterogeneous sedimentation, in-situ microbial methanogenesis, the influence of pore size contrast on solubility gradients, and the impact of salt exclusion from the hydrate phase on dissolved methane equilibrium in pore water. Using environmental parameters from Walker Ridge in the Gulf of Mexico, we first simulate hydrate formation in and around a thin, dipping, planar sand stratum surrounded by clay lithology as it is buried to 295mbsf. We find that with sufficient methane being supplied by organic methanogenesis in the clays, a 200x pore size contrast between clays and sands allows for a strong enough concentration gradient to significantly drop the concentration of methane hydrate in clays immediately surrounding a thin sand layer, a phenomenon that is observed in well log data. Building upon previous work, our simulations account for the increase in sand-clay solubility contrast with depth from about 1.6% near the top of the sediment column to 8.6% at depth, which leads to a progressive strengthening of the diffusive flux of methane with time. By including an exponentially decaying organic methanogenesis input to the clay lithology with depth, we see a decrease in the aqueous methane supplied to the clays surrounding the sand layer with time, which works to further enhance the contrast in hydrate saturation between the sand and surrounding clays. Significant diffusive methane transport is observed in a clay interval of about 11m above the sand layer and about 4m below it, which matches well log observations. The clay-sand pore size contrast alone is not enough to completely eliminate hydrate (as observed in logs), because the diffusive flux of aqueous methane due to a contrast in pore size occurs slower than the rate at which methane is supplied via organic methanogenesis. Therefore, it is likely that additional mechanisms are at play, notably bound water activity reduction in clays. Three-dimensionality allows for inclusion of lithologic heterogeneities, which focus fluid flow and subsequently allow for heterogeneity in the methane migration mechanisms that dominate in marine sediments at a local scale. Incorporating recently acquired 3D seismic data from Walker Ridge to inform the lithologic structure of our modeled reservoir, we show that even with deep adjective sourcing of methane along highly permeable pathways, local hydrate accumulations can be sourced either by diffusive or advective methane flux; advectively-sourced hydrates accumulate evenly in highly permeable strata, while diffusively-sourced hydrates are characterized by thin strata-bound intervals with high clay-sand pore size contrasts.

Nole, Michael↗

DOE OSTI2016

Magnetic method for stimulating transport in fluids

A method for producing mass and heat transport in fluids, wherein the method does not rely on conventional convection, that is, it does not require gravity, a thermal gradient, or a magnetic field gradient. This method gives rise to a unique class of vigorous, field-controllable flow patterns termed advection lattices. The advection lattices can be used to transport heat and/or mass in any desired direction using only magnetic fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Nonsteady Load Responses to Daytime Atmospheric Turbulence Eddies on the DOE 1.5 MW Wind Turbine at NREL

Field data collected from the NREL/GE 1.5MW wind turbine (WT) and met tower (MetT) at the NREL Wind Technology Center near Boulder, CO June-October 2018 were analyzed to quantify the impacts of turbulence eddies on the load responses measured from sensors on the main shaft, blade and tower. The passage of individual mountain-generated eddies from the met tower to the WT were critically determined by correlating the optimal time shifts in signal between MetT and nacelle anemometers with mean advection time. Loading responses from mountain eddy passage were compared with atmospheric eddies from the north/south, unimpeded by the mountains, and found to be similar. Whereas time variations in torque were highly correlated with time changes in horizontal eddy velocity, the out-of-plane bending moments on the main shaft (directly forcing the main bearing) were uncorrelated with horizontal eddy velocity. This result is consistent with a previous LES study indicating that the main bearing is forced by asymmetrical interactions between the WT rotor and turbulence eddies, while power fluctuations respond primarily to advective eddy velocity. Surprisingly, the nacelle anemometer produced statistics very similar to the MetT.

ENGINEERING,WIND ENERGY↗

DOE OSTI2023

Nonsteady Load Responses to Mountain-Generated Turbulence Eddies on the DOE 1.5 MW Wind Turbine at the National Wind Technology Center

Field data collected from the NREL/GE 1.5MW wind turbine and met tower at the NREL Wind Technology Center near Boulder, Colorado from June-October 2018 were analyzed to quantify the impacts of turbulence eddies on the load responses measured from sensors on the main shaft, blade and tower. The passage of individual mountain-generated eddies from the met tower to the wind turbine were critically determined by correlating the optimal time shifts in signal between met tower and nacelle anemometers with mean advection time. Loading responses from mountain eddy passage were compared with atmospheric eddies from the north/south, unimpeded by the mountains, and found to be similar. Whereas time variations in torque were highly correlated with time changes in horizontal eddy velocity, the out-of-plane bending moments on the main shaft (directly forcing the main bearing) were uncorrelated with horizontal eddy velocity. This result is consistent with a previous LES study indicating that the main bearing is forced by asymmetrical interactions between the rotor and turbulence eddies, while power fluctuations respond primarily to advective eddy velocity. Surprisingly, the nacelle anemometer produced statistics very similar to the met tower.

ENGINEERING,WIND ENERGY↗

DOE OSTI2020

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

The Influence of Aqueous Dispersivity on Transport of Radiological Material at LANL MDA G - 20450

Aqueous dispersivity is a phenomenon by which a solute moving in groundwater mixes with water that does not contain the solute and becomes diluted. Solutes take paths of varying lengths through pore networks or experience slower velocities due to interactions with pore walls. Aqueous dispersivity leads to breakthrough curves with earlier arrival than expected in purely advective systems and persistent 'tails,' or late-arriving, dilute concentrations of the solute. Dispersivity can be longitudinal (in the direction of flow) or transverse (orthogonal to flow). LANL has disposed of radioactive material at Material Disposal Area (MDA) G since the late 1950's. Most of the disposed materials at MDA G are low-level waste (LLW), but other types of radioactive wastes are disposed or stored at the site. For example, there is a collection of remote-handled transuranic waste (RH-TRU) waste in 33 shafts, which were intended to be stored retrievably. Aqueous-phase radionuclides included in models of contaminant transport at MDA G are subject to dispersive processes in porous media of the vadose zone and saturated zone. To model dispersivity at MDA G probabilistically, we created site-specific distributions of aqueous dispersivity. First, we performed a comprehensive review of peer-reviewed publications, reports from other performance assessments of MDA G, and documentation from sites that were considered similar based on expert opinion, such as Yucca Mountain. Data were included from a wide range of sources, so that distributions developed based on the data reflect the current state of knowledge in properties of the subsurface environment. Data were grouped based on value type (e.g., determined via experiment or used in a different model), rock type, zone type (saturated or unsaturated), and site, and were weighted in the distributions according to their relevance to the site and zone type. We combined Monte Carlo simulation and bootstrapping techniques to produce approximate distributions of mean aqueous dispersivity. Distributions of dispersivity for saturated and unsaturated zone transport are differentiated by data weight and the length scales over which transport occurs. We find that underestimating or disregarding transverse dispersivity may lead to overstating simulated well capture in the saturated zone at MDA G. Neglecting transverse dispersivity is a common assumption due in part to use of one-dimensional transport models that do not consider solute spreading orthogonal to flow. Also, we find that variability in longitudinal dispersivity is associated with substantial travel time alterations in the MDA G unsaturated zone and in simulated contaminant breakthrough times at a downstream observation plane in the MDA G saturated zone. Finally, at large scales (>1000 m), a linear trend between the dispersivity coefficient and length over which transport occurs may overestimate the dispersivity coefficient. Future modeling work with application to solute spreading could include interacting advective and diffusive zones in the highly fractured basalt below MDA G. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

A remapping scheme for history-dependent material state variables

Here, we present an interpolation scheme for use with material models involving history-dependent state variables. The interpolation scheme is based on the solution of a constrained optimization problem, to identify the smoothest monotone curve which produces an admissible solution. Unlike classical schemes used for transport equations, the proposed methodology accounts for the interdependence of the history variables and produces admissible interpolations of the material's state variables. After introducing a general solution scheme for the problem, a simplified and more computationally efficient method is derived. Finally, as a validation problem for the proposed methods, we consider the advection of internal state variables belonging to material models with highly nonlinear yield surfaces.

42 ENGINEERING↗