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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2022

A non-hierarchical correlation discrete variable representation

The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. Here, the non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH 3 demonstrate the accuracy of the non-hierarchical CDVR approach.

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

Bio-inspired redox mediated electrolyte for high performance flexible supercapacitor applications over broad temperature domain

For practical supercapacitors, recent advances in electrolytes have diminished long-standing problems such chemical stability, low ionic conductivity and limited electrochemical performance. Herewith, we communicate a facile method to produce anhydrous redox-active Bio-inspired electrolytes for flexible supercapacitors with application over a broad temperature range. The glycerol (Gly) gel electrolyte consists of various fractions of potassium hydroxide (KOH) and redox active unit, i.e., ammonium molybdate, (NH 4 ) 2 MoO 4 (abbreviated as Mo), in the anhydrous state. The doped glycerol matrix acts as a framework that provides a stable gel electrolyte, forming hierarchical pathway for ion diffusion. The developed gel network with Gly-KOH-MoX (X = 3, 5, 10, and 15%, w/w) provides an effective route to carry the ions onto the pores of the electrodes. The supercapacitor with Gly-KOH5-Mo10 gel provides a specific capacitance of 328 F g –1 , and the same device achieved an energy density of 45.6 Wh kg –1 at a power density of 497 W kg –1 . In addition, the supercapacitor maintained extraordinary cyclic robustness and wide temperature tolerance. Thus, the low cost and environmentally benign, cyclic durability of Gly-KOH-Mo makes these gel electrolytes a potential candidate for use in various energy storage devices.

25 ENERGY STORAGE↗
NASA NTRS2016

Improved Convergence and Robustness of USM3D Solutions on Mixed-Element Grids

Several improvements to the mixed-elementUSM3Ddiscretization and defect-correction schemes have been made. A new methodology for nonlinear iterations, called the Hierarchical Adaptive Nonlinear Iteration Method, has been developed and implemented. The Hierarchical Adaptive Nonlinear Iteration Method provides two additional hierarchies around a simple and approximate preconditioner of USM3D. The hierarchies are a matrix-free linear solver for the exact linearization of Reynolds-averaged Navier-Stokes equations and a nonlinear control of the solution update. Two variants of the Hierarchical Adaptive Nonlinear Iteration Method are assessed on four benchmark cases, namely, a zero-pressure-gradient flat plate, a bump-in-channel configuration, the NACA 0012 airfoil, and a NASA Common Research Model configuration. The new methodology provides a convergence acceleration factor of 1.4 to 13 over the preconditioner-alone method representing the baseline solver technology.

Pandya, Mohagna J.↗
NASA NTRS2016

Improved Convergence and Robustness of USM3D Solutions on Mixed-Element Grids

Several improvements to the mixed-element USM3D discretization and defect-correction schemes have been made. A new methodology for nonlinear iterations, called the Hierarchical Adaptive Nonlinear Iteration Method, has been developed and implemented. The Hierarchical Adaptive Nonlinear Iteration Method provides two additional hierarchies around a simple and approximate preconditioner of USM3D. The hierarchies are a matrix-free linear solver for the exact linearization of Reynolds-averaged Navier-Stokes equations and a nonlinear control of the solution update. Two variants of the Hierarchical Adaptive Nonlinear Iteration Method are assessed on four benchmark cases, namely, a zero-pressure-gradient flat plate, a bump-in-channel configuration, the NACA 0012 airfoil, and a NASA Common Research Model configuration. The new methodology provides a convergence acceleration factor of 1.4 to 13 over the preconditioner-alone method representing the baseline solver technology.

Pandya, Mohagna J.↗
NASA NTRSDate not supplied

Hierarchical Multiscale Process Modeling of A Textile Composite Y-Joint for the Aurora D8 Aircraft

Woven polymer matrix composites (PMCs) are widely used in aerospace applications for their favorable specific properties. However, the behavior of the inter- and intra-tow matrix during manufacturing of woven PMCs is difficult to isolate experimentally due to geometric and processing complexities. Classical closed-form homogenization approaches allow for composite thermos-chemo-mechanical properties to be characterized at any fiber volume fraction. In this work, virtual curing is applied to the inter- and intra-tow matrix materials of a plain weave repeating unit cell (RUC). Thermo-chemo-mechanical property evolution within the tows is calculated using the Concentric Cylinder Model (CCM) and the Rule of Mixtures (ROM). Curing is implemented through user-written subroutines within the commercial finite element method software Abaqus. Virtual curing results are compared with those from a linear-elastic thermal cool down simulation in Abaqus. Results indicate that tow property evolution during curing has a non-negligible effect on the final stress state.

curing↗
DOE OSTI2021

Hierarchical quark and lepton masses with sequential U(1) gauge symmetry

Instead of right-handed neutrino singlets, the standard model is extended to include lepton triplets (Σ + , Σ 0 , Σ - ). Each quark and lepton family may now transform under an anomaly-free U(1) X gauge symmetry, known already for many years. A new sequential application is presented, using just the one Higgs doublet of the standard model, together with two U(1) X Higgs singlets. The resulting structure has hierarchical quark and lepton masses, as well as a viable seesaw neutrino mass matrix.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Butterfly-Accelerated Volume Integral Equation Solver for Broad Permittivity and Large-Scale Electromagnetic Analysis

In this work, a butterfly-accelerated volume integral equation (VIE) solver is proposed for fast and accurate electromagnetic (EM) analysis of scattering from heterogeneous objects. The proposed solver leverages the hierarchical off-diagonal butterfly (HOD-BF) scheme to construct the system matrix and obtain its approximate inverse, used as a preconditioner. Complexity analysis and numerical experiments validate the O(N log 2 N) construction cost of the HOD-BF-compressed system matrix and O(N log 1.5 N) inversion cost for the preconditioner, where N is the number of unknowns in the high-frequency EM scattering problem. For many practical scenarios, the proposed VIE solver requires less memory and computational time to construct the system matrix and obtain its approximate inverse compared to a H matrix-accelerated VIE solver. The accuracy and efficiency of the proposed solver have been demonstrated via its application to the EM analysis of large-scale canonical and real-world structures comprising of broad permittivity values and involving millions of unknowns.

42 ENGINEERING↗
NASA NTRS2000

GENOA-PFA: Progressive Fracture in Composites Simulated Computationally

GENOA-PFA is a commercial version of the Composite Durability Structural Analysis (CODSTRAN) computer program that simulates the progression of damage ultimately leading to fracture in polymer-matrix-composite (PMC) material structures under various loading and environmental conditions. GENOA-PFA offers several capabilities not available in other programs developed for this purpose, making it preferable for use in analyzing the durability and damage tolerance of complex PMC structures in which the fiber reinforcements occur in two- and three-dimensional weaves and braids. GENOA-PFA implements a progressive-fracture methodology based on the idea that a structure fails when flaws that may initially be small (even microscopic) grow and/or coalesce to a critical dimension where the structure no longer has an adequate safety margin to avoid catastrophic global fracture. Damage is considered to progress through five stages: (1) initiation, (2) growth, (3) accumulation (coalescence of propagating flaws), (4) stable propagation (up to the critical dimension), and (5) unstable or very rapid propagation (beyond the critical dimension) to catastrophic failure. The computational simulation of progressive failure involves formal procedures for identifying the five different stages of damage and for relating the amount of damage at each stage to the overall behavior of the deteriorating structure. In GENOA-PFA, mathematical modeling of the composite physical behavior involves an integration of simulations at multiple, hierarchical scales ranging from the macroscopic (lamina, laminate, and structure) to the microscopic (fiber, matrix, and fiber/matrix interface), as shown in the figure. The code includes algorithms to simulate the progression of damage from various source defects, including (1) through-the-thickness cracks and (2) voids with edge, pocket, internal, or mixed-mode delaminations.

Murthy, Pappu L. N.↗
NASA NTRS1990

Efficient parallel solution of parabolic equations - Implicit methods on the Cedar multicluster

A class of implicit methods for the parallel solution of linear parabolic differential equations based on Pade and Chebyshev rational approximations to the matrix exponential are presented. It is pointed out that this approach incorporates both natural hierarchical parallelism, improved intrinsic efficiency, and fewer timesteps. These advantages lead to an extremely fast family of methods for the solution of certain time-dependent problems. These techniques are illustrated with numerical experiments on the University of Illinois Cedar multicluster architecture. The experiments indicate that implicit methods of very high degree offer great promise for the solution of certain parabolic problems when in computational environment with parallel resources. Hierarchically organized parallel computers, such as the Cedar multicluster, are found to be especially attractive for these schemes.

Gallopoulos, E.↗
NASA NTRS1993

Hierarchical Simulation of Hot Composite Structures

Computational procedures are described to simulate the thermal and mechanical behavior of high temperature metal matrix composites (HT-MMC) in the following three broad areas: (1) Behavior of HT-MMC's from micromechanics to laminate via Metal Matrix Composite Analyzer (METCAN), (2) tailoring of HT-MMC behavior for optimum specific performance via Metal Matrix Laminate Tailoring (MMLT), and (3) HT-MMC structural response for hot structural components via High Temperature Composite Analyzer (HITCAN). Representative results from each area are presented to illustrate the effectiveness of computational simulation procedures. The sample case results show that METCAN can be used to simulate material behavior such as strength, stress-strain response, and cyclic life in HTMMC's; MMLT can be used to tailor the fabrication process for optimum performance such as that for in-service load carrying capacity of HT-MMC's; and HITCAN can be used to evaluate static fracture and fatigue life of hot pressurized metal matrix composite rings.

Chamis, C. C.↗
DOE OSTI2020

Hierarchical morphologies in co-sputter deposited thin films

Co-depositions of immiscible alloy films at specific processing conditions have yielded hierarchical microstructures which consist of distinct features at multiple length scales, often agglomerates and concentration modulations on the sub-micrometer-scale and fine nanoprecipitates in a matrix on the nanoscale. The present work examined a series of immiscible alloy systems: Cu-Mo, Cu-Ag, Cu-Fe, Cu-Ta, Mo-Ag, Cu-Mo-Ag, to determine the kinetic conditions favorable for hierarchical organization and the formation mechanism of such structures. Thin films of six immiscible systems were sputter co-deposited over a range of deposition rates from 0.12 to 2 nm/s and various temperatures from 400 to 800°C. The resulting microstructures indicate that hierarchical structures form with sufficient disparity in kinetic energy between the constituent atoms, one species being highly mobile (A) and the other relatively immobile (B). This condition arises typically at elevated deposition temperatures and reduced deposition rates but is also alloy dependent. The hierarchical structures form during deposition via phase separation and self-organization processes across the multiple length scales. The adatoms diffuse on the film surface with the highly mobile species swiftly agglomerating into A-rich domains within which B-rich nanoprecipitates form, often self-organizing into periodic arrays. The smallest B-rich nanoprecipitates in the A-rich domains are found to be coherent and in a metastable crystal structure (B taking the structure of A), but coarser precipitates that exhibit the equilibrium structure of B element. The A-rich domains are surrounded by a B-rich matrix that phase-separates into a concentration modulated structure. In conclusion, the observations are interpreted via a model incorporating material properties and process parameters.

36 MATERIALS SCIENCE↗