Engineering Papers⌕ Search

SEARCH · Engineering Papers

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↗

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)↗

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↗

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↗