Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistency”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

Self-consistent implementation of locally scaled self-interaction-correction method

Recently proposed local self-interaction correction (LSIC) method is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange–correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew–Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizsäcker and Kohn–Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-Löwdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of the AE6 dataset better than some of the most widely used generalized-gradient-approximation (GGA) functional and barrier heights of the BH6 database better than some of the most widely used hybrid functionals. The LSIC method [a mean absolute error (MAE) of 0.008 Å] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 Å) and LSDA (0.011 Å). Furthermore, this work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗

A Consistent Phylogenetic Backbone for the Fungi

The kingdom of fungi provides model organisms for biotechnology, cell biology, genetics, and life sciences in general. Only when their phylogenetic relationships are stably resolved, can individual results from fungal research be integrated into a holistic picture of biology. However, and despite recent progress, many deep relationships within the fungi remain unclear. Here, we present the first phylogenomic study of an entire eukaryotic kingdom that uses a consistency criterion to strengthen phylogenetic conclusions. We reason that branches (splits) recovered with independent data and different tree reconstruction methods are likely to reflect true evolutionary relationships. Two complementary phylogenomic data sets based on 99 fungal genomes and 109 fungal expressed sequence tag (EST) sets analyzed with four different tree reconstruction methods shed light from different angles on the fungal tree of life. Eleven additional data sets address specifically the phylogenetic position of Blastocladiomycota, Ustilaginomycotina, and Dothideomycetes, respectively. The combined evidence from the resulting trees supports the deep-level stability of the fungal groups toward a comprehensive natural system of the fungi. In addition, our analysis reveals methodologically interesting aspects. Enrichment for EST encoded data—a common practice in phylogenomic analyses—introduces a strong bias toward slowly evolving and functionally correlated genes. Consequently, the generalization of phylogenomic data sets as collections of randomly selected genes cannot be taken for granted. A thorough characterization of the data to assess possible influences on the tree reconstruction should therefore become a standard in phylogenomic analyses.

59 BASIC BIOLOGICAL SCIENCES↗

Elliptic Preconditioner for Accelerating the Self-Consistent Field Iteration in Kohn--Sham Density Functional Theory

We discuss techniques for accelerating the self consistent field (SCF) iteration for solving the Kohn-Sham equations. These techniques are all based on constructing approximations to the inverse of the Jacobian associated with a fixed point map satisfied by the total potential. They can be viewed as preconditioners for a fixed point iteration. We point out different requirements for constructing preconditioners for insulating and metallic systems respectively, and discuss how to construct preconditioners to keep the convergence rate of the fixed point iteration independent of the size of the atomistic system. We propose a new preconditioner that can treat insulating and metallic system in a unified way. The new preconditioner, which we call an elliptic preconditioner, is constructed by solving an elliptic partial differential equation. The elliptic preconditioner is shown to be more effective in accelerating the convergence of a fixed point iteration than the existing approaches for large inhomogeneous systems at low temperature.

97 MATHEMATICS AND COMPUTING↗

A Shift Selection Strategy for Parallel Shift-invert Spectrum Slicing in Symmetric Self-consistent Eigenvalue Computation

The central importance of large-scale eigenvalue problems in scientific computation necessitates the development of massively parallel algorithms for their solution. Recent advances in dense numerical linear algebra have enabled the routine treatment of eigenvalue problems with dimensions on the order of hundreds of thousands on the world’s largest supercomputers. In cases where dense treatments are not feasible, Krylov subspace methods offer an attractive alternative due to the fact that they do not require storage of the problem matrices. However, demonstration of scalability of either of these classes of eigenvalue algorithms on computing architectures capable of expressing massive parallelism is non-trivial due to communication requirements and serial bottlenecks, respectively. In this work, we introduce the SISLICE method: a parallel shift-invert algorithm for the solution of the symmetric self-consistent field (SCF) eigenvalue problem. The SISLICE method drastically reduces the communication requirement of current parallel shift-invert eigenvalue algorithms through various shift selection and migration techniques based on density of states estimation and k-means clustering, respectively. This work demonstrates the robustness and parallel performance of the SISLICE method on a representative set of SCF eigenvalue problems and outlines research directions that will be explored in future work.

97 MATHEMATICS AND COMPUTING↗

RPP-PRT-58489, Revision 1, One Systems Consistent Safety Analysis Methodologies Report. 24590-WTP-RPT-MGT-15-014

In 2012, One System Nuclear Safety performed a comparison of the safety bases for the Tank Farms Operations Contractor (TOC) and Hanford Tank Waste Treatment and Immobilization Plant (WTP) (RPP-RPT-53222 / 24590-WTP-RPT-MGT-12-018, “One System Report of Comparative Evaluation of Safety Bases for Hanford Waste Treatment and Immobilization Plant Project and Tank Operations Contract”), and identified 25 recommendations that required further evaluation for consensus disposition. This report documents ten NSSC approved consistent methodologies and guides and the results of the additional evaluation process using a new set of evaluation criteria developed for the evaluation of the new methodologies.

61 RADIATION PROTECTION AND DOSIMETRY↗

Confirming a pyrolitic lower mantle using self-consistent pressure scales and new constraints on CaSiO 3 perovskite

In this study, we have examined the lower mantle composition and mineralogy by modeling the density ( ρ ), bulk sound velocity ( V Φ ), and dln ρ /dln V Φ profiles of candidate lower mantle minerals using literature and new experimental equation of state (EoS) results. Additionally, for CaSiO 3 perovskite, complimentary synchrotron X-ray diffraction measurements in a laser-heated diamond anvil cell were conducted up to 156 GPa between 1200 K and 2600 K to provide more reliable constraints on the thermal EoS parameters. These new experimental results as well as literature P-V-T data sets are systematically analyzed using an internally self-consistent pressure scale. We have modeled ρ , V Φ , and dln ρ /dln V Φ profiles of the lower mantle with representative pyrolitic and chondritic compositional models in which the effect of Fe spin transition in ferropericlase is also taken into account. Our modeling results show that a pyrolitic lower mantle with an aggregate mineralogy of 75 vol % bridgmanite, 17 vol % ferropericlase, and 8 vol % CaSiO 3 perovskite produces ρ and V Φ profiles in better agreement with preliminary reference Earth model than a lower mantle with a chondritic composition. The modeled ρ , V Φ , and dln ρ /dln V Φ are mainly affected by the relative ratio of bridgmanite and ferropericlase but are not sensitive to the variation of the CaSiO 3 perovskite content. In addition, the spin crossover of Fe in ferropericlase can greatly raise the value of dln ρ /dln V Φ in the middle lower mantle, which is useful in detecting the presence of ferropericlase in the region. Based on these new mineral physical constraints and radial seismic structure, our study suggests the lower mantle is pyrolitic, which is chemically indistinguishable from the upper mantle.

58 GEOSCIENCES↗

Tundra is a consistent source of CO 2 at a site with progressive permafrost thaw during 6 years of chamber and eddy covariance measurements: Tundra CO 2 Fluxes

Current and future warming of high-latitude ecosystems will play an important role in climate change through feedbacks to the global carbon cycle. This study compares 6 years of CO 2 flux measurements in moist acidic tundra using autochambers and eddy covariance (Tower) approaches. Here, we found that the tundra was an annual source of CO 2 to the atmosphere as indicated by net ecosystem exchange using both methods with a combined mean of 105 ± 17 g CO 2 C m-2 y-1 across methods and years (Tower 87 ± 17 and Autochamber 123 ± 14). Furthermore, the difference between methods was largest early in the observation period, with Autochambers indicated a greater CO 2 source to the atmosphere. This discrepancy diminished through time, and in the final year the Autochambers measured a greater sink strength than tower. Active layer thickness was a significant driver of net ecosystem carbon exchange, gross ecosystem primary productivity, and Reco and could account for differences between Autochamber and Tower. The stronger source initially attributed lower summer season gross primary production (GPP) during the first 3 years, coupled with lower ecosystem respiration (Reco) during the first year. The combined suppression of GPP and Reco in the first year of Autochamber measurements could be the result of the experimental setup. Root damage associated with Autochamber soil collar installation may have lowered the plant community's capacity to fix C, but recovered within 3 years. And while this ecosystem was a consistent CO 2 sink during the summer, CO 2 emissions during the nonsummer months offset summer CO 2 uptake each year.

58 GEOSCIENCES↗

A decade of curtailment studies demonstrates a consistent and effective strategy to reduce bat fatalities at wind turbines in North America

Abstract There is a rapid, global push for wind energy installation. However, large numbers of bats are killed by turbines each year, raising concerns about the impacts of wind energy expansion on bat populations. Preventing turbine blades from spinning at low wind speeds, referred to as curtailment, is a method to reduce bat fatalities, but drawing consistent inference across studies has been challenging. We compiled publicly available studies that evaluated curtailment at six wind energy facilities in North America across 10 years. We used meta‐regression of 29 implemented treatments to determine fatality reduction efficacy as well as sources of variation influencing efficacy. We also estimated species‐specific fatality reduction for three species that comprise most fatalities in North America: hoary bat ( Lasiurus cinereus ), eastern red bat ( Lasiurus borealis ) and silver‐haired bat ( Lasionycteris noctivagans ). We found that curtailment reduced total bat fatalities by 33% with every 1.0 ms −1 increase in curtailment wind speed. Estimates of the efficacy for the three target species were similar (hoary bats: 28% per ms −1 , 95% CI: 0.4%–48%, eastern red bats: 32% per ms −1 , 95% CI: 13%–47% and silver‐haired bats: 32% per ms −1 , 95% CI: 3%–53%). Across multiple facilities and years, a 5.0 ms −1 cut‐in speed was estimated to reduce total bat fatalities by an average of 62% (95% CI: 54%–69%). Mortality reductions at individual facilities in any given year were estimated to fall between 33%–79% (95% prediction interval). Inter‐annual differences rather than inter‐site or turbine characteristics accounted for most of the variation in efficacy rates. Species‐specific average mortality reduction at 5.0 ms −1 curtailment wind speed was 48% (95% CI: 24%–64%) for hoary bats, 61% (95% CI: 42%–74%) for eastern red bats and 52% (95% CI: 30%–66%) for silver‐haired bats. Practical implication . curtailment reduced bat mortality at wind turbines in this North American study. Efficacy increased proportionally as curtailment speed is raised, and patterns and rates of efficacy were similar across species. This indicates that curtailment is an effective strategy to reduce bat fatalities at wind energy facilities, but exploration of further refinements could both minimize bat mortality and maximize energy production.

17 WIND ENERGY↗

High‐Energy LiNiO 2 Li Metal Batteries Enabled by Hybrid Electrolyte Consisting of Ionic Liquid and Weakly Solvating Fluorinated Ether

Abstract In pursuit of the highest possible energy density, researchers shift their focus to the ultimate anode material, lithium metal (Li 0 ), and high‐capacity cathode materials with high nickel content (Ni > 80%). The combination of these aggressive electrodes presents unprecedented challenges to the electrolyte. Here, we report a hybrid electrolyte consisting of a highly fluorinated ionic liquid and a weakly solvating fluorinated ether, whose hybridization structure enables the reversible operation of a battery chemistry based on Li 0 and LiNiO 2 (Ni = 100%), delivering nearly theoretical capacity of the latter (up to 249 mAh g −1 ) for >300 cycles with retention of 78.6% and in absence of unwanted morphological changes in both electrodes. Extensive characterization assisted by molecular dynamic simulation and density functional theory calculations reveals the function of the fluorinated ether to be far more profound than simple dilution and viscosity reduction. Instead, it induces drastic changes in Li + ‐solvation environment, the consequence of which engenders simultaneous stabilization of electrode/electrolyte and interfacing via formation of respective interfacial chemistries. This study further unlocks fundamental knowledge underneath the prevailing “diluent strategy” that is extensively applied by the electrolyte researchers and opens more design space for the next‐generation electrolytes and interphases for these coveted battery chemistries.

25 ENERGY STORAGE↗

Efficient, Thermally Stable, and Mechanically Robust All-Polymer Solar Cells Consisting of the Same Benzodithiophene Unit-Based Polymer Acceptor and Donor with High Molecular Compatibility

All-polymer solar cells (all-PSCs) are a highly attractive class of photovoltaics for wearable and portable electronics due to their excellent morphological and mechanical stabilities. Recently, new types of polymer acceptors (P A s) consisting of non-fullerene small molecule acceptors (NFSMAs) with strong light absorption have been proposed to enhance the power conversion efficiency (PCE) of all-PSCs. However, polymerization of NFSMAs often reduces entropy of mixing in PSC blends and prevents the formation of intermixed blend domains required for efficient charge generation and morphological stability. One approach to increase compatibility in these systems is to design P A s that contain the same building blocks as their polymer donor (P D ) counterparts. Here, a series of NFSMA-based P A s [P(BDT2BOY5-X), (X = H, F, Cl)] are reported, by copolymerizing NFSMA (Y5-2BO) with benzodithiophene (BDT), a common donating unit in high-performance P D s such as PBDB-T. All-PSC blends composed of PBDB-T PD and P(BDT2BOY5-X) P A show enhanced molecular compatibility, resulting in excellent morphological and electronic properties. Specifically, PBDB-T:P(BDT2BOY5-Cl) all-PSC has a PCE of 11.12%, which is significantly higher than previous PBDB-T:Y5-2BO (7.02%) and PBDB-T:P(NDI2OD-T2) (6.00%) PSCs. Additionally, the increased compatibility of these all-PSCs greatly improves their thermal stability and mechanical robustness. For example, the crack onset strain (COS) and toughness of the PBDB-T:P(BDT2BOY5-Cl) blend are 15.9% and 3.24 MJ m -3 , respectively, in comparison to the PBDB-T:Y5-2BO blends at 2.21% and 0.32 MJ m -3 .

14 SOLAR ENERGY↗

Numerical simulation of a thermodynamically consistent four‐species tumor growth model

We develop a thermodynamically consistent four‐species model of tumor growth on the basis of the continuum theory of mixtures. Unique to this model is the incorporation of nutrient within the mixture as opposed to being modeled with an auxiliary reaction‐diffusion equation. The formulation involves systems of highly nonlinear partial differential equations of surface effects through diffuse‐interface models. A mixed finite element spatial discretization is developed and implemented to provide numerical results demonstrating the range of solutions this model can produce. A time‐stepping algorithm is then presented for this system, which is shown to be first order accurate and energy gradient stable. The results of an array of numerical experiments are presented, which demonstrate a wide range of solutions produced by various choices of model parameters.

97 MATHEMATICS AND COMPUTING↗

Self‐Consistent Convolutional Density Functional Approximations: Application to Adsorption at Metal Surfaces

The exchange-correlation (XC) functional in density functional theory is used to approximate multi-electron interactions. A plethora of different functionals are available, but nearly all are based on the hierarchy of inputs commonly referred to as “Jacob's ladder.” This paper introduces an approach to construct XC functionals with inputs from convolutions of arbitrary kernels with the electron density, providing a route to move beyond Jacob's ladder. We derive the variational derivative of these functionals, showing consistency with the generalized gradient approximation (GGA), and provide equations for variational derivatives based on multipole features from convolutional kernels. A proof-of-concept functional, PBEq, which generalizes the PBEα framework with mathematical equation being a spatially-resolved function of the monopole of the electron density, is presented and implemented. It allows a single functional to use different GGAs at different spatial points in a system, while obeying PBE constraints. Analysis of the results underlines the importance of error cancellation and the XC potential in data-driven functional design. After testing on small molecules, bulk metals, and surface catalysts, the results indicate that this approach is a promising route to simultaneously optimize multiple properties of interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Continuum kinetic modelling of cross-separatrix plasma transport in a tokamak edge including self-consistent electric fields

Eulerian kinetic calculations are presented here for the axisymmetric cross-separatrix plasma transport at the edge of a tokamak. The simulations are performed with a high-order finite volume code COGENT, which solves the long-wavelength limit of the full-f ion gyrokinetic equation including the non-linear Fokker–Plank collision model. Self-consistent 2D electrostatic potential variations are obtained from the quasi-neutrality (∇·j = 0) equation coupled to an isothermal fluid electron response. Finally, illustrative solutions presented for parameters characteristic of the DIII-D tokamak show qualitative agreement with experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗