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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Apparatus and method for safety analysis evaluation with data-driven workflow

An apparatus and method for system safety analysis evaluation is provided, the apparatus including processing circuitry configured for generating a calculation matrix for a system, generating a plurality of models based on the calculation matrix, performing a benchmarking or convolution analysis of the plurality of models, identifying a design envelope based on the benchmarking or convolution analysis, deriving uncertainty models from the benchmarking or convolution analysis, deriving an assessment judgment based on the uncertainty models and acceptance criteria, defining one or more limiting scenarios based on the design envelope, and determining a safety margin in at least one figure-of-merit for the system based on the design envelope and the acceptance criteria.

Martin, Robert P.↗
DOE OSTI2023

Apparatus and method for safety analysis evaluation with data-driven workflow

An apparatus and method for system safety analysis evaluation is provided, the apparatus including processing circuitry configured for generating a calculation matrix for a system, generating a plurality of models based on the calculation matrix, performing a benchmarking or convolution analysis of the plurality of models, identifying a design envelope based on the benchmarking or convolution analysis, deriving uncertainty models from the benchmarking or convolution analysis, deriving an assessment judgment based on the uncertainty models and acceptance criteria, defining one or more limiting scenarios based on the design envelope, and determining a safety margin in at least one figure-of-merit for the system based on the design envelope and the acceptance criteria.

Martin, Robert P.↗
DOE OSTI2019

Transverse Young's modulus of carbon/glass hybrid fiber composites

Hybrid fiber composites offer designers a means of tailoring the stress–strain behavior of lightweight materials used in high-performance structures. While the longitudinal stress–strain behavior of unidirectional hybrid fiber composites has been thoroughly evaluated experimentally and analytically, relatively little information is available on the transverse behavior. The objective of the current investigation is to present data on the transverse modulus of elasticity of unidirectional composites with five different ratios of carbon and glass fiber and to compare the data with predictive and fitted models. The transverse modulus increases monotonically with the proportion of glass fiber in the composite. Finite element analysis was used to evaluate different ways to model voids in the matrix and allowed the unknown transverse properties of the carbon fibers to be backed out using experimental data from the all-carbon composite. The finite element results show that the transverse modulus can be accurately modeled if voids are modeled explicitly in the matrix region and if modulus is calculated based on stress applied along the minimum interfiber distance path between adjacent fibers arranged in a rectangular array. The transverse modulus was under-predicted by the iso-stress model and was well predicted by a modified iso-stress model and a modified Halpin–Tsai model.

Materials Science↗
DOE OSTI2025

Predicting Flow in Fracture Networks With Quantum Algorithms

Uncertainty quantification plays a crucial role in the modeling of subsurface flow. For instance, uncertainties in the properties of geologic fracture networks significantly impact flow, requiring numerous simulations to accurately estimate quantities of interest. However, each simulation is computationally expensive because it requires solving a large linear system to capture features that involve both small and large fractures. An example is in percolation, where the interaction of many small fractures (which cumulatively can have a large surface area) with the rock matrix must be modeled precisely. Quantum computing is an emerging tool with the potential to address this issue. Quantum algorithms offer a significant speedup in solving linear systems, achieving efficiencies that are challenging to match with classical approaches. These classical approaches include direct solvers, such as LU decomposition, and iterative methods, notably preconditioned conjugate gradient, commonly used in subsurface modeling to solve large sparse systems. However, applying quantum algorithms to geologic fracture flow requires careful attention to algorithmic and problem-specific constraints to fully realize this quantum advantage. In this work we describe a quantum algorithm for generalized Monte Carlo applications with a quadratic speedup over the classical approaches which can be combined with the quantum speedup, currently under investigation, for solving quantum linear systems for subsurface flow. We show that for quantum algorithms the computational cost of estimating a quantity of interest for a statistical ensemble of networks is roughly the same as that of a single realization, essentially implying that one can get uncertainty quantification for free.

58 GEOSCIENCES↗
DOE OSTI2015

Mueller-matrix mapping of biological tissues in differential diagnosis of optical anisotropy mechanisms of protein networks

We report a model of Mueller-matrix description of optical anisotropy of protein networks in biological tissues with allowance for the linear birefringence and dichroism. The model is used to construct the reconstruction algorithms of coordinate distributions of phase shifts and the linear dichroism coefficient. In the statistical analysis of such distributions, we have found the objective criteria of differentiation between benign and malignant tissues of the female reproductive system. From the standpoint of evidence-based medicine, we have determined the operating characteristics (sensitivity, specificity and accuracy) of the Mueller-matrix reconstruction method of optical anisotropy parameters and demonstrated its effectiveness in the differentiation of benign and malignant tumours. (laser applications and other topics in quantum electronics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Coupled-cavity terahertz quantum cascade lasers for single mode operation

We demonstrate the operation of coupled-cavity terahertz frequency quantum-cascade lasers composed of two sub-cavities separated by an air gap realized by optical lithography and dry etching. This geometry allows stable, single mode operation with typical side mode suppression ratios in the 30–40 dB range. We employ a transfer matrix method to model the mode selection mechanism. The obtained results are in good agreement with the measurements and allow prediction of the operating frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Strong thermal leptogenesis and the absolute neutrino mass scale

We show that successful strong thermal leptogenesis, where the final asymmetry is independent of the initial conditions and in particular a large pre-existing asymmetry is efficiently washed-out, favours values of the lightest neutrino mass m{sub 1}∼>10 meV for normal ordering (NO) and m{sub 1}∼>3 meV for inverted ordering (IO) for models with orthogonal matrix entries respecting |Ω{sub ij}{sup 2}|∼<2. We show analytically why lower values of m{sub 1} require a higher level of fine tuning in the seesaw formula and/or in the flavoured decay parameters (in the electronic for NO, in the muonic for IO). We also show how this constraint exists thanks to the measured values of the neutrino mixing angles and could be tightened by a future determination of the Dirac phase. Our analysis also allows us to place a more stringent constraint for a specific model or class of models, such as SO(10)-inspired models, and shows that some models cannot realise strong thermal leptogenesis for any value of m{sub 1}. A scatter plot analysis fully supports the analytical results. We also briefly discuss the interplay with absolute neutrino mass scale experiments concluding that they will be able in the coming years to either corner strong thermal leptogenesis or find positive signals pointing to a non-vanishing m{sub 1}. Since the constraint is much stronger for NO than for IO, it is very important that new data from planned neutrino oscillation experiments will be able to solve the ambiguity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The Sylvester Graphical Lasso (SyGlasso)

This paper introduces the Sylvester graphical lasso (SyGlasso) that captures multiway dependencies present in tensor-valued data. The model is based on the Sylvester equation that defines a generative model. The proposed model complements the tensor graphical lasso that imposes a Kronecker sum model for the inverse covariance matrix by providing an alternative Kronecker sum model that is generative and interpretable. A nodewise regression approach is adopted for estimating the conditional independence relationships among variables. The statistical convergence of the method is established, and empirical studies are provided to demonstrate the recovery of meaningful conditional dependency graphs. Here we apply the SyGlasso to an electroencephalography (EEG) study to compare the brain connectivity of alcoholic and nonalcoholic subjects. Here, we demonstrate that our model can simultaneously estimate both the brain connectivity and its temporal dependencies.

Wang, Yu↗
DOE OSTI2015

Density matrix embedding in an antisymmetrized geminal power bath

Density matrix embedding theory (DMET) has emerged as a powerful tool for performing wave function-in-wave function embedding for strongly correlated systems. In traditional DMET, an accurate calculation is performed on a small impurity embedded in a mean field bath. Here, we extend the original DMET equations to account for correlation in the bath via an antisymmetrized geminal power (AGP) wave function. The resulting formalism has a number of advantages. First, it allows one to properly treat the weak correlation limit of independent pairs, which DMET is unable to do with a mean-field bath. Second, it associates a size extensive correlation energy with a given density matrix (for the models tested), which AGP by itself is incapable of providing. Third, it provides a reasonable description of charge redistribution in strongly correlated but non-periodic systems. Thus, AGP-DMET appears to be a good starting point for describing electron correlation in molecules, which are aperiodic and possess both strong and weak electron correlation.

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

Phase identification using co‐association matrix ensemble clustering

Calibrating distribution system models to aid in the accuracy of simulations such as hosting capacity analysis is increasingly important in the pursuit of the goal of integrating more distributed energy resources. The recent availability of smart meter data is enabling the use of machine learning tools to automatically achieve model calibration tasks. This research focuses on applying machine learning to the phase identification task, using a co‐association matrix‐based, ensemble spectral clustering approach. The proposed method leverages voltage time series from smart meters and does not require existing or accurate phase labels. This work demonstrates the success of the proposed method on both synthetic and real data, surpassing the accuracy of other phase identification research.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2013

Nuclear matrix elements for 0νβ{sup −}β{sup −} decays: Comparative analysis of the QRPA, shell model and IBM predictions

In this work we report on general properties of the nuclear matrix elements involved in the neutrinoless double β{sup −} decays (0νβ{sup −}β{sup −} decays) of several nuclei. A summary of the values of the NMEs calculated along the years by the Jyväskylä-La Plata collaboration is presented. These NMEs, calculated in the framework of the quasiparticle random phase approximation (QRPA), are compared with those of the other available calculations, like the Shell Model (ISM) and the interacting boson model (IBA-2)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Thermal aging modeling and validation on the Mo-containing Fe-Cr-Ni alloys

Thermodynamics of intermetallic phases in Fe-rich Fe-Cr-Ni-Mo alloys is critical knowledge to understand thermal aging effect on the phase stability of Mo-containing austenitic steels, which subsequently facilitates alloy design/improvement and degradation mitigation of these materials for reactor applications. Among the intermetallic phases, Chi (χ), Laves, and Sigma (σ) are often of concern because of their tendency to cause embrittlement of the materials. The focus of this study is thermal stability of the Chi and Laves phases as they were less studied compared to the Sigma phase. Coupled with thermodynamic modeling, thermal stability of intermetallic phases in Mo containing Fe-Cr-Ni alloys was investigated at 1000, 850 and 700°C for different annealing times. The morphologies, compositions and crystal structures of the precipitates of the intermetallic phases were carefully examined by scanning electron microscopy, electron probe microanalysis, X-ray diffraction, and transmission electron microscopy. Three key findings resulted from this study. First, the Chi phase is stable at high temperature, and with decreasing temperature it transforms into the Laves phase that is stable at low temperature. Secondly, Cr, Mo, Ni are soluble in both the Chi and Laves phases, with the solubility of Mo playing a major role in the relative stability of the intermetallic phases. Thirdly, in situ transformation from Chi phase to Laves phase was directly observed, which increased the local strain field, generated dislocations in the intermetallic phases, and altered the precipitate phase orientation relationship with the austenitic matrix. The thermodynamic models that were developed and validated were then applied to evaluating the effect of Mo on the thermal stability of intermetallic phases in type 316 and NF709 stainless steels.

36 MATERIALS SCIENCE↗
DOE OSTI2014

WE-E-18A-08: Towards a Next-Generation Electronic Portal Device for Simultaneous Imaging and Dose Verification in Radiotherapy

Purpose: This work forms part of an ongoing study to develop a next-generation electronic portal imaging device (EPID) for simultaneous imaging and dose verification in radiotherapy. Monte Carlo (MC) simulations were used to characterize the imaging performance of a novel EPID that has previously been demonstrated to exhibit a water-equivalent response. The EPID ' s response was quantified in several configurations and model parameters were empirically validated against experimental measurements. Methods: A MC model of a novel a-Si EPID incorporating an array of plastic scintillating fibers was developed. Square BCF-99-06A scintillator fibers with PMMA cladding (Saint-Gobain Crystals) were modelled in a matrix with total area measuring 150×150 mm{sup 2}. The standard electromagnetic and optical physics Geant4 classes were used to simulate radiation transport from an angled slit source (6 MV energy spectrum) through the EPID and optical photons reaching the photodiodes were scored. The prototype's modulation transfer function (MTF) was simulated and validated against experimental measurements. Several optical transport parameters, fiber lengths and thicknesses of an air gap between the scintillator and photodiodes were investigated to quantify their effects on the prototype's detection efficiency, sensitivity and MTF. Results: Simulated EPID response was more sensitive to variations in geometry than in the optical parameters studied. The MTF was particularly sensitive to the introduction of a 0.5–1.0 mm air gap between the scintillator and photodiodes, which lowered the MTF relative to that simulated without the gap. As expected, increasing the fiber length increased the detector efficiency and sensitivity while decreasing the MTF. Conclusion: A model of a novel water-equivalent EPID has been developed and benchmarked against measurements using a physical prototype. We have demonstrated the feasibility of this new device and are continuing to optimize the design to achieve an imaging response that warrants the development of a next-generation prototype.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Quantum dynamics simulations of the 2D spectroscopy for exciton polaritons

We develop an accurate and numerically efficient non-adiabatic path-integral approach to simulate the non-linear spectroscopy of exciton–polariton systems. This approach is based on the partial linearized density matrix approach to model the exciton dynamics with explicit propagation of the phonon bath environment, combined with a stochastic Lindblad dynamics approach to model the cavity loss dynamics. Through simulating both linear and polariton two-dimensional electronic spectra, we systematically investigate how light–matter coupling strength and cavity loss rate influence the optical response signal. Our results confirm the polaron decoupling effect, which is the reduced exciton–phonon coupling among polariton states due to the strong light–matter interactions. We further demonstrate that the polariton coherence time can be significantly prolonged compared to the electronic coherence outside the cavity.

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

Kaon mixing beyond the standard model with physical masses

We present nonperturbative results for beyond the standard model kaon mixing matrix elements in the isospin symmetric limit ( m u = m d ) of QCD, including a complete estimate of all dominant sources of systematic error. Our results are obtained from numerical simulations of lattice QCD with N f = 2 + 1 flavors of dynamical domain wall fermions. For the first time, these quantities are simulated directly at the physical pion mass m π ∼ 139 MeV for two different lattice spacings. We include data at three lattice spacings in the range a = 0.11 – 0.07 fm and with pion masses ranging from the physical value up to 450 MeV. Compared to our earlier work, we have added both direct calculations at physical quark masses and a third lattice spacing making the removal of discretization effects significantly more precise and eliminating the need for any significant mass extrapolation beyond the range of simulated data. We renormalize the lattice operators nonperturbatively using RI-SMOM off-shell schemes. These schemes eliminate the need to model and subtract nonperturbative pion poles that arises in the RI-MOM scheme and, since the calculations are performed with domain wall fermions, the unphysical mixing between chirality sectors is suppressed. Our results for the bag parameters in the MS ¯ scheme at 3 GeV are B K ≡ B 1 = 0.5240 ( 17 ) ( 54 ) , B 2 = 0.4794 ( 25 ) ( 35 ) , B 3 = 0.746 ( 13 ) ( 17 ) , B 4 = 0.897 ( 02 ) ( 10 ) and B 5 = 0.6882 ( 78 ) ( 94 ) , where the first error is from lattice uncertainties and the second is the uncertainty due to the perturbative matching to MS ¯ . Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

A revised partiality model and post-refinement algorithm for X-ray free-electron laser data

An updated partiality model and post-refinement algorithm for XFEL snapshot diffraction data is presented and confirmed by observing anomalous density for S atoms at an X-ray wavelength of 1.3 Å. Research towards using X-ray free-electron laser (XFEL) data to solve structures using experimental phasing methods such as sulfur single-wavelength anomalous dispersion (SAD) has been hampered by shortcomings in the diffraction models for X-ray diffraction from FELs. Owing to errors in the orientation matrix and overly simple partiality models, researchers have required large numbers of images to converge to reliable estimates for the structure-factor amplitudes, which may not be feasible for all biological systems. Here, data for cytoplasmic polyhedrosis virus type 17 (CPV17) collected at 1.3 Å wavelength at the Linac Coherent Light Source (LCLS) are revisited. A previously published definition of a partiality model for reflections illuminated by self-amplified spontaneous emission (SASE) pulses is built upon, which defines a fraction between 0 and 1 based on the intersection of a reflection with a spread of Ewald spheres modelled by a super-Gaussian wavelength distribution in the X-ray beam. A method of post-refinement to refine the parameters of this model is suggested. This has generated a merged data set with an overall discrepancy (by calculating the R{sub split} value) of 3.15% to 1.46 Å resolution from a 7225-image data set. The atomic numbers of C, N and O atoms in the structure are distinguishable in the electron-density map. There are 13 S atoms within the 237 residues of CPV17, excluding the initial disordered methionine. These only possess 0.42 anomalous scattering electrons each at 1.3 Å wavelength, but the 12 that have single predominant positions are easily detectable in the anomalous difference Fourier map. It is hoped that these improvements will lead towards XFEL experimental phase determination and structure determination by sulfur SAD and will generally increase the utility of the method for difficult cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Neutrinoless double-β decay of 124 Sn, 130 Te, and 136 Xe in the Hamiltonian-based generator-coordinate method

We present a generator-coordinate method for realistic shell-model Hamiltonians that closely approximates the full shell model calculations of the matrix elements for the neutrinoless double-β decay of 124 Sn, 130 Te, and 136 Xe. Here, we treat not only quadrupole deformations but also the proton-neutron pairing amplitudes as generator coordinates. We validate this method by calculating and comparing spectroscopic quantities with the exact shell model results and experimental data. Our Hamiltonian-based generator-coordinate method produces 0νββ matrix elements much closer to the shell model ones, compared to the existing energy-density-functional-based generator-coordinate approaches. The remaining overestimation of 0νββ nuclear matrix element suggests that additional correlations may be needed to be taken into account for 124 Sn, 130 Te, and 136 Xe when calculating with the Hamiltonian-based generator-coordinate method. The validation of this method may open the possibility of calculating 0νββ matrix element of 150 Nd in a large shell-model space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗