SEARCH · Engineering Papers
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.
We study spectra output from spatially proximal DPRK underground nuclear explosions from data challenged by several constraints.
SUMMARY The spatial correlation of earthquake ground motion intensity can be measured from strong motion data; however, the data used in past studies is sparsely sampled in space, and only the interstation distance was considered as a correlation variable. These limitations mean that we have only weak constraints on the true correlation structure of ground motion and that potentially important aspects of spatial correlation are unconstrained. In this study, we combine a large-N seismic array and graph analytics to explore this issue at a local scale using small local and regional earthquakes. Our result suggests site conditions, and how they interact with the incident seismic wavefield, strongly condition the spatial correlation of ground motion. Future progress in characterizing ground motion spatial variability will require dense wavefield measurements, either through nodal deployments, or perhaps distributed acoustic sensing measurements, of seismic wavefields. Aftershock sequences of major earthquakes would provide particularly data-rich targets of opportunity.
As one of the major agricultural production areas in the world, the United States (U.S.) Midwest plays a vital role in the global food supply and agricultural ecosystem services. Although significant efforts have been made in modeling the carbon cycle dynamics over this area, large uncertainty still exists in the previous simulations in terms of reproducing individual components of the carbon cycle and their responses to environmental variability. Here we evaluated the performance of an advanced agroecosystem model, ecosys, in simulating carbon budgets over the U.S. Midwest, considering both the magnitude of carbon flux/yield and its response to environmental (climate and soil) variability. We conducted model simulations and evaluations at 7 cropland eddy-covariance sites as well as over 293 counties of Illinois, Indiana, and Iowa in the U.S. Midwest. The site-level simulations showed that ecosys captured both the magnitude and seasonal patterns of carbon fluxes (i.e., net ecosystem carbon exchange (NEE), ecosystem gross primary production (GPP), and ecosystem respiration (Reco)), leaf area index (LAI), and dynamic plant carbon allocation processes, with R2 equal to 0.92, 0.87, 0.87, and 0.78 for GPP, NEE, Reco, and LAI, respectively across all the sites compared with the observations. For regional scale simulations, ecosys reproduced the spatial distribution and interannual variability of corn and soybean yields with the constraints of observed yields and a new remotely sensed GPP product, with R2 of multi-year averaged simulated and observed yield equal 0.83 and 0.80 for corn and soybean, respectively. The simulated responses of carbon cycle dynamics to environmental variability were consistent with that from the empirical observations at both site and regional scales. Furthermore, our results demonstrated the applicability of ecosys in simulating the carbon cycle and soil carbon dynamics of the U.S. Midwestern agroecosystems under different climate and soil conditions.
Throughout many fields, practitioners often rely on differential equations to model systems. Yet, for many applications, the theoretical derivation of such equations and/or the accurate resolution of their solutions may be intractable. Instead, recently developed methods, including those based on parameter estimation, operator subset selection, and neural networks, allow for the data-driven discovery of both ordinary and partial differential equations (PDEs), on a spectrum of interpretability. The success of these strategies is often contingent upon the correct identification of representative equations from noisy observations of state variables and, as importantly and intertwined with that, the mathematical strategies utilized to enforce those equations. Specifically, the latter has been commonly addressed via unconstrained optimization strategies. Representing the PDE as a neural network, we propose to discover the PDE (or the associated operator) by solving a constrained optimization problem and using an intermediate state representation similar to a physics-informed neural network (PINN). The objective function of this constrained optimization problem promotes matching the data, while the constraints require that the discovered PDE is satisfied at a number of spatial collocation points. We present a penalty method and a widely used trust-region barrier method to solve this constrained optimization problem, and we compare these methods on numerical examples. Our results on several example problems demonstrate that the latter constrained method outperforms the penalty method, particularly for higher noise levels or fewer collocation points. This work motivates further exploration into using sophisticated constrained optimization methods in scientific machine learning, as opposed to their commonly used, penalty-method or unconstrained counterparts. For both of these methods, we solve these discovered neural network PDEs with classical methods, such as finite difference methods, as opposed to PINNs-type methods relying on automatic differentiation. Here, we briefly highlight how simultaneously fitting the data while discovering the PDE improves the robustness to noise and other small, yet crucial, implementation details.
Polydimethylsiloxane is an important polymeric material with a wide range of applications. However, environmental effects like low temperature can induce crystallization in this material with resulting changes in its structural and dynamic properties. The incorporation of phenyl-siloxane components, e.g., as in a poly(dimethyl-co-diphenyl)siloxane random copolymer, is known to suppress such crystallization. Molecular dynamics (MD) simulations can be a powerful tool to understand such effects in atomistic detail. Unfortunately, all-atomistic molecular dynamics (AAMD) is limited in both spatial dimensions and simulation times it can probe. Here, to overcome such constraints and to extend to more useful length- and time-scales, we systematically develop a coarse-grained molecular dynamics (CGMD) model for the poly(dimethyl-co-diphenyl)siloxane system with bonded and non-bonded interactions determined from all-atomistic simulations by the iterative Boltzmann inversion (IBI) method. Additionally, we propose a lever rule that can be useful to generate non-bonded potentials for such systems without reference to the all-atomistic ground truth. Our model captures the structural and dynamic properties of the copolymer material with quantitative accuracy and is useful to study long-time dynamics of highly-entangled systems, sequence-dependent properties, phase behaviour, etc.
Radially resolved hydrogenic isotope fraction measurement capabilities have been developed for DIII-D using the main-ion charge exchange recombination (MICER) spectroscopy system in preparation for mixed hydrogen and deuterium experiments. Constraints on the hydrogenic ion temperatures and velocities based on measurements of the impurity ion properties are required to accurately fit the spectrum. Corrections for cross sectional distortions, spatial smearing due to the halo, and a neoclassical offset between the impurity and hydrogenic toroidal rotation are applied to the constraints prior to fitting the MICER spectrum. Extensive atomic physics calculations have been performed using the FIDASIM code, which has recently been improved to allow simulations using mixtures of hydrogenic species. Overall, these results demonstrate that for the same plasma parameters, the Dα emission is 20%–30% brighter than Hα due to differences in rate coefficients associated with the different ion thermal velocities for the same temperature and therefore must be taken into consideration when calculating absolute densities. However, despite these differences, the absolute error when estimating the hydrogen isotope fraction [n H /(n H + n D )] by using the Hα radiance fraction [L Hα /(L Hα + L Dα )] is typically less than 5% due to the way the fraction is formed, making the radiance fraction a reasonably accurate estimate of the isotope fraction for most cases.
The ENIGMA science focus area (SFA) is a multi-disciplinary, multi-institutional research effort focused on addressing foundational knowledge gaps in environmental microbial communities by studying groundwater and sediment microbiomes in the shallow subsurface at the contaminated Oak Ridge Field Research Site (FRC). We seek to discover and characterize the reciprocal interactions between the microbial communities and the geochemical and geophysical parameters of the shallow subsurface within the contamination plume. The goal of the Fields lab is to help manage and conduct laboratory experiments informed by field observations using ecological and physiological approaches to study microbial populations in situ as well as ex situ (in the laboratory). In particular, we aim to study and characterize the impact of hydrological constraints on free-living and biofilm biomass and activity with increasing spatial and temporal resolution under static and flow conditions. We also aim to help characterize novel microbial groups that are present and active under relevant field conditions, including the development of molecular techniques for in situ detection as well as metabolic interactions that underlie pertinent physiology and ecology.
Aerosols mediate the radiative fluxes in clear and cloudy skies and dominate the uncertainty in the radiative forcing of climate. Relatively dense networks of aerosol optical depth measurements have been used to effectively constrain simulated aerosol optical properties, but model diversity of aerosol microphysical properties is much larger (Mann et al. 2014, Myhre et al. 2009). A better understanding of aerosol microphysics is essential as they are the fundamental pieces of information required to convert aerosol emissions information to radiative and cloud-nucleating properties that drive radiative effects and forcing of climate. Model representation of aerosol microphysical properties has lagged in part due to their inherently greater spatial variability, but also due to a lack of available observations. This project – the Printed Optical Particle Spectrometer Network (POPSnet)-Southern Great Plains (SGP) Pilot – launched the first spatially dense network of aerosol size distribution measurements over an area the size of a global model grid cell and demonstrated its use in providing valuable information regarding the spatial variability of aerosol microphysical properties and its drivers for improving model constraints.
Abstract In this work, we present a binary assembly model that can predict the co-assembly structure and spatial frequency spectra of monodispersed nanoparticles with two different particle sizes. The approach relies on an iterative algorithm based on geometric constraints, which can simulate the assembly patterns of particles with two distinct diameters, size distributions, and at various mixture ratios on a planar surface. The two-dimensional spatial-frequency spectra of the modeled assembles can be analyzed using fast Fourier transform analysis to examine their frequency content. The simulated co-assembly structures and spectra are compared with assembled nanoparticles fabricated using transfer coating method are in qualitative agreement with the experimental results. The co-assembly model can also be used to predict the peak spatial frequency and the full-width at half-maximum bandwidth, which can lead to the design of the structure spectra by selection of different monodispersed particles. This work can find applications in fabrication of non-periodic nanostructures for functional surfaces, light extraction structures, and broadband nanophotonics.
ABSTRACT We compare the constraints from two (2019 and 2021) compilations of H ii starburst galaxy (H iiG) data and test the model independence of quasar (QSO) angular size data using six spatially flat and non-flat cosmological models. We find that the new 2021 compilation of H iiG data generally provides tighter constraints and prefers lower values of cosmological parameters than those from the 2019 H iiG data. QSO data by themselves give relatively model-independent constraints on the characteristic linear size, lm, of the QSOs within the sample. We also use Hubble parameter [H(z)], baryon acoustic oscillation (BAO), Pantheon Type Ia supernova (SN Ia) apparent magnitude (SN-Pantheon), and DES-3 yr binned SN Ia apparent magnitude (SN-DES) measurements to perform joint analyses with H iiG and QSO angular size data, since their constraints are not mutually inconsistent within the six cosmological models we study. A joint analysis of H(z), BAO, SN-Pantheon, SN-DES, QSO, and the newest compilation of H iiG data provides almost model-independent summary estimates of the Hubble constant, $H_0=69.7\pm 1.2\ \rm {km\,s^{-1}\,Mpc^{-1}}$, the non-relativistic matter density parameter, $\Omega _{\rm m_0}=0.293\pm 0.021$, and lm = 10.93 ± 0.25 pc.
Prediction and control of chemical mixing are vital for many scientific areas such as subsurface reactive transport, climate modeling, combustion, epidemiology, and pharmacology. Due to the complex nature of mixing in heterogeneous and anisotropic media, the mathematical models related to this phenomenon are not analytically tractable. Numerical simulations often provide a viable route to predict chemical mixing accurately. However, contemporary modeling approaches for mixing cannot utilize available spatial-temporal data to improve the accuracy of the future prediction and can be compute-intensive, especially when the spatial domain is large and for long-term temporal predictions. To address this knowledge gap, in this work we will present in this paper a deep learning (DL) modeling framework applied to predict the progress of chemical mixing under fast bimolecular reactions. This framework uses convolutional neural networks (CNN) for capturing spatial patterns and long short-term memory (LSTM) networks for forecasting temporal variations in mixing. By careful design of the framework—placement of non-negative constraint on the weights of the CNN and the selection of activation function, the framework ensures non-negativity of the chemical species at all spatial points and for all times. Our DL-based framework is fast, accurate, and requires minimal data for training. The time needed to obtain a forecast using the model is a fraction (≈ O(-6)) of the time needed to obtain the result using a high-fidelity simulation. To achieve an error of 10% (measured using the infinity norm) for capturing local-scale mixing features such as interfacial mixing, only 24% to 32% of the sequence data for model training is required. To achieve the same level of accuracy for capturing global-scale mixing features, the sequence data required for model training is 64% to 70% of the total spatial-temporal data. Hence, the proposed approach—a fast and accurate way to forecast long-time spatial-temporal mixing patterns in heterogeneous and anisotropic media—will be a valuable tool for modeling reactive-transport in a wide range of applications.
We report translating the surging interest in neuromorphic electronic components, such as those based on nonlinearities near Mott transitions, into large-scale commercial deployment faces steep challenges in the current lack of means to identify and design key material parameters. These issues are exemplified by the difficulties in connecting measurable material properties to device behavior via circuit element models. Here, the principle of local activity is used to build a model of VO 2 /SiN Mott threshold switches by sequentially accounting for constraints from a minimal set of quasistatic and dynamic electrical and high-spatial-resolution thermal data obtained via in situ thermoreflectance mapping. By combining independent data sets for devices with varying dimensions, the model is distilled to measurable material properties, and device scaling laws are established. The model can accurately predict electrical and thermal conductivities and capacitances and locally active dynamics (especially persistent spiking self-oscillations). The systematic procedure by which this model is developed has been a missing link in predictively connecting neuromorphic device behavior with their underlying material properties, and should enable rapid screening of material candidates before employing expensive manufacturing processes and testing procedures.
Abstract This work presents a novel numerical procedure for reconstructing volumetric density and velocity fields from planar laser-induced fluorescence (PLIF) and stereoscopic particle image velocimetry (SPIV) data. This new method is theoretically and practically demonstrated to provide more accurate 3D vortical structures and density fields in high shear flows than reconstruction methods based on the mean convective velocity. While Taylor’s hypothesis of frozen turbulence is commonly applied by using the local mean streamwise velocity, the proposed algorithm uses the measured local instantaneous velocity for data convection. It consists of a step-by-step reconstruction based on a mixed Lagrangian–Eulerian solver that includes the 3D interpolation of scattered flow data and that relaxes the Taylor’s hypothesis by iterative enforcement of the incompressibility constraint on the velocity field. This methodology provides 3D fields with temporal resolution, spatial resolution, and accuracy comparable to that of real 3D snapshots, thus providing a practical alternative to tomographic measurements. The procedure is validated using numerical data of the constant-density channel flow available on the Johns Hopkins University Turbulence Database (JHTDB), showing the accurate reconstruction of the 3D velocity field. The algorithm is applied to an experimental dataset of PLIF and SPIV measurements of a variable-density jet flow, demonstrating its capability to provide 3D velocity and density fields that are more consistent with the Navier–Stokes equations compared to the mean flow convective method. Graphic abstract
Here, we use updated Type Ia Pantheon+ supernova, baryon acoustic oscillation, and Hubble parameter (now also accounting for correlations) data, as well as new reverberation-measured C $\tiny{IV}$ quasar data, and quasar angular size, H $\tiny{II}$ starburst galaxy, reverberation-measured Mg $\tiny{II}$ quasar, and Amati correlated gamma-ray burst data to constrain cosmological parameters. We show that these data sets result in mutually consistent constraints and jointly use them to constrain cosmological parameters in six different spatially-flat and non-flat cosmological models. Our analysis provides summary model independent determinations of two key cosmological parameters: the Hubble constant, H 0 = 69.8 ± 1.3 km s -1 Mpc -1 , and the current non-relativistic matter density parameter, Ω m0 = 0.288 ± 0.017. Our summary error bars are 2.4 and 2.3 times those obtained using the flat ΛCDM model and Planck TT,TE,EE + lowE + lensing cosmic microwave background (CMB) anisotropy data. Our H 0 value is very consistent with that from the local expansion rate based on the Tip of the Red Giant Branch and Type Ia supernova (SN Ia) data, is 2σ lower than that from the local expansion rate based on Cepheid and SN Ia data, and is 2σ higher than that in the flat ΛCDM model based on Planck TT,TE,EE + lowE + lensing CMB data. Our data compilation shows at most mild evidence for non-flat spatial hypersurfaces, but more significant evidence for dark energy dynamics, 2σ or larger in the spatially-flat dynamical dark energy models we study.
We present a joint analysis of the cosmic microwave background (CMB) lensing power spectra measured from the Data Release 6 of the Atacama Cosmology Telescope (ACT) and Planck PR4, cross-correlations between the ACT and Planck lensing reconstruction and galaxy clustering from unWISE, and the unWISE clustering auto-spectrum. We obtain 1.5% constraints on the matter density fluctuations at late times parametrized by the best constrained parameter combination 𝑆$^{3x2pt}_{8}$ ≡ 𝜎 8 (Ω 𝑚 /0.3) 0.4 = 0.815 ± 0.012. The commonly used 𝑆 8 ≡ 𝜎 8 (Ω 𝑚 /0.3) 0.5 parameter is constrained to 𝑆 8 = 0.816 ± 0.015. In combination with baryon acoustic oscillation (BAO) measurements we find 𝜎 8 = 0.815 ± 0.012. We also present sound-horizon-independent estimates of the present day Hubble rate of 𝐻 0 = 66.4$^{+3.2}_{−3.7}$ km s −1 Mpc −1 from our large scale structure data alone and 𝐻 0 = 64.3$^{+2.1}_{−2.4}$ km s −1 Mpc −1 in combination with uncalibrated supernovae from Pantheon+. Using parametric estimates of the evolution of matter density fluctuations, we place constraints on cosmic structure in a range of high redshifts typically inaccessible with cross-correlation analyses. Combining lensing cross- and autocorrelations, we derive a 3.3% constraint on the integrated matter density fluctuations above 𝑧 = 2.4, one of the tightest constraints in this redshift range and fully consistent with a Λ cold dark matter (ΛCDM) model fit to the primary CMB from Planck. Finally, combining with primary CMB observations and using the extended low redshift coverage of these combined datasets we derive constraints on a variety of extensions to the ΛCDM model including massive neutrinos, spatial curvature, and dark energy. We find in flat ΛCDM ∑𝑚 𝜈 < 0.12 eV at 95% confidence using the large scale structure data, BAO measurements from Sloan Digital Sky Survey, and primary CMB observations.
African swine fever virus (ASFv) is a virulent pathogen that threatens domestic swine industries globally and persists in wild boar populations in some countries. Persistence in wild boar can challenge elimination and prevent disease-free status, making it necessary to address wild swine in proactive response plans. In the United States, invasive wild pigs are abundant and found across a wide range of ecological conditions that could drive different epidemiological dynamics among populations. Information on the size of the control areas required to rapidly eliminate the ASFv in wild pigs and how this area should change with management constraints and local ecology is needed to optimize response planning. We developed a spatially explicit disease transmission model contrasting wild pig movement and contact ecology in two ecosystems in Southeastern United States. We simulated ASFv spread and determined the optimal response area (reported as the radius of a circle) for eliminating ASFv rapidly over a range of detection times (when ASFv was detected relative to the true date of introduction), culling capacities (proportion of wild pigs in the culling zone removed weekly) and wild pig densities. Large radii for response areas (14 km) were needed under most conditions but could be shortened with early detection (≤ 8 weeks) and high culling capacities (≥ 15% weekly). Under most conditions, the ASFv was eliminated in less than 22 weeks using optimal control radii, although ecological conditions with high rates of wild pig movement required higher culling capacities (≥ 10% weekly) for elimination within 1 year. The results highlight the importance of adjusting response plans based on local ecology and show that wild pig movement is a better predictor of the optimal response area than the number of ASFv cases early in the outbreak trajectory. Furthermore, our framework provides a tool for determining optimal control plans in different areas, guiding expectations of response impacts, and planning resources needed for rapid elimination.
Abstract Unlike streamflow, which can be sampled in aggregate at the catchment outlet, evapotranspiration (ET) is spatially dispersed, challenging large‐scale age estimation. Here, we introduce an approach for constraining the age of ET via mass balance and present the minimum flux‐weighted age of ET across the continental United States using distributed, publicly available water flux data sets. The lower‐bound constraint on ET age can be calculated by assuming that ET is preferentially sourced from the most recent precipitation through a last‐in, first‐out algorithm. From 2012 to 2017, ET was at least several months old across large areas of the western continental United States, including in Mediterranean and (semi‐)arid climate zones and shrub and evergreen needleleaf plant communities. The primary limitation of this approach is that it provides only a minimum flux‐weighted average age to satisfy the mass balance of outgoing fluxes; true ET fluxes are composed of distributions of ages and may be composed of much older water. The primary advantage of the approach is that flux time series of precipitation and ET are sufficient to constrain ET age, and model parameterization is unnecessary. ET ages can be used to validate tracer‐aided and modeling approaches and inform studies of biogeochemistry, water‐rock interactions, and plant water sourcing under drought.
Abstract The delensing of cosmic microwave background (CMB) maps will be increasingly valuable for extracting as much information as possible from future CMB surveys. Delensing provides many general benefits, including sharpening of the acoustic peaks, more accurate recovery of the damping tail, and reduction of lensing-inducedB-mode power. In this paper we present several applications of delensing focused on testing theories of early-universe inflation with observations of the CMB. We find that delensing the CMB results in improved parameter constraints for reconstructing the spectrum of primordial curvature fluctuations, probing oscillatory features in the primordial curvature spectrum, measuring the spatial curvature of the universe, and constraining several different models of isocurvature perturbations. In some cases we find that delensing can recover almost all of the constraining power contained in unlensed spectra, and it will be a particularly valuable analysis technique to achieve further improvements in constraints for model parameters whose measurements are not expected to improve significantly when utilizing only lensed CMB maps from next-generation CMB surveys. We also quantify the prospects of testing the single-field inflation tensor consistency condition using delensed CMB data; we find it to be out of reach of current and proposed experimental technology and advocate for alternative detection methods.