Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regularity”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Spectral Threshold for Extremal Cyclic Edge-Connectivity

In this report, the cyclic edge-connectivity of a graph G is the least k such that there exists a set of k edges whose removal disconnects G into components where every component contains a cycle. We show that for graphs of minimum degree at least 3 and girth g at least 4, the cyclic edge-connectivity is bounded above by (Δ-2)g where Δ is the maximum degree. We then prove that if the second eigenvalue of the adjacency matrix of a d-regular graph of girth g ≥ 4 is sufficiently small, then the cyclic edge-connectivity is (d-2)g, providing a spectral condition for when this upper bound on cyclic edge-connectivity is tight.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Controls on the frequency content of near-source infrasound at open-vent volcanoes: a case study from Volcán Villarrica, Chile

Infrasound signals from open-vent volcanoes can contain information related to that volcano’s eruption dynamics and future activity. Studying a specific volcano’s acoustic signature during periods of relative quiescence is thus important for understanding potential warning signs of impending eruptions. Volcán Villarrica, located in southern Chile, has a long-lived active lava lake that produces continuous infrasound with spectral peaks centered around 1 Hz that vary by + / - ~ 0.2 Hz over day-to-week-long time scales. The infrasound frequency content has been shown to illuminate key volcanic properties such as eruption style and crater shape. Leading up to Villarrica’s most recent paroxysm in 2015, for instance, infrasound spectral changes coincided with a rise in Villarrica’s lava lake level. Quantifying and understanding the regular fluctuation in recorded infrasonic frequencies from Villarrica and other open-vent volcanoes is thus imperative during periods of relative calm. Our experiment entailed a week-long period of crater rim infrasound observations associated with stable, open-vent activity and revealed two independent source processes: spatter bursting events (relatively broadband infrasound) and lava lake-induced tremor (sustained signals peaked at 1 Hz). A comparison of these acoustic signals with results from 3-D finite-difference time-domain wave propagation models shows that while sound speed variability can influence Villarrica’s spectrum, it cannot explain the full extent of the observed frequency excursions. We conclude that source spectrum variability is primarily responsible for the frequency excursions. In conclusion, our work highlights the utility of data collected from a dense distribution of twenty infrasound sensors operating at the summit coupled with numerical modeling of sound radiation.

58 GEOSCIENCES↗
DOE OSTI2025

On the Sampling-Based Computation of Nash Equilibria Under Uncertainty via the Nikaido–Isoda Function

We consider the computation of an equilibrium of a stochastic Nash equilibrium problem, where the player objectives are assumed to be L 0 -Lipschitz continuous and convex, given rival decisions with convex and closed player-specific feasibility sets. To address this problem, we consider minimizing a suitably defined value function defined using the Nikaido–Isoda function. Such an avenue does not necessitate either monotonicity properties of the concatenated gradient map or potentiality requirements on the game but does require a suitable regularity requirement under which a stationary point is a Nash equilibrium. We design and analyze a sampling-enabled projected-gradient-response method, reliant on inexact resolution of a player-level best-response subproblem. Here, by deriving suitable Lipschitzian guarantees on the value function, we derive both asymptotic guarantees for the sequence of generated iterates as well as rate and complexity guarantees for computing a stationary point by appropriate choices of the sampling rate and inexactness sequence.

Nikaido-Isoda function↗
DOE OSTI2021

A primal–dual algorithm for risk minimization

In this paper, we develop an algorithm to efficiently solve risk-averse optimization problems posed in reflexive Banach space. Such problems often arise in many practical applications as, e.g., optimization problems constrained by partial differential equations with uncertain inputs. Unfortunately, for many popular risk models including the coherent risk measures, the resulting risk-averse objective function is nonsmooth. Here, this lack of differentiability complicates the numerical approximation of the objective function as well as the numerical solution of the optimization problem. To address these challenges, we propose a primal–dual algorithm for solving large-scale nonsmooth risk-averse optimization problems. This algorithm is motivated by the classical method of multipliers and by epigraphical regularization of risk measures. As a result, the algorithm solves a sequence of smooth optimization problems using derivative-based methods. We prove convergence of the algorithm even when the subproblems are solved inexactly and conclude with numerical examples demonstrating the efficiency of our method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A proximal trust-region method for nonsmooth optimization with inexact function and gradient evaluations

Many applications require minimizing the sum of smooth and nonsmooth functions. For example, basis pursuit denoising problems in data science require minimizing a measure of data misfit plus an $\ell^1$-regularizer. Similar problems arise in the optimal control of partial differential equations (PDEs) when sparsity of the control is desired. Here, we develop a novel trust-region method to minimize the sum of a smooth nonconvex function and a nonsmooth convex function. Our method is unique in that it permits and systematically controls the use of inexact objective function and derivative evaluations. When using a quadratic Taylor model for the trust-region subproblem, our algorithm is an inexact, matrix-free proximal Newton-type method that permits indefinite Hessians. We prove global convergence of our method in Hilbert space and demonstrate its efficacy on three examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Differential effects of redox conditions on the decomposition of litter and soil organic matter

Soil redox conditions exert substantial influence on biogeochemical processes in terrestrial ecosystems. Humid tropical forests are often characterized by fluctuating redox dynamics, yet how redox fluctuations affect patterns in soil versus litter decomposition and associated CO2 fluxes are not well understood in these ecosystems. We used a 13C-labeled litter addition to explicitly follow the decomposition of litter-derived vs. native soil-derived organic matter in response to four different soil redox regimes—static oxic or anoxic, and two oscillating treatments—in soil from the Luquillo Experimental Forest, Puerto Rico. We coupled this incubation experiment with high-resolution mass spectrometry analysis to characterize the preferential decomposition of specific classes of organic molecules. CO2 production from litter and soil organic matter (SOM) showed distinctly different responses to redox manipulation. The cumulative production of SOM-derived CO2 was positively correlated with the length of soil exposure to an oxic headspace (r = 0.89, n = 20), whereas cumulative 13C-litter-derived CO2 production was not linked to oxygen availability. The CO2 production rate from litter was highest under static anoxic conditions in the first half of the incubation period, and later dropped to the lowest among all redox treatments. In the consistently anoxic soils, we observed the depletion of more oxidized water-extractable organic matter (especially amino sugars, carbohydrates, and proteins) over time, suggesting that under anaerobic conditions, microbes preferentially used more oxidized litter-derived compounds during the early stages of decomposition. Results from kinetic modeling showed that more frequent anoxic exposure limited the decomposition of a slow-cycling C pool, but not a fast-cycling pool. Overall, our results demonstrate that substrate source—freshly added litter vs. native organic matter—plays an important role in the redox sensitivity of organic matter decomposition. In soil environments that regularly experience redox fluctuations, anaerobic heterotrophs can be surprisingly effective in degrading fresh plant litter.

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

Stochastic projective splitting

Here, we present a new, stochastic variant of the projective splitting (PS) family of algorithms for inclusion problems involving the sum of any finite number of maximal monotone operators. This new variant uses a stochastic oracle to evaluate one of the operators, which is assumed to be Lipschitz continuous, and (deterministic) resolvents to process the remaining operators. Our proposal is the first version of PS with such stochastic capabilities. We envision the primary application being machine learning (ML) problems, with the method’s stochastic features facilitating “mini-batch” sampling of datasets. Since it uses a monotone operator formulation, the method can handle not only Lipschitz-smooth loss minimization, but also min–max and noncooperative game formulations, with better convergence properties than the gradient descent-ascent methods commonly applied in such settings. The proposed method can handle any number of constraints and nonsmooth regularizers via projection and proximal operators. We prove almost-sure convergence of the iterates to a solution and a convergence rate result for the expected residual, and close with numerical experiments on a distributionally robust sparse logistic regression problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Thermodynamic stabilization of nanocrystalline aluminum

Nanocrystalline metals are generally unstable due to a large volume fraction of high-energy grain boundaries associated with a small grain size. Preferential dopant segregation to the high-energy grain boundaries is observed to enhance the stability of the material’s microstructure by minimizing its energy. Nanocrystalline aluminum-dopant systems were evaluated for thermodynamic stability against grain growth and phase precipitation via the mechanism of grain boundary segregation according to a modified regular nanocrystalline solution model. Fifty-one potential dopant elements have been evaluated for their efficacy in stabilizing nanostructures with three potential candidates, magnesium, lanthanum, and silicon, identified possessing the characteristics to promote grain boundary segregation and a state of thermodynamic stability in aluminum’s nanocrystalline regime. Here, the minimum dopant content required to achieve nanocrystalline microstructure stability is identified for each of the three candidate elements. Beyond this minimum content, further addition of the dopant elements decreased the final microstructure’s stability with no effects on the existence of a stable nanocrystalline state.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Haar-Like Wavelets on Hierarchical Trees

Here, discrete wavelet methods, originally formulated in the setting of regularly sampled signals, can be adapted to data defined on a point cloud if some multiresolution structure is imposed on the cloud. A wide variety of hierarchical clustering algorithms can be used for this purpose, and the multiresolution structure obtained can be encoded by a hierarchical tree of subsets of the cloud. Prior work introduced the use of Haar-like bases defined with respect to such trees for approximation and learning tasks on unstructured data. This paper builds on that work in two directions. First, we present an algorithm for constructing Haar-like bases on general discrete hierarchical trees. Second, with an eye towards data compression, we present thresholding techniques for data defined on a point cloud with error controlled in the $L$ $\infty$ norm and in a Hölder-type norm. In a concluding trio of numerical examples, we apply our methods to compress a point cloud dataset, study the tightness of the $L$ $\infty$ error bound, and use thresholding to identify MNIST classifiers with good generalizability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

On the Analysis of the Tin-Inside-H 3 S Mössbauer Experiment

Here, a simple analysis is presented of the particular experiment used to prove the bulk nature of very-high-Tc superconductivity in H 3 S compound under ultra-high pressure. In the experiment, an internal magnetic field was sensed by the synchrotron Mössbauer spectroscopy in tin placed inside the H 3 S sample. The experiment showed peculiar anisotropy with respect to the direction of the applied field at first sight. By considering actual experimental geometries and parameters of the experiment, we show that this particular observation is consistent with the expectations for a regular type-II superconductor with Meissner expulsion and pinning.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

On the closedness and geometry of tensor network state sets

Tensor network states (TNS) are a powerful approach for the study of strongly correlated quantum matter. The curse of dimensionality is addressed by parametrizing the many-body state in terms of a network of partially contracted tensors. These tensors form a substantially reduced set of effective degrees of freedom. In practical algorithms, functionals like energy expectation values or overlaps are optimized over certain sets of TNS. Concerning algorithmic stability, it is important whether the considered sets are closed because, otherwise, the algorithms may approach a boundary point that is outside the TNS set and tensor elements diverge. Here we discuss the closedness and geometries of TNS sets, and we propose regularizations for optimization problems on non-closed TNS sets. We show that sets of matrix product states (MPS) with open boundary conditions, tree tensor network states, and the multiscale entanglement renormalization ansatz are always closed, whereas sets of translation-invariant MPS with periodic boundary conditions (PBC), heterogeneous MPS with PBC, and projected entangled pair states are generally not closed. The latter is done using explicit examples like the W state, states that we call two-domain states, and fine-grained versions thereof.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Residual Stress Measurements in Extreme Environments for Hazardous, Layered Specimens

In nuclear fuel plates of low-enriched U-10Mo (LEU) clad with aluminum by hot isostatic pressing (HIP), post-irradiation stresses arising during reactor shutdown are a major concern for safe reactor operations. Measurement of those residual stresses has not previously been possible because the high radioactivity of the plates requires handling only by remote manipulation in a hot cell. The incremental slitting method for measuring through-thickness stress profiles was modified, and a system for automated, remote operation was built and tested. Experimental modifications consisted of replacing electric-discharge machining (EDM) with a small end mill and strain-gauge measurements with cantilever displacement measurements. The inverse method used to calculate stresses was the pulse-regularization method modified to allow discontinuities across material interfaces. The new system was validated by comparing with conventional slitting on a depleted U-10Mo (DU) fuel plate. The new system was applied to two measurements each on six as-fabricated (pre-irradiation) LEU miniature fuel plates. Variations between the measurements at two locations in the same plate were strongly correlated with measured geometrical heterogeneity in the plate—a tilt in the fuel foil. Compressive stresses in the U-10Mo were shown to increase from 20 to 250 MPa as the ratio of aluminum thickness to U-10Mo thickness increased causing increased constraint during cooling. Faster cooling rates during processing also increased stress magnitudes. The measurements trends agreed with data in the literature from similar plates made with DU, which further validates the method. Because other methods are impractical in a hot cell, the modified slitting method is now poised for the first measurements of post-irradiation stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Combining shallow-water and analytical wake models for tidal array micro-siting

For tidal-stream energy to become a competitive renewable energy source, clustering multiple turbines into arrays is paramount. Array optimisation is thus critical for achieving maximum power performance and reducing cost of energy. However, ascertaining an optimal array layout is a complex problem, subject to specific site hydrodynamics and multiple inter-disciplinary constraints. In this work, we present a novel optimisation approach that combines an analytical-based wake model, FLORIS, with an ocean model, Thetis. The approach is demonstrated through applications of increasing complexity. By utilising the method of analytical wake superposition, the addition or alteration of turbine position does not require re-calculation of the entire flow field, thus allowing the use of simple heuristic techniques to perform optimisation at a fraction of the computational cost of more sophisticated methods. Using a custom condition-based placement algorithm, this methodology is applied to the Pentland Firth for arrays with turbines of 3.05 m/s rated speed, demonstrating practical implications whilst considering the temporal variability of the tide. For a 24-turbine array case, micro-siting using this technique delivered an array 15.8% more productive on average than a staggered layout, despite flow speeds regularly exceeding the rated value. Performance was evaluated through assessment of the optimised layout within the ocean model that treats turbines through a discrete turbine representation. Used iteratively, this methodology could deliver improved array configurations in a manner that accounts for local hydrodynamic effects.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Effect of Heave Plate Hydrodynamic Force Parameterization on a Two-Body Wave Energy Converter

Heave plates are one approach to generating the reaction force necessary to harvest energy from ocean waves. In a Morison equation description of the hydrodynamic force, the components of drag and added mass depend primarily on the heave plate oscillation. These terms may be parameterized in three ways: (1) as a single coefficient invariant across sea state, most accurate at the reference sea state, (2) coefficients dependent on the oscillation amplitude, but invariant in phase, that are most accurate for relatively small amplitude motions, and (3) coefficients dependent on both oscillation amplitude and phase, which are accurate for all oscillation amplitudes. We validate a MATLAB model for a two-body point absorber wave energy converter against field data and a dynamical model constructed in ProteusDS. We then use the MATLAB model to evaluate the effect of these parameterizations on estimates of heave plate motion, tension between the float and heave plate, and wave energy converter electrical power output. We find that power predictions using amplitude-dependent coefficients differ by up to 30% from models using invariant coefficients for regular waves ranging in height from 0.5 to 1.9 m. Amplitude- and phase-dependent coefficients, however, yield less than a 5% change when compared with coefficients dependent on amplitude only. This suggests that amplitude-dependent coefficients can be important for accurate wave energy converter modeling, but the added complexity of phase-dependent coefficients yields little further benefit. We show similar, though less pronounced, trends in maximum tether tension, but note that heave plate motion has only a weak dependence on coefficient fidelity. Finally, we emphasize the importance of using experimentally derived added mass over that calculated from boundary element methods, which can lead to substantial under-prediction of power output and peak tether tension.

dynamical model↗
DOE OSTI2024

Recent developments in mathematical aspects of relativistic fluids

Abstract We review some recent developments in mathematical aspects of relativistic fluids. The goal is to provide a quick entry point to some research topics of current interest that is accessible to graduate students and researchers from adjacent fields, as well as to researches working on broader aspects of relativistic fluid dynamics interested in its mathematical formalism. Instead of complete proofs, which can be found in the published literature, here we focus on the proofs’ main ideas and key concepts. After an introduction to the relativistic Euler equations, we cover the following topics: a new wave-transport formulation of the relativistic Euler equations tailored to applications; the problem of shock formation for relativistic Euler; rough (i.e., low-regularity) solutions to the relativistic Euler equations; the relativistic Euler equations with a physical vacuum boundary; relativistic fluids with viscosity. We finish with a discussion of open problems and future directions of research.

Disconzi, Marcelo (ORCID:0000000234497778)↗
DOE OSTI2021

Seasonal Disorder in Urban Traffic Patterns: A Low Rank Analysis

This article proposes several advances to sparse nonnegative matrix factorization (SNMF) as a way to identify large-scale patterns in urban traffic data. The input to our model is traffic counts organized by time and location. Nonnegative matrix factorization additively decomposes this information, organized as a matrix, into a linear sum of temporal signatures. Penalty terms encourage this factorization to concentrate on only a few temporal signatures, with weights which are not too large. Our interest here is to quantify and compare the regularity of traffic behavior, particularly across different broad temporal windows. In addition to the rank and error, we adapt a measure introduced by Hoyer to quantify sparsity in the representation. Combining these, we construct several curves which quantify error as a function of rank (the number of possible signatures) and sparsity; as rank goes up and sparsity goes down, the approximation can be better and the error should decreases. Plots of several such curves corresponding to different time windows leads to a way to compare disorder/order at different time scalewindows. In this paper, we apply our algorithms and procedures to study a taxi traffic dataset from New York City. In this dataset, we find weekly periodicity in the signatures, which allows us an extra framework for identifying outliers as significant deviations from weekly medians. We then apply our seasonal disorder analysis to the New York City traffic data and seasonal (spring, summer, winter, fall) time windows. We do find seasonal differences in traffic order.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Multifunctional Fibroblasts Enhanced via Thermal and Freeze-Drying Post-treatments of Aligned Electrospun Nanofiber Membranes

Parallel fibrous scaffolds play a critical role in controlling the morphology of cells to be more natural and biologically inspired. Among popular tissue engineering materials, poly(2-hydroxyethyl methacrylate) (pHEMA) has been widely investigated in conventional forms due to its biocompatibility, low toxicity, and hydrophilicity. However, the swelling of pHEMA in water remains a major concern. To address this issue, randomly oriented and aligned as-spun pHEMA nanofibrous scaffolds were first fabricated at speeds of 300 and 2000 rpm in this study, which were then post-treated using either a thermal or a freeze-drying method. In cell assays, human dermal fibroblasts (HDFs) adhered to the freeze-drying treated substrates at a significantly faster rate, whereas they had a higher cell growth rate on thermally-treated substrates. Results indicated that the structural properties of pHEMA nanofibrous scaffolds and subsequent cellular behaviors were largely dependent on post-treatment methods. Moreover, this study suggests that aligned pHEMA nanofibrous substrates tended to induce regular fibroblast orientation and unidirectionally oriented actin cytoskeletons over random pHEMA nanofibrous substrates. Such information has predictive power and provides insights into promising post-treatment methods for improving the properties of aligned pHEMA scaffolds for numerous tissue engineering applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Chapter 7 - CRISPR-Based Tools for Microbial Cell Factories

The development of microbial chassis for the production of a variety of biochemicals and biofuels is a growing area of research. How to efficiently manipulate genetic information to achieve optimal production of these compounds is a key area of focus in the field. In recent years, clustered regularly interspaced palindromic repeats (CRISPR) and its associated proteins (Cas) have become a popular strategy for gene editing and regulation in many organisms due to its versatility and efficacy. Here, we describe methods developed utilizing CRISPR-Cas systems for engineering microbial cell factories (e.g., bacteria and yeast) at the single gene to genome scale for mutagenesis and transcriptional regulation of target genes. Finally, we provide a perspective on the challenges and opportunities for the applications of advanced CRISPR-Cas-based tools for engineering microbial cell factories.

BIOMASS FUELS↗