Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic State Estimation”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quantifying uncertainty in analysis of shockless dynamic compression experiments on platinum. II. Bayesian model calibration

Dynamic shockless compression experiments provide the ability to explore material behavior at extreme pressures but relatively low temperatures. Typically, the data from these types of experiments are interpreted through an analytic method called Lagrangian analysis. Here, in this work, alternative analysis methods are explored using modern statistical methods. Specifically, Bayesian model calibration is applied to a new set of platinum data shocklessly compressed to 570 GPa. Several platinum equation-of-state models are evaluated, including traditional parametric forms as well as a novel non-parametric model concept. The results are compared to those in Paper I obtained by inverse Lagrangian analysis. The comparisons suggest that Bayesian calibration is not only a viable framework for precise quantification of the compression path, but also reveals insights pertaining to trade-offs surrounding model form selection, sensitivities of the relevant experimental uncertainties, and assumptions and limitations within Lagrangian analysis. The non-parametric model method, in particular, is found to give precise unbiased results and is expected to be useful over a wide range of applications. The calibration results in estimates of the platinum principal isentrope over the full range of experimental pressures to a standard error of 1.6%, which extends the results from Paper I while maintaining the high precision required for the platinum pressure standard.

Brown, Justin Lee↗
DOE OSTI2025

Proton Dynamics Scenarios in the Integrable Optics Test Accelerator (IOTA) at Fermilab

The Integrable Optics Test Accelerator (IOTA) at Fermilab provides a versatile platform for studying the interplay of space-charge, impedance, and non-linear optics in high-intensity hadron beams within synchrotrons and storage rings. This report examines the parameters and dynamics of 2.5~MeV proton beam operations in two configurations of the bare IOTA lattice: one for demonstrating Non-linear Integrable Optics with the Danilov-Nagaitsev magnet, and the other for use with electron cooling. We offer order-of-magnitude estimates of the transverse emittance growth rate as a function of beam intensity, highlighting contributions from residual gas scattering, intra-beam scattering, and space-charge effects. Under nominal conditions, the beam lifetime is projected to be less than 7~minutes at low intensity with the current vacuum quality, and fewer than 100,000~turns at high intensity due to strong space-charge effects. The calculations presented here will guide strategies to mitigate emittance growth and inform future IOTA experiments.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

A New Global Storage-Area-Depth Data Set for Modeling Reservoirs in Land Surface and Earth System Models

Reservoir storage-area-depth relationships are the most important factors controlling thermal stratification in reservoirs and, more broadly, the water, energy and biogeochemical dynamics in the reservoirs and subsequently their impacts on downstream rivers. However, most land surface or earth system models do not account for the gradual changes of reservoir surface area and storage with the changing depth, inhibiting a consistent and accurate representation of mass, energy and biogeochemical balances in reservoirs. Here we present a physically coherent parameterization of reservoir storage-area-depth dataset at the global scale. For each reservoir, the storage-area-depth relationships were derived from an optimal geometric shape selected iteratively from five possible regular geometric shapes that minimizes the error of total storage and surface area estimation. We applied this algorithm to 6,800 reservoirs included in the Global Reservoir and Dam database (GRanD). The relative error between the estimated and observed total storage is no more than 5% and 50% for 66% and 99% of all GRanD reservoirs, respectively. More importantly, the storage-depth profiles derived from the approximated reservoir geometry compared well with remote sensing based estimation at 40 major reservoirs from previous studies, and ground-truth measurements for 34 reservoirs in the United States and China. Furthermore, the new global reservoir storage-area-depth dataset is critical for advancing future modeling and understanding of reservoir processes and subsequent effects on the terrestrial hydrological, ecological and biogeochemical cycles at the regional and global scales.

13 HYDRO ENERGY↗
DOE OSTI2021

Electron thermalization and relaxation in laser-heated nickel by few-femtosecond core-level transient absorption spectroscopy

Direct measurements of photoexcited carrier dynamics in nickel are made using few-femtosecond extreme ultraviolet (XUV) transient absorption spectroscopy at the nickel M 2,3 edge. It is observed that the core-level absorption line shape of photoexcited nickel can be described by a Gaussian broadening (σ) and a red shift (ω s ) of the ground-state absorption spectrum. Theory predicts and the experimental results verify that after initial rapid carrier thermalization, the electron temperature increase (ΔT ) is linearly proportional to the Gaussian broadening factor σ, providing quantitative real-time tracking of the relaxation of the electron temperature. Measurements reveal an electron cooling time for 50 nm thick polycrystalline nickel films of 640 ± 80 fs. With hot thermalized carriers, the spectral red shift exhibits a power-law relationship with the change in electron temperature of ω s ∝ ΔT 1.5 . Rapid electron thermalization via carrier-carrier scattering accompanies and follows the nominal 4-fs photoexcitation pulse until the carriers reach a quasithermal equilibrium. Entwined with a <6 fs instrument response function, carrier thermalization times ranging from 34 fs to 13 fs are estimated from experimental data acquired at different pump fluences and it is observed that the electron thermalization time decreases with increasing pump fluence. Here, the study provides an initial example of measuring electron temperature and thermalization in metals in real time with XUV light, and it lays a foundation for further investigation of photoinduced phase transitions and carrier transport in metals with core-level absorption spectroscopy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Sticking Coefficients of Fusion Reactor Impurities from Molecular Dynamics Simulations for the Design of Cryopumps

A cryopump can be utilized as an impurity removal component of a direct internal recirculation (DIR) system for the fusion fuel cycle. The DIR facilitates a low fuel inventory by continuously pumping unburnt fuel while removing impurities from the fusion exhaust stream. A cryopump can target multiple impurity species by maintaining a temperature lower than the gas triple-point temperature that promotes desublimation. The desublimation/condensation of gases in cryopumps can be characterized by the sticking coefficient, which is defined as the probability for a gas particle to stick to a (cryo-)surface upon collision. The sticking coefficient is one of the important design/operation parameters for cryopumps, and it depends on a variety of surface and gas properties. Here, in this study, molecular dynamics simulations were utilized to estimate the sticking coefficients of typical fusion gas impurity species N 2 , CO 2 , and CH 4 over a Cu surface for a range of gas temperatures and surface coverages. The molecular dynamics study showed that the sticking coefficients for gases decrease with an increase in gas temperature. The presence of a single full monolayer of condensate on the metallic surface showed an adverse effect on the sticking of gases; however the sticking improved with two full monolayers of condensate on the surface. The sticking of gases over the mixed condensate on a surface was more favorable than the condensate of the same species for N 2 and CH 4 , with an exception for CO 2 , which showed a decrease in sticking over the mixed condensate.

cryopump↗
DOE OSTI2024

a priori uncertainty quantification of reacting turbulence closure models using Bayesian neural networks

While many physics-based closure model forms have been posited for the sub-filter scale (SFS) in large eddy simulation (LES), vast amounts of data available from direct numerical simulations (DNS) create opportunities to leverage data-driven modeling techniques. Albeit flexible, data-driven models still depend on the dataset and the functional form of the model chosen. Increased adoption of such models requires reliable uncertainty estimates both in the data-informed and out-of-distribution regimes. Here, in this work, we employ Bayesian neural networks (BNNs) to capture both epistemic and aleatoric uncertainties in a reacting flow model. In particular, we model the filtered progress variable scalar dissipation rate which plays a key role in the dynamics of turbulent premixed flames. We demonstrate that BNN models can provide unique insights about the structure of uncertainty of the data-driven closure models. We also propose a method for the incorporation of out-of-distribution information in a BNN, which can be used for out-of-distribution query detection. The efficacy of the model is demonstrated by a priori evaluation on a dataset consisting of a variety of flame conditions and fuels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Critical Role of the Steric Factor in the Viscoelasticity of Vitrimers

Dynamic covalent networks (DCNs) are a promising solution to mitigate plastic-waste-related issues through improved recyclability enabled by dynamic bonds. However, our understanding of the mechanisms controlling their viscoelasticity, especially in vitrimers where bond exchange relies on associative reactions, remains limited. Here, in this study, we investigate the dynamics in model DCNs with boric ester functionalities, and the analysis of the temperature dependence of their terminal relaxation times revealed a puzzling result: extremely large Arrhenius prefactors associated with their dynamic bond rearrangement times. We ascribe this observation to the often-overlooked chemical steric factor that slows down chemical reactions, therefore decreasing the vitrimers’ bond exchange rate by many orders. The estimated steric factor of the bond exchange in our model DCNs is comparable to those observed in boronic ester exchange reactions between small molecules. Additional analysis of literature data revealed an overall low steric factor also for imine bond exchange, thus highlighting the role of this parameter in tremendously slowing down bond exchange in DCNs despite low activation energy barriers. We propose a general approach for designing vitrimers with desired viscoelastic and creep properties considering the critical role of the steric factor in bond rearrangement mechanisms, in addition to the traditionally considered activation energy and matrix properties.

Carden, Gregory P. [Univ. of Tennessee, Knoxville,↗
DOE OSTI2019

The Fire and Smoke Model Evaluation Experiment—A Plan for Integrated, Large Fire–Atmosphere Field Campaigns

The Fire and Smoke Model Evaluation Experiment (FASMEE) is designed to collect integrated observations from large wildland fires and provide evaluation datasets for new models and operational systems. Wildland fire, smoke dispersion, and atmospheric chemistry models have become more sophisticated, and next-generation operational models will require evaluation datasets that are coordinated and comprehensive for their evaluation and advancement. Integrated measurements are required, including ground-based observations of fuels and fire behavior, estimates of fire-emitted heat and emissions fluxes, and observations of near-source micrometeorology, plume properties, smoke dispersion, and atmospheric chemistry. To address these requirements the FASMEE campaign design includes a study plan to guide the suite of required measurements in forested sites representative of many prescribed burning programs in the southeastern United States and increasingly common high-intensity fires in the western United States. Here we provide an overview of the proposed experiment and recommendations for key measurements. The FASMEE study provides a template for additional large-scale experimental campaigns to advance fire science and operational fire and smoke models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

The Arrow of Time in the Collapse of Collisionless Self-gravitating Systems: Non-validity of the Vlasov–Poisson Equation during Violent Relaxation

The collapse of a collisionless self-gravitating system, with the fast achievement of a quasi-stationary state, is driven by violent relaxation, with a typical particle interacting with the time-changing collective potential. It is traditionally assumed that this evolution is governed by the Vlasov–Poisson equation, in which case entropy must be conserved. We run N-body simulations of isolated self-gravitating systems, using three simulation codes, NBODY-6 (direct summation without softening), NBODY-2 (direct summation with softening), and GADGET-2 (tree code with softening), for different numbers of particles and initial conditions. At each snapshot, we estimate the Shannon entropy of the distribution function with three different techniques: Kernel, Nearest Neighbor, and EnBiD. For all simulation codes and estimators, the entropy evolution converges to the same limit as N increases. During violent relaxation, the entropy has a fast increase followed by damping oscillations, indicating that violent relaxation must be described by a kinetic equation other than the Vlasov–Poisson equation, even for N as large as that of astronomical structures. This indicates that violent relaxation cannot be described by a time-reversible equation, shedding some light on the so-called “fundamental paradox of stellar dynamics.” The long-term evolution is well-described by the orbit-averaged Fokker–Planck model, with Coulomb logarithm values in the expected range 10--12. By means of NBODY-2, we also study the dependence of the two-body relaxation timescale on the softening length. The approach presented in the current work can potentially provide a general method for testing any kinetic equation intended to describe the macroscopic evolution of N-body systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Theory of Ion-Mediated Segmental Localization, Activated Structural Relaxation, and the Glass Transition in Polymerized Ionic Liquids

Here we combine polymer integral equation theory for structural correlations with a force-level theory for activated segmental relaxation to construct a microscopic theory for how ions modify the onset of segmental localization, the polymer alpha relaxation time, and glass transition temperature, T g , in dry polymerized ionic liquids (PolyILs). Our most important findings concerning the behavior of T g are (i) a weak dependence on anion–cation Coulomb attraction strength due to strong decoupling of ion and segmental dynamics, (ii) strong reduction with increasing mobile ion size due to both plasticization effects and weakening of ion-modified monomer caging, and (iii) increase with backbone rigidity and intrachain dynamic cooperativity. We quantitatively estimate T g variations as ion size and polymer static and dynamic rigidity are widely varied. The experimental data of 19 PolyILs are organized into a specific pattern consistent with the theoretical predictions with glass transition temperatures varying over a very wide range (~250 K) as observed experimentally. Based on the theoretical results, suggestions are made concerning the search for high conductivity Li- or Na-based PolyILs based on further lowering of T g via new polymer synthesis. Calculations of the alpha time as a function of packing fraction reveal the dynamic fragility increases for larger mobile ions, a more flexible backbone, and/or a smaller degree of intrachain dynamic cooperativity, as a consequence of the increasing importance of collective elasticity in determining the activation barrier. An analysis of transient segmental localization and a crossover to activated dynamics reveals similar, but not identical, trends as predicted at the laboratory vitrification timescale.

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

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Preliminary Look at the LBEG & MPEG Beam Transmissions between HPSim and Operation Data

This report summarizes recent work on estimating the transmission of LANSCE H- beams from the end of the present DTL to both the PSR stripper foil (LBEG) and WNR target 4 (MPEG). The LBEG beam might be considered a more typical LINAC beam and lends itself to continuous monitoring of transmission through the various stages of the accelerator. For the MPEG, however, the widely-spaced micropulses and lack of the requisite sensitivity current monitors throughout the accelerator make these measurements extremely difficult and more uncertain. Therefore, to assist in estimating the MPEG beam transmission, beam-dynamics simulations using HPSim were employed. These beam transmission estimates from the end of the Drift Tube Linac (DTL) to Target 4 (MPEG) and the PSR stripper foil (LBEG) are vital to determine the charge requirements for LAMP’s front end. In this technote, we demonstrate three major efforts in determining the transmission: (1) Convert the WNR beamline lattice from TRANSPORT to HPSim for use in the simulation; (2) Simulate the optimization process in the Central Control Room (CCR) that brings down the losses by up to a factor of 4000 between the end of the initial physics-based tuneup phase and production beam operation; (3) Analyze operational data to deduce measured transmissions. Table 1 shows the estimated losses with HPSim and operational analysis. Finally, a better measurement and other improvements to refine the results are also proposed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Markov state models from short non-equilibrium simulations—Analysis and correction of estimation bias

To begin, many state-of-the-art methods for the thermodynamic and kinetic characterization of large and complex biomolecular systems by simulation rely on ensemble approaches, where data from large numbers of relatively short trajectories are integrated. In this context, Markov state models (MSMs) are extremely popular because they can be used to compute stationary quantities and long-time kinetics from ensembles of short simulations, provided that these short simulations are in “local equilibrium” within the MSM states. However, over the last 15 years since the inception of MSMs, it has been controversially discussed and not yet been answered how deviations from local equilibrium can be detected, whether these deviations induce a practical bias in MSM estimation, and how to correct for them. In this paper, we address these issues: We systematically analyze the estimation of MSMs from short non-equilibrium simulations, and we provide an expression for the error between unbiased transition probabilities and the expected estimate from many short simulations. We show that the unbiased MSM estimate can be obtained even from relatively short non-equilibrium simulations in the limit of long lag times and good discretization. Further, we exploit observable operator model (OOM) theory to derive an unbiased estimator for the MSM transition matrix that corrects for the effect of starting out of equilibrium, even when short lag times are used. Finally, we show how the OOM framework can be used to estimate the exact eigenvalues or relaxation time scales of the system without estimating an MSM transition matrix, which allows us to practically assess the discretization quality of the MSM. Applications to model systems and molecular dynamics simulation data of alanine dipeptide are included for illustration. The improved MSM estimator is implemented in PyEMMA of version 2.3.

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

High Performance Computing to Enable Next-Generation Low-Temperature Waste Heat Recovery

US manufacturing industry fails to recover an estimated 900 trillion BTUs of low-temperature waste heat from its processes each year. A grand research challenge has been to develop waste heat recovery technology that can be applied to industrial manufacturing processes and vehicle operations. The Oak Ridge National Laboratory (ORNL)/Eaton research team performed computational research and development to design an innovative direct-contact heat exchanger (DCHE) technology to deliver low-cost, compact, longer-lifetime, high-efficiency waste heat recovery that is optimized for a low-temperature organic Rankine cycle. ORNL resources and expertise in high performance computing (HPC) and multiphase flows were utilized to realize this goal while advancing the fundamental understanding of two-phase, two-immiscible-fluid turbulent flow heat transfer for heat exchangers. In the first stage of the project, a computational fluid dynamics (CFD) model of a DCHE was developed. Verification and limited validation of the numerical solution were performed using experimental data from published literature. To support efforts in verification and validation, two-dimensional CFD models were also developed with the objective of investigating different boundary conditions to reach converged, steady-state solutions.Based on benchmarking efforts, a baseline industry-grade design of a DCHE was developed, which consists of two horizontal pipes of different cross-sectional areas – joined by a converging-diverging nozzle in the middle and relies on a three-phase system with co-current flows for operation. Evaporation of the cold liquid (n-pentane or methanol) by the hot liquid water occurs in the first pipe. The first pipe contains two inlets for liquid water and liquid pentane (or liquid methanol). The second pipe with the larger cross-sectional area serves as a gravity-driven phase-separator. It contains two outlets that are staggered to avoid entrainment of the liquid water phase by the gas phase. This baseline design resulted from several trial-and-error evaluations using two-dimensional CFD models. Confirmatory three-dimensional CFD models were also run for different mass flow rates and inlet temperatures (inlet water temperature ranging from 50°C to 90°C), and efficiency plots were produced to fully characterize the proposed industrial-grade design of the DCHE.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Anthromes and forest carbon responses to global change

Human effects on ecosystems date back thousands of years, and anthropogenic biomes—anthromes—broadly incorporate the effects of human population density and land use on ecosystems. Forests are integral to the global carbon cycle, containing large biomass carbon stocks, yet their responses to land use and climate change are uncertain but critical to informing climate change mitigation strategies, ecosystem management, and Earth system modeling. Using an anthromes perspective and the site locations from the Global Forest Carbon (ForC) Database, we compare intensively used, cultured, and wildland forest lands in tropical and extratropical regions. We summarize recent past (1900-present) patterns of land use intensification, and we use a feedback analysis of Earth system models from the Coupled Model Intercomparison Project Phase 6 to estimate the sensitivity of forest carbon stocks to CO 2 and temperature change for different anthromes among regions. Modeled global forest carbon stock responses are positive for CO 2 increase but neutral to negative for temperature increase. Across anthromes (intensively used, cultured, and wildland forest areas), modeled forest carbon stock responses of temperate and boreal forests are less variable than those of tropical forests. Tropical wildland forest areas appear especially sensitive to CO 2 and temperature change, with the negative temperature response highlighting the potential vulnerability of the globally significant carbon stock in tropical forests. The net effect of anthropogenic activities—including land-use intensification and environmental change and their interactions with natural forest dynamics—will shape future forest carbon stock changes. These interactive effects will likely be strongest in tropical wildlands.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Rate/Concentration Kinetic Petals: A Transient Method to Examine the Interplay of Surface Reaction Processes

Transient pulse response experiments are used to construct rate/concentration (RC) kinetic dependencies—RC petals—and provide a new method to distinguish the timing and interplay of adsorption, surface reaction, and product formation on complex (industrial) materials. A petal shape arises as the dynamic “reaction-diffusion” experiment forces the concentration and reaction rate to return to zero. In contrast to the typical steady-state “Langmuir-type” RC dependence, RC petals have two branches, which arise as a result of decoupled gas and surface concentrations in the non-steady-state regime. To demonstrate this approach, the characteristics of petal shapes using ammonia decomposition as a probe reaction are presented. Ammonia, hydrogen, and nitrogen transformation rates are compared on three simple materials: Fe, Co, and a bimetallic CoFe preparation when ammonia is pulsed at 550°C in a low-pressure diffusion reactor. All materials demonstrate a two-branch kinetic RC dependence for ammonia adsorption, and rate constants are quantified in the low-coverage regime. We found that H 2 and N 2 product formation was dependent on the concentration of surface intermediates for all materials with one exception: for Co, an additional fast hydrogen generation process was observed; the rate of which coincided with ammonia adsorption. Nitrogen generation was only significant for CoFe and Co and on the CoFe catalyst, a self-inhibition property was observed. A method for estimating the number of active sites based on the RC petals is presented and was applied to the Fe and CoFe samples. The surface coverage and rate of formation/conversion of surface intermediates are interpreted from the examination of shape characteristics of the RC petals for each material.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Physics-Informed Active Learning With Simultaneous Weak-Form Latent Space Dynamics Identification

The parametric greedy latent space dynamics identification (gLaSDI) framework has demonstrated promising potential for accurate and efficient modeling of high-dimensional nonlinear physical systems. However, it remains challenging to handle noisy data. Here, to enhance robustness against noise, we incorporate the weak-form estimation of nonlinear dynamics (WENDy) into gLaSDI. In the proposed weak-form gLaSDI (WgLaSDI) framework, an autoencoder and WENDy are trained simultaneously to discover intrinsic nonlinear latent-space dynamics of high-dimensional data. Compared with the standard sparse identification of nonlinear dynamics (SINDy) employed in gLaSDI, WENDy enables variance reduction and robust latent space discovery, therefore leading to more accurate and efficient reduced-order modeling. Furthermore, the greedy physics-informed active learning in WgLaSDI enables adaptive sampling of optimal training data on the fly for enhanced modeling accuracy. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including viscous and inviscid Burgers' equations, time-dependent radial advection, and the Vlasov equation for plasma physics. With data that contains 5%–10% Gaussian white noise, WgLaSDI outperforms gLaSDI by orders of magnitude, achieving 1%–7% relative errors. Compared with the high-fidelity models, WgLaSDI achieves 121 to 1779x speed-up.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Gaussian integral method for void fraction

Here, a novel method, the Gaussian Integral Method (GIM), is presented for calculating void fractions in Computational Fluid Dynamics–Discrete Element Method (CFD-DEM) simulations. GIM is versatile and applicable to various grid types, including structured and unstructured polyhedral meshes, without requiring special boundary treatments. An optimization technique is introduced to make GIM independent of grid resolution and type. The method is validated against experimental data from a fluidized bed, demonstrating that GIM produces realistic simulations closely resembling experimental observations. Additionally, unstructured polyhedral grids using GIM outperform structured grids of equivalent resolution, yielding results more aligned with experimental data. The gradient of the void fraction is computed in the CFD solver and utilized in the DEM solver for precise estimation at particle locations. Overall, GIM provides an effective solution for void fraction calculations in particulate media simulations with complex geometries, enhancing the accuracy and applicability of CFD-DEM simulations for industrial processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗