Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Predictive Control”

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 415 records · Page 23

Cluster Dynamics Modeling Needs for the Advanced Materials and Manufacturing Technologies Program

This milestone report aims to identify and assess the cluster dynamics (CD) modeling requirements within the Department of Energy's Office of Nuclear Energy (DOE-NE) Advanced Materials and Manufacturing Technologies (AMMT) program and to communicate these needs to the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The goal is to ensure NEAMS is well-informed about the CD modeling requirements to support AMMT's mission of accelerating the development, qualification, demonstration, and deployment of advanced structural materials and manufacturing for nuclear energy applications. CD modeling is an essential tool for predicting the degradation of structural materials under irradiation, which is a key component of AMMT's accelerated qualification process. The AMMT program focuses on both additively manufactured and wrought structural alloys, such as laser powder-bed fusion 316H austenitic stainless steel, alloy 709, Haynes 244, and alloy 617. These materials require a generalized CD modeling framework to facilitate rapid model development and computational simulation. A flexible, generalized CD software, similar to the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework, would enable modeling of various cluster types, including defect clusters, defect-solute clusters, and multicomponent clusters, incorporating thermodynamics and kinetics parameters. Radiation effects, microstructural feature evolution, and multi-dimensional modeling are critical considerations for the CD model. The usability of the CD code should allow for easy modification and coupling with MOOSE-based simulations. Additionally, the software should adhere to Nuclear Quality Assurance-1 standards, include a testing suite for verification and validation, and be version-controlled within a national laboratory-managed Git repository. Benchmark problems are needed to assess code predictions and performance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Satellite Reentry Predictions During Sudden Stratospheric Warmings

Prediction of uncontrolled satellite reentry time and trajectory are essential to avoid damages to human being and properties over land. Reentry time and trajectory are largely controlled by aerodynamic drag and mass neutral density below 200 km. In this paper, using a numerical model of the whole atmosphere, WACCM‐X, and a Precision Orbital Determination software, we demonstrate the major sudden stratospheric warming (SSW) in January 2009, during the deep solar minimum of 2009, may cause 0.2–0.3 day of variation in reentry time. The SSW can alter the mass density and composition mixing ratio in the lower thermosphere below 200 km through changed wave forcing and general circulation in the lower thermosphere. This is comparable to the effect of a 2‐SFU (Solar Flux Unit) solar radiance variation on the neutral density. 0.1 day of reentry time uncertainty corresponds to a global shift in impact location, that is, thousands of kilometers distance. For comparison, a major SSW during solar maximum contributes little to reentry time variation, which is predominantly controlled by the solar and geomagnetic forcing. We suggest a few proxies such as O/N 2 in the lower thermosphere or mean meridional wind to account for the lower atmosphere impacts in an operational model. This needs to be investigated in rigorous future work.

54 ENVIRONMENTAL SCIENCES↗

Topography, surface water distribution and subsurface structure in 2023 across an Arctic coastal tundra site near Utqiagvik, Alaska

Subsurface electrical resistivity tomography (ERT), active layer thickness measurements, photogrammetry, and topographic data were collected in September 2023 along a 475 m long, 20 m wide corridor that traverses various polygon types within the Barrow Environmental Observatory (BEO) on the Alaskan Arctic Coastal Plain, approximately 4 miles from the Beaufort Sea near Utqiaġvik, Alaska. These measurements were designed to assess decadal changes in surface water distribution, topography, and subsurface structure across this dynamic landscape. This archive contains the datasets acquired in 2023 and references to the datasets acquired previously at the same location. The ERT survey was conducted along the 475 m transect using 0.5 m electrode spacing and a roll-along acquisition strategy. Thaw layer thicknesses were measured with a tile probe along the same transect. Photogrammetry data were acquired using an unoccupied aerial vehicle (UAV) and were used to generate a digital elevation model and an RGB mosaic. A real-time kinematic (RTK) GPS was used to survey the ERT electrodes and the ground control points for the aerial imagery. The dataset contains 5 *.csv data files, 6 *.csv metadata files, and 6 *.tif files.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Advances in 3D transient plasma dynamics and control through MHD and hybrid fluid-kinetic simulations with JOREK

Transient phenomena and their control are of high relevance in magnetic confinement fusion plasmas to guarantee a stable and safe plasma operation. Interpretative simulations can maximize the insights gained from experiments on present machines and predictive simulations can help in the preparation of design, mitigation techniques and operational scenarios for future devices. In this article, we provide an overview of recent advances and novel scientific results obtained with the 3D non-linear hybrid fluid-kinetic code JOREK, covering physics of plasma transients from the core to the scrape-off layer (SOL) both for tokamak and stellarator devices. Substantial progress was made in the physics understanding, model validation with experiments and experiment interpretation, thus, giving confidence for predictions to devices like DTT, ITER and DEMO. The topics addressed comprise a wide range: the edge physics of new operation scenarios and edge localized mode suppression; major disruptions with a focus on runaway electrons and vertical displacement events as well as disruption mitigation by shattered pellet injection; the physics mechanisms and operational limits of the flux pumping regime for sawtooth control; MHD limits of stellarators and work towards incorporating advanced edge/SOL/exhaust dynamics; continuing improvements of the code for more efficient hybrid simulations on conventional and accelerated high performance computing architectures.

disruptions↗

Online thermal profile prediction for large format additive manufacturing: A hybrid CNN-LSTM based approach

Large format additive manufacturing (LFAM) is an advanced 3D printing technique that efficiently fabricates large-scale components through a layer-by-layer extrusion and deposition process. Accurate surface layer temperature monitoring is essential to prevent manufacturing failures and ensure final product quality. Traditional physics-based offline approaches for simulating thermal behavior are often inefficient and complex, posing challenges on real-time, in-situ monitoring. Here, to address this, we propose a data-driven hybrid CNN-LSTM model to predict sequential thermal images of arbitrary length using real-time infrared thermal imaging. In this approach, a Convolutional Neural Networks (CNN) is trained offline to capture spatial features, reduce dimensional complexity, and enhance time efficiency, while a stacked Long Short-Term Memory (LSTM) is applied online to capture temporal information for improved prediction of future thermal behavior in subsequent printing layers. Model performance is evaluated using MSE, SSIM, and PSNR metrics and is benchmarked against stacked LSTM and convolutional LSTM models, demonstrating superior accuracy and applicability. Additionally, to mitigate noise from moving extruders and gantry backgrounds in thermal images, a fine-tuned semantic segmentation model is implemented offline to extract printing geometry, enabling precise temperature tracking along the tool path for further thermal analysis. The frameworks developed in this study significantly advance temperature monitoring, thermal analysis, and in-situ manufacturing control for LFAM, bridging the gap between theoretical modeling and practical application.

Geometry extraction↗

Online LIBS–ML Framework for Dynamic Characterization of Heterogeneous Waste-Derived Gasification Feedstocks

LIBS−ML framework for real time feedstock characterization during continuous conveyor transport Heterogeneous waste derived feedstocks (e.g., waste coal, biomass and blends) introduce rapid variability in heating value and ash chemistry that affect gasifier operation, yet conventional laboratory characterization techniques are too slow to support proactive control. To address this gap, this study reports on an online, in situ, dynamic characterization framework that couple’s laser-induced breakdown spectroscopy (LIBS) with leakage safe machine learning (ML) regression to deliver real time, decision quality predictions of gasifier relevant properties. A controlled sample matrix spanning two different waste coals, two different biomasses, and engineered blends under two particle size conditions were constructed and benchmarked using standardized laboratory analyses for proximate/ultimate properties and ash composition. LIBS spectra were acquired dynamically as material flowed on a conveyor belt, using high energy 1064 nm laser ablation and shot averaging to improve repeatability and precision. Supervised regression models (multi layer perceptron (MLP) /artificial neural network (ANN), random forest (RF), and support vector regression (SVR)) and an optimized weighted ensemble were trained on emission line feature sets using nested cross validation with Bayesian hyperparameter tuning and validated against an independent hold out set. The proposed LIBS−ML workflow achieves near laboratory predictive fidelity across parametric targets (including higher heating value (HHV), ash content, fixed carbon, sulfur, major ash forming oxides, and initial deformation temperature (IDT)), with the weighted ensemble providing a robust default predictor under dynamic measurement conditions. These results demonstrate a practical pathway for real time feedstock characterization that can enable feedforward adjustments and more resilient gasifier operation for variable quality waste derived fuels.

Biomass↗

Towards a More Predictive Framework for Laser-Driven Particle Sources through Experimental Data-Informed Models

Laser-driven particle acceleration (LDPA) has emerged as a critical technology for high-energydensity physics applications since its discovery at Lawrence Livermore National Laboratory twenty years ago. However, realizing the full potential of these particle sources requires understanding the fundamental acceleration mechanisms and developing enhanced target designs for improved performance. This research addressed the need for controllable, high-performance laser-driven proton sources through two complementary approaches: experimentally investigating sheath field dynamics in multi-picosecond laser regimes and developing novel three-dimensional printed microstructured targets to achieve enhanced particle acceleration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Identifying rolling resistance and air resistance simultaneously for an electric truck

Accurately estimating rolling and air resistance is essential for predicting the energy consumption of vehicles. This study presents a field-based approach using a rolldown test to simultaneously determine rolling and air resistance coefficients. Unlike prior methods that frequently used simulations or models, we employ a goal programming methodology to improve precision and evaluate the actual vehicle and environmental conditions. Our methodology was tested using a Class 8 Freightliner eCascadia on a surveyed road section, ensuring controlled conditions for data collection. By analyzing the time–velocity relationship across multiple test runs, we derived resistance coefficients for both loaded and unloaded conditions. The study confirms that rolling resistance is largely independent of velocity at low speeds but exhibits a nonlinear dependency at higher speeds. Additionally, road surface conditions, tire condition, axle configuration, aerodynamic properties, and weather conditions significantly impact resistance values, emphasizing the need for real-world testing rather than relying solely on standardized projections. Our results align with existing literature while demonstrating the efficacy of the goal programming approach in refining resistance estimates. This work contributes to improved vehicle energy modeling, offering practical insights for fleet operators and policymakers seeking accurate energy consumption predictions for electric trucks operating under varying environmental conditions.

47 OTHER INSTRUMENTATION↗

The use of a benign fast-growing cyanobacterial species to control microcystin synthesis from Microcystis aeruginosa

Introduction Microcystis aeruginosa(M. aeruginosa), one of the most abundant blue-green algae in aquatic environments, produces microcystin by causing harmful algal blooms (HABs). This study investigated the combined effects of nutrients and competition among cyanobacterial subpopulations on the synthesis of microcystin-LR. Methods Under varying nitrogen and phosphorus concentrations, cyanobacterial coculture, and the presence of algicidal DCMU, the growth was monitored by optical density analysis or microscopic counting, and the microcystin production was analyzed using high-performance liquid chromatography-UV. Furthermore, growth and toxin production were predicted using a kinetic model. Results and discussion First, coculture with the fast-growing cyanobacteriumSynechococcus elongatusUTEX 2973 (S. elongatus) reducedM. aeruginosabiomass and microcystin production at 30°C. Under high nitrogen and low phosphorus conditions,S. elongatuswas most effective, limitingM. aeruginosagrowth and toxin synthesis by up to 94.7% and 92.4%, respectively. Second, this biological strategy became less effective at 23°C, whereS. elongatusgrew more slowly. Third, the photosynthesis inhibitor DCMU (3-(3,4-dichlorophenyl)-1,1-dimethylurea) inhibitedM. aeruginosagrowth (at 0.1 mg/L) and microcystin production (at 0.02 mg/L). DCMU was also effective in controlling microcystin production inS. elongatus–M. aeruginosacocultures. Based on the experimental results, a multi-substrate, multi-species kinetic model was built to describe coculture growth and population interactions. Conclusion Microcystin from representative toxin-producingM. aeruginosacan be controlled by coculturing fast-growing benign cyanobacteria, which can be made even more efficient if appropriate algicide is applied. This study improved the understanding of the biological control of microcystin production under complex environmental conditions.

Microbiology↗

Predicting multi-nodal in-nozzle particle interactions in high-viscosity fluid mediums for acoustophoretic direct-ink writing of line-patterned composites

Patterned functional materials offer improved properties (electrical, thermal, etc.) over their bulk counterparts in many applications, including energy storage, flexible electronics, and sensors. However, manufacturing approaches for patterning materials over large areas with features on the order of hundreds of microns or less are limited. Acoustophoresis, which uses acoustic forces to control particle arrangement in a fluid medium, is a pathway to address this challenge. This process is dependent on particle and fluid properties and enables patterning of a broad range of materials. Herein, a model with experimental validation is presented to demonstrate that acoustophoresis can be combined with direct-ink writing (DIW) to fabricate line patterns over large cm-scale areas. An in-nozzle particle interaction model was developed to investigate the impact of processing conditions on multi-nodal acoustophoretic DIW. The model predicts patterned line widths within a factor of two relative to experimental results for a high viscosity case study. Here, the model was used to investigate the impact of frequency, particle loading, particle radius, and acoustic pressure on line width and patterning time, providing critical feedback regarding the processing conditions suitable for a target application. Model results illustrate that frequency has the greatest impact on line patterns: increasing from 1 to 3 MHz resulted in a greater than 65% reduction in line width and a greater than 85% reduction in patterning time. Additionally, experiments were conducted with an alumina-epoxy ink and a ~21 cm 2 area pattern was rastered in ~5.5 minutes, demonstrating a path towards large-area line-patterned composite fabrication.

25 ENERGY STORAGE↗

Inverse Thermodynamics: Designing Interactions for Targeted Phase Behavior

The traditional goal of inverse self-assembly is to design interactions that drive particles toward a desired target structure. However, achieving successful self-assembly also requires tuning the thermodynamic conditions under which the structure is stable. In this work, we extend the inverse design paradigm to explicitly address this challenge by developing a framework for inverse thermodynamics, i.e., the design of interaction potentials that realize specific thermodynamic behavior. As a step in this direction, using patchy particle mixtures as a model system, we demonstrate how precise control over both bonding topology and bond energetics enables the programming of targeted phase behavior. In particular, we establish design principles for azeotropic demixing and show how to create mixtures that exhibit azeotropy at any prescribed composition. Our predictions are validated through Gibbs-ensemble simulations [Panagiotopoulos, Mol. Phys. 1987, 61, 813−826]. These results highlight the necessity of coupling structural design with thermodynamic engineering, and provide a blueprint for controlling complex phase behavior in multicomponent systems.

Azeotropes↗

Influence of Shear Strength Assumptions on BISON Debonding Simulations

Accurately predicting the thermomechanical response of buffer–IPyC debonding in TRISO fuel particles requires reliable mechanical property inputs for each coating layer, particularly the normal and shear strengths that influence interlayer delamination and stress concentrations. Micro tensile testing of AGR-2 fuel particles provided experimentally measured normal strengths for the buffer, IPyC, and buffer–IPyC interface; however, shear strength was not measured. As a result, BISON simulations of interface debonding must rely on assumed shear strength values, typically estimated as 20–40% of the measured ultimate tensile strength. This study evaluates how these assumed shear strength values influence cohesive zone model (CZM) predictions of buffer–IPyC separation in AGR 2 TRISO particles. Using micro tensile data from three AGR 2 compacts (2 1 3, 5 1 3, and 6 3 3), BISON simulations were performed with multiple shear strength assumptions to quantify their effect on radial and tangential stress evolution, debonding, and gap propagation. The results show that shear strength is a high sensitivity parameter: increasing the assumed shear strength significantly alters the stress distribution at the buffer–IPyC junction, shifts the predicted debonding location, and changes the extent of partial gap formation. While normal strength controls the initiation of interface separation, shear strength strongly influences the mode mixity of the failure process and the resulting stress concentrations transmitted to the IPyC and SiC layers. These findings highlight a critical gap in current TRISO mechanical characterization. Without experimentally measured shear strength, BISON simulations must rely on approximations that introduce uncertainty into predictions of coating layer integrity and fission product barrier performance. Future fuel qualification campaigns should therefore consider measurement of shear strength at the interlayer interfaces to reduce model uncertainty and improve the fidelity of TRISO fuel performance simulations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Physics-Guided Machine Learning (PGML) for Improved Aerostructure Manufacturing

The Physics-Guided Machine Learning (PGML) for Improved Aerostructure Manufacturing project objective is to automatically tune machining parameter predictions from physics-based models using process data and Bayesian machine learning. The intent is to enable a step change in aerospace manufacturing by combining machine learning, physics-based process models, and sensors/data in a comprehensive digital environment that simultaneously considers the computer numerically controlled (CNC) machining center capabilities, the workpiece material and geometry, and the workpiece support (fixturing). The project hypothesis is that this combination will enable improved performance in machining operations.

42 ENGINEERING↗

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗

Wind farm structural response and wake dynamics for an evolving stable boundary layer: computational and experimental comparisons

Abstract. The wind turbine design process requires performing thousands of simulations for a wide range of inflow and control conditions, which necessitates computationally efficient yet time-accurate models, especially when considering wind farm settings. To this end, FAST.Farm is a dynamic-wake-meandering-based mid-fidelity engineering tool developed by the National Renewable Energy Laboratory targeted at accurately and efficiently predicting wind turbine power production and structural loading in wind farm settings, including wake interactions between turbines. This work is an extension of a study that addressed constructing a diurnal cycle evolution based on experimental data (Quon, 2024). Here, this inflow is used to validate the turbine structural and wake-meandering response between experimental data, FAST.Farm simulation results, and high-fidelity large-eddy simulation results from the coupled Simulator fOr Wind Farm Applications (SOWFA)–OpenFAST tool. The validation occurs within the nocturnal stable boundary layer when corresponding meteorological and turbine data are available. To this end, we compared the load results from FAST.Farm and SOWFA–OpenFAST to multi-turbine measurements from a subset of a full-scale wind farm. Computational predictions of blade-root and tower-base bending loads are compared to 10 min statistics of strain gauge measurements during 3.5 h of the evolving stable boundary layer, generally with good agreement. This time period coincided with an active wake-steering campaign of an upstream turbine, resulting in time-varying yaw positions of all turbines. Wake meandering was also compared between the computational solutions, generally with excellent agreement. Simulations were based on a high-fidelity precursor constructed from inflow measurements and using state-of-the-art mesoscale-to-microscale coupling.

17 WIND ENERGY↗

Harnessing ionic complexity: A modeling approach for hierarchical ionic circuit design

Since the 1950s, soft ionic devices have evolved from individual components to an expanding library of sensors, actuators, signal transmitters, and processors. However, integrating these components into complex, multifunctional systems remains challenging due to the nonintuitive and nonlinear interactions between ionic elements. In this work, we address these fundamental challenges by developing a lumped element model that enables interrogation of the physics that governs ionic circuits, as well as rapid design and optimization. Our model captures features specific to ionic charge carriers, while preserving the hierarchical design flexibility and computational efficiency of traditional circuit modeling. We demonstrate that our model can not only fit individual device behavior but also accurately predict the behavior of larger circuits formed by combining those devices. Additionally, we show how our tool utilizes the intrinsic nonlinearities of ionic systems to enable extended functionality, revealing how factors such as ion enrichment, ion leakage, and polymer charge density influence performance. Lastly, we present a fully ionic power supply, sensor, control system, and actuator for a soft robot that adapts its motion in response to environmental salt, illustrating the tool’s potential to accelerate advancements in chemical sensing, biointerfacing, biomimetic systems, and adaptive materials.

42 ENGINEERING↗

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗