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 OSTI2020

Fracture resistance of hierarchical Cu–Mo nanocomposite thin films

We report the fracture behavior of high-temperature co-sputtered Cu/Mo nanocomposite thin films was investigated through in situ 3-point bend testing of pre-notched microbeams in SEM. The as-synthesized nanocomposites present hierarchical nano/microstructures composed of a matrix of phase-separated Cu–Mo with nanoscale ligaments dispersed with sub-micron scale Cu-rich islands containing Mo nano-precipitates. Results show a significant crack growth resistance in the hierarchical nanocomposite, several times higher than that measured in Cu/Mo nanoscale multilayers. Based on electron microscopy characterization, three mechanisms of crack growth resistance in the hierarchical structures are proposed: crack bridging by the Cu-rich layer, crack deflection via shear along the Cu/Mo interface, and multiple cracking. This work demonstrates an approach to increase toughness in high strength nanocomposites through interface micro-structure design.

36 MATERIALS SCIENCE↗
NASA NTRS1984

On the development of hierarchical solution strategies for nonlinear finite element formulations

This paper develops a hierarchical type solution scheme which can handle the field equations associated with nonlinear finite element simulations. The overall procedure possesses various levels of application namely degree of freedom, nodal, elemental, substructural as well as global. In particular iteration, updating, assembly and solution control occurs at the various hierarchical levels. Due to the manner of formulation, the degree of matrix inversion depends on the size of the various hierarchical partitioned groups. In this context, degree of freedom partitioning requires no inversion. To benchmark the overall scheme, the results of several numerical examples are presented.

Padovan, J.↗
DOE OSTI2025

Hierarchically Structured Vitrimer Biocomposites for Sustainable Manufacturing

Polymers containing dynamic covalent bonds (DCBs) exhibit thermoplastic-like flow above their topology freezing temperature (T v ) while maintaining thermoset-like properties below it, making them promising for sustainable manufacturing. However, their large-scale adoption remains limited due to challenges in accurately determining T v and achieving efficient fiber-matrix bonding in composite applications. Here, hierarchically structured epoxy-anhydride-based polyester vitrimer composites reinforced with cellulosic filaments is demonstrated, where hydroxyl groups on fiber surfaces participate directly in transesterification with the matrix. Further, this dynamic interfacial bonding delivers exceptional mechanical properties, including ≈70 MPa shear strength and >10% strain-to-failure, while enabling thermal malleability. Using a combination of nuclear magnetic resonance and nano-infrared spectroscopies, direct evidence is provided that chemical bond exchange begins well below the conventionally measured T v , supporting the hypothesis that rheologically determined T v reflects a combination of chemical exchange and frictional dynamics rather than a discrete transition. The composites demonstrate excellent processability through vacuum-assisted resin transfer molding and maintain >90% of their mechanical properties after multiple thermal reforming cycles. These findings advance both the fundamental understanding of vitrimeric transitions and the practical development of sustainable, high-performance composite materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
NASA NTRS2004

Fossil Signatures Using Elemental Abundance Distributions and Bayesian Probabilistic Classification

Elemental abundances (C6, N7, O8, Na11, Mg12, Al3, P15, S16, Cl17, K19, Ca20, Ti22, Mn25, Fe26, and Ni28) were obtained for a set of terrestrial fossils and the rock matrix surrounding them. Principal Component Analysis extracted five factors accounting for the 92.5% of the data variance, i.e. information content, of the elemental abundance data. Hierarchical Cluster Analysis provided unsupervised sample classification distinguishing fossil from matrix samples on the basis of either raw abundances or PCA input that agreed strongly with visual classification. A stochastic, non-linear Artificial Neural Network produced a Bayesian probability of correct sample classification. The results provide a quantitative probabilistic methodology for discriminating terrestrial fossils from the surrounding rock matrix using chemical information. To demonstrate the applicability of these techniques to the assessment of meteoritic samples or in situ extraterrestrial exploration, we present preliminary data on samples of the Orgueil meteorite. In both systems an elemental signature produces target classification decisions remarkably consistent with morphological classification by a human expert using only structural (visual) information. We discuss the possibility of implementing a complexity analysis metric capable of automating certain image analysis and pattern recognition abilities of the human eye using low magnification optical microscopy images and discuss the extension of this technique across multiple scales.

Hoover, Richard B.↗
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↗
NASA NTRS1999

Development of GENOA Progressive Failure Parallel Processing Software Systems

A capability consisting of software development and experimental techniques has been developed and is described. The capability is integrated into GENOA-PFA to model polymer matrix composite (PMC) structures. The capability considers the physics and mechanics of composite materials and structure by integration of a hierarchical multilevel macro-scale (lamina, laminate, and structure) and micro scale (fiber, matrix, and interface) simulation analyses. The modeling involves (1) ply layering methodology utilizing FEM elements with through-the-thickness representation, (2) simulation of effects of material defects and conditions (e.g., voids, fiber waviness, and residual stress) on global static and cyclic fatigue strengths, (3) including material nonlinearities (by updating properties periodically) and geometrical nonlinearities (by Lagrangian updating), (4) simulating crack initiation. and growth to failure under static, cyclic, creep, and impact loads. (5) progressive fracture analysis to determine durability and damage tolerance. (6) identifying the percent contribution of various possible composite failure modes involved in critical damage events. and (7) determining sensitivities of failure modes to design parameters (e.g., fiber volume fraction, ply thickness, fiber orientation. and adhesive-bond thickness). GENOA-PFA progressive failure analysis is now ready for use to investigate the effects on structural responses to PMC material degradation from damage induced by static, cyclic (fatigue). creep, and impact loading in 2D/3D PMC structures subjected to hygrothermal environments. Its use will significantly facilitate targeting design parameter changes that will be most effective in reducing the probability of a given failure mode occurring.

Abdi, Frank↗
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↗
NASA NTRS1997

Computational Infrastructure for Engine Structural Performance Simulation

Select computer codes developed over the years to simulate specific aspects of engine structures are described. These codes include blade impact integrated multidisciplinary analysis and optimization, progressive structural fracture, quantification of uncertainties for structural reliability and risk, benefits estimation of new technology insertion and hierarchical simulation of engine structures made from metal matrix and ceramic matrix composites. Collectively these codes constitute a unique infrastructure readiness to credibly evaluate new and future engine structural concepts throughout the development cycle from initial concept, to design and fabrication, to service performance and maintenance and repairs, and to retirement for cause and even to possible recycling. Stated differently, they provide 'virtual' concurrent engineering for engine structures total-life-cycle-cost.

Chamis, Christos C.↗
NASA NTRS1995

Sensitivity analysis for large-deflection and postbuckling responses on distributed-memory computers

A computational strategy is presented for calculating sensitivity coefficients for the nonlinear large-deflection and postbuckling responses of laminated composite structures on distributed-memory parallel computers. The strategy is applicable to any message-passing distributed computational environment. The key elements of the proposed strategy are: (1) a multiple-parameter reduced basis technique; (2) a parallel sparse equation solver based on a nested dissection (or multilevel substructuring) node ordering scheme; and (3) a multilevel parallel procedure for evaluating hierarchical sensitivity coefficients. The hierarchical sensitivity coefficients measure the sensitivity of the composite structure response to variations in three sets of interrelated parameters; namely, laminate, layer and micromechanical (fiber, matrix, and interface/interphase) parameters. The effectiveness of the strategy is assessed by performing hierarchical sensitivity analysis for the large-deflection and postbuckling responses of stiffened composite panels with cutouts on three distributed-memory computers. The panels are subjected to combined mechanical and thermal loads. The numerical studies presented demonstrate the advantages of the reduced basis technique for hierarchical sensitivity analysis on distributed-memory machines.

Watson, Brian C.↗
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↗