Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak grids”

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 163 records · Page 9

A cohesive zone treatment for the material point method involving problems of large deformation and damage

A new algorithm is described that permits the use of cohesive zones in the material point method for problems involving large deformation and fracture. In contrast to previous cohesive zone implementations, this method does not utilize massless surface-element particles. Instead, cohesive tractions are computed using the shape function mappings from a reference grid configuration in combination with explicitly defined particle surface normals and surface positions. These normals and relative surface positions are updated each time step according to particle deformation. The tractions are converted to cohesive forces using the nodal areas and mapped back to particles using the same reference shape function mappings. These forces are then remapped by conventional particle-to-grid interpolation as external forces using the current-configuration shape-function mappings. This allows highly compliant cohesive zones to function over jump displacements larger than a grid cell. Upon damage, these interfaces can revert to conventional multi-field contact surfaces. This approach is general and readily applies to two and three dimensions as well as being compatible with damage-field gradient partitioning offering exceptional computational flexibility. The framework for this method enables other capabilities, such as improved contact precision using explicitly defined surface normals and positions, and a method to mitigate spurious material damage at weak discontinuities between stiff brittle materials and soft or compliant materials.

Cohesive zone↗

Evaluation of Mesoscale Convective Systems in Climate Simulations: Methodological Development and Results from MPAS-CAM over the U.S.

In this study, a process-oriented approach is developed to evaluate warm-season mesoscale convective system (MCS) precipitation and their favorable large-scale meteorological patterns (FLSMPs) over the U.S. This approach features a novel observation-driven MCS-tracking algorithm using infrared brightness temperature and precipitation feature at 12, 25 and 50 km resolution and metrics to evaluate the model large-scale environment favorable for MCS initiation. The tracking algorithm successfully reproduces the observed MCS statistics from a reference 4-km radar MCS database. To demonstrate the utility of the new methodologies in evaluating MCS in climate simulations with mesoscale resolution, the process-oriented approach is applied to two climate simulations produced by the Variable-Resolution Model for Prediction Across Scales coupled to the Community Atmosphere Model physics, with refined horizontal grid spacing at 50 km and 25 km over North America. With the tracking algorithm applied to simulations and observations at equivalent resolutions, the simulated number of MCS and associated precipitation amount, frequency and intensity are found to be consistently underestimated in the Central U.S., particularly from May to August. The simulated MCS precipitation shows little diurnal variation and lasts too long, while MCS precipitation area is too large and intensity is too weak. The model is able to simulate four types of observed FLSMP associated with frontal systems and low-level jets (LLJ) in spring, but the frequencies are underestimated because of low-level dry bias and weaker LLJ. Precipitation simulated under different FLSMPs peak during daytime, in contrast to the observed nocturnal peak. Implications of these findings for future model development and diagnostics are discussed.

54 ENVIRONMENTAL SCIENCES↗

Feasibility and strategic implications of deploying nuclear power reactors in Africa

This report assesses the feasibility and strategic implications of deploying nuclear power reactors, including large-scale plants, advanced small modular reactors (SMRs), and microreactors, in African countries. Case studies focus on South Africa, Egypt, Kenya, Ghana, and Nigeria, examining nuclear energy’s role in Africa’s rapidly evolving energy landscape, marked by fast-growing demand, significant electricity access gaps, increasing renewable penetration, and strong policy commitments to industrialization and energy security. Several U.S. reactor technologies and designs are considered based on their development status and readiness for deployment. The analysis finds that nuclear power can provide reliable, clean baseload and flexible generation, as well as high-temperature process heat for desalination, hydrogen production, and industrial applications. However, suitability is highly country-specific, depending on grid size and stability, transmission capacity, cooling water availability, regulatory readiness, and fuel supply chains. Near-term deployment opportunities are strongest for light-water reactors (such as NuScale, BWRX-300, AP300, and SMR-300) that use low-enriched uranium and build on proven technology. More advanced concepts, including gas-cooled, sodium-cooled, molten-salt cooled reactors, and microreactors, will likely be relevant for African deployment in the 2030s or later, contingent on demonstration projects, high-assay low-enriched uranium (HALEU) fuel availability, and mature international licensing frameworks. Economic analysis shows that SMRs are capital-intensive, with projected overnight costs for 300 MWe units in 2025 ranging from approximately 1.4 to 2.6 billion USD per module. The levelized cost of electricity (LCOE) is highly sensitive to the weighted average cost of capital (WACC). Given typically higher financing costs and utility balance-sheet weaknesses in many African countries, bankable project structures will require sovereign guarantees, robust offtake arrangements, and layered financing from export credit agencies, development finance institutions, and vendor nations. Comparisons with recent large nuclear projects in the United Arab Emirates (UAE) and Egypt underscore the central role of state-backed loans, long tenors, and concessional terms. Country case studies illustrate a spectrum of readiness and opportunity. South Africa operates two 920 MWe pressurized light water reactors (totaling 1,840 MWe) at Koeberg and has the most mature regulatory and industrial base, positioning it as a prime candidate for both large reactors and SMRs to replace coal, support desalination, and anchor industrial hubs. Egypt is constructing four VVER-1200 units at El Dabaa with strong state leadership and could later complement this fleet with SMRs for coastal and industrial applications. Kenya and Ghana are advancing through IAEA Milestones with growing institutional capacity and clear interest in SMRs that match their smaller grids and industrialization plans. Nigeria has the largest demand potential but faces acute constraints in grid reliability, project bankability, and regulatory capacity; targeted deployments of large reactors and SMRs near coastal or industrial sites could have high impact if accompanied by major grid upgrades and institutional reforms. The report identifies cross-cutting challenges such as financing, political continuity, public acceptance, nonproliferation and security, waste and back-end management, regulatory capacity, grid adequacy, and long deployment timelines for first-of-a-kind designs, and ANL/NSE-26/3 ii proposes broad directions for resolution. These include stronger multifaceted financing for nuclear, long-term national energy strategies that transcend electoral cycles, proactive stakeholder engagement, strengthened regional and national regulators, and systematic workforce development through centers of excellence and expanded training. The United States should develop partnerships with African countries and offer end-to-end nuclear package similar to those used effectively by competitors: coordinated project development, state-backed financing, long-term fuel services, and durable in-country support through regional offices and sustained workforce/regulatory training. With timely planning, sustained political commitment, and appropriate financing and institutional support, nuclear energy, both large reactors and advanced SMRs, can become a meaningful, though not dominant, pillar of Africa’s future power mix, enhancing energy security, enabling industrial growth, and supporting climate goals.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The BTSbot-nearby Discovery of SN 2024jlf: Rapid, Autonomous Follow-up Probes Interaction in an 18.5 Mpc Type IIP Supernova

We present observations of the Type IIP supernova (SN) SN 2024jlf, including spectroscopy beginning just 0.7 days (∼17 hr) after first light. Rapid follow-up was enabled by the new BTSbot-nearby program, which involves autonomously triggering target-of-opportunity requests for new transients in Zwicky Transient Facility data that are coincident with nearby (D < 60 Mpc) galaxies and identified by the BTSbot machine learning model. Early photometry and nondetections shortly prior to first light show that SN 2024jlf initially brightened by >4 mag day −1 , quicker than ∼90% of Type II SNe. Early spectra reveal weak flash ionization features: narrow, short-lived (1.3 < τ[days] < 1.8) emission lines of Hα, He II , and C IV . Assuming a wind velocity of v w = 50 km s −1 , these properties indicate that the red supergiant progenitor exhibited enhanced mass loss in the last year before explosion. We constrain the mass-loss rate to $1{0}^{-4}\lt \dot{M}\,[{M}_{\odot }\,{\mathrm{yr}}^{-1}]\lt 1{0}^{-3}$ by matching observations to model grids from two independent radiative hydrodynamics codes. BTSbot-nearby automation minimizes spectroscopic follow-up latency, enabling the observation of ephemeral early-time phenomena exhibited by transients.

core-collapse supernovae↗

Improvement of Numerical Methods in Pronghorn

During the fiscal year of 2021, the finite-volume method (FVM) was deployed in Pronghorn to provide improved efficiency, stability, and accuracy for coarse-mesh, thermal-hydraulics problems. While the main goals of fiscal year 2021 were met, several issues emerged from the early deployment of the finite-volume method in Pronghorn. These issues were: The compressible and incompressible formulations that were implemented are inadequate for many nuclear reactor flow problems; Omission of terms accounting for the porosity and Darcy-Forchheimer body force discontinuities in the Rhie-Chow interpolation lead to oscillations in pressure and velocity at these discontinuities; The FVM lacks a correction for non-orthogonal grids for computing accurate pressure gradient leading to loss of accuracy in regions with skewed elements; The FVM currently uses a monolithic solver. Monolithic solvers have issues dealing with the saddle-point nature of the discretized fluid equations. The result are bad convergence if direct factorization is not used and large memory consumption when direct factorization is used. These four issues are addressed in this report. In particular, we report the completion of the following task: Implementation of a weakly compressible formulation in the MOOSE Navier-Stokes module; Implementation of Moukalled’s method for including body forces in the RCI. Additionally, we identified the need to smooth the porosity using Moukalled’s face-cell smoothing operator; Implementation of a non-orthogonal correction for the Green-Gauss gradient computation; Preliminary implementation of a SIMPLE segregated solver.

97 MATHEMATICS AND COMPUTING↗

Case Study: Applying the Idaho National Laboratory Resilience Framework to St. Mary’s, Alaska

The Idaho National Laboratory (INL) resilience framework has been developed to broadly apply to EEDS so that all elements of systems that contain distributed wind can be part of the resilience evaluation. The users or audience for this framework can include any stakeholders associated with the EEDS. Not all electrical energy systems have the same stakeholders; customers, owners, and operators are generally present but have different interests. Considering the broad electrical grid, customers, regulators, investors, utility planners, engineers, and operators each have an interest in system resilience driven from different motivating factors. This document focuses on the planning stage of the framework. In this document, each step is explained briefly before demonstrating its application to the St. Mary’s-Mt. Village system. The framework can be used for many types of resilience planning. It can be used to evaluate current overall resilience, or the resilience of certain subsystems. It can be used to explore existing resilience weak points and propose mitigations. It can also be used to evaluate the resilience benefits of a new investment. We use the latter application for this case study. Although the wind turbine in St. Mary’s has already been installed, the resilience benefits that the turbine provided were not well defined. It was installed with the main objective to generate electric power from a renewable resource in an effort to reduce the local dependency on fuel oil as the sole source of electric power generation, which is a resilience goal on its own, but there are other ways in which the turbine can add resilience to the system, as well as scenarios of interest to analyze how resilient the wind turbine itself is against different hazards. In this case study, we analyze the operation of the St. Mary’s power system both with the wind installed and without the wind installed during different resilience hazards of interest. This allows us to compare the performance with wind and without wind and to quantify the resilience benefits provided by wind. Our MIRACL partners at PNNL will then take the resilience benefits and assign value to the resilience provided by wind based on costs and costs avoided in the different scenarios.

17 WIND ENERGY↗

Towards a Unified Setup to Simulate Mid‐Latitude and Tropical Mesoscale Convective Systems at Kilometer‐Scales

Abstract Mesoscale convective systems (MCSs) are the main source of precipitation in the tropics and parts of the mid‐latitudes and are responsible for high‐impact weather worldwide. Studies showed that deficiencies in simulating mid‐latitude MCSs in state‐of‐the‐art climate models can be alleviated by kilometer‐scale models. However, whether these models can also improve tropical MCSs and whether we can find model settings that perform well in both regions is understudied. We take advantage of high‐quality MCS observations collected over the Atmospheric Radiation Measurement (ARM) facilities in the US Southern Great Plains (SGP) and the Amazon basin near Manaus (MAO) to evaluate a perturbed physics ensemble of simulated MCSs with 4 km horizontal grid spacing. A new model evaluation method is developed that enables to distinguish biases stemming from spatiotemporal displacements of MCSs from biases in their reflectivity and cloud shield. Amazon MCSs are similarly well simulated across these evaluation metrics than SGP MCSs despite the challenges anticipated from weaker large‐scale forcing in the tropics. Generally, SGP MCSs are more sensitive to the choice of model microphysics, while Amazon cases are more sensitive to the planetary boundary layer (PBL) scheme. Although our tested model physics combinations had strengths and weaknesses, combinations that performed well for SGP simulations result in worse results in the Amazon basin and vice versa. However, we identified model settings that perform well at both locations, which include the Thompson and Morrison microphysics coupled with the Yonsei University (YSU) PBL scheme and the Thompson scheme coupled with the Mellor‐Yamada‐Janjic PBL scheme.

54 ENVIRONMENTAL SCIENCES↗

Combined Cycle Integrated Thermal Energy Storage “CiTES” (Final Scientific/Technical Report)

The Phase I of this project confirmed the technical feasibility of a Combined Cycle integrated Thermal Energy Storage “CiTES” system, calculated the key performance parameters like power efficiency and costs, and proved its commercial value with full-year simulations for several US electricity markets with high degree of variable renewable generation and volatile hourly electricity prices. The core element of this project is the Electro Thermal Energy Storage (ETES) technology from Siemens Gamesa Renewable Energy GmbH, using thermally stable and inexpensive volcanic rocks as storage material and air as heat transfer medium. This technology is backed by more than 10 years of experience and a 440MMBTU (130MWh-th) pilot plant in Hamburg, Germany, which is in operation since 2019. The integration of this thermal storage in an existing combined cycle power plant (CCPP) is typical power plant technology without any major technology risks. It allows the storage of inexpensive renewable energy during times of surplus renewable generation and the discharge of this energy in times of high energy demand when the fossil plant is in operation. This supplements the fossil power generation with CO2-emission-free energy. The secondary effect of the CiTES system is that a small part of the stored thermal energy is used to keep the heat recovery steam generator (HRSG) and steam turbine (ST) of the combined cycle power plant in hot and ready-to-start condition. This enables the plant to start rapidly when fossil generation is required to satisfy demand as soon variable generation drops off in the evenings or during cloud cover and calm wind periods. Without pre-warming of the HRSG and ST, the CCPP would need several hours for a cold or warm start, burn a lot of gas and release high NOx emissions during start and wouldn’t be able to use the short times of high energy prices in an efficient or economical manner. The economic parameters of CiTES were determined by a full year “8760” simulation using a data set calculation for each of the hours of the year, and historical electricity and gas prices. For consistency, the simulations were focused on the pre-COVID year 2019. The financially most attractive markets were in the Energy Reliability Council of Texas (ERCOT) region, which allowed substantial value generation with arbitrage (charge with cheap energy during renewable surplus times and discharge when energy is needed and expensive). The improvement of flexibility with the CiTES system by pre-heating and warm-keeping of the CCPP allowed for additional power generation during short time periods when demand is high but renewable generation is down; when the hourly energy prices are highest in these markets. The simulations are based on 2019 data, when ERCOT had 27GW of installed photovoltaic (PV) and wind generation. They showed that the created revenue with the prototypically sized CiTES system of 1,000MMBTU (300MWh-th) falls a little bit short of what is expected from a commercially viable investment. The system has optimization opportunities for cost reduction and increased effectiveness which will be realized during a potential Phase II Pre-FEED study following this project. Furthermore, it is safe to assume that a lot of renewable generation capacity will be added all over the US in the coming years. As an example, ERCOT is predicting to more than double its renewable generation from 27GW in 2019 to a forecasted 63GW in 2023. This will increase the amount of renewable overproduction exponentially. This rapid increase of local overproduction and the need to curtail renewable generation is well documented by the California ISO (www.CAISO.com / managing oversupply). However, the simulations also revealed a weakness in the structure of the electricity markets in the US. More specifically, when electricity prices are very low and approaching negative levels, the owners of Variable Renewable Energy (VRE) will curtail a part of their facility to stabilize the price by reducing supply. This results in a situation in which storage facilities, which are integrated in existing fossil assets and don’t have the behind-the-meter benefit of a VRE, won’t be able to purchase low cost – otherwise curtailed – renewable energy off the grid. A special tariff, which motivates VRE owners to sell otherwise curtailed renewable energy to storage facilities (Hydrogen, thermal, pumped hydro, etc.) can solve this issue. The implementation of such a regulating tariff by Independent System Operators, thus avoiding renewable curtailment, is a pre condition for successful commercialization for renewable energy storage technologies. With this advancement of design and technology and improvements in the market environment, it can be expected that the Combined Cycle integrated Thermal Energy Storage proves itself as an important innovation to keep highly efficient, natural gas-based power generation economically successful and relevant for the power industry in the United States of America.

Wolf, Thorsten↗

Representation of atmosphere-induced heterogeneity in land–atmosphere interactions in E3SM–MMFv2

In the Energy Exascale Earth System Model (E3SM) Multi-scale Modeling Framework (MMF), where parameterizations of convection and turbulence are replaced by a 2-D cloud-resolving model (CRM), there are multiple options to represent land–atmosphere interactions. Here, we propose three different coupling strategies, namely the (1) coupling of a single land surface model to the global grid (MMF), (2) coupling a single land copy directly to the embedded CRM (SFLX2CRM), and (3) coupling a single copy of land model to each column of the CRM grid (MAML). In the MAML (Multi-Atmosphere Multi-Land) framework, a land model is coupled to CRM at the CRM-grid scale by coupling an individual copy of a land model to each CRM grid. Therefore, we can represent intra-CRM heterogeneity in the land–atmosphere interaction processes. There are 5-year global simulations run using these three coupling strategies, and we find some regional differences but overall small changes with respect to whether a land model is coupled to CRM or a global atmosphere. In MAML, the spatial heterogeneity within CRM induces stronger turbulence, which leads to the changes in soil moisture, surface heat fluxes, and precipitation. However, the differences in the MAML from the other two cases are rather weak, suggesting that the impact of using MAML does not justify the increase in cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnetic fields in the formation of the first stars – II. Results

ABSTRACT Beginning with cosmological initial conditions at z = 100, we simulate the effects of magnetic fields on the formation of Population III stars and compare our results with the predictions of Paper I. We use gadget-2 to follow the evolution of the system while the field is weak. We introduce a new method for treating kinematic fields by tracking the evolution of the deformation tensor. The growth rate in this stage of the simulation is lower than expected for diffuse astrophysical plasmas, which have a very low resistivity (high magnetic Prandtl number); we attribute this to the large numerical resistivity in simulations, corresponding to a magnetic Prandtl number of order unity. When the magnetic field begins to be dynamically significant in the core of the minihalo at z = 27, we map it on to a uniform grid and follow the evolution in an adaptive mesh refinement, MHD simulation in orion2. The non-linear evolution of the field in the orion2 simulation violates flux-freezing and is consistent with the theory proposed by Xu & Lazarian. The fields approach equipartition with kinetic energy at densities ∼1010–1012 cm−3. When the same calculation is carried out in orion2 with no magnetic fields, several protostars form, ranging in mass from ∼1 to 30 M⊙; with magnetic fields, only a single ∼30 M⊙ protostar forms by the end of the simulation. Magnetic fields thus suppress the formation of low-mass Pop III stars, yielding a top-heavy Pop III initial mass function and contributing to the absence of observed Pop III stars.

79 ASTRONOMY AND ASTROPHYSICS↗

A protocol and analysis of year-long simulations of global storm-resolving models and beyond

We propose a protocol to evaluate and analyze year-long simulations of global storm-resolving models (GSRMs). The proposed protocol complements an earlier 40-day simulation protocol under the DYAMOND (DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains) project to allow the analysis of the seasonal cycle and associated climatic relevant phenomena. This intercomparison aims to reveal how GSRMs, which can simulate mesoscale convective systems (MCSs) in the global domain, reproduce atmospheric large-scale structures related to convection beyond month-long simulations. The intercomparison for one-year simulations is conducted by either atmosphere-only models or atmosphere–ocean coupled models with atmospheric horizontal mesh sizes less than 5 km. We recommend the continuous four seasons from March 2020 to February 2021 as a target period for the intercomparison but with options for many groups to join more flexibly. The output variables are collected at 0.25° resolution, and archives of a small set of native grid variables are encouraged to analyze tropical cyclones and MCSs. Through the proposed global storm-resolving simulation, we will evaluate the climatological distributions of the atmospheric large-scale circulations, such as the Intertropical Convergence Zone (ITCZ), monsoon, midlatitude jets, their time evolution, and the upscale impacts on them. We present sample analyses from a one-year simulation using the 3.5 km mesh Nonhydrostatic Icosahedral Atmospheric Model (NICAM), revealing the realistic zonal contrast of tropical precipitation, no double ITCZ structure, the reasonable midlatitude jet position and intensity but a weak bias of storm track activities, and a warm bias over the Eurasia during boreal winter. We also clarify the cross-scale interaction, such as the effects of cold pools on mean precipitation over the Maritime Continent through the precipitation diurnal cycle and the effects of resolved gravity waves on midlatitude mean flows. The proposed one-year simulation protocol is referred to as the “Sendai Protocol.” This protocol is not unique or definite for evaluating GSRMs; we prospect a hierarchical set of experiments from short-term to multi-year simulations as GSRM intercomparisons.

54 ENVIRONMENTAL SCIENCES↗

Department of Energy Project DE-SC0012198: Probing Electron-Molecule Scattering with HHG: Theory and Experiment (Final Report)

The study of the motion of atoms in molecules is important to understanding many areas of physical and life sciences. Such motion occurs on many different times scales, with electronic motion occurring on a sub-femtosecond time scale, and simple vibrational motion in tens to hundreds of femtoseconds. One way to follow such processes in real time is by the use of shortpulsed lasers, and in particular by studying time-resolved photoionization and the related process of high-harmonic generation (HHG). Thus there has been much effort to develop the tools necessary to probe molecular systems using short pulse lasers and understanding the sensitivity of the different possible probes to the time dependent geometric structure as well as the electronic structure of molecules. Our research has particularly focused on the connection between high-field processes and the more traditional weak field photoionization processes. Strong field and weak field processes can be connected through models that involve the same matrix elements. We have demonstrated in our study of HHG from SF 6 that the spectrum is sensitive to the interplay between the angular dependence of the ionization step and recombination step. In our study of rescattering spectroscopy, we have shown that with a combination of experiment and theory, we can use this high-field spectroscopy to determine molecular structure in molecules such as C 2 H 4 . We have also developed new computational tools based on overset grids to enable studies on larger molecular systems which use much more robust numerical approaches so that the resulting code can be a tool that non-specialists can use to study related systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport

OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Alternating Direction Decomposition with Strong Bounding and Convexification (ADDSBC) for Solving Security Constrained AC Unit Commitment Problems

This project aims to develop efficient and robust computational methods for solving the security-constrained unit commitment and alternating current optimal power flow problem (SC-UC-ACOPF). The SC-UC-ACOPF problem is at the center of the short-term operation of the U.S. Power Grid. It is solved every week, every day, and every 10 minutes to plan for the optimal action of electricity generation and consumption by minimizing the generation cost and maintaining power system reliability against potential disruptions of equipment failures. In mathematical terms, SC-UC-ACOPF is a challenging large-scale mixed-integer nonlinear optimization model. This means that the decisions involve both discrete variables, e.g. the turning on and off of generators and switching of transmission lines and transformers, and continuous decisions, e.g. the amount of energy generated by each generator and the power flows in the power grid. The physics of the power flow is described by nonlinear equations involving real and reactive power and bus voltages. Another key feature is the large number of contingencies, i.e. the system needs to stay reliable in face of failure of any one equipment, such as transmission lines and generators. The U.S. power grids are extremely complicated and large scale with more than 5,000 generators, 50,000 buses, and 100,000 high-voltage transmission lines, making the SC-UC-ACOPF a very large-scale computation challenge. The research developed in this project aims to solve the SC-UC-ACOPF problems in the three timescales, i.e. weekly, daily, and every 10-min. The proposed computational methods are built on a principled algorithmic approach of decomposition and penalization. More specifically, the algorithm develops spatial and temporal decomposition by exploiting the strong temporal coupling and weak spatial coupling of the UC problem and the complementary feature, i.e. weak temporal coupling and strong spatial coupling of the ACOPF problem. The algorithm also leverages recent progresses in strong convex relaxation of ACOPF. A unique feature of the proposed approach is that it generates a valid, global upper bound on the optimal maximum profit. In this way, a global optimality gap is available to measure the quality of the solution. To further speed up computation, the research team has developed a plethora of effective heuristics to strengthen the iterative penalty-based decomposition framework. For instance, a heuristic is developed to construct inner approximations of the time coupling constraints within the time decoupled problems. Contingencies are pre-screened and low-rank matrix computation is exploited to find the almost unique solution to each contingency. A novel heuristic for line switching is proposed and tested with positive impacts on instances where line switching is beneficial. Taking a systematic approach and carefully handling every detail of the problem pays off. The TIM-GO’s performance throughout the trials and the final event was stellar. TIM-GO garnered the second highest total prize money and is ranked in the top three positions across all categories of comparison.

97 MATHEMATICS AND COMPUTING↗

Study of few-electron backgrounds in the LUX-ZEPLIN detector

The LUX-ZEPLIN (LZ) experiment aims to detect rare interactions between dark matter particles and xenon. Although the detector is designed to be the most sensitive to GeV/𝑐 2 –TeV/𝑐 2 weakly interacting massive particles (WIMPs), it is also capable of measuring low-energy ionization signals down to a single electron that may be produced by scatters of sub-GeV/𝑐 2 dark matter. The major challenge in exploiting this sensitivity is to understand and suppress the ionization background in the few-electron regime. We report a characterization of the delayed electron backgrounds following energy depositions in the LZ detector under different detector conditions. In addition, we quantify the probability for photons to be emitted in coincidence with electron emission from the high voltage grids. We then demonstrate that spontaneous grid electron emission can be identified and rejected with a high efficiency using a coincident photon tag, which provides a tool to improve the sensitivity of future dark matter searches.

Akerib, D. S. [SLAC National Accelerator Laborator↗

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗

Highly Reliable Multi-Port Smart Inverter Modules for PV-Based Energy Systems

Generally, large electrolytic capacitors (E-caps) are required to decouple the double-line frequency power fluctuation inherent in all single-phase dc-ac or ac-dc converters. These E-caps pose the weak link in the photovoltaic (PV) inverter because they have poor reliability as a result of short life expectancy and higher equivalent series resistance (ESR) which leads to lower system efficiency. In this paper, multiport smart dual-inverter modules are proposed for residential PV inverter systems with balanced outputs to eliminate the requirement of large decoupling capacitors, thus leading to more reliable and highly efficient inverter modules. The balanced dual outputs are controlled in phase quadrature such that the double frequency current is completely eliminated. One phase supports the standalone residential loads, and the other phase interfaces with the grid, which is synchronized to the utility grid. Also, the proposed modular configuration is scalable to accommodate higher power rating by stacking and paralleling multiple modules. Simulation results for the dual-inverter modules rated for 1 kW each are presented to validate the proposed concept. Experimental results from a controller hardware-in-the-loop (CHIL) platform also verify its feasibility.

14 SOLAR ENERGY↗