Search NASA⌕ Search

SEARCH · Search NASA

Results for “Total variation decomposition”

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.

Domain Decomposition for Integer Optimal Control with Total Variation Regularization

Total variation integer optimal control problems admit solutions and necessary optimality conditions via geometric variational analysis. In spite of the existence of said solutions, algorithms which solve the discretized objective suffer from high numerical cost associated with the combinatorial nature of integer programming. Hence, such methods are often limited to small and medium-sized problems. We propose a globally convergent, coordinate descent–inspired algorithm that allows tractable subproblem solutions restricted to a partition of the domain. Our decomposition method solves relatively small trust-region subproblems that modify the control variable on a subdomain only. Given nontrivial subdomain overlap, we prove that a global first-order necessary optimality condition is equivalent to a first-order necessary optimality condition per subdomain. We additionally show that a sufficient decrease is achieved on a single subdomain by way of a trust-region subproblem solver using geometric measure–theoretic arguments, which we integrate with a greedy patch selection to prove convergence of our algorithm. In conclusion, we demonstrate the practicality of our algorithm on a benchmark large-scale, PDE-constrained integer optimal control problem and find that our method is faster than the state of the art.

domain decomposition↗

Facile alignment estimation in carbon nanotube films using image processing

Whether a macroscopic assembly of carbon nanotubes can exhibit the one-dimensional properties expected from individual nanotubes critically depends on how well the nanotubes are aligned inside the assembly. Therefore, a simple and accurate method for assessing the degree of alignment is desired for the rapid characterization of carbon nanotube films and fibers. Here, we present an end-to-end solution for determining the global and local spatial orientation of carbon nanotubes in films within a short amount of time using a fast, precise, and economical approach based on an image processing method applied to scanning electron microscopy images. Further, we first use Laplacian edge enhancement filtering for improving the appearance of edge regions, which is followed by image partitioning into multiple blocks to capture the nanoscale orientation characteristics and total variation-based image decomposition of these image blocks. We then perform a 2D-fast Fourier transform on the image decomposed textural components of these edge-enhanced image blocks to determine the orientation distribution, which is utilized to estimate the 2D nematic order parameter. To show the effectiveness of our method, we corroborated our results against results obtained with other state-of-the-art image processing and experimental techniques.

2D-FFT↗

Experimental quantum learning of a spectral decomposition

Currently available quantum hardware allows for small-scale implementations of quantum machine learning algorithms. Such experiments aid the search for applications of quantum computers by benchmarking the near-term feasibility of candidate algorithms. Here we demonstrate the quantum learning of a two-qubit unitary by a sequence of three parameterized quantum circuits containing a total of 21 variational parameters. Moreover, we variationally diagonalize the unitary to learn its spectral decomposition, i.e., its eigenvalues and eigenvectors. We illustrate how this can be used as a subroutine to compress the depth of dynamical quantum simulations. One can view our implementation as a demonstration of entanglement-enhanced machine learning, as only a single (entangled) training data pair is required to learn a 4 × 4 unitary matrix.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Point-scale organic-matter decomposition in streambeds is weakly associated with reach-scale respiration

Stream and river ecosystems play a central role in the movement and decomposition of particulate organic matter, serving as a conduit between terrestrial hillslopes and coastal environments. Microbial-catalyzed decomposition generates simpler organic molecules that fuel respiration, often in the sediments of these ecosystems. However, the degree of connection between sediment-associated respiration (ER??d) and organic-matter decomposition remains poorly understood. How that relationship compares to decomposition’s relationship with whole ecosystem (ERtot) and water column (ERwc) respiration is also not clear. We examined the link between particulate organic matter decomposition—using cellulose-based cotton strips as a standardized substrate—and all three components of respiration across 48 sites in the environmentally diverse Yakima River Basin (Washington State, USA). We hypothesized that decomposition within sediments would be most strongly related to ERsed, but decomposition rates were more closely associated with ERtot, with little connection to ERsed or ERwc. This suggests that particulate organic matter decomposition within stream/river sediments reflects integrated system respiration rather than processes confined to sediments or the water column alone. Further, across the basin, decomposition rates nearly spanned the previously reported global range for streams and rivers and were best explained by total dissolved nitrogen (TDN), sediment grain size, and aridity of the upstream drainage area. These results highlight the strong influence of land cover and basin-scale biophysical variation on sediment-associated decomposition processes and indicate that mechanistic models of organic matter decomposition in streams/rivers should account for coupled sediment–water–land interactions.

Stegen, James C. (ORCID:0000000191357424)↗

Fungal and bacterial growth variation due to drought and nitrogen addition experimental treatments. Loma Ridge Experimental Project. 2010-2012

Terrestrial ecosystem models assume that microbial communities respond instantaneously, or are immediately resilient, to environmental change. Here we tested this assumption by quantifying the resilience of a leaf litter community to changes in precipitation or nitrogen availability. By manipulating composition within a global change experiment, we decoupled the legacies of abiotic parameters versus that of the microbial community itself. After one rainy season, more variation in fungal composition could be explained by the original microbial inoculum than the litterbag environment (18% versus 5.5% of total variation). This compositional legacy persisted for 3 years, when 6% of the variability in fungal composition was still explained by the microbial origin. In contrast, bacterial composition was generally more resilient than fungal composition. Microbial functioning (measured as decomposition rate) was not immediately resilient to the global change manipulations; decomposition depended on both the contemporary environment and rainfall the year prior. Finally, using metagenomic sequencing, we showed that changes in precipitation, but not nitrogen availability, altered the potential for bacterial carbohydrate degradation, suggesting why the functional consequences of the two experiments may have differed. Predictions of how terrestrial ecosystem processes respond to environmental change may thus be improved by considering the legacies of microbial communities. This data package includes ten csv files (five data files and their corresponding data dictionaries) and one file-level metadata excel file. Data files contains information about which plots were exposed to treatments related to drought and nitrogen, information about litter bags reciprocal transplants manipulation for water input and nitrogen, detail information about water addition, precipitation records, and litter variables collected. Data dictionary files include detail explanation for each column in the data files. The file-level metadata file describes each file mentioned above. All the analyses were done using the R software.

54 ENVIRONMENTAL SCIENCES↗

Disruption event characterization and forecasting in tokamaks

Disruption prediction and avoidance is a critical need for next-step tokamaks, such as ITER. Disruption Event Characterization and Forecasting (DECAF) research fully automates analysis of tokamak data to determine chains of events that lead to disruptions and to forecast their evolution allowing sufficient time for mitigation or complete avoidance of the disruption. Disruption event chains related to local rotating or global magnetohydrodynamic (MHD) modes and vertical instability are examined with warnings issued for many off-normal physics events, including density limits, plasma dynamics, confinement transitions, and profile variations. Along with Greenwald density limit evaluation, a local radiative island power balance theory is evaluated and compared to the observation of island growth. Automated decomposition and analysis of rotating tearing modes produce physical event chains leading to disruptions. A total MHD state warning model comprised of 15 separate criteria produces a disruption forecast about 180 ms before a standard locked mode detector warning. Single DECAF event analyses have begun on KSTAR, MAST, and NSTX/-U databases with thousands of shot seconds of device operation using from 0.5 to 1 × 10 6 tested sample times per device. An initial multi-device database comparison illustrates a highly important result that plasma disruptivity does not need to increase as βN increases. Global MHD instabilities, such as resistive wall modes (RWMs), can give the briefest time period of warning before disruption compared to other physics events. In an NSTX database with unstable RWMs, the mode onset, loss of boundary and current control, and disruption event warnings are found in all cases and vertical displacement events are found in 91% of cases. An initial time-dependent reduced physics model of kinetic RWM stabilization created to forecast the disruption chain predicts instability 84% of the time for experimentally unstable cases with a relatively low false positive rate. Instances of the disruption event chain analysis illustrate dynamics including H–L back transitions for rotating MHD and global RWM triggering events. Furthermore, disruption warnings are issued with sufficient time before the disruption (on transport timescales) to potentially allow active profile control for disruption avoidance, active mode control, or mitigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulating measurable ecosystem carbon and nitrogen dynamics with the mechanistically defined MEMS 2.0 model

Abstract. For decades, predominant soil biogeochemical models have used conceptual soil organic matter (SOM) pools and only simulated them to a shallow depth in soil. Efforts to overcome these limitations have prompted the development of the new generation SOM models, including MEMS 1.0, which represents measurable biophysical SOM fractions, over the entire root zone, and embodies recent understanding of the processes that govern SOM dynamics. Here we present the result of continued development of the MEMS model, version 2.0. MEMS 2.0 is a full ecosystem model with modules simulating plant growth with above- and belowground inputs, soil water and temperature by layer, decomposition of plant inputs and SOM, and mineralization and immobilization of nitrogen (N). The model simulates two commonly measured SOM pools – particulate and mineral-associated organic matter (POM and MAOM, respectively). We present results of calibration and validation of the model with several grassland sites in the US. MEMS 2.0 generally captured the soil carbon (C) stocks (R2 of 0.89 and 0.6 for calibration and validation, respectively) and their distributions between POM and MAOM throughout the entire soil profile. The simulated soil N matches measurements but with lower accuracy (R2 of 0.73 and 0.31 for calibration and validation of total N in SOM, respectively) than for soil C. Simulated soil water and temperature were compared with measurements, and the accuracy is comparable to the other commonly used models. The seasonal variation in gross primary production (GPP; R2 = 0.83), ecosystem respiration (ER; R2 = 0.89), net ecosystem exchange (NEE; R2 = 0.67), and evapotranspiration (ET; R2 = 0.71) was well captured by the model. We will further develop the model to represent forest and agricultural systems and improve it to incorporate new understanding of SOM decomposition.

54 ENVIRONMENTAL SCIENCES↗

On the Existence of Fast Modes in Compressible Magnetohydrodynamic Turbulence

Abstract We study the existence and properties of fast magnetosonic modes in 3D compressible MHD turbulence by carrying out a number of simulations with compressible and incompressible driving conditions. We use two approaches to determine the presence of fast modes: mode decomposition based on spatial variations only and spatio-temporal 4D fast Fourier transform (4D FFT) analysis of all fluctuations. The latter method enables us to quantify fluctuations that satisfy the dispersion relation of fast modes with finite frequency. Overall, we find that the fraction of fast modes identified via the spatio-temporal 4D FFT approach in total fluctuation power is either tiny with nearly incompressible driving or ∼2% with highly compressible driving. We discuss the implications of our results for understanding the compressible fluctuations in space and astrophysical plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

Ecohydrological controls on root and microbial respiration in the East River watershed of Colorado

The main objective of this project was to conduct exploratory work to quantify how snow and rain water inputs influence the CO2 coming from the soil surface (soil CO2 flux), and its plant and microbial sources, in the East River watershed, near Crested Butte, Colorado. Knowledge gained from this effort laid the groundwork for a more comprehensive (ongoing) follow-on grant that is using a combination of experiments, field observations, machine learning and simulation modeling to fully disentangle these relationships. New field measurements were made at four locations along an elevational transect on Snodgrass Mountain at Rocky Mountain Biological Laboratory. These sites were chosen to differ in total snowpack and water table depth, and included both deciduous (aspen) and evergreen (spruce/fir) forest types. We used automated measurements of soil CO2 concentrations to quantify the total soil CO2 flux, and the vertical CO2 production within the soil profile at each site. Isotope (radiocarbon, 14C) measurements determined how much of the CO2 emitted from the soil surface came from plant respiration (root metabolism) versus microbial respiration (decomposition of soil organic matter) sources. Supporting data on plant phenology and microbial dynamics provided context for the observed variation in respiration sources. This work was motivated by our overarching hypothesis that quantifying belowground plant and microbial processes separately, and how they are influenced by snow and rain inputs, is necessary for understanding and predicting how the East River watershed ecosystems will respond to future environmental change.

54 ENVIRONMENTAL SCIENCES↗

Towards robust surrogate models: Benchmarking machine learning approaches to expediting phase field simulations of brittle fracture

Data-driven approaches have the potential to make modeling complex, nonlinear physical phenomena significantly more computationally tractable. For example, computational modeling of fracture is a core challenge where machine learning techniques have the potential to provide a much needed speedup that would enable progress in areas such as multi-scale modeling and uncertainty quantification. Currently, phase field modeling (PFM) of fracture is one such approach that offers a convenient variational formulation to model crack nucleation, branching and propagation. To date, machine learning techniques have shown promise in approximating PFM simulations. While standard fracture benchmarks represent realistic scenarios frequently observed in practice, they typically do not provide sufficiently challenging tests for data-driven methods. Here, to address this gap, we introduce a challenging dataset based on PFM simulations designed to benchmark and advance ML methods for fracture modeling. This dataset includes three energy decomposition methods, two boundary conditions, and 1000 random initial crack configurations for a total of 6000 simulations. Each sample contains 100 time steps capturing the temporal evolution of the crack field. Alongside this dataset, we also implement and evaluate Physics Informed Neural Networks (PINN), Fourier Neural Operators (FNO), and UNet models as baselines, and explore the impact of ensembling strategies on prediction accuracy. With this combination of our dataset and baseline models drawn from the literature we aim to provide a standardized and challenging benchmark for evaluating machine learning approaches to solid mechanics. Our results highlight both the promise and limitations of popular current models, and demonstrate the utility of this dataset as a testbed for advancing machine learning in fracture mechanics research.

Benchmark dataset↗

Multiwavelength Bulge–Disk Decomposition for the Galaxy M81 (NGC 3031). I. Morphology

A panchromatic investigation of morphology for the early-type spiral galaxy M81 is presented in this paper. We perform bulge–disk decomposition in M81 images at a total of 20 wave bands from far-UV to near-IR (NIR) obtained with GALEX, Swift, Sloan Digital Sky Survey, WIYN, Two Micron All Sky Survey, Wide-field Infrared Survey Explorer, and Spitzer. Morphological parameters such as Sérsic index, effective radius, position angle, and axis ratio for the bulge and the disk are thus derived at all of the wave bands, which enables quantifying the morphological K -correction for M81 and makes it possible to reproduce images for the bulge and the disk in the galaxy at any wave band. The morphology as a function of wavelength appears as a variable-slope trend of the Sérsic index and the effective radius, in which the variations are steep at UV–optical and shallow at optical–NIR bands; the position angle and the axis ratio keep invariable at least at optical–NIR bands. It is worth noting that the Sérsic index for the bulge reaches ~4–5 at optical and NIR bands, but drops to ~1 at UV bands. This difference brings forward a caveat that a classical bulge is likely misidentified for a pseudobulge or no bulge at high redshifts where galaxies are observed through rest-frame UV channels with optical telescopes. The next work of this series is planned to study spatially resolved spectral energy distributions for the bulge and the disk, respectively, and thereby explore stellar population properties and star formation/quenching history for the galaxy composed of the subsystems.

79 ASTRONOMY AND ASTROPHYSICS↗

The Transcriptional Response of Soil Bacteria to Long-Term Warming and Short-Term Seasonal Fluctuations in a Terrestrial Forest

Terrestrial ecosystems are an important carbon store, and this carbon is vulnerable to microbial degradation with climate warming. After 30 years of experimental warming, carbon stocks in a temperate mixed deciduous forest were observed to be reduced by 30% in the heated plots relative to the controls. In addition, soil respiration was seasonal, as was the warming treatment effect. We therefore hypothesized that long-term warming will have higher expressions of genes related to carbohydrate and lipid metabolism due to increased utilization of recalcitrant carbon pools compared to controls. Because of the seasonal effect of soil respiration and the warming treatment, we further hypothesized that these patterns will be seasonal. We used RNA sequencing to show how the microbial community responds to long-term warming (~30 years) in Harvard Forest, MA. Total RNA was extracted from mineral and organic soil types from two treatment plots (+5°C heated and ambient control), at two time points (June and October) and sequenced using Illumina NextSeq technology. Treatment had a larger effect size on KEGG annotated transcripts than on CAZymes, while soil types more strongly affected CAZymes than KEGG annotated transcripts, though effect sizes overall were small. Although, warming showed a small effect on overall CAZymes expression, several carbohydrate-associated enzymes showed increased expression in heated soils (~68% of all differentially expressed transcripts). Further, exploratory analysis using an unconstrained method showed increased abundances of enzymes related to polysaccharide and lipid metabolism and decomposition in heated soils. Compared to long-term warming, we detected a relatively small effect of seasonal variation on community gene expression. Together, these results indicate that the higher carbohydrate degrading potential of bacteria in heated plots can possibly accelerate a self-reinforcing carbon cycle-temperature feedback in a warming climate.

54 ENVIRONMENTAL SCIENCES↗