Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A two-fluid single-column model of turbulent shallow convection, Part III: Results and parameter sensitivity

Two-fluid modelling has recently emerged as a promising approach to representing cumulus convection in weather and climate models. This study applies the two-fluid model described in Part II Thuburn et al., 2021b to a shallow cumulus convection case study over land (ARM). Large eddy simulation data is used to tune the majority of the closures that determine the properties of entrained and detrained air. The two-fluid model is generally able to reproduce the profiles of the mean and turbulent quantities over all stages of the diurnal cycle. As such, the initiation of shallow convection and the evolution of the cloud layer are well captured. The robustness of the two-fluid model is further verified using a steady-state test case (BOMEX), in which the cloud properties are also well modelled.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The role of thermal stratification on the co‐spectral properties of momentum transport above an Amazonian forest

The influence of thermal stratification on the turbulent kinetic energy balance has been widely studied; however, its influence on the turbulent stress remains less explored in the presence of tall vegetated canopies and less ideal micrometeorological conditions. Here, the impact of thermal stratification on turbulent momentum flux is considered in the roughness sublayer (RSL) and the atmospheric surface layer (ASL) using the Amazon Tall Tower Observatory (ATTO) in Brazil. A scalewise co‐spectral budget (CSB) model is developed using standard closure schemes for the pressure–velocity decorrelation. The CSB revealed that the co‐spectrum ${F}_{wu}\left({k}_x\right)$ between longitudinal (u') and vertical (w') velocity fluctuations is impacted by the energy spectrum of the vertical velocity ${E}_{ww}\left({k}_x\right)$ and the much less studied longitudinal heat‐flux co‐spectrum ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$, where ${\theta}_{\mathrm{v}}^{\prime }$ are temperature fluctuations and ${k}_x$ is the longitudinal wavenumber. Under stable, very stable, and dynamic–convective conditions, the scaling exponent ${F}_{wu}\left({k}_x\right)$ in for the inertial subrange (ISR) scales is dominated by ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$ instead of ${E}_{ww}\left({k}_x\right)$. A near ${k}_x^{-7/3}$scaling in ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$ robust to large variations in thermal stratification is found, whereas the Kolmogorov ISR scaling for ${E}_{ww}\left({k}_x\right)\sim {k}_x^{-5/3}$ is not found. The scale‐dependent decorrelation time between u' and w' is dominated by ${\epsilon}^{-1/3}{k}_x^{-2/3}$ in the ISR, but is nearly constant for eddies larger than the vertical velocity integral scale, regardless of stability. Implications of these findings for generalized stability correction functions that are based on the turbulent stress budget instead of the turbulent kinetic energy budget are discussed.

canopy turbulence↗
DOE OSTI2023

High-throughput mass spectrometry analysis using immediate drop-on-demand technology coupled with an open port sampling interface

Rationale: The sampling throughput of immediate drop-on-demand technology (I.DOT) coupled with an open port sampling interface (OPSI) is limited by software communication. To enable much-needed high-throughput mass spectrometry (MS) analysis capabilities, a novel software was developed that allows for flexible sample selection from a 96-well plate and for maximized analysis throughput using I.DOT/OPSI-MS coupling. Methods: Wells of a 96-well I.DOT plate were filled with propranolol solution and were used to test maximum sampling throughput strategies to minimize analysis time. Demonstration of chemical reaction monitoring was done using acid-catalyzed ring closure of 2,3-diaminonaphthalene (DAN) with nitrite to form 2,3-naphthotriazole (NAT). Analytes were detected in positive electrospray ionization mode using selected reaction monitoring. Results: A maximum throughput of 1.54 s/sample (7.41 min/96-well plate with three technical replicates) was achieved, and it was limited by the peak width of the MS signal resulting in an occasional slight overlap between the peaks. Further, relative standard deviation was 10 ± 1% with all tested sampling strategies. Chemical reaction monitoring of DAN to NAT using nitrite was successfully accomplished with 2 s/sample throughout showing almost complete transformation in 10 min with no signal overlap. Conclusions: This work illustrates the development of a noncontact, automated I.DOT/OPSI-MS system with improved throughput achieved through an optimized software interface. Its achievable analysis time and precision make it a viable approach for drug discovery and in situ reaction monitoring studies.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Heterogeneous Mixtures of Dictionary Functions to Approximate Subspace Invariance in Koopman Operators: Why Deep Koopman Operators Work

Abstract Koopman operators model nonlinear dynamics as a linear dynamic system acting on a nonlinear function as the state. This nonstandard state is often called a Koopman observable and is usually approximated numerically by a superposition of functions drawn from a dictionary . In a widely used algorithm, extended dynamic mode decomposition (EDMD), the dictionary functions are drawn from a fixed class of functions. Deep learning combined with EDMD has been used to learn novel dictionary functions in an algorithm called deep dynamic mode decomposition (deepDMD). The learned representation both (1) accurately models and (2) scales well with the dimension of the original nonlinear system. In this paper, we analyze the learned dictionaries from deepDMD and explore the theoretical basis for their strong performance. We explore State-Inclusive Logistic Lifting (SILL) dictionary functions to approximate Koopman observables. Error analysis of these dictionary functions show they satisfy a property of subspace approximation, which we define as uniform finite approximate closure. Typically, a Koopman dictionary’s nonlinear functions are homogeneous. In this paper, we discover that structured mixing of heterogeneous dictionary functions drawn from different classes of nonlinear functions achieve the same accuracy and dimensional scaling as the deep-learning-based deepDMD algorithm Yeung et al. ( In: 2019 American Control Conference (ACC), 2019). We specifically show this by building a heterogeneous dictionary comprised of SILL functions and conjunctive radial basis functions (RBFs). This mixed dictionary achieves similar accuracy and dimensional scaling to deepDMD with an order of magnitude reduction in parameters, while maintaining geometric interpretability. These results strengthen the viability of dictionary-based Koopman models to solving high-dimensional nonlinear learning problems.

Johnson, Charles A.↗
DOE OSTI2025

Hydraulic constraints to stomatal conductance in flooded trees

Stomatal closure is a pervasive response among trees exposed to flooded soil. We tested whether this response is caused by reduced hydraulic conductance in the soil-to-leaf hydraulic continuum (k total ), and particularly by reduced root hydraulic conductance (k root ), which has been widely hypothesized. We tracked stomatal conductance at the leaf level (g s ) and canopy scale (G s ) along with physiological conditions in two temperate tree species, Magnolia grandiflora and Quercus virginiana, that were subjected to flood and control conditions in a greenhouse experiment. Flooding reduced g s , G s , k root and k total . Path analysis showed strong support for direct effects of k total on g s and for flood duration on k total , but not k root on k total . A process-based model that accounted for the k total reduction predicted the timeseries of G s in flood and control treatment trees reasonably well (predicted versus observed G s R 2 = 0.80 and 0.51 for M. grandiflora and Q. virginiana, respectively). However, accounting only for k root reduction in flooded trees was insufficient for predicting observed G s reduction. Together, these results suggest that hydraulic constraints were not limited to roots and highlight the need to account for flooding effects on k total when projecting forest ecosystem function using process-based models.

Plant stress↗
DOE OSTI2021

A Time-Dependent Directional Damage Theory for Brittle Rocks Considering the Kinetics of Microcrack Growth

In this work, a novel microcrack damage theory for describing the time-dependent behavior of brittle rocks is proposed. Instead of using scalar or tensorial variables to approximate the effect of microcracks, the concept of directional damage is introduced by upscaling the directional distribution density of microcracks through the noninteracting homogenization scheme. The contribution of an individual crack on global strain is first derived from the crack-opening displacements including sliding, closure, and dilation. The model is then cast in the free energy form, and the range of parameters ensuring compliance with the second law of thermodynamics is derived rigorously. Considering that the subcritical propagation of microcracks is the dominant mechanism behind the time-dependent deformation and failure of brittle rocks, the driving force for directional damage is defined as the corresponding energy release rate, and the effect of time is introduced through a damage evolution law inspired from the kinetics of subcritical crack growth. The model is applied to simulate the behavior of basalt rock under various loading conditions. In addition to predicting the stress–strain–time response, the model offers enhanced resolution in representing the anisotropic characteristics of microcrack-induced damage in brittle rocks and their time-dependent evolutions.

58 GEOSCIENCES↗
DOE OSTI2023

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗
DOE OSTI2024

Revisiting a Drag Partition Model For Canopy-Like Roughness Elements

Turbulent flows over a large surface area (S) covered by n obstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter area A and a shelter volume V. The individual sizes of A and V without any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element to n/S elements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducing A or V with each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction in A and V only when the ratio of n/S is small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughout S. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.

canopy turbulence↗
DOE OSTI2026

Characterizing Turbulence at a Forest Edge: Comparing Sub-Filter Scale Turbulence Models in Simulations of Flow over a Canopy

In wildfires, atmospheric turbulence plays a major role in the transfer of turbulent kinetic energy. Understanding how turbulence feeds back into a dynamical system is important, down to the varying small scales of fuel structures (i.e. pine needles, grass). Large eddy simulations (LES) are a common way of numerically representing turbulence. The Smagorinsky model (1963) serves as one of the most studied sub-grid scale representations in LES. In this investigation, the Smagorinsky model was implemented in HIGRAD/FIRETEC, LANL’s coupled fire-atmosphere model. This study was motivated by the need to quantitatively investigate the vorticity budget equation in HIGRAD/FIRETEC. The Smagorinsky turbulent kinetic energy (TKE) was compared to FIRETEC’s 1.5-order TKE eddy-viscosity subgrid-scale model, known as the Linn turbulence model. This was done in simulations of flow over flat terrain with a homogeneous, cuboidal canopy in the center of the domain. Examinations of the modeled vertical TKE profile and turbulent statistics at the leading edge, and throughout the canopy, show that the Smagorinsky model provides comparable results to that of the original closure model posed in FIRETEC.

58 GEOSCIENCES↗
DOE OSTI2021

Advanced mass balance characterization and fractionation of algal biomass composition

Abstract Opportunities associated with biomass production and bioproduct isolation from algae-derived feedstocks are plentiful and promising; however, there are challenges associated with realizing these applications. One of the most important, and often overlooked, challenges is the lack of availability of a strong foundation of compositional analysis methods validated on microalgal biomass. Currently, compositional analysis in algae is dominated by the use of interference-prone methods, a lack of full mass balance accounting, and the use of top-down approaches that bin all unaccounted-for mass into a single category, such as carbohydrates. We present here an approach based on a bottom-up algal biomass characterization aimed at moving towards, and highlighting the importance of, full and accurate mass closure to achieve the maximum economic potential from a sustainable and renewable feedstock. Algal biomass representing three genera, Nannochloropsis , Scenedesmus , and Monoraphidium , was subjected to a cell rupture and fractionation process, followed by detailed characterization of each fraction to determine the partitioning of measured and unknown components. The goal of this work is to identify where the missing components partition, and develop a strategy to close the mass balance or identify the unknowns, while utilizing a rigorous characterization approach for characterizing algal biomass. Although only 75–80% of the biomass was accounted for, the fractionation approach utilized here provides key insight into possible chemical components for future investigations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Progress of Divertor Heat and Particle Flux Control in EAST for Advanced Steady-State Operation in the Last 10 Years

Active control of the excessively high heat and particle fluxes on the divertor target plates is of fundamental importance to the steady state operation of tokamaks, especially for fusion reactors. A series of experiments have been carried out on this critical issue to relieve the plasma-wall interactions in the experimental advanced superconducting tokamak (EAST) in the last ten years, not only contributing to the long pulse operation of EAST itself, but also providing physical understandings and potential techniques to the next-generation devices like ITER. Here we have characterized the power deposition pattern and broadened the divertor footprint width effectively. The plasma-wetted area is actively handled using either 3-dimentional edge magnetic topology or advanced plasma equilibrium, thereby peak heat flux around the strike point is reduced. Active control of detachment or radiation compatible with core plasma performance has progressed significantly in very recent years, with a series of active feedback control modules developed and utilized successfully, based on the divertor physics advances with both experiments and simulation. The upper divertor of EAST was upgraded from graphite to active water-cooling ITER-like tungsten in 2014, exhibiting much enhanced heat removal capability. As for the particle exhaust including both fueling and impurity particles, in addition to wall conditioning and impurity source control, the efficiency of particle flux exhaust is optimized by making full use of the divertor closure and the plasma drifts in both scrape-off layer and divertor volume. These heat and particle exhaust advances contribute greatly to a series of EAST achievements like H-mode operation over 100 s. A brief near-term plan on the integrated control of divertor plasma-wall interactions in long-time scale will also be introduced, aiming to provide favorable divertor operation solution for ITER and CFETR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Drought conditioning of rhizosphere microbiome influences maize water use traits

Background and Aims: Beneficial plant–microbe interactions can improve plant performance under drought; however, we know less about how drought-induced shifts in microbial communities affect plant traits. Methods: We cultivated Zea mays in fritted clay with soil microbiomes originating from contrasting environments (agriculture or forest) under two irrigation treatments (well-watered or water limited). Using this design, we investigated whether water conditioning was carried forward through the microbiome to affect a subsequent plant cohort that was subjected to either a well-watered or water limited treatment. Results: Regardless of the microbiome-origin, plants inoculated with a microbiome from a water limited legacy had traits that allowed them to avoid stress but conserve water. They produced longer roots to explore soil, generated greater soil dissolved organic carbon, potentially stimulating the microbiome, and slower soil water content loss during drought. A well-watered legacy resulted in plants that delayed permanent stomatal closure and higher photosynthetic nitrogen use efficiency. In plants with a forest-originated microbiome, a well-watered legacy and water treatment also resulted in higher rates of photosynthesis and stomatal conductance. Conclusion: These results demonstrate that soil microbiomes can be developed to influence plant drought performance, impacting crop resilience, using short-term microbial conditioning.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

The Influence of Residual Stress on Fatigue Crack Growth Rates in Stainless Steel Processed by Different Additive Manufacturing Methods

The properties and microstructure of Type 304L stainless steel produced by two additive manufacturing (AM) methods—directed energy deposition (DED) and powder bed fusion (PBF)—are evaluated and compared. Localized heating and steep temperature gradients of AM processes lead to significant residual stress and distinctive microstructures, which may be process-specific and influence mechanical behavior. Test data show that materials produced by DED and PDF have small differences in tensile strengths but clear differences in residual stress and microstructural features. Measured fatigue crack growth rates (FCGRs) for cracks propagating parallel to and perpendicular to the build directions differ between the two AM materials. To separate the influences of residual stress and microstructure, K-control test procedures with decreasing and constant stress intensity factor ranges are used to measure FCGRs in the near-threshold regime (crack growth rates ≤ 1 × 10 −8 m/cycle). Residual stress is quantified by the residual stress intensity factor, K res , measured by the online crack compliance method. Correcting the FCGR data for differences in K res brings results for specimens of the two AM materials into agreement with each other and with results for wrought specimens, when the latter are corrected for crack closure. Differences in microstructure and tensile strength have an insignificant influence on FCGRs in these tests.

additive manufacturing↗
DOE OSTI2024

Evaluation of sintered bauxite proppant for binary enhanced geothermal systems

Abstract Solid granular proppant particles are widely used in oil and gas development to sustain permeability through fractures after hydraulic stimulation. Similar proppants are of interest for geothermal applications where the goal of sustaining permeability is the same, but the harsh geothermal environment risks rapid proppant degradation that will reduce fracture permeability. Here, we present proppant conductivity experiments using saw cut granite, tensile fractured granite, and aluminum control specimens packed with sintered bauxite ceramic proppants at concentrations of 0.0, 0.1, and 1.0 kg/m 2 . Simulated geothermal conditions included temperatures up to 130 °C and normal closure stresses up to 60 MPa. Compared to unpropped specimens, peak fracture conductivity enhancement was up to 6 orders of magnitude. At simulated geothermal conditions, proppants were able to sustain fracture conductivity over 60 h, but chemical dissolution and decreasing permeability over time were evident. Irreversible conductivity reductions with crushing and embedment of proppants during loading stages were also observed. Overall, sintered bauxite proppant remains a promising option for low-temperature binary-cycle enhanced geothermal systems.

58 GEOSCIENCES↗
DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2022

Sensitivity analysis of generic deep geologic repository with focus on spatial heterogeneity induced by stochastic fracture network generation

Geologic Disposal Safety Assessment Framework is a state-of-the-art simulation software toolkit for probabilistic post-closure performance assessment of systems for deep geologic disposal of nuclear waste developed by the United States Department of Energy. This paper presents a generic reference case and shows how it is being used to develop and demonstrate performance assessment methods within the Geologic Disposal Safety Assessment Framework that mitigate some of the challenges posed by high uncertainty and limited computational resources. Variance-based global sensitivity analysis is applied to assess the effects of spatial heterogeneity using graph-based summary measures for scalar and time-varying quantities of interest. Behavior of the system with respect to spatial heterogeneity is further investigated using ratios of water fluxes. This analysis shows that spatial heterogeneity is a dominant uncertainty in predictions of repository performance which can be identified in global sensitivity analysis using proxy variables derived from graph descriptions of discrete fracture networks. New quantities of interest defined using water fluxes proved useful for better understanding overall system behavior.

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

An Unmanned Aerial System (UAS) for concurrent measurements of solar-induced chlorophyll fluorescence and hyperspectral reflectance toward improving crop monitoring

Unmanned aerial system (UAS)-based remote sensing can provide high spatial- and temporal-resolution crop monitoring for precision management. Existing crop monitoring UAS primarily use conventional techniques such as multi-spectral broadband vegetation indices (VIs) to track canopy structure/biomass. Recently developed lightweight hyperspectral sensors can track crop physiology and performance via complementary signals such as solar-induced chlorophyll fluorescence (SIF) and photochemical reflectance index (PRI) that can be derived from hyperspectral reflectance. However, few existing UAS can acquire high-quality SIF and hyperspectral reflectance. Here, we designed a novel UAS that simultaneously captures far-red SIF and hyperspectral reflectance, using a single bifurcated fiber and motorized arm to measure downwelling and upwelling irradiance. The UAS was tested over (1) a heterogeneous corn field with spatially varying crop yield potential and (2) a set of soybean and corn plots under differing nutrient treatments. Our UAS can maintain stability under low-to-moderate winds (2–3 m/s) with <0.4 m geolocation accuracy, <0.1 m altitude drift, and <1.5° error of the fiber optic from nadir position during scanning. Weekly seasonal campaigns reveal that SIF outperforms conventional VIs such as the normalized difference vegetation index (NDVI) for distinguishing plots with different crop yield potential, especially after canopy closure when NDVI tends to saturate. The UAS can capture the diurnal dynamics of SIF and PRI for both corn and soybean. At the seasonal scale, SIF acquired from the UAS correlates well with that from a fixed SIF tower system ( R 2 = 0.81), which measured the same target albeit with different footprints, demonstrating the capability of UAS in characterizing the seasonal progression of crop activity. In conclusion, our newly developed UAS provides high-quality SIF and hyperspectral reflectance, facilitating mechanistic understanding of the physiological control on photosynthesis dynamics. Our findings also imply that this UAS could enable extension of fixed tower-based systems to monitor multiple targets from diurnal to seasonal scales and improve monitoring of heterogeneous fields.

59 BASIC BIOLOGICAL SCIENCES↗