Search NASA⌕ Search

SEARCH · Search NASA

Results for “separability approach”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Neural network potentials with effective charge separation for non-equilibrium dynamics of ionic solids: a ZnO case study

Developing neural network potentials (NNPs) accurate under non-equilibrium dynamics is challenging, as such systems require extensive sampling beyond equilibrium phases. Here we construct high-fidelity NNPs for zinc oxide (ZnO), a polymorphic ionic solid, using density functional theory (DFT) reference data. To efficiently capture transitional configurations, we combine enhanced-sampling molecular dynamics with empirical potentials, data distillation, and pretraining on short-range atomic energies (A-Train), followed by transfer learning with DFT-relabeled datasets. This hierarchical approach improves transferability across polymorphs and stress states. We further introduce effective charge separation, treating long-range Coulombic terms analytically while short-range residual interactions are learned by the NNP. The optimal effective charges fall in the range 0.5–1.0 q e , consistent with dielectric-screened values derived from formal charges but distinct from Bader estimates. Motivated by this observation, we propose a simple data-driven protocol in which effective charges are optimized by comparing DFT reference energies with explicit Coulomb calculations, without additional NNP training. This strategy improves accuracy and transferability in DFT-level predictions of energies, forces, and stress. Together, these results provide a practical charge-selection framework for robust NNP development in ionic solids, enabling reliable simulation of polymorphic phase transformations and non-equilibrium dynamics.

Chemistry↗
DOE OSTI2021

New approach to precisely measure γ-ray intensities for long-lived fission products, with results for the decay of 95 Zr

For many fission products, the rays emitted following decay provide an easily-detectable signature that can be used to identify their quantities and distributions in a sample. As a result, γ-ray spectroscopy is often exploited to study fission-product yields, provided sufficiently accurate information on the γ-ray intensity is available. However, in many cases, the uncertainties in the existing nuclear data are large enough that they compromise the precision achievable for modern experiments and applications. To address this need, we have developed a new experimental method that is well suited to precisely measure absolute γ-ray intensities in the β decay of long-lived fission products. The approach involves the production of a radiopure sample by implantation of a mass-separated ion beam from the CAlifornium Rare Isotope Breeder Upgrade (CARIBU) facility on a thin carbon foil. The emitted β-decay radiation is detected with a 4π gas proportional counter and a meticulously efficiency-calibrated high-purity germanium (HPGe) detector. As a first measurement to demonstrate the approach, we studied the absolute γ-ray intensities of the strongest transitions following the β decay of 95 Zr and its decay-daughter 95 Nb, and determined them to fractional precisions of better than 1–2%. In addition, with a larger sample of activity produced through neutron irradiation of an isotopically-enriched Zr foil, we performed a high-precision measurement of the relative γ-ray intensities following the decay of 95 Zr with just the HPGe detector. Here, the sample-production method at CARIBU and the coincidence detection approach demonstrated here can be applied to study fission products with half-lives longer than a day, which includes isotopes important not only for nuclear-energy and national-security applications, but also for medical-isotope research and environmental monitoring.

4π↗
DOE OSTI2022

Development of Carbon Molecular Sieves Hollow Fiber Membranes based on Polybenzimidazole Doped with Polyprotic Acids with Superior H 2 /CO 2 Separation Properties (Final Report)

The goal of this project was to develop a highly efficient membrane-based process to capture CO 2 from coal-derived syngas with 95% CO 2 purity, achieving the cost of electricity (COE) 30% below the baseline capture approaches (i.e., Selexol process) when coupled with the advancement in other areas of the power generation facility. Our core approach is based on high-permeance hollow fiber membranes (HFMs) with superior H 2 /CO 2 separation properties at the syngas process conditions, which can then be further utilized to design membrane reactors for process intensification of H 2 production and purification and CO 2 capture. Three organizations with complementary skills collaborated to achieve the goal, including the University at Buffalo (UB), Los Alamos National Laboratory (LANL), and Trimeric Corporation (Trimeric). We formulated logical steps to bring the membrane technology from Technology Readiness Level (TRL) 3 (Experimental proof of concept) to TRL 4 (Laboratory scale validation in relevant environment). During the budget period (BP) 1, we screened various polymeric materials and identified polybenzimidazole doped with inorganic polyprotic acids as the desirable platform. The acid doping increases the H 2 /CO 2 selectivity, and the sequential carbonization increases H 2 permeability while retaining the high selectivity. By manipulating the acid type and doping level and the carbonization temperature, we demonstrated advanced carbon molecular sieving (CMS) materials with H 2 permeability of above 200 Barrer (1 Barrer = 10 -10 cm 3 (STP) cm cm -2 s -2 cmHg -1 ) and H 2 /CO 2 selectivity of above 40 at 200-300°C with simulated syngas containing CO and water vapor. For example, the PBI-(H 3 PO 4 ) 0.11 carbonized at 700 °C exhibits H 2 permeability of 200 Barrer and H 2 /CO 2 selectivity of 60 at 200 °C, and H 2 permeability of 240 Barrer and H 2 /CO 2 selectivity of 54 at 225 °C, which meets the targeted properties and surpasses Robeson’s upper bound. During the BP2, we focused on the conversion of the advanced CMS materials to stable HFMs. Membranes with H 2 permeance of 1,090 GPU (1 GPU = 10 -6 cm 3 (STP) cm -2 s -2 cmHg -1 ) and H 2 /CO 2 selectivity of 57 at 300 °C were successfully fabricated. The effects of temperature, gas compositions, pressure, and time on the separation properties were systematically investigated. Pencil modules were continuously evaluated for 219 h (dry pure gas) and 669 h (dry simulated syngas) and showed initial decline in permeance and increased H 2 /CO 2 selectivity over time, ultimately achieving a steady state stable value, indicating that the ageing phenomena in the nanoporous structures of the membranes during the test. Membrane reactors were fabricated based on the CMS membranes and evaluated for water-gas shift (WGS) reaction. The use of membranes slightly improves the conversion of the CO. However, parametric tests of the membrane reactors at various temperatures and flow rates need to be conducted, as well as the membranes with improved separation performance. We performed a sensitivity analysis on the impact of H 2 /CO 2 selectivity on the COE based on a hybrid process of a membrane unit and cryogenic unit developed by Membrane Technology and Research, Inc. (MTR). Three H 2 /CO 2 selectivity (40, 60, and 15) cases were developed and compared with the baseline capture process (Case B5B) provided by the DOE report. Increasing the membrane H 2 /CO 2 selectivity reduces COE, but the rate of the decrease of COE also diminishes. The COE values for H 2 /CO 2 selectivities of 40 and 60 were nearly the same. As the H 2 /CO 2 selectivity increases, the inert recycling decreases, leading to smaller equipment, less auxiliary power requirements, and less heating, cooling, and refrigeration duty. The refrigeration system used to liquefy the CO 2 is the most expensive piece of equipment and consumes the most electricity within the CO 2 capture process. Increasing the CO 2 concentration in the recycle stream would improve the economics of the process by reducing the refrigeration duty requirement of the unit and also allow for higher liquefaction temperatures. The high H 2 -selective membrane developed by our team may be applicable in other separation processes where lower pressure H 2 retains value. Typically, hydrogen retains its pressure when it is separated from syngas components. Residual components may be used as low-quality fuel and then vented to the atmosphere. Applications might include control of H 2 /CO ratios or mitigation of the water gas shift reaction by CO 2 recycling to the feed of a gasifier or steam methane reformer. To summarize, we have developed industrial HFMs with the best H 2 /CO 2 separation performance reported in the literature. The membranes demonstrate stability with simulated syngas and show great potential for membrane reactors for WGS reactions, lowering the cost of blue H 2 production.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

REBOUND: Reverse Engineering Bidirectional Outflow Under Non-Equilibrium Diffusion

Rare-earth elements (REEs) are essential for electronics, renewable energy, and defense technologies. However, the current supply of REEs relies on mining concentrated in a few countries and energy-intensive separations. DOE’s Basic Energy Sciences (BES) program has launched a grand challenge which aims to ensure a sustainable supply of critical REEs by developing innovative and environmentally friendly separation methods. As an alternative to costly and harmful traditional methods, the Non-Equilibrium Transport Driven Separations (NETS) initiative has created a microfluidic Y-channel co-flow method that applies external fields to exploit magneto- and electrohydrodynamic effects for separating dilute REE ions from complex feedstocks. Computational fluid dynamics (CFD) studies have identified a few operating conditions with promising ion selectivity and separation efficiency. However, challenges remain regarding Y-channel versatility across feedstocks and accurate incorporation of physical phenomena into CFD models. In this work, we develop a multi-fidelity modelling approach which integrates experimental results with CFD simulation to build a surrogate model for the dependence of separation efficiency to variation of design parameters. The surrogate model enables a reinforcement learning (RL) method to adaptively launch CFD and experimental runs, improving model fidelity around optimal Y-channel parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Entropy–Preserving and Entropy–Stable Relaxation IMEX and Multirate Time–Stepping Methods

In this work, we propose entropy-preserving and entropy-stable partitioned Runge–Kutta (RK) methods. In particular, we extend the explicit relaxation Runge–Kutta methods to IMEX–RK methods and a class of explicit second-order multirate methods for stiff problems arising from scale-separable or grid-induced stiffness in a system. The proposed approaches not only mitigate system stiffness but also fully support entropy-preserving and entropy-stability properties at a discrete level. The key idea of the relaxation approach is to adjust the step completion with a relaxation parameter so that the time-adjusted solution satisfies the entropy condition at a discrete level. The relaxation parameter is computed by solving a scalar nonlinear equation at each timestep in general; however, as for a quadratic entropy function, we theoretically derive the explicit form of the relaxation parameter and numerically confirm that the relaxation parameter works the Burgers equation. Several numerical results for ordinary differential equations and the Burgers equation are presented to demonstrate the entropy-conserving/stable behavior of these methods. We also compare the relaxation approach and the incremental direction technique for the Burgers equation with and without a limiter in the presence of shocks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Single-size and cluster dynamics modeling of intra-granular fission gas bubbles in UO 2

For this work, we perform simulations of intra-granular fission gas bubble evolution in UO 2 using both a relatively simple, computationally inexpensive single-size model and a detailed cluster dynamics model. Simulations encompass 36 experimental cases from 4 different databases, covering various temperature and burnup levels. We systematically compare results from the two models to each other and to post-irradiation experimental data of bubble average size and number density. Overall, the model-to-model comparisons reveal an excellent agreement across the set of simulations. This outcome indicates that, in spite of the underlying assumptions, the single-size model provides a good approximation of the complex physical behavior that is more rigorously described by the cluster dynamics model. Qualitatively, both models reproduce the trends of the experimental data with temperature and burnup correctly. Quantitatively, calculated results are either in good agreement with the data or within errors that appear consistent with the inherent uncertainties. Moreover, for the single-size model, we demonstrate and assess a multiscale approach whereby values for the fission gas atom diffusion coefficient from separate atomistic calculations are used. Systematic comparisons to experimental data point out a credible accuracy of the multiscale model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Thermal equilibrium of collisional non-neutral plasma in a magnetic dipole trap

This paper discusses thermal equilibrium states of single-species plasmas, such as pure electron plasmas and pure positron plasmas, that are confined in a dipole trap. Thermal equilibrium states for such plasmas are routinely realized in the homogeneous magnetic field of Penning–Malmberg traps. We generalize the theory of these states to include inhomogeneous magnetic dipole fields. The approach to thermal equilibrium takes place in two stages with well separated time scales. On the collision time scale, thermal equilibrium is established along each magnetic field line. On the much longer transport time scale, heat conduction and viscosity bring the plasmas on different flux contours into thermal equilibrium, we call this a state of global thermal equilibrium. We present numerical results for local and global thermal equilibria. These results agree with the analytic predictions for charge collections that are large compared with the Debye length. There is, in principle, no limit to the confinement time of a single-species plasma in a global thermal equilibrium state. Experiments with hot electron–ion plasmas performed in the LDX and RT1 devices give us confidence that, in contrast to a Penning–Malmberg trap, a magnetic dipole field can also confine cold quasi-neutral electron–positron pair plasmas on the time scale of the phenomena of interest. Such pair plasmas are assumed to form in the magnetosphere of neutron stars but have so far not been realized in a laboratory. The creation of an electron–positron pair plasma is the main goal of the APEX collaboration.

Steinbrunner, P. (ORCID:0000000317693686)↗
DOE OSTI2026

Multi-Spectroscopic Determination of Exchange Coupling, Zero-Field Splitting, and g-Matrices in Radical-Bridged Dinuclear Fe(III) Complexes

When the energy gap, Δ, between the lowest-lying spin manifolds within a spin-exchange coupled molecule approaches Δ/k B ≈ 300 K, the traditional temperature-dependence (T < 400 K) of the molar magnetic susceptibility is not always a reliable way to obtain a good estimate of intramolecular exchange couplings. We develop a spectroscopic approach capable of accurately parametrizing complex magnetic Hamiltonians by exploiting the separation of the anisotropy and exchange energy scales in strongly coupled magnetic molecules. Specifically, we combine inelastic neutron scattering, high-frequency electron paramagnetic resonance, far-infrared magneto-spectroscopy and magnetometry, and obtain detailed information about the magnetic properties of a series of diiron complexes derived from [[Fe(cth)] 2 (dxbq)] 3+ (H 2 dxbq: 2,5-dihydroxy-1,4-benzoquinone (x = h) or 3,6-dichloro-2,5-dihydroxy-1,4-benzoquinone (x = c), cth: 5,5,7,12,12,14-hexamethyl-1,4,8,11-tetraazacyclotetradecane). Well-isolated S = 9/2 ground states emerge due to strong direct antiferromagnetic exchange between the Fe 3+ centers (S = 5/2) and the radical bridging benzoquinone ligand (S = 1/2). The specific sensitivities and transition selection rules of the applied methods allow us to determine the parameters of the microscopic Hamiltonian including exchange coupling, fourth-order Stevens operators and g-factors. Our methodology is directly portable to other strongly coupled molecular compounds.

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

Field-Driven Simulations to Probe the Impact of Ionic Correlations on Solution Transport Coefficients in Binary, Ternary, and Reciprocal Quaternary Aqueous Electrolytes

Ionic correlations play a critical role in governing transport properties of mixed-salt aqueous electrolyte solutions, yet their quantitative characterization remains challenging, particularly for multicomponent aqueous systems. Here, we develop a general nonequilibrium molecular dynamics framework to efficiently compute Onsager transport coefficients and ionic correlations in mixed-salt aqueous solutions. Using field-driven simulations, we obtain accurate Onsager matrices for LiCl/ KCl, KCl/KBr, and LiBr/KCl electrolyte solutions with significantly reduced computational cost relative to equilibrium Green−Kubo methods. The framework enables direct assessment of how attractive cation−anion and repulsive like-ion interactions contribute to conductivity and salt diffusivity across compositions. While static ion pairing exhibits strong composition dependence, dynamic ion correlations remain nearly invariant, leading to constant deviations from Nernst−Einstein predictions. These results highlight the disconnect between static ion association and dynamic transport correlations, and they establish a transferable approach for analyzing ion transport in complex electrolyte environments relevant to separation processes and electrochemical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Water-Fed, Photovoltaic-Driven Anion-Exchange Membrane Water Electrolysis for Solar Hydrogen Production

Direct coupled photovoltaic (PV)-electrolysis is a promising approach for low-cost hydrogen production, avoiding the need for separate electricity generation. Without solar concentration, the operating current densities of a PV-electrolyzer are small, and low capital costs are needed to reach hydrogen production cost targets (<$2/kg). Anion-exchange membrane (AEM) electrolyzers could be well-suited for this application due to their ability to use platinum group metal (PGM)-free catalysts and operate without supporting electrolytes, but a water-fed PV-AEM system has not yet been demonstrated. In this work, the performance of two AEM electrolyzer designs under pure-water, low-temperature, and diurnal-cycling conditions was evaluated. A simple PV-electrolyzer system design with direct electric coupling to a commercial 84 cm2 Si mini module and passive heating and water flow to the electrolyzer was used for on-sun testing in October 2025. The best-performing PV-AEM system achieved an average solar-to-hydrogen (STH) efficiency of 6.6% and a production rate of 15 mg/kWh/m2 of solar irradiance over 11 days. Minimal electrolyzer corrosion was observed, with no loss in efficiency over the diurnal cycles. While highlighting areas for improved electrolyzer and system design, this work is a proof of concept for distributed hydrogen production using inexpensive and abundant materials.

08 HYDROGEN↗
DOE OSTI2022

Predicting solid state material platforms for quantum technologies

Semiconductor materials provide a compelling platform for quantum technologies (QT). However, identifying promising material hosts among the plethora of candidates is a major challenge. Therefore, we have developed a framework for the automated discovery of semiconductor platforms for QT using material informatics and machine learning methods. Different approaches were implemented to label data for training the supervised machine learning (ML) algorithms logistic regression, decision trees, random forests and gradient boosting. We find that an empirical approach relying exclusively on findings from the literature yields a clear separation between predicted suitable and unsuitable candidates. In contrast to expectations from the literature focusing on band gap and ionic character as important properties for QT compatibility, the ML methods highlight features related to symmetry and crystal structure, including bond length, orientation and radial distribution, as influential when predicting a material as suitable for QT.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Impact of beam coupling impedance on crab cavity noise induced emittance growth

Crab cavities will be deployed as a part of the High Luminosity Large Hadron Collider (HL-LHC) upgrade to mitigate the luminosity reduction induced by the crossing angle at the main experiments (ATLAS and CMS). Two prototype crab cavities have been installed in the CERN Super Proton Synchrotron (SPS) in 2018 for studies with proton beams. An issue of concern is the transverse emittance growth induced by noise in the crab cavity radio frequency (rf) system, which is anticipated to limit the performance of the HL-LHC. In measurements conducted in the SPS in 2018, the crab cavity noise-induced emittance growth was measured to be a factor of 4 lower than predicted from the existing analytical models. In this paper, it is shown that the observed discrepancy is explained by damping effects from the beam coupling impedance, which were not included in the models up to now. Using the van Kampen mode approach, a new theory is developed, suggesting that the impedance can separate the coherent tune from the incoherent spectrum leading to an effective reduction of the crab cavity rf noise-induced emittance growth. This mechanism is validated in tracking simulations using the SPS impedance model as well as in dedicated experimental measurements conducted in the SPS in 2022. The implications for the HL-LHC project are discussed. Published by the American Physical Society 2024

Physics↗
DOE OSTI2024

A General Framework for Error-controlled Unstructured Scientific Data Compression

Data compression plays a key role in reducing storage and I/O costs. Traditional lossy methods primarily target data on rectilinear grids and cannot leverage the spatial coherence in unstructured mesh data, leading to suboptimal compression ratios. We present a multi-component, error-bounded compression framework designed to enhance the compression of floating-point unstructured mesh data, which is common in scientific applications. Our approach involves interpolating mesh data onto a rectilinear grid and then separately compressing the grid interpolation and the interpolation residuals. This method is general, independent of mesh types and typologies, and can be seamlessly integrated with existing lossy compressors for improved performance. We evaluated our framework across twelve variables from two synthetic datasets and two real-world simulation datasets. The results indicate that the multi-component framework consistently outperforms state-of-the-art lossy compressors on unstructured data, achieving, on average, a 2.3 − 3.5× improvement in compression ratios, with error bounds ranging from 1 × 10 the −6 to 1×10−2. We further investigate impact of hyperparameters, such as grid spacing and error allocation, to deliver optimal compression ratios in diverse datasets.

Gong, Qian↗
DOE OSTI2024

Streamlined Loads Analysis of Floating Wind Turbines With Fiber Rope Mooring Lines

This paper presents an approach for more accurate yet relatively streamlined accounting for the nonlinear characteristics of synthetic fiber rope mooring lines for floating wind turbines. First, we select a minimal set of parameters that can efficiently approximate a mooring line material's quasi-static and dynamic stiffness characteristics. We also develop a set of baseline coefficients for different rope materials based on published product information and research papers. We then expand a quasi-static mooring model to include dynamic mooring line stiffness terms in a way that allows the nonlinear stiffness behavior of fiber ropes to be considered even in quasi-static analyses. For dynamic analysis, we have updated the model MoorDyn, coupled with OpenFAST, to work with the new dynamic mooring line stiffness terms. This includes the implementation of a new viscoelastic approach that allows the tension-strain relation of each mooring line segment to vary between two stiffness values, depending on the strain rate. After presenting the formulation of the approaches for modeling synthetic ropes, they are demonstrated on a floating wind turbine mooring system with taut polyester rope mooring lines. The results are compared with those of other approaches of similar fidelity, including the static-dynamic method and separate simulations with static and dynamic stiffness values. Comparing the results shows the ability of the new method to match the results of the previous methods in a more streamlined manner.

elasticity↗
DOE OSTI2020

Critical Materials Capabilities at LANL [Slides]

Critical materials are a recognized problem. In addition to being it’s own cross-cutting topic area, it is called out in Advanced Energy Storage Initiative, Transportation Sector Priorities, and Energy Efficiency Sector Priorities. The need for domestic battery technology is a priority. Domestic supply, separations and processing technologies are required to reduce dependence on foreign capabilities. LANL maintains many capabilities that are applicable to REEs and critical materials: Actinide processing capability for defense programs and extensive separation capabilities – trace analysis up to pilot scale. Development of new approaches for reprocessing technologies are often tested first on lanthanides.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Developing a sustainability science approach for water systems

We convened a workshop to enable scientists who study water systems from both social science and physical science perspectives to develop a shared language. This shared language is necessary to bridge a divide between these disciplines’ different conceptual frameworks. As a result of this workshop, we argue that we should view socio-hydrological systems as structurally co-constituted of social, engineered, and natural elements and study the “characteristic management challenges” that emerge from this structure and reoccur across time, space, and socioeconomic contexts. This approach is in contrast to theories that view these systems as separately conceptualized natural and social domains connected by bi-directional feedbacks, as is prevalent in much of the water systems research arising from the physical sciences. A focus on emergent characteristic management challenges encourages us to go beyond searching for evidence of feedbacks and instead ask questions such as: What types of innovations have successfully been used to address these challenges? What structural components of the system affect its resilience to hydrological events and through what mechanisms? Are there differences between successful and unsuccessful strategies to solve one of the characteristic management challenges? If so, how are these differences affected by institutional structure and ecological and economic contexts? To answer these questions, social processes must now take center stage in the study and practice of water management. We also argue that water systems are an important class of coupled systems with relevance for sustainability science because they are particularly amenable to the kinds of systematic comparisons that allow knowledge to accumulate. Indeed, the characteristic management challenges we identify are few in number and recur over most of human history and in most geographical locations. This recurrence should allow us to accumulate knowledge to answer the above questions by studying the long historical record of institutional innovations to manage water systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Streamlined Loads Analysis of Floating Wind Turbines With Fiber Rope Mooring Lines: Preprint

This paper presents an approach for more accurate yet relatively streamlined accounting for the nonlinear characteristics of synthetic fiber rope mooring lines for floating wind turbines. First, we select a minimal set of parameters that can efficiently approximate a mooring line material's quasi-static and dynamic stiffness characteristics. We also develop a set of baseline coefficients for different rope materials based on published product information and research papers. We then expand a quasi-static mooring model to include dynamic mooring line stiffness terms in a way that allows the nonlinear stiffness behavior of fiber ropes to be considered even in quasi-static analyses. For dynamic analysis, we have updated the model MoorDyn, coupled with OpenFAST, to work with the new dynamic mooring line stiffness terms. This includes the implementation of a new viscoelastic approach that allows the tension-strain relation of each mooring line segment to vary between two stiffness values depending on the strain rate. After presenting the formulation of the approaches for modeling synthetic ropes, they are demonstrated on a floating wind turbine mooring system with taut polyester rope mooring lines. The results are compared with those of other approaches of similar fidelity including the static-dynamic method and separate simulations with static and dynamic stiffness values. Comparing the results shows the ability of the new method to match the results of the previous methods in a more streamlined manner.

elasticity↗
DOE OSTI2021

An Orthogonal Recursive Bisection (ORB) Based Time Advancement Algorithm for CFD-DEM Solvers

The time integration of the granular phase in coupled computational fluid dynamics (CFD) – discrete element method (DEM) simulations presents a unique computational challenge brought about by the large variations in particle collisional time scales. Particles in the dilute regions of the computational domain can be advanced with large time steps while dense regions require much smaller time increments. However, the time step size in most solvers is globally set as the limit for accuracy and stability imposed by the collisions and is typically orders of magnitude less than that required away from collisions. This work addresses this precise issue and provides a strategy to avoid the use of a global conservative small time step size for the entire set of particles.A novel time stepping algorithm for CFD-DEM solvers using a partitioning approach using orthogonal recursive bisection (ORB) that allows for variable time steps among particles is described and its computational performance is compared against baseline explicit methods, typically used in several CFD-DEM solvers. ORB has advantages of being relatively quick and easy to update incrementally and has the required heuristic behavior (i.e., it will split the region in half with a cluster on each side) when groups of particles are well separated (clustered). The algorithm presented in this work uses a local time stepping approach to resolve collisional time scales for subsets of particles that are present at the leaves of the ORB, thereby resulting in substantial reduction of computational cost. The parallel implementation of this method where a ``knapsack” algorithm is used in tandem with ORB for effective load-balancing is also presented, where a best possible partitioning is obtained based on number of particles and local time-stepping costs. The algorithm is tested against benchmark problems with varying particle distributions that include fluidized bed and riser flow scenarios. Preliminary results indicate that the approach is 2-3X faster than traditional explicit methods for problems that involve both dense and dilute regions, while maintaining the same level of accuracy.

adaptive timestepping↗