Stellar Interpretation of Meteoritic Data and PLotting for Everyone (SIMPLE): Isotope Mixing Lines for Six Sets of Core-collapse Supernova Models
Explore the source record for details and available documents.
Physics↗
SEARCH · Search NASA
Results for “Thermodynamics of mixing”
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.
Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling
Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .
Close-clearance impellers↗
The impact of aerosol mixing state on immersion freezing: insights from classical nucleation theory and particle-resolved simulations
Immersion freezing, initiated by ice-nucleating particles (INPs) in supercooled aqueous droplets, plays an important role in the formation of ice crystals within clouds. The efficiency of immersion freezing depends strongly on INP composition and, crucially, on the mixing state – how chemical species are distributed across the particle population. Here, we quantify the impact of aerosol mixing state on immersion freezing using a combined theoretical and particle-resolved modeling approach. We derive analytical expressions for the frozen fraction of internally and externally mixed INP populations based on classical nucleation theory, showing that the frozen fraction is sensitive to whether ice-active species are present in all particles or only in a subset of the population. We introduce a multi-species immersion freezing scheme into the particle-resolved model PartMC, using the water activity-based immersion freezing model (ABIFM) to compute freezing probabilities for mixed-composition particles. To improve computational efficiency, we implement a Binned Tau-Leaping algorithm and demonstrate an order-of-magnitude speedup with minimal accuracy loss. Simulations reproduce the analytical trends in limiting cases and extend the analysis to more general aerosol populations, where mixing state continues to exert a substantial control on frozen fraction. Sensitivity analyses across particle size, species type, and cooling condition reveal that the mixing state effect is most pronounced when small amounts of highly efficient INPs are mixed with less efficient materials. These findings underscore the need to represent aerosol mixing state explicitly in models of heterogeneous ice nucleation to reduce uncertainty in cloud-phase partitioning.
54 ENVIRONMENTAL SCIENCES↗
Interface-induced fast Li+ transport in mixed ionic–electronic conductors
Interfacial instability between lithium metal and solid-state electrolytes limits the performance of all-solid-state lithium metal batteries (ASSLBs), leading to parasitic reactions, non-uniform Li+ flux, and dendrite growth. Here, we develop a composite interlayer composed of the anti-perovskite Li2OHCl0.75Br0.25 (AP) and carbon nanotubes (CNTs) to enhance both interfacial stability and ionic transport. The AP–CNT interlayer exhibits enhanced Li+ conductivity arising from interfacial electron transfer from AP to CNTs, which generates a built-in electric field that facilitates Li+ migration. Lithium symmetric cells incorporating this interlayer achieve a high critical current density of 2.4 mA cm−2 at 55 °C. This design integrates chemical robustness with coupled ion–electron transport, offering a generalizable strategy for safe, dendrite-free, and high-performance ASSLBs.
He, Chenche↗
The Reaction of Atomic Hydrogen with an Imine: Direct Kinetic Measurements and Computational Studies
Imines are of interest in the combustion modeling of nitrogen-based fuels such as ammonia, as intermediates in atmospheric oxidation of amines, and may offer a new class of environmentally benign substitutes for perhalogenated reagents. Elementary gas-phase rate constants for the consumption of H atoms by 1,1,1,3,3,3-hexafluoro-2-propanimine were measured over 294–736 K using the laser flash photolysis/atomic resonance fluorescence technique. The results are summarized as 3.0 × 10 –11 exp(−14.5 kJ mol –1 /RT) cm 3 molecule –1 s –1 with a confidence interval of ±13%. These are the first temperature-dependent kinetic experiments on an imine of the form R 2 C═NH and thereby provide the only validation of transition-state modeling for such systems. Several pathways are feasible; comparison of theory and experiment suggests the dominant channel is addition of H atoms to the nitrogen atom. On this basis, the observations are rationalized quantitatively. Computations also yield the thermochemistry of (CF 3 ) 2 CNH, with Δ f H 298 = −1232 kJ mol –1 .
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Computation of three-dimensional mixed convective boundary layer flow
The paper presents the numerical solution of heat and mass transfer during cross-flow (orthogonal) mixed convection. In this class of flow, a buoyancy-driven transport in the vertical direction and a forced convective flow in the horizontal direction results in a three-dimensional boundary layer structure adjacent to the plate. The rates of heat and mass transfer are determined by a combined influence of the two transport processes. The equations for the conservation of mass, momentum, energy, and species concentration were solved along with appropriate boundary conditions to determine the distributions of velocity components, temperature, and concentration across the thickness of the boundary layer at different locations on the plate. Results were expressed in dimensionless form using Reynolds number, Richardson number for heat transfer, Richardson number for mass transfer, Prandtl number, and Schmidt number as parameters. It was found that the transport is dominated by buoyancy at smaller vertical locations and at larger distances away from the forced convection leading edge. Effects of forced convection appeared to be very strong at smaller horizontal distances from the leading edge. The cross stream forced convection enhanced the rate of heat and mass transfer by a very significant amount.
Prashandt Gadepalli↗
Methanol mixing−controlled compression ignition enabled via homogeneous charge compression ignition of dimethyl ether through catalytic decomposition of methanol
Methanol is a potentially attractive fuel for marine and off−road engines owing to its availability at bunkering and global distribution locations. Although methanol is well−distributed worldwide, its fuel chemistry and ignition properties make it poorly suited as a direct drop−in replacement for diesel fuel in compression−ignition engines. However, industrial processes are regularly used to convert methanol, via catalytic dehydration, to dimethyl ether (DME) over nonprecious metal catalysts. This chemical conversion can occur at relatively low pressures, temperatures, and catalyst space velocities, highlighting a potential opportunity to generate DME via onboard catalytic dehydration of methanol. DME’s fuel kinetic and ignition properties for compression ignition are much more favorable than those of methanol or even diesel fuel, but DME is more challenging than diesel fuel or methanol to pump, store, and deliver through conventional diesel fueling injection hardware. Thus, a potential opportunity exists to use the ignition and kinetic properties of DME, with the transportation and delivery advantages of methanol, in a methanol−fueled mixing−controlled compression−ignition engine. The present work explores performance, combustion behavior, and emissions reduction opportunities for methanol mixing−controlled combustion, enabled by a HCCI of DME that represents a small fraction of the total fuel energy that can be generated onboard via catalytic dehydration of methanol.
Splitter, Derek [ORNL] (ORCID:0000000174044047)↗
Cosmological probes of Dark Radiation from Neutrino Mixing
Abstract Models of stepped dark radiation have recently been found to have an important impact on the anisotropies of the cosmic microwave background, aiding in easing the Hubble tension. In this work, we study models with a sector of dark radiation with a step in its abundance, which thermalizes after big bang nucleosynthesis by mixing with the standard model neutrinos. For this, we extend an earlier work which has focused on the background evolution only until the dark sector thermalizes by deriving the full background and perturbation equations of the model and implementing them in an Einstein-Boltzmann solving code. We expound on the behavior of this model, discussing the wide range of parameters that result in interesting and viable cosmologies that dynamically generate dark radiation during a range of epochs. We find that for the strongly self-coupled regime, there is no large cosmological impact for a tight prior on the mass, whereas larger mass ranges allow a smooth interpolation between a behavior close to the ΛCDM cosmological standard model and close to an additional component of strongly self-interacting dark radiation. In the weakly self-coupled regime we find that we can accommodate a parameter space relevant for the neutrino anomalies as well as one relevant to easing the Hubble tension.
Astronomy & Astrophysics↗
Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection
This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data
Navier Stokes Equations↗
Characterization of Flashback and Flame-Holding in a Jet-in-Crossflow Mixing Configuration With Methane–Hydrogen Fuel Blends
An approach that combines experimental and numerical analyses has been implemented to characterize the fundamentals of flashback events and flame-holding phenomena during high-hydrogen combustion in a jet-in-crossflow (JICF) configuration. Such flame dynamics are visualized experimentally using nanosecond (ns)-based hydroxyl planar laser-induced fluorescence (OH-PLIF) and chemiluminescence diagnostics techniques. The JICF burner has an optically accessible premixing tube allowing the optical diagnostics. The testing was conducted for varied premixer velocities (V) and equivalence ratio (ϕ) for 90–100% (H 2 , by mole) H 2 /CH 4 reactant mixtures at atmospheric temperature and pressure conditions. Two distinct flashback events were identified – conventional rich flashback and lean flashback, recorded while increasing ϕ and decreasing ϕ, respectively. The cause of lean flashback was attributed to lower momentum flux ratio which bends the jet sharply, closer toward the injection plane. The mean OH-PLIF images characterized the flame-holding behavior, where the flame was found to be stabilized on the leeward side only or on both windward and leeward sides as a lifted flame near the fuel port. A large eddy simulation (LES) with detailed chemistry combustion modeling approach was implemented along with an OH* submechanism, and it showed qualitative agreement with the integrated line-of-sight chemiluminescence results as well as the planar OH-PLIF measurements.
chemiluminescence diagnostic↗
Control of an Axisymmetric Jet Using Vortex Generators
The results of an experimental investigation on the effect of vortex generators, in the form of small tabs at the nozzle exit on the evolution of a jet, are reported in this paper. Primarily tabs of triangular shape are considered, and the effect is studied up to a jet Mach number of 1.8. Each tab is found to produce a dominant pair of counter-rotating streamwise vortices having a sense of rotation opposite to that expected from the wrapping of the boundary layer. This results in an inward indentation of the mixing layer into the core of the jet. A triangular-shaped tab with its apex leaning downstream, referred to as a delta tab, is found to be the most effective in producing such vortices, with a consequential large influence on the overall jet evolution. Two delta tabs, spaced 180 deg apart, completely bifurcate the jet. Four delta tabs stretch the mixing layer into four 'fingers,' resulting in a significant increase in the jet mixing downstream. For six delta tabs the mixing layer distortion settles back to a three finger configuration through an interaction of the streamwise vortices. The tabs are found to be equally effective in jets with turbulent or laminar initial boundary layers. Two sources of streamwise vorticity are postulated for the flow under consideration. One is the upstream 'pressure hill,' generated by the tab, which constitutes the main contributor of vorticity to the dominant pair. Another is due to vortex filaments shed from the sides of the tab and reoriented downstream by the mean shear of the mixing layer. Depending on the orientation of the tab, the latter source can produce a vortex pair having a sense of rotation opposite to that of the dominant pair. In the case of the delta tab, vorticity from the two sources add, explaining the strong effect in that configuration.
Streamwise Vorticity↗
Investigation of Lunar-Inspired Geopolymer Concrete Formulations Mixed and Cured in Microgravity on the International Space Station (ISS)
The research outlined in this presentation investigates the use of various lunar regolith simulants in geopolymer lunar concrete mixes mixed and cured on the International Space Station (ISS). The motivation for this work is to study the effects of gravity on the microstructure of alkali-activated materials cured with heat, and to develop materials for the construction of long-term infrastructure on the lunar surface with in-situ resource utilization (ISRU). ISRU for construction materials reduces the cost and mass of payloads related to lunar construction. The advantage of geopolymer concrete as opposed to traditional portland cement concrete is that water acts as a medium for the polymerization reaction and leaves the system throughout the process, reducing its demand. Twelve samples of lunar regolith simulant and a solution composed of sodium hydroxide and sodium silicate were sent to the ISS. The three simulants were OPRH2N, OPRL2N, and JSC-1AF, using only particles less than 53 µm in diameter to increase reactivity of the simulant. Simulant to solution ratios were determined by workability while mixing. The simulant and solution were sealed Burst Pouches® along with 2 other sealed bags to prevent material from leaking. Crew member F-14 conducted testing on the ISS by introducing the solution to the simulant in the Burst Pouch®, mixing the sample with a spatula, and then clamping the specimen in the fresh state to prevent flow inside the Burst Pouch®. These specimens were then put in a thermos heated to 80C via sealed drinking water bags to cure for 24 hours with a temperature logger. The cured specimens remained in microgravity for at least 28 days and were returned from the ISS in February 2025. The specimens were then brought to the NASA Marshall Space Flight Center (MSFC) to analyze. Material characterization consisted of conducting Micro-CT tests of entire samples in their sealed apparatus to a resolution of 25µm. 2D image slices were saved in each orthogonal direction of each specimen at a 0.03 mm step size from the 3D model to conduct analytical porosity calculations. Representative samples from each specimen were sampled to perform helium gas pycnometery and were then mounted in resin for SEM imaging, EDS, and nanoindentation. Porosity was analyzed analytically using micromechanics modelling with the assistance of the NASA Multiscale Analysis Tool (NASMAT), as well as the NASA Advanced Supercomputing (NAS) servers (V. Saseendran & N. Yamamoto, 2024). Density was measured using helium gas pycnometery and was then compared to the theoretical density for experimental porosity calculation. Due to the samples’ non-uniform shape being cured in a pouch, traditional compression and tensile strength testing could not be performed. Nanoindentation was conducted at Clarkson University to determine the microhardness and reduced modulus of elasticity. Results from flight samples can be compared to ground samples currently in DLR’s possession to determine the effect on microstructure from being mixed and cured in microgravity. This study gives further insight and understanding of geopolymer lunar concrete and its viability as a lunar construction material with ISRU.
Adam Johnson↗
Effect of Fe on Co-Based SiO2Al2O3 Mixed Support Catalyst for Fischer–Tropsch Synthesis in 3D-Printed SS Microchannel Microreactor
This research explores the effect of a composite support of SiO2 and Al2O3 with Fe and Co incorporated as catalysts for Fischer–Tropsch synthesis (FTS) using a 3D-printed stainless steel (SS) microchannel microreactor. Two mesoporous catalysts, FeCo/SiO2Al2O3 and Co/SiO2Al2O3, were synthesized via a one-pot (OP) method and extensively characterized using N2 physisorption, XRD, SEM, TEM, H2-TPR, TGA-DSC, FTIR, and XPS. H2-TPR results revealed that the synthesis method significantly affected the reducibility of metal oxides, thereby influencing the formation of active FTS sites. SEM-EDS and TEM further revealed a well-defined hexagonal matrix with a porous surface morphology and uniform metal ion distribution. FTS reactions, carried out in the 200–350 °C temperature range at 20 bar with a H2/CO molar ratio of 2:1, exhibited the highest activity for FeCo/SiO2Al2O3, with up to 80% CO conversion. Long-term stability was evaluated by monitoring the catalyst performance for 30 h on stream at 320 °C under identical reaction conditions. The catalyst was initially active for the methanation reaction for up to 15 h, after which the selectivity for CH4 declined. Correspondingly, the C4+ selectivity increased after 15 h of time-on-stream, indicating a shift in the product distribution toward longer-chain hydrocarbons. This trend suggests that the catalyst undergoes gradual activation or restructuring under reaction conditions, which enhances chain growth over time. The increase in C4+ products may be attributed to the stabilization of the active sites and suppression of methane or light hydrocarbon formation.
Biochemistry & Molecular Biology↗
Laboratory Development of a Pot Calcination Process for Converting Liquid Wastes to Solids
Simulated Darex, Purex, Thorex, and TBP-25 wastes were batch-evaporated and calcined to solids. The weight of residue was decreased to ~ 8% of the original waste for Darex, to ~ 4% for Purex, to ~ 2% for Thorex, and ~ 7% for TBP-25 by calcination to 1200°C. Semicontinuous evaporation and calcination to solids at 900°C of Darex and Purex wastes gave volume reduction factors of ~ 8 and 7 to 10, respectively. The nitrate contents of the residues from batch calcination to 800°C were from ~ 0.07 to 0.5 wt % while the nitrate contents of the residues from semicontinuous calcination to 900°C were from ~ 0.1 to 0.7 wt %. Sodium, calcium, and magnesium additives to Purex waste decreased the percent of sulfate volatilized into the condensate from ~ 30% to <0.5%. Replacement of the atmosphere above the waste with nitric oxide decreased ruthenium volatility from 50 to 60% of that originally present in Purex waste to 0.5 to 3.5%. The thermal conductivities for the calcined wastes, measured in situ, were all >0.1 Btu/hr ft °F at ~ 400°F and increased almost linearly with increasing temperature to >0.3 Btu/hr ft °F at ~ 1600°F in all cases studied.
052001* -- Nuclear Fuels-- Waste Processing↗
Anion-dependent phase behavior of methylimidazolium-based ionic liquids mixed with water: Correlation between local molecular structure and mesoscale behaviors
The mesoscopic phase behavior of decylmethylimidazolium (C 10 mim) ionic liquids (ILs) bearing three monovalent anions—thiocyanate (SCN − ), nitrate (NO 3 − ), and chloride (Cl − )—mixed with water at relatively high IL contents (50–95 wt%) was investigated. Small-angle and wide-angle X-ray scattering (SAXS/WAXS) were employed to follow the evolution of both local and mesoscale structures across this composition range. In the absence of water, C 10 mimSCN and C 10 mimCl behaved as disordered liquids, whereas C 10 mimNO 3 spontaneously formed a hexagonally ordered cylindrical mesophase and displayed a sticky-solid macroscopic appearance. Upon addition of water, C 10 mimSCN remained a viscous liquid and only weakly ordered lamellar domains were observed. This limited ordering is attributed to the lack of hydrogen bonding and weak interaction energy between SCN − anions and C 10 mim + cations. In contrast, the trigonal-planar NO 3 − and point-like Cl − anions promoted the formation of well-defined hexagonal mesophases up to 35–45 wt% water. The formation of hydrogen bonding of the two anions with imidazolium ring protons likely enabled the creation of compact ion clusters that effectively exclude water molecules from the immediate vicinity of the IL aggregates. These findings demonstrate that the interaction energy between ion pairs dictate IL-water interactions and therefore control the transition from disordered liquids to ordered mesophases in IL/water mixtures with high IL contents. In conclusion, the combined SAXS/WAXS analysis reveals a correlation between local intermolecular structure and the emergence of mesoscopic order, providing a systematic framework for tailoring mesoscale structures in alkylimidazolium-based IL/water systems.
36 MATERIALS SCIENCE↗
Influences of Introduced Yttrium Oxide Particles on Superalloy 718
Additive manufacturing has proven useful for introducing oxide particles into superalloys, to provide oxide dispersion strengthening (ODS) at high temperatures. This study screened how such an approach can influence microstructure and failure modes for superalloy 718. The objective of this study was to compare tensile and creep failure modes for additively manufactured superalloy 718 having introduced oxide particles. Fine yttrium-oxide (yttria) particles were introduced into 718 powder using two different powder mixing methods. Additive manufacturing by laser powder bed fusion was then used to prepare specimen blanks for each case, oriented parallel and transverse to the building direction. These were subsequently given consistent stress relief, hot isostatic pressurization, and final heat treatments, then subjected to tensile and creep tests at varied temperatures. Additively manufactured 718 without ODS (“no-ODS”) had nearly equiaxed grains at 50 μm in width. Roll-mixed ODS material had a bimodal grain size distribution, with fine grains at 9 μm and coarse grains at 66 μm in average width. The fine grains were elongated about 5x in the building direction. Acoustic-mixed ODS materials had grains at 10 μm in width that were elongated about 10x in the build direction. Compared to no-ODS, roll-mixed and acoustic mixed ODS materials did not show improved tensile strengths at room temperature, 760 °C, and 1093 °C. The tensile failure strains parallel to the building direction of roll-mixed ODS material were lowest in these tests. Flattened clumps of yttria had formed transverse to the building direction in roll-mixed material, which were easily cracked. Creep rupture response at 760 °C in the building direction was highest for acoustic-mixed material test. Creep tests at 1093 °C showed varied results, with low rupture lives associated with enhanced cavitation at grain boundaries and surface oxidation. Transverse to the building direction, all three materials had low failure strains in tests at 760 °C. The elongated grain boundaries failed readily in all transverse specimens tested at 760 °C and 1093 °C.
oxide dispersion strengthening↗
The critical role of intrinsic defects and many-body interactions on the stability of MnBi2Te4
Intrinsic antisite defects pose a major challenge to understanding and predicting the exotic properties of the layered topological magnetic insulator MnBi2Te4 (MBT). In this work, we study the origin of the abundance of intrinsic defects in MBT, including many-body defect–defect interactions and many-body electronic correlations. Until now, ab initio methods have struggled to explain thermodynamic stability and properties influenced by defect behavior in MBT. We model native Mn–Bi antisite defects in MBT at finite temperatures using a cluster expansion that includes defect–defect interactions. To overcome the limitations of conventional density functional theory (DFT), we introduce a hybrid approach that incorporates high-accuracy quantum Monte Carlo (QMC) calculations, introducing missing correlations. This strategy allows for accurate estimation of defect energetics and finite-temperature properties. We compute the configurational free energy, defect concentration, and configurational heat capacity, revealing a second-order order–disorder phase transition near the experimental synthesis temperature. Our study provides the first theoretical insight into the thermodynamics of intrinsic defects in MBT. The negative free energy relative to pristine MBT at synthesis temperatures indicates that Mn–Bi antisite formation is thermodynamically spontaneous. We also present a broadly applicable general framework for correcting low-level theoretical theories using highly accurate many-body corrections from QMC.
Ghaffar, Abdul [ORNL] (ORCID:0000000241190168)↗
Solvent-Mediated Control of Nanocellulose Dispersion: An Integrated Computational and Experimental Investigation
Fibrillated cellulose derived from forestry feedstocks represents a renewable and high-strength materials platform for circular bioeconomies. However, its practical implementation is hindered by the irreversible aggregation of nanocellulose architectures, including cellulose nanofibers (CNFs). Solvent-based dispersion offers a simple and practical route to prevent CNF aggregation. Here, in this work, we integrate classical and enhanced sampling molecular dynamics (MD) simulations with experimental suspension rheology and atomic force microscopy (AFM) to elucidate how solvent environments tune CNF–CNF interactions and dispersion stability. CNF–CNF contact free energies computed from MD simulations reveal reduced aggregation in acetone/water, γ-valerolactone (GVL)/water, and tetrahydrofuran (THF)/water and pure acetone compared with pure water, reflecting stronger CNF-solvent relative to inter-CNF interactions. Correspondingly, CNF-solvent suspensions in these solvent systems exhibit stronger inter-fibril network structures and enhanced recovery compared to water, indicating improved CNF-solvent affinity. Liquid cell AFM imaging in acetone–water mixtures and in pure acetone further confirm the presence of well-dispersed CNFs. By combining multiscale computation with targeted experiments, this study establishes a rational framework for solvent design to achieve stable nanocellulose dispersions for high-strength biobased materials and efficient bioenergy conversion.
cellulose↗