SEARCH · Engineering Papers

Results for “multiple time stepping”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Distributed model predictive control for coordinated, grid-interactive buildings

Continued focus on reducing carbon emissions and improving energy efficiency requires buildings to become grid-interactive and not just behave as static consumers. A distributed model predictive control (DMPC) algorithm known as Limited-Communication (LC) DMPC is modified to enable grid-interactive buildings. A grid-aggregator subsystem is added that allows for a bulk grid power reference signal to be followed while the individual building subsystems also achieve their local comfort objectives. The LC-DMPC algorithm is applied for the first time to systems with multiple buildings. Adequate power tracking is shown for different simulation scenarios involving heterogeneous buildings, and next steps are discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

COVID-19 dynamics across the US: A deep learning study of human mobility and social behavior

This paper presents a deep learning framework for epidemiology system identification from noisy and sparse observations with quantified uncertainty. The proposed approach employs an ensemble of deep neural networks to infer the time-dependent reproduction number of an infectious disease by formulating a tensor-based multi-step loss function that allows us to efficiently calibrate the model on multiple observed trajectories. The method is applied to a mobility and social behavior-based SEIR model of COVID-19 spread. The model is trained on Google and Unacast mobility data spanning a period of 66 days, and is able to yield accurate future forecasts of COVID-19 spread in 203 US counties within a time-window of 15 days. Interestingly, a sensitivity analysis that assesses the importance of different mobility and social behavior parameters reveals that attendance of close places, including workplaces, residential, and retail and recreational locations, has the largest impact on the effective reproduction number. Furthermore, the model enables us to rapidly probe and quantify the effects of government interventions, such as lock-down and re-opening strategies. Taken together, the proposed framework provides a robust workflow for data-driven epidemiology model discovery under uncertainty and produces probabilistic forecasts for the evolution of a pandemic that can judiciously provide information for policy and decision making. All codes and data accompanying this manuscript are available at https://github.com/PredictiveIntelligenceLab/DeepCOVID19.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Multimodal Analysis of Reaction Pathways of Cathode Materials for Lithium Ion Batteries

Conversion mechanism in lithium ion batteries provides higher capacity than intercalation mechanism since multiple numbers of electrons and lithium ions are associated. However, poor cycling stability, large voltage hysteresis, and low energy efficiency have been great challenges of conversion reaction. To address those issues, conversion-type electrode materials have been reformed via doping or substituting other elements. For example, iron fluorides (FeF 2 , FeF 3 ) have modified as iron oxyfluorides (FeF 1-x O x ), showing enhanced long-term stability. Furthermore, co-substituted (both anion and cation substituted) Fe 0.9 Co 0.1 OF (FeCoOF) was demonstrated excellent cycling stability. (capacity of 350 mAh g -1 at a current of 500 mA g -1 for 1000 cycles). Substituting anion and cation in iron fluoride has been suggested as an effective method to achieve better reversibility but understanding of lithiation reactions in co-substituted FeCoOF is not clear. This work takes advantage of ex-situ/ in-situ synchrotron X-ray based techniques and transmission electron microscopy to elucidate structural changes with lithium ion, which may provide fundamental insights into modifying conversion-type materials. Figure 1 presents discharge-charge curves and pair distribution function patterns acquired at each potential. As lithium ions were inserted, structural changes were noticed both at short-range and long-range. However, long-range ordering was nearly maintained even at 1 V, indicating absence of conversion reaction. Figure 2 shows lithiation induced structural evolution of a single FeCoOF nanorod observed in real time. As lithiation proceeds, the width of nanorods shows a stepwise increase, particularly A in figure 2d, which may indicate multiple steps of lithiation occur. Considering that conversion reaction takes place around 2 V in FeF 3 , co-substitution Co and O into iron fluoride may change thermodynamic features of lithiation reactions by lowering the initiation potential for conversion reaction. Instead, phase transformations occur at long-range order, which may help maintaining structural integrity during operation, eventually, achieving cycling stability.

25 ENERGY STORAGE↗

DOE OSTI2026

Postirradiation Examination of WIRE-21 Experiment Irradiated in the High Flux Isotope Reactor

Westinghouse Electric Company is developing wireless sensors to monitor the centerline temperature and internal pressure of commercial light-water reactor fuel rods during irradiation. Oak Ridge National Laboratory and Westinghouse Electric Company developed the Wireless Instrumented RB Experiment 2021 (WIRE-21) to test wireless temperature and pressure sensor technologies in a removable beryllium position in the High Flux Isotope Reactor. The experiment was irradiated for a total of 75 days, at temperatures ranging from approximately 150°C to 400°C, resulting in a peak fast (energy > 0.1 MeV) neutron fluence of about 3 × 10 21 n/cm 2 . The temperature was intentionally cycled multiple times to compare the response of the wireless temperature sensor to collocated thermocouples. Similarly, the pressure sensor was actuated in multiple steps to compare the response of the wireless measurement to excore pressure transducers (Petrie et al., 2023). After irradiation, the experiment was disassembled in the Irradiated Fuels Examination Laboratory (IFEL) hot cell at Oak Ridge National Laboratory with the intent to recover sensor and dosimetry components, document the as-irradiated condition of the hardware, and investigate possible causes of the sensor behavior observed during irradiation. The postirradiation examination successfully recovered and preserved key WIRE-21 components. After the housing was removed using a milling machine, the internal experiment sections were examined. All eight fiber-optic sensors were recovered, cut, and stored. The silicon carbide thermometry, temperature sensor, pressure sensor, lower spacer, and selected pressure and temperature cable sections were also removed and stored. The metal bellows of the pressure sensor was found to be plastically deformed, indicating that it had been properly pressurized during irradiation and generally behaved as expected. X-ray diffraction analysis of a section of one of the irradiated inductor cores within the pressure sensor was performed and confirmed the presence of phase-pure alpha ferrite (i.e., no unexpected phase transformations). Inductance testing was performed on the irradiated pressure sensor cores using an unirradiated test coil, and DC resistance measurements of the transceiver coils were also performed. The measurements with irradiated inductor cores assembled inside unirradiated coils showed slightly reduced inductance compared to measurements made with unirradiated cores, but the difference was not sufficient to explain the more significant reductions in inductance that were observed in-pile. Therefore, it is suspected that degradation of the inductor coils (specifically the wire wrapping) is responsible for the reduced inductance observed in-pile.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

PARA: A New Platform for the Rapid Assembly of gRNA Arrays for Multiplexed CRISPR Technologies

Multiplexed CRISPR technologies have great potential for pathway engineering and genome editing. However, their applications are constrained by complex, laborious and time-consuming cloning steps. In this research, we developed a novel method, PARA, which allows for the one-step assembly of multiple guide RNAs (gRNAs) into a CRISPR vector with up to 18 gRNAs. Here, we demonstrate that PARA is capable of the efficient assembly of transfer RNA/Csy4/ribozyme-based gRNA arrays. To aid in this process and to streamline vector construction, we developed a user-friendly PARAweb tool for designing PCR primers and component DNA parts and simulating assembled gRNA arrays and vector sequences.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance

This study focuses on the performance of the elastic–viscous–plastic (EVP) dynamical solver within the sea ice model, CICE v6.5.1. The study has been conducted in two steps. First, the standard EVP solver was extracted from CICE for experiments with refactored versions, which are used for performance testing. Second, one refactored version was integrated and tested in the full CICE model to demonstrate that the new algorithms do not significantly impact the physical results. The study reveals two dominant bottlenecks, namely (1) the number of Message Parsing Interface (MPI) and Open Multi-Processing (OpenMP) synchronization points required for halo exchanges during each time step combined with the irregular domain of active sea ice points and (2) the lack of single-instruction, multiple-data (SIMD) code generation. The standard EVP solver has been refactored based on two generic patterns. The first pattern exposes how general finite differences on masked multi-dimensional arrays can be expressed in order to produce significantly better code generation by changing the memory access pattern from random access to direct access. The second pattern takes an alternative approach to handle static grid properties. The measured single-core performance improvement is more than a factor of 5 compared to the standard implementation. The refactored implementation of strong scales on the Intel® Xeon® Scalable Processors series node until the available bandwidth of the node is used. For the Intel® Xeon® CPU Max series, there is sufficient bandwidth to allow the strong scaling to continue for all the cores on the node, resulting in a single-node improvement factor of 35 over the standard implementation. This study also demonstrates improved performance on GPU processors.

58 GEOSCIENCES↗

DOE OSTI2023

Writing the Programs of Programmable Catalysis

It has long been known that non-steady state and periodic catalytic reactor operation in terms of temperature, pressure, and composition can lead to higher overall productivity and/or product selectivity than the best steady operation. Recently, the emergence of catalysts whose intrinsic properties can be made to oscillate with time, introduces advanced forcing capabilities that can be “programmed” into the catalysts to broaden the scope and applicability of periodic operation to surface chemistry. In this work, an algorithmic approach is implemented to significantly accelerate the discovery and optimization of periodic steady states of catalytic reactors. Decomposition of complex dynamics into fundamental mechanistic fast–slow steps is seen to improve conceptual understanding of the relationship between binding energy oscillation protocols and overall catalytic rates. Finding structured forcing protocols, optimally tailored to the multiple time scales of a given individual mechanism, requires an efficient search of high-dimensional parameter spaces. Here, this is enabled here through active learning (Bayesian optimization, enhanced by our proposed Bayesian continuation). Implementation of these methods is shown to accelerate the evaluation of catalyst programs by up to several orders of magnitude. Faster screening of programmable catalysts to discover periodic steady states enables the optimization of catalytic operating protocols and thus opens the possibility for catalyst engineering based on optimal forcing programs to control rate and product selectivity, even for complex multistep catalytic mechanisms.

catalysis↗

DOE OSTI2025

Macro-level mechanical interlocking: A rapid joining approach for additively manufactured compression molded composite panels

Composite joining typically involves multiple steps, such as drilling and surface treatment, as part of the manufacturing process, which leads to low throughput and long cycle times. In the present study, we demonstrated a macro-level mechanical interlocking (MI) based, rapid joining technique to assemble additively manufactured compression molded (AMCM) panels, enabling the production of parts larger than the mold dimensions. Composite panels made of 20 wt% short carbon fiber reinforced acrylonitrile butadiene styrene (CF/ABS) were joined using MI features of various geometries, namely tree (TR), dovetail (Dov), rectangle 2 (Rect2), and rectangle 1 (Rect1), and their in-plane strength was evaluated. The resultant strength of the tested MI joints reached up to 74 % of the baseline tensile strength (i.e., the ‘no joint’ case). Observations from optical and scanning electron microscopy revealed inadequate polymer diffusion between the adherends, indicating that the joint strength was primarily derived from mechanical interlocking. Additionally, the fracture surfaces exhibited stress-whitening marks, which were characterized using differential scanning calorimetry (DSC). The increase in melting enthalpy suggested local stretching of polymer chains due to MI. Finite element analysis (FEA) indicated that the Rect1 MI feature, which generated the lowest stress concentration, outperformed the others in terms of joint strength, achieving 42 MPa. As a demonstration of the MI joining method, a battery box tray measuring 108 cm × 34 cm using a mold with an effective dimension of 36 cm × 34 cm successfully manufactured, resulting in a part with an area three times larger than the mold. In conclusion, this study presents a promising approach to improving composite joining techniques while minimizing production complexities.

In-plane joining↗

DOE OSTI2021

High-Quality Revision of the Israeli Seismic Bulletin

Seismic bulletins, with trustworthy phase picks, origin times, and source locations are key for regional seismic studies, such as travel-time (TT) tomography, attenuation tomography, and anisotropy studies. To lay the groundwork for such studies in Israel, we revised the seismic bulletin of Israel and the surrounding area and obtained a trustworthy TT data set. From the earthquake and explosion bulletins of the Geophysical Institute of Israel, we compiled a starting data set of about 123,000 earthquakes and explosions that occurred during the past 40 yr. After screening out the poorly recorded events, we were left with a data set of ~38,000 well-recorded events. We then revised the remaining data set in two consecutive steps. In the first, we reviewed and updated station metadata, including changes in station metadata parameters over time. In the second step, we jointly relocated a list of selected seismic events, using the Bayesian hierarchical location software package (BayesLoc) of Myers et al. (2007) that performs joint relocation of multiple events. We observed striking dissimilarities between the spatial distributions of the newly relocated catalog and the initial locations. Although the depth distribution of the starting catalog is trimodal with peaks at 0, 5, and 10 km, the distribution in this study is unimodal, with a broad peak between 7.5 and 12.5 km. By differencing the observed arrival times and the origin times obtained through relocation with BayesLoc, we obtained a revised TT database that consists of 261,336 Pg, 132,876 Pn, 114,816 Sg, and 60,394 Sn arrivals, from a set of 30,458 jointly relocated seismic sources. In this work, we compared prerevision and postrevision TTs as a function of epicentral distance and concluded that the revised data set contains far fewer outliers and inconsistencies than the original data set. The revised TT data set may be used for seismic studies, such as TT tomography, attenuation tomography, and anisotropy studies.

58 GEOSCIENCES↗

DOE OSTI2020

Microscale mechanical modeling of deformable geomaterials with dynamic contacts based on the numerical manifold method

Abstract Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.

58 GEOSCIENCES↗

DOE OSTI2021

Study of the niobium oxide structure and microscopic effect of plasma processing on the Nb surface

A study of the niobium oxide structure is presented here, focusing on the niobium suboxides. Multiple steps of argon sputtering and XPS measurements were carried out until the metal surface was exposed. Subsequently, the sample was exposed to air for different time intervals and the oxide regrowth was studied. In addition, three Nb samples prepared with different surface treatments were studied before and after being subjected to plasma processing. The scope is investigating the microscopic effect that the reactive oxygen contained in the glow discharge may have on the niobium surface. This study suggests that the Nb 2 O 5 thickness may increase. Nevertheless, since the Nb 2 O 5 is dielectric, its thickening would not negatively affect the cavity performance.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Study of the Niobium Oxide Structure and Microscopic Effect of Plasma Processing on the Nb Surface

36 MATERIALS SCIENCE↗

DOE OSTI2025

National Energy Water Treatment & Speciation (NEWTS): A Water & Critical Mineral Database and Dashboard

The scarcity of water resources, the need for beneficial water reuse, and the challenges of wastewater treatment are becoming increasingly pressing in economic, social, and environmental domains. Addressing these concerns requires effective treatment strategies to manage wastewater streams and tackle environmental and economic issues. Furthermore, the recovery of critical minerals from the waste streams associated with energy production holds the promise of offsetting treatment costs and securing local sources of valuable minerals. However, relevant data on these waste streams are dispersed and challenging to locate. The process of ingesting such data into modeling software often involves multiple steps, requiring data restructuring to meet software-input requirements. The non-standardized reporting of water data makes data aggregation and reformatting a time-consuming process. Additionally, essential attributes necessary for modeling water treatment and mineral scale formation are frequently missing. Moreover, data gaps vary depending on the region of interest. Consequently, there is a pressing need for high-quality energy-water composition data that can be easily imported into water chemistry modeling software. To address this need, the National Energy Technology Laboratory has created the National Energy Water Treatment and Speciation (NEWTS) Database and Dashboard—a free online tool catering to community leaders and water researchers. NEWTS facilitates a comprehensive understanding of the composition of energy-related wastewater streams in the United States. The datasets provide detailed concentrations and speciation of major and minor aqueous compounds in energy-related wastewater streams, including power plant leachate, acid mine drainage, brackish water, and oil and gas produced water across the United States. Many of the aqueous species are critical minerals (Li, REEs) in high demand to modernize the world’s energy infrastructure. Many of the datasets also contain volumetric flow-rates needed to model the treatment and reuse scenarios in advanced aqueous chemistry software programs. The NEWTS Database and Dashboard offer public access to hitherto challenging-to-access datasets, presented in a standardized format that is tailored for easy input into aqueous chemistry modeling software. By performing the work needed to transform dispersed, disparate data sources into unified, model-ready datasets, NEWTS serves as an essential resource in advancing water treatment research and sustainable water resource management.

produced water management↗

DOE OSTI2021

BETO 2021 Peer Review - WBS 4.1.2.32: Bioeconomy Scenario Analysis

The Bioeconomy Scenario Analysis project uses systems thinking and analysis to assess current and/or prospective techno-economics, research and development, deployment strategies, policy, and market conditions and their impact on the potential development trajectories of the bioenergy industry over time. Results from this project include identification of opportunities and constraints to industrial development, quantification of multiple metrics (energy, economic, environmental) and informing researchers, decision makers, and industry of the steps needed for a sustainable, nationwide biofuels industry. Analyses from this project enable the creation of a bioenergy industry by (1) inciting policy-makers to explore scenarios for nationwide biofuels production, identifying policy actions consistent with pathways for growth; (2) improving industry’s understanding of industry growth potential under different technology and investment conditions, better targeting their development efforts; and (3) providing universities and other interested stakeholders with tools and analyses that can be adapted to meet research and teaching objectives, connecting students with careers that build the industry. One of the many modeling tools used in this project, the Biomass Scenario Model (BSM) is a publicly-available, unique, validated, state-of-the-art, award-winning, fourth-generation model of the domestic biofuels supply chain which explicitly focuses on how and under what conditions biofuel technologies might be deployed to contribute to the U.S. transportation energy sector. We use models like the BSM to examine the implications of policies and incentives as well as their potential side-effects. The BSM uses a system-dynamics simulation to model dynamic interactions and transitions across the supply chain; it tracks the deployment of biofuels given industrial learning and the reaction of the investment community in the context of land availability, projected oil markets, consumer demand for biofuels, and government policies over time. Under expected market conditions, analyses using the BSM suggest that the biofuels industry may require significant external actions in the early years to thrive. Interventions that accelerate the industrial learning process (e.g. operation of pre-commercial and commercial facilities) have been identified as having strong influence in starting the growth of a commercial biofuel industry. Policies which are coordinated across the whole supply chain in BSM foster the growth of the biofuels industry and production of tens of billions of gallons of biofuels may occur under sufficiently favorable conditions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2022

QuYBE - An Algebraic Compiler for Quantum Circuit Compression

QuYBE is an open-source algebraic compiler for the compression of quantum circuits. It has been applied for the efficient simulation of the Heisenberg Hamiltonian on quantum computers. Currently, it can simulate the time dynamics of one-dimensional chains. It includes modules to generate the quantum circuits for the above as well as produce the compressed circuits, which are independent of the time step. It utilizes the Yang-Baxter equation (YBE) to perform the compression. QuYBE enables users to seamlessly design, execute, and analyze the time dynamics of the Heisenberg Hamiltonian on quantum computers. QuYBE is the first step toward making the YBE technique available to a broader community of scientists from multiple domains. The QuYBE compiler is available at https://github.com/ZichangHe/QuYBE.

Gulania, Sahil↗

DOE OSTI2024

GCAM–GLORY v1.0: representing global reservoir water storage in a multi-sector human–Earth system model

Abstract. Reservoirs play a significant role in modifying the spatiotemporal availability of surface water to meet multi-sector human demands, despite representing a relatively small fraction of the global water budget. Yet the integrated modeling frameworks that explore the interactions among climate, land, energy, water, and socioeconomic systems at a global scale often contain limited representations of water storage dynamics that incorporate feedbacks from other systems. In this study, we implement a representation of water storage in the Global Change Analysis Model (GCAM) to enable the exploration of the future role (e.g., expansion) of reservoir water storage globally in meeting demands for, and evolving in response to interactions with, the climate, land, and energy systems. GCAM represents 235 global water basins, operates at 5-year time steps, and uses supply curves to capture economic competition among renewable water (now including reservoirs), non-renewable groundwater, and desalination. Our approach consists of developing the GLObal Reservoir Yield (GLORY) model, which uses a linear programming (LP)-based optimization algorithm and dynamically linking GLORY with GCAM. The new coupled GCAM–GLORY approach improves the representation of reservoir water storage in GCAM in several ways. First, the GLORY model identifies the cost of supplying increasing levels of water supply from reservoir storage by considering regional physical and economic factors, such as evolving monthly reservoir inflows and demands, and the leveled cost of constructing additional reservoir storage capacity. Second, by passing those costs to GCAM, GLORY enables the exploration of future regional reservoir expansion pathways and their response to climate and socioeconomic drivers. To guide the model toward reasonable reservoir expansion pathways, GLORY applies a diverse array of feasibility constraints related to protected land, population, water sources, and cropland. Finally, the GLORY–GCAM feedback loop allows evolving water demands from GCAM to inform GLORY, resulting in an updated supply curve at each time step, thus enabling GCAM to establish a more meaningful economic value of water. This study improves our understanding of the sensitivity of reservoir water supply to multiple physical and economic dimensions, such as sub-annual variations in climate conditions and human water demands, especially for basins experiencing socioeconomic droughts.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Hamiltonian learning using machine-learning models trained with continuous measurements

Here, we build upon recent work on the use of machine-learning models to estimate Hamiltonian parameters using continuous weak measurement of qubits as input. We consider two settings for the training of our model: (1) supervised learning, where the weak-measurement training record can be labeled with known Hamiltonian parameters, and (2) unsupervised learning, where no labels are available. The first has the advantage of not requiring an explicit representation of the quantum state, thus potentially scaling very favorably to a larger number of qubits. The second requires the implementation of a physical model to map the Hamiltonian parameters to a measurement record, which we implement using an integrator of the physical model with a recurrent neural network to provide a model-free correction at every time step to account for small effects not captured by the physical model. We test our construction on a system of two qubits and demonstrate accurate prediction of multiple physical parameters in both the supervised context and the unsupervised context. We demonstrate that the model benefits from larger training sets, establishing that it is “learning,” and we show robustness regarding errors in the assumed physical model by achieving accurate parameter estimation in the presence of unanticipated single-particle relaxation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Poplar lignin structural changes during extraction in γ-valerolactone (GVL)

In this paper, we describe an approach for producing both high quality and high quantity of lignin through studying the structural change of lignin during treatment of poplar wood in γ-valerolactone (GVL) for a range of temperatures (from 80 to 120 °C) and reaction time at temperature (from 1 to 24 h). Throughout the study, various techniques, including nuclear magnetic resonance (NMR) spectroscopies (solution- and gel-state 1 H –13 C 2D HSQC and 31 P) and gel-permeation chromatography (GPC) were applied to characterize the lignin structures. As the GVL-extracted lignin yield increases, the level of β-ether units decreases and the level of condensation products increases. The β-ether content, the aliphatic hydroxyl group content, and the molecular weight of the GVL-extracted lignin fractions were close to the poplar lignin from other preparation methods (e.g., enzyme lignin). A two-step hydrolytic process (120 °C, 2 × 15 min) gave a higher lignin yield (56.5% vs. 54.8%) with three times higher β-ether content (31.9% vs. 10.6%) than lignin extracted from a single-step process at 120 °C for 1 h. The results demonstrate that multiple-step cycling of cosolvent-assisted hydrolysis can help preserve more of the virgin ether-bond structures of GVL-extracted poplar lignin. Such a strategy can also be applied to a fully continuous-flow reactor system in future research to further improve both the productivity and quality of GVL-extracted lignin.

2D HSQC NMR↗