SEARCH · Engineering Papers

Results for “dynamic matrix completion”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1989

A review of some problems in global-local stress analysis

The various types of local-global finite-element problems point out the need to develop a new generation of software. First, this new software needs to have a complete analysis capability, encompassing linear and nonlinear analysis of 1-, 2-, and 3-dimensional finite-element models, as well as mixed dimensional models. The software must be capable of treating static and dynamic (vibration and transient response) problems, including the stability effects of initial stress, and the software should be able to treat both elastic and elasto-plastic materials. The software should carry a set of optional diagnostics to assist the program user during model generation in order to help avoid obvious structural modeling errors. In addition, the program software should be well documented so the user has a complete technical reference for each type of element contained in the program library, including information on such topics as the type of numerical integration, use of underintegration, and inclusion of incompatible modes, etc. Some packaged information should also be available to assist the user in building mixed-dimensional models. An important advancement in finite-element software should be in the development of program modularity, so that the user can select from a menu various basic operations in matrix structural analysis.

Nelson, Richard B.↗

NASA NTRS1991

Solar dynamic heat receiver technology

A full-size, solar dynamic heat receiver was designed to meet the requirements specified for electrical power modules on the U.S. Space Station, Freedom. The heat receiver supplies thermal energy to power a heat engine in a closed Brayton cycle using a mixture of helium-xenon gas as the working fluid. The electrical power output of the engine, 25 kW, requires a 100 kW thermal input throughout a 90 minute orbit, including when the spacecraft is eclipsed for up to 36 minutes from the sun. The heat receiver employs an integral thermal energy storage system utilizing the latent heat available through the phase change of a high-temperature salt mixture. A near eutectic mixture of lithium fluoride and calcium difluoride is used as the phase change material. The salt is contained within a felt metal matrix which enhances heat transfer and controls the salt void distribution during solidification. Fabrication of the receiver is complete and it was delivered to NASA for verification testing in a simulated low-Earth-orbit environment. This document reviews the receiver design and describes its fabrication history. The major elements required to operate the receiver during testing are also described.

Sedgwick, Leigh M.↗

DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Conformal collider physics meets LHC data

The remarkably high energies of the Large Hadron Collider (LHC) have allowed for the first measurements of the shapes and scalings of multipoint correlators of energy flow operators, ⟨ Ψ | E ( n → 1 ) E ( n → 2 ) ⋯ E ( n → k ) | Ψ ⟩ , providing new insights into the Lorentzian dynamics of quantum chromodynamics (QCD). In this letter, we use recent advances in effective field theory to derive a rigorous factorization theorem for the light-ray density matrix, ρ = | Ψ ⟩ ⟨ Ψ | , inside high transverse momentum jets at the LHC. Using the light-ray operator product expansion, the scaling behavior of multipoint correlators can be computed from the expectation value of the twist-2 spin- J light-ray operators, O [ J ] , in this state, Tr [ ρ O [ J ] ] . We compute the light-ray density matrix at next-to-leading order, and combine this with results for the next-to-leading logarithmic scaling behavior of the correlators up to six-points, comparing with CMS open data. This theoretical accuracy allows us to resolve the quantum scaling dimensions of QCD light-ray operators inside jets at the LHC. Our factorization theorem for the light-ray density matrix at the LHC completes the link between recent developments in the study of energy correlators and LHC phenomenology, opening the door to a wide variety of precision jet substructure studies. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Shock compaction heating and collisional processes in the production of type 3 ordinary chondrites: Lessons from the (nearly) unique L3 chondrite melt breccia Northwest Africa 8709

Abstract Northwest Africa (NWA) 8709 is a rare example of a type 3 ordinary chondrite melt breccia and provides critical information for the shock compaction histories of chondrites. An L3 protolith for NWA 8709 is inferred on the basis of oxygen isotope composition, elemental composition, diverse mineral chemistry, and overall texture. NWA 8709 is among the most strongly shocked type 3 chondrites known, and experienced complete melting of the matrix and partial melting of chondrules. Unmelted phases underwent FeO reduction and partial homogenization, with reduction possibly occurring by reaction of olivine and low‐Ca pyroxene with an S‐bearing gas that was produced by vaporization. Chondrules and metal grains became foliated by uniaxial compaction, and during compression, chondrules and fragments became attached to form larger clumps. This process, and possibly also melt incorporation into chondrules to cause “inflation,” may have contributed to anomalously large chondrule sizes in NWA 8709. The melt breccia character is attributed to strong shock affecting a porous precursor. Data‐model comparisons suggest that a precursor with 23% porosity that was impacted by a 3 km/s projectile could have produced the meteorite. The rarity of other type 3 ordinary chondrite melt breccias implies that the immediate precursors to such chondrites were lower in porosity than the NWA 8709 precursor, or experienced weaker shocks. Altogether, the data imply a predominantly “quiet” dynamical environment to form most type 3 ordinary chondrites, with compaction occurring in a series of relatively weak shock events.

Ruzicka, Alex M.↗

DOE OSTI2022

Self-consistency in GWΓ formalism leading to quasiparticle-quasiparticle couplings

Within many-body perturbation theory, Hedin's formalism offers a systematic way to iteratively compute the self-energy Σ of any dynamically correlated interacting system, provided one can evaluate the interaction vertex Γ exactly. This is, however, impossible, in general, for it involves the functional derivative of Σ with respect to the Green's function. Here, we analyze the structure of this derivative, splitting it into four contributions and outlining the type of quasiparticle interactions that each of them generate. Moreover, we show how, in the implementation of self-consistency, the action of these contributions can be classified into two: A quantitative renormalization of previously included interaction terms and the inclusion of qualitatively distinct interaction terms through successive functional derivatives of Γ itself. Implementing this latter type of self-consistency can extend the validity of perturbative approximations based on Hedin's equations toward the high interaction limit, as we show in the example of the Hubbard dimer. Furthermore, our analysis also provides a unifying perspective on the perturbation theory landscape, showing how the T-matrix approach is completely contained in Hedin's formalism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2017

Incorporation of Failure Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in various coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

Impact↗

NASA NTRS2001

Seismic Study of the Solar Interior: Inferences from SOI/MDI Observations During Solar Activity

We have continued in collaboration with Dr. Eff-Darwich (University of La Laguna, Tenerife, Spain) the study of the structure, asphericity and dynamics of the solar interior from p-mode frequencies and frequency splittings. In March 2001, Dr. Eff-Darwich came for 3 weeks visit to CfA. During this visit we completed our work on the inversion of the internal solar rotation rate, and submitted a paper describing this work to the Astrophysical Journal. This paper has been recently revised in response to the referee comments and I expect that it will be accepted for publication very soon. We also have analyzed helioseismic data looking for temporal variations of the solar stratification near the base of the convection zone. We have expanded on the initial work that was presented at the SOHO-10/GONG-2000 meeting (October 2000, Tenerife), and are in the process of writing this up. Substantial progress towards the characterization of high-degree p-modes has been achieved. Indeed, in collaboration Dr. Rabello-Soares (Stanford University), we have gained a clear conceptual understanding of the various elements that affect the leakage matrix of the SOI/MDI instrument. This was presented in an invited talk at the SOHO-10/GONG-2000 meeting (October 2000, Tenerife). Once we will have successfully migrated from a qualitative to a quantitative assessment of these effects, we should be able to generate high-degree p-modes frequencies so crucial in the diagnostic of the layers just below solar surface.

Korzennik, Sylvain G.↗

NASA NTRS2000

Monitoring Astronaut Health at the Nanoscale Cellular Level Through the Eye

A user friendly goggles-like head-mounted device equipped with a suite of instruments for several non-invasive and quantitative medical evaluation of the eye, skin, and brain is desired for monitoring the health of astronauts during space travel and exploration of neighboring and distant planets. Real-time non-invasive evaluation of the different structures within the above organs can provide indices of the health of not just these organs, but the entire body. The techniques such as dynamic light scattering (for the early detection of uveitis, cholesterol levels, cataract, changes in the vitreous and possibly Alzheimer's disease), corneal autofluorescence (to assess extracellular matrix biology e.g., in diabetes), optical activity measurements (of anterior ocular fluid to evaluate blood-glucose levels), laser Doppler velocimetry (to assess retinal, optic nerve, and choroidal blood flow), reflectometry/oximetry (for assessing ocular and central nervous system oxygen metabolism), optical coherence tomography (to determine retinal tissue microstructure) and possibly scanning laser technology (for intraocular tissue imaging and scanning) will he integrated into this compact device. Skin sensors will also be mounted on the portion of the device in contact with the periocular region. This will enable monitoring of body temperature, EEG, and electrolyte status. This device will monitor astronaut health during long-duration space travel by detecting aberrations from pre-established "nonns", enabling prompt diagnosis and possibly the initiation of early preventative/curative therapy. The non-invasive nature of the device technologies permits frequent repetition of tests, enabling real-time complete crew health monitoring. This device may ultimately be useful in tele-medicine to bring modern healthcare to under-served areas on Earth as well as in so-called "advanced" care settings (e.g. diabetes in the USA).

Ansari, Rafat R.↗

DOE OSTI2024

The accuracies of effective interactions in downfolding coupled-cluster approaches for small-dimensionality active spaces

Here, this paper evaluates the accuracy of the Hermitian form of the downfolding procedure using the double unitary coupled cluster (DUCC) ansatz on the benchmark systems of linear chains of hydrogen atoms, H6 and H8. The computational infrastructure employs the occupation-number-representation codes to construct the matrix representation of arbitrary second-quantized operators, allowing for the exact representation of exponentials of various operators. The tests demonstrate that external amplitudes from standard single-reference coupled cluster methods that sufficiently describe external (out-of-active-space) correlations reliably parameterize the Hermitian downfolded effective Hamiltonians in the DUCC formalism. The results show that this approach can overcome the problems associated with losing the variational character of corresponding energies in the corresponding SR-CC theories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2005

A Fundamental Mathematical Model of a Microbial Predenitrification System

Space flight beyond Low Earth Orbit requires sophisticated systems to support all aspects of the mission (life support, real-time communications, etc.). A common concern that cuts across all these systems is the selection of information technology (IT) methodology, software and hardware architectures to provide robust monitoring, diagnosis, and control support. Another dimension of the problem space is that different systems must be integrated seamlessly so that communication speed and data handling appear as a continuum (un-interrupted). One such team investigating this problem is the Advanced Integration Matrix (AIM) team whose role is to define the critical requirements expected of software and hardware to support an integrated approach to the command and control of Advanced Life Support (ALS) for future long-duration human space missions, including permanent human presence on the Moon and Mars. A goal of the AIM team is to set the foundation for testing criteria that will assist in specifying tasks, control schemes and test scenarios to validate and verify systems capabilities. This project is to contribute to the goals of the AIM team by assisting with controls planning for ALS. Control for ALS is an enormous problem it involves air revitalization, water recovery, food production, solids processing and crew. In more general terms, these systems can be characterized as involving both continuous and discrete processes, dynamic interactions among the sub-systems, nonlinear behavior due to the complex operations, and a large number of multivariable interactions due to the dimension of the state space. It is imperative that a baseline approach from which to measure performance is established especially when the expectation for the control system is complete autonomous control.

Hoo, Karlene A.↗

DOE OSTI2022

The Effect of Geometry, Spin, and Orbital Optimization in Achieving Accurate, Correlated Results for Iron–Sulfur Cubanes

Iron–sulfur clusters comprise an important functional motif in the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from correlated, large active space calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. Here, to provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of the [Fe(II/III) 4 S 4 (SMe) 4 ] -2 clusters using two complementary, correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (<1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the singlet and the spin state 2S + 1 = 19, corresponding to a maximal number of Fe-d electrons being unpaired and of parallel spin, is strongly geometry dependent, changing by a factor of 3 upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap compared to typical mean-field orbital approaches. Our results show the need for performing large active space calculations to unveil the challenging electronic structure of these complex catalytic centers and should serve as accurate starting points for fully correlated treatments upon inclusion of dynamical correlation outside the active space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Semi-Annual Report for Horizontal Compact High Temperature Gas Reactor (HC-HTGR) Development during Performance Period April 2023 – September 2023

Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed HC-HTGR in 3 years and support its commercialization as a safe, low-cost HTGR. Argonne National Laboratory (Argonne) is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. Additionally, Argonne is providing analysis of the primary coolant system to ensure temperatures within the core remain below safety margins during steady-state and potential accident scenarios. This fourth semi-annual report summarized the progress made at Argonne on the two tasks during the second half of FY23. As a part of the RCCS design task, recent efforts have been made to complete a conceptual design of the RCCS for the HC-HTGR, including the design update of the water panel and system configuration favorable in point of view of fabrication and system operation. Design calculations were conducted under various heat load conditions to validate the system design. Transient simulations using RELAP5-3D were conducted to investigate system dynamics under transients of interest and to evaluate the system performance in the design condition. The results demonstrated the overall system feasibility that the RCCS design maintains structures temperatures lower than maximum allowable temperature with sufficient system inventory without any active heat sink. In the primary system thermal hydraulics task, the preliminary analysis was performed for a long term pressurized conduction cooldown (PCC) transient. This analysis used a combination of a fully resolved and coarse homogenized mesh to predict the temperature distribution for the steady-state initial condition and the PCC transient. The steady-state initial condition was determined using a fully resolved full core model with 3D solid to 1D fluid coupling. The fully resolved mesh was also used to model the first 20 seconds of the PCC. The temperature difference between fuel pins and the graphite matrix becomes minimal and the dominant heat transfer shifts to a larger scale radially towards the RCCS. After 20 seconds, a homogenized coarse mesh is used, greatly reducing the computational costs of the model. These results demonstrate that the core is designed to passively remove enough decay heat in a protected loss of primary coolant flow to prevent an unsafe rise of core temperatures.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

On the synthesis of resonance lines in dynamical models of structured hot-star winds

We examine basic issues involved in synthesizing resonance-line profiles from 1-D, dynamical models of highly structured hot-star winds. Although these models exhibit extensive variations in density as well as velocity, the density scale length is still typically much greater than the Sobolev length. The line transfer is thus treated using a Sobolev approach, as generalized by Rybicki & Hummer (1978) to take proper account of the multiple Sobolev resonances arising from the nonmonotonic velocity field. The resulting reduced-lambda-matrix equation describing nonlocal coupling of the source function is solved by iteration, and line profiles are then derived from formal solution integration using this source function. Two more approximate methods that instead use either a stationary or a structured, local source function yield qualitatively similar line-profiles, but are found to violate photon conservation by 10% or more. The full results suggest that such models may indeed be able to reproduce naturally some of the qualitative properties long noted in observed UV line profiles, such as discrete absorption components in unsaturated lines, or the blue-edge variability in saturated lines. However, these particular models do not yet produce the black absorption troughs commonly observed in saturated lines, and it seems that this and other important discrepancies (e.g., in acceleration time scale of absorption components) may require development of more complete models that include rotation and other 2-D and/or 3-D effects.

Puls, J.↗

NASA NTRS1993

Fluctuations at the blue edge of saturated wind lines in IUE spectra of O-type stars

We examine basic issues involved in synthesizing resonance-line profiles from 1-D, dynamical models of highly structured hot-star winds. Although these models exhibit extensive variations in density as well as velocity, the density scale length is still typically much greater than the Sobolev length. The line transfer is thus treated using a Sobolev approach, as generalized by Rybicki & Hummer (1978) to take proper account of the multiple Sobolev resonances arising from the nonmonotonic velocity field. The resulting reduced-Lambda-matrix equation describing nonlocal coupling of the source function is solved by iteration, and line profiles and then derived from formal solution integration using this source function. The more appropriate methods that instead use either a stationary or a structured, local source function yield qualitatively similar line-profiles, but are found to violate photon conservation by 10 percent or more. The full results suggest that such models may indeed be able to reproduce naturally some of the qualitative properties long noted in observed UV line profiles, such as discrete absorption components in unsaturated lines, or the blue-edge variability in saturated lines. However, these particular models do not yet produce the black absorption troughs commonly observed in saturated lines, and it seems that this and other important discrepancies (e.g., in acceleration time scale of absorption components) may require development of more complete models that include rotation and other 2-D and/or 3-D effects.

Owocki, Stanley P.↗

NASA NTRS2008

Post Test Evaluation of HSCT Nozzle Acoustic Liner Subcomponents Subjected to a Hot Acoustic Durability Test

The acoustic liner system designed for use in the High Speed Civil Transport (HSCT) was tested in a thermal-acoustic environment. Five ceramic matrix composite (CMC) acoustic tile configurations, five bulk acoustic absorbers, and one thermal protection system design were tested. The CMC acoustic tiles were subjected to two 2 3/4 hr ambient temperature acoustic exposures to measure their dynamic response. One exposure was conducted on the tiles alone and the second exposure included the tiles and the T-foam bulk absorber. The measured tile RMS strains were small. With or without the T-foam absorber, the dynamic strains were below strain levels that would cause damage during fatigue loading. After the ambient exposure, a 75-hr durability test of the entire acoustic liner system was conducted using a thermal-acoustic cycle that approximated the anticipated service cycle. Acoustic loads up to 139 dB/Hz and temperatures up to 1670 F (910 C) were employed during this 60 cycle test. During the durability test, the CMC tiles were exposed to temperatures up to 1780 F and a transient through thickness gradient up to 490 F. The TPS peak temperatures on the hot side of the panels ranged from 750 to 1000 F during the 60 cycles. The through thickness delta T ranged from 450 to 650 F, varying with TPS location and cycle number. No damage, such as cracks or chipping, was observed in the CMC tiles after completion of the testing. However, on tile warped during the durability test and was replaced after 43 or 60 cycles. No externally observed damage was found in this tile. No failure of the CMC fasteners occurred, but damage was observed. Cracks and missing material occurred, only in the fastener head region. No indication of damage was observed in the T-foam acoustic absorbers. The SiC foam acoustic absorber experienced damage after about 43 cycles. Cracking in the TPS occurred around the attachment holes and under a vent. In spite of the development of damage, the TPS maintained its insulative capability throughout the durability test. The durability test results demonstrate damage-tolerant CMC tile, CMC fastener, TPS, and T-foam absorber designs for the combined thermal and acoustic engine nozzle environment.

Verrilli, Michael J.↗

DOE OSTI2023

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗