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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

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

This work discusses 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 are critical in many power systems 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.

ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Annual report for DOE VTO

Carbon fiber (CF)/polymer composites are a transformative class of high-performance, lightweight material, where high aspect-ratio CFs reinforce a polymer matrix and exceed the strength of steel alloys at a fraction of the density. Despite the advantages of such a class of material, the broader implementation of CF composites in a range of automotive, aerospace, and energy applications is hindered by limitations of current manufacturing methods. These current techniques (e.g., hand lay-up, wet filament winding) are costly and impose severe limitations on fiber placement, orientation, and angle, and thus a composite’s ultimate properties. Today’s CF composites are expensive to manufacture, limited in form factor, and utilize costly and sub-optimal continuous filament CF. Advanced additive manufacturing (AM) processes, combined with computational design optimization and new approaches to resin development, offer alternative design and manufacturing paradigms that have the realistic potential to lift these constraints. Such integrated AM approaches could thus help to realize the full potential of CF composite materials. One relevant application of CF composite materials where manufacturing constraints limit the cost-benefit ratio is in the manufacture of high-performance composite pressure vessels for onboard compressed natural gas (CNG) storage. Current CNG storage vessels (Types 3–5) are made from load-bearing filament-wound carbon-fiber composite and are ~3.5 times as expensive as an all-metallic Type-1 vessel. This cost is invariably tied to the complex and labor-intensive nature of conventional filament winding processes and the large volumes of expensive high tensile-strength CF tow feedstock required in manufacture. Our proposed approach to CNG storage vessel manufacture is based on a combination of AM technologies for CF composite printing and design optimization tools that were pioneered at Lawrence Livermore National Laboratory (LLNL), with advances in resin/composite formulation enabled by chemical and nano-material modification. Through the successful development of this technology, LLNL seeks to demonstrate the capability for advanced CNG storage vessel manufacture at reduced cost with no reduction in performance versus the most advanced, extant Type-5 designs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
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↗
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↗
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↗
DOE OSTI2020

Transformational Challenge Reactor preconceptual core design studies

In the nuclear industry, a manufacturing-informed design approach has the potential to yield the most benefit from advanced manufacturing. By leveraging advanced materials, data science, and rapid testing and deployment, manufacturing-informed design can drive down costs and development times, ultimately improving future commercial viability. This approach is being demonstrated in the US Department of Energy Office of Nuclear Energy (DOE-NE) Transformational Challenge Reactor (TCR) program. Preconceptual design activities for TCR have been focused on analyzing and maturing four reactor core design concepts: two fast-spectrum and two thermal-spectrum systems. Furthermore, the designs were iteratively modified and analyzed, and subcomponents were manufactured in parallel over weeks instead of months or years. To meet key program initiatives (e.g., timeline and material use), several constraints—including fissile material availability, component availability, materials compatibility, and additive manufacturing capabilities—were factored into the design effort, yielding small cores less than one cubic meter in volume with near-term viability. Additionally, the TCR program has made significant progress on development of advanced moderator materials such as yttrium hydride, advancing the feasibility of gas-cooled thermal spectrum systems using less than 250 kg of high-assay low enriched uranium (HALEU) and occupying less than 1 m3. Each of the two resulting thermal designs uses a different fuel form: traditional UO2 ceramic fuel and tristructural isotropic (advanced TRISO) fuel particles embedded inside a SiC matrix. Core neutronics and thermal performance for these systems were assessed and summarized. Evaluation of the performance metrics for these two moderated designs has yielded the downselected TCR design: a TRISO-fueled and yttrium hydride moderated gas-cooled reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Disentangling the Influences of Storm-Relative Flow and Horizontal Streamwise Vorticity on Low-Level Mesocyclones in Supercells

Sufficient low-level storm-relative flow is a necessary ingredient for sustained supercell thunderstorms and is connected to supercell updraft width. Assuming a supercell exists, the role of low-level storm-relative flow in regulating supercells’ low-level mesocyclone intensity is less clear. One possibility considered in this article is that storm-relative flow controls mesocyclone and tornado width via its modulation of overall updraft extent. This hypothesis relies on a previously postulated positive correspondence between updraft width, mesocyclone width, and tornado width. An alternative hypothesis is that mesocyclone characteristics are primarily regulated by horizontal streamwise vorticity irrespective of storm-relative flow. A matrix of supercell simulations was analyzed to address the aforementioned hypotheses, wherein horizontal streamwise vorticity and storm-relative flow were independently varied. Among these simulations, mesocyclone width and intensity were strongly correlated with horizontal streamwise vorticity, and comparatively weakly correlated with storm-relative flow, supporting the second hypothesis. Accompanying theory and trajectory analysis offers the physical explanation that, when storm-relative flow is large and updrafts are wide, vertically tilted streamwise vorticity is projected over a wider area but with a lesser average magnitude than when these parameters are small. These factors partially offset one another, degrading the correspondence of storm-relative flow with updraft circulation and rotational velocity, which are the mesocyclone attributes most closely tied to tornadoes. In conclusion, these results refute the previously purported connections between updraft width, mesocyclone width, and tornado width, and emphasize horizontal streamwise vorticity as the primary control on low-level mesocyclones in sustained supercells.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Easy-plane anisotropic-exchange magnets on a honeycomb lattice: Quantum effects and dealing with them

We provide analytical and numerical insights into the phase diagram and other properties of the extended Kitaev-Heisenberg model on the honeycomb lattice in the easy-plane limit, in which interactions are only between spin components that belong to the plane of magnetic ions. This parameter subspace allows for a much-needed systematic quantitative investigation of spin excitations in the ordered phases and of their generic features. Specifically, we demonstrate that in this limit one can consistently take into account magnon interactions in both zero-field zigzag and field-polarized phases. For the nominally polarized phase, we propose a regularization of the unphysical divergences that occur at the critical field and are plaguing the 1/S approximation in this class of models. For the explored parameter subspace, all symmetry-allowed terms of the standard parametrization of the extended Kitaev-Heisenberg model, such as K, J, and Γ, are significant, making the offered consideration relevant to a much wider parameter space. Furthermore, the dynamical structure factor near the paramagnetic critical point illustrates this relevance by showing features that are reminiscent of the ones observed in α–RuCl 3 , underscoring that they are not unique and should be common to a wide range of parameters of the model and, by extension, to other materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗