SEARCH · Engineering Papers

Results for “transportation simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

DOE OSTI2026

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)↗

DOE OSTI2026

An International Round-Robin Study on Thermoelectric Module Testing and Development of Standard Power Generation Modules

An international round-robin study on thermoelectric power generation modules was conducted with nine participating laboratories. Two types of commercially available bismuth telluride modules, 30 mm × 30 mm and 40 mm × 40 mm, were used. A test protocol was followed with five temperature set points from 50°C to 150°C. Graphite sheets were used as thermal interface materials with test pressure at 100 psi (0.69 MPa). The results showed large lab-to-lab variations and the key source of uncertainty for module efficiency was identified as the heat flux measurement. In the meantime, significant uncertainty was also found in maximum electrical power (P max ) measurements. As a result of the round-robin, a “standard module” with 4 × 4 legs on a 20 mm × 20 mm platform was suggested. A skutterudite module and a half-Heusler module were produced with identical geometry and 4 mm × 4 mm × 8 mm legs. All transport properties to calculate the figure-of-merit, zT, were measured from ambient temperature to 500°C. Module performance was measured by two laboratories. Two finite-element-analysis (FEA)-based models were developed independently to simulate and predict the module performance. In conclusion, the standard modules eliminated significant test uncertainties and are aimed at assisting device design and achieving more accurate performance predictions.

Round-robin↗

DOE OSTI2026

MAPS: the MFEM Anisotropic Plasma Solver

Simulating magnetically confined fusion plasmas presents a uniquely challenging problem due to the nonlinear anisotropic heat conduction. We introduce the MAPS (MFEM Anisotropic Plasma Solver) tool, which uses a high-order finite element method to compute transport solutions on unstructured meshes. We show results for a set of three 2-D verification tests, two of which demonstrate the expected convergence properties for various mesh resolutions and polynomial degrees. We then discuss the convergence rate for the third test.

Barnett, Rhea [ORNL] (ORCID:0000000317527979)↗

DOE OSTI2024

Magnetically Arrested Circumbinary Accretion Flows

Abstract Binary systems with comparable masses and a surrounding accretion disk can accrete gas through spiral accretion streams penetrating the central cavity formed by tidal interactions. Using three-dimensional Newtonian magnetohydrodynamics simulations, we investigate the possibility of a magnetically arrested accretion flow through the cavity. Rather than solely continuously feeding the binary through spiral accretion streams, the accretion is regulated by the strong magnetic field inside the cavity. Transport of mass and angular momentum onto the binary then proceeds largely periodically in magnetic flux eruption episodes. The ejected flux tubes carry angular momentum outward and away from the binary, inject hot plasma into the disk, and can launch flares. This likely intermittent scenario could have potential implications for the emission signatures of supermassive black hole binaries and shed light onto the role magnetic fields play in the binary’s orbital evolution.

Astronomy & Astrophysics↗

DOE OSTI2026

Stochastic Thermo-Hydro Modeling and Neural Network Surrogate Development for Thermal Resource Assessment of the Galleries-to-Calories Geobattery

The Galleries-to-Calories Geobattery concept explores the use of abandoned coal mine workings for large-scale thermal energy transport and storage. The system involves injecting waste heat from a supercomputing facility into flooded mine galleries, where groundwater flow can store and transport thermal energy for potential recovery in downgradient district heating and cooling applications. To evaluate the feasibility and performance of the Geobattery under geological and operational uncertainty, we developed a suite of stochastic thermo-hydrological (TH) simulations using Monte Carlo sampling of key uncertain parameters (e.g., permeability, porosity, thermal conductivity, specific heat capacity) and operating conditions (e.g., injection rate, injection temperature). Results identified injection rate and temperature as the most influential parameters governing thermal front propagation, while the geometry of the room-and-pillar structure played a critical role in directing the extent and orientation of thermal advancement. Optimal combinations of material properties for maximizing heat recovery were also determined. To address the high computational cost of coupled-process stochastic modeling, we trained a neural network surrogate model on 24,000 physics-based realizations, achieving an R² > 0.99 and MAE < 0.1 for temperature predictions at monitoring locations. This surrogate enabled an additional 100,000 realizations for global sensitivity analysis and probabilistic thermal resource assessment. The integrated stochastic physics–surrogate modeling framework offers a computationally efficient tool for quantifying uncertainty, identifying key drivers, and informing early-stage design decisions for Geobattery systems.

15 - GEOTHERMAL ENERGY↗

DOE OSTI2025

Community-Informed Urban Flood Modeling for Impact Mitigation

The intensification of the hydrologic cycle due to climate change poses a threat to aging and under-designed water infrastructure systems which cannot adequately manage intense storm events. Developing a comprehensive plan for managing rain-driven flooding events is challenging due to uncertainties in the magnitude and frequency of future storm events and conflicting stakeholder objectives. In the City of Baltimore, Maryland, stormwater infrastructure is struggling to keep up with rainfall-driven (pluvial) flooding events, which regularly damage housing and disrupt transportation for residents. In this study, a hybrid of community engagement, numerical modeling, and artificial intelligence techniques are employed to explore prospective urban flooding adaptations. Community engagement drives the development of an urban flooding model (EPA Storm Water Management Model) for the Baltimore Harbor watershed. The model integrates complex surface and subsurface stormwater infrastructure data from the City, high-resolution spatial data, insights from local public works experts, and the lived experiences of City residents. This co-developed model simulates adaptations of interest to stakeholders in the city, including green and grey infrastructure and operational management strategies. Stormwater management scenarios focused on inlet cleaning and spatially concentrated green infrastructure are found to be the most effective in reducing flood depths in community priority locations. Together, these adaptations can reduce the duration of intersection inundation by more than twenty minutes, allowing for quicker emergency response and restoration of typical transportation systems. Future work will combine this community engaged flooding model with the Deep Uncertainties Pathways framework to explore tradeoffs between adaptations and develop dynamic adaptations which align with community objectives, enhance climate resilience in Baltimore, and can be adjusted in response to changing future conditions.

Ava, Spangler [Pennsylvania State University]↗

DOE OSTI2025

Enabling Low-Temperature (LTP) Ignition Technologies for Multi-Mode Engines through the Development of a Validated High-Fidelity LTP Model for Predicative Simulations Tools

The goal of multi-mode engine architectures is to extend current lean-burn dilution limits with renewable fuels, which requires spark plugs to deposit high energies (hundreds of mJ) in order to initiate ignition and complete combustion. At elevated energy deposition rates, spark plugs experience increased electrode erosion and thermal losses, which ultimately shortens the spark-plug lifetime and lowers ignition efficiency. As such, in order to safeguard the efficiency gains of multi-mode concepts, new and improved ignition technologies are required. Recently, non-equilibrium low-temperature plasmas (LTP) have been shown to promote energy-efficient ignition via quenching and transport of electronically excited atoms and molecules, selective radical production and fast heating of hydrocarbon/air mixtures [1-2]. Thus, LTP is seen as a technology that can potentially improve the energy extraction efficiency of fuels, while enabling kinetically controlled combustion modes towards fuel leaner conditions to realize current DOE VTO goals of improving the sustainability of future mobility [3]. Although many previous studies have demonstrated the efficacy of plasma-assisted ignition to enhance combustion, the detailed enhancement mechanisms remain largely unknown, especially for oxygenated fuels and at elevated pressures that are most relevant to practical engine conditions. These barriers hinder the development of accurate and comprehensive numerical models that seek to describe LTP-based ignition in existing engine design software tools and methods. Current state-of-the-art simulation capabilities for LTP ignition systems are in need of improvements since they deliver qualitative results only due to important limitations of existing approaches. Firstly, validated kinetic models with elementary steps for plasma discharges in oxygenated fuel/air mixtures of relevance to the transportation sector are required. Such kinetic models do not exist at present and will be developed and validated within this project. Secondly, plasma discharges and reactive mixture ignition are multi-scale, unsteady processes requiring high-performance numerical methods and software that execute efficiently on DOE supercomputers. Such software does not exist at present and will be developed and applied to practical LTP ignition scenarios as part of this project. Thirdly, experimental databases that are tailored to serve as benchmark in support of the development of predictive computational models of LTP ignition do not exist and will be part of this project.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Magnetically tuned metal-insulator transition in LaAlO 3 /SrTiO 3 nanowire arrays

A wide family of two-dimensional (2D) systems, including stripe-phase superconductors, sliding Luttinger liquids, and anisotropic 2D materials, can be modeled by an array of coupled one-dimensional (1D) electron channels or nanowire arrays. Here we report experiments in arrays of conducting nanowires with gate and field tunable interwire coupling, that are programmed at the LaAlO 3 /SrTiO 3 interface. We find a magnetically tuned metal-to-insulator transition in which the transverse resistance of the nanowire array increases by up to four orders of magnitude, which can be further tuned with a gate voltage. To explain this behavior, we develop a minimal model of a coupled two-wire system where a Wenzel-Kramers-Brillouin-based approach is used to estimate the transverse tunneling conductance. We demonstrate the existence of distinct conductance features and highlight the crucial role played by the field dependence of the interwire potential barrier on transport properties. Since our model makes minimal assumptions, we expect our predictions to hold for a wide class of coupled 1D systems. The nanowire arrays can serve as model systems to understand the origin of exotic behavior in correlated materials via analog quantum simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗

DOE OSTI2025

Long-lived neutron-star remnants from asymmetric binary neutron star mergers: element formation, kilonova signals and gravitational waves

We present 3D general-relativistic neutrino-radiation hydrodynamics simulations of two asymmetric binary neutron star mergers producing long-lived neutron stars remnants and spanning a fraction of their cooling time scale. The mergers are characterized by significant tidal disruption with neutron rich material forming a massive disc around the remnant. The latter develops one-armed dynamics that is imprinted in the emitted kilo-Hertz gravitational waves. Angular momentum transport to the disc is initially driven by spiral-density waves and enhanced by turbulent viscosity and neutrino heating on longer timescales. The mass outflows are composed by neutron-rich dynamical ejecta of mass ∼10 −3 –10 −2 M ⊙ followed by a persistent spiral-wave/neutrino-driven wind of ≳ 10 −2 M ⊙ with material spanning a wide range of electron fractions, ∼0.1–0.55. Dynamical ejecta (winds) have fast velocity tails up to ∼0.8 (∼0.4) c. The outflows are further evolved to days timescale using 2D ray-by-ray radiation-hydrodynamics simulations that include an online nuclear network. We find complete r-process yields and identify the production of 56 Ni and the subsequent decay chain to 56 Co and 56 Fe. Synthetic kilonova light curves predict an extended (near-) infrared peak a few days postmerger originating from r-process in the neutron-rich/high-opacity ejecta and UV/optical peaks at a few hours (ten minutes) postmerger originating from weak r-process (free-neutron decay) in the faster ejecta components. Additionally, the fast tail of tidal origin generates kilonova afterglows potentially detectable in radio and X band on a few to ten years time scale. Quantitative effects originating from the tidal disruption merger dynamics are reflected in the multimessenger emissions.

abundances↗

DOE OSTI2026

Aerosol-deep convection interaction based on joint cell-thermal tracking in Large Eddy Simulations during the TRACER campaign

In cumulus clouds, aerosol concentrations control cloud droplet concentrations, modifying cloud radiative properties, precipitation processes, and cloud electrification. However, mechanisms of aerosol-deep convection interactions are not well understood due to complex cloud dynamics and microphysics. We investigate the interaction of aerosols with isolated deep convection using Large Eddy Simulations of two cases during the TRacking Aerosol Convection interactions ExpeRiment (TRACER) near Houston, Texas, using a joint cell-thermal tracking algorithm. Cumulus thermals are droplet generators, since supersaturation and droplet nucleation coincide with thermal centers, where the strongest updrafts occur. Primary ice crystal formation does not take place inside thermals, but at layers where previous thermals detrained moisture. As subsequent thermals containing supercooled droplets penetrate these layers, hail and graupel form at or near these thermals. Higher aerosol concentrations result in higher droplet concentrations that suppress drizzle, delay warm rain processes, and transport more moisture aloft. This increases snow and ice amount, as well as graupel and hail, leading to more lightning. Polluted thermals initiate at slightly higher altitudes, and are slightly larger and faster, suggesting a weak invigoration. We also find more thermals per cell, but fewer isolated cells, since convection is more aggregated and intense, especially near the end of the 24 h simulation. Non-linear mesoscale feedback likely triggered by temperature and moisture responses to aerosol-thermal interactions causes the aggregation. Time-lagged aerosol-reinitialization experiments show that the mesoscale response is the predominant forcing for the invigoration. These changes happen within one day, on a smaller scale than previously suggested.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Numerical investigation of two-dimensional buoyancy-driven eddies in liquid metal magnetohydrodynamic flows in breeding blankets *

Abstract Lead–lithium flows are key features in the design of tokamak breeding blanket concepts such as the dual-coolant lead–lithium (DCLL). Since they flow under magnetic fields, they are affected by magnetohydrodynamic (MHD) effects. The neutron flux originating in the tokamak plasma heats the breeding blanket channels in a non-uniform manner, inducing buoyancy forces in the liquid metal. Buoyancy may become a source of quasi-two-dimensional (Q2D) turbulence, and the appearance of eddies may affect the transport of heat and tritium across the blanket. Blankets characterized by high-speed liquid metal flows (such as DCLL) will need ceramic insulating walls to reduce the MHD-related pressure drop in the channels. In our simulations, we have used the Q2D model proposed by Sommeria and Moreau (SM82) which is especially suitable for modeling electrically insulating channel flows. On top of that, we have modelled buoyancy forces in the momentum equation using the Oberbeck–Boussinesq approximation. In this work, we include a validation of the implemented Q2D model in buoyancy-driven cases and identify a thermohydraulic configuration that promotes the generation of eddies and the accumulation of tritium. We also show the results provided by our post-processing tool based on the bi-dimensional fast Fourier transform for eddy detection and characterization. We complete our investigation by performing an initial assessment of how eddies can accumulate tritium in the breeding blanket. We conclude by discussing the relevance of preventing accumulation of tritium within the blanket and suggest a possible solution.

Suarez, Daniel (ORCID:0000000251740100)↗

DOE OSTI2024

A theoretical/computational framework to measure SiO2 and MgO viscosity at high pressure

The convection of the mantle of Earth and super-Earths is important for many terrestrial phenomena, from plate tectonics to outgassing. Rheological properties, such as viscosity, regulate the transport of thermal energy and mass. However, the viscosity of mantle-relevant materials, such as MgO, at relevant pressures (>120 GPa) are not well constrained. The objective of this work was to develop a computational platform to simulate novel experiments aiming to measure the viscosity of MgO at high pressures. Experiments performed by our collaborators at Johns Hopkins University and Lawrence Livermore National Laboratory use the OMEGA EP laser facility to shock a corrugated MgO interface to 170 GPa, with resulting velocity evolution governed by the viscous Richtmyer-Meshkov instability. We used an in-house hydrocode to simulate this process and thus provide a bound on viscosity by comparing our simulations results to these experiments, as well as to examine the physical processes at play. We simulated the experiments with different values of MgO viscosity (from inviscid to 10,000 Pa∙s) while taking into account the unsteady laser pulse, the material's equation of state, and the rate-dependent constitutive relation for MgO, and the materials' equations of state. Our results suggest that MgO at these conditions has a viscosity within an order of magnitude of 5000 Pa∙s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

First-principles elucidation of defect-mediated Li transport in hexagonal boron nitride

Hexagonal boron nitride (hBN) is a promising candidate as a protective membrane or separator in Li-ion and Li–S batteries, given its excellent chemical stability, mechanical robustness, and high thermal conductivity. In addition, hBN can be functionalized by introducing defects and dopants, or be directly integrated into other active components of batteries, which further augments its appeal to the field. Here, we use first-principles simulations to evaluate the role of atomic defects in hBN in regulating the Li-ion diffusion mechanism and associated kinetics. Specifically, the following four distinct types of vacancy defects are considered: isolated single B and N vacancies, a B–N vacancy pair, and a B 3 N vacancy cluster. It is found that these defect sites generally favor Li intercalation and out-of-plane diffusion but slow down in-plane Li-ion diffusion due to a strong Li trapping effect at the defect sites. Such a trapping effect is, however, highly local such that it does not necessarily affect the overall Li-ion conductivity in defected hBN layers. The present systematic evaluation of the impact of atomic defects on Li ion migration and accompanied charge analysis of hBN lattice in response to Li-ion diffusion provide a mechanistic understanding of Li-ion transport behavior in defected hBN and highlight the potential of defect engineering to achieve optimal material performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Interactions between Per- and Polyfluoroalkyl Substances (PFAS) at the Water–Air Interface

Per- and polyfluoroalkyl substances (PFAS)--so-called “forever chemicals”--contaminate the drinking water of about 100 million people in the U.S. alone and are inefficiently removed by standard treatment techniques. A key property of these compounds that underlies their fate and transport and the efficacy of several promising remediation approaches is that they accumulate at the water−air interface. This phenomenon remains incompletely understood, particularly under conditions relevant to natural and treatment systems where water−air interfaces often carry significant loads of other organic contaminants or natural organic matter. To understand the impact of organic loading on PFAS adsorption, we carried out molecular dynamics simulations of PFAS at varying interfacial densities. We find that adsorbed PFAS form strong mutual interactions (attraction between perfluoroalkyl chains and electrostatic interactions among charged head groups) that give rise to ordered interfacial coatings. Furthermore, these interactions often involve near-cancellation of hydrophobic attraction and Coulomb repulsion. Our findings explain an apparent paradox whereby PFAS adsorption isotherms often suggest minimal mutual interactions while simultaneously displaying a high sensitivity to the composition and density of interfacial coatings. Consideration of the compounds present with PFAS at the interface has the potential to allow for more accurate predictions of fate and transport and the design of more efficient remediation approaches.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Edge ML for CAN bus intrusion detection in AVs

Autonomous Vehicles (AVs) are revolutionizing transportation, but their reliance on interconnected cyber-physical systems exposes them to unprecedented cybersecurity risks. This study addresses the critical challenge of detecting real-time cyber intrusions in self-driving vehicles by leveraging a dataset from the Udacity self-driving car project. We simulate four high-impact attack vectors, Denial of Service (DoS), spoofing, replay, and fuzzy attacks, by injecting noise into spatial features (e.g., bounding box coordinates) to replicate adversarial scenarios. We develop and evaluate two lightweight neural network architectures (NN-1 and NN-2) alongside a logistic regression baseline (LG-1) for intrusion detection. The models achieve exceptional performance, with NN-2 attaining an AUC score of 93.15% and 93.15% accuracy, demonstrating their suitability for edge deployment in AV environments. Through explainable AI techniques, we uncover unique forensic fingerprints of each attack type, such as spatial corruption in fuzzy attacks and temporal anomalies in replay attacks, offering actionable insights for feature engineering and proactive defense. Visual analytics, including confusion matrices, ROC curves, and feature importance plots, validate the models' robustness and interpretability. This research sets a new benchmark for AV cybersecurity, delivering a scalable, field-ready toolkit for Original Equipment Manufacturers (OEMs) and policymakers. By aligning intrusion fingerprints with SAE J3061 automotive security standards, we provide a pathway for integrating machine learning into safety-critical AV systems. Our findings underscore the urgent need for security-by-design AI, ensuring that AVs not only drive autonomously but also defend autonomously. This work bridges the gap between theoretical cybersecurity and life-preserving engineering, offering a leap toward safer, more secure autonomous transportation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Evaluating the Transport Mechanism for Lignocellulosic Inhibitors in Zymomonas mobilis

The conversion of lignocellulosic biomass to biofuels and bioproducts is limited by small inhibitory molecules generated during the pretreatment process. The diverse compounds within the cell wall hydrolysate including aldehydes, carboxylic acids, phenolics, and alcohols inhibit microbial growth by disrupting microbial membranes. Although Zymomonas mobilis combines high ethanol productivity with hopanoid-rich membranes that confer solvent tolerance, hydrolysate compounds still inhibit microbial growth. One proposed mitigation strategy is to pump these molecules out through active transport processes, reducing direct toxicity. However, if passive permeation rates are high, active transport would create a futile cycle where the exported molecules diffuse back in, creating a net drag on fitness. To investigate this hypothesis, we measure the passive permeation for thirty-three lignocellulose derived compounds using atomistic molecular dynamics simulations with enhanced sampling. Permeability coefficients determined spanned more than eight orders of magnitude driven largely by the balance between hydrophobicity and polarity. Hydrophobic inhibitors such as aromatic acids, aldehydes, and short chain alcohols readily partitioned into the membrane and showed high permeability consistent with rapid passive influx and strong toxicity. Most polar molecules have lower but still appreciable permeation rates, with only the most polar sugars having low permeabilities. By integrating these predictions together with hydrolysate composition information, we find that most of these molecules are so permeable that futile cycles where passive permeation overwhelms efflux is highly likely. These findings provide a predictive framework for strain engineering, highlighting that, for highly permeable inhibitors, reducing extracellular inhibitor concentrations through hydrolysate detoxification, pretreatment optimization, or dilution may be necessary to mitigate futile cycling and improve microbial productivity in industrial biorefineries.

Zymomonas mobilis↗

DOE OSTI2025

Co‐Evolution of Specific Stiffness and Permeability of Rock Fractures Offset in Shear

Abstract Fractures and faults represent planes of weakness and compliance in rock masses that serve as focal points for both microearthquakes and fluid transport, with seismicity and permeability evolution closely linked. Contact stiffness is highly stress‐sensitive and directly influences permeability. We explore the co‐evolution of specific stiffness and permeability of rough fractures under normal stress and shear offset using numerical simulations. Individual rough fractures are represented by variable amplitude (Root mean square) and wavelength ( λ ) using a granular mechanics model. Contacting rough surfaces are mated, offset in shear, and then compacted in displacement mode. The compacting fractures generate stress‐dependent changes in contact porosity, which govern both permeability and stiffness evolution. We establish a universal dimensionless relationship linking specific stiffness and permeability that inherently incorporates the effects of surface roughness, shear offset, and microcracking. The observed cracking effect—where local stress redistribution and pressure‐driven microcrack propagation dynamically alter the aperture field—introduces a nonlinear permeability response at high stress. Increased roughness amplitude and larger shear offsets reduce stiffness while dampening permeability sensitivity to stress, demonstrating a strong interplay between surface texture and hydro‐mechanical behavior. While the model captures this behavior effectively, deviations emerge at very low porosities due to extreme aperture sensitivity in this limit.

He, Xinxin [Energy and Mineral Engineering Pennsyl↗