Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale 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.

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

At least 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Unveiling the catalytic mechanism of GTP hydrolysis in microtubules

Microtubules (MTs) are large cytoskeletal polymers, composed of αβ-tubulin heterodimers, capable of stochastically converting from polymerizing to depolymerizing states and vice versa. Depolymerization is coupled with hydrolysis of guanosine triphosphate (GTP) within β-tubulin. Hydrolysis is favored in the MT lattice compared to a free heterodimer with an experimentally observed rate increase of 500- to 700-fold, corresponding to an energetic barrier lowering of 3.8 to 4.0 kcal/mol. Mutagenesis studies have implicated α-tubulin residues, α:E254 and α:D251, as catalytic residues completing the β-tubulin active site of the lower heterodimer in the MT lattice. The mechanism for GTP hydrolysis in the free heterodimer, however, is not understood. Additionally, there has been debate concerning whether the GTP-state lattice is expanded or compacted relative to the GDP state and whether a “compacted” GDP-state lattice is required for hydrolysis. In this work, extensive quantum mechanics/molecular mechanics simulations with transition-tempered metadynamics free-energy sampling of compacted and expanded interdimer complexes, as well as a free heterodimer, have been carried out to provide clear insight into the GTP hydrolysis mechanism. α:E254 was found to be the catalytic residue in a compacted lattice, while in the expanded lattice, disruption of a key salt bridge interaction renders α:E254 less effective. The simulations reveal a barrier decrease of 3.8 ± 0.5 kcal/mol for the compacted lattice compared to a free heterodimer, in good agreement with experimental kinetic measurements. Additionally, the expanded lattice barrier was found to be 6.3 ± 0.5 kcal/mol higher than compacted, demonstrating that GTP hydrolysis is variable with lattice state and slower at the MT tip.

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

Machine-learning modeling of magnetization dynamics in quasi-equilibrium and driven metallic spin systems

Here, we present a perspective on recent progress in machine-learning (ML) force-field approaches for large-scale Landau–Lifshitz–Gilbert (LLG) simulations of metallic spin systems. Building on a generalization of the Behler–Parrinello (BP) architecture originally developed for quantum molecular dynamics, we develop scalable and transferable ML models that faithfully capture the complex, environment-dependent electron-mediated exchange fields characteristic of itinerant magnets. A central ingredient of this framework is the implementation of symmetry-aware magnetic descriptors based on group-theoretical bispectrum formalisms. Leveraging these ML force fields, LLG simulations faithfully reproduce hallmark non-collinear magnetic orders—such as the 120° and tetrahedral states—on the triangular lattice, and successfully capture the complex spin textures emerging in the mixed-phase states of a square-lattice double-exchange model under thermal quench. We further discuss a generalized potential theory that extends the BP formalism to incorporate both conservative and nonconservative electronic torques, thereby enabling ML models to learn nonequilibrium exchange fields from computationally demanding microscopic approaches such as nonequilibrium Green’s-function techniques. This extension yields quantitatively accurate predictions of voltage-driven domain-wall motion and establishes a foundation for quantum-accurate, multiscale modeling of nonequilibrium spin dynamics and spintronic functionalities.

Descriptors↗
DOE OSTI2025

Influences of shock imprinting on mix in a 3D-printed porous media

Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.

2 Photon Polymerization↗
DOE OSTI2024

Driving macro-scale transformations in three-dimensional-printed biopolymers through controlled induction of molecular anisotropy at the nanoscale

Motivated by the need to harness the properties of renewable and biodegradable polymers for the design and manufacturing of multi-scale structures with complex geometries, we have employed our additive manufacturing platform that leverages molecular self-assembly for the production of metre-scale structures characterized by complex geometries and heterogeneous material composition. As a precursor material, we used chitosan, a chemically modified form of chitin, an abundant and sustainable structural polysaccharide. We demonstrate the ability to control concentration-dependent crystallization as well as the induction of the preferred orientation of the polymer chains through the combination of extrusion-based robotic fabrication and directional toolpathing. Anisotropy is demonstrated and assessed through high-resolution micro-X-ray diffraction in conjunction with finite element simulations. Using this approach, we can leverage controlled and user-defined small-scale propagation of residual stresses to induce large-scale folding of the resulting structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Determining the Length Scale of Transport Impedances in Li-Ion Electrodes: Li(Ni 0.33 Mn 0.33 Co 0.33 )O 2

Battery cathodes are complex multiscale, multifunctional materials. The length scale at which the dominant impedance arises may be difficult to determine even with the most advanced experimental characterization efforts, and thus modeling can play an important role in analysis. Discharge and voltage relaxation curves, interrogated with theory, are used to distinguish between transport impedance that arise on the scale of the active crystal and on the scale of agglomerates (secondary particles) comprised of nanoscale crystals. Model-selection algorithms are applied to determine that the agglomerate scale is dominant in the Li Ni 0.33 Mn 0.33 Co 0.33 O 2 electrode studied here. Furthermore, conditions where the agglomerate and crystal-scale models yield distinct simulation results are demonstrated, providing approaches that can be applied to other systems.

25 ENERGY STORAGE↗
DOE OSTI2023

The Role of Topography in Controlling Evapotranspiration Age

Abstract Evapotranspiration (ET) age is a key metric of water sustainability but a major unknown partly due to the extreme difficulty in modeling it. Groundwater is found to be important in ET age variations in small‐scale studies, yet our understanding is insufficient because groundwater systems are nested across scales. Here, we conducted GPU‐accelerated particle tracking with integrated hydrologic modeling to quantify the variations in ET age at a regional scale of ∼0.4 M km 2 . Simulation results reveal topography‐driven flow paths shaping the spatial and temporal patterns of ET age variations. On ridges, where root zone decoupling with deep subsurface storage, ET age is generally young, with seasonal variations dominated by meteorological conditions. In the valley bottom, ET age is generally old, with significant subseasonal variations caused by the convergence of subsurface flow paths. On hillslopes with water table depths ranging from 1 to 10 m, ET age shows strong seasonal variations caused by the connections with lateral groundwater regulated by ET demand. Our modeling approach provides insights into the basic linkages between ET age and topography at large scale. Our work highlights the perspective of multiscale studies of ET age, suggesting new field experiments to test these process connections and to determine if such linkages warrant inclusion in Earth System Models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A novel conditional formulation of the Vlasov–Ampère equations: a conservative, positivity, asymptotic and Gauss law preserving scheme

We propose a novel reformulation of the Vlasov–Ampère equations for plasmas that reveals discrete symmetries that enables simultaneous conservation of mass, momentum and energy; preservation of Gauss’s law; positivity of the distribution function; and consistency with quasi-neutral asymptotics. The approach employs variable and coordinate transformations to yield a coupled system comprising a modified Vlasov equation and associated moment–field equations. The modified Vlasov equation advances a conditional distribution function that excludes mass, momentum and energy densities, which are instead evolved through moment equations enforcing the relevant symmetries, conservation laws and involution constraints. This reformulation aligns naturally with a recent slow-manifold reduction technique, which separates fast electron time scales and simplifies the treatment of the quasi-neutral limit within the reduced moment–field subsystem. Using this framework, we develop a numerical method for the reduced 1D1V subsystem that, for the first time in the literature, satisfies all key physical constraints while maintaining a quasi-neutral asymptotic behaviour. The advantages of the method are demonstrated on canonical electrostatic test problems, including the multiscale ion acoustic shock wave.

1D1V↗
DOE OSTI2023

Model America - Arizona extract from ORNL's AutoBEM v1.1

Oak Ridge National Laboratory (ORNL) has developed the Automatic Building Energy Modeling (AutoBEM) software suite to process multiple types of data, extract building-specific descriptors, generate building energy models, and simulate them on High Performance Computing (HPC) resources. For more information, see AutoBEM-related publications (bit.ly/AutoBEM).Two sets of sample data are provided for 2,555,152 buildings located within the boundary of Arizona in the United States:Data (846.3MB *.csv) - minimalist list of each building (rows) for the following fields (columns) • ID - unique building ID • Centroid - building center location in latitude/longitude (from Footprint2D) • Footprint2D - building polygon of 2D footprint (lat1/lon1_lat2/lon2_...) • State_abbr - state name • Area - estimate of total conditioned floor area (ft2) • Area2D - footprint area (ft2) • Height - building height (ft) • NumFloors - number of floors (above-grade) • WWR_surfaces - percent of each facade (pair of points from Footprint2D) covered by fenestration/windows (average 14.5% for residential, 40% for commercial buildings) • CZ - ASHRAE Climate Zone designation • BuildingType - DOE prototype building designation (IECC=residential) as implemented by OpenStudio-standards • Standard - building vintage • Sample Models (114GB*.zip by county) - OpenStudio and EnergyPlus building energy models named according to IDThis data is made free and openly available in hopes of stimulating any simulation-informed use case. Data is provided as-is with no warranties, express or implied, regarding fitness for a particular purpose. We wish to thank our sponsors which include Oak Ridge National Laboratory (ORNL), U.S. Dept. of Energy’s (DOE) Building Technologies Office (BTO), Office of Electricity (OE), and Biological and Environmental Research (BER).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

SBML Level 3: an extensible format for the exchange and reuse of biological models

Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction-based models and packages that extend the core with features suited to other model types including constraint-based models, reaction-diffusion models, logical network models, and rule-based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single-cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Automated Multiscale Approach To Predict Self-Diffusion from a Potential Energy Field

For large-scale screening studies there is a need to estimate the diffusion of gas molecules in nanoporous materials more efficiently than (brute force) molecular dynamics. In particular for systems with low diffusion coefficients molecular dynamics can be prohibitively expensive. An alternative is to compute the hopping rates between adsorption sites using transition state theory. For large-scale screening this requires the automatic detection of the transition states between the adsorption sites along the different diffusion paths. Here an algorithm is presented that analyzes energy grids for the moving particles. It detects the energies at which diffusion paths are formed, together with their directions. This allows for easy identification of nondiffusive systems. For diffusive systems, it partitions the grid coordinates assigned to energy basins and transitions states, permitting a transition state theory based analysis of the diffusion. We test our method on CH 4 diffusion in zeolites, using a standard kinetic Monte Carlo simulation based on the output of our grid analysis. We find that it is accurate, fast, and rigorous without limitations to the geometries of the diffusion tunnels or transition states.

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

xSDK: Building an ecosystem of highly efficient math libraries for exascale

Current efforts to build increasingly powerful computer architectures are opening up new avenues for more complex and higher fidelity simulations coupled with data analytics and learning, leading to new scientific insights and deeper understanding. At one extreme, exascale computers will be much faster than previous computer generations (performing 10 18 operations per second—that is, 1,000 times faster than petascale). To achieve these performance improvements, computer architectures are becoming increasingly complex, with deep memory hierarchies, very high node and core counts, and heterogeneous features such as graphics processing units (GPUs). Such architectural changes impact the full breadth of computing scales, as heterogeneity pervades even current-generation laptops, workstations, and moderate-sized clusters. While emerging advanced architectures provide unprecedented opportunities, they also present significant challenges for developers of scientific applications, such as multiphysics and multiscale codes, who must adapt their software to handle disruptive changes in architectures and new programming models that have not yet stabilized. Developers must consider increasing concurrency while reducing communication and synchronization, and other complexities such as the potential for using mixed precision to leverage the compute power available in low-precision tensor cores. On one hand, developers must implement new scientific capabilities, which in turn increase code complexity. On the other hand, the codes must be ported to new architectures, requiring the inclusion of new programming models and the restructuring of code to achieve good performance. Addressing these issues is beyond the capability of any single person or team—leading to the need for collaboration among many teams, who encapsulate their expertise in reusable software and work together to create sustainable software ecosystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Chemo-mechanical-microstructural coupling in the tarsus exoskeleton of the scorpion Scorpio palmatus

The multiscale structure of biomaterials enables their exceptional mechanical robustness, yet the impact of each constituent at their relevant length scale remains elusive. Here, we used SAXD analysis to expose the intact chitin-fiber architecture within the exoskeleton on a scorpion's claw, revealing varying orientations, including Bouligand and unidirectional regions different from other arthropod species. We uncovered the contribution of individual components’ constituent behavior to its mechanical properties from the micro- to the nanoscale. At the microscale, in-situ micromechanical experiments were used to determine site-specific stiffness, strength, and failure of the biocomposite due to fiber orientation, while metal-crosslinking of proteins is characterized via fluorescence maps. At the constituent level, combined with FEA simulations, we uncovered the behavior of fiber-matrix deformation with fiber diameter <53.7 nm and protein modulus in the range 1.4–11 MPa. The unveiled microstructure-mechanics relationship sheds light on the evolved structural functionalities and constituents' interactions within the scorpion cuticle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Anisotropic yield models for lattice unit cell structures exploiting orthotropic symmetry

We report numerical homogenization enables efficient computational analysis and design of multiscale structures made of micro-architected materials including lattice unit cells. However, predicting yield is nontrivial because it requires accurate and efficient predictions of the maximum stress inside the homogenized unit cells. To address this challenge, we develop a macroscale anisotropic yield function for micro-architected materials. The yield function depends on the three-dimensional macroscale stress state and the parameters describing a family of micro-architectures, such as the radii of the struts in a lattice unit cell. To ensure accuracy, we determine the maximum stress using three-dimensional continuum finite-element analysis. To ensure efficiency, we construct surrogate models from the aforementioned high-fidelity results for yield prediction. To reduce simulation costs and surrogate modeling complexity, we leverage orthotropic symmetry commonly found in lattice unit cells. In this paper, we provide a thorough presentation of group representations and the systematic procedure to exploit orthotropic domain symmetry in homogenization and surrogate modeling. We illustrate the use of linear homogenization results to predict yield in specific unit cells without further simulations. We furthermore show that surrogate modeling presents a viable option for anisotropic yield prediction in a continuously parametrized family of micro-architectures. More specifically, despite the dimensionality and degree of nonlinearity in the maximum micro von Mises stress within the unit cell, the surrogate models can predict it with less than 5% error at least 90% of the time. Moreover, the largest under-prediction error, which is more critical than the over-prediction error, is typically less than 10%.

42 ENGINEERING↗
DOE OSTI2019

ATF material model development and validation for priority fuel concepts

Development, implementation, and validation of material and behavior models for accident tolerant fuel (ATF) concepts in the Bison fuel performance code began in 2014 in response to the events that occurred at the Fukushima Daichii nuclear power plant in March 2011. Early on the focus was on U 3 Si 2 fuel and FeCrAl cladding as part of a high impact problem through the Nuclear Energy and Advanced Modeling Simulation (NEAMS) program. Then, developments for Cr 2 O 3 -doped UO 2 fuel, and SiC-SiC and Cr-coated zirconium-based claddings began based upon industry interests. In late fiscal year 2018 the Consortium for Advanced Simulation of Light Water Reactors (CASL) took over further ATF work in Bison in support of the Nuclear Regulatory Commission (NRC) engagement. Discussions with the NRC identified their list of priority fuel and cladding concepts, which included Cr 2 O 3 -doped UO 2 and U 3 Si 2 fuels, and Cr-coated zirconium-based and FeCrAl claddings. In particular, the NRC suggested that reports similar in form to NUREG/CR-7024 [1] that was developed for traditional LWR materials UO 2 and zirconium-based claddings (i.e., Zircaloy-4, M5®, ZIRLO) be created for the priority ATF concepts. The approach to ATF capability development in Bison since the beginning has been two-fold: (1) empirical correlations and (2) multiscale model development. Both approaches have uncertainty inherent to them. Uncertainty in empirical correlations is bounded by the experimental data upon which with the correlation was developed. Models developed through a multiscale approach have uncertainty associated with the lower length scale calculations and input parameters that must be propagated to the engineering scale model in Bison. In this report, the recommended models, their range of applicability (e.g., temperature, burnup), and associated uncertainty for the NRC priority fuel concepts Cr 2 O 3 -doped UO 2 and U 3 Si 2 are presented in a manner similar to the aforementioned NUREG. In addition, the Cr 2 O 3 -doped UO 2 models are validated to the Halden Reactor tess IFA-677.1 rods 1 and 5 and IFA-716.1 rod 1. For U 3 Si 2 the models are validated to the recent post-irradiation examination (PIE) data ATF-13 and ATF-15 rods that were irradiated in the Advanced Test Reactor (ATR) at Idaho National Laboratory. Finally, an uncertainty quantification (UQ) and sensitivity analysis (SA) is completed on a single rod for each fuel concept that takes into account the defined uncertainty in select models. The results of the UQ and SA analyses indicate that given the large uncertainty in some of the input ATF models that the uncertainty on the Bison predictions of fuel performance metrics of interest bound the available experimental data. Correlation coefficients are also reported that identified the uncertain inputs with the strongest correlation (positive or negative) with the outputs of interest. Fuel performance metrics investigated include fuel elongation, rod internal pressure, fuel centerline temperature, and fission gas release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

ATF material model development and validation for priority cladding concepts

Development, implementation, and validation of material and behavior models for accident tolerant fuel (ATF) concepts in the Bison fuel performance code began in 2014 in response to the events that occurred at the Fukushima Daichii nuclear power plant in March 2011. Early on the focus was on U 3 Si 2 fuel and FeCrAl cladding as part of a high impact problem through the Nuclear Energy and Advanced Modeling Simulation (NEAMS) program. Then, developments for Cr 2 O 3 -doped UO 2 fuel, and SiC-SiC and Cr-coated zirconium-based claddings began based upon industry interests. In late fiscal year 2018 the Consortium for Advanced Simulation of Light Water Reactors (CASL) took over further ATF work in Bison in support of the Nuclear Regulatory Commission (NRC) engagement. Discussions with the NRC identified their list of priority fuel and cladding concepts, which included Cr 2 O 3 -doped UO 2 and U 3 Si 2 fuels, and Cr-coated zirconium-based and FeCrAl claddings. In particular, the NRC suggested that reports similar in form to NUREG/CR-7024 that was developed for traditional LWR materials UO 2 and zirconium-based claddings (i.e., Zircaloy-4, M5 ® , ZIRLO TM ) be created for the priority ATF concepts. The approach to ATF capability development in Bison since the beginning has been two-fold: (1) empirical correlations and (2) multiscale model development. Both approaches have uncertainty inherent to them. Uncertainty in empirical correlations is bounded by the experimental data upon which with the correlation was developed. Models developed through a multiscale approach have uncertainty associated with the lower length scale calculations and input parameters that must be propagated to the engineering scale model in Bison. In this report, the recommended models, their range of applicability (e.g., temperature, fluence), and associated uncertainty for the NRC priority clad concepts Cr-coated zirconium based clad and FeCrAl are presented in a manner similar to the aforementioned NUREG. In addition, in the absence of experimental data for validation, demonstration cases are presented that investigate the ballooning and burst behavior of Cr-coated and FeCrAl clads as well as tritium permeation in FeCrAl. A portion of the study by Gamble et al. on FeCrAl failure is also revisited with the latest updated models. First, the nominal results are shown using the as implemented models. Then a sensitivity analysis and uncertainty quantification is performed to assess the impact of the uncertainties in the models on the parameters of interest. For the Cr-clad ballooning case the parameter of interest is the time to burst. For the revisited FeCrAl study, the parameter of interest is the pressure at burst.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Demonstrate Improved Ag Diffusion and Describe the Basis for Pd Penetration Modeling in SiC

In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Solving the master equation without kinetic Monte Carlo: Tensor train approximations for a CO oxidation model

In multiscale modeling of heterogeneous catalytic processes, one crucial point is the solution of a Markovian master equation describing the stochastic reaction kinetics. Usually, this is too high-dimensional to be solved with standard numerical techniques and one has to rely on sampling approaches based on the kinetic Monte Carlo method. In this study we break the curse of dimensionality for the direct solution of the Markovian master equation by exploiting the Tensor Train Format for this purpose. The performance of the approach is demonstrated on a first principles based, reduced model for the CO oxidation on the RuO{sub 2}(110) surface. We investigate the complexity for increasing system size and for various reaction conditions. The advantage over the stochastic simulation approach is illustrated by a problem with increased stiffness.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Deconstruction of biomass enabled by local demixing of cosolvents at cellulose and lignin surfaces

A particularly promising approach to deconstructing and fractionating lignocellulosic biomass to produce green renewable fuels and high-value chemicals pretreats the biomass with organic solvents in aqueous solution. Here, neutron scattering and molecular-dynamics simulations reveal the temperature-dependent morphological changes in poplar wood biomass during tetrahydrofuran (THF):water pretreatment and provide a mechanism by which the solvent components drive efficient biomass breakdown. Whereas lignin dissociates over a wide temperature range (>25 °C) cellulose disruption occurs only above 150 °C. Neutron scattering with contrast variation provides direct evidence for the formation of THF-rich nanoclusters (R g ∼ 0.5 nm) on the nonpolar cellulose surfaces and on hydrophobic lignin, and equivalent water-rich nanoclusters on polar cellulose surfaces. The disassembly of the amphiphilic biomass is thus enabled through the local demixing of highly functional cosolvents, THF and water, which preferentially solvate specific biomass surfaces so as to match the local solute polarity. A multiscale description of the efficiency of THF:water pretreatment is provided: matching polarity at the atomic scale prevents lignin aggregation and disrupts cellulose, leading to improvements in deconstruction at the macroscopic scale.

09 BIOMASS FUELS↗