Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix factorization”

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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Enhancing The Thermal Resistivity of Rigid Polyisocyanurate Foam Insulation

The development of rigid polyurethane foam insulation has garnered considerable attention because of its promising applications in the buildings and construction industry. Its low thermal conductivity makes it an attractive choice for improving energy performance in buildings. Current Rigid Polyurethane foams have thermal resistivity (R-value/in.) 5.5 to 6.5 h.ft2F/BTU/in.· that could be further improved by diminishing the heat transfer through the foam matrix. Nevertheless, minimizing both conduction (through gas and solid) and radiation simultaneously in porous solids is a significant challenge due to the trade-off between these two mechanisms. This study focuses on improving the R-value of the insulation foams via several strategies: such as type and the amount of the surfactants, blowing agent content and precooling and premixing polyol mixture. These methods optimized thermal properties of the PIR foams, achieving R/in. as high as 8.3. This excellent R/in. is anticipated to be a critical factor in significantly advancing the thermal insulation performance of rigid polyurethane cellular foams, thereby enhancing their efficacy in energy-efficient building applications.

Wanasinghe Mudiyanselage Pahala Gedara, Shiwanka V↗
DOE OSTI2024

Improved Evaluation of Large Network Matrices for Linear Power Flow Within Optimization Problems

This work presents methods for evaluating the Power Transfer Distribution Factor (PTDF) and Line Outage Distribution Factor (LODF) matrices by employing sparse linear algebra for large-scale computing applications. These matrices play a critical role in many power system applications, such as the Unit Commitment Problem (UC), pre- and post-contingency power flow analysis, and transmission expansion. These matrices are typically dense, which means they require a significant amount of time and memory to be computed for large networks. However, by analyzing the structure of the matrices and their computation method, it is possible to use reduced memory methods based on sparse matrix operations. This paper shows that sparse linear algebra algorithms are faster and require less memory and time than traditional dense approaches. Additionally, we explore the effect of matrix sparsification by eliminating trailing digits on power flow calculations.

large scale↗
DOE OSTI2020

Transparent window film with embedded nano-shades for thermoregulation

Every year, more than 20% of the energy consumption in the United States and more than 10% of the global consumption is used towards HVAC (Heating, Ventilating, Air-Conditioning) systems in buildings. Although continued efforts in advancing renewable energy, efficient appliances, and smart building systems are desired, the most critical factors that cause extensive indoor energy consumption are thermal leakage and thermal waste generation. We demonstrate a smart window film that allows natural lighting with a clear view while blocking oblique incident sunlight to make the building interior cooler without using electricity or generating heat. Building upon the concept of, often overlooked, traditional window shades and privacy window films, we emphasize and elevate the application and impact of using magnetically arranged nanoscale material for energy conservation purpose in buildings. By embedding vertically aligned Ni flakes (VANF) in a polymer matrix film, nanoscale mirror array serve as shades to block the high-angled incident sunlight, blocking approximately 73% of total solar irradiance and bringing down the indoor temperature by 1.92 °C during the daytime in a model house while providing nearly 90% transparency. Global application was simulated by applying the respective solar angle and intensity at selected latitudes, indicating their potential. This passively energy-conserving smart film is a low cost, maintenance free, and simple product for users and easy-to-manufacture product for the industry.

36 MATERIALS SCIENCE↗
NASA NTRS1983

Numerical studies of major and minor stratospheric warmings caused by orographic forcing

A primitive equation spectral model using spherical harmonics is formulated to study dynamic interactions between the troposphere and stratosphere in association with sudden stratospheric warmings. Using sigma coordinates for five tropospheric layers and log-pressure coordinates for 26 stratospheric and mesospheric layers, separate model equations for each system are combined to form single matrix governing equations. The gradual introduction of large scale topography to balanced initial states representative of observed mean winter conditions in the Northern Hemisphere is used for the generation of planetary waves during 40-day time integrations. Results of these integrations indicate that stratospheric warmings can be simulated by this orographic forcing and that mean momentum flux divergence due to zonal mean motion appears to be an essential mechanism of these simulated sudden warmings. It was found that the strength of the polar night jet can be a determining factor whether a warming becomes 'major' or 'minor'.

Koermer, J. P.↗
DOE OSTI2023

Measurements and computational analysis of the natural decay of 176 Lu

Background: Mainly because of its long half-life and despite its scientific relevance, spectroscopic measurements of 176 Lu forbidden β decays are very limited and lack formulation of shape factors. A direct precise measurement of its Q value is also presently unreported. In addition, the description of forbidden decays provides interesting challenges for nuclear theory. The comparison of precise experimental results with theoretical calculations for these decays can help to test underlying models and can aid the interpretation of data from other experiments. Purpose: Perform the first precision measurements of 176 Lu β-decay spectra and attempt the observation of its electron capture decays, as well as perform the first precision direct measurement of the 176 Lu β-decay Q value. Compare the shape of the precisely determined experimental β spectra to theoretical calculations, and compare the end point energy to that obtained from an independent Q value measurement. Method: Here, the 176 Lu β-decay spectra measurements and the search for electron capture decays were performed with an experimental setup that employed lutetium-containing scintillator crystals and a NaI(Tl) spectrometer for coincidence counting. The β decay Q value was determined via high-precision Penning trap mass spectrometry (PTMS) with the LEBIT facility at the National Superconducting Cyclotron Laboratory. The β-spectrum calculations were performed within the Fermi theory formalism with nuclear structure effects calculated using a shell model approach. Results: Both β transitions of 176 Lu were experimentally observed and corresponding shape factors formulated in their entire energy ranges. The search for electron capture decay branches led to an experimental upper limit of 6.3×10 –6 relative to its β decays. The 176 Lu β-decay and electron capture Q values were measured using PTMS to be 1193.0(6) and 108.9(8) keV, respectively. This enabled precise β end point energies of 596.2(6) and 195.3(6) keV to be determined for the primary and secondary β decays, respectively. The conserved vector current hypothesis was applied to calculate the relativistic vector matrix elements. The β-spectrum shape was shown to significantly depend on the Coulomb displacement energy and on the value of the axial vector coupling constant gA, which was extracted according to different assumptions. Conclusion: The implemented self-scintillation method has provided unmatched observations of 176 Lu, independently validated by the first direct measurements of its β-decay Q value by Penning trap mass spectrometry. Theoretical study of the main β transition led to the extraction of very different effective gA and log10f values, showing that a high-precision description of this transition would require a realistic nuclear structure with nucleus deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Lithium-Ion Battery Life Model with Electrode Cracking and Early-Life Break-in Processes

This paper develops a physically justified reduced-order capacity fade model from accelerated calendar- and cycle-aging data for 32 lithium-ion (Li-ion) graphite/nickel-manganese-cobalt (NMC) cells. The large data set reveals temperature-, charge C-rate-, depth-of-discharge-, and state of charge (SOC)-dependent degradation patterns that would be unobserved in a smaller test matrix. Model structure is informed by incremental capacity analysis that shows loss of lithium inventory and cathode-material loss as the dominant capacity fade mechanisms. The model includes terms attributable to solid-electrolyte interface (SEI) growth, electrode cracking, cycling-driven acceleration of SEI growth, and "break-in" mechanisms that slightly decrease or increase available Li inventory early in life. The study explores what mathematical couplings of these mechanisms best describe calendar aging, cycle aging, and mixed calendar/cycle aging. Various approaches are discussed for extracting relevant stress factors from complex cycling profiles to predict lifetime during real-world battery loads using models trained on constant-current laboratory test results. The complexity of the present human-driven model identification process motivates future work in machine learning to more widely search and statistically discern the optimal model that correctly extrapolates capacity fade based on physical knowledge.

25 ENERGY STORAGE↗
DOE OSTI2024

Luttinger liquid phase in the Aubry-André Hubbard chain

We study the interplay between an on-site Hubbard repulsion and quasiperiodic potential in one-dimensional fermion chains using the density matrix renormalization group. We find that, at half-filling, the quasiperiodic potential can destroy the Mott gap, leading to a metallic Luttinger liquid phase between the gapped Mott insulator at strong repulsion and localized gapless Aubry- André insulator at strong quasiperiodic potential. Away from half-filing, the metallic phase of the interacting model persists to larger critical strengths of the potential than in the non-interacting case, suggesting interaction-stabilized delocalization at finite doping. Here, we characterize the Luttinger liquid through its charge and spin correlations, structure factors, and entanglement entropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

At-source Recovery of Rare Earth Elements from Coal Mine Drainage

Conventional acid mine drainage (AMD) treatment consists of pH adjustment, oxidation, and separation of solid precipitates from discharge quality water. Alkaline addition increases the pH to the range 7 to 9 to precipitate the regulated metals Al and Mn ions as hydroxides. Mechanical or chemical oxidation then converts the reduced iron species Fe 2+ to Fe 3+ which precipitates at a low pH value as ferric hydroxide (Fe(OH) 3 ). Co-precipitation of Rare Earth Elements (REEs) with hydroxides creates a gangue-rich matrix from which REEs must be separated. Fe(OH)3 is generally the most abundant gangue-forming metal. Current research by the project team has shown raw, untreated AMD has an average total rare earth element content (TREE) of about 287 μg/L (0.287 ppm). Operators of AMD-producing facilities are obliged to treat it to neutralize acidity and remove regulated metals (e.g., iron (Fe), aluminum (Al), and manganese (Mn)) prior to discharge. Ongoing research indicates conventional AMD treatment concentrates TREE in resulting precipitates (AMD sludge) by an average factor of 2,635x to an average dry weight, elemental concentration of 708 g/t (ppm). While the economics of recovering REEs from AMD sludge are appealing, much of the cost of REE recovery would be incurred during separation of REEs from the gangue metals: Fe, Al, and Mn. Significant gains can be made by precipitating REEs when most of these gangue metals are still in solution. Significant improvements in REE extraction efficiency can be obtained through separation of REEs from the aqueous phase AMD, upstream of conventional AMD treatment by: 1) creating an enriched REE feedstock, 2) producing a more consistent feedstock, 3) reducing transportation costs to an REE refinery, 4) reducing acid consumption in the acid leaching step, and 5) reducing the volume of waste produced at the ALSX plant. This project is exploring two novel approaches for extracting REEs upstream of the conventional AMD treatment plants to create a purified REE feedstock while leaving the bulk of the Fe, Al, and Mn in solution for subsequent treatment. While the average TREE concentration of raw AMD is low and individual sites range between 10 and 2,200 μg/L, successful at-source REE separation will generate a superior feedstock to a conventional ALSX process and significantly improve the economics and environmental performance of REE extraction from AMD. AMD can be classified into two types: type A is net acid and iron is partly oxidized while type B is net alkaline, and iron occurs in the reduced, ferrous state. In Case A, the pH is raised slightly, which will precipitate Fe 3+ and Al 3+ but not REEs. The other important type of AMD is net alkaline water found in flooded, anoxic deep mines. These generally have a pH > 6.0 and are net alkaline. REE concentrations are also much lower than those found in acidic AMD while the volumes and flux are much higher. Most of these are pumped discharges with almost no Al. Case B AMD is strongly reduced and, like Case A water, Fe and Mn ions are similarly reduced and therefore soluble at pH <9.0. REEs and an array of transition ions will be the only trivalent cations in this type of AMD. In Case B, both upward pH adjustment under reducing conditions and application of an electrochemically stimulated supported liquid membrane strategy to separate REEs from ferrous ion will be explored.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

LLNL measurements of thermally irradiated HEU sample and saltwater-matrixed HEU sample

Two scintillation vials, kindly provided by the PNNL team, were received in the LLNL radiochemistry building on June 10, 2024; one containing approximately 2 mL dissolved thermally-irradiated HEU, and the other containing dried salt from irradiated seawater. A quick, semiquantitative screening count of the salt showed that the activation product activities were very low, with the highest activity being ~1 Bq 24 Na. The dissolved HEU solution was transferred to a Prindle vial (LLNL standard counting geometry) and weighed. Approximately 2 mL 3 M HNO 3 was added so that the solution completely covered the bottom of the vial and Gamanal software would be able to accurately generate an efficiency curve. This sample, called “PNNLFP24” is the irradiated HEU solution received, gravimetrically diluted by a factor of 1.7356 ± 0.0003, and was sent for quick gamma counting prior to further modification. Meanwhile the salt was quantitatively transferred to a 250 mL polyethylene bottle and dissolved in 160 mL 3 M HNO 3 . After counting, 50 µL PNNLFP24 was aliquoted for Resonance Ionization Mass Spectrometry (RIMS) analysis, and then 10.0 mL of the salt solution (6.03% of the total salt, gravimetrically) was added to the remaining PNNLFP24 solution to make “PNNLFPSW24”. After weighing, two 100 µL aliquots were removed for RIMS and DES, and 3 mL was aliquoted for microfluidic chemistry. The remaining “PNNLFPSW24” solution was weighed and proceeded for singles and coincidence counting on the MCBOS system.

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

Simple and Robust Carrier-Based PWM Technique for Single-Stage Three-phase Rectifier Indirect Matrix Converter

Indirect matrix converters (IMC) can offer high efficiency, high power density, and bidirectional power flow capability with a single-stage three-to single-phase AC-AC power conversion integrating a conventional three-phase active rectifier with a primary side full-bridge inverter of a dual active bridge converter. However, the IMCs generally require space vector modulation techniques, leading to complex sectors and dwell times calculations. A practical implementation of the modulation methods with a commonly used digital signal processor becomes even more challenging when considering power factor correction, harmonics reduction, and deadtime between sector transitions. This paper proposes a simple and robust carrier-based pulse width modulation (PWM) technique determined through a voltage-source inverter (VSI) to current-source inverter (CSI) to IMC conversion. This translation enables the adoption of all conventional VSI control methods with active and reactive power control in the IMCs through simple carrier and reference waveform comparison. The VSI vectors are first one-on-one matched with CSI vectors and translated to IMC modulation signals without calculating sectors and dwell times. Nine different carrier waveforms are proposed, and the converter performance is analyzed. Grid current THD and efficiency are estimated for each modulation method. The symmetric A and B methods had the lowest THD of the grid current, 2.80%, and the efficiency remained at approximately 98.3%, making it the most suitable carrier-based method proposed in this paper. The space vector PWM (SVPWM) modulation method was implemented to the DSP and generated the expected twelve gate signals.

Benson, Mikayla↗
DOE OSTI2022

Emergence of many-body quantum chaos via spontaneous breaking of unitarity

It is suggested that many-body quantum chaos appears as the spontaneous symmetry breaking of unitarity in interacting quantum many-body systems. It has been shown that many-body level statistics, probed by the spectral form factor (SFF) defined as K(η,t)=$\langle$|Tr exp(-ηH+itH)| 2 $\rangle$, is dominated by a diffuson-type mode in a field theory analysis. The key finding of this Letter is that the “unitary” η = 0 case is different from the η → 0 ± limit, with the latter leading to a finite mass of these modes due to interactions. This mass suppresses a rapid exponential ramp in the SFF, which is responsible for the fast emergence of Poisson statistics in the noninteracting case, and gives rise to a nontrivial random matrix structure of many-body levels. The interaction-induced mass in the SFF shares similarities with the dephasing rate in the theory of weak localization and the Lyapunov exponent of the out-of-time-ordered correlators.

36 MATERIALS SCIENCE↗
NASA NTRS1974

The two-dimensional tunnel structures of K3Sb5O14 and K2Sb4O11

The structures of K3Sb5O14 and K2Sb4O11 have been solved by the single-crystal X-ray direct method and the heavy-atom method, respectively. The structure of K3Sb5O14 is orthorhombic, with space group Pbam and cell parameters a = 24.247 (4), b = 7.157 (2), c = 7.334 (2) A, Z = 4. The structure of K2Sb4O11 is monoclinic, with space group C2/m and cell parameters a = 19.473 (4), b = 7.542 (1), c = 7.198 (1) A, beta = 94.82 (2) deg, Z = 4. A full-matrix least-squares refinement gave R = 0.072 and R = 0.067, respectively. In both structures, oxygen atoms form an octahedron around each Sb atom and an irregular polyhedron around each K atom. By sharing corners and edges, the octahedra form a skeleton network having intersecting b-axis and c-axis tunnels. The K(+) ions, which have more than ten oxygen near neighbors, are located in these tunnels. Evidence for K(+)-ion transport within and between tunnels comes from ion exchange of the alkali ions in molten salts and anisotropic temperature factors that are anomalously large in the direction of the tunnels.

Hong, H. Y.-P.↗
DOE OSTI2020

Theory for transport in magnetized plasmas

We present a two fluid model for anomalous thermal, particle and momentum transport which has been derived from kinetic theory using a small parameter ε~10-2. This small parameter can come from either the ratio of wave frequency to cyclotron frequency or the ratio of potential to thermal energies. This fluid model is strongly nonlinear and has an exact closure within the accuracy ε. It can be used to simulate the full radius of a tokamak including internal and edge transport barriers at the same time without introducing any local assumptions for physics model or grid size. It involves a full transport matrix, thus giving possibility for both density and heat pinches. The model can be derived from the Vlasov equation, adding also close collisions and includes particle pinches of the type found in the magnetosphere and in the Massachusetts Institute of Technology elevated dipole machine. Pinches are also in agreement for the heat pinch experiments on Doublet III D including simultaneous particle and heat pinches. For tokamaks it gives the right scaling of confinement time versus heating power. Finally, we have also recovered the Turbulent Equipartition result that the density scales as inverse safety factor q in typical simulations. Thus this model is both fundamental and of high practical use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Vibration localization in mono- and bi-coupled bladed disks - A transfer matrix approach

A transfer matrix approach to the analysis of the dynamics of mistuned bladed disks is presented. The study focuses on mono-coupled systems, in which each blade is coupled to its two neighboring blades, and bi-coupled systems, where each blade is coupled to its four nearest neighbors. Transfer matrices yield the free dynamics, both the characteristic free wave and the normal modes - in closed form for the tuned assemblies. Mistuned assemblies are represented by random transfer matrices and an examination of the effect of mistuning on harmonic wave propagation yields the localization factor - the average rate of spatial wave amplitude decay per blade - in the mono-coupled assembly. Based on a comparison of the wave propagation characteristics of the mono- and bi-coupled assemblies, important conclusions are drawn about the effect of the additional coupling coordinate on the sensitivity to mistuning and the strength of mode localization predicted by a mono-coupled analysis.

Ottarsson, Gisli↗
NASA NTRS1997

The Path Resistance Method for Bounding the Smallest Nontrivial Eigenvalue of a Laplacian

We introduce the path resistance method for lower bounds on the smallest nontrivial eigenvalue of the Laplacian matrix of a graph. The method is based on viewing the graph in terms of electrical circuits; it uses clique embeddings to produce lower bounds on lambda(sub 2) and star embeddings to produce lower bounds on the smallest Rayleigh quotient when there is a zero Dirichlet boundary condition. The method assigns priorities to the paths in the embedding; we show that, for an unweighted tree T, using uniform priorities for a clique embedding produces a lower bound on lambda(sub 2) that is off by at most an 0(log diameter(T)) factor. We show that the best bounds this method can produce for clique embeddings are the same as for a related method that uses clique embeddings and edge lengths to produce bounds.

Guattery, Stephen↗
NASA NTRS1987

A simplified computer solution for the flexibility matrix of contacting teeth for spiral bevel gears

A computer code, FLEXM, was developed to calculate the flexibility matrices of contacting teeth for spiral bevel gears using a simplified analysis based on the elementary beam theory for the deformation of gear and shaft. The simplified theory requires a computer time at least one order of magnitude less than that needed for the complete finite element method analysis reported earlier by H. Chao, and it is much easier to apply for different gear and shaft geometries. Results were obtained for a set of spiral bevel gears. The teeth deflections due to torsion, bending moment, shearing strain and axial force were found to be in the order 10(-5), 10(-6), 10(-7), and 10(-8) respectively. Thus, the torsional deformation was the most predominant factor. In the analysis of dynamic load, response frequencies were found to be larger when the mass or moment of inertia was smaller or the stiffness was larger. The change in damping coefficient had little influence on the resonance frequency, but has a marked influence on the dynamic load at the resonant frequencies.

Hsu, C. Y.↗
NASA NTRS2000

A Matrix Approach to Software Process Definition

The Software Engineering Laboratory (SEL) is currently engaged in a Methodology and Metrics program for the Information Systems Center (ISC) at Goddard Space Flight Center (GSFC). This paper addresses the Methodology portion of the program. The purpose of the Methodology effort is to assist a software team lead in selecting and tailoring a software development or maintenance process for a specific GSFC project. It is intended that this process will also be compliant with both ISO 9001 and the Software Engineering Institute's Capability Maturity Model (CMM). Under the Methodology program, we have defined four standard ISO-compliant software processes for the ISC, and three tailoring criteria that team leads can use to categorize their projects. The team lead would select a process and appropriate tailoring factors, from which a software process tailored to the specific project could be generated. Our objective in the Methodology program is to present software process information in a structured fashion, to make it easy for a team lead to characterize the type of software engineering to be performed, and to apply tailoring parameters to search for an appropriate software process description. This will enable the team lead to follow a proven, effective software process and also satisfy NASA's requirement for compliance with ISO 9001 and the anticipated requirement for CMM assessment. This work is also intended to support the deployment of sound software processes across the ISC.

Schultz, David↗