Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Step”

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

Monitoring Based Commissioning: Step-by-Step Approach to Implementation in Chiller Plants

Chilled water systems consume a significant portion of the energy in manufacturing, and inefficient operation could significantly influence utility bills and the environmental footprint of a facility. One way to ensure long-term energy-efficient plant operation is to ensure that chilled water system equipment, such as chillers, cooling towers, and pumps, are running at their best. Periodic operational assessments help to meet efficiency objectives for cooling plants. However, occasional set point and operational changes resulting from short-term fluctuations cause the plant operation to deviate from optimal. These changes accumulate over time, and operational efficiency deteriorates. Significant energy efficiency improvement potential in chilled water systems can be achieved by operational changes alone. Recent developments in information technology and the industrial internet of things data processing ability have overcome traditional problems and given rise to a new paradigm: monitoring-based commissioning (MBCx). Commissioning platforms and fault detection and diagnostics capabilities have evolved over time and presented more sophisticated abilities to efficiently control process cooling plants. The typical architectures of the MBCx process and a detailed approach to implement the process for chilled water systems are presented. A case study highlighting the energy-and cost-savings potential of MBCx systems is also presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

ZTF SN Ia DR2 follow-up: Exploring the origin of the Type Ia supernova host galaxy step through Si II velocities

The relation between Type Ia supernovae (SNe Ia) and the stellar masses of their host galaxy is well documented. In particular, Hubble residuals display a distinct luminosity shift based on host mass. This is known as the mass step. This effect is widely used as an additional correction factor in the standardisation of SN Ia luminosities. We investigate the Hubble residuals and the mass step of normal SNe Ia in the context of Si IIλ6355 velocities based on 277 normal SNe Ia that are near their peak in the second data release (DR2) of the Zwicky Transient Facility (ZTF). We divided the sample into high-velocity (HV) and normal-velocity (NV) SNe Ia, separated at 12,000 km s −1 . This produced a sample of 70 HV and 207 NV objects. We then explored potential environment- and/or progenitor-related effects by investigating the Si IIλ6355 velocities with parameters such as the light-curve stretch x 1 , the colour c, and the host galaxy properties. Although we only find a marginal difference between the Hubble residuals of HV and NV SNe Ia, the NV mass step is 0.149 ± 0.024 mag (6.3σ). The HV mass step is smaller, 0.046 ± 0.041 mag (1.1σ), and is consistent with zero. The difference between the NV and HV mass steps is modest, at ∼2.2σ. Moreover, the clearest subtype difference appears for SNe in central regions (d DLR < 1), where NV SNe Ia show a large mass step, whereas HV SNe Ia are consistent with no step, yielding a difference of 3.1–3.6σ between NV and HV SNe Ia. We observe a host-colour step for both subtypes. NV SNe Ia show a step of 0.142 ± 0.024 mag (5.9σ), while HV SNe Ia show a step of 0.158 ± 0.042 mag (3.8σ), where the HV SNe Ia step appears to be larger, but the significance is lower because the sample size is smaller. Overall, the NV and HV colour steps are statistically consistent. HV SNe Ia also show modest (∼2.5–3σ) steps in certain subsets, such as those in outer regions (d DLR > 1), whereas NV SNe display stronger environmental trends. Our results indicate that NV SNe Ia appear to be more environmentally sensitive, particularly in central likely metal-rich and older regions, while HV SNe Ia show weaker and subset-dependent trends. This suggests that applying a universal mass-step correction might introduce biases, and that incorporating refined classifications and/or environment-dependent factors, such as the location within the host, might improve future cosmological analyses beyond the standard x 1 and c cuts.

supernovae: general↗
DOE OSTI2022

Structure Matters: Asymmetric CO Oxidation at Rh Steps with Different Atomic Packing

Curved crystals are a simple but powerful approach to bridge the gap between single crystal surfaces and nanoparticle catalysts, by allowing a rational assessment of the role of active step sites in gas-surface reactions. Using a curved Rh(111) crystal, here, we investigate the effect of A-type (square geometry) and B-type (triangular geometry) atomic packing of steps on the catalytic CO oxidation on Rh at millibar pressures. Imaging the crystal during reaction ignition with laser-induced CO 2 fluorescence demonstrates a two-step process, where B-steps ignite at lower temperature than A-steps. Such fundamental dissimilarity is explained in ambient pressure X-ray photoemission (AP-XPS) experiments, which reveal partial CO desorption and oxygen buildup only at B-steps. AP-XPS also proves that A-B step asymmetries extend to the active stage: at A-steps, low-active O–Rh–O trilayers buildup immediately after ignition, while highly active chemisorbed O is the dominant species on B-type steps. We conclude that B-steps are more efficient than A-steps for the CO oxidation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development and application of two-step uncertainty propagation and sensitivity analysis methodology for fast reactor safety analysis

Uncertainty quantification (UQ) in nuclear reactors for transients is directly linked with safety assessment through the cross-sections uncertainties, provided as a covariance matrix, which are propagated through the reactor system to output of interest pertaining to reactor safety, such as peak temperatures in fuel/clad/coolant. Using a two-step approach, uncertainties are first quantified and propagated from basic input variables (such as reaction cross-sections) to intermediate quantities (such as reactivity feedback coefficients) through lattice level calculations. Uncertainties of intermediate quantities (from the first step) are then propagated through the system transient calculations, in the second step, to obtain uncertainties on reactor safety output parameters of interest. The scope of this work consists of Uncertainty Quantification & Propagation of nuclear data uncertainties that are highly correlated through unprotected transient overpower and unprotected loss of flow to assess their impact on core safety parameters. This two-step approach in the presence of covariance renders the sensitivity analysis very challenging. In fact, usually the sensitivity analysis is restricted to each step, which limits its application since the sensitivities between the system output quantities and the basic input variables are difficult to obtain. Here, in this work, we address this issue by proposing a simple, general methodology to combine the sensitivity indices obtained in each step by assuming the model behavior being linear. For the first step Generalized Perturbation theory based indices are used while in the second step the recently studied Johnson indices. The uncertainty quantification and sensitivity methodologies discussed here are demonstrated on a generic LFR design which is based on the 500 MWth demonstration Lead-cooled fast reactor (DLFR) using oxide fuel, developed by Westinghouse Electric Company (WEC).

42 - ENGINEERING↗
DOE OSTI2020

Experimental studies of the two-step scheme with an intense radioactive 132 Sn beam for next-generation production of very neutron-rich nuclei

The study of nuclei in the yet-unexplored, very neutron-rich, mid- to heavy-mass region of the nuclear chart, comprising perhaps half of all nuclei predicted to exist, is leading to new and/or upgraded facilities, as well as to research into the underlying production processes. In the present study, the usefulness of the two-step scheme, an alternate method to produce very neutron-rich nuclei, by a combination of an isotope-separation online (ISOL) system as a first step, and in-beam fragmentation of re-accelerated radioactive isotopes (RIs) as a second step, was investigated with a 132 Sn beam. In this work, very neutron-rich RIs around the neutron-rich neutron number N = 82 region were produced from the 278-MeV/nucleon 132 Sn beam impinging on a 5.97-mm Be target, and their production cross sections were measured. Yields were then estimated for the two-step scheme with the 132 Sn beam relative to the ones by a one-step scheme, in-flight fission of a 238 U beam, for 1-MW proton and 238 U beams at respective RI-beam facilities. This comparison suggests that the two-step scheme with the 132 Sn beam provides yields >40-times higher than those with the one-step scheme for the very neutron-rich N = 82 region. Moreover, by using various RI beams over the nuclear chart from ISOL, certain kinds of very neutron-rich RIs around the supernova r-process path can be produced with greater yields than by the one-step approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Quantifying and attributing time step sensitivities in present-day climate simulations conducted with EAMv1

Abstract. This study assesses the relative importance of time integration error in present-day climate simulations conducted with the atmosphere component of the Energy Exascale Earth System Model version 1 (EAMv1) at 1∘ horizontal resolution. We show that a factor-of-6 reduction of time step size in all major parts of the model leads to significant changes in the long-term mean climate. Examples of changes in 10-year mean zonal averages include the following: up to 0.5 K of warming in the lower troposphere and cooling in the tropical and subtropical upper troposphere, 1 %–10 % decreases in relative humidity throughout the troposphere, and 10 %–20 % decreases in cloud fraction in the upper troposphere and decreases exceeding 20 % in the subtropical lower troposphere. In terms of the 10-year mean geographical distribution, systematic decreases of 20 %–50 % are seen in total cloud cover and cloud radiative effects in the subtropics. These changes imply that the reduction of temporal truncation errors leads to a notable although unsurprising degradation of agreement between the simulated and observed present-day climate; to regain optimal climate fidelity in the absence of those truncation errors, the model would require retuning. A coarse-grained attribution of the time step sensitivities is carried out by shortening time steps used in various components of EAM or by revising the numerical coupling between some processes. Our analysis leads to the finding that the marked decreases in the subtropical low-cloud fraction and total cloud radiative effect are caused not by the step size used for the collectively subcycled turbulence, shallow convection, and stratiform cloud macrophysics and microphysics parameterizations but rather by the step sizes used outside those subcycles. Further analysis suggests that the coupling frequency between the subcycles and the rest of EAM significantly affects the subtropical marine stratocumulus decks, while deep convection has significant impacts on trade cumulus. The step size of the cloud macrophysics and microphysics subcycle itself appears to have a primary impact on cloud fraction in the upper troposphere and also in the midlatitude near-surface layers. Impacts of step sizes used by the dynamical core and the radiation parameterization appear to be relatively small. These results provide useful clues for future studies aiming at understanding and addressing the root causes of sensitivities to time step sizes and process coupling frequencies in EAM. While this study focuses on EAMv1 and the conclusions are likely model-specific, the presented experimentation strategy has general value for weather and climate model development, as the methodology can help researchers identify and understand sources of time integration error in sophisticated multi-component models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

High-order multirate explicit time-stepping schemes for the baroclinic-barotropic split dynamics in primitive equations

In order to treat the multiple time scales of ocean dynamics in an efficient manner, the baroclinic-barotropic splitting technique has been widely used for solving the primitive equations for ocean modeling. Based on the framework of strong stability-preserving Runge-Kutta approach, we propose two high-order multirate explicit time-stepping schemes (SSPRK2-SE and SSPRK3-SE) for the resulting split system in this paper. The proposed schemes allow for a large time step to be used for the three-dimensional baroclinic (slow) mode and a small time step for the two-dimensional barotropic (fast) mode, in which each of the two mode solves just need to satisfy their respective CFL conditions for numerical stability. Specifically, at each time step, the baroclinic velocity is first computed by advancing the baroclinic mode and fluid thickness of the system with the large time-step and the assistance of some intermediate approximations of the barotropicmode obtained by substepping with the small time step; then the barotropic velocity is corrected by using the small time step to re-advance the barotropic mode un-der an improved barotropic forcing produced by interpolation of the forcing terms from the preceding baroclinic mode solves; lastly, the fluid thickness is updated by coupling the baroclinic and barotropic velocities. Additionally, numerical inconsistencies on the discretized sea surface height caused by the mode splitting are relieved via a reconciliation process with carefully calculated flux deficits. Here, two benchmark tests from the “MPAS-Ocean” platform are carried out to numerically demonstrate the performance and parallel scalability of the proposed SSPRK-SE schemes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mechanistic Pathways for the Molecular Step Growth of Calcium Oxalate Monohydrate Crystal Revealed by In Situ Liquid-Phase Atomic Force Microscopy

Calcium oxalate monohydrate (COM) crystal is the most common crystalline component of human kidney stones. The molecular-scale inhibitory mechanisms of COM crystal growth by urinary biomolecules such as citrate and osteopontin adsorbed onto the crystal surface are now well understood. However, the pathways by which dissolved calcium and oxalate ions are incorporated into the molecular step of the COM crystal surface, leading to COM crystal growth—a prerequisite to be elucidated for developing effective therapeutics to inhibit COM stones—remain unknown. Here, using in situ liquid-phase atomic microscopy along with a step kinetic model, we reveal the pathways of the calcium and oxalate ions into the COM molecular step via the growth speed analysis of the molecular steps with respect to their step width at the nanoscale. Our results show that, primarily, the ions are adsorbed onto the terrace of the crystal surface from the solution—the rate-controlling stage for the molecular step growth, i.e., COM crystal growth—and then diffuse over it and are eventually incorporated into the steps. This primary pathway of the ions is unaffected by the model peptide D-Asp6 adsorbed on the COM crystal surface, suggesting that urinary biomolecules will not alter the pathway. These new findings rendering an essential understanding of the fundamental growth mechanism of COM crystal at the nanoscale provide crucial insights beneficial to the development of effective therapeutics for COM kidney stones.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Explicit Monotone Stable Super-Time-stepping Methods for Finite Time Singularities

We explore a novel way to numerically resolve the scaling behavior of finite-time singularities in solutions of nonlinear parabolic PDEs. The Runge–Kutta–Legendre (RKL) and Runge–Kutta–Gegenbauer (RKG) super-time-stepping methods were originally developed for nonlinear complex physics problems with diffusion. These are multistage single step second-order, forward-in-time methods with no implicit solves. The advantage is that the time-step size for stability scales with stage number 𝑠 as $\mathcal{O}$⁡(𝑠 2 ). Many interesting nonlinear PDEs have finite-time singularities, and the presence of diffusion often limits one to using implicit or semi-implicit time-step methods for stability constraints. Finite-time singularities are particularly challenging due to the large range of scales that one desires to resolve, often with adaptive spatial grids and adaptive time steps. Here, in this study, we show two examples of nonlinear PDEs for which the self-similar singularity structure has time and space scales that are resolvable using the RKL and RKG methods, without forcing even smaller time steps. Compared to commonly used implicit numerical methods, we achieve a significantly smaller run time while maintaining comparable accuracy. We also prove numerical monotonicity for both the RKL and RKG methods under their linear stability conditions for the constant coefficient heat equation, in the case of infinite domain and periodic boundary condition, leading to a theoretical guarantee of the superiority of the RKL and RKG methods over traditional super-time-stepping methods, such as the Runge-Kutta-Chebyshev and the orthogonal Runge-Kutta-Chebyshev methods. Code can be found at https://github.com/ZT220501/SRK-Singularity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Mixed precision s –step Lanczos and conjugate gradient algorithms

Compared to the classical Lanczos algorithm, the s-step Lanczos variant has the potential to improve performance by asymptotically decreasing the synchronization cost per iteration. However, this comes at a price; despite being mathematically equivalent, the s-step variant may behave quite differently in finite precision, potentially exhibiting greater loss of accuracy and slower convergence relative to the classical algorithm. It has previously been shown that the errors in the s-step version follow the same structure as the errors in the classical algorithm, but are amplified by a factor depending on the square of the condition number of the O(s)-dimensional Krylov bases computed in each outer loop. As the condition number of these s-step bases grows (in some cases very quickly) with s, this limits the s values that can be chosen and thus can limit the attainable performance. In this work, we show that if a select few computations in s-step Lanczos are performed in double the working precision, the error terms then depend only linearly on the conditioning of the s-step bases. This has the potential for drastically improving the numerical behavior of the algorithm with little impact on per-iteration performance. Our numerical experiments demonstrate the improved numerical behavior possible with the mixed precision approach, and also show that this improved behavior extends to mixed precision s-step CG. Here, we present preliminary performance results on NVIDIA V100 GPUs that show that the overhead of extra precision is minimal if one uses precisions implemented in hardware.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Shapiro steps and stability of skyrmions interacting with alternating anisotropy under the influence of ac and dc drives

Here we use atomistic simulations to examine the sliding dynamics of a skyrmion in a two-dimensional system containing a periodic one-dimensional stripe pattern of variations between low and high values of the perpendicular magnetic anisotropy. The skyrmion changes in size as it crosses the interface between two anisotropy regions. On applying combined dc and ac driving in either parallel or perpendicular directions, we observe a wide variety of Shapiro steps, Shapiro spikes, and phase-locking phenomena. The phase-locked orbits have two-dimensional dynamics due to the gyrotropic or Magnus dynamics of the skyrmions and are distinct from the phase-locked orbits found for strictly overdamped systems. Along a given Shapiro step when the ac drive is perpendicular to the dc drive, the velocity parallel to the ac drive is locked while the velocity in the perpendicular direction increases with increasing drive to form Shapiro spikes. At the transition between adjacent Shapiro steps, the parallel velocity jumps up to the next step value, and the perpendicular velocity drops. The skyrmion Hall angle shows a series of spikes as a function of increasing dc drive, where the jumps correspond to the transition between different phase-locked steps. At high drives, the Shapiro steps and Shapiro spikes are lost. When both the ac and dc drives are parallel to the stripe periodicity direction, Shapiro steps appear, while if the dc drive is parallel to the stripe periodicity direction and the ac drive is perpendicular to the stripe periodicity, then there are only two locked phases, and the skyrmion motion consists of a combination of sliding along the interfaces between the two anisotropy values and jumping across the interfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

First Solar Thermal Energy Planner (STEP 1) and Nationwide Industrial Heat and Power Analysis

The First Solar Thermal Energy Planner (STEP 1) and Nationwide Industrial Heat and Power Analysis (aka the STEP 1 Project or the Project) aimed to (1) developed a brand-new web tool that could provide decision support through free, rapid techno-economic analysis of behind-the-meter solar+storage systems for industrial process heat and (2) conduct high-level analyses of the cost-competitiveness of the same systems across the US in key sectors. The STEP 1 web tool collects key location, land availability, thermal load profile, proccess heat temperature and media, and other key parameters through an easy-to-use user interface (UI). The UI was designed to meet the user at their level of understanding by minimizing the number of required inputs as much as possible while including options for more nuanced inputs if the user desires. STEP 1 advises users on which solar thermal tehcnologies that could fit their needs based on the inputs provided (namely process media and temperature). The tool can model a wide range of solar thermal technologies including flat plate collectors, evacuated tubes, parabolic troughs, linear Fresnel, and molten salt towers all with corresponding thermal energy storage (TES) - solar PV with resistive heating and TES is also included. Once the parameters are collected, a nominal thermal energy production profile for the facility's location is generated using NREL's System Advisory Model (SAM) and then passed to a modified version of NREL's REopt platform to optimize the size (capacity) and dispatch of the solar+storage system to minimize lifecycle costs subject to energy balance, fuel and electricity rates, emission reductions goals, and other constraints. This entire process takes less than 20 minutes, is completely free, requires zero coding skills, and provides the user with a high-level assessment on the techno-economic feasability of deploying solar+storage systems for their energy needs. In addition to the development of the STEP 1 web tool, the project completed two complementary analyses focused on leveraging the backend code of STEP 1, public industrial facility locations and fuel consumption data, and key sector information to assess the economic opportunity of reducing fuel costs at various levels of capacity factor around the US. Both analyses found that there are key markets, locations, industrial sectors (namely food and beverage), and degrees of offset where solar thermal technologies could be cost-effectively deployed, highlighting a key market entry point for these technologies, and assessing deployment potential. In summary, the STEP 1 project improved the opportunities for solar+storage systems to expand into the industrial process heat market through breaking down barriers to assessing these technologies.

14 SOLAR ENERGY↗
DOE OSTI2021

A Two-Step PCR Protocol Enabling Flexible Primer Choice and High Sequencing Yield for Illumina MiSeq Meta-Barcoding

High-throughput amplicon sequencing that primarily targets the 16S ribosomal DNA (rDNA) (for bacteria and archaea) and the Internal Transcribed Spacer rDNA (for fungi) have facilitated microbial community discovery across diverse environments. A three-step PCR that utilizes flexible primer choices to construct the library for Illumina amplicon sequencing has been applied to several studies in forest and agricultural systems. The three-step PCR protocol, while producing high-quality reads, often yields a large number (up to 46%) of reads that are unable to be assigned to a specific sample according to its barcode. Here, we improve this technique through an optimized two-step PCR protocol. We tested and compared the improved two-step PCR meta-barcoding protocol against the three-step PCR protocol using four different primer pairs (fungal ITS: ITS1F-ITS2 and ITS1F-ITS4, and bacterial 16S: 515F-806R and 341F-806R). We demonstrate that the sequence quantity and recovery rate were significantly improved with the two-step PCR approach (fourfold more read counts per sample; determined reads ≈90% per run) while retaining high read quality (Q30 > 80%). Given that synthetic barcodes are incorporated independently from any specific primers, this two-step PCR protocol can be broadly adapted to different genomic regions and organisms of scientific interest.

16S rDNA↗
DOE OSTI2022

LEAKAGE CHARACTERISTICS OF A STEPPED BALANCE PISTON LABYRINTH SEAL IN SUPERCRITICAL CO2

In an integrally geared compressor, there is a net thrust load on the pinion shaft from gear loads and aerodynamic forces. The net thrust on the pinion is resolved with either a thrust bearing or a thrust collar (rider ring). As the pressure increases, this increases the aerodynamic force proportionally. In highpressure systems this can overwhelm the thrust bearing or thrust collar if thrust is mismanaged. In a recent test program for an sCO2 integrally geared compander (compressor/turbine), a balance piston was created on the opposite end of the turbine to provide an equalizing thrust; while this is not novel, the balance piston provided a test apparatus for measuring the leakage on a 20 tooth stepped labyrinth seal operating at pressures in excess of 250 bar (3600 psi). The stepped balance piston seal is a 20 tooth tooth-onrotor balance piston seal. This type of seal was selected since it was predicted to have substantially less leakage than a comparably sized straight through labyrinth seal. Comparisons between the two designs will be made, showing that the leakage of the stepped balance piston seal is considerably lower. The tooth clearance on the seal was selected to be 30 mils radial. While this seems large, a stack-up tolerance assessment showed that this magnitude of clearance was appropriate to prevent rotor stator contact. During the testing campaign, data was recorded showing the leakage across this stepped balance piston seal at pressures in excess of 250 bar (3600 psi) with a pressure ratio of 2. This leakage data was measured, and will be compared to leakage data predicted with CFD. Two measurement methods are presented for the stepped balance piston seal: an orifice plate and valve. The predicted seal leakage from a straight through seal was four times that of the predicted stepped labyrinth seal. The measured flow from the stepped seal was on the same order as the CFD prediction and offered adequate rotordynamic performance

wilkes, jason↗
DOE OSTI2024

Models implemented in the methodological approach to design the initial STEP first wall contour

The official Spherical Tokamak for Energy Production mission aims to demonstrate the ability to generate net electricity from fusion with the STEP Prototype Power plant. One of the key technological and engineering challenges in fusion power plants is managing the loads on the first wall within acceptable limits. Therefore, the conceptual design development of the STEP Prototype Power plant needs to be based on load estimates derived using legitimate plasma physics assumptions through dynamic and flexible tools. The current design foresees the STEP main chamber first wall to withstand steady-state heat loads of up to ~1 MW/m 2 , excluding critical regions expected to receive higher heat loads such as the baffle regions approaching the divertors. These critical areas will require ad hoc assessments and will be designed with the presence of limiters. This article focuses on the models and methodology adopted for designing the 2-D poloidal contour of the STEP first wall, based on the anticipated charged particle and radiation heat loads during normal operation. Firstly, the models adopted for calculating the charged particle and radiation heat loads are introduced. The first model is validated through benchmarking against the particle tracing code SMARDDA, while the second model is verified by comparing it with data from the MAST-U experiment. Secondly, the model used to design the 2-D first wall contour according to the heat loads is explained. We acknowledge that this preliminary design stage assumes certain simplifications, notably an axisymmetric geometry, for computational efficiency and clarity in presentation. It is understood that subsequent design phases will address the complexities of real-world engineering, including non-axisymmetric effects, transient plasma scenarios, and the impact of disruptions on the first wall design. Finally, an automatic procedure based on these models is presented for defining the 2-D poloidal contour of the STEP first wall to minimize heat loads, taking into account the need to radiate most of the alpha-particle and auxiliary heating power. Here, by providing an overview of the models, methodology, and an automatic procedure, this paper contributes to the design process of the STEP first wall, addressing the engineering challenges associated with fusion power plant development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Storm Surge Modeling as an Application of Local Time-Stepping in MPAS-Ocean

This paper presents the first practical application of local time-stepping (LTS) schemes in the Model for Prediction Across Scales-Ocean (MPAS-O). We use LTS schemes in a single-layer, global ocean model that predicts the storm surge around the eastern coast of the United States during Hurricane Sandy. The variable-resolution meshes used are of unprecedentedly high resolution in MPAS-O, containing cells as small as 125 m wide in Delaware Bay. It is shown that a particular, third-order LTS scheme (LTS3) produces sea-surface height solutions that are of comparable quality to solutions produced by the classical four-stage, fourth-order Runge-Kutta method (RK4) with a uniform time step on the same meshes. Furthermore, LTS3 is up to 35% faster in the best cases considered, where the number of cells using the coarse time-step relative to those using the fine time-step is as low as 1:1. This shows that LTS schemes are viable for use in MPAS-O with the added benefit of substantially less computational cost. The results of these performance experiments inform us of the requirements for efficient mesh design and configuration of LTS regions for LTS schemes. In particular, we see that for LTS to be efficient on a given mesh, it is important to have enough cells using the coarse time-step relative to those using the fine time-step, typically at least 1:5 to see an increase in performance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Scalable Two–Step Production of High–Efficiency Perovskite Solar Cells and Modules

Perovskite solar cells (PSCs) fabricated in laboratories have already achieved power conversion efficiency (PCE) comparable to market-dominant crystalline silicon solar cells. However, this promising photovoltaic technology suffers from severe loss of PCE during scaling up, limiting its progress toward commercialization. One critical question is to develop scalable, low-cost, high throughput, and well-controlled production methods to deposit high-quality perovskite films. Among various approaches, two-step sequential deposition methods have their unique advantages but have been long overlooked. This review provides an overview of two-step methods for fabricating efficient and stable perovskite solar modules (PSMs). We first discuss the mechanisms of two-step perovskite conversion and advanced engineering approaches to modulate the perovskite formation process. We survey the progress of efficient PSCs prepared by different two-step methods and compare the advantages and disadvantages of each method for scalable production of PSMs. Particularly, we highlight that the vapor-based two-step methods are promising for high-throughput in-line production of PSMs. Finally, we provide insights into the challenges and outlook of the industrialization of two-step processes for producing PSMs.

14 SOLAR ENERGY↗
DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗