Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “factorization methods”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Engineering chirality at wafer scale with ordered carbon nanotube architectures

Creating artificial matter with controllable chirality in a simple and scalable manner brings new opportunities to diverse areas. Here we show two such methods based on controlled vacuum filtration - twist stacking and mechanical rotation - for fabricating wafer-scale chiral architectures of ordered carbon nanotubes (CNTs) with tunable and large circular dichroism (CD). By controlling the stacking angle and handedness in the twist-stacking approach, we maximize the CD response and achieve a high deep-ultraviolet ellipticity of 40 ± 1 mdeg nm -1 . Our theoretical simulations using the transfer matrix method reproduce the experimentally observed CD spectra and further predict that an optimized film of twist-stacked CNTs can exhibit an ellipticity as high as 150 mdeg nm -1 , corresponding to a g factor of 0.22. Furthermore, the mechanical rotation method not only accelerates the fabrication of twisted structures but also produces both chiralities simultaneously in a single sample, in a single run, and in a controllable manner. The created wafer-scale objects represent an alternative type of synthetic chiral matter consisting of ordered quantum wires whose macroscopic properties are governed by nanoscopic electronic signatures and can be used to explore chiral phenomena and develop chiral photonic and optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Acceleration Factors for Combined‐Accelerated Stress Testing of Photovoltaic Modules

Combined‐accelerated stress testing (C‐AST) is developed to establish the durability of photovoltaic (PV) products, including for degradation modes that are not a priori known or examined in standardized tests. C‐AST aims to comprehensively represent the sample, stress factors, and their combinations using levels at the statistical tails of the natural environment. Acceleration factors for relevant climate sequences within the C‐AST cycle with respect to the Florida USA climate are estimated for selected degradation mechanisms. It is found that for degradation of the outer backsheet polymer layer, the acceleration factor of the tropical climate sequence (the longest of the climate sequences) is f ( T , G ) = 17.3 with ultraviolet photodegradation; for polyethylene terephthalate hydrolysis (backsheets), f ( T , RH ) = 426; for electrochemical corrosion (PV cell), f ( I ) = 14.1; and for PbSn solder fatigue f (Δ T , r ( T )) = 17.3. Here, T is the module temperature, G is the broadband spectrum irradiance on the plane of array of the module, RH is the relative humidity on the module surface, I is the leakage current through the module packaging, and r ( T ), the number of temperature reversals. The methods discussed herein are generally applicable for evaluating acceleration factors in other accelerated test methods.

14 SOLAR ENERGY↗
DOE OSTI2020

Prediction of Climate-Specific Degradation Rate for Photovoltaic Encapsulant Discoloration

Encapsulant discoloration is a well-known field degradation mode of crystalline-silicon photovoltaic modules, particularly in the hot climate zones. The discoloration rate is influenced by several weathering factors, such as UV light, module temperature, and humidity, as well as the permeability of oxygen into the module. Here, a rate dependence model employing the modified Arrhenius equations to predict the degradation rate for encapsulant discoloration in different climates is presented. Two modeling approaches are introduced, which utilize the field and accelerated UV testing degradation data in conjunction with the field meteorological data to determine the acceleration factor for encapsulant browning. A novel method of accelerated UV stress testing at three simultaneous module temperatures in a single environmental chamber test run is implemented to estimate the activation energy for browning. The test was performed on three field-retrieved modules to capture the wear-out failure mechanism. The degradation in short-circuit current Isc rather than maximum power is used as a decisive parameter for the discoloration analysis. Furthermore, the developed model has been used to predict the Isc degradation rate for the Arizona field characterized by a hot and dry climate and is validated against the field-measured value. It has also been applied to other climate types, e.g., the cold and dry climate of New York.

14 SOLAR ENERGY↗
DOE OSTI2023

A Decision-Relevant Factor-Fixing Framework: Application to Uncertainty Analysis of a High-Dimensional Water Quality Model

Factor Fixing (FF) is a common method for reducing the number of model parameters to lower computational cost. FF typically starts with distinguishing the insensitive parameters from the sensitive and pursues uncertainty quantification (UQ) on the resulting reduced-order model, fixing each insensitive parameter at a fixed value. There is a need, however, to expand such a common approach to consider the effects of decision choices in the FF-UQ procedure on metrics of interest. Therefore, to guide the use of FF and increase confidence in the resulting dimension-reduced model, we propose a new adaptive framework consisting of four principles: (a) re-parameterize the model first to reduce obvious non-identifiable parameter combinations, (b) focus on decision relevance especially with respect to errors in quantities of interest (QoI), (c) conduct adaptive evaluation and robustness assessment of errors in the QoI across FF choices as sample size increases, and (d) reconsider whether fixing is warranted. The framework is demonstrated on a spatially-distributed water quality model. The error in estimates of QoI caused by FF can be estimated using a Polynomial Chaos Expansion (PCE) surrogate model. Built with 70 model runs, the surrogate is computationally inexpensive to evaluate and can provide global sensitivity indices for free. For the selected catchment, just two factors may provide an acceptably accurate estimate of model uncertainty in the average annual load of Total Suspended Solids (TSS), suggesting that reducing the uncertainty in these two parameters is a priority for future work before undertaking further formal uncertainty quantification.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A polarimetry-based field-deployable non-interruptive mirror soiling detection method

The soiling level of heliostat mirrors in Concentrated Solar Power (CSP) fields is one of the key factors that significantly influences optical efficiency. State-of-the-art methods of monitoring heliostats soiling levels still face various challenges, including slow speed, labor-intensive operations, resolution and accuracy constraints or interruptions to solar field operations. Here, we present a rapid, cost-effective, and non-intrusive method for mirror soiling detection based on polarimetric imaging, referred to as Polarimetric Imaging-based Mirror Soiling (PIMS). The compact PIMS device is designed for integration with unmanned aerial vehicles (UAVs), enabling rapid, large-area assessments of heliostat mirrors for efficient soiling detection. Our method utilizes the correlation between the Degree of Linear Polarization (DoLP) and surface soiling level based on Mie scattering theory and Monte Carlo simulations. Field deployment of the PIMS method requires minimal device installation, and its UAV-based operation allows for soiling detection without interrupting plant activities. The PIMS method holds the potential for mirror soiling detection across various concentrated solar power (CSP) plants and can be further adapted for other types of solar fields, such as parabolic trough systems.

CSP Field↗
DOE OSTI2022

Spectral dispersion of the linewidth enhancement factor and four wave mixing conversion efficiency of an InAs/GaAs multimode quantum dot laser

The spectral dependence of the linewidth enhancement factor (αH-factor) of a multimode InAs/GaAs quantum dot laser is analyzed. Amplified spontaneous and high-frequency modulation methods are used to experimentally retrieve the αH-factor of each longitudinal mode below and above the threshold. A dispersion of the αH-factor is unlocked across the entire optical spectrum, which is further illustrated in the context of four wave mixing experiments. The results show that the induced conversion efficiency is increased at lasing wavelengths where the linewidth enhancement is lower. These results highlight the importance of carefully monitoring the linewidth enhancement factor in quantum dot lasers especially for frequency combs and mode-locking applications in future optical communication systems.

Ding, Shihao (ORCID:0000000346416944)↗
DOE OSTI2023

An adaptive synchronous extraction (ASE) method for estimating intensity and footprint of surface urban heat islands: A case study of 254 North American cities

The urban heat island (UHI) effect has attracted great attention due to its potential impacts on rapidly growing urban areas. Using remotely sensed estimates of land surface temperature (LST), a large number of studies have focused on the surface UHI (SUHI) effect, which can be characterized by its two fundamental properties: intensity and footprint. The SUHI intensity reflects the LST difference between the urban area and the background reference area (BRA), and the SUHI footprint indicates the spatial extent influenced by the heat island. Currently, numerous methods have been developed to estimate the SUHI intensity and footprint, but are still greatly challenged by three main issues. Namely, the discrepancy in BRA selection criterion brings great uncertainty to the estimated SUHI intensity, the estimation of SUHI footprint is largely constrained by the predefined models, and the quantification of SUHI effect is potentially influenced by several confounding factors. Here, we proposed an adaptive synchronous extraction (ASE) method, which is capable of adaptively selecting the most optimal BRA while removing the influence of confounding factors, and achieving synchronous estimation of SUHI intensity and footprint. We applied the ASE method to 254 North American cities and conducted an in-depth comparative analysis to discuss its applicability and benefits. The main results include: (1) The ASE method avoids the limitations of existing methods in BRA selection and model presetting, and shows resilience to parameter variations. This makes the ASE method highly applicable to quantify the SUHI intensity and footprint in cities with various thermal characteristics. (2) The ASE method can better highlight the spatial, seasonal and day-night contrasts in the estimated SUHI intensity. This superiority is particularly evident when comparing it to methods based on the equal-area buffer or the simplified urban-extent algorithm. (3) Confounding factors pose non-negligible impacts on the quantification of the SUHI effect. Typically, ignoring the influence of topographic relief or missing LST data can lead to an overall overestimation of the SUHI intensity, while not removing surrounding urban areas will cause some underestimation of the SUHI intensity. In conclusion, overall, the proposed ASE method provides a new generalizable tool for quantifying the SUHI effect, which has great potentials for future studies and urban climate assessments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through two methods, thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine the thermal quadrant powers. The N-16 system utilizes a beta chamber detector that outputs reactivity levels to calculate lobe power through an algorithm called lobe power calculation and indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. This is crucial for maintaining safety limits. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have assisted system engineers at INL with the objective to reach a more dynamic system. In return the system demonstrates an increase in the accuracy of power reading while maintaining safety margins. The new proposed coefficient and multiple method implements a cycle specific, and a physics-based model intended for changing coefficients during operation. This calibration experiment focuses on power splitting, outer shim and neck shim, as well as fuel burning into the RDAS weighting factors. Results show that the Physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. This proved true for both constrained and unconstrained testing. This is most likely due to the physics data fitting approach that resulted in favorable coefficients and multipliers for 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Validating and Comparing Energy Estimation Methods at Water Resource Recovery Facilities

Water resource recovery facilities play a crucial role in the water-energy nexus, consuming a substantial amount of energy in the United States. Growing treatment volumes and more stringent water quality standards are expected to increase the amount of energy needed to treat wastewater, but accurately estimating energy consumption and potential remains challenging due to variability in scale, treatment methods, and effluent treatment standards. In this study, we used publicly available data to evaluate the accuracy of methods for estimating energy consumption and generation, then quantified uncertainty based on key factors like flow rate, treatment level, and geographic location. To validate methods, we estimated energy consumption and generation at the facility-level, then compared estimates to self-reported data from utilities in major U.S. cities. We found that process models of treatment trains under best practice configurations were accurate relative to other methods for estimating electricity use, total energy use, and electricity generation from biogas utilization, and less complex methods based on effluent treatment level and prime movers also performed well for estimating electricity consumption and generation, respectively. Applying the evaluated methods to a national inventory of treatment facilities, we estimate that annual energy consumption ranged from 56.3 x 10^3 to 82.5 x 10^3 TJ in 2012 and 83.6 x 10^3 to 127 x 10^3 TJ in 2042. Our results indicate that not all estimation methods are suited for every use case, so we recommend that researchers and practitioners select an estimation method based on data availability and desired computational intensity.

Hodson, Abigayle↗
DOE OSTI2021

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through various methods, two of them being thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine thermal quadrant powers. The N-16 system utilizes a beta detector that outputs Nitrogen activation levels to calculate lobe power through an algorithm called Lobe Power Calculation and Indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. These correct indications are crucial for maintaining safety limits because operators rely on this information for decision making. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have increased the accuracy of power readings by making the system more dynamic. The new proposed coefficient and multiplier method implements a cycle specific and physics-based model intended for changing coefficients during operation. This calibration experiment focused on power splitting, outer shim and neck shim, as well as fuel burning into the reactor digital acquisition system (RDAS) weighting factors. Results demonstrated that the physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Cost-efficient copper-nickel alloy for active cooling applications

It is recently proposed that active (Peltier) and passive (conductive) cooling could be combined to design highly-efficient thermal conductors as active heat sinks. Here we study the potential of low-cost copper-nickel alloys for active cooling applications. We show that for this application, copper-nickel alloys prepared using industrial-grade powders to lower the manufacturing cost, are as good as copper-nickel alloys made out of high purity metals. Furthermore, directed energy deposition, a type of additive manufacturing, can produce copper-nickel alloys which have about 10 percent lower thermoelectric power factor compared to mechanical-alloying methods but have larger thermal conductivity values, hence, possessing similar effective thermal conductivity values, reaching 403 W m -1 K -1 at 390K. Much larger values are expected at higher temperatures. The usage of additive manufacturing enables complex geometries and large-scale production.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Towards green synthesis of Mn 4+ -doped fluoride phosphors: a review

Mn 4+ -doped fluoride phosphors are attractive for application in white solid-state lighting devices because their distinct red emission can improve the color rendering index, while the earth-abundance of Mn 4+ ions can reduce the production cost. However, highly hazardous HF is used as the fluoride source in the synthesis processes, which requires extreme caution. In this review, we introduce various methods for synthesizing Mn 4+ -doped fluoride phosphors, categorized according to the presence or absence of HF in the synthesis. The underlying optical properties of the phosphors fabricated using the various methods are discussed. In addition, in the discussion of each method, we emphasize on the key factors and strategies for developing phosphors. Finally, this review highlights perspectives on and future directions for the safe synthesis of Mn 4+ -doped fluoride phosphors to achieve solid-state lighting devices with superior optical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗
DOE OSTI2021

Insights into the emergence of mass from studies of pion and kaon structure

There are two mass generating mechanisms in the standard model of particle physics (SM). One is related to the Higgs boson and fairly well understood. The other is embedded in quantum chromodynamics (QCD), the SM’s strong interaction piece; and although responsible for emergence of the roughly 1 GeV mass scale that characterises the proton and hence all observable matter, the source and impacts of this emergent hadronic mass (EHM) remain puzzling. As bound states seeded by a valence-quark and -antiquark, pseudoscalar mesons present a simpler problem in quantum field theory than that associated with the nucleon. Consequently, there is a large array of robust predictions for pion and kaon properties whose empirical validation will provide a clear window onto many effects of both mass generating mechanisms and the constructive interference between them. This has now become significant because new-era experimental facilities, in operation, construction, or planning, are capable of conducting such tests and thereby contributing greatly to resolving the puzzles of EHM. These aspects of experiment, phenomenology, and theory, along with contemporary successes and challenges, are reviewed herein. In addition to providing an overview of the experimental status, we focus on recent progress made using continuum Schwinger function methods and lattice-regularised QCD. Further, advances made using other theoretical tools are also sketched. Our primary goal is to highlight the potential gains that can accrue from a coherent effort aimed at finally reaching an understanding of the character and structure of Nature’s Nambu–Goldstone modes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A flexible class of priors for orthonormal matrices with basis function-specific structure

Statistical modeling of high-dimensional matrix-valued data motivates the use of a low-rank representation that simultaneously summarizes key characteristics of the data and enables dimension reduction. Low-rank representations commonly factor the original data into the product of orthonormal basis functions and weights, where each basis function represents an independent feature of the data. However, the basis functions in these factorizations are typically computed using algorithmic methods that cannot quantify uncertainty or account for basis function correlation structure a priori. While there exist Bayesian methods that allow for a common correlation structure across basis functions, empirical examples motivate the need for basis function-specific dependence structure. We propose a prior distribution for orthonormal matrices that can explicitly model basis function-specific structure. The prior is used within a general probabilistic model for singular value decomposition to conduct posterior inference on the basis functions while accounting for measurement error and fixed effects. We discuss how the prior specification can be used for various scenarios and demonstrate favorable model properties through synthetic data examples. Finally, we apply our method to two-meter air temperature data from the Pacific Northwest, enhancing our understanding of the Earth system’s internal variability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Telecommunication-wavelength two-dimensional photonic crystal cavities in a thin single-crystal diamond membrane

We demonstrate two-dimensional photonic crystal cavities operating at telecommunication wavelengths in a single-crystal diamond membrane. We use a high-optical-quality and thin (~300 nm) diamond membrane, supported by a polycrystalline diamond frame, to realize fully suspended two-dimensional photonic crystal cavities with a high theoretical quality factor of ~8 × 10 6 and a relatively small mode volume of ~2(λ/n) 3 . The cavities are fabricated in the membrane using electron-beam lithography and vertical dry etching. We observe cavity resonances over a wide wavelength range spanning the telecommunication O- and S-bands (1360–1470 nm) with Q factors of up to ~1800. Our method paves the way for on-chip diamond nanophotonic applications in the telecommunication-wavelength range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Energy-conserving coupled trajectory mixed quantum–classical dynamics

The coupled-trajectory mixed quantum–classical method (CTMQC), derived from the exact factorization approach, has successfully predicted photo-chemical dynamics in a number of interesting molecules, capturing population transfer and decoherence from first principles. Furthermore, due to the approximations made, CTMQC does not guarantee energy conservation. We propose a modified algorithm, CTMQC-E, which redefines the integrated force in the coupled-trajectory term so to restore energy conservation, and demonstrate its accuracy on scattering in Tully’s extended coupling region model and photoisomerization in a retinal chromophore model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗