Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Simulation”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

C3MechLite: An integrated component library of compact kinetic mechanisms for low-carbon, carbon neutral and zero-carbon fuels

Based on our latest detailed chemical reaction mechanism, C3MechV4.0, we have developed two reduced reaction mechanisms—C3MechLite and C3MechCore—targeting C 0 –C 3 chemical species including NH 3 . C3MechLite (61 species), contains a number of species comparable to GRI-Mech (53 species), that can accurately predict the combustion characteristics of hydrogen, carbon monoxide, ammonia, methane, natural gas, nitrogen oxides, and their mixtures for a wide range of conditions. C3MechCore (118 species) targets a more comprehensive range of C 0 –C 3 fuels, including ammonia, methanol, ethanol, and dimethyl ether. Both mechanisms demonstrate predictive accuracy comparable to C3MechV4.0 for the combustion characteristics of the target fuels. C3MechLite is designed with a component library structure, enabling further reduction in mechanism size depending on the fuel(s) of interest for 2D/3D numerical simulations. Various combinations of component libraries were validated, and the average prediction error remains within 1 % compared to C3MechLite. Furthermore, the mechanism was applied to 3D LES simulations of H 2 lifted flames and was confirmed to reproduce flame characteristics with high accuracy. C3MechLite and its component library structure enable high-fidelity and computationally efficient chemical kinetic mechanisms, paving the way for application in more complex combustion simulations.

Ammonia↗
DOE OSTI2025

Turbulent burning velocity of lean premixed hydrogen/air flames at engine conditions: Effects of turbulence intensity and length scale

For turbulent lean premixed hydrogen flames with strong thermodiffusively instabilities, most previous studies have focused on the influence of turbulence intensity, whereas the role of turbulence length scale is less well understood. Here, this study addresses this gap by conducting direct numerical simulations (DNS) of statistically planar turbulent premixed flames for a lean (ϕ=0.35) hydrogen/air mixture under independently varied turbulence intensity (u') and length scale (l T ) at engine-relevant thermodynamics conditions. Results show that as u' increases, the flame front becomes increasingly wrinkled, forming smaller cellular structures. In contrast, l T variations do not significantly alter the size of these structures. For the turbulent burning velocity (S T ), the normalized S T (i.e., S T /S L , where S L is the laminar flame speed) increases linearly with u', driven by both enhanced flame surface wrinkling (i.e., increased A T /A L ) and enhanced local burning rate (i.e., increased I 0 ). However, increasing l T reduces I 0 , despite a continued increase in A T /A L , resulting in only a marginal increase in S T /S L . To reveal the underlying mechanisms, especially the decreasing trend of I 0 with l T , local flame dynamics analyses are performed. It is found that as l T increases, the interaction between thermodiffusive effects and turbulence weakens due to the reduced tangential strain rate, while the flame curvature remains largely unchanged. This suppresses local reactivity enhancement and thus decreases I 0 , In contrast, an increase in u' enhances the interaction by amplifying both curvature fluctuation and tangential strain rate, leading to increased local reactivity (increased I 0 ). Finally, based on the DNS data, several new scaling models are proposed for the three global properties, S T /S L , A T /A L , and I 0 , and showed improvements compared to existing models. These findings provide new insights into the flame-turbulence interactions in thermodiffusively unstable hydrogen flames. The DNS dataset is also useful for the development of turbulent combustion models applicable to practical engine simulations.

Engine-relevant condition↗
DOE OSTI2024

Nanocrystal programmable assembly beyond hard spheres (or shapes) and other (simple) potentials

Ligands are the key to almost any strategy in the assembly of programmable nanocrystals (or nanoparticles) and must be accurately considered in any predictive model. Hard Spheres (or Shapes) provide the simplest and yet quite successful approach to assembly, with remarkable sophisticated predictions verified in experiments. There are, however, many situations where hard spheres/shapes predictions fail. This prompts three important questions: In what situations should hard spheres/shapes models be expected to work? and when they do not work, Is there a general model that successfully corrects hard sphere/shape predictions? and given other successful models where ligands are included explicitly, and of course, numerical simulations, can we unify hard sphere/shape models, explicit ligand models and all atom simulations? . Further, the Orbifold Topological Model (OTM) provides a positive answer to these three questions. In this paper, I give a detailed review of OTM, describing the concept of ligand vortices and how it leads to spontaneous valence and nanoparticle “eigenshapes” while providing a prediction of the lattice structure, without fitting parameters, which accounts for many body effects not captured by (two-body) potentials. I present a thorough survey of experiments and simulations and show that, to this date, they are in full agreement with the OTM predictions. I conclude with a discussion on whether NC superlattices are equilibrium structures and some significant challenges in structure prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Comparative analysis of thermal management systems in electric vehicles at extreme weather conditions: Case study on Nissan Leaf 2019 Plus, Chevrolet Bolt 2020 and Tesla Model 3 2020

With the surge in electric vehicle (EV) adoption and the need for extended driving ranges, optimizing energy efficiency, particularly through thermal management, is critical, especially in extreme weather. Managing the substantial energy needed for cabin climate control and battery temperature regulation can increase energy demands by over 50 %, severely limiting range. This study conducts a comparative analysis of thermal management systems (TMS) in three popular EV vehicles, 2020 Chevrolet Bolt, 2019 Nissan Leaf Plus, and 2020 Tesla Model 3, evaluating their distinct TMS configurations and performance under varied weather conditions. Using both numerical simulations and experimental data collected on a controlled test bench at Argonne National Laboratory, we assess how TMS architecture and operational modes influence energy consumption and range. A comprehensive TMS model was developed, integrating cabin and battery thermal sub-models in the Autonomie software platform, to simulate temperature fluctuations and range impacts. Cabin climate was modeled using a mono-zonal approach, while battery cell temperature distribution was estimated through a 2D nodal structure. Each vehicle's distinct TMS setup was evaluated: the Chevrolet Bolt and Tesla Model 3 use a dual evaporator vapor compression cycle with a PTC heater for the cabin and a coolant loop for battery thermal management; the Nissan Leaf Plus employs a heat pump with a PTC heater for the cabin and air-cooling for the battery. Tests conducted at ambient temperatures of 35°C, 22°C, -7°C, and -18°C reveal significant differences in energy use and range reduction across both configurations and conditions. At 35°C, the Tesla Model 3, Chevrolet Bolt, and Nissan Leaf Plus have a range reduction of 8%, 9%, and 13%, respectively, due to air conditioning. In winter, heating technology is paramount; at -7°C, the Nissan Leaf's heat pump configuration achieves a lower range reduction (19.3%) compared to the Tesla and Chevrolet Bolt PTC heaters, which reduce range by 28.3% and 31%, respectively. Further, this study provides valuable insights for automotive engineers, EV technology researchers, and thermal management system designers aiming to enhance electric vehicle performance by understanding how different weather conditions and TMS architectures impact energy consumption and driving range.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

An update to the Sandia method for creating Typical Meteorological Years from a limited pool of calendar years

Typical Meteorological Years (TMYs) are essential for the efficient evaluation of energy system performance. Ideally, 30 years of weather data are required to generate TMYs, but significantly fewer years are typically available due to practical limitations. To address this issue, an update to the Sandia method was developed, referred to as the Argonne method, to create TMYs from a limited number of years. Furthermore, this method enhances candidate diversity by systematically shifting original candidate months forward or backward by specific days, creating an expanded pool of candidates. The effectiveness of the Argonne method was validated through statistical testing, comparison of monthly average weather parameters, and numerical simulations. The results demonstrate a high probability of identifying at least one shifted month whose cumulative distribution functions of weather parameters closely align with long-term distributions. In 67 % of all comparisons, the monthly average weather parameters in TMYs generated using the Argonne method exhibit better agreement with long-term averages than TMY3. Moreover, in 74 % of the 318 building simulation cases, the Argonne method outperforms TMY3 in estimating long-term average building heating and cooling demands. Therefore, the Argonne method effectively diversifies the candidate pool and produces typical years that provide more accurate estimations of long-term averages compared to TMY3 when only a limited pool of calendar years (10 years or fewer) is available.

Building energy modeling↗
DOE OSTI2025

Computationally efficient models for aqueous organic redox flow batteries

The rising usage of intermittent energy has garnered the need for large scale energy storage systems. Redox flow batteries (RFB) based energy storage system shows promising potential. Numerical simulations and machine learning approaches have been widely used to study RFB performance. The development of autonomous material discovery framework and digital twin of energy storage system usually needs to query cell performance through fast response models. In this study, two computationally efficient models are introduced: a physics-based analytical flow battery model (EZBattery), and a machine learning operator model (Deep Operator Network, denoted by DeepONet). Both models can provide cell performance near instantly, and prediction accuracy was systematically examined on an application of evaluating the performances of a 780 cm 2 aqueous organic redox flow battery (AORFB), using potential anolyte candidates in dihydroxyphenazine (DHP)-based family of organic materials. A validated computationally expansive 3-dimensional multi-physics finite element model by COMSOL was used as the ground truth and provided the training data set for the DeepONet. 1280 samples were generated with 10 properties to mimic the different possible anolyte candidates, and the cell performances were evaluated under 10 different combined operating conditions. The accuracy comparisons for the two computationally efficient models show that both models can provide comparable accuracy in predicting cell charging/discharging voltage curves. DeepONet can provide slightly higher overall accuracy than EZBattery with faster calculation speed, but highly relies on the training dataset. EZBattery does not need a training dataset and can provide interpretable physics-based explanations of the results, while being more flexible to adjust to adapt any different cell designs, flow battery architectures, and electrolyte materials.

Analytical model↗
DOE OSTI2026

Examination of entrainment, combustibility, and heat transfer due to cell venting inside simplified battery energy storage enclosures

Battery energy storage systems (ESS) assembled of modules containing lithium-ion batteries can pose fire and explosion hazards during thermal runaway because vented gases may accumulate, combust, or escape enclosures. Here, we propose a prototypical configuration composed of the main components of an ESS rack, where modules are represented as rectangles within an enclosure and vent gas is injected steadily through the top of a module. This idealized geometry allows for exploration of how rack-level characteristics can influence smaller single-cell or module-level scales. Through this geometry, we conduct numerical simulations to understand how variations in the vent velocity, temperature, and vertical location affect entrainment rates, fraction of uncombusted fuel, rich-gas coverage, and heat flux to other modules. At lower vent velocities and vertical positions, buoyancy drives enough entrainment and mixing to react with most of the vent gases on the converging plenum. As jet momentum increases and/or buoyancy decreases via higher vent position in the rack, vent gases can propagate into and burn in the divergence plenum. Heat flux to the impinged module scales strongly with vent momentum and temperature, while heat transfer to upper modules peaks when flames in the diverging plenum reach the surface and declines when combustion is suppressed.

Computational fluid dynamics↗
DOE OSTI2025

A new tubular ironless linear induction motor using non-interrupted winding method

The ironless linear induction motor device can be used for any application which requires a large motion span or rotation restrictedly by the ironless feature due to the strong magnetic field environment. One of the possible applications is ITER diagnostics shutters. The tokamak reactor systems have a severe radiation issue during the operation. The traditional Kapton-insulated wires have significantly lower radiation resistance than an expected large radiation dose at ITER Motion Stark Effect (MSE) shutter location. Furthermore, the traditional single-sided linear induction motor winds the wire to make the N-turn pitch coils independently and then connect them in series to make a one-phase winding. Therefore, it is a big challenge to make so many radiation-protection connectors in such a limited space as inside the vacuum vessel of fusion reactors. The approved mineral-insulated kinds of cables need a large bend radius for the end turns, so the traditional single-sided linear motor with two-layer three-phase windings need to be upgraded to the tubular linear motors. The new proposed tubular ironless linear induction motor can not only overcome the large bend radius issue with the mineral-insulated cables but also allow the motor coils wound continuously around the tubular-shaped stator, forward and backward N times to make a N-turn winding for each phase, in order to resolve the radiation protection issue with the winding connection. The 3D transient Maxwell electromagnetic model with a tubular three-phase linear induction motor has been analyzed. Furthermore, this paper will present the detailed design and winding method of ironless linear induction motor actuators, and results of the 3D transient numerical simulation of the actuator performance.

Ironless↗
DOE OSTI2025

NdPO 4 solubility and aqueous Neodymium speciation in supercritical fluids: An experimental study at 500–700 °C and 1.7 kbar

A key aspect in the formation of rare earth elements (REE) deposits is the role of REE transport as aqueous REE complexes in supercritical hydrothermal solutions, where the nature of the aqueous complex is controlled by solution composition, temperature and pressure. Despite chloride being considered as one of the most abundant transporting ligands in magmatic-hydrothermal fluids, experimental investigations on the stability of aqueous REE chloride complexes are scarce above 300 °C. In this study, synthetic NdPO 4 crystals were reacted with non-saline and saline (0, 0.05 and 0.5 mNaCl), acidic (0.01 mHCl) aqueous solutions in a series of solubility experiments conducted at 500–700 °C and 1.7 kbar, where the solubilities were determined using a stable Nd isotope ( 145 Nd isotope spike) dilution technique. NdPO 4 solubility ranges between 28 ppm and 10,858 ppm, where solubility increases with both temperature and salinity. At 500 °C, log mNdPO 4 increases from –3.93 to –1.60 and there is a strong correlation between NdPO 4 solubility and NaCl concentrations (slope of 1.2 ± 0.3), indicating stabilization of the Nd chloride aqueous complexes with a stoichiometry corresponding to NdCl 2+ . At 600 °C, this correlation is weaker (slope of 0.4, log mNdPO 4 increases from –2.63 to –1.88) indicating the stabilization of both Nd chloride and hydroxyl species controlling solubility. At 700 °C, NdPO 4 solubility is largely independent of NaCl concentration indicating that solubility is controlled by Nd hydroxyl complexes, where stoichiometry suggests the neutral Nd(OH) 3 0 species is dominant. The solubility product (Ksp) of NdPO4 is derived from experimental data with the relation: log K sp = -41.81 – 0.057T – 20987/T, with T temperature in Kelvin. Comparison of the measured Nd phosphate solubility to thermodynamic predictions using the available Helgeson-Kirkham-Flowers equation of state parameters for aqueous Nd complexes indicate that predictions are up to three orders of magnitude lower compared to experimental observations. This discrepancy is most pronounced in saline solutions, suggesting that thermodynamic properties of the REE chloride species in supercritical fluids require revision. Numerical simulations of fluid-rock interaction between acidic, saline fluids and a Strange Lake felsic mineral assemblage demonstrates that NdPO 4 solubility predictions from models are four to six orders of magnitude lower than those calculated based on empirical fits from experiments, which suggests that acidic, saline fluids may play an important role in mobilizing large amounts of light REE from 450 to 700 °C.

58 GEOSCIENCES↗
DOE OSTI2024

Experimental and numerical investigation of fracture conductivity between non-smooth rock surfaces with and without proppant

The enhancement of fracture conductivity is vital for the efficient recovery of subsurface resources, such as geothermal energy and petroleum hydrocarbons. Proppants, granular materials injected into hydraulic fractures to maintain their conductivity, have been studied primarily in the context of smooth fractures (i.e., fractures between smooth rock surfaces). However, non-smooth fractures (i.e., fractures between rough rock surfaces) are common in geoenergy reservoirs and thus require further investigations. In this study, we conducted laboratory measurements of fracture conductivity on shale slabs with non-smooth surfaces and carried out numerical simulation using the lattice Boltzmann (LB) method, which aimed to investigate the conductivity of non-smooth fractures with and without proppants placement. When ceramic proppant with an areal concentration of 2 lb/ft 2 was placed in the fracture, the conductivity was enhanced by roughly 3-8 times compared to fractures without proppant. In fractures with proppant, gas-measured conductivity was higher than that measured with water due to proppant embedment caused by water. The experiments demonstrate the advantages of using proppant in fractures, even if the rock surface roughness can provide certain fracture conductivity via the self-propping mechanism. For fractures without proppants, high rock surface roughness is not necessarily favorable for enhancing fracture conductivity because the self-propping mechanism requires shear slip along the fracture surface. If there is no shear slip, high rock surface roughness can cause a detrimental effect on the fracture conductivity due to the interlocking effect. Utilizing advanced experimental equipment and LB modeling, this research explores the interplays between proppant placement, fracture geometry, and stress conditions to develop a comprehensive understanding of the productivity in non-smooth fractures. Further, the outcomes of this investigation indicate the importance of creating fractures with surface roughness during hydraulic fracturing and will contribute to the development of more efficient stimulation techniques for subsurface energy extraction.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Exploring capabilities of Micro-Fabricated 3D-printed capsules for studying effects of material mix on thermonuclear burn

The introduction of a high-precision, 3D-printing method, known as two-photon polymerization (2PP), has created a new means of producing novel targets to study inertial confinement fusion (ICF). Here, this work aims to explore if 2PP fabricated capsules can meet ICF requirements such as size, uniformity, deuteration degree, and ability to maintain gas fill. Two types of capsules with 3D-printed lattices were evaluated on the Omega-60 Laser Facility: (1) carbon-deuterium-oxygen (CDO) lattice with a H 2 gas fill and (2) carbon-hydrogen-oxygen (CHO) lattice with a D 2 gas fill. Experimental results and numerical simulations, which assumed complete mixing of capsule materials, reasonably agreed for both target types. A further study into the effects of capsule preheat or 2PP geometry on material mix is underway. Our work demonstrates that the unique capabilities of customizable, 3D-printed capsules present promising opportunities for further investigation into more sophisticated mix and burn experiments in ICF.

2PP↗
DOE OSTI2024

A Review: Indirect Optimal Control of Wave Energy Converters

Wave energy conversion has been the subject of interest in the past several years. While there are several concepts for converting wave power into electric power, the cost of the electric power harvested from ocean waves remains high. One of the main challenges, though it receives less attention, is the control of the wave energy converter (WEC). This paper presents a treatment for the WEC control problem within the context of optimal control theory. The result is systematic development for an explicit expression for a control that maximizes the harvested energy while meeting operational constraints such as the maximum device stroke and the maximum control force. The control presented here can also be adjusted to meet device design constraints such as a limitation on the amount of reactive power available from the power take-off (PTO) unit; this feature enables a control co-design for the PTO unit. Numerical simulations are presented in this paper for demonstration.

control design↗
DOE OSTI2026

Theoretical model of dielectrophoretic water vapor condensation during multiphase flow of electrically charged water droplets in air

Water droplet growth in air, influenced by vapor diffusion and electrical charging, has significant implications for indoor air quality, cloud formation, and aerosol dynamics. Here, this study presents a comprehensive theoretical model that describes the growth dynamics of electrically charged water droplets under external electric fields, incorporating diffusion, electrostatic, and dielectrophoretic forces to enhance understanding and practical applications in air dehumidification and environmental control. A theoretical model was established to analyze the growth dynamics of the charged water droplets. The approach extended the thermodynamic framework of the Modified Kelvin-Thomson equation by integrating electrostatic and dielectrophoretic effects, enabling a comprehensive evaluation of vapor concentration near the surface of charged droplets. A sensitivity analysis was performed to assess the roles of droplet size, electric field intensity, surface charge, and vapor concentration on condensation rates. The study progressed beyond the classical diffusion-limited growth of neutral droplets, addressing the behavior of charged droplets in electrohydrodynamic environments. The model highlighted the significance of dielectrophoretic drift flux for dipolar vapor molecules, especially in micron-sized droplets ranging from 10 to 100 µm. Findings demonstrated that, for smaller droplets subjected to elevated electrical charge, average growth rates varied from 0.3 to 1.6 µm/s if the air relative humidity increased from 50 to 100 %. Analytical results are corroborated by numerical simulations, which exhibited strong concordance, thereby validating the proposed theoretical approach. This research introduced an innovative and scalable analytical methodology for electric field-enhanced droplet growth, offering valuable insights into the advancement of energy-efficient air dehumidification, water desalination, and atmospheric water harvesting technologies.

42 ENGINEERING↗
DOE OSTI2025

Directed energy deposition of radially functionally graded nickel based superalloys IN625 and IN738: Microstructure and mechanical properties

This study investigates the microstructure and mechanical properties of radially functionally graded (RFG) materials fabricated from nickel-based superalloys IN625 and IN738 using directed energy deposition (DED). RFG cylindrical deposits with radial gradients of IN625-IN738 and IN738-IN625 were successfully fabricated and confirmed to be crack-free through tailored composition grading, scan strategy and process parameters that minimized thermal gradients and solidification velocities. Microstructural analysis revealed distinct variations across the gradient, characterized by differences in dendrite arm spacing, microsegregation, and the formation of distinct γ' precipitates and carbide/boride phases (MC, M 6 C, M 2 C, and MB 2 ) in as-deposited (AD) conditions, which evolved into M 23 C 6 and M 2 B upon subsequent heat treatment (HT). EBSD showed predominant columnar grains with a [100] texture, while microhardness profiles demonstrated consistent variation across the gradients. Mechanical testing highlighted the synergy of the graded structure, combining IN738's strength with IN625's ductility, resulting in enhanced properties such as a yield strength of 1131 MPa, ultimate tensile strength of 1382 MPa, and elongation of 12 % for HT IN625-IN738 grading. Numerical simulations employing CALPHAD thermodynamic calculations and finite element thermal modeling are used to establish phase evolution and microstructural feature correlation with thermal profiles, aiding in the optimization of process-structure-property relationship.

Additive manufacturing↗
DOE OSTI2024

Fundamentals of wildlife dosimetry and lessons learned from a decade of measuring external dose rates in the field

Methods for determining the radiation dose received by exposed biota require major improvements to reduce uncertainties and increase precision. We share our experiences in attempting to quantify external dose rates to free-ranging wildlife using GPS-coupled dosimetry methods. The manuscript is a primer on fundamental concepts in wildlife dosimetry in which the complexities of quantifying dose rates are highlighted, and lessons learned are presented based on research with wild boar and snakes at Fukushima, wolves at Chornobyl, and reindeer in Norway. GPS-coupled dosimeters produced empirical data to which numerical simulations of external dose using computer software were compared. Our data did not support a standing paradigm in risk analyses: Using averaged soil contaminant levels to model external dose rates conservatively overestimate the dose to individuals within a population. Following this paradigm will likely lead to misguided recommendations for risk management. The GPS-dosimetry data also demonstrated the critical importance of how modeled external dose rates are impacted by the scale at which contaminants are mapped. When contaminant mapping scales are coarse even detailed knowledge about each animal’s home range was inadequate to accurately predict external dose rates. Importantly, modeled external dose rates based on a single measurement at a trap site did not correlate to actual dose rates measured on free ranging animals. These findings provide empirical data to support published concerns about inadequate dosimetry in much of the published Chernobyl and Fukushima dose-effects research. Furthermore, our data indicate that a huge portion of that literature should be challenged, and that improper dosimetry remains a significant source of controversy in radiation dose-effect research.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2026

A finite viscoelastic constitutive model for low to high strain rate response of elastomers with application of strain rate-induced glass transition

Amorphous elastomers exhibit significant rate-stiffening and unique viscous flow characteristics across a wide range of strain rates, often undergoing glass transition above a strain rate threshold. We have developed a thermodynamically-consistent and micromechanically-inspired constitutive model for soft elastomeric materials to capture the rate-dependent stress-strain behavior and hysteresis when subjected to low to high strain rates. Here, our proposed constitutive model encapsulates the viscous flow of materials through molecular motion at low strain rates and local rearrangement and alignment of the molecules trying to overcome the intermolecular resistance at high strain rates, essentially covering the glass transition. We applied our constitutive model to uniaxial compression experiments performed at low and high strain rates for polyborosiloxane (PBS) to identify the material parameters, and subsequently, performed numerical simulations of single and multi-cycle compression, stress relaxation, and small amplitude oscillatory tension-compression. Our analyses indicate that the model predicts higher total energy dissipation with increasing strain rate; however, dissipation associated with molecular relaxation decreases (forming a cusp) because, beyond a crossover strain rate, intermolecular rearrangement and alignment become dominant, which is consistent with the onset of the glass transition. For cyclic loading-unloading, we observed that dissipation over a cycle remains constant at low strain rates but decreases non-monotonically at high strain rates before becoming constant, with the peak stress over the cycle becoming higher, which can be interpreted as more loading being carried elastically by the polymer network as the intermolecular rearrangement process occurs. Additionally, our model was able to predict the qualitative nature of the storage modulus and loss modulus in the limit of small strain over a wide range of frequency sweeps.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Angular-spatial hp -adaptivity for radiative transfer with discontinuous Galerkin spectral element methods

Radiative transfer is important for many science and engineering applications, and numerical simulations of radiative transfer can be challenging. For instance, the radiation field is seven-dimensional – three spatial, two angular, one wavelength, and one temporal – and often features steep gradients. Therefore, memory usage is a key issue. To reduce memory, some past work has investigated the use of adaptive mesh refinement (AMR), typically for either the spatial or angular coordinate, and typically for only h -adaptivity. Here, we propose the use of AMR for the spatial and angular coordinates together, and the use of h - and p -adaptivity together as hp -AMR for the potential for further memory savings. We implemented the proposed method for several test cases in two spatial and one angular dimension, with the discontinuous Galerkin spectral element method. These test cases featured highly anisotropic angular radiation, with or without steep spatial gradients. Our primary findings from these test cases were: (1) Angular hp -adaptivity can deliver the radiation solution with the same accuracy as, and with much less computational memory than, uniform angular h - or p -refinements, or angular h -adaptivity alone. This is most obvious when the incoming radiation is highly anisotropic, in which case the savings can be orders of magnitude. (2) Full spatial-angular hp -adaptivity is more efficient in solution representation, compared to solely spatial or solely angular -adaptivity. This is most evident when steep gradients are present in both the spatial and angular distribution. These results suggest that adaptive spatial- hp angular-refinement may perform well in large-scale seven-dimensional applications.

Adaptive refinement↗
DOE OSTI2026

Omega-Limit Sets and Input-To-State Stability in Power Grids with Switching Equilibria

This paper studies a power transmission system with both conventional generators (CGs) and distributed energy assets (DEAs) providing frequency control. We consider an operating condition with demand aggregating two dynamic components: one that switches between different values on a finite set, and one that varies smoothly over time. Such dynamic operating conditions may result from protection scheme activations, external cyber-attacks, or due to the integration of dynamic loads, such as data centers. Mathematically, the dynamics of the resulting system are captured by a system that switches between a finite number of vector fields - or modes -, with each mode having a distinct equilibrium point induced by the demand aggregation. To analyze the stability properties of the resulting switching system, we leverage tools from hybrid dynamic inclusions and the concept of ..omega.. -limit sets from sets. Specifically, we characterize a compact set that is semi-globally practically asymptotically stable under the assumption that the switching frequency and load variation rate are sufficiently slow. For arbitrarily fast variations of the load, we use a level-set argument with multiple Lyapunov functions to establish input-to-state stability of a larger set and with respect to the rate of change of the loads. The theoretical results are illustrated via numerical simulations on the IEEE 39-bus test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗