Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix factorization”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of Fe atoms on hardening of a nickel matrix: Nanoindentation experiments and atom-scale numerical modeling

Significant hardening effect due to Fe concentrations in Ni-based alloys with face-centered-cubic structure has been studied by using a combined experimental and atomistic-based computational approach via nanoindentation tests. The obtained experimental load–displacement data for the [0 0 1] crystal orientation reached a qualitative good agreement with molecular dynamics simulations results, leading to strong evidence that the main strengthening factors are associated to sluggish dislocation diffusion, reduced defect sizes and the nucleation of tetrahedral stacking faults. Here, interstitial type prismatic dislocation loops mainly formed by $\frac{1}{6}$$\langle$1 1 2$\rangle$ Shockley dislocations are nucleated during the loading process, where their interaction leads to the formation of pyramidal shaped stacking fault which are mainly created by $\frac{1}{3}$$\langle$1 0 0$\rangle$ Hirth dislocations lines. Observing both types of defects coexisting in the same plastic deformation zone by both approaches. Reported mechanical data, measured experimentally and interpreted numerically, are also in accordance with microstructural SEM and TEM investigations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamic Control of Photocatalytic Proton Reduction through the Mechanical Actuation of a Hydrogel Host Matrix

This paper describes a photocatalytic hydrogen evolution system that is dynamically and reversibly responsive to the pH of the surrounding solution through actuation of a micro-hydrogel (microgel) matrix that hosts the photocatalysts (CdSe/CdS nanorods). Here, the reversible actuation occurs within 0.7 s (swelling) and 1.7 s (contraction). A ΔpH of 0.01 relative to the pK a of the tertiary amine on the microgel polymer (7.27) results in a reversible change in the average diameter of the microgel hosts by a factor of 2.5 and a change in the photocatalytic turnover frequency (TOF) by a factor of five. Kinetic isotope effect and photoluminescence quenching experiments reveal that the scavenging of the photoexcited hole by sulfite ions is the rate-limiting step and leads to the observed response of the TOF to pH through actuation of the microgel. Molecular dynamics simulations quantify a greater local concentration of sulfite hole scavengers for pH < pK a .

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

The changing mouse embryo transcriptome at whole tissue and single-cell resolution

During mammalian embryogenesis, differential gene expression gradually builds the identity and complexity of each tissue and organ system 1 . Here we systematically quantified mouse polyA-RNA from day 10.5 of embryonic development to birth, sampling 17 tissues and organs. The resulting developmental transcriptome is globally structured by dynamic cytodifferentiation, body-axis and cell-proliferation gene sets that were further characterized by the transcription factor motif codes of their promoters. We decomposed the tissue-level transcriptome using single-cell RNA-seq (sequencing of RNA reverse transcribed into cDNA) and found that neurogenesis and haematopoiesis dominate at both the gene and cellular levels, jointly accounting for one-third of differential gene expression and more than 40% of identified cell types. By integrating promoter sequence motifs with companion ENCODE epigenomic profiles, we identified a prominent promoter de-repression mechanism in neuronal expression clusters that was attributable to known and novel repressors. Focusing on the developing limb, single-cell RNA data identified 25 candidate cell types that included progenitor and differentiating states with computationally inferred lineage relationships. We extracted cell-type transcription factor networks and complementary sets of candidate enhancer elements by using single-cell RNA-seq to decompose integrative cis -element (IDEAS) models that were derived from whole-tissue epigenome chromatin data. These ENCODE reference data, computed network components and IDEAS chromatin segmentations are companion resources to the matching epigenomic developmental matrix, and are available for researchers to further mine and integrate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Matrix-free subcell residual distribution for Bernstein finite elements: Monolithic limiting

This paper is focused on the aspects of limiting in residual distribution (RD) schemes for high-order finite element approximations to advection problems. Both continuous and discontinuous Galerkin methods are considered in this work. Discrete maximum principles are enforced using algebraic manipulations of element contributions to the global nonlinear system. The required modifications can be carried out without calculating the element matrices and assembling their global counterparts. The components of element vectors associated with the standard Galerkin discretization are manipulated directly using localized subcell weights to achieve optimal accuracy. Low-order nonlinear RD schemes of this kind were originally developed to calculate local extremum diminishing predictors for flux-corrected transport (FCT) algorithms. In the present paper, we incorporate limiters directly into the residual distribution procedure, which makes it applicable to stationary problems and leads to well-posed nonlinear discrete problems. To circumvent the second-order accuracy barrier, the correction factors of monolithic limiting approaches and FCT schemes are adjusted using smoothness sensors based on second derivatives. Finally, the convergence behavior of presented methods is illustrated by numerical studies for two-dimensional test problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Convergence Analysis of the Alternating Anderson–Picard Method for Nonlinear Fixed-Point Problems

Anderson acceleration (AA) has been widely used to solve nonlinear fixed-point problems due to its rapid convergence. This work focuses on a variant of AA in which multiple Picard iterations are performed between each AA step, referred to as the Alternating Anderson–Picard (AAP) method. Furthermore, despite introducing more “slow” Picard iterations, this method has been shown to be efficient and even more robust in both linear and nonlinear cases. However, there is a lack of theoretical analysis for AAP in the nonlinear case. In this paper, we address this gap by establishing the equivalence between AAP and a multisecant-GMRES method that uses GMRES to solve a multisecant linear system at each iteration. From this perspective, we show that AAP “converges” to the Newton-GMRES method. Specifically, as the residual approaches zero, the multisecant matrix, the approximate Jacobian inverse, the search direction, and the optimization gain of AAP converge to their counterparts in the Newton-GMRES method. These connections provide insights for analyzing the asymptotic convergence properties of AAP. Consequently, we show that AAP is locally 𝑞-linear convergent and provide an upper bound for the convergence factor of AAP. To validate the theoretical results, numerical examples are provided.

Anderson acceleration↗
DOE OSTI2021

Cellular mRNA triggers structural transformation of Ebola virus matrix protein VP40 to its essential regulatory form

The Ebola virus matrix protein VP40 forms distinct structures linked to distinct functions in the virus life cycle. Dimeric VP40 is a structural protein associated with virus assembly, while octameric, ring-shaped VP40 is associated with transcriptional control. In this study, we show that suitable nucleic acid is sufficient to trigger a dynamic transformation of VP40 dimer into the octameric ring. Deep sequencing reveals a binding preference of the VP40 ring for the 3' untranslated region of cellular mRNA and a guanine- and adenine-rich binding motif. Complementary analyses of the nucleic-acid-induced VP40 ring by native mass spectrometry, electron microscopy, and X-ray crystal structures at 1.8 and 1.4 Å resolution reveal the stoichiometry of RNA binding, as well as an interface involving a key guanine nucleotide. The host factor-induced structural transformation of protein structure in response to specific RNA triggers in the Ebola virus life cycle presents unique opportunities for therapeutic inhibition.

3′ UTR↗
DOE OSTI2024

Defect-engineered monolayer MoS 2 with enhanced memristive and synaptic functionality for neuromorphic computing

Two-dimensional transition metal dichalcogenides (TMDs)-based memristors are promising candidates for realizing artificial synapses in next-generation computing. However, practical implementation faces several challenges, such as high non-linearity and asymmetry in synaptic weight updates, limited dynamic range, and cycle-to-cycle variability. Here, utilizing optimal-power argon plasma treatment, we significantly enhance the performance matrix of memristors fabricated from monolayer MoS 2 . Our approach not only improves linearity and symmetry in synaptic weight updates but also increases the number of available synaptic weight updates and enhances Spike-Time Dependent Plasticity. Notably, it broadens the switching ratio by two orders, minimizes cycle-to-cycle variability, reduces non-linear factors, and achieves an energy consumption of ~30 fJ per synaptic event. Implementation of these enhancements is demonstrated through Artificial Neural Network simulations, yielding a learning accuracy of ~97% on the MNIST hand-written digits dataset. Our findings underscore the significance of defect engineering as a powerful tool in advancing the synaptic functionality of memristors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Accelerating $\mathrm{LHC}$ event generation with simplified pilot runs and fast $\mathrm{PDF}$s

Poor computing efficiency of precision event generators for LHC physics has become a bottleneck for Monte-Carlo event simulation campaigns. We provide solutions to this problem by focusing on two major components of general-purpose event generators: The PDF evaluator and the matrix-element generator. For a typical production setup in the ATLAS experiment, we show that the two can consume about 80% of the total runtime. Using NLO simulations of pp→ℓ + ℓ - +jets and pp→$t\bar{t}$+jets as an example, we demonstrate that the computing footprint of LHAPDF and SHERPA can be reduced by factors of order 10, while maintaining the formal accuracy of the event sample. The improved codes are made publicly available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

HNRNPU promotes the progression of hepatocellular carcinoma by enhancing CDK2 transcription

Highlights: • The expression of HNRNPU is upregulated in hepatocellular carcinoma. • HNRNPU knockdown significantly inhibits HCC cell proliferation. • HNRNPU knockdown induces cell cycle arrest. • HNRNPU regulates CDK2 transcription by directly binding to the CDK2 gene locus. The nuclear matrix-associated protein Heterogeneous Nuclear Ribonucleoprotein U (HNRNPU), also known as SAF-A, is known to maintain active chromatin structure in mouse hepatocytes. However, the functional roles and molecular mechanisms of HNRNPU in the development of hepatocellular carcinoma (HCC) remain largely unknown. Herein, we found that HNRNPU was upregulated in HCC, and the proliferation of HCC cells was inhibited in vitro and in vivo upon HNRNPU knockdown. Moreover, the upregulation of HNRNPU was correlated with poor prognosis in HCC. Mechanistically, HNRNPU bound to the CDK2 gene locus, a key factor in cell cycle regulation, where it was enriched with H3K27 acetylation (H3K27ac), H3K9 acetylation (H3K9ac), and H3K4 mono-methylation (H3K4me1). Furthermore, HNRNPU knockdown reduced the levels of H3K27ac and H3K9ac at the binding site, where the levels of H3K27 tri-methylation (H3K27me3) were increased, eventually leading to the downregulation of CDK2. Collectively, our results provide a new mechanism whereby HNRNPU promotes HCC development by enhancing the transcription of CDK2.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Experimental and numerical investigation of proppant embedment and conductivity reduction within a fracture in the Caney Shale, Southern Oklahoma, USA

The current worldwide energy supply is insufficient to meet the rising demand. As a result, the energy prices are expected to keep soaring despite the recent increases in a variety of renewable energy resources. Although not renewable, shale oil and gas — “unconventional” hydrocarbon resources are relatively clean forms of energy resources, which still hold a vast share of the energy market. For many oil and gas companies, meeting profitable production goals from shale reservoirs is sometimes challenging, due to the loss of fracture conductivity and premature declines in the production. In this paper we investigate the stress-dependent changes in the hydraulic conductivity of proppant-filled fractures and mechanical fracture–proppant interactions in Caney Shale, a calcareous, organic-rich mudrock, through laboratory experiments and numerical modeling. American Petroleum Institute (API) fracture conductivity tests were conducted using 2% KCl on five locations within the Caney Shale that consisted of selecting three brittle (reservoir) zones and two ductile zones. Confining pressures ranged from 1,000 psi to 12,000 psi at 210 °F. Conductivity, permeability as well as embedment were measured during the test. Also, an additional, laboratory in-situ visualization test was conducted to examine the detailed proppant-shale matrix interaction under elevated stress (3,920 psi effective stress) and temperature (252 °F), with a synthetic reservoir fluid. Our experimental results have confirmed that improved fracture conductivity is attributed to proppant size, and that the increase in porosity of the proppant pack, closure pressure changes and the reduction in fracture conductivity are a function of many factors such as fracture closure stress.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Structural origin of recovered ferroelectricity in BaTi O 3 nanoparticles

Nanoscale BaTiO 3 particles (≈10 nm) prepared by ball-milling a mixture of oleic acid and heptane have been reported to have an electric polarization several times larger than that for bulk BaTiO 3 . In this work, detailed local, intermediate, and long-range structural studies are combined with spectroscopic measurements to develop a model structure of these materials. The X-ray spectroscopic measurements reveal large Ti off-centering as the key factor producing the large spontaneous polarization in the nanoparticles. Temperature-dependent lattice parameter changes reveal the sharpening of the structural phase transitions in these BaTiO 3 nanoparticles compared to the pure nanoparticle systems. In conclusion, sharp crystalline-type peaks in the barium oleate Raman spectra suggest that this component in the composite core-shell matrix, a product of mechanochemical synthesis, stabilizes an enhanced polar structural phase of the BaTiO 3 core nanoparticles.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Subsurface mechanical damage of fused silica glass during grinding by various sub-aperture tools with and without ultrasonics

The subsurface mechanical damage (SSD) depth after grinding fused silica glass with a comprehensive set of sub-aperture fixed abrasive grinding tools [cup, wheel, belt, pad, and rotary face mill (with and without ultrasonics)] and process parameters has been statistically measured using the taper wedge technique and evaluated. Consistent with a previously reported grinding model [J. Non-Cryst. Solids 352, 5601–5617 (2006) Crossref , Materials Science and Technology of Optical Fabrication (Wiley & Sons, 2018)], based on the sliding indentation fracture by sliding particles or asperities where the normal load per particle determines the depth of the fracture (and ultimately the overall SSD depth), the dominant factor controlling SSD depth was found to be the abrasive size regardless of the tool type and process conditions. Compared to full aperture grinding methods, the overall SSD depth was higher using the sub-aperture tools, likely due to the higher effective pressure and higher load per particle distribution. Here, in addition to abrasive size, a significant reduction in SSD depth was achieved by: (1) reducing the load distribution on the abrasive particles via increase in contact area and/or decrease in mechanical loading; (2) using a more compliant host tool medium; and (3) in what we believe is a more novel way, using ultrasonics. Combining low abrasive size, larger contact area, and a compliant host, the 6 µm diamond in a resin matrix (Trizact) on a foam pad led to very low SSD depth (~ 4.6 µm), relatively fast grinding rate (186 mm 3 /h), and little or no figure degradation. This grinding tool/process is an attractive choice for final grind, resulting in significantly reduced polish out (i.e., “grey out”) times. With the rotary face mill tool, the use of ultrasonics consistently led to a SSD depth reduction (ranging from 17%–34%). A new fracture mechanics-based model, to the best of our knowledge, where the relevant normal load is parallel to the feed direction, has been developed to explain how ultrasonics leads to lower SSD depth. The key factors, supported by finite element stress analysis and load measurements, are (1) the initiation of fractures at higher z heights during the tool’s ultrasonic vertical oscillations, thus propagating less deep into workpiece; (2) reduction in load (and therefore reduction in fracture propagation distance) due to smaller tool-workpiece feed direction contact area (again caused by higher heights relative to depth of cut); (3) upward movement of the tool during oscillation leads to fracturing toward the surface instead into the depth; and finally (4) at tool’s lowest point of oscillation cycle, there may not be enough time for the fracture to propagate to its full length.

Optics and optical instruments↗
DOE OSTI2024

Bragg-Primakoff axion photoconversion in crystal detectors

Abstract Axions and axion-like pseudoscalar particles with dimension-5 couplings to photons exhibit coherent Primakoff scattering with ordered crystals at keV energy scales, making for a natural detection technique in searches for solar axions. We find that there are large suppressive corrections, potentially greater than a factor of$$ \mathcal{O} $$ O (10 3 ), to the coherent enhancement when taking into account absorption of the final state photon. This effect has already been accounted for in light-shining-through-wall experiments through the language of Darwin classical diffraction, but is missing from the literature in the context of solar axion searches that use a matrix element approach. We extend the treatment of the event rate with a heuristic description of absorption effects to bridge the gap between these two languages. Furthermore, we explore the Borrmann effect of anomalous absorption in lifting some of the event rate suppression by increasing the coherence length of the conversion. We study this phenomenon in Ge, NaI, and CsI crystal experiments and its impact on the projected sensitivities of SuperCDMS, LEGEND, and SABRE to the solar axion parameter space. Lastly, we comment on the reach of multi-tonne scale crystal detectors and strategies to maximize the discovery potential of experimental efforts in this vein.

Physics↗
DOE OSTI2021

Enhanced relaxed physical factorization preconditioner for coupled poromechanics

The relaxed physical factorization (RPF) preconditioner is a recent algorithm allowing for the efficient and robust solution to the block linear systems arising from the three-field displacement-velocity-pressure formulation of coupled poromechanics. For its application, however, it is necessary to invert blocks with the algebraic form C^ = (C + βFF T ), where C is a symmetric positive definite matrix, FF T a rank-deficient term, and β a real non-negative coefficient. The inversion of C^, performed in an inexact way, can become unstable for large values of β, as it usually occurs at some stages of a full poromechanical simulation. In this work, we propose a family of algebraic techniques to stabilize the inexact solve with C^. This strategy can prove useful in other problems as well where such an issue might arise, such as augmented Lagrangian preconditioning techniques for Navier-Stokes or incompressible elasticity. First, we introduce an iterative scheme obtained by a natural splitting of matrix C^. Second, we develop a technique based on the use of a proper projection operator annihilating the near-kernel modes of C^. Both approaches give rise to a novel class of preconditioners denoted as Enhanced RPF (ERPF). Furthermore, effectiveness and robustness of the proposed algorithms are demonstrated in both theoretical benchmarks and real-world large-size applications, outperforming the native RPF preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Represent precipitation-induced geological hazards in Earth system models using artificial intelligence

Precipitation-induced geological hazards, such as debris flow, landslides, mudflow and rockfalls (hereafter referred to as landslides), pose serious threats to public safety in many areas through the world. As residential properties and infrastructure in the US have increasingly expanded into landslide-prone areas and the drivers of landslides (e.g., wildfires and hurricanes) are predicted to intensify under climate warming, losses and fatalities from landslides are likely to increase in the future. Although our understanding of geoenvironmental factors and mechanisms contributing to landslides has greatly improved, only moderate progress has been made in predicting landslides out of well-studied watersheds. Furthermore, explicitly representing various landslide-related processes in Earth system models (ESMs), from the buckling of local bearing elements in granular materials, to frictional sliding between grains, formation of microcracks in the soil matrix, rupture of capillary bridges, or breakage of plant roots3 is still very unlikely within the next decade, even with the help of exascale computers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Feedstock Variability: Differences Among Biomass Anatomical Fractions and Tissues

Corn stover is a versatile feedstock for producing biofuels and other sustainable products. However, its diverse chemical and physical attributes can hinder the processes that convert it into these compounds. The composition of corn stover varies according to factors including corn plant variety, cultivation conditions and maturity at harvest. In essence, corn stover consists of the polymers cellulose, hemicellulose and lignin. The structural core of plant cell walls consists of elongated chains of glucose molecules arranged in bundles of cellulose. Additionally, lower concentrations of hemicellulose, a complex carbohydrate composed of various sugars, form part of the matrix surrounding cellulose bundles. Lignin is a complex aromatic polymer that reinforces plant cell walls. Corn stover generally comprises 40%-50% cellulose, 20%-30% hemicellulose, and 20%-30% lignin. The ratios of these components vary based on the specific type of corn plant and the growth conditions. For example, corn plants cultivated in drought conditions exhibit higher lignin levels than those grown under optimal conditions. The compositional diversity resulting from growth conditions is most pronounced among the different anatomical fractions of the corn plant. Understanding and accommodating the variation among the anatomical fractions of corn stover can enhance the effective utilization of this resource.

09 BIOMASS FUELS↗
DOE OSTI2021

Island for gravitationally prepared state and pseudo entanglement wedge

We consider spacetime initiated by a finite-sized initial boundary as a generalization of the Hartle-Hawking no-boundary state. We study entanglement entropy of matter state prepared by such spacetime. We find that the entanglement entropy for large subregion is given either by the initial state entanglement or the entanglement island, preventing the entropy to grow arbitrarily large. Consequently, the entanglement entropy is always bounded from above by the boundary area of the island, leading to an entropy bound in terms of the island. The island I is located in the analytically continued spacetime, either at the bra or the ket part of the spacetime in Schwinger-Keldysh formalism. The entanglement entropy is given by an average of complex pseudo generalized entropy for each entanglement island. We find a necessary condition of the initial state to be consistent with the strong sub-additivity, which requires that any probe degrees of freedom are thermally entangled with the rest of the system. We then find a large parameter region where the spacetime with finite-sized initial boundary, which does not have the factorization puzzle at leading order, dominates over the Hartle-Hawking no-boundary state or the bra ket wormhole. Due to the absence of a moment of time reflection symmetry, the island in our setup is a generalization of the entanglement wedge, called pseudo entanglement wedge. In pseudo entanglement wedge reconstruction, we consider reconstructing the bulk matter transition matrix on A ∪ I, from a fine-grained state on A. The bulk transition matrix is given by a thermofield double state with a projection by the initial state. We also provide an AdS/BCFT model by considering EOW branes with corners. We also find the exponential hardness of such reconstruction task using a generalization of Python’s lunch conjecture to pseudo generalized entropy.

2D gravity↗
DOE OSTI2025

Quantum chaos, integrability, and late times in the Krylov basis

Quantum chaotic systems are conjectured to display a spectrum whose fine-grained features (gaps and correlations) are well described by random matrix theory (RMT). We propose and develop a complementary version of this conjecture: quantum chaotic systems display a Lanczos spectrum whose local means and covariances are well described by RMT. To support this proposal, we first demonstrate its validity in examples of chaotic and integrable systems. We then show that for Haar-random initial states in RMTs the mean and covariance of the Lanczos spectrum suffice to produce the full long-time behavior of general survival probabilities including the spectral form factor, as well as the spread complexity. In addition, for initial states with continuous overlap with energy eigenstates, we analytically find the long-time averages of the probabilities of Krylov basis elements in terms of the mean Lanczos spectrum. This analysis suggests a notion of eigenstate complexity, the statistics of which differentiate integrable systems and classes of quantum chaos. Lastly, we clarify the relation between spread complexity and the universality classes of RMT by exploring various values of the Dyson index and Poisson distributed spectra.

Combinatorics↗