Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Global optimization”

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

Assessing the life-cycle environmental impacts of the wood pallet sector in the United States

Wood pallets play a critical role in the movement and storage of goods worldwide. They are an important component in the complex global supply chain and used by almost every industry. It is therefore important to assess the environmental implications of the wood pallet supply chain and identify optimization strategies that can be implemented. In this paper, primary 2018 annual production data collected from U.S. pallet manufacturers were used to develop the first industry-average life-cycle inventory (LCI). A new functional unit was proposed to perform a more refined and accurate environmental life-cycle assessment of the wood pallet supply chain. Using the LCI data developed, a cradle-to-grave industry-average life-cycle impact assessment was performed. This novel approach quantifies environmental impacts of a generic multi-use pallet, including repair and remanufacturing. The total global warming impact was 10.4 kg CO 2 e per 45.4 t of pallet loads of product delivered using wood pallets. The manufacturing stage contributed the most, about 35%, followed by the raw material supply stage. About 41% of total primary energy consumption was from renewable sources, with most sourced from biomass. Fossil fuels comprised about 52% of the total (225 MJ per functional unit) primary energy consumption. Total environmental impact was significantly affected by two main parameters: reference service life and load-bearing capacity. Pallet repair was also found to be an important component of the wood pallet supply chain, which has a low environmental footprint compared with the overall impact of a pallet and enables mitigation of overall impact by extending the reference service life. Finally, at end-of-life, common industry practices demonstrated substantial potential environmental benefits that can minimize overall environmental impact.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Optimization of Refrigerant Compositions for Low-GWP Refrigerant Mixtures Using Segment-by-segment Heat Exchanger and Detailed System Models

The recently introduced hydrofluoroolefin (HFO) refrigerants, including R1234yf and R1234ze(E), have significantly lower global warming potentials (GWPs) than traditional hydrofluorocarbon (HFC) refrigerants like R410A. However, prior tests show that direct drop-in of pure R1234yf or R1234ze(E) into equipment designed for R410A results in a decrease in heat exchanger capacity and the system coefficient of performance. The primary reason is the lower in-tube heat transfer performance of R1234yf and R1234ze(E) compared with that of R410A. To address this issue, previous studies have mixed the mildly flammable HFC R32 with HFOs to improve system performance, with HFC R125 also added to suppress flammability. Previous studies selected compositions based on simple cycle analyses and did not consider modifications of the heat exchanger circuitry configuration to adapt to the new refrigerants. This study presents a novel multi-objective optimization approach to design a refrigerant composition that maximizes energy efficiency within flammability and GWP limits. The approach in this work simultaneously optimizes mixture composition and heat exchanger circuitry configuration. A case study on a rooftop unit indicates that, compared with mixture-only optimization, simultaneous optimization of mixture and heat exchanger circuitry yields a 5.9% improvement in cycle efficiency and a 48.6% reduction in refrigerant flammability with a GWP of 268. Circuitry optimization using refrigerants with different temperature glides shows that the larger the temperature glide is, the larger EER improvement is obtained. The results show that zeotropic blends with a large temperature glide are more sensitive to the refrigerant circuitry than pure refrigerants and may suffer significant performance degradation with subpar heat exchanger circuitry design. The proposed optimization approach is generally applicable to mixtures with any number of components. Using this approach to design a HVAC system can yield higher system efficiency within flammability and GWP constraints.

Li, Zhenning↗
DOE OSTI2022

Prevalence and scalable control of localized networks

The ability to control network dynamics is essential for ensuring desirable functionality of many technological, biological, and social systems. Such systems often consist of a large number of network elements, and controlling large-scale networks remains challenging because the computation and communication requirements increase prohibitively fast with network size. Here, we introduce a notion of network locality that can be exploited to make the control of networks scalable, even when the dynamics are nonlinear. We show that network locality is captured by an information metric and is almost universally observed across real and model networks. In localized networks, the optimal control actions and system responses are both shown to be necessarily concentrated in small neighborhoods induced by the information metric. This allows us to develop localized algorithms for determining network controllability and optimizing the placement of driver nodes. This also allows us to develop a localized algorithm for designing local feedback controllers that approach the performance of the corresponding best global controllers, while incurring a computational cost orders-of-magnitude lower. Here, we validate the locality, performance, and efficiency of the algorithms in Kuramoto oscillator networks, as well as three large empirical networks: synchronization dynamics in the Eastern US power grid, epidemic spreading mediated by the global air-transportation network, and Alzheimer’s disease dynamics in a human brain network. Taken together, our results establish that large networks can be controlled with computation and communication costs comparable to those for small networks.

42 ENGINEERING↗
DOE OSTI2021

Overcoming Metastable CO 2 Adsorption in a Bulky Diamine-Appended Metal-Organic Framework

We report carbon capture at fossil fuel-fired power plants is a critical strategy to mitigate anthropogenic contributions to global warming, but widespread deployment of this technology is hindered by a lack of energy-efficient materials that can be optimized for CO 2 capture from a specific flue gas. As a result of their tunable, step-shaped CO 2 adsorption profiles, diamine-functionalized metal-organic frameworks (MOFs) of the form diamine-Mg 2 (dobpdc) (dobpdc 4- = 4,4'-dioxidobiphenyl-3,3'-dicarboxylate) are among the most promising materials for carbon capture applications. Here, we present a detailed investigation of dmen-Mg 2 (dobpdc) (dmen = 1,2-diamino-2-methylpropane), one of only two MOFs with an adsorption step near the optimal pressure for CO 2 capture from coal flue gas. While prior characterization suggested that this material only adsorbs CO 2 to half capacity (0.5 CO 2 per diamine) at 1 bar, we show that the half-capacity state is actually a metastable intermediate. Under appropriate conditions, the MOF adsorbs CO 2 to full capacity, but conversion from the half-capacity structure happens on a very slow time scale, rendering it inaccessible in traditional adsorption measurements. Data from solid-state magic angle spinning nuclear magnetic resonance spectroscopy, coupled with van der Waals-corrected density functional theory, indicate that ammonium carbamate chains formed at half capacity and full capacity adopt opposing configurations, and the need to convert between these states likely dictates the sluggish post-half-capacity uptake. By use of the more symmetric parent framework Mg 2 (pc-dobpdc) (pc-dobpdc 4- = 3,3'-dioxidobiphenyl-4,4'-dicarboxylate), the metastable trap can be avoided and the full CO 2 capacity of dmen-Mg 2 (pc-dobpdc) accessed under conditions relevant for carbon capture from coal-fired power plants.

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

Optimizing accuracy and efficacy in data-driven materials discovery for the solar production of hydrogen

The production of hydrogen fuels, via water splitting, is of practical relevance for meeting global energy needs and mitigating the environmental consequences of fossil-fuel-based transportation. Water photoelectrolysis has been proposed as a viable approach for generating hydrogen, provided that stable and inexpensive photocatalysts with conversion efficiencies over 10% can be discovered, synthesized at scale, and successfully deployed. While a number of first-principles studies have focused on the data-driven discovery of photocatalysts, in the absence of systematic experimental validation, the success rate of these predictions may be limited. We address this problem by developing a screening procedure with co-validation between experiment and theory to expedite the synthesis, characterization, and testing of the computationally predicted, most desirable materials. Here, starting with 70 150 compounds in the Materials Project database, the proposed protocol yielded 71 candidate photocatalysts, 11 of which were synthesized as single-phase materials. Experiments confirmed hydrogen generation and favorable band alignment for 6 of the 11 compounds, with the most promising ones belonging to the families of alkali and alkaline-earth indates and orthoplumbates. This study shows the accuracy of a nonempirical, Hubbard-corrected density-functional theory method to predict band gaps and band offsets at a fraction of the computational cost of hybrid functionals, and outlines an effective strategy to identify photocatalysts for solar hydrogen generation.

08 HYDROGEN↗
DOE OSTI2020

Multi-fidelity machine-learning with uncertainty quantification and Bayesian optimization for materials design: Application to ternary random alloys

We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used density functional theory as high-fidelity prediction, while a ML interatomic potential is used as low-fidelity prediction. Practical materials’ design efficiency is demonstrated by reproducing the ternary composition dependence of a quantity of interest (bulk modulus) across the full aluminum–niobium–titanium ternary random alloy composition space. The MFGP is then coupled to a Bayesian optimization procedure, and the computational efficiency of this approach is demonstrated by performing an on-the-fly search for the global optimum of bulk modulus in the ternary composition space. The framework presented in this manuscript is the first application of MFGP to atomistic materials simulations fusing predictions between density functional theory and classical interatomic potential calculations.

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

Clinical practice gaps and challenges in non‐alcoholic steatohepatitis care: An international physician needs assessment

Abstract Background and aims Even as several pharmacological treatments for non‐alcoholic steatohepatitis (NASH) are in development, the incidence of NASH is increasing on an international scale. We aim to assess clinical practice gaps and challenges of hepatologists and endocrinologists when managing patients with NASH in four countries (Germany/Italy/United Kingdom/United States) to inform educational interventions. Methods A sequential mixed‐method design was used: qualitative semi‐structured interviews followed by quantitative online surveys. Participants were hepatologists and endocrinologists practising in one of the targeted countries. Interview data underwent thematic analysis and survey data were analysed with chi‐square and Kruskal‐Wallis tests. Results Most interviewees ( n = 24) and surveyed participants (89% of n = 224) agreed that primary care must be involved in screening for NASH, yet many faced challenges involving and collaborating with them. Endocrinologists reported low knowledge of which blood markers to use when suspecting NASH (56%), when to order an MRI (65%) or ultrasound/FibroScan® (46%), and reported sub‐optimal skills interpreting alanine aminotransferase (ALT, 37%) and aspartate aminotransferase (AST, 38%) blood marker test results, causing difficulty during diagnosis. Participants believed that more evidence is needed for upcoming therapeutic agents; yet, they reported sub‐optimal knowledge of eligibility criteria for clinical trials. Knowledge and skill gaps when managing comorbidities, as well as skill gaps facilitating patient lifestyle changes were reported. Conclusions Educational interventions are needed to address the knowledge and skill gaps identified and to develop strategies to optimize patient care, which include implementing relevant care pathways, encouraging referrals and testing, and multidisciplinary collaboration, as suggested by the recent Global Consensus statement on NAFLD.

Lazure, Patrice↗
DOE OSTI2023

Efficient QAOA Optimization using Directed Restarts and Graph Lookup

Variational Quantum Algorithms (VQA) aim to enhance the capabilities of Noisy Intermediate-Scale Quantum (NISQ) devices. These algorithms utilize parameterized circuits and classical optimizers to iteratively execute circuits with varying parameters. However, VQA faces computational overheads due to repeated iterations and random restarts. Prior work suggests using basic sub-graphs to transfer parameters for the input graph, reducing optimizer overheads but limiting applicability to structured regular graphs. In real-world applications, random irregular graphs are common, and existing methods are not scalable or practical for such graphs. This paper presents a framework that aims to improve random irregular graphs in VQA. The framework uses graph similarity and important features like total edge counts, average edge counts, and variance. It follows an iterative process to choose basis sub-graphs from a small database and adjust parameters accordingly. Classical optimizers then utilize these parameters to determine when to restart and perform gradient descent. This approach increases the chances of reaching global maximum points.

Wang, Meng↗
DOE OSTI2021

AI-Batt (Autonomous Identification of Battery Life Models) [SWR 21-36]

Autonomous Identification of Battery Life Models (AI-Batt) AI-Batt is a MATLAB code base for developing lifetime models for batteries from accelerated aging data. The code base provides many functions for processing, visualizing, and modeling battery aging data, making the data processing, exploration, and modeling workflow substantially faster. These tools are tailored for working with battery aging data sets, which usually consist of many separate time-series for each cell, with many test conditions and possible replicates at each condition, which makes it difficult to simply process or visualize the data set. Complex modeling tasks, such as cross-validation, sensitivity analysis, and uncertainty quantification have been implemented to enable thorough statistical investigation of model predictions. Additionally, several machine-learning algorithms are implemented to autonomously identify suitable models via symbolic regression. Data processing functions automatically cast data from the struct data type, which is commonly used to store experimental data, but is not an acceptable input for most algorithms, to the table data type, which can be easily used as input to any optimization algorithm. Also, the data can be separated into time-invariant and time-variant data tables, which is helpful for exploring the data set as well as developing separate models for time-variant and time-invariant aging mechanisms. For example, in aging tests with constant temperature, temperature is a time-invariant experimental condition. Visualization tools enable plotting of data, model fits, and model simulations possible with single-line function calls, empowering data exploration of complex data sets with both time-varying and time-invariant trends. Plots can be automatically generated for the whole data set, or separated by data group (groups of test replicates) or individual data series. Data points or data series can be automatically colored by the value of a variable with a variety of color maps, and model predictions can also be colored by the value of a fit statistic. Comparisons between data sets and the predictions/simulations of different models on the same data set can be easily plotted as well. Distributions of parameter values from bootstrap resampling can be plotted to visualize the reliability of parameter estimation, or determine any correlations between parameters. Modeling tools handle the complex task of creating and parsing symbolic equations for modeling battery lifetime. Equations are parsed to grab relevant data variables, parameter values, or specified sub-models for input into optimization, evaluation, or simulation functions. Models can be optimized locally (one set of parameters for each data series), bi-level (some parameters shared across the data set), or globally (single set of parameters for all data). Functions implementing symbolic regression algorithms help users to discover effective model equations, even in poorly sampled, high-dimensional data.

Smith, Kandler↗
DOE OSTI2023

Tropical Red Macroalgae Cultivation with a Focus on Compositional Analysis

To create carbon efficient sources of bioenergy feedstocks and feedstuff for aquaculture and terrestrial livestock, it is critical to develop and commercialize the most efficient seaweed cultivation approach with a sustainable nutrient input supply. Here, we present data for a novel, onshore tropical macroalgae cultivation system, based on influent deep seawater as the nutrient and carbon sources. Two red algal species were selected, Agardhiella subulata and Halymenia hawaiiana, as the basis for growth optimization. Highest productivity in small-scale cultivation was demonstrated with A. subulata in the 10% deep seawater (64.7 µg N L−1) treatment, growing at up to 26% specific growth rate day−1 with highest yields observed at 247.5 g m−2 day−1 fresh weight. The highest yields for H. hawaiiana were measured with the addition of 10% deep seawater up to 8.8% specific growth rate day−1 and yields at 63.3 g fresh weight m−2 day−1 equivalent. Biomass should be culled weekly or biweekly to avoid density limitations, which likely contributed to a decrease in SGR over time. With a measured 30–40% carbon content of the ash-free dry weight (20–30% of the dry weight) biomass, this translates to an almost 1:1 CO2 capture to biomass ratio. The compositional fingerprint of the high carbohydrate content of both Agardhiella and Halymenia makes for an attractive feedstock for downstream biorefinery applications. By focusing on scaling and optimizing seaweed farming technologies for large-scale onshore farms, the opportunities for yield potential, adaptability to cultivation conditions, and meeting global sustainability goals through novel, carbon-negative biomass sources such as seaweed can be realized.

09 BIOMASS FUELS↗
DOE OSTI2022

Calibration of cloud and aerosol related parameters for solar irradiance forecasts in WRF-solar

Model parameters are a major source of uncertainty in numerical weather prediction. Recently, the Weather Research and Forecasting model with Solar extensions (WRF-Solar) has been upgraded by enhancing the treatment of sub-grid scale cloud and aerosols with augmentations of a sub-grid scale cloud scheme (CLD3) and an upgraded aerosol-aware Thompson-Eidhammer scheme (TE14). However, the value of model parameters associated with these parameterizations are assigned based on limited measurements or theoretical calculations. Calibrating the most sensitive parameters has the potential to improve solar irradiance predictions. Here, we adopted a multiobjective surrogate-based optimization (SBO) framework to calibrate nine parameters used in CLD3 and TE14 that lead to the largest sensitivity in simulated irradiance. The normalized mean-absolute-error (NMAE) of global horizontal irradiance (GHI) and direct normal irradiance (DNI) are minimized by calibrating WRF-Solar over two regions including the Southern Great Plains (SGP) and Central California, in order to focus on parameter calibration under cloudy conditions with different aerosol loading. The results show that generalized linear model (GLM)-based surrogate models approximate physical models well, particularly when the third order and three-way interaction terms are considered. The SBO framework efficiently searches the parameter space for optimal solutions with less computational costs than directly calibrating the physical model. We first calibrate CLD3 parameters over the less-polluted SGP region. Optimized CLD3 parameters alone result in NMAE reduction by 14% for the site-mean and up to 33% for individual cases over the SGP region. With further calibration of TE14 parameters over the Central California during active fire periods, the optimized parameters lead to over 20% reductions of NMAE. Our investigation reveals, however, that optimizing TE14 has a limited impact on irradiance simulations under less-polluted conditions in the SGP.

14 SOLAR ENERGY↗
DOE OSTI2020

Sensitivity Studies, Gap Analysis, and Benchmark Experiment Optimization for Reactor Physics and Criticality Safety Applications

Many new reactor designs, such as advanced reactors and micro reactors, have materials that lack nuclear data validation. This is also true for many other applications in criticality safety and global security. Both differential and integral experiments are needed to validate cross-section data. Without this, a user cannot have confidence in the predicted results of a radiation-transport code. This work describes an approach called ARCHIMEDES (Application Relevant Critical/Subcritical HEU/Pu-based Integral Measurements for Enhancing Data and Evaluating Sensitivities) to design new criticality experiments that have similar k eff cross-section sensitivities to an application of interest. This process involves simulations to generate cross-section sensitivities to a parameter of interest (such as k eff ), a gap analysis to determine which existing benchmarks are most similar to the application, and an experiment optimization. Recently, there has been a great deal of interest in the reactor physics community on advanced reactors, micro reactors, and accelerator driven systems (ADS). This work will apply the described method to specific examples in this area. The focus of this work will be on the sensitivity study and gap analysis, while future work will include experiment design.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Nitrogen and phosphorous acquisition strategies drive coexistence patterns among archaeal lineages in soil

Soil represents the largest reservoir of Archaea on Earth. Present-day archaeal diversity in soils globally is dominated by members of the class Nitrososphaeria. The evolutionary radiation of this class is thought to reflect adaptations to a wide range of temperatures, pH, and other environmental conditions. However, the mechanisms that govern competition and coexistence among Nitrososphaeria lineages in soil remain poorly understood. Here we show that predominant soil Nitrososphaeria lineages compose a patchwork of gene inventory and expression profiles for ammonia, urea, and phosphate utilization. In contrast, carbon fixation, respiration, and ATP synthesis genes are conserved and expressed consistently among predominant phylotypes across 12 major evolutionary lineages commonly found in soil. In situ gene expression profiles closely resemble pure culture reference strains under optimal growth conditions. Together, these results reveal resource-based coexistence patterns among Nitrososphaeria lineages and suggest complementary ecophysiological niches associated with differential nutrient acquisition strategies among globally predominant archaeal lineages in soil.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Forest residue harvest optimization: spanning the bridge between plant biology and biorefinery performance

Forestry residues have immense potential as alternative feedstocks to petroleum, yet their inherent complexity remains a major challenge to widespread use. Pairing the temporal rhythms of plant biology with biorefinery performance is critical to industrial-scale biorefinery development. Here, we provide the first report of a techno-economic analysis (TEA) and life cycle assessment (LCA) for a model integrated reductive catalytic fractionation (RCF)–molten salt hydrolysis process for forestry residues varying in tree part, species, and phenophase. All forestry residues resulted in net-negative greenhouse gas (GHG) emissions vs. comparable petroleum feedstocks, with GHG emissions potentially reduced >4.0× through composition-based feedstock selection (e.g., harvesting American beech bark in spring vs. summer). Moreover, American beech twigs/branchlets and bark in leafed and emergence phenophases, respectively, had 7.9× lower predicted phenolic minimum selling prices (MSPs) vs. other feedstocks and MSPs within the current global phenolic market range. Hemicellulose content and RCF yield emerged as key parameters impacting GHG emissions and biorefinery revenue, identifying hardwood twigs/branchlets in the leafed phenophase as optimal biofeedstocks. Biorefinery expenses were dominated by purchased equipment, raw materials, and utility costs, highlighting essential areas for future study. Notably, RCF reactor pressures drove 85–90% of equipment costs, but sensitivity analysis revealed that decreasing the pressure 20% could reduce the phenol MSP 4-fold. Structural carbohydrate dynamics were also investigated using a two-step acid hydrolysis method to resolve tissue- and species-level patterns in biomass composition throughout the year to enable harvest optimization based on TEA/LCA findings. Ultimately, elucidating the impact of biofeedstock dynamics on biorefinery performance enables harvest optimization, informed engineering design, and progress towards an expanded bioeconomy.

Shapiro, Alison J. [University of Delaware, Newark↗
DOE OSTI2021

Levenberg–Marquardt multi-classification using hinge loss function

Incorporating higher-order optimization functions, such as Levenberg-Marquardt (LM) have revealed better generalizable solutions for deep learning problems. However, these higher-order optimization functions suffer from very large processing time and training complexity especially as training datasets become large, such as in multi-view classification problems, where finding global optima is a very costly problem. To solve this issue, we develop a solution for LM-enabled classification with, to the best of knowledge first-time implementation of hinge loss, for multiview classification. Hinge loss allows the neural network to converge faster and perform better than other loss functions such as logistic or square loss rates. Here we prove our method by experimenting with various multiclass classification challenges of varying complexity and training data size. The empirical results show the training time and accuracy rates achieved, highlighting how our method outperforms in all cases, especially when training time is limited. Our paper presents important results in the relationship between optimization and loss functions and how these can impact deep learning problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Piston Bowl Geometry Effects on Gasoline Compression Ignition in a Heavy-Duty Diesel Engine

A design optimization campaign was conducted to search for improved combustion profiles that enhance gasoline compression ignition in a heavy-duty diesel engine with a geometric compression ratio of 17.3. A large-scale design of experiments approach was used for the optimization, employing three-dimensional computational fluid dynamics simulations. The main parameters explored include geometric features, injector specifications, and swirl motion. Both stepped-lip and re-entrant bowls were included in order to assess their respective performance implications. A total of 256 design candidates were prepared using the software package CAESES for automated and simultaneous geometry generation and combustion recipe perturbation. The design optimization was conducted for three engine loads representing light to medium load conditions. The design candidates were evaluated for fuel efficiency, emissions, fuel-air mixing, and global combustion behavior. Simulation results showed that the optimum designs were all stepped-lip bowls, due to improvements in fuel-air mixing, as well as reduced heat loss and emissions formation. Improvements in indicated specific fuel consumption of up to 3.2% were achieved while meeting engine-out NOx emission targets of 1-1.5 g/kW center dot h. Re-entrant bowls performed worse compared to the baseline design, and significant performance variations occurred across the load points. Specifically, the re-entrant bowls were on par with the stepped-lip bowls under light load conditions, but significant deteriorations occurred under higher load conditions. As a final task, selected optimized designs were then evaluated under full-load conditions.

CFD↗
DOE OSTI2020

Matrix-free subcell residual distribution for Bernstein finite elements: Monolithic limiting

This paper is focused on the aspects of limiting in residual distribution (RD) schemes for high-order finite element approximations to advection problems. Both continuous and discontinuous Galerkin methods are considered in this work. Discrete maximum principles are enforced using algebraic manipulations of element contributions to the global nonlinear system. The required modifications can be carried out without calculating the element matrices and assembling their global counterparts. The components of element vectors associated with the standard Galerkin discretization are manipulated directly using localized subcell weights to achieve optimal accuracy. Low-order nonlinear RD schemes of this kind were originally developed to calculate local extremum diminishing predictors for flux-corrected transport (FCT) algorithms. In the present paper, we incorporate limiters directly into the residual distribution procedure, which makes it applicable to stationary problems and leads to well-posed nonlinear discrete problems. To circumvent the second-order accuracy barrier, the correction factors of monolithic limiting approaches and FCT schemes are adjusted using smoothness sensors based on second derivatives. Finally, the convergence behavior of presented methods is illustrated by numerical studies for two-dimensional test problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Superconductivity, charge density waves, and bipolarons in the Holstein model

The electron-phonon ( e-ph ) interaction remains of great interest in condensed matter physics and plays a vital role in realizing superconductors, charge density waves (CDW), and polarons. We study the two-dimensional Holstein model for e-ph coupling using determinant quantum Monte Carlo across a wide range of its phase diagram as a function of temperature, electron density, dimensionless e-ph coupling strength, and the adiabatic ratio of the phonon frequency to the Fermi energy. We describe the behavior of the CDW correlations, the competition between superconducting and CDW orders and polaron formation, the optimal conditions for superconductivity, and the transition from the weak-coupling regime to the strong-coupling regime. Superconductivity is optimized at intermediate e-ph coupling strength and intermediate electron density, and the superconducting correlations increase monotonically with phonon frequency. The global maximum for superconductivity in the Holstein model occurs at large phonon frequency, the limit where an attractive Hubbard model effectively describes the physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗