Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density 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 1,513 records · Page 84

Benchmarking Mixed Mode Failure in Progressive Damage and Failure Analysis Methods

The verification and validation of progressive-damage-analysis finite element methods are difficult but critical tasks to undertake during their development. Verification exercises assess whether a predictive analysis tool produces results that are consistent with the fundamental concepts and assumptions of the tool under evaluation. Ideally, closed-form analytical solutions can be derived for which method verification results can be compared. Problems selected for computational tool verification are often simple and isolate individual features of the tool. In the case of progressive damage finite element methods, verifications should be performed to evaluate the ability of the model to predict the initiation of damage and its growth through the finite element mesh under a variety of conditions. Mabson et al. proposed a test case of a unidirectional, fiber-reinforced plate with a center crack subjected to tensile loads to evaluate matrix crack propagation predictions. The problem was modeled using the Abaqus Hashin continuum damage mechanics (CDM) model for fiber-reinforced composites. Different combinations of matrix strengths and element sizes were used in the simulations, and the results were compared to a closed-form solution based on linear elastic fracture mechanics (LEFM). It was determined that the Abaqus CDM model could predict the LEFM solution of Mode I cracks only when the finite element mesh density met specific requirements based on the material properties. This paper presents closed-form LEFM solutions for a center notch mixed mode (CNMM) verification problem. Parametric finite element analyses were developed using progressive damage analysis methods of both the Discrete Damage Mechanics (DDM) and CDM classes. The progressive damage analysis methods applied in the analyses of the CNMM problem include CompDam and the Floating Node Method. Analyses were conducted with various mode mixities and element sizes to verify that the damage models were working as intended and to identify any limits of applicability.

Leone, Frank↗

Inverse Calculation of Burden Distribution Matrix Using B-spline Model Based PDF control in Blast Furnace Burden Charging Process

The inverse calculation of burden distribution matrix (BDM) is one of the most important challenges in the blast furnace operation in iron-making processes. In general, blast furnace consumes 65% of the total energy for the whole steel-making. Focusing on this practical challenge, this article proposes a new burden distribution spatial model in calculating burden charging process, and develops a B-spline approximation-based probability density function (PDF) control algorithm to assign the expected thickness distribution of burden layer and, thus, develops a new method for the required inverse calculation of BDM. First, a novel method for the thickness distribution of burden layer is given using B-spline model to produce an expected distribution shape subjected to a desired tracking within a specific spatial constraint. Then, according to the coexistence of continuous and bounded discrete variables in BDM, a novel hybrid optimization control method by combining integer programming and PDF tracking is further established for the effective inverse calculation of BDM. Finally, the proposed PDF-based iterative inverse calculation of BDM using B-spline models are tested using various data from industrial examples. Furthermore, the simulation results show that the proposed method is well suited to solve the BDM inverse calculation problem in practice.

42 ENGINEERING↗

Developing Reaction Chemistry Models from Reactive Molecular Dynamics: TATB

Reactive Molecular Dynamic (RMD) are used to simulate the cook-off chemistry of TATB at a variety of fixed density and fixed temperature conditions. The chemical transformations are monitored using a Coordination Geometry Analysis (CGA) approach which tracks which atom types are bonded to each specific atom. This particularly identifies oxidation state changes that occur during the transformations. Correlations between these different chemical changes are identified using a Non-negative Matrix Factorization (NMF) approach. These identify reduced order chemistry models for the TATB system which contains six components whose concentration profiles are a function of both the temperature and density/pressure. The time histories of these transformations appear to show exponential growth/decay properties that could be fit with Arrhenius rates. These components should form the basis of deflagration rate models for these materials which could then be used in mesoscale simulations to analyze accidental initiation, shock-to-detonation and detonation propagation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of thermal treatments on electrochemical behavior of binder jetted 17-4 PH stainless steel

Binder jetted 17–4 PH stainless steel was post-processed to relate heat treatment, microstructure, and corrosion in 3.5 wt% NaCl. Specimens were sintered at 1380 or 1400 °C, solution-annealed at 1055 °C for 1 h, and aged at 482 °C for 1 h. Here, as-sintered parts showed α′-martensitic matrix with a δ-ferrite network and Cu-rich precipitates in ferrite; inclusions (MnS, NbC) promoted localized attack. solutionizing redistributed elements and reduced ferrite, while aging generated coherent Cu nano-precipitates. Corrosion resistance was highly sensitive to post-processing in which aged specimens exhibited the lowest corrosion current density and formed a thicker, stable Cr 2 O 3 -rich passive film, whereas sintered specimens degraded most. Pitting potential depended on sintering temperature and microstructure, with 1400 °C sintering yielding more positive pitting potentials and the best overall performance after aging. Although the 1400-solutionized condition showed a relatively noble pitting response versus 1380-solutionized, it displayed unstable corrosion kinetics attributed to an imperfect passive film linked to higher NbC density. XPS depth profiles corroborated these trends, showing thicker, more continuous Cr 2 O 3 in aged states and discontinuous/thinning oxides in less resistant conditions. Practically, high-temperature sintering (∼1400 °C) followed by solutionizing and aging is recommended, with further gains expected from reducing NbC/MnS populations and porosity via powder and process control.

36 MATERIALS SCIENCE↗

ElasTool v3.0: Efficient computational and visualization toolkit for elastic and mechanical properties of materials

Efficient computation and visualization of elastic and mechanical properties are crucial in the selection of materials and the design of new materials. Here, the ElasTool v3.0 toolkit marks a significant advancement in the computational analysis and visualization of elastic and mechanical properties of materials, essential in material selection and design. This enhanced version extends beyond standard calculations like elastic tensor, Young's modulus, bulk modulus, and Poisson's ratio. It introduces capabilities for computing minimum thermal conductivity, linear compressibility, rendering the Christoffel equation, and elastic energy density. Notably, it integrates advanced visualization tools, including compatibility with Plotly and Elate web platforms for interactive web-based property exploration. A key feature of ElasTool v3.0 is the implementation of second-order elastic constants (SOECs) for tubular 2D-based nanostructures and nanotubes. Leveraging high-efficiency strain-matrix sets (OHESS), the toolkit now facilitates efficient computation of elastic constants and mechanical properties at both zero and finite temperatures for 1D, 2D, and 3D dimensions. ElasTool is openly accessible on GitHub: https://github.com/gmp007/elastool.

1D, 2D, 3D, and tubular 2D nanostructure and nanot↗

Nanoscale Ion Transport Enhances Conductivity in Solid Polymer-Ceramic Lithium Electrolytes

The predictive design of flexible and solvent-free polymer electrolytes for solid-state batteries requires an understanding of the fundamental principles governing the ion transport. In this work, we establish a correlation among the composite structures, polymer segmental dynamics, and lithium ion (Li + ) transport in a ceramic-polymer composite. Elucidating this structure–property relationship will allow tailoring of the Li + conductivity by optimizing the macroscopic electrochemical stability of the electrolyte. The ion dissociation from the slow polymer segmental dynamics was found to be enhanced by controlling the morphology and functionality of the polymer/ceramic interface. The chemical structure of the Li + salt in the composite electrolyte was correlated with the size of the ionic cluster domains, the conductivity mechanism, and the electrochemical stability of the electrolyte. Polyethylene oxide (PEO) filled with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) or lithium bis(fluorosulfonyl) imide (LiFSI) salts was used as a matrix. A garnet electrolyte, aluminum substituted lithium lanthanum zirconium oxide (Al-LLZO) with a planar geometry, was used for the ceramic nanoparticle moieties. Further, the dynamics of the strongly bound and highly mobile Li + were investigated using dielectric relaxation spectroscopy. The incorporation of the Al-LLZO platelets increased the number density of more mobile Li + . The structure of the nanoscale ion-agglomeration was investigated by small-angle X-ray scattering, while molecular dynamics (MD) simulation studies were conducted to obtain the fundamental mechanism of the decorrelation of the Li + in the LiTFSI and LiFSI salts from the long PEO chain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vertical Architecture Solution-Processed Quantum Dot Photodetectors with Amorphous Selenium Hole Transport Layer

Colloidal quantum dots (CQDs) provide wide spectral tunability and high absorption coefficients owing to quantum confinement and large oscillator strengths, which along with solution processability, allow a facile, low-cost, and room-temperature deposition technique for the fabrication of photonic devices. However, many solution-processed CQD photodetector devices demonstrate low specific-detectivity and slow temporal response. To achieve improved photodetector characteristics, limiting carrier recombination and enhancing photogenerated carrier separation are crucial. In this study, we develop and present an alternate vertical-stack photodetector wherein we use a solution-processed quantum dot photoconversion layer coupled to an amorphous selenium (a-Se) wide-bandgap charge transport layer that is capable of exhibiting single-carrier hole impact ionization and is compatible with active-matrix readout circuitry. This a-Se chalcogenide transport layer enables the fabrication of high-performance and reliable solution-processed quantum dot photodetectors, with enhanced charge extraction capabilities, high specific detectivity (D* ~ 0.5–5 × 10 12 Jones), fast 3 dB electrical bandwidth (3 dB BW ~ 22 MHz), low dark current density (J D ~ 5–10 pA/cm 2 ), low noise current (i n ~ 20–25 fW/Hz 1/2 ), and high linear dynamic range (LDR ~130–150 dB) across the measured visible electromagnetic spectrum (~405–656 nm).

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Photo-actuators via epitaxial growth of microcrystal arrays in polymer membranes

Photomechanical crystals composed of three-dimensionally ordered and densely packed photochromes hold promise for high-performance photochemical actuators. However, bulk crystals with high structural ordering are severely limited in their flexibility, resulting in poor processibility and a tendency to fragment upon light exposure, while previous nano- or microcrystalline composites have lacked global alignment. Here we demonstrate a photon-fuelled macroscopic actuator consisting of diarylethene microcrystals in a polyethylene terephthalate host matrix. These microcrystals survive large deformations and show a high degree of three-dimensional ordering dictated by the anisotropic polyethylene terephthalate, which critically also has a similar stiffness. Overall, these ordered and compliant composites exhibit rapid response times, sustain a performance of over at least hundreds of cycles and generate work densities exceeding those of single crystals. In conclusion, our composites represent the state-of-the-art for photochemical actuators and enable properties unattainable by single crystals, such as controllable, reversible and abrupt jumping (photosalient behaviour).

36 MATERIALS SCIENCE↗

Impact of CH 4 addition on the electron properties and electric field dynamics in a Ar nanosecond-pulsed dielectric barrier discharge

Non-equilibrium plasmas derive their low temperature reactivity from producing and driving energetic electrons and active species under large electric fields. Therefore, the impact of reactants on the plasma properties including electron number density, electric field, and electron temperature is critical for applications such as plasma CH 4 reforming. Due to experimental complexity, electron properties and the electric field are rarely measured together in the same discharge. Here, we combine time-resolved Thomson scattering and electric field induced second harmonic generation (EFISH) to probe electron temperature, electron density, and electric field strength in a 60 Torr CH 4 /Ar nanosecond-pulsed dielectric barrier discharge (ns-DBD) while varying the CH 4 mole fraction from 0 to 8%. These measurements are compared to a 1-D numerical model to benchmark its predictions and identify areas of uncertainty. Nonlinear coupling between CH 4 addition, electron temperature, electron density, and the electric field was directly observed. Contrary to previous measurements in He, the electron temperature increased with CH 4 mole fraction. This rise in electron temperature is identified as electron heating by residual electric fields that increased with larger CH 4 mole fraction. Moreover, the electron number density has been found to decrease rapidly with the increase of methane mole fraction. Comparison of these measurements with the model yielded better agreement at higher CH 4 mole fractions and with the usage of ab initio calculated Ar electron-impact cross-sections from the B-spline R-matrix (BSR) database. Furthermore, the calculated plasma properties are shown to be sensitive to the residual surface charge implanted on the quartz dielectric surfaces. Without considering surface charge in the simulations, the calculated electric field profiles agreed well with the measurements, but the electron properties were underpredicted by more than a factor of three. Therefore, measurements of either the electric field or electron properties measurements alone are insufficient to fully validate modelling predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Report on partially complete post-irradiation-examination of the INL samples

Ferritic-martensitic steels (NF616 and T91), and austenitic stainless steel 800H and its Grain Boundary Engineering (GBE)-treated version 800H-TMP (ThermoMechanical Processing) were irradiated in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) and their microstructures and mechanical properties were characterized. Ferritic-martensitic steels were irradiated up to 8.16 dpa in the irradiation temperatures ranging from 241 °C to 447.5 °C. Selected 800H and 800H-TMP samples were also irradiated up to 9.12 dpa at 359°C to 431°C. Vickers hardness measurement, fractography, and microstructural characterization were performed on the selected samples in the LAMDA lab. Radiation-induced hardening of NF616 and T91 was observed in the hardness measurements which show ~37% to ~65% increases, depending on irradiation dose and irradiation temperature. Within the irradiation conditions studied in this work, samples with a higher dose had a higher hardness. All the tested alloy 800H and 800H-TMP in this work showed irradiation-hardening by ~96±7% to ~152±10%. Alloy 800HTMP tended to have slightly smaller irradiation-hardening than alloy 800H. The fractography of NF616 samples D2 (2.96 dpa at 291.5°C), D4 (5.91 dpa at 359°C), and D6 (8.16 dpa at 431°C) indicated loss of ductility with negligible necking for the sample D2, yet obvious ductile failure for samples D4 and D6. Fractography of alloy 800H samples N4 (7.27 dpa at 359°C), N5 (3.9 dpa at 451.5°C), and N6 (9.01 dpa at 431°C), and alloy 800H-TMP samples P4 (7.36 dpa at 359°C), P5 (3.95 dpa at 451.5°C), and P6 (9.12 dpa at 431°C) indicated that all samples failed in a ductile mode, with obvious necking. Dimples were observed where large Ti-containing particles were observed for the samples. TEM characterization was performed on alloy 800H (N4, N5, N6) and 800H-TMP (P4, P5 and P6). Accumulation of M23C6 precipitates at grain boundaries was observed in all samples, and the presence of Ti(C,N) precipitates at grain boundaries and in the matrix was observed in alloy 800H-TMP. M23C6 precipitates were found to have the cube-on-cube orientation relationship with the matrix. In addition, nanoscale Si-rich clusters were observed in the matrix of all samples. Atom probe tomography was conducted on the same samples, and their results are being analyzed to be integrated with the TEM results for a confident description of the clustering behavior. Due to the low irradiation temperature (up to 451.5°C), dislocation loops formed in all samples, with the density and the average size of dislocation loops in the order of 1022 – 1023 m-3 and 11.7 – 15.9 nm, respectively. The loop densities in alloy 800H were higher than that in alloy 800H-TMP under the same irradiation conditions. Further systematic data analyses, together with some complementary experiments, will be pursued for these samples to foster peer-reviewed journal article publications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Carbon-Carbon Piston Architectures

An improved structure for carbon-carbon composite piston architectures is disclosed. The improvement consists of replacing the knitted fiber, three-dimensional piston preform architecture described in U.S. Pat.No. 4,909,133 (Taylor et al.) with a two-dimensional lay-up or molding of carbon fiber fabric or tape. Initially, the carbon fabric of tape layers are prepregged with carbonaceous organic resins and/or pitches and are laid up or molded about a mandrel, to form a carbon-fiber reinforced organic-matrix composite part shaped like a "U" channel, a "T"-bar, or a combination of the two. The molded carbon-fiber reinforced organic-matrix composite part is then pyrolized in an inert atmosphere, to convert the organic matrix materials to carbon. At this point, cylindrical piston blanks are cored from the "U"-channel, "T"-bar, or combination part. These blanks are then densified by reimpregnation with resins or pitches which are subsequently carbonized. Densification is also accomplished by direct infiltration with carbon by vapor deposition processes. Once the desired density has been achieved, the piston billets are machined to final piston dimensions; coated with oxidation sealants; and/or coated with a catalyst. When compared to conventional steel or aluminum alloy pistons, the use of carbon-carbon composite pistons reduces the overall weight of the engine; allows for operation at higher temperatures without a loss of strength; allows for quieter operation; reduces the heat loss; and reduces the level of hydrocarbon emissions.

Rivers, H. Kevin↗

Development and Characterization of Embedded Sensory Particles Using Multi-Scale 3D Digital Image Correlation

A method for detecting fatigue cracks has been explored at NASA Langley Research Center. Microscopic NiTi shape memory alloy (sensory) particles were embedded in a 7050 aluminum alloy matrix to detect the presence of fatigue cracks. Cracks exhibit an elevated stress field near their tip inducing a martensitic phase transformation in nearby sensory particles. Detectable levels of acoustic energy are emitted upon particle phase transformation such that the existence and location of fatigue cracks can be detected. To test this concept, a fatigue crack was grown in a mode-I single-edge notch fatigue crack growth specimen containing sensory particles. As the crack approached the sensory particles, measurements of particle strain, matrix-particle debonding, and phase transformation behavior of the sensory particles were performed. Full-field deformation measurements were performed using a novel multi-scale optical 3D digital image correlation (DIC) system. This information will be used in a finite element-based study to determine optimal sensory material behavior and density.

Cornell, Stephen R.↗

Scalable, Infiltration-Free Ceramic Matrix Composite (SIF-CMC) manufacturing for Molten Salt Receiver: Investigation of particle-loaded preceramic resin

Due to their unique combination of high temperature, corrosion resistance, and toughness, Ceramic Matrix Composites (CMC) are the optimal materials to be used in the next generation of Concentrated Solar Power (CSP), which will operate at working temperatures up to 800 °C. However, today, CMC is too costly to be used in CSP applications due to its complex manufacturing process based on repeated infiltration-pyrolysis cycles. To solve this issue, PARC is developing a scalable, infiltration-free method (SIF) to manufacture CMC, which can reduce the production cost, increase batch consistency, and leverage all of the intrinsic advantages of current CMCs. The key to the SIF method is using a flowable preceramic resin with high loading of functional components. The formulation impregnates the fiber during curing and leads to high-density, low-shrinkage parts upon pyrolysis. In this paper, we demonstrate a proof-of-concept of the SIF process and report the relationships between key process parameters in the manufacturing process: functional material loading in the resin vs. resin’s flowability and morphology changes upon pyrolysis.

ceramic matrix composite, concentrated solar power↗

Quantitative assessment of fitting errors associated with streak camera noise in Thomson scattering data analysis

Thomson scattering measurements in high energy density experiments are often recorded using optical streak cameras. In the low-signal regime, noise introduced by the streak camera can become an important and sometimes the dominant source of measurement uncertainty. In this paper, we present a formal method of accounting for the presence of streak camera noise in our measurements. We present a phenomenological description of the noise generation mechanisms and present a statistical model that may be used to construct the covariance matrix associated with a given measurement. This model is benchmarked against simulations of streak camera images. We demonstrate how this covariance may then be used to weight fitting of the data and provide quantitative assessments of the uncertainty in the fitting parameters determined by the best fit to the data and build confidence in the ability to make statistically significant measurements in the low-signal regime, where spatial correlations in the noise become apparent. These methods will have general applicability to other measurements made using optical streak cameras.

47 OTHER INSTRUMENTATION↗

Experimental evaluation of localized creep deformation in grade 91 steel weldments

Spatially resolved measurement of localized creep deformation in heterogeneous creep resistant steel weldments is crucial but challenging for lifetime assessments of critical steam components in power plants. Here, experimental approaches were established to quantitatively evaluate commonly observed localized creep deformation in multi-pass Grade 91 steel weldments. An in-situ digital image correlation (DIC) system was utilized with a creep testing frame to monitor and measure both full-field strain and localized strain accumulation across the weldments during long-term creep testing at elevated temperatures. The in situ DIC method measured not only the creep deformation behavior of the weld metal, heat affected zone (HAZ), and base metal, but also creep strain evolution for each sub-region within the HAZ itself, including coarse-grained HAZ, fine-grained HAZ (FGHAZ), and intercritical HAZ (ICHAZ). The DIC results revealed that local creep strain in the ICHAZ reached up to 90% strain before final rupture, whereas nominal creep strain measured by the standard extensometer of the tested cross-weld specimen was below 10%, indicative of Type IV cracking of the Grade 91 weld. Microstructural analyses revealed that the faster creep degradation/deformation in the HAZ was caused mainly by accelerated matrix grain recrystallization/growth and a reduced pinning effect from the segregated and coarsened precipitates in the FGHAZ and ICHAZ. The ultimate creep rupture occurred in the ICHAZ owing to its lowest creep resistance induced by the largest recrystallized grain size, the lowest fraction of coincidence site lattice, and the lowest local strain energies/dislocation densities.

36 MATERIALS SCIENCE↗

Data assimilation in operator algebras

We develop an algebraic framework for sequential data assimilation of partially observed dynamical systems. In this framework, Bayesian data assimilation is embedded in a nonabelian operator algebra, which provides a representation of observables by multiplication operators and probability densities by density operators (quantum states). In the algebraic approach, the forecast step of data assimilation is represented by a quantum operation induced by the Koopman operator of the dynamical system. Moreover, the analysis step is described by a quantum effect, which generalizes the Bayesian observational update rule. Projecting this formulation to finite-dimensional matrix algebras leads to computational schemes that are i) automatically positivity-preserving and ii) amenable to consistent data-driven approximation using kernel methods for machine learning. Moreover, these methods are natural candidates for implementation on quantum computers. Applications to the Lorenz 96 multiscale system and the El Niño Southern Oscillation in a climate model show promising results in terms of forecast skill and uncertainty quantification.

97 MATHEMATICS AND COMPUTING↗

Entanglement and many-body effects in collective neutrino oscillations

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings, such as supernovae, where the neutrino density is large. In this regime, neutrino-neutrino interactions are important and simulations in the mean-field approximation show evidence for collective oscillations occurring at timescales much shorter than those associated with vacuum oscillations. In this work, we study the out-of-equilibrium dynamics of a corresponding spin model using matrix product states and show how collective bipolar oscillations can be triggered by many-body correlations if appropriate initial conditions are present. We find entanglement entropies scaling at most logarithmically in the system size suggesting that classical tensor network methods could be efficient in describing collective neutrino dynamics more generally. These observation provide a clear path forward, not only to increase the accuracy of current simulations, but also to elucidate the mechanism behind collective flavor oscillations without resorting to the mean-field approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mechanical properties of monolithic and particulate composites of L12 forms of Al3Ti

Preliminary studies of simple compressive properties of particulate composites with L1(2)-modified Al3Ti's as the matrices were performed as functions of matrix composition, amount of TiB2, processing technique, temperature, and strain rate. At and below about 900 K the introduction of TiB2 particles overwhelms minor chemical effects, and extremely high-yield strengths (about 2000 MPa at the 20 vol pct TiB2 level) were obtained. As L1(2)-modified Al3Ti has about half the density of the Ni-base superalloys, such values translate to strength/density ratios approximately four times those for superalloys. This advantage is not maintained by the Al3Ti-based composites above about 1000 K, particularly under slow strain rate plastic flow conditions, and it is proposed that this behavior could be due to the dissolution of intermetallic second phases as well as grain boundary creep mechanisms.

Whittenberger, J. D.↗