Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical approximations”

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

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

At least 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Multiscale Ecosystem for solving Maxwell-Schrodinger equations of open quantum systems (OpenMS)

Light-matter interactions play an important role in many branches of physics, chemistry, energy, and materials science. In the strong coupling regime, light-matter interactions are able to tune the materials properties via the formation of new quasiparticles, such as plasmons and polaritons. However, theoretical and numerical modeling of the light-matter interaction-mediated processes remain a big challenge because light-matter interactions are fundamentally multiscale and multiphysics problems involving multiple interactions between electrons, nuclei, and photons at different time/length scales. In the current community, light-matter interactions were treated at different levels of theoretical complexity in quantum chemistry and quantum optics. In quantum optics or quantum photonics, the matter is usually simplified as a few-level system, and the light is treated quantum-mechanically. On the other hand, quantum chemistry explores first-principles methods, including both single-particle and many-body-based techniques, to describe the electronic properties of matter in detail. However, light is usually prescribed as a classical electromagnetic field, and the light-matter interaction is taken into account as an external potential via classical approximations. This software is designed to fill current modeling shortcomings by delivering a first-ever scalable multiscale platform for simulating light-matter interactions in realistic electromagnetic environments. The software solves Maxwell and Schrodinger equations self-consistent on the heterogeneous platforms. It implements HF/DFT, TDDFT, and coupled-cluster counterparts for light-matter interactions and adopts modular programming to offload massively parallel algorithms on a large number of CPU/GPUs.

Zhang, Yu↗
DOE OSTI2021

Large-Scale Welding Process Simulation by GPU Parallelized Computing

The computational design of industrially relevant welded structures is extremely time consuming due to coupled physics and high nonlinearity. Previously, most welding distortion and residual stress simulations have been limited to small coupons and reduced order (from three-dimensional [3D] to two-dimensional [2D]), or inherent strain approximations were used for large structures. In this current study, an explicit finite element code based on a graphics processing unit was utilized to perform 3D transient thermomechanical simulation of structural components during welding. Laser brazing of aluminum alloy panels as representative of automotive manufacturing scenarios was simulated to predict out-of-plane distortion under different clamping conditions. The predicted deformation pattern and magnitude were validated by laser scanning data of physical assemblies. In addition, the code was used to investigate residual stresses developed during multipass arc welding of a nuclear industry pressurizer surge nozzle and subsequent welding repair where a 3D simulation was necessary. Taking the experimental data as reference, the 3D model predicted better residual stress distribution than a typical 2D asymmetrical model. Stress evolution in welding repair was also presented and discussed in this study. Furthermore, the efficient numerical model made it feasible to use integrated computational welding engineering to simulate welding processes for large-scale structures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Sparse low rank approximation of potential energy surfaces with applications in estimation of anharmonic zero point energies and frequencies

We introduce a method that exploits sparse representation of potential energy surfaces (PES) on a polynomial basis set selected by compressed sensing . The method is useful for studies involving large numbers of PES evaluations, such as the search for local minima, transition states, or integration. We apply this method for estimating zero point energies and frequencies of molecules using a three step approach. In the first step, we interpret the PES as a sparse tensor on polynomial basis and determine its entries by a compressed sensing based algorithm using only a few PES evaluations. Then, we establish a rank reduction strategy to compress this tensor in a suitable low-rank canonical tensor format using standard tensor compression tools. This allows representing a high dimensional PES as a small sum of products of one dimensional functions. Finally, a low dimensional Gauss–Hermite quadrature rule is used to integrate the product of sparse canonical low-rank representation of PES and Green’s function in the second-order diagrammatic vibrational many-body Green’s function theory (XVH2) for estimation of zero-point energies and frequencies. Numerical tests on molecules considered in this work indicate a more efficient scaling of computational cost with molecular size as compared to other methods.

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

Reynolds-Averaged Turbulence Modeling Using Deep Learning with Local Flow Features: An Empirical Approach

Reynolds-Averaged Navier-Stoke (RANS) models offer an alternative avenue in predicting flow characteristics when the corresponding experiments are difficult to achieve due to geometry complexity, limited budget, or knowledge. RANS models require the knowledge of subgrid scale physics to solve conservation equations for mass, energy, and momentum. Mechanistic turbulence models, such as k-epsilon, are generally evaluated and calibrated for specific flow conditions with various degrees of uncertainty. These models have limited capability to assimilate a substantial amount of data due to model form constraints. Meanwhile, deep learning (DL) has been proven to be universal approximators with the potential to assimilate available, relevant, and adequately evaluated data. Moreover, deep neural networks (DNNs) can create surrogate models without knowing function forms. Such a data-driven approach can be used in updating fluid models based on observations as opposed to hard-wiring models with precalibrated correlations. The paper presents progress in applying DNNs to model Reynolds stress using two machine learning (ML) frameworks. A novel flow feature coverage mapping is proposed to quantify the physics coverage of DL-based closures. It can be used to examine the sufficiency of training data and input flow features for data-driven turbulence models. The case of a backward-facing step is formulated to demonstrate that not only can DNNs discover underlying correlation behind fluid data but also they can be implemented in RANS to predict flow characteristics without numerical stability issues. Finally, the presented research is a crucial stepping-stone toward the data-driven turbulence modeling, which potentially benefits the design of data-driven experiments that can be used to validate fluid models with ML-based fluid closures.

42 ENGINEERING↗
DOE OSTI2023

Raman Wavelength Conversion in Ionic Liquids

Here, we explore the use of room-temperature ionic liquids (ILs) as Raman wavelength converters. ILs provide an engineerable framework to design suitable liquids for wavelength conversion over a broad spectral range, through careful selection of the molecular structures of the IL anions and cations so that specific characteristics can be obtained, such as a desirable Raman shift, low Brillouin scattering, and good optical transmission in the pump and Stokes wavelengths. Applying such criteria, we demonstrate that 1-ethyl-3-methylimidazolium dicyanamide (EMIM DCA) is an effective medium for conversion of 532-nm pulses from a Q-switched Nd :YAG laser to 603 nm. This corresponds to an approximate 2200 cm –1 shift, which can be used to generate mid-infrared radiation through subsequent difference frequency generation for optical pumping of CO 2 lasers. Threefold-higher Raman conversion efficiency is obtained in EMIM DCA compared with water under identical conditions in a proof-of-principle single-pass conversion setup, resulting in an efficient generation of multimillijoule, <6 ns duration, high-quality orange laser pulses in a wavelength region that is difficult to access at high energies. Consequently, we examine ILs representing two other classes of Raman-active functional groups and obtain conversion up to the fifth-order Stokes shift and first anti-Stokes shift. Through the tunable selection of their components and their useful dynamical properties, ILs provide a platform for efficient, simple, and alignment-tolerant high-energy Raman shifting with numerous industrial and technological applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Spectral properties and enhanced superconductivity in renormalized Migdal-Eliashberg theory

Migdal-Eliashberg theory describes the properties of the normal and superconducting states of electron-phonon-mediated superconductors based on a perturbative treatment of the electron-phonon interactions. It is necessary to include both electron and phonon self-energies self-consistently in Migdal-Eliashberg theory in order to match numerically exact results from determinantal quantum Monte Carlo in the adiabatic limit. Here in this work we provide a method to obtain the real-axis solutions of the Migdal-Eliashberg equations with electron and phonon self-energies calculated self-consistently. Our method avoids the typical challenge of computing cumbersome singular integrals on the real axis and is numerically stable and exhibits fast convergence. Analyzing the resulting real-frequency spectra and self-energies of the two-dimensional Holstein model, we find that self-consistently including the lowest-order correction to the phonon self-energy significantly affects the solution of the Migdal-Eliashberg equations. The calculation captures the broadness of the spectral function, renormalization of the phonon dispersion, enhanced effective electron-phonon coupling strength, minimal increase in the electron effective mass, and the enhancement of superconductivity which manifests as a superconducting ground state despite strong competition with charge-density-wave order. We discuss surprising differences in two common definitions of the electron-phonon coupling strength derived from the electron mass and the density of states, quantities which are accessible through experiments such as angle-resolved photoemission spectroscopy and electron tunneling. An approximate upper bound on 2Δ/T c for conventional superconductors mediated by retarded electron-phonon interactions is proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Triangular lattice Hubbard model physics at intermediate temperatures

Moiré systems offer an exciting playground to study the many-body effects of strongly correlated electrons in regimes that are not easily accessible in conventional material settings. Motivated by a recent experiment on WSe 2 /WS 2 moiré bilayers [Y. Tang et al., Nature (London) 579, 353 (2020)], which realized a triangular superlattice with a small hopping (of approximately 10 Kelvin), with tunable density of holes, we explore the Hubbard model on the triangular lattice for intermediate temperatures t ≲ T < U. Employing finite-temperature Lanczos calculations and closely following the fitting protocols used in the experiment, we recover the observed trends in the reported Curie-Weiss temperature Θ with filling, using the reported interaction strength U/t=20. Here, we focus on the large increase of |Θ| on decreasing the density below half filling and the sign change of Θ at higher fillings, which signals the onset of ferromagnetism. The increase in |Θ| is also seen in the t–J model (the low-energy limit of the Hubbard model) in the intermediate temperature range, which we clarify is opposite to the trend in its high temperature limit. Differences between the low, intermediate, and high temperature regimes are discussed. Our numerical calculations also capture the crossover between short-range antiferromagnetic to ferromagnetic order in the intermediate temperature regime, a result broadly consistent with the experimental findings. We find that this behavior is a finite-temperature remnant of the underlying zero-temperature phase transition, which we explore with ground-state density matrix renormalization group calculations. We provide evidence of ferromagnetism characterized by weak (but robust) correlations that explain the small Θ seen in the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Infrared thermography NDT for in-situ defect detection in sandwich composite panel manufacturing

Composite manufacturing presents numerous challenges, as defects can arise from various sources throughout the process. In sandwich composite structures, the integration of a foam core introduces additional complexity and increases the likelihood of defect formation like delamination. To mitigate these issues and reduce the risk of future structural failures, in-situ monitoring during manufacturing is essential. This study investigates infrared (IR) thermography as a non-destructive technique for detecting manufacturing defects in foam-core sandwich composite panels under thermally excited conditions representative of in-situ processing. A stationary FLIR A8590 IR camera (640 × 512 pixels, 30Hz, 17mm lens, 9 ft stand-off distance) was used to monitor prefabricated panels subjected to controlled external heating simulating compression molding and resin cure exotherm. Interlaminar delamination defects with characteristic sizes ranging from 0.25 × 0.25in² to 5 × 5in² produced measurable surface temperature depressions of approximately 4–10°C during transient cooling, exceeding the effective noise floor of the camera by more than two standard deviations. Thicker laminates exhibited prolonged defect detectability windows due to increased thermal diffusion time. In contrast, embedded Teflon inclusions generated weak thermal contrasts of ≤ 3°C, approaching the measurement noise floor, due to limited thermal property contrast with the surrounding glass fiber composite. These results establish quantitative detectability limits for stationary thermographic inspection of sandwich composite panels under manufacturing-representative thermal cycles.

Barakat, Abdallah [ORNL] (ORCID:0000000296141398)↗
DOE OSTI2019

Sharp Drop in the Mobility of Holes with a Decrease in Their Two-Dimensional Concentration by an External Voltage in Boron δ-Doped Diamond Layers

It is shown, both analytically and numerically, that the hole mobility in boron δ-doped (i.e., with a thickness on the order of several lattice constants) diamond layers drops with a decrease in their two-dimensional concentration during depletion of the δ-doped layer by an external voltage. This drop is most pronounced for the maximum initial two-dimensional hole concentrations ~ 3 × 10{sup 13} cm{sup –2} (limited from above by the condition for the possibility of their significant reduction without the electrical breakdown of diamond). This is explained by a decrease in the degree of screening of scattering Coulomb potentials of ionized boron atoms and an increase in the scattering efficiency of degenerate holes from these atoms due to a drop in the kinetic energy of holes. The corresponding calculations of the transport scattering cross section of holes are performed beyond the Born approximation (i.e., perturbation theory), which is inapplicable for boron δ-doped diamond layers. The predicted effect can be used to increase the efficiency of source-to-drain current modulation by a gate voltage in diamond field-effect transistors (FETs) with δ-doped conduction channels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Atomic Mixing of Initially Separated Mixtures

We consider the atomic mixing of two materials, initially separated in space, at fixed temperature– pressure conditions [T, p], and zero mass-mean velocity. Each material is composed of a specified set of atomic species whose density corresponds to the specieswise real Equation of State (EOS) $\tilde{&#961}$ i [T, p]. Because there are no spatial gradients in the mean [T, p] state, the only mechanism for mixing is a specieswise nonequilibrium thermal velocity distribution. The recipe for nonequilibrium species mass fluxes is taken directly from Hirschfelder, Curtiss, and Bird, who express the diffusion velocities in terms of gradients in the molar concentrations and pairwise coefficients of diffusion. EOS data from the SESAME library is used for all species. Data for the specieswise self diffusion coefficients is adopted from the Orbital Free Molecular Dynamics (OFMD) simulations of Ticknor et al.; the so-called Darken approximation is used to express the mutual diffusion coefficients of a mixture in the manner used by White et al. A one-dimensional spatial domain is used; with the left half occupied by one multispecies material and the right half occupied by the other, and with zero–gradient boundary conditions at the ends of the domain. Hence the total mass in the domain is a constant. At time t = 0 the materials are separated. Accordingly, spatial gradients at the separation point are infinite and the species fluxes are also infinite. As t → $\textbf{&#x221E;}$ the spatial gradients approach zero and a uniform set of species concentrations is obtained; at this point the system is said to be at “thermodynamic equilibrium.” Analytic solution of the problem is not possible in general due to the nonlinear concentration dependence of the species fluxes. Hence we utilize the numerical solution apparatus in Mathematica for obtaining accurate time-dependent solutions for the mixing of various initial species masses for conditions of interest in laser-driven Inertial Confinement Fusion. Concentration profiles and the speed of propagation of the mixing fronts are discussed in light of the nonlinearity of the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Stage-local partitioned two-step runge-kutta methods for large systems of ordinary differential equations

We introduce stage-local partitioned two-step Runge-Kutta methods are an extension of standard two-step Runge-Kutta methods, which are an alternative to the standard additive two-step Runge-Kutta methods currently existing in the literature. Furthermore, these new schemes are designed with an eye towards truly N-partitioned systems and leverage local stage approximations to make several computationally interesting approximations viable. Specifically, the focus on local stage approximations makes possible the construction of truly asynchronous schemes, in the parallel sense, possible. In addition, we show that an implicit-explicit approach to these schemes can lead to methods that require the inversion of only local nonlinear systems.

Applied Dynamical Systems↗
DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Sensitivity of the Shallow‐To‐Deep Convective Transition to Moisture and Wind Shear in the Amazon

Deep convection is the primary influence on weather and climate in tropical regions. However, understanding and simulating the shallow‐to‐deep (STD) convective transition has long been challenging. Here, we conduct high‐resolution numerical simulations to assess the environmental controls on the evolution of isolated convection in the Amazon during the wet season. The large‐scale forcing derived through a constrained variational analysis approach for the GoAmazon2014/5 Experiment is used in the simulations. Through sensitivity experiments, we examine the relative importance of moisture and wind shear in controlling the shallow‐to‐deep convective transition for isolated convective events. Convection exhibits the greatest sensitivity to humidity within the lowest 1.5 km, where a 4 mm reduction in column water vapor nearly suppresses ice water formation on deep convective days. In contrast, a reduction in column water vapor in the free troposphere by a factor of two or more is necessary to produce a comparable impact on convection. Increasing low‐level wind speed from 6 to 9 m s -1 enhances afternoon deep convection, raising the cloud ice mixing ratio by approximately 25%. Conversely, upper‐level wind shear reveals the weakest correlation with daytime convection in our simulations. Our results help characterize the role of moisture and wind shear on the STD transition and our understanding of the underlying mechanisms.

environmental controls↗
DOE OSTI2024

Towards dynamical low-rank approximation for neutrino kinetic equations. Part I: Analysis of an idealized relaxation model

Dynamical low-rank approximation (DLRA) is an emerging tool for reducing computational costs and provides memory savings when solving high-dimensional problems. Here, in this work, we propose and analyze a semi-implicit dynamical low-rank discontinuous Galerkin (DLR-DG) method for the space homogeneous kinetic equation with a relaxation operator, modeling the emission and absorption of particles by a background medium. Both DLRA and the discontinuous Galerkin (DG) scheme can be formulated as Galerkin equations. To ensure their consistency, a weighted DLRA is introduced so that the resulting DLR-DG solution is a solution to the fully discrete DG scheme in a subspace of the standard DG solution space. Similar to the standard DG method, we show that the proposed DLR-DG method is well-posed. We also identify conditions such that the DLR-DG solution converges to the equilibrium. Numerical results are presented to demonstrate the theoretical findings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Poisson–Nernst–Planck Model of Ion Transport and Interface Segregation in Metal–Insulator–Semiconductor Structures and Solar Cells

A numerical model that describes the transport of mobile ionic species in metal–insulator–semiconductor (MIS) and photovoltaic (PV) devices subject to temperature and voltage stress is presented. The finite element method (FEM) is used to solve the Nernst–Planck equation while imposing Poisson's equation self-consistently as a restriction for the electrostatic potential. This allows the contribution of the ionic species to the potential to be taken into account. Using a variational formulation eases the implementation of diverse boundary conditions, including the incorporation of segregation kinetics at the device interfaces. Segregation across the dielectric–semiconductor interface is relevant to modeling the electronic device degradation in systems where contamination reaches the semiconductor. The model in closed systems with no-flux boundary conditions is validated first. In the limiting case of low contamination levels with respect to the gate bias, the FEM solution matches analytically derived approximations. Then, the implementation is broadened to include an open boundary at the dielectric–semiconductor interface to account for leakage of ions. The predicted time dependence of the flatband voltage in Na-contaminated MIS test structures agrees well with measurements. This model successfully captures the role of long-range ion transport at concentrations of relevance to electronic and PV device instability and neuromorphic computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Asymptotically compatible reproducing kernel collocation and meshfree integration for the peridynamic Navier equation

Here, we study reproducing kernel (RK) collocation method for peridynamic Navier equation. In the first part, we apply a linear RK approximation to both displacement and dilatation, and then back-substitute dilatation and solve the peridynamic Navier equation in a pure displacement form. The RK collocation scheme converges to the nonlocal limit for a fixed nonlocal interaction length and also to the local limit as nonlocal interactions vanish. The stability is shown by comparing the collocation scheme with the standard Galerkin scheme using Fourier analysis. In the second part, we apply the RK collocation to the quasi-discrete peridynamic Navier equation and show its convergence to the correct local limit when the ratio between the nonlocal length scale and the discretization parameter is fixed. The analysis is carried out on a special family of rectilinear Cartesian grids for the RK collocation method with a designated kernel with finite support. We assume the Lamé parameters satisfy λ ≥ μ to avoid extra assumptions on the nonlocal kernel. Finally, numerical experiments are conducted to validate the theoretical results.

42 ENGINEERING↗
DOE OSTI2021

Transient cycling of nitrogen, organic carbon and oxygen within the free-flowing Columbia River corridor: Linking exposure time dependent biogeochemical reactions to river stage fluctuations (Final project report)

The objective of this project was to develop next-generation techniques for representing the transformations of complex reaction networks like those impacting transient river-corridors, then apply them to develop upscaling tools. The project contributed significantly to the development of Lagrangian “mass transfer particle tracking” (MTPT) tools that are the first numerical framework that explicitly separates mixing and spreading processes; a crucial distinction for accurate representations of reaction rates across scales. The MTPT approach uses an operator splitting scheme where physical transport processes are simulated using classical random walk methods and each particle is treated as a “container” that carries with it any number of chemical species. Mixing between particles is simulated using a colocation-probability based mass transfer kernel and reactions are evaluated on each particle after mixing. The method is stable, accurate, and also allows the explicit modeling of residence time distributions of the individual containers of mass. The MTPT scheme allowed us to explore several research questions related to reactive transport in river corridors and complex reaction networks. We found that it is possible to quantify the exposure time distributions (ETDs) of reactants and that this does suggest some pathways to upscaling. However, comparable approximations can be obtained using simplified (reduced dimensionality) MTPT simulations in less time with similar uncertainty, so development of extensive ETD-based methods was not productive. The major findings are that i) the MTPT schemes are robust and highly accurate across a range of arbitrarily complex reactions, ii) the methods can be efficiently parallelized and the parallel performance characteristics of MTPT are predictable, and iii) the MTPT tools allow simultaneous tracking of residence time. The development of these research tools into efficient software packages continues. The methods are already available to the community because open-source, working examples have been included with all publications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Uptake and elimination of 137Cs in aquatic biota inhabiting a contaminated effluent canal

Radiocesium (137Cs) is among the most prevalent and concerning radionuclide released from anthropogenic activities due to its relatively slow decay time and high bioavailability. Numerous studies have reported activity concentrations of 137Cs in exposed biota, but few have assessed detailed toxicokinetics. The objective of our study was to quantify and compare the uptake and elimination of 137Cs in aquatic taxa from multiple classes with different ecological roles inhabiting a canal contaminated from nuclear production activities on the Savannah River Site in South Carolina, USA. We tracked accumulation and elimination of 137Cs in four understudied aquatic species: one aquatic macrophyte [American white water lily (Nymphae odorata), two benthic associated omnivorous species (American bullfrog (Lithobates catesbeianus) and red swamp crayfish (Procambarus clarkii), and one pelagic carnivorous fish (Eastern mosquitofish (Gambusia holbrooki)]. We monitored accumulation of 137Cs from initial exposure until equilibrium for biota restricted within the canal and subsequently tracked elimination of 137Cs over time by moving organisms from the canal to a nearby reference pond. Results indicated bullfrog tadpoles reached 95 % equilibrium in approximately 9 days, the fastest among the species evaluated. After bullfrog tadpoles, the next fastest time to equilibrium was achieved by red swamp crayfish (51 days) followed by mosquitofish (86 days). At equilibrium, bullfrog tadpoles were found to have similar 137Cs activity concentrations to mosquitofish (3.16 and 3.57 Bq/g, respectively), and both were significantly higher than red swamp crayfish (1.77 Bq/g). For elimination, we found biological half-lives to range between 8 and 43 days with bullfrog tadpoles eliminating 137Cs the fastest followed by water lilies, mosquitofish, and red swamp crayfish. Overall, we found the uptake and elimination of 137Cs to vary widely among our studied species, but benthic organisms are not at an increased risk to accumulate 137Cs compared to organisms in less frequent contact with the sediment.

Benthic organism Bioaccumulation Bioindicator Elim↗