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

Compact representations of structured BFGS matrices

For general large-scale optimization problems compact representations exist in which recursive quasi-Newton update formulas are represented as compact matrix factorizations. For problems in which the objective function contains additional structure, recent structured quasi-Newton methods exploit available second-derivative information and approximate unavailable second derivatives. Here, this article develops the compact representations of two structured Broyden-Fletcher-Goldfarb-Shanno update formulas. The compact representations enable efficient limited memory and initialization strategies. Two limited memory line search algorithms are described for which extensive numerical results demonstrate the efficacy of the algorithms, including comparisons to IPOPT on large machine learning problems, and to L-BFGS on a real world large scale ptychographic imaging application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Radionuclide surrogate aerosolization, resuspension and suppression in hazardous situations

Radioactive fire following nuclear accidents causes resuspension of radioactive contamination that can spread over large areas. Although many papers focus on particle deposition and resuspension, a limited number of studies model its suppression. Even less data is available for the comparison of emissions across particle size, surface, and mode of resuspension. In this study, cerium dioxide surrogate for internal hazards such as radionuclide particle matter was aerosolized into a large chamber and subsequently resuspended after sedimentation. Resuspension was induced by common modes of human movement including walking, wind, and driving on common surfaces including polyvinyl chloride (PVC), artificial grass (turf), and concrete in a 5.8 m long, 3.7 m wide and 3 m high tent. The resuspension was experimentally measured for two particle sizes (1 and 10 μm) at three different heights (0.3 m, 1.2 m, and 2.4 m) and compared across size, surface, and mode of resuspension. Particle size distribution was measured using an aerodynamic particle sizer and the resuspended particles were collected with the wetted wall cyclone aerosol collector and quantitated using inductively coupled plasma mass spectrometry. The testing was repeated with the F-500 fire retardant application to study its effect on resuspension. In this study, artificial grass flooring was found to generate the most resuspension, with walking on turf for the 1 μm particle testing yielding the highest factor overall. Concrete had the lowest resuspension, showing variations in the size distribution of aerosol as a function of height from the source resuspension factors. The application of fire retardant was found to significantly reduce surrogate resuspension, regardless of resuspension method or testing surface. Here, the results provide impactful information regarding the different modes of resuspension including human transport (walking and driving) and wind for the resuspension of radionuclide surrogate particles deposited onto common surfaces and the effect of fire retardant application on their resuspension, leading to higher levels of safety for nuclear energy plants and facilities.

42 ENGINEERING↗
DOE OSTI2023

Optimization of operational strategies for rich gas enhanced oil recovery based on a pilot test in the Bakken tight oil reservoir

Horizontal well drilling and multistage hydraulic fracturing have been demonstrated as effective approaches for stimulating oil production in the Bakken tight oil reservoir. However, after multiple years of production, primary oil recovery in the Bakken is generally less than 10% of the estimated original oil in place. Gas huff ‘n’ puff (HnP) has been tested in the Bakken Formation as an enhanced oil recovery (EOR) method; however, most field pilot test results showed no significant incremental oil production. One of the factors affecting HnP EOR performance is premature gas breakthrough, which is one of the most critical issues observed in the field because of the presence of interwell fractures. Consequently, injected gas rapidly reaches adjacent production wells without contacting reservoir rock and increasing oil recovery. Proper conformance control is therefore needed to avoid early gas breakthrough and improve EOR performance. In this study, a rich gas EOR pilot in the Bakken was carefully analyzed to collect the essential reservoir and operational data. A simulation model with 16 wells was then developed to reproduce the production history and predict the EOR performance with and without conformance control. EOR operational strategies, including single- and multiple-well HnP, with different gas injection constraints were investigated. The simulation results of single-well HnP without conformance control showed that a rich gas injection rate of at least 10 MMscfd was needed to yield meaningful incremental oil production. The strategy of conformance control via water injection could significantly improve oil production in the HnP well, but injecting an excessive amount of water also leads to water breakthrough and loss of oil production in the offset wells. By analyzing the production performance of the wells individually, the arrangement of wells was optimized for multiple-well HnP EOR. The multiwell results showed that rich gas EOR could improve oil production up to 7.4% by employing conformance control strategies. Furthermore, replacing rich gas with propane as the injection gas could result in 14% of incremental oil production.

02 PETROLEUM↗
DOE OSTI2021

An investigation of hard-disk drive circularity accounting for socio-technical dynamics and data uncertainty

The installed data storage capacity in the U.S. will reach 2.2 Zettabytes by 2025, generating about 50 million units of end-of-life (EOL) hard-disk drives (HDDs) per year. Due to data security concerns, most EOL HDDs are currently shredded (even when still functioning), representing an economic loss. Moreover, raw material extraction linked to the increased demand for storage causes environmental impacts. Besides mitigating the threat posed by sudden restrictions of raw materials, the circular economy (CE) offers to maximize value retention in the economy and reduce the environmental impacts of human activities. Common CE strategies are reusing and recycling products. However, the reuse of hard disk drives is currently burdened by the lack of trust HDD end-users have toward other non-physical means of data removal than shredding. Here, an agent-based modeling (ABM) approach is proposed to explore how techno-economic and social factors affect end-users' decisions to adopt EOL management practices other than shredding. The proposed method also accounts for data uncertainty by applying a semi-quantitative approach. Results demonstrate how increased green procurement and more robust standards could spur end-users' trust toward data-wiping technologies. Even when accounting for uncertainty, HDDs' reuse brings better environmental and economic benefits than HDD shredding followed by material recovery. The semi-quantitative approach proposed in this study could be more universally applied in future ABM, especially given the often-stochastic nature of such models. The developed ABM is also the first to represent several HDD industry stakeholders and demonstrate how the HDD shredding lock-in situation could be resolved.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Polishing the Gold Standard: The Role of Orbital Choice in CCSD(T) Vibrational Frequency Prediction

While CCSD(T) with spin-restricted Hartree-Fock (RHF) orbitals has long been lauded for its ability to accurately describe closed-shell interactions, the performance of CCSD(T) on open-shell species is much more erratic, especially when using a spin-unrestricted HF (UHF) reference. Previous studies have shown improved treatment of open-shell systems when a non-HF set of molecular orbitals, like Brueckner or Kohn-Sham density functional theory (DFT) orbitals, is used as a reference. Inspired by the success of regularized orbital-optimized second-order Møller-Plesset perturbation theory (κ-OOMP2) orbitals as reference orbitals for MP3, we investigate the use of κ-OOMP2 orbitals and various DFT orbitals as reference orbitals for CCSD(T) calculations of the corrected ground-state harmonic vibrational frequencies of a set of 36 closed-shell (29 neutrals, 6 cations, 1 anion) and 59 open-shell diatomic species (38 neutrals, 15 cations, 6 anions). The aug-cc-pwCVTZ basis set is used for all calculations. The use of κ-OOMP2 orbitals in this context alleviates difficult cases observed for both UHF orbitals and OOMP2 orbitals. Removing two multireference systems and 12 systems with ambiguous experimental data leaves a pruned data set. Overall performance on the pruned data set highlights CCSD(T) with a B97 orbital reference (CCSD(T):B97), CCSD(T) with a κ-OOMP2 orbital reference (CCSD(T):κ-OOMP2), and CCSD(T) with a B97M-rV orbital reference (CCSD(T):B97M-rV) with RMSDs of 8.48 cm -1 , and 8.50 cm -1 , and 8.75 cm -1 respectively, outperforming CCSD(T):UHF by nearly a factor of 5. Moreover, the performance on the closed- and open-shell subsets shows these methods are able to treat open-shell and closed-shell systems with comparable accuracy and robustness. CCSD(T) with RHF orbitals is seen to improve upon UHF for the closed-shell species, while spatial symmetry breaking in a number of restricted open-shell HF (ROHF) references leads CCSD(T) with ROHF reference orbitals to exhibit the poorest statistical performance of all methods surveyed for open-shell species. The use of κ-OOMP2 orbitals has also proven useful in diagnosing multireference character that can hinder the reliability of CCSD(T).

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

Connectivity-Dependent Exciton–Phonon Coupling in Cesium Bismuth Halide Quantum Dots

Metal halide octahedra form the fundamental functional building blocks of metal halide perovskites, dictating their structures, optical properties, electronic structures, and dynamics. Here, in this study, we show that the connectivity of bismuth halide octahedra in Cs 3 Bi 2 Br 9 and Cs 3 Bi 2 I 9 quantum dots (QDs) changes with different halide elements. We use first-principles calculations to reveal the key role of the connectivity of bismuth halide octahedra on the wave function symmetry, Huang-Rhys factor, and exciton-phonon interaction strength. Following QD synthesis via a ligand-mediated transport method, the effect of connectivity is verified with transient absorption spectroscopy, where we contrast Cs 3 Bi 2 Br 9 and Cs 3 Bi 2 I 9 QD exciton dynamics. In photoexcited Cs 3 Bi 2 I 9 QDs, phonons related to the vibrational motions of face-sharing [BiI 6 ] 3- bioctahedra couple strongly to the electronic state and drive rapid carrier relaxation. Equivalent signals are not observed for photoexcited Cs 3 Bi 2 Br 9 QDs, implying a lack of phonon involvement in band-edge absorption and subsequent exciton relaxation. Our findings suggest that structural engineering can effectively tune the exciton-phonon coupling and therefore influence exciton relaxation and recombination in perovskite nanomaterials.

TDDFT↗
DOE OSTI2020

Controlled Assembly of Upconverting Nanoparticles for Low-Threshold Microlasers and Their Imaging in Scattering Media

Micron-sized lasers fabricated from upconverting nanoparticles (UCNP) coupled to whispering gallery mode (WGM) microresonators can exhibit continuous-wave anti-Stokes lasing useful for tracking cells, environmental sensing, and coherent stimulation of biological activity. The integration of these microlasers into organisms and microelectronics requires even smaller diameters, however, which raises threshold pump powers beyond practical limits for biological applications. To meet the need for low lasing thresholds and high fidelity fabrication methods, we use correlative optical and electron microscopy to uncover the nanoparticle assembly process and structural factors that determine efficient upconverted lasing. We show that 5 μm microspheres with controlled submonolayer UCNP coatings exhibit, on average, 25-fold lower laser thresholds (1.7 ± 0.7 kW/cm 2 ) compared to the mean values of the lowest threshold UCNP lasers, and variability is reduced 30-fold. WGMs are observed in the upconversion spectra for TiO 2 -coated microspheres as small as 3 μm, a size at which optical losses had previously prevented such observations. Finally, we demonstrate that the WGM signatures of these upconverting microlasers can be imaged and distinguished through tissue-mimicking phantoms. These advances will enable the fabrication of more efficient upconverting lasers for imaging, sensing, and actuation in optically complex environments.

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

Estimating the CO 2 Fertilization Effect on Extratropical Forest Productivity From Flux‐Tower Observations

Abstract The land sink of anthropogenic carbon emissions, a crucial component of mitigating climate change, is primarily attributed to the CO 2 fertilization effect on global gross primary productivity (GPP). However, direct observational evidence of this effect remains scarce, hampered by challenges in disentangling the CO 2 fertilization effect from other long‐term confounding drivers, particularly climatic changes. Here, we introduce a novel statistical approach to separate the CO 2 fertilization effect on photosynthetic carbon uptake using eddy covariance (EC) records across 38 extratropical forest sites. We find the median stimulation rate of GPP to be 3.2 ± 0.9 gC m −2 yr −1 ppm −1 (or 16.4 ± 4.2% per 100 ppm) under increasing atmospheric CO 2 across these sites, respectively. To validate the robustness of our findings, we test our statistical method using factorial simulations of an ensemble of process‐based land surface models. We address additional factors, including nitrogen deposition and land management, that may impact plant productivity, potentially confounding the attribution to the CO 2 fertilization effect. Assuming these site‐specific effects offset to some extent across sites as random factors, the estimated median value still reflects the strength of the CO 2 fertilization effect. However, disentanglement of these long‐term effects, often inseparable by timescale, requires further causal research. Our study provides direct evidence that the photosynthetic stimulation is maintained under long‐term CO 2 fertilization across multiple EC sites. Such observation‐based quantification is key to constraining the long‐standing uncertainties in the land carbon cycle under rising CO 2 concentrations.

Environmental Sciences & Ecology↗
DOE OSTI2020

Structure retrieval in liquid-phase electron scattering

Electron scattering on liquid samples has been enabled recently by the development of ultrathin liquid sheet technologies. The data treatment of liquid-phase electron scattering has been mostly reliant on methodologies developed for gas electron diffraction, in which theoretical inputs and empirical fittings are often needed to account for the atomic form factor and remove the inelastic scattering background. In this work, we present an alternative data treatment method that is able to retrieve the radial distribution of all the charged particle pairs without the need of either theoretical inputs or empirical fittings. Here, the merits of this new method are illustrated through the retrieval of real-space molecular structure from experimental electron scattering patterns of liquid water, carbon tetrachloride, chloroform, and dichloromethane.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Monte Carlo evaluation of the semiclassical multi-step direct reaction series

We show how the multi-step direct reaction series can be evaluated using Monte Carlo methods. Nucleon-nucleon collisions occur according to the random selection of a nucleon’s attenuation factor along its classical trajectory. The particles and hole excited in a collision are selected from a local Fermi distribution. We assume that the particles continue to propagate and possibly collide again before leaving the nucleus. We assume that holes collide in place, to possibly produce other particles and holes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Experimental and computational studies of the kinetics of the reaction of hydrogen peroxide with the amidogen radical

In this work, the pulsed-laser photolysis/laser-induced fluorescence method is used to study the kinetics of the reaction of NH 2 with H 2 O 2 to yield a second-order rate constant of (2.42 ± 0.55) × 10 –14 cm 3 molecule –1 s –1 at 412 K in 10–22 mbar of Ar bath gas. There are no prior measurements for comparison. To check this value and enable reliable extrapolation to other temperatures, we also compute thermal rate constants for this process over the temperature range 298–3000 K via multi-structural canonical variational transition-state theory with small-curvature multidimensional tunneling (MS-CVT/SCT). The CVT/SCT rate constants are derived using a dual-level direct dynamics approach utilizing single-point CCSD(T)-F12b/cc-pVQZ-F12 energies—corrected for core-valence and scalar relativistic effects—and M06-2X/MG3S geometries, gradients, and Hessians—for all stationary and non-stationary points along the reaction path. The multistructural method with torsional anharmonicity, based on a coupled torsional potential, is then employed to calculate correction factors for the rate constants, accounting for the comprehensive effects of torsional anharmonicity on the kinetics of this reaction system. The final MS-CVT/SCT rate constants are found to be in good agreement with our measurements and can be expressed in modified Arrhenius form as 2.13 × 10 –15 (T/298 K) 4.02 exp(–513 K/T) cm 3 molecule –1 s –1 over the temperature range of 298–3000 K.

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

An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe–Salpeter equation calculations of molecules

The ab initio GW plus Bethe–Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn–Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm–Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn–Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel’s set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.

McKeon, Caroline A. (ORCID:0000000217373503)↗
DOE OSTI2021

The inflationary wavefunction from analyticity and factorization

Here, we study the analytic properties of tree-level wavefunction coefficients in quasi-de Sitter space. We focus on theories which spontaneously break dS boost symmetries and can produce significant non-Gaussianities. The corresponding inflationary correlators are (approximately) scale invariant, but are not invariant under the full conformal group. We derive cutting rules and dispersion formulas for the late-time wavefunction coefficients by using factorization and analyticity properties of the dS bulk-to-bulk propagator. This gives a unitarity method which is valid at tree-level for general n-point functions and for fields of arbitrary mass. Using the cutting rules and dispersion formulas, we are able to compute n-point functions by gluing together lower-point functions. As an application, we study general four-point, scalar exchange diagrams in the EFT of inflation. We show that exchange diagrams constructed from boost-breaking interactions can be written as a finite sum over residues. Finally, we explain how the dS identities used in this work are related by analytic continuation to analogous identities in Anti-de Sitter space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

NASMDR: a framework for miRNA-drug resistance prediction using efficient neural architecture search and graph isomorphism networks

Abstract As a frontier field of individualized therapy, microRNA (miRNA) pharmacogenomics facilitates the understanding of different individual responses to certain drugs and provides a reasonable reference for clinical treatment. However, the known drug resistance-associated miRNAs are not yet sufficient to support precision medicine. Although existing methods are effective, they all focus on modelling miRNA-drug resistance interaction graphs, making their performance bounded by the interaction density. In this study, we propose a framework for miRNA-drug resistance prediction through efficient neural architecture search and graph isomorphism networks (NASMDR). NASMDR uses attribute information instead of the commonly used interactive graph information. In the cross-validation experiment, the proposed framework can achieve an AUC of 0.9468 on the ncDR dataset, which is 2.29% higher than the state-of-the-art method. In addition, we propose a novel sequence characterization approach, k-mer Sparse Nonnegative Matrix Factorization (KSNMF). The results show that NASMDR provides novel insights for integrating efficient neural architecture search and graph isomorphic networks into a unified framework to predict drug resistance-related miRNAs. The codes for NASMDR are available at https://github.com/kaizheng-academic/NASMDR.

Zheng, Kai↗
DOE OSTI2022

Genetic-algorithm retrieval of the molecular alignment distribution with high-order harmonics generated from transiently aligned CO 2 molecules

Using genetic algorithm, we propose a method to retrieve the alignment distribution of transiently aligned CO 2 molecules from the high-order harmonic generation (HHG) spectra. The retrieval method is based on the quantitative rescattering (QRS) model where averaged photorecombination transition dipole can be factored out from the parallel (or perpendicular) harmonic spectra after the propagation of the harmonic fields in the gas medium. We examine how the retrieved alignment distributions are affected by uncertainty in alignment dependent ionization probability and on multiple orbital contribution to the HHG. We further confirm that alignment distribution is more accurately retrieved by using the minima in the HHG spectra driven by a long-wavelength laser. In addition, we show that earlier experimental data on the ratios between the perpendicular and the parallel HHG components of aligned CO 2 molecules are in better agreement with the QRS model if the macroscopic propagation and multiple orbital interference are included in the theoretical calculation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Predicting band gaps and band-edge positions of oxide perovskites using density functional theory and machine learning

Density functional theory (DFT) within the local or semilocal density approximations, i.e., the local density approximation (LDA) or generalized gradient approximation (GGA), has become a workhorse in the electronic structure theory of solids, being extremely fast and reliable for energetics and structural properties, yet remaining highly inaccurate for predicting band gaps of semiconductors and insulators. The accurate prediction of band gaps using first-principles methods is time consuming, requiring hybrid functionals, quasiparticle GW, or quantum Monte Carlo methods. Efficiently correcting DFT-LDA/GGA band gaps and unveiling the main chemical and structural factors involved in this correction is desirable for discovering novel materials in high-throughput calculations. In this direction, we, in this study, use DFT and machine learning techniques to correct band gaps and band-edge positions of a representative subset of ABO 3 perovskite oxides. Relying on the results of HSE06 hybrid functional calculations as target values of band gaps, we find a systematic band-gap correction of ~1.5 eV for this class of materials, where ~1 eV comes from downward shifting the valence band and ~0.5 eV from uplifting the conduction band. The main chemical and structural factors determining the band-gap correction are determined through a feature selection procedure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Bifurcation Buckling Analysis and Non-Linear Collapse Analysis of Teardrop Shaped Vacuum Chamber

The Spallation Neutron Source (SNS) accelerator is being upgraded to increase the beam power from 1.4MW at 1GeV to 2.8MW at 1.3GeV. The currents in the middle two injection chicane magnets cannot simply be scaled up to accommodate the increased injection energy of 1.3GeV due to potential excessive H- stripping; the magnets must be replaced with longer, lower-field magnets and the associated vacuum chambers need to be redesigned. A new teardrop-shaped vacuum chamber was initially designed to accommodate the new magnets and the updated beam paths and instrumentation. This paper focuses on the structural stability study of the teardrop shape vacuum chamber based on buckling analysis. Protection against collapse from buckling according to the ASME BPVC requirement has been evaluated in depth. First, a Type-1 bifurcation buckling analysis using a linear eigenvalue solution to determine the critical load factor was performed. Subsequently, a Type-3 nonlinear collapse analysis was conducted using the static Riks method with elastic-plastic material properties and imperfections explicitly considered in the model geometry. The critical buckling load for the teardrop shape vacuum chamber was confidently estimated based upon this two-stage approach.

Jiang, Hao↗