Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Operator approximation”

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

Effects of finite element discretization and model simplification on calculations of ductile failure initiation

The finite element method is a scheme to discretize the infinite number of degrees of freedom in continuum-level problems down to a finite number of degrees of freedom. This discretization is done in conjunction with methods that also reduce the field differential equations to sets of algebraic ones that can be solved by arithmetical operations. Therefore, solutions attained by finite element models are approximations to the exact solutions of the field equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Jet Fuel Production at the Pittsburgh Airport: GTL via Fischer-Tropsch Synthesis

The Pittsburgh International Airport (PIT)—with the Allegheny County Airport Authority (which manages PIT)—has established itself as a leader in resiliency by becoming the first major United States (U.S.) airport to have a self-sustaining microgrid, providing electricity, heating, and cooling for airport operations. The microgrid is powered by natural gas and solar power produced on the airport property and was completed in Summer 2021. This study examines the feasibility of producing jet fuel at the airport to provide a secure supply of aviation fuel, furthering PIT’s ability to weather supply disruptions and operate self-sufficiently. Gas-to-liquids (GTL) is a commercially available technology that converts natural gas to liquid hydrocarbons, including synthetic jet fuel. A GTL facility at PIT could convert natural gas from onsite wells to jet fuel, effectively doubling the onsite fuel stores in the event of a supply disruption. Moreover, GTL provides a pathway to renewable jet fuel production and reduced greenhouse gas (GHG) emissions from the aviation sector, particularly if renewable natural gas (RNG) is used as a feedstock or other renewable energy sources are used for energy inputs. This study has found that it would be technically feasible to construct and operate a GTL facility on PIT’s property. The approximately 6,000-barrel per day (BPD) facility evaluated would produce nearly 70 million (MM) gallons (gal) of synthetic jet fuel per year, which could supplant nearly all (85 percent) current jet fuel consumption at PIT. Given the current blend limitation of 50 percent Fischer-Tropsch fuels by volume, the plant would have excess production capacity available for the United States Air Force (USAF) Pittsburgh Air Reserve Station and the USAF 171st Air Refueling Wing co-located at the airport.

03 NATURAL GAS↗
DOE OSTI2022

Integration of Pumped Heat Energy Storage with Fossil-Fired Power Plant (Final Report)

The project team of Southwest Research Institute ® (SwRI ® ), Malta Inc. (Malta), and Luminant Generation Company LLC (Luminant) completed a feasibility study for the integration of a 100-MW, 10-hour (1000-MWh) Malta Pumped Heat Energy Storage (MPHES) system with multiple full-sized fossil-fired electricity generation units (EGU) in Luminant portfolio. MPHES is a long-duration, molten-salt-energy storage technology that uses turbomachinery and heat exchangers to transfer energy to a thermal storage media when charging, and removes the heat in a similar fashion when discharging. With high round trip efficiency (60-65%) and long lifespan (30+ years), MPHES provides economic benefits to the fossil-asset owners that can be scaled to integrate with assets across their portfolio. This technology uses hardware components, workforce personnel, and skillsets similar to those used by fossil EGUs, allowing for synergy when co-locating the two technologies. Luminant has approximately 39,000 megawatts of generation across 12 states, operating in six of the seven competitive markets in the U.S. and powered by a diverse portfolio of natural gas, nuclear, coal, and solar facilities. The DeCordova plant in Granbury, Texas, a simple cycle natural gas peaker power plant, was used as the fossil-fired asset in this project. The local market in Granbury, Texas has many influences, including several nearby power plants, a Luminant-owned nuclear plant (Comanche Peak), and substantial wind energy, which causes both negative pricing at night and high market volatility. Reducing false starts of the DeCordova plant and better responding to market volatility would be an economic advantage. Luminant is currently integrating battery storage plant on site to begin addressing these challenges. Integrating long-duration storage, like MPHES, would expand this capability beyond one hour of storage and have the potential to greatly reduce the total number of gas turbine starts. The MPHES charging requirement could help offset the overnight operating costs of Comanche Peak, which cannot load follow, and the nearby Luminant-owned Wise County combined cycle plant that cycles too often. Following the assessment of Luminant’s ERCOT-based portfolio for integration compatibility with Malta’s PHES system and the project tasks of conceptual study, technoeconomic analysis, technology gap assessment, and commercialization plan, the project team effort has resulted in several key outcomes: (1) Identification of market trends in a high-wind penetration market outside a major metropolitan area; (2) Creation of a dispatching model for the MPHES system and the pairing of Li-ion battery with a gas turbine in a real time market; (3) Revenue and cost estimations for operating MPHES alongside a gas peaker plant with real dispatching considerations and comparison with variations in the Malta implementation, including doubling the storage capacity and using two discharge drivetrains; (4) Potential carbon emission reductions possible by replacing gas turbine operation with Malta PHES operation; and (5) Summary of literature-based future market predictions for Texas.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Evaluation of the Radioactive Material Release in the Harborview Research and Training Building and Implications for Emergency Response

On May 2, 2019, during the 137 Cs source recovery operation, a source capsule in a research irradiator containing approximately 77.1 TBq was breached. Based on a geometric reconstruction analysis of the damage to the capsule, approximately 46.3 GBq (0.04%) was impacted by the chop saw (grinder) inside a Mobile Hot Cell (MHC) on the loading dock at the University of Washington Harborview Research and Training (HRT) Building. A very small fraction of the material impacted, less than 1%) was released from the Mobile Hot Cell and then to the rest of the HRT Building. The objectives of this project were to assess the accidental release of 137 CsCl and its implications related to emergency response methods and the ramifications of 137 CsCl transport. The phenomenology of this event was also compared with past alkali halide dispersal events. The vast number of measurements and samples collected by the remediation contractors, the Department of Energy's Nuclear Emergency Support Team, and the small number of retrospective samples collected by the authors informed the analysis. The techniques included (1) autoradiography and electron microscopy of samples collected from the HRT Building and the irradiator, (2) 3D visualization of deposition on surfaces and within the ventilation system, and (3) a study of the damage to the source capsule to evaluate the Cs particle size and particle composition due to the grinding accident. Subsequently, the cesium contaminant transport through the numerous pathways in the building was reconstructed to assess the deposition on surfaces as a function of particle size. The implications for emergency response are relevant to data quality and management. A Data Quality Objective (DQO) guides data collection methods so that they have appropriate accuracy and precision for the intended application. Recommendations were made with respect to the sample collection protocols and sample archival.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Magnesium production by molten salt electrolysis with liquid tin cathode and multiple effect distillation

Low-cost clean primary production of magnesium metal is important for its use in many applications, from light-weight structural components to energy technologies. This work describes new experiments and cost and emissions analysis for a magnesium metal production process. The process combines molten salt electrolysis of MgO using MgF₂-CaF₂ electrolyte and a reactive liquid tin cathode, with gravity-driven multiple effect thermal system (G-METS) distillation to separate out the magnesium product, and re-use of the tin. Electrolysis experiments with carbon anodes showed current yield above 90%, while a yttria-stabilized zirconia solid oxide membrane (SOM) anode experiment showed 84% current yield. G-METS distillation is an important component of the envisioned process. It can potentially lower costs and energy use considerably compared with conventional magnesium distillation. Techno-economic analysis including detailed mass and energy balances shows that this electrolyte composition could lower costs by utilizing CaO, which is the primary impurity in MgO, as the Hall-Héroult process uses the sodium impurity in alumina. Analysis options include: raw material types (magnesite rock vs. brine or seawater), drying and calcining using electricity vs. natural gas, and carbon vs. SOM anode type. Using SOM inert anodes results in a cost premium around 10%–15%, mostly due to higher electrical energy usage resulting from membrane resistance, and reduces GHG emissions by approximately 1 kg CO₂/kg Mg product. Capital and operating cost estimates, and cradle to gate greenhouse gas (GHG) emissions analysis under several raw material and process technology scenarios, show comparable costs and emissions to those of aluminum production.

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

Efficient Implementation of Artificial Neural Networks for Sensor Data Analysis Based on a Genetic Algorithm

The reliability of many industrial processes depends on the sensor system. However, these sensors can be affected by noise, perturbations and failures. Hence, sensor monitoring and diagnosis are fundamental to guarantee the quality of an industrial process. Nowadays, artificial neural networks (ANN) are widely used in sensor signal processing and diagnosis. However, those ANNs usually require many artificial neurons, being difficult to implement in software and hardware due to their high computational costs. This paper presents an optimized implementation of artificial neurons in ANNs for sensor data analysis using a Genetic Algorithm (GA). The objective of GA is to find an adequate segmentation to reduce the activation function approximation error. One of the advantages of the proposed approach is that the cost function used in GA considers the effect of factors such as the ANN architecture or the number of bits used in arithmetic operations. The proposed ANN implementation technique aims to get the best possible approximation for a specific ANN architecture, making easier its implementation in software and hardware. Simulation and experimental results using FPGA (Field Programmable Gate Array) prove the advantages of the proposed approach for implementing sensor data analysis systems based on ANNs.

D estefani, André↗
DOE OSTI2023

Godiva IV Central Cavity Neutron Environment Characterization with Threshold Neutron Detectors

Godiva IV is a cylindrical fast burst reactor comprised of approximately 65 kg of highly enriched uranium that is operated by Los Alamos National Laboratory and sited at the National Criticality Experiments Research Center at the Nevada National Security Site in Nevada in the United States. Godiva IV is typically operated at delayed critical and in the regime spanning from sub-prompt to super-prompt bursts. Godiva IV is used for sample irradiations, criticality safety demonstrations, dosimetry studies, and for studying super-prompt behavior. In preparation for both an upcoming experiment to reduce uncertainties in the prompt fission spectrum for 235 U using threshold neutron detectors, and for future research using Godiva IV, it was desired to exercise the process of the selection of threshold neutron detectors/activation foils, radiation metrology, and the subsequent adjustment of the neutron spectrum. For this exercise, nine high purity threshold neutron detectors/activation foils were irradiated in a Godiva IV burst. The foils were then analyzed using a high purity germanium detector in the NCERC counting laboratory to determine end of irradiation specific activities for available IRDFF-II reactions. This work summarizes the Godiva IV foil irradiation, radiation metrology results, and adjusted neutron spectrum. The results of this exercise ultimately characterized the neutron environment inside the sample irradiation cavity inside Godiva IV to a higher degree than previously performed, informed decisions for the upcoming larger scale experiment, and will inform future neutron spectrum characterizations at NCERC.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Network-Level Optimization for Unbalanced Power Distribution System: Approximation and Relaxation

The nonlinear programming (NLP) problem to solve distribution-level optimal power flow (D-OPF) poses convergence issues and does not scale well for unbalanced distribution systems. The existing scalable D-OPF algorithms either use approximations that are not valid for an unbalanced power distribution system, or apply relaxation techniques to the nonlinear power flow equations that do not guarantee a feasible power flow solution. In this paper, we propose scalable D-OPF algorithms that simultaneously achieve optimal and feasible solutions by solving multiple iterations of approximate, or relaxed, D-OPF subproblems of low complexity. The first algorithm is based on a successive linear approximation of the nonlinear power flow equations around the current operating point, where the D-OPF solution is obtained by solving multiple iterations of a linear programming (LP) problem. The second algorithm is based on the relaxation of the nonlinear power flow equations as conic constraints together with directional constraints, which achieves optimal and feasible solutions over multiple iterations of a second-order cone programming (SOCP) problem. Finally, it is demonstrated that the proposed algorithms are able to reach an optimal and feasible solution while significantly reducing the computation time as compared to an equivalent NLPD-OPF model for the same distribution system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Inverter Reliability Estimation for Advanced Inverter Functionality

In the near future, grid operators are expected to regularly use advanced distributed energy resource (DER) functions, defined in IEEE 1547-2018, to perform a range of grid-support operations. Many of these functions adjust the active and reactive power of the device through commanded or autonomous modes, which will produce new stresses on the grid-interfacing power electronics components, such as DC/AC inverters. In previous work, multiple DER devices were instrumented to evaluate additional component stress under multiple reactive power setpoints. We utilize quasi-static time-series simulations to determine voltage-reactive power mode (volt-var) mission profile of inverters in an active power system. Mission profiles and loss estimates are then combined to estimate the reduction of the useful life of inverters from different reactive power profiles. It was found that the average lifetime reduction was approximately 0.15% for an inverter between standard unity power factor operation and the IEEE 1547 default volt-var curve based on thermal damage due to switching in the power transistors. For an inverter with an expected 20-year lifetime, the 1547 volt-var curve would reduce the expected life of the device by 12 days. This framework for determining an inverter's useful life from experimental and modeling data can be applied to any failure mechanism and advanced inverter operation.

component degradation↗
DOE OSTI2025

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

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

Understanding carbon sourcing and transport originating from the helicon antenna surfaces during high-power helicon discharge in DIII-D Tokamak

The high-power helicon wave system in the DIII-D tokamak could potentially introduce new plasma–material interaction (PMI) challenges owing to rectified RF sheath potentials that develop near the antenna and surrounding plasma-facing components. We present the first application of the STRIPE (Simulated Transport of RF Impurity Production and Emission) framework to helicon-induced PMIs, extending previous STRIPE studies of ICRH antennas by incorporating net erosion, local re-deposition, and three-dimensional global impurity transport. The integrated workflow couples SOLPS, COMSOL, RustBCA, GITR, and GITRm to simulate carbon erosion, re-deposition, and impurity transport for two experimentally constrained DIII-D H-mode helicon operating scenarios with different antenna–plasma gaps, coupled RF powers, and edge plasma conditions. COMSOL predicts rectified RF sheath potentials of 1–5 kV localized near the lower portion of the antenna, where the magnetic field intersects the surface at grazing incidence. Carbon self-sputtering dominates the erosion source, whereas RF-accelerated D + ions contribute approximately 1% of the total gross erosion. The smaller-gap operating scenario exhibits substantially stronger gross erosion, enhanced local re-deposition (∼12%), and a larger confined carbon inventory owing to increased plasma accessibility and broader RF sheath coverage. Comparison with available DIII-D measurements shows no distinct change in the global carbon signal that correlates with the helicon RF pulse, consistent with the simulations indicating that the helicon-generated carbon source remains small compared with the existing background carbon inventory under the present graphite-wall operating conditions. These results demonstrate the capability of STRIPE to integrate RF sheath modeling, plasma transport, surface interaction physics, and three-dimensional impurity transport for the interpretation of helicon-induced PMIs. The study further identifies the principal sources of modeling uncertainty, including grazing-angle RF sheath physics, slow-wave resolution, plasma-background extrapolation, and trace-impurity transport assumptions, providing a framework for future validation and model development.

Kumar, Atul [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2025

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗
DOE OSTI2021

Proactive Operations and Investment Planning via Stochastic Optimization to Enhance Power Systems’ Extreme Weather Resilience

We present scalable stochastic optimization approaches for improving power systems’ resilience to extreme weather events. We consider both proactive redispatch and transmission line hardening as alternatives for mitigating expected load shed due to extreme weather, resulting in large-scale stochastic linear programs (LPs) and mixed-integer linear programs (MILPs). We solve these stochastic optimization problems with progressive hedging (PH), a parallel, scenario-based decomposition algorithm. Our computational experiments indicate that our proposed method for enhancing power system resilience can provide high-quality solutions efficiently. With up to 128 scenarios on a 2,000-bus network, the operations (redispatch) and investment (hardening) resilience problems can be solved in approximately 6 min and 2 h of wall-clock time, respectively. Additionally, we solve the investment problems with up to 512 scenarios, demonstrating that the approach scales very well with the number of scenarios. Moreover, the method produces high quality solutions that result in statistically significant reductions in expected load shed. Our proposed approach can be augmented to incorporate a variety of other operational and investment resilience strategies, or a combination of such strategies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

INTEGRATION OF FLEX EQUIPMENT AND OPERATOR ACTIONS IN PLANT FORCE-ON-FORCE MODELS WITH DYNAMIC RISK ASSESSMENT

The overall operation and maintenance cost to protect nuclear power plants accounts for approximately 7% of the total cost of power generation, with labor accounting for half of this cost. In the current research, from interaction with utilities and other stakeholders, it was determined that physical security forces account for nearly 20% of the entire workforce at several nuclear power plants. Labor costs continue to rise in the U.S., so any measures to reduce the cost of operating a nuclear power plant will need to include a reduction in labor. The physical security pathway within the DOE’s Light Water Reactor Sustainability program aims to lower the cost of physical security through directed research into modeling and simulation, application of advanced sensors or deployment of advanced weapons. This report presents a modeling and simulation framework for integrating Diverse and Flexible Mitigation Capability (FLEX) portable equipment performance with Force on Force models of a plant’s physical security posture. The generic framework is described in detail, followed by a case study of modeling an adversarial attack aimed at causing a radiological release by sabotaging the plant’s power supply and its ultimate heat sink capabilities at a hypothetical nuclear power plant. Two different FLEX deployment strategies, series and parallel, are modeled with distinct timelines. The results of the adversarial attack modeled in a commercial Force on Force tool are integrated with the FLEX deployment model in INL’s dynamic modeling tool EMRALD. Monte Carlo simulation is used to model the distribution of the timeline in FLEX deployment strategies. The results demonstrate that, even in the extreme case of a successful adversarial attack, deployment of FLEX equipment can result in a significantly high likelihood of preventing radiological release. The modeling and simulation framework integrating FLEX equipment with Force on Force models enables the nuclear power plants to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Quantifying Uncertainty of Deep Reinforcement Learning Based Decision Making for Operations and Maintenance of Nuclear Power Plant

This paper summarizes research that integrates condition monitoring and prognostics with decision making for nuclear power plant operations and maintenance. As part of this research, we have developed an online asset management tool to help reduce life-cycle maintenance and repair costs. Using the latest advancements in condition monitoring, supply chain analytics, and deep reinforcement learning, we have created a predictive maintenance tool that can optimize the maintenance and spare-part management of a repairable nuclear system. To demonstrate these methods, preliminary studies were conducted on a simple, representative maintenance system undergoing a stochastic degradation process that requires repairs or replacement to continue operation. Through Monte Carlo simulations, we were able to reduce maintenance spending by approximately 50% compared to optimized, time-based maintenance strategies. Not only does the decision maker reduce the average life-cycle costs, it also minimizes the chance of high cost scenarios, lowering the variance of the expected cost distributions, and reducing overall financial risk. Furthermore, this work also studies the ability of the decision maker to handle various levels of noise from observation uncertainty. By introducing uncertainty into the decision-making process, we have quantified the robustness and resiliency of the decision maker, as well as identified necessary levels of observability to demonstrate cost effectiveness.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Development of Real-Time System Identification to Detect Abnormal Operations in a Gas Turbine Cycle

Here, we present a novel online system identification methodology for monitoring the performance of power systems. This methodology was demonstrated in a gas turbine recuperated power plant designed for a hybrid configuration. A 120-kW Garrett microturbine modified to test dynamic control strategies for hybrid power systems designed at the National Energy Technology Laboratory (NETL) was used to implement and validate this online system identification methodology. The main component of this methodology consists of an empirical transfer function model implemented in parallel to the turbine speed operation and the fuel control valve, which can monitor the process response of the gas turbine system while it is operating. During fully closed-loop operations or automated control, the output of the controller, fuel valve position, and the turbine speed measurements were fed for a given period of time to a recursive algorithm that determined the transfer function parameters during the nominal condition. After the new parameters were calculated, they were fed into the transfer function model for online prediction. The turbine speed measurement was compared against the transfer function prediction, and a control logic was implemented to capture when the system operated at nominal or abnormal conditions. To validate the ability to detect abnormal conditions during dynamic operations, drifting in the performance of the gas turbine system was evaluated. A leak in the turbomachinery working fluid was emulated by bleeding 10% of the airflow from the compressor discharge to the atmosphere, and electrical load steps were performed before and after the leak. This tool could detect the leak 7 s after it had occurred, which accounted for a fuel flow increase of approximately 15.8% to maintain the same load and constant turbine speed operations.

algorithms↗
DOE OSTI2023

Exploiting Kronecker structure in exponential integrators: Fast approximation of the action of φ $-$functions of matrices via quadrature

Here, in this article, we propose an algorithm for approximating the action of φ $-$ functions of matrices against vectors, which is a key operation in exponential time integrators. In particular, we consider matrices with Kronecker sum structure, which arise from problems admitting a tensor product representation. The method is based on quadrature approximations of the integral form of the φ $-$ functions combined with a scaling and modified squaring method. Owing to the Kronecker sum representation, only actions of 1D matrix exponentials are needed at each quadrature node and assembly of the full matrix can be avoided. Additionally, we derive a priori bounds for the quadrature error, which show that, as expected by classical theory, the rate of convergence of our method is supergeometric. Guided by our analysis, we construct a fast and robust method for estimating the optimal scaling factor and number of quadrature nodes that minimizes the total cost for a prescribed error tolerance. We investigate the performance of our algorithm by solving several linear and semilinear time-dependent problems in 2D and 3D. The results show that our method is accurate and orders of magnitude faster than the current state-of-the-art.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Compositional texture engineering for highly stable wide-bandgap perovskite solar cells

The development of highly stable and efficient wide-bandgap (WBG) perovskite solar cells (PSCs) based on bromine-iodine (Br–I) mixed-halide perovskite (with Br greater than 20%) is critical to create tandem solar cells. However, issues with Br–I phase segregation under solar cell operational conditions (such as light and heat) limit the device voltage and operational stability. This challenge is often exacerbated by the ready defect formation associated with the rapid crystallization of Br-rich perovskite chemistry with antisolvent processes. We combined the rapid Br crystallization with a gentle gas-quench method to prepare highly textured columnar 1.75–electron volt Br–I mixed WBG perovskite films with reduced defect density. Here, with this approach, we obtained 1.75–electron volt WBG PSCs with greater than 20% power conversion efficiency, approximately 1.33-volt open-circuit voltage (V oc ), and excellent operational stability (less than 5% degradation over 1100 hours of operation under 1.2 sun at 65°C). When further integrated with 1.25–electron volt narrow-bandgap PSC, we obtained a 27.1% efficient, all-perovskite, two-terminal tandem device with a high V oc of 2.2 volts.

14 SOLAR ENERGY↗