Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hierarchical matrix”

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

In Situ Growth of Crystalline and Polymer-Incorporated Amorphous ZIFs in Polybenzimidazole Achieving Hierarchical Nanostructures for Carbon Capture

We report mixed matrix materials (MMMs) hold great potential for membrane gas separations by merging nanofillers with unique nanostructures and polymers with excellent processability. In situ growth of the nanofillers is adapted to mitigate interfacial incompatibility to avoid the selectivity loss. Surprisingly, functional polymers have not been exploited to co-grow the nanofillers for membrane applications. Herein, in situ synergistic growth of crystalline zeolite imidazole framework-8 (ZIF-8) in polybenzimidazole (PBI), creating highly porous structures with high gas permeability, is demonstrated. More importantly, PBI contains benzimidazole groups (similar to the precursor for ZIF-8, i.e., 2-methylimidazole) and induces the formation of amorphous ZIFs, enhancing interfacial compatibility and creating highly size-discriminating bottlenecks. For instance, the formation of 15 mass% ZIF-8 in PBI improves H- 2 permeability and H- 2 /CO 2 selectivity by approximate to 100% at 35 degrees C, breaking the permeability/selectivity tradeoff. This work unveils a new platform of MMMs comprising functional polymer-incorporated amorphous ZIFs with hierarchical nanostructures for various applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Accelerating structural dynamics simulations with localised phenomena through matrix compression and projection‐based model order reduction

In this work, a novel approach is introduced for accelerating the solution of structural dynamics problems in the presence of localised phenomena, such as cracks. For this category of problems, conventional projection-based Model Order Reduction (MOR) methods are either limited with respect to the range of system configurations that can be represented or require frequent solutions of the Full Order Model (FOM) to update the low-dimensional spaces, in which solutions are represented. In the proposed approach, low-dimensional spaces, constructed for the healthy structure, are enriched with appropriately selected columns of the flexibility matrix of the system. It can be shown that these spaces contain the solution to the original problem for the static case, while their dimension is much smaller. In order to allow their online construction for arbitrary localised features, the full flexibility matrix of the system should be available. To this end, a hierarchical representation is used for the matrices involved, allowing to compute the flexibility matrix efficiently and with reduced memory requirements. The resulting method offers significant speedups, without sacrificing the flexibility and accuracy of the full order model. The performance and limitations of the approach are studied through a series of examples in structural dynamics.

fracture mechanics↗
DOE OSTI2023

Structure-informed clustering for population stratification in association studies

Background: Identifying variants associated with complex traits is a challenging task in genetic association studies due to linkage disequilibrium (LD) between genetic variants and population stratification, unrelated to the disease risk. Existing methods of population structure correction use principal component analysis or linear mixed models with a random effect when modeling associations between a trait of interest and genetic markers. However, due to stringent significance thresholds and latent interactions between the markers, these methods often fail to detect genuinely associated variants. Results: To overcome this, we propose CluStrat, which corrects for complex arbitrarily structured populations while leveraging the linkage disequilibrium induced distances between genetic markers. It performs an agglomerative hierarchical clustering using the Mahalanobis distance covariance matrix of the markers. In simulation studies, we show that our method outperforms existing methods in detecting true causal variants. Applying CluStrat on WTCCC2 and UK Biobank cohorts, we found biologically relevant associations in Schizophrenia and Myocardial Infarction. CluStrat was also able to correct for population structure in polygenic adaptation of height in Europeans. Conclusions: CluStrat highlights the advantages of biologically relevant distance metrics, such as the Mahalanobis distance, which captures the cryptic interactions within populations in the presence of LD better than the Euclidean distance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Hierarchical Bayesian Thermonuclear Rate for the 7 Be(n,p) 7 Li Big Bang Nucleosynthesis Reaction

Big Bang nucleosynthesis provides the earliest probe of standard model physics, at a time when the universe was less than 1000 seconds old. It determines the abundances of the lightest nuclides, which give rise to the subsequent history of the visible matter in the universe. This work derives new 7 Be(n,p) 7 Li thermonuclear reaction rates based on all available experimental information. This reaction sensitively impacts the primordial abundances of 7 Be and 7 Li during big bang nucleosynthesis. We critically evaluate all available data and disregard experimental results that are questionable. For the nuclear model, we adopt an incoherent sum of single-level, two-channel, R -matrix approximation expressions, which are implemented into a hierarchical Bayesian model, to analyze the remaining six data sets we deem most reliable. In the fitting of the data, we consistently model all known sources of uncertainty, including discrepant absolute normalizations of different data sets, and also take the variation of the neutron and proton channel radii into account, hence providing less biased estimates of the 7 Be(n,p) 7 Li thermonuclear rates. From the resulting posteriors, we extract R -matrix parameters ($E_r, γ^2_n, γ^2_p$) and derive excitation energies and partial and total widths. Our fit is sensitive to the contributions of the first three levels above the neutron threshold. Reaction rates were computed by integrating 10,000 samples of the reduced cross section. Our 7 Be(n,p) 7 Li thermonuclear rates have uncertainties between 1.5% and 2.0% at temperatures of ≤1 GK. Finally, we compare our rates to previous results and find that the 7 Be(n,p) 7 Li rates most commonly used in big bang simulations have uncertainties that are too optimistic.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initialization

Obtaining lightweight and accurate approximations of discretized objective functional Hessians in inverse problems governed by partial differential equations (PDEs) is essential to make both deterministic and Bayesian statistical large-scale inverse problems computationally tractable. The cubic computational complexity of dense linear algebraic tasks, such as Cholesky factorization, that provide a means to sample Gaussian distributions and determine solutions of Newton linear systems is a computational bottleneck at large-scale. These tasks can be reduced to log-linear complexity by utilizing hierarchical off-diagonal low-rank (HODLR) matrix approximations. In this work, we show that a class of Hessians that arise from inverse problems governed by PDEs are well approximated by the HODLR matrix format. In particular, we study inverse problems governed by PDEs that model the instantaneous viscous flow of ice sheets. In these problems, we seek a spatially distributed basal sliding parameter field such that the flow predicted by the ice sheet model is consistent with ice sheet surface velocity observations. Here, we demonstrate the use of HODLR Hessian approximation to efficiently sample the Laplace approximation of the posterior distribution with covariance further approximated by HODLR matrix compression. Computational studies are performed which illustrate ice sheet problem regimes for which the Gauss–Newton data-misfit Hessian is more efficiently approximated by the HODLR matrix format than the low-rank (LR) format. We then demonstrate that HODLR approximations can be favorable, when compared to global LR approximations, for large-scale problems by studying the data-misfit Hessian associated with inverse problems governed by the first-order Stokes flow model on the Humboldt glacier and Greenland ice sheet.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Additive manufacturing enabled synergetic strengthening of bimodal reinforcing particles for aluminum matrix composites

An additive manufactured TiB 2 /Al-Cu-Mg-Ni composite with a minor amount of Sc was fabricated by laser powder bed fusion (LPBF). The composite shows a yield strength of ~370 MPa, almost doubling the strength of its wrought matrix counterpart, and an elongation of ~7 %. The superior mechanical properties are attributed to a unique micro-nano hierarchical microstructure, consisting of nanoscale and microscale TiB 2 particles dispersed in a matrix of fine aluminum grains (3 µm) together with intragranular intermetallic nanoparticles and nano-cellular networks (cell size 30 nm). The formation of the TiB 2 nanoparticles is the consequence of partial dissolution of the TiB 2 particulates and the enrichment of Sc in the newly formed nanoparticles. Remarkable strengthening effects are achieved by the bimodal TiB 2 particles, intermetallic nanoparticles and intragranular nano-cellular networks. Finally, this study provides new insights into the role of additive manufacturing in tailoring the microstructure of particulate reinforced metal matrix composites (MMCs) with advanced properties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Evaluating the efficiency of touch-spun scaffolds in producing dense cell cultures for tissue engineering applications

Developing efficient scaffolds for long-term cell cultivation remains a challenge in tissue engineering. Biomimetic approaches aim to create a three-dimensional (3D) extracellular matrix (ECM)-like fiber network with a tunable hierarchical structure to promote sufficient cell attachment, differentiation, and overall viability. Among the fiber fabrication methods documented in the literature, mechanical fiber drawing techniques, such as touch-spinning, have garnered significant research interest. This is due to the simplicity of the equipment, the ability to control fiber diameter and interfiber spacing at the nanoscale without the need for external fields, and the absence of specific requirements for material dielectric properties. Despite the advantages of mechanically drawn scaffolds in biomedical research, the methodologies for cell culturing and analysis for these materials have not been adequately addressed. In this study, we assess the potential of touch-spun scaffolds in promoting NIH/3T3-GFP fibroblast cell growth for tissue engineering applications. Polycaprolactone/polyethylene oxide (PCL/PEO)-based fiber arrays with a controlled interfiber spacing of 91.9 ± 25.0 μm (N = 50) were fabricated using a modified touch-spinning apparatus and then assembled into 2D and 3D scaffolds through additive manufacturing technology. A comparative cell analysis conducted for single- and multi-layered structures showed that the 3D touch-spun scaffolds support healthy growth of up to 6.5 million fibroblast cells within 21 days and offer enhanced cell viability compared to conventional 2D fiber scaffolds, as confirmed by the Presto Blue assay. Furthermore, the metabolic activity of fibroblasts on 3D scaffolds assessed by the MTT test is approximately four times higher than that of the positive control, making the 3D touch-spun materials ideal for long-term cell culture applications.

Peranidze, Kristina [University of Georgia, Athens↗
DOE OSTI2022

Particle physics and cosmology of the string derived no-scale flipped SU(5)

In a recent paper, we identified a cosmological sector of a flipped SU(5) model derived in the free fermionic formulation of the heterotic superstring, containing the inflaton and the goldstino superfields with a superpotential leading to Starobinsky type inflation, while SU(5) × U(1) is still unbroken. Here, we study the properties and phenomenology of the vacuum after the end of inflation, where the gauge group is broken to the Standard Model. We identify a set of vacuum expectation values, triggered by the breaking of an anomalous U(1) A gauge symmetry at roughly an order of magnitude below the string scale, that solve the F and D-flatness supersymmetric conditions up to 6th order in the superpotential which is explicitly computed, leading to a successful particle phenomenology. In particular, all extra colour triplets become superheavy guaranteeing observable proton stability, while the Higgs doublet mass matrix has a massless pair eigenstate with realistic hierarchical Yukawa couplings to quarks and leptons. The supersymmetry breaking scale is constrained to be high, consistent with the non observation of supersymmetry signals at the LHC.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A parallel hierarchical blocked adaptive cross approximation algorithm

This article presents a low-rank decomposition algorithm based on subsampling of matrix entries. The proposed algorithm first computes rank-revealing decompositions of submatrices with a blocked adaptive cross approximation (BACA) algorithm, and then applies a hierarchical merge operation via truncated singular value decompositions (H-BACA). The proposed algorithm significantly improves the convergence of the baseline ACA algorithm and achieves reduced computational complexity compared to the traditional decompositions such as rank-revealing QR. Numerical results demonstrate the efficiency, accuracy, and parallel scalability of the proposed algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

PETSc/TAO Users Manual (Rev. 3.19)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than packages such as MATLAB or a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or lldb, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

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

Using Crystallization to Control Filler Dispersion and Vice Versa in Polymer Nanocomposites

Our overall goal was to exploit our ability to control the spatial dispersion of spherical nanoparticles in a semicrystalline polymer matrix to tune the properties of the resulting nanocomposite. In particular, our two tasks were: Hierarchically organize spherical nanoparticles into unique morphologies by varying the ratio of the polymer crystallization rate to the filler diffusion rate during matrix solidification. Organize nanoparticles in the matrix melt into superstructrues as a means to tailor polymer crystallization rates, and potentially the crystal structure and/or morphology of the semicrystalline matrix.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Harnessing peptide–cellulose interactions to tailor the performance of self-assembled, injectable hydrogels

Taking inspiration from natural systems, such as spider silk and mollusk nacre, that employ hierarchical assembly to attain robust material performance, we leveraged matrix–filler interactions within reinforced polymer–peptide hybrids to create self-assembled hydrogels with enhanced properties. Specifically, cellulose nanocrystals (CNCs) were incorporated into peptide–polyurea (PPU) hybrid matrices to tailor key hydrogel features through matrix–filler interactions. Herein, we examined the impact of peptide repeat length and CNC loading on hydrogelation, morphology, mechanics, and thermal behavior of PPU/CNC composite hydrogels. The addition of CNCs into PPU hydrogels resulted in increased gel stiffness; however, the extent of reinforcement of the nanocomposite gels upon nanofiller inclusion also was driven by PPU architecture. Temperature-promoted stiffening transitions observed in nanocomposite PPU hydrogels were dictated by peptide segment length. Analysis of the peptide secondary structure confirmed shifts in the conformation of peptidic domains (α-helices or β-sheets) upon CNC loading. Finally, PPU/CNC hydrogels were probed for their injectability characteristics, demonstrating that nanofiller–matrix interactions were shown to aid rapid network reformation (∼10 s) upon cessation of high shear forces. Overall, this research showcases the potential of modulating matrix–filler interactions within PPU/CNC hydrogels through strategic system design, enabling the tuning of functional hydrogel characteristics for diverse applications.

42 ENGINEERING↗
DOE OSTI2020

PETSc Users Manual (Rev. 3.13)

This manual describes the use of PETSc for the numerical solution of partial differential equations and related problems on high-performance computers. The Portable, Extensible Toolkit for Scientific Computation (PETSc) is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all message-passing communication. PETSc includes an expanding suite of parallel linear solvers, nonlinear solvers, and time integrators that may be used in application codes written in Fortran, C, C++, and Python. PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or dbx, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than \rolling them" yourself.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Pyrolysis in Molten Salts Converts Plastics into a Mesoporous Electrocatalyst with a High Density of Atomic Fe–N–C Active Site

Upcycling plastic waste into value-added products is a promising strategy for both economic and environmental sustainability. However, achieving control over structure–property correlations during the upcycling process remains challenging. Here, in this study, we investigate the modulation of active site formation during pyrolytic conversion of plastics to carbon by composition-controlled molten salts. Using melamine formaldehyde (MF) foams as a precursor, we demonstrate that their pyrolysis in eutectic chloride molten salts with Fe 2+ , Na + , and K + ions directs the carbonization process toward the formation of mesoporous functional carbon enriched with atomically dispersed Fe–N–C sites, which are well-known for their catalytic activity in oxygen reduction reaction (ORR). Compared to pyrolysis in the FeCl 2 -only salt, the incorporation of alkali metal ions (Na + , K + ) in the eutectic mixture facilitates nitrogen retention in the carbon matrix and modulates iron speciation. This synergistic environment promotes a higher density and more uniform dispersion of Fe–N–C moieties within the carbon matrix. Consequently, the resulting catalyst exhibits a high surface area, hierarchical porosity, and improved electrochemical properties. When tested as an ORR catalyst, the Fe–N–C-enriched catalyst delivers a high half-wave potential (E 1/2 = 0.875 V vs RHE) with nearly exclusive 4e – pathway. These properties are comparable to those of commercial Pt/C catalysts, along with improved stability.

Fe−N−C active sites↗
DOE OSTI2025

ThinCurr: An open-source 3D thin-wall eddy current modeling code for the analysis of large-scale systems of conducting structures

In this paper we present a new thin-wall eddy current modeling code, ThinCurr, for studying inductively-coupled currents in 3D conducting structures -- with primary application focused on the interaction between currents flowing in coils, plasma, and conducting structures of magnetically-confined plasma devices. The code utilizes a boundary finite element method on an unstructured, triangular grid to accurately capture device structures. The new code, part of the broader Open FUSION Toolkit, is open-source and designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. Scalability to large models is enabled through use of hierarchical off-diagonal low-rank compression of the inductance matrix, which is otherwise dense. Ease of handling large models of complicated geometry is further supported by automatic determination of supplemental elements through a greedy homology approach. Here, a detailed description of the numerical methods of the code and verification of the implementation of those methods using cross-code comparisons against the VALEN code and Ansys commercial analysis software is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Dynamics of photosynthetic light harvesting systems interacting with N-photon Fock states

Here, we develop a method to simulate the excitonic dynamics of realistic photosynthetic light harvesting systems, including non-Markovian coupling to phonon degrees of freedom, under excitation by N-photon Fock state pulses. This method combines the input–output and the hierarchical equations of motion formalisms into a double hierarchy of density matrix equations. We show analytically that under weak field excitation relevant to natural photosynthesis conditions, an N-photon Fock state input and a corresponding coherent state input give rise to equal density matrices in the excited manifold. However, an N-photon Fock state input induces no off-diagonal coherence between the ground and excited subspaces, in contrast with the coherences created by a coherent state input. We derive expressions for the probability to absorb a single Fock state photon with or without the influence of phonons. For short pulses (or, equivalently, wide bandwidth pulses), we show that the absorption probability has a universal behavior that depends only upon a system-dependent effective energy spread parameter Δ and an exciton–light coupling constant Γ. This holds for a broad range of chromophore systems and for a variety of pulse shapes. We also analyze the absorption probability in the opposite long pulse (narrow bandwidth) regime. We then derive an expression for the long time emission rate in the presence of phonons and use it to study the difference between collective vs independent emission. Finally, we present a numerical simulation for the LHCII monomer (14-mer) system under single photon excitation that illustrates the use of the double hierarchy equations.

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

Evaluation of a pair-based, joint-likelihood association approach for regional infrasound event identification

SUMMARY A Bayesian framework for the association of infrasonic detections is presented and evaluated for analysis at regional propagation scales. A pair-based, joint-likelihood association approach is developed that identifies events by computing the probability that individual detection pairs are attributable to a hypothetical common source and applying hierarchical clustering to identify events from the pair-based analysis. The framework is based on a Bayesian formulation introduced for infrasonic source localization and utilizes the propagation models developed for that application with modifications to improve the numerical efficiency of the analysis. Clustering analysis is completed using hierarchical analysis via weighted linkage for a non-Euclidean distance matrix defined by the negative log-joint-likelihood values. The method is evaluated using regional synthetic data with propagation distances of hundreds of kilometres in order to study the sensitivity of the method to uncertainties and errors in backazimuth and time of arrival. The method is found to be robust and stable for typical uncertainties, able to effectively distinguish noise detections within the data set from those in events, and can be made numerically efficient due to its ease of parallelization.

58 GEOSCIENCES↗