Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix models”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Interactive Reliability Model for Whisker-toughened Ceramics

Wider use of ceramic matrix composites (CMC) will require the development of advanced structural analysis technologies. The use of an interactive model to predict the time-independent reliability of a component subjected to multiaxial loads is discussed. The deterministic, three-parameter Willam-Warnke failure criterion serves as the theoretical basis for the reliability model. The strength parameters defining the model are assumed to be random variables, thereby transforming the deterministic failure criterion into a probabilistic criterion. The ability of the model to account for multiaxial stress states with the same unified theory is an improvement over existing models. The new model was coupled with a public-domain finite element program through an integrated design program. This allows a design engineer to predict the probability of failure of a component. A simple structural problem is analyzed using the new model, and the results are compared to existing models.

Palko, Joseph L.↗
NASA NTRS1993

Establishing the relationship between manufacturing and component performance in stretch formed thermoplastic composites

Flexible manufacturing methods are needed to reduce the cost of using advanced composites in structural applications. One method that allows for this is the stretch forming of long discontinuous fiber materials with thermoplastic matrices. In order to exploit this flexibility in an economical way, a thorough understanding of the relationship between manufacturing and component performance must be developed. This paper reviews some of the recent work geared toward establishing this understanding. Micromechanics models have been developed to predict the formability of the material during processing. The latest improvement of these models includes the viscoelastic nature of the matrix and comparison with experimental data. A finite element scheme is described which can be used to model the forming process. This model uses equivalent anisotropic viscosities from the micromechanics models and predicts the microstructure in the formed part. In addition, structural models have been built to account for the material property gradients that can result from the manufacturing procedures. Recent developments in this area include the analysis of stress concentrations and a failure model each accounting for the heterogeneous material fields.

Santare, Michael H.↗
NASA NTRS1985

A General-applications Direct Global Matrix Algorithm for Rapid Seismo-acoustic Wavefield Computations

A new matrix method for rapid wave propagation modeling in generalized stratified media, which has recently been applied to numerical simulations in diverse areas of underwater acoustics, solid earth seismology, and nondestructive ultrasonic scattering is explained and illustrated. A portion of recent efforts jointly undertaken at NATOSACLANT and NORDA Numerical Modeling groups in developing, implementing, and testing a new fast general-applications wave propagation algorithm, SAFARI, formulated at SACLANT is summarized. The present general-applications SAFARI program uses a Direct Global Matrix Approach to multilayer Green's function calculation. A rapid and unconditionally stable solution is readily obtained via simple Gaussian ellimination on the resulting sparsely banded block system, precisely analogous to that arising in the Finite Element Method. The resulting gains in accuracy and computational speed allow consideration of much larger multilayered air/ocean/Earth/engineering material media models, for many more source-receiver configurations than previously possible. The validity and versatility of the SAFARI-DGM method is demonstrated by reviewing three practical examples of engineering interest, drawn from ocean acoustics, engineering seismology and ultrasonic scattering.

Schmidt, H.↗
NASA NTRS1997

The Influence of Correlated Crustal Signals in Modelling the Main Geomagnetic Field

Algorithms used in geomagnetic main-field modelling have for the most part treated the noise in the field measurements as if it were white. A major component of the noise consists of the field due to magnetization in the crust and it has been realized for some time that such signals are highly correlated at satellite altitude. Hence approximation by white noise, while of undoubted utility, is of unknown validity. In this paper we study two plausible statistical models for the crustal magnetization, in which the magnetization is a realization of a stationary, isotropic, random process. At a typical satellite altitude the associated fields exhibit significant correlation over ranges as great as 15 deg. or more, which introduces off-diagonal elements into the covariance matrix, elements that have usually been neglected in modelling procedures. Dealing with a full covariance matrix for a large data set would present a formidable computational challenge, but fortunately most of the entries in the covariance matrix are so small that they can be replaced by zeros. The resultant matrix comprises only about 3 per cent non-zero entries and thus we can take advantage of efficient sparse matrix techniques to solve the numerical system. We construct several main-field models based on vertical-component data from a selected 5 deg. by 5 deg. data set derived from the Magsat mission. Models with and without off-diagonal terms are compared.

Rygaard-Hjalsted, C.↗
DOE OSTI2025

Linear Solver for Electromagnetic Simulation of General Distribution Feeders

High-fidelity electromagnetic transient (EMT) modeling is required for accurate simulation and analysis of power system dynamics in modern distribution feeders. However, the high-fidelity of EMT models often leads to significant computational challenges, particularly in terms of computational resources and simulation time. This paper investigates the development and application of a detailed EMT model for general distribution feeders, with a focus on improving computational efficiency. A direct linear solver is proposed for a bordered block diagonal (BBD) matrix structure commonly encountered in a EMT model of distribution feeders. The solver integrates the Schur complement method with the block tridiagonal matrix algorithm to enhance the computational performance. The proposed solver is validated using the primary feeder of the IEEE 342-node test system, demonstrating its accuracy and efficiency in EMT simulations. Furthermore, the solver’s performance is benchmarked against MATLAB’s built-in linear solvers, showing significant improvements in computation time while maintaining high fidelity and accuracy in simulation results.

Choi, Jongchan [ORNL] (ORCID:000000025952455X)↗
DOE OSTI2026

Reservoir-scale model of geologic hydrogen production from serpentinization: Cyclic injection in a dual-permeability fracture-matrix system

Geologic hydrogen (GeoH 2 ) from serpentinization is a promising low-carbon resource, but its reservoir-scale behavior remains poorly understood. We develop a dual-permeability reactive transport model for an injector–producer well pair in ultramafic rock that couples multiphase flow, heat transfer, geochemistry, and porosity–permeability evolution. Here, the model is calibrated to olivine flow-through experiments and upscaled to two-year long simulations with continuous injection and cyclic injection with shut-in-to-injection ratios (SIR = 1, 0.5, 0.1). Olivine reacts along high-flux pathways to form lizardite and magnetite, increasing pH; H 2 (aq) and H 2 (g) peak early and then decline as exsolution and advective export outpace local generation. Continuous injection yields the highest cumulative H 2 but the lowest water-use efficiency (1.138 x 10 –6 mol/kgw). Cyclic injection increases this ratio to 1.35 x 10 –6 , 1.377 x 10 –6 , and 1.182 x 10 –6 mol/kgw for SIR = 1, 0.5, and 0.1, respectively; SIR = 0.5 provides a ~20% improvement over continuous injection and the best compromise for pilot design.

08 HYDROGEN↗
NASA NTRS1982

Micromechanical predictions of crack initiation, propagation and crack growth resistance in boron/aluminum composites

An elastoplastic, axisymmetric finite element model has been used to predict the initiation and propagation of a crack in a composite model consisting of a single broken boron fiber embedded in an annular sheath of aluminum matrix. The accuracy of the axisymmetric finite element model for crack problems has been established by solving the classical problem of a penny-shaped crack in a thick cylindrical rod under axial tension. Also, the stress intensity factors predicted by the present numerical model are compared with continuum results. A constant displacement boundary condition applied during an increment of crack growth permits a substantial amount of stable crack growth in the matrix material. The concept of Crack Growth Resistance Curves (KR-curves) has been used to determine the point of crack instability

Mahishi, J. M.↗
NASA NTRS2019

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↗
DOE OSTI2020

A quantitative phase-field model for gas bubble evolution in UO2

Due to the large formation energy of vacancies and noble gas atoms in the form of interstitials or substitutional atoms in nuclear fuel (UO2), the thermodynamic equilibrium concentrations of these species are very low in the nuclear fuel matrix even at very high temperatures, which imposes difficulties upon the quantitative study of bubble evolution via the phase-field method. In this study, a quantitative phase-field model is proposed to deal with this problem. The system’s free energy density is derived according to the principles of thermodynamics, with consideration of the elastic effect and with the use of real material parameters from experiments. The model is useful for the study of the kinetics of gas bubble growth with very dilute concentrations of vacancy and gas atoms in the matrix. This model is applied to study single bubble growth and multiple bubble growth under various concentrations of vacancy and gas atoms and at various temperatures. The elastic effect and the effects of the generation rate of vacancies and gas atoms on bubble growth are analyzed.

low defect concentration, nuclear fuel, gas bubble↗
DOE OSTI2020

Design of High Efficient Mid-Wavelength Infrared Polarizer on ORMOCHALC Polymer

While an organically modified chalcogenide (ORMOCHALC) can be used to fabricate a polymeric mid-wavelength infrared (MWIR) polarizer with competitive extinction ratio compared to the commercial wire-grid polarizers, which are composed of fragile inorganic materials, there is still a knowledge gap regarding the systematic design process to obtain high transmission efficiency and extinction ratio. As such, a computational parameter study for design optimization is conducted with the geometric parameters of the bilayer grating ORMOCHALC polarizer. The computational study shows that the Fabry–Pérot cavity is the primary mechanism that determines the transmission behaviors and the extinction ratio. A bilayer grating design, guided by the parameter study, is realized through the thermal nanoimprint and metal deposition processes. The extinction ratios measured with the Fourier-transform infrared are 245, 304, and 351 at the wavelength of 3, 4, and 5 μm, respectively. Compared to the state-of-the-art of the polymeric MWIR linear polarizers, the extinction ratio is improved by 1.4 times, and the transmission efficiency is increased by 2.5 times. Theoretical analysis with the multiple-layer model based on the transfer matrix method predicts a matched transmission behavior with the experiment and a full-wave electromagnetic simulation.

42 ENGINEERING↗
DOE OSTI2024

LSZ in action: extracting form factors from correlators nonperturbatively in 2d ϕ4 theory

Abstract In this paper, we compute multiparticle form factors of local operators in 2dϕ 4 theory using a recently proposed method [1] for efficiently implementing the LSZ prescription with Hamiltonian Truncation methods, and we adopt Lightcone Conformal Truncation (LCT) in particular for our calculations. We perform various checks of our results at weak and strong coupling, and elucidate the parametric behavior of truncation errors. This opens up the possibility to compute S-matrix in various strongly coupled models using the LSZ method in LCT.

Physics↗
DOE OSTI2020

Measurement of the cross section for $\text{t}\bar{\text{t}}$ production with additional jets and b jets in pp collisions at $\sqrt{s}=$ 13 TeV

Measurements of the cross section for the production of top quark pairs in association with a pair of jets from bottom quarks $ \left({\sigma}_{\mathrm{t}\overline{\mathrm{t}}\mathrm{b}\overline{\mathrm{b}}}\right) $ and in association with a pair of jets from quarks of any flavor or gluons $ \left({\sigma}_{\mathrm{t}\overline{\mathrm{t}}\mathrm{jj}}\right) $ and their ratio are presented. The data were collected in proton-proton collisions at a center-of-mass energy of 13 TeV by the CMS experiment at the LHC in 2016 and correspond to an integrated luminosity of 35.9 fb$^{−1}$. The measurements are performed in a fiducial phase space and extrapolated to the full phase space, separately for the dilepton and lepton+jets channels, where lepton corresponds to either an electron or a muon. The results of the measurements in the fiducial phase space for the dilepton and lepton+jets channels, respectively, are $ {\sigma}_{\mathrm{t}\overline{\mathrm{t}}\mathrm{jj}} $ = 2.36±0.02 (stat)±0.20 (syst) pb and 31.0±0.2 (stat)±2.9 (syst) pb, and for the cross section ratio 0.017 ± 0.001 (stat) ± 0.001 (syst) and 0.020 ± 0.001 (stat) ± 0.001 (syst). The values of $ {\sigma}_{\mathrm{t}\overline{\mathrm{t}}\mathrm{b}\overline{\mathrm{b}}} $ are determined from the product of the $ {\sigma}_{\mathrm{t}\overline{\mathrm{t}}\mathrm{jj}} $ and the cross section ratio, obtaining, respectively, 0.040±0.002 (stat)±0.005 (syst) pb and 0.62±0.03 (stat)±0.07 (syst) pb. These measurements are the most precise to date and are consistent, within the uncertainties, with the standard model expectations obtained using a matrix element calculation at next-to-leading order in quantum chromodynamics matched to a parton shower.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Unsupervised Clustering of Microseismic Events and Focal Mechanism Analysis at the CO 2 Injection Site in Decatur, Illinois

Characterization of induced microseismicity at a carbon dioxide (CO 2 ) storage site is critical for preserving reservoir integrity and mitigating seismic hazards. We apply a multilevel machine learning (ML) approach that combines the nonnegative matrix factorization and hidden Markov model to extract spectral representations of microseismic events and cluster them to identify seismic patterns at the Illinois Basin-Decatur Project. Unlike traditional waveform correlation methods, this approach leverages spectral characteristics of first arrivals to improve event classification and detect previously undetected planes of weakness. By integrating ML-based clustering with focal mechanism analysis, we resolve small-scale fault structures that are below the detection limits of conventional seismic imaging. Our findings reveal temporal bursts of microseismicity associated with brittle failure, providing insights into the spatio-temporal evolution of fault reactivation during CO 2 injection. This approach enhances seismic monitoring capabilities at CO 2 injection sites by improving fault characterization beyond the resolution of standard geophysical surveys.

Willis, Rachel Marie [Sandia National Laboratories↗
DOE OSTI2021

Two-color, intracavity pump–probe, cavity ringdown spectroscopy

Here, we report a proof-of-principle demonstration of intracavity pump–probe, cavity ringdown (CRD) detection in a three-mirror, traveling-wave cavity. With cavity-enhanced pump power and probe absorption path length, the technique is a generally applicable, high-sensitivity, high-selectivity detection method. In our experiments, the pump radiation is switched off during every other probe ringdown, which allows uncorrelated measurements of analyte and background cavity decay rates. The net, two-color signal from the difference between the pump-on and pump-off decay rates is immune to empty-CRD drifts and spectral overlaps from non-target molecular transitions. The immunity to the ringdown drifts allows longer signal-averaging and, thus, higher detection sensitivity. The ability to compensate for the background absorption enhances the detection selectivity in spectrally congested regions. Our technique is well-suited for trace-detection in the mid-IR region, where pump–probe schemes based on strong rovibrational transitions can be applied. In this work, two-color CRD detection is implemented on a ladder-type, three-level system based on the N 2 O, ν 3 = 1 ← 0, P(19) (pump) and ν 3 = 2 ← 1, R(18) (probe), rovibrational transitions. By frequency-locking two-quantum cascade lasers to the p-polarization (pump, Finesse = 5280) and s-polarization (probe, Finesse = 67 700) cavity modes, we achieve high intracavity pump power (36 W) and high probe ringdown rates (>2 kHz). The observed two-color spectra are simulated by a density-matrix, three-level system model that is solved under the constraints of the cavity resonance conditions. In addition to its background compensation capability, experimental flexibility in the selection of pump–probe schemes and signal insensitivity to intracavity laser power are further features that enhance the utility of our technique for mid-IR trace-detection.

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

Weighted FFT estimators for 1D and 3D correlations of the Lyman- α forest

Correlations in the Lyman-α (Lyα) forest, both as a function of line of sight separation (1D) and 3D separation, provide a unique window to the distribution of matter at redshifts not accessible by current galaxy surveys. While optimal quadratic estimators have been used to measure 1D correlations, they are computationally expensive and difficult to extend to 3D analyses. On the other hand, estimators based on the Fast Fourier Transform (FFT) are significantly faster, but are affected by missing data in the spectra (masked pixels) and so far have not used pixel weights to reduce the uncertainties in the measurement. In this publication we describe how to compute the window matrix that enables forward-modelling the impact of masked pixels and weights on the FFT-based estimators. Here, we use Gaussian and hydrodynamical simulations with artificially masked pixels to validate the method on the measurement of 1D correlations. Finally, we show that the formalism can be extended to model the impact on 3D correlations, in particular on the cross-spectrum, the correlation of 1D Fourier modes as a function of transverse separation. This work will enable more precise clustering measurements with the Lyα forest dataset recently collected by the Dark Energy Spectroscopic Instrument (DESI).

Lokken, Martine [Univ. Autonoma de Barcelona (Spai↗
DOE OSTI2020

Prediction of exotic magnetic states in the alkali-metal quasi-one-dimensional iron selenide compound Na 2 FeSe 2

The magnetic and electronic phase diagram of a model for the quasi-one-dimensional alkali-metal iron selenide compound Na 2 FeSe 2 is presented. The novelty of this material is that the valence of iron is Fe 2+ , contrary to most other iron-chain compounds with valence Fe 3+ . Using first-principles techniques, we developed a three-orbital tight-binding model that reproduces the ab initio band structure near the Fermi level. Including Hubbard and Hund couplings and studying the model via the density-matrix renormalization group and Lanczos methods, we constructed the ground-state phase diagram. A robust region where the block state ↑↑↓↓↑↑↓↓ is stabilized was unveiled. The analog state in iron ladders, employing 2×2 ferromagnetic blocks, is by now well established, but in chains a block magnetic order has not been observed yet in real materials. The phase diagram also contains a large region of canonical staggered spin order ↑↓↑↓↑↓↑ at very large Hubbard repulsion. At the block-to-staggered transition region, an exotic phase is stabilized with a mixture of both states: an inhomogeneous orbital-selective charge density wave with the exotic spin configuration ↑↑↓↑↓↓↑↓. Our predictions for Na 2 FeSe 2 may guide crystal growers and neutron-scattering experimentalists towards the realization of block states in one-dimensional iron selenide chain materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Peierls transition, ferroelectricity, and spin-singlet formation in monolayer VOI 2

Using ab initio density functional theory and single-orbital Hubbard model calculations via the density matrix renormalization group method, we systematically studied the monolayer VOI 2 with a 3d 1 electronic configuration. Our phonon calculations indicate that the orthorhombic Pmm2 FE-II phase is the most likely ground state, involving a ferroelectric (FE) distortion along the a axis and V-V dimerization along the b axis. Specifically, the “pseudo Jahn-Teller” effect caused by the coupling between empty V (d xz/yz and d 3z 2 -r 2 ) and O 2p states is proposed as the mechanism that stabilizes the FE distortion from the paraelectric phase. Moreover, the half-filled metallic d xy band displays a Peierls instability along the b axis, inducing a V-V dimerization. Additionally, we found very short-range antiferromagnetic coupling along the V-V chain due to the formation of nearly decoupled spin singlets in the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

High-precision monitoring of outgassing species in model extreme ultraviolet photoresists with a cavity ring-down spectrometer

Background Extreme ultraviolet (EUV) photoresists play a pivotal role in advancing nanopatterning technologies by balancing image quality and sensitivity. The outgassing behavior of photoresist thin films under EUV and deep ultraviolet (DUV) exposure reveals chemical details relevant to their performance. Aim Here, we focus on utilizing an analytical technique not previously used in photolithography, the tabletop cavity ring-down spectrometer, to investigate outgassing dynamics in EUV photoresists, enabling precise chemical identification and deeper insights into resist processing. Approach The spectrometer’s enhanced laser path length (∼ 20 km) and broadband absorption capabilities in the C–H overtone region allow for sensitive and temporally resolved detection of mixtures of outgassed species. Using a model resist comprising a polymer matrix with a photoacid generator and quencher, we analyzed the influence of time delays between exposure and post-exposure bake (PEB) as well as storage under varying environmental conditions. Results Suppression of isobutylene outgassing and thickness loss was observed with extended delays between exposure and PEB, potentially linked to water absorption and acid deactivation. The technique proved highly effective in distinguishing subtle chemical differences between processing stages. Conclusions Delay times and their environmental conditions, particularly humidity, reduce outgassing and thickness loss of photoresists during PEB, suggesting decreased acid-driven deprotection. This can potentially impact sensitivity, defectivity, and roughness of resist patterns, necessitating precise monitoring and control.

Extreme ultraviolet photoresists↗