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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Elastic scattering of low-energy electrons by 1,4-dioxane

We report calculated cross sections for elastic collisions of low-energy-electrons with 1,4-dioxane. Our calculations employed the Schwinger multichannel method with pseudopotentials and were carried out in the static-exchange and static-exchange plus polarization approximations for energies up to 30 eV. Our results show the presence of three shape resonances belonging to the B{sub u}, A{sub u}, and B{sub g} symmetries and located at 7.0 eV, 8.4 eV, and 9.8 eV, respectively. We also report the presence of a Ramsauer-Townsend minimum located at around 0.05 eV. We compare our calculated cross sections with experimental data and R-matrix and independent atom model along with the additivity rule corrected by using screening coefficients theoretical results for 1,4-dioxane obtained by Palihawadana et al. [J. Chem. Phys. 139, 014308 (2013)]. The agreement between the present and the R-matrix theoretical calculations of Palihawadana et al. is relatively good at energies below 10 eV. Our calculated differential cross sections agree well with the experimental data, showing only some discrepancies at higher energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Robust Matrix Completion State Estimation in Distribution Systems

Due to the insufficient measurements in the distribution system state estimation (DSSE), full observability and redundant measurements are difficult to achieve without using the pseudo measurements. The matrix completion state estimation (MCSE) combines the matrix completion and power system model to estimate voltage by exploring the low-rank characteristics of the matrix. This paper proposes a robust matrix completion state estimation (RMCSE) to estimate the voltage in a distribution system under a low-observability condition. Tradition state estimation weighted least squares (WLS) method requires full observability to calculate the states and needs redundant measurements to proceed a bad data detection. The proposed method improves the robustness of the MCSE to bad data by minimizing the rank of the matrix and measurements residual with different weights. It can estimate the system state in a low-observability system and has robust estimates without the bad data detection process in the face of multiple bad data. The method is numerically evaluated on the IEEE 33-node radial distribution system. The estimation performance and robustness of RMCSE are compared with the WLS with the largest normalized residual bad data identification (WLS-LNR), and the MCSE.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Multiresolution spatiotemporal mechanical model of the heart as a prior to constrain the solution for 4D models of the heart

In several nuclear cardiac imaging applications (SPECT and PET), images are formed by reconstructing tomographic data using an iterative reconstruction algorithm with corrections for physical factors involved in the imaging detection process and with corrections for cardiac and respiratory motion. The physical factors are modeled as coefficients in the matrix of a system of linear equations and include attenuation, scatter, and spatially varying geometric response. The solution to the tomographic problem involves solving the inverse of this system matrix. This requires the design of an iterative reconstruction algorithm with a statistical model that best fits the data acquisition. The most appropriate model is based on a Poisson distribution. Using Bayes Theorem, an iterative reconstruction algorithm is designed to determine the maximum a posteriori estimate of the reconstructed image with constraints that maximizes the Bayesian likelihood function for the Poisson statistical model. The a priori distribution is formulated as the joint entropy (JE) to measure the similarity between the gated cardiac PET image and the cardiac MRI cine image modeled as a FE mechanical model. In conclusion, the developed algorithm shows the potential of using a FE mechanical model of the heart derived from a cardiac MRI cine scan to constrain solutions of gated cardiac PET images.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Modeling of laser interactions with composite materials

In this study, we develop models of laser interactions with composite materials consisting of fibers embedded within a matrix. A ray-trace model is shown to determine the absorptivity, absorption depth, and optical power enhancement within the material, as well as the angular distribution of the reflected light. We also develop a macroscopic model, which provides physical insight and overall results. We show that the parameters in this model can be determined from the ray trace model.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Robust Matrix Completion State Estimation in Distribution Systems: Preprint

Due to the insufficient measurements in the distribution system state estimation (DSSE), full observability and redundant measurements are difficult to achieve without using the pseudo measurements. The matrix completion state estimation (MCSE) combines the matrix completion and power system model to estimate voltage by exploring the low-rank characteristics of the matrix. This paper proposes a robust matrix completion state estimation (RMCSE) to estimate the voltage in a distribution system under a low-observability condition. Tradition state estimation weighted least squares (WLS) method requires full observability to calculate the states and needs redundant measurements to proceed a bad data detection. The proposed method improves the robustness of the MCSE to bad data by minimizing the rank of the matrix and measurements residual with different weights. It can estimate the system state in a low-observability system and has robust estimates without the bad data detection process in the face of multiple bad data. The method is numerically evaluated on the IEEE 33-node radial distribution system. The estimation performance and robustness of RMCSE are compared with the WLS with the largest normalized residual bad data identification (WLS-LNR), and the MCSE.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Harmonic Modeling, Data Generation and Analysis of Power Electronics-Interfaced Residential Loads

The share of electronics-based residential load is expected to rise as devices such as variable frequency drives (VFDs), electric vehicle chargers, and inverter-based distributed energy resources (DERs), e.g., photovoltaic (PV) systems become more common. These loads may introduce significant harmonics into power networks that need to be closely studied in order to perform accurate load modeling and forecasting. However, it can be difficult to obtain harmonic-rich voltage and current data - necessary for identifying accurate load models - for residential electrical loads. Recognizing this need, we identify and model a set of electronics-based end-use loads and DERs in an electromagnetic transients program (EMTP) tool for a residence with a single- phase split-phase supply. Further, a procedure is developed to model harmonic interactions between end-use loads connected to the same non-ideal supply voltage in a residential setting. Finally, a harmonic-rich dataset produced via the proposed procedure is utilized to identify frequency coupling matrix (FCM) based load model. Numerical results demonstrate the accuracy of the model, and explore model identifiability with limited data points.

harmonics, power quality, load modeling↗
DOE OSTI2015

Brief announcement: Hypergraph parititioning for parallel sparse matrix-matrix multiplication

The performance of parallel algorithms for sparse matrix-matrix multiplication is typically determined by the amount of interprocessor communication performed, which in turn depends on the nonzero structure of the input matrices. In this paper, we characterize the communication cost of a sparse matrix-matrix multiplication algorithm in terms of the size of a cut of an associated hypergraph that encodes the computation for a given input nonzero structure. Obtaining an optimal algorithm corresponds to solving a hypergraph partitioning problem. Furthermore, our hypergraph model generalizes several existing models for sparse matrix-vector multiplication, and we can leverage hypergraph partitioners developed for that computation to improve application-specific algorithms for multiplying sparse matrices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Structural conserved moiety splitting of a stoichiometric matrix

Characterising biochemical reaction network structure in mathematical terms enables the inference of functional biochemical consequences from network structure with existing mathematical techniques and spurs the development of new mathematics that exploits the peculiarities of biochemical network structure. The structure of a biochemical network may be specified by reaction stoichiometry, that is, the relative quantities of each molecule produced and consumed in each reaction of the network. A biochemical network may also be specified at a higher level of resolution in terms of the internal structure of each molecule and how molecular structures are transformed by each reaction in a network. The stoichiometry for a set of reactions can be compiled into a stoichiometric matrix N ϵ Z mxn where each row corresponds to a molecule and each column corresponds to a reaction. We demonstrate that a stoichiometric matrix may be split into the sum of m - rank(N) moiety transition matrices, each of which corresponds to a subnetwork accessible to a structurally identifiable conserved moiety. The existence of this moiety matrix splitting is a property that distinguishes a stoichiometric matrix from an arbitrary rectangular matrix.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Spectral fluctuations in the Sachdev-Ye-Kitaev model

We present a detailed quantitative analysis of spectral correlations in the Sachdev-Ye-Kitaev (SYK) model. We find that the deviations from universal Random Matrix Theory (RMT) behavior are due to a small number of long-wavelength fluctuations (of the order of the number of Majorana fermions N ) from one realization of the ensemble to the next one. These modes can be parameterized effectively in terms of Q-Hermite orthogonal polynomials, the main contribution being due to scale fluctuations for which we give a simple analytical estimate. Our numerical results for N = 32 show that only the lowest eight polynomials are needed to eliminate the nonuniversal part of the spectral fluctuations. The covariance matrix of the coefficients of this expansion can be obtained analytically from low-order double-trace moments. We evaluate the covariance matrix of the first six moments and find that it agrees with the numerics. We also analyze the spectral correlations in terms of a nonlinear σ-model, which is derived through a Fierz transformation, and evaluate the one and two-point spectral correlation functions to two-loop order. The wide correlator is given by the sum of the universal RMT result and corrections whose lowest-order term corresponds to scale fluctuations. However, the loop expansion of the σ-model results in an ill-behaved expansion of the resolvent, and it gives universal RMT fluctuations not only for q= 4 or higher even q-body interactions, but also for the q= 2 SYK model albeit with a much smaller Thouless energy while the correct result in this case should have been Poisson statistics. In our numerical studies we analyze the number variance and spectral form factor for N = 32 and q = 4. We show that the quadratic deviation of the number variance for large energies appears as a peak for small times in the spectral form factor. After eliminating the long-wavelength fluctuations, we find quantitative agreement with RMT up to an exponentially large number of level spacings for the number variance or exponentially short times in the case of the spectral form factor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Multiscale Modeling of Silicon Carbide Cladding for Nuclear Applications: Thermal Performance Modeling

The complex multiscale and anisotropic nature of silicon carbide (SiC) ceramic matrix composite (CMC) makes it difficult to accurately model its performance in nuclear applications. The existing models for nuclear grade composite SiC do not account for the microstructural features and how these features can affect the thermal and structural behavior of the cladding and its anisotropic properties. In addition to the microstructural features, the properties of individual constituents of the composites and fiber tow architecture determine the bulk properties. Models for determining the relationship between the individual constituents’ properties and the bulk properties of SiC composites for nuclear applications are absent, although empirical relationships exist in the literature. Here, a hierarchical multiscale modeling approach was presented to address this challenge. This modular approach addressed this difficulty by dividing the various aspects of the composite material into separate models at different length scales, with the evaluated property from the lower-length-scale model serving as an input to the higher-length-scale model. The multiscale model considered the properties of various individual constituents of the composite material (fiber, matrix, and interphase), the porosity in the matrix, the fiber volume fraction, the composite architecture, the tow thickness, etc. By considering inhomogeneous and anisotropic contributions intrinsically, our bottom-up multiscale modeling strategy is naturally physics-informed, bridging constitutive law from micromechanics to meso-mechanics and structural mechanics. The effects that these various physical attributes and thermo-physical properties have on the composite’s bulk thermal properties were easily evaluated and demonstrated through the various analyses presented herein. Since silicon carbide fiber-reinforced SiC CMCs are also promising thermal–structural materials with a broad range of high-end technology applications beyond nuclear applications, we envision that the multiscale modeling method we present here may prove helpful in future efforts to develop and construct reinforced CMCs and other advanced composite nuclear materials, such as MAX phase materials, that can service under harsh environments of ultrahigh temperatures, oxidation, corrosion, and/or irradiation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Objective identification of informative wavelength regions in galaxy spectra

Understanding the diversity in spectra is the key to determining the physical parameters of galaxies. The optical spectra of galaxies are highly convoluted with continuum and lines that are potentially sensitive to different physical parameters. Defining the wavelength regions of interest is therefore an important question. In this work, we identify informative wavelength regions in a single-burst stellar population model using the CUR Matrix Decomposition. Simulating the Lick/IDS spectrograph configuration, we recover the widely used D {sub n}(4000), Hβ, and Hδ {sub A} to be most informative. Simulating the Sloan Digital Sky Survey spectrograph configuration with a wavelength range 3450-8350 Å and a model-limited spectral resolution of 3 Å, the most informative regions are: first region—the 4000 Å break and the Hδ line; second region—the Fe-like indices; third region—the Hβ line; and fourth region—the G band and the Hγ line. A principal component analysis on the first region shows that the first eigenspectrum tells primarily the stellar age, the second eigenspectrum is related to the age-metallicity degeneracy, and the third eigenspectrum shows an anti-correlation between the strengths of the Balmer and the Ca K and H absorptions. The regions can be used to determine the stellar age and metallicity in early-type galaxies that have solar abundance ratios, no dust, and a single-burst star formation history. The region identification method can be applied to any set of spectra of the user's interest, so that we eliminate the need for a common, fixed-resolution index system. We discuss future directions in extending the current analysis to late-type galaxies. ASCII formatted tables of the regional eigenspectra are available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Leptonic Decays of Charged Pseudoscalar Mesons - 2015

We review the physics of purely leptonic decays of $\pi^\pm$, $K^\pm$, $D^{\pm}$, $D_s^\pm$, and $B^\pm$ pseudoscalar mesons. The measured decay rates are related to the product of the relevant weak-interaction-based CKM matrix element of the constituent quarks and a strong interaction parameter related to the overlap of the quark and antiquark wave-functions in the meson, called the decay constant $f_P$. The leptonic decay constants for $\pi^\pm$, $K^\pm$, $D^{\pm}$, $D_s^\pm$, and $B^\pm$ mesons can be obtained with controlled theoretical uncertainties and high precision from {\it ab initio} lattice-QCD simulations. The combination of experimental leptonic decay-rate measurements and theoretical decay-constant calculations enables the determination of several elements of the CKM matrix within the standard model. These determinations are competitive with those obtained from semileptonic decays, and also complementary because they are sensitive to different quark flavor-changing currents. They can also be used to test the unitarity of the first and second rows of the CKM matrix. Conversely, taking the CKM elements predicted by unitarity, one can infer "experimental" values for $f_P$ that can be compared with theory. These provide tests of lattice-QCD methods, provided new-physics contributions to leptonic decays are negligible at the current level of precision. This review is the basis of the article in the Particle Data Group's 2016 edition, updating the versions in Refs. [1-3].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

NLOX, a one-loop provider for Standard Model processes

NLOX is a computer program for calculations in high-energy particle physics. It provides fully renormalized scattering matrix elements in the Standard Model of particle physics, up to one-loop accuracy for all possible coupling-power combinations in the strong and electroweak couplings, and for processes with up to six external particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Optical Design and Validation of an Infrared Transmissive Spectrum Splitting Concentrator Photovoltaic Module

A new modular, hybrid solar power system is designed to generate both electrical and thermal energy by utilizing the full solar spectrum. The key element, an infrared-transparent concentrator photovoltaic (CPV) module, acts as a spectrum splitter, dividing solar radiation into two parts. The ultraviolet and visible light (“in-band”) are converted to electricity with high efficiency in CPV cells, while the infrared light (“out-of-band”) is transmitted directly to a thermal receiver, where thermal power may be converted to electricity by a suitable heat engine or used directly for industrial process heat applications whenever needed. Here, we describe the optical design, modeling, fabrication, and performance validation of this novel spectrum splitting CPV module. A transfer matrix style approach, cumulative transmission model, is built to study the reflection, absorption, and transmission in each layer of the CPV module. To optimize the optical performance, different materials for module superstrate/substrate, encapsulant, cell substrate, and cooling fluids are compared in order to enhance the transmission of out-of-band light through the CPV module by minimizing absorption. Six antireflection coatings along with front and backside electrical contact grids are designed to maximize transmittance of in-band light to the cell and out-of-band light to the thermal receiver. Here, the final design, currently being prototyped, predicts out-of-band light transmission to the thermal receiver of 74.1% (for the passively cooled version) and 65.3% (for the actively cooled version). When epitaxial liftoff technology is applied, the transmission will change to 80.8% (passively cooled) and 71.9% (actively cooled). Experimental prototypes show good agreement with modeled optical performance.

14 SOLAR ENERGY↗
DOE OSTI2019

Ultralight Door Design, Development and Demonstration Project (Final Scientific Report)

The project objectives associated with the Ultralight Door Project were to develop a lightweight driver's side door, 42.5% lighter than a current model year high-volume baseline door, at an incremental cost not to exceed $ \$ $ 5 per pound saved, while maintaining equivalent functionality, safety, corrosion, durability and sound transmission performance of the baseline vehicle. The project approach included the selection of a door architecture, development of a key performance target matrix and incremental cost model. Several unique door structures were developed and evaluated using computer-aided engineering (CAE) and cost modeling tools. The aluminum-intensive door structure was the only design capable of meeting both the performance and cost criteria objectives. The door module, glass, interior trim, latch and electrical subsystems were functionally integrated with the aluminum door structure, enabling additional mass reduction beyond material substitution. The lightweight door design was prototyped, integrated with baseline vehicles and tested to validate equivalent functionality and performance. The project results demonstrate a mass reduction of 15.2 kg (40%) at an incremental cost of $ \$ $2.81 per pound saved, while maintaining functionality and performance of the baseline door. The environmental benefit associated with the Ultralight Door enables a 6g CO 2 eq. per kilometer environmental benefit, for a 4-door vehicle powered by a downsized fuel injected gasoline combustion engine. The project results were presented and published in the 2018 SAE International Conference proceedings and in the Journal of Life Cycle Assessment.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Comparison of Coarse-Grained Approaches in Predicting Polymer Nanocomposite Phase Behavior

Because of the considerable parameter space, efficient theoretical and simulation methods are required to predict the morphology and guide experiments in polymer nanocomposites (PNCs). Unfortunately, theoretical and simulation methods are restricted in their ability to accurately map to experiments based on necessary approximations and numerical limitations. In this study, we provide direct comparisons of two recently developed coarse-grained approaches for modeling polymer nanocomposites (PNCs): polymer nanocomposite field theory (PNC-FT) and dynamic mean-field theory (DMFT). These methods are uniquely suited to efficiently capture mesoscale phase behavior of PNCs in comparison to other theoretical and simulation frameworks. We demonstrate the ability of both methods to capture macrophase separation and describe the thermodynamics of PNCs. We systematically test how the nanoparticle morphology in PNCs is affected by a uniform probability distribution of grafting sites, common in field-based methods, versus random discrete grafting sites on the nanoparticle surface. We also analyze the accuracy of the mean-field approximation in capturing the phase behavior of PNCs. Moreover, the DMFT method introduces the ability to describe nonequilibrium phase behavior while the PNC-FT method is strictly an equilibrium method. With the DMFT method we are able to show the evolution of nonequilibrium states toward their equilibrium state and a qualitative assessment of the dynamics in these systems. These simulations are compared to experiments consisting of polystyrene grafted gold nanorods in a poly(methyl methacrylate) matrix to ensure the model gives results that qualitatively agree with the experiments. This study reveals that nanoparticles in a relatively high matrix molecular weight are trapped in a nonequilibrium state and demonstrates the utility of the DMFT framework in capturing nonequilibrium phase behavior of PNCs. In conclusion, both the PNC-FT and DMFT framework are promising methods to describe the thermodynamic and nonequilibrium phase behavior of PNCs.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Hanle-Zeeman Scattering Matrix for Magnetic Dipole Transitions

The polarization of the light that is scattered by the coronal ions is influenced by the anisotropic illumination from the photosphere and the magnetic field structuring in the solar corona. The properties of the coronal magnetic fields can be well studied by understanding the polarization properties of coronal forbidden emission lines that arise from magnetic dipole ( M 1) transitions in the highly ionized atoms that are present in the corona. We present the classical scattering theory of the forbidden lines for a more general case of arbitrary-strength magnetic fields. We derive the scattering matrix for M 1 transitions using the classical magnetic dipole model of Casini and Lin and applying the scattering matrix approach of Stenflo. We consider a two-level atom model and neglect collisional effects. The scattering matrix so derived is used to study the Stokes profiles formed in coronal conditions in those regions where the radiative excitations dominate collisional excitations. To this end, we take into account the integration over a cone of an unpolarized radiation from the solar disk incident on the scattering atoms. Furthermore, we also integrate along the line of sight to calculate the emerging polarized line profiles. We consider radial and dipole magnetic field configurations and spherically symmetric density distributions. For our studies we adopt the atomic parameters corresponding to the [Fe xiii] 10747 Å coronal forbidden line. We also discuss the nature of the scattering matrix for M 1 transitions and compare it with that for the electric dipole ( E 1) transitions.

79 ASTRONOMY AND ASTROPHYSICS↗