Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Small Reactor”

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 973 records · Page 54

Performance assessment of long-legged tightly-baffled divertor geometries in the ARC reactor concept

Extremely intense power exhaust channels are projected for tokamak-based fusion power reactors; a means to handle them remains to be demonstrated. Advanced divertor configurations have been proposed as potential solutions. Recent modelling of tightly baffled, long-legged divertor geometries for the divertor test tokamak concept, ADX, has shown that these concepts may access passively stable, fully detached regimes over a broad range of parameters. The question remains as to how such divertors may perform in a reactor setting. To explore this, numerical simulations are performed with UEDGE for the long-legged divertor geometry proposed for the ARC pilot plant conceptual design - a device with projected heat flux power width (λq||) of 0.4 mm and power exhaust of 93 MW - first for a simplified Super-X divertor configuration (SXD) and then for the actual X-point target divertor (XPTD) being proposed. It is found that the SXD, combined with 0.5% fixed-fraction neon impurity concentration, can produce passively stable, detached divertor regimes for power exhausts in the range of 80-108 MW - fully accommodating ARC's power exhaust. The XPTD configuration is found to reduce the strike-point temperature by a factor of ~10 compared to the SXD for small separations (~1.4λq||) between main and divertor X-point magnetic flux surfaces. Even greater potential reductions are identified for reducing separations to ~1λq|| or less. The power handling response is found to be insensitive to the level of cross-field convective or diffusive transport assumed in the divertor leg. By raising the separatrix density by a factor of 1.5, stable fully detached divertor solutions are obtained that fully accommodate the ARC exhaust power without impurity seeding. To our knowledge, this is the first time an impurity-free divertor power handling scenario has been obtained in edge modelling for a tokamak fusion power reactor with λq|| of 0.4 mm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Temperature-dependent model of helium bubble expansion and bursting in tungsten

Continuum-scale models that can reliably predict the behavior of helium in tungsten are of interest to the fusion community due to the projected impact of these materials on fusion reactor operation. In this work, we perform molecular dynamics simulations of spherical helium bubbles of various sizes in tungsten at different temperatures and depths with the goal of determining a mathematical model of the pressure and density at which the bubbles initially expand or burst as a function of depth, size, temperature, and surface orientation. The bubbles are small enough that their loop-punching pressures cannot be accurately predicted with continuum mechanics, and their expansion behavior is important, as it appears to cause many of the features observed on helium-irradiated tungsten surfaces. We vary the temperature, bubble size, bubble depth, and surface orientation in each case, recording the bubble pressure and density that result in bubble expansion. An exponential function with three adjustable parameters is found to fit the results well; the parameters that best fit our results are provided.

36 MATERIALS SCIENCE↗

Characterization of the ELM-free negative triangularity edge on DIII-D

Tokamak plasmas with strong negative triangularity (NT) shaping typically exhibit fundamentally different edge behavior than conventional L-mode or H-mode plasmas. On DIII-D, every plasma with sufficiently negative triangularity (δ < δ crit ≃ -0.12) is found to be inherently free of edge localized modes (ELMs), even at injected powers well above the predicted L-H power threshold. It is also possible to access an ELM free state at weaker average triangularities provided that at least one of the two x-points is still sufficiently negative. Access to the ELM-free NT scenario is found to coincide with the closure of the second stability region for infinite-n ballooning modes, suggesting that ballooning stability may play a role in limiting the accessible pressure gradient in NT plasmas. Despite this, NT plasmas are able to support small pedestals and are typically characterized by an enhancement of edge pressure gradients beyond those found in traditional L mode plasmas. Further, the pressure gradient inside of this small pedestal is unusually steep, allowing access to high core performance that is competitive with other ELM-free regimes previously achieved on DIII-D. Since ELM-free operation in NT is linked directly to the magnetic geometry, NT fusion pilot plants are predicted to maintain advantageous edge conditions even in burning plasma regimes, potentially eliminating reactor core-integration issues caused by ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reactor System Facility Modification to Detect Compromised Human Machine Interfaces

This study focuses on a multi-layered Industrial Control System (ICS)/Operational Technology (OT) security architecture to aid in the discovery and mitigation of compromised Human Machine Interface (HMI)/Instrumentation & Control (I&C) based systems for modifying a prototypical reactor condition test facility called the Flowing Autoclave System (FAS) at Idaho National Laboratory (INL). This is achieved through a three-layered combination of network security solutions, hash-based algorithms, and blockchain technologies. Hash algorithms are mathematical functions used to generate a predetermined set of fixed-length values. They are widely used in computer security to verify the integrity of system information and data, both on a local network and the wider internet. Even small amounts of unauthorized system modification will cause the hash algorithm to output a set of characters that deviate significantly from its original value. Assisting secure hash functions, blockchain technology is a secure and distributed technology used to provide an immutable set of records replicated on all devices within a decentralized network. Blockchain offers a cost-effective solution to detect system compromise by providing a traceable breadcrumb trail of all network activity and data modification happening on a system. If both are used in conjunction with network monitoring tools, the integration of this three-pronged approach can become an asset in detecting suspected system compromises before any real damage can occur.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

2019 accomplishments: degradation of mechanical properties in structural metals and welds for FTS reservoirs

This report documents work in 2019 at SRNL in support of the Aging and Lifetimes program. Specimens for fracture toughness testing were prepared in prior years by thermally pre-charging samples, provided by Sandia National Laboratories, of welded types 304L and 21-6-9 stainless steels. These samples were pre-charged in 2017. The fracture toughness of the first aged specimens was measured and compared to historical data. The measurements of the 304L and 21-6-9 weld and particularly HAZ represent some of the first fracture data relevant to modern GTS reservoirs. This work represents an on-going collaboration between SRNL and SNL to understand tritium embrittlement of structural metals in Gas Transfer System reservoirs which informs lifetime assessments. Continued testing was conducted on ongoing historical studies, including the effects of very long-term aging, and the relative performance of differently forged stockpile relevant 304L and 316L samples. New samples were pre-charged with tritium and placed in a freezer for an aging study to allow helium to build into the material. The precharged samples consist of smooth and notched tensile samples of 304L stainless steel, as well as tube tensile samples of different thermal treatments supplied by Sandia National Laboratories. Their charging marks the beginning of a new sequence of studies to explore fundamental deformation mechanisms underpinning tritium embrittlement. At predesignated times, samples will be removed from aging and tested to provide data which will help develop an understanding of material degradation with increasing aging time. Also documented is progress in employing various advanced characterization techniques including high resolution electron microscopy, thermally programed desorption, and small angle x-ray scattering, as well as facilities renovations and equipment procurements in support of sustained tritium impacts on materials program at SRNL.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Understanding the formation of a low-pressure pedestal in the presence of a strong internal transport barrier in DIII-D high β p plasmas

As a promising scenario for fusion reactors, the high poloidal-beta (β P ) scenario is characterized by a strong large radius internal transport barrier (ITB), which significantly enhances the overall confinement quality and the bootstrap current fraction for fully non-inductive operation. It is frequently observed that in the presence of a strong ITB, the pedestal height is lower and is accompanied by small edge localized modes (ELMs), which further improves the compatibility of a high performance core with an edge solution. A mechanism for the formation of the low pedestal is proposed in this paper. It is found that the strong ITB creates an off-axis bootstrap current to clamp the local safety factor q, and thus the magnetic shear in the outer core/pedestal region is increased. Gyrokinetic simulations with the CGYRO code show that the higher magnetic shear brings the experimental profiles into the range where the growth rate of drift-wave instabilities and thus transport is higher, and therefore a lower pedestal gradient is expected. Here, the combination of low pedestal and high magnetic shear further enhances the turbulent transport across the whole pedestal, consistent with power balance analysis. Such a positive feedback mechanism ultimately results in a lower pressure pedestal as observed in experiments. Under such a low pedestal, linear simulations with BOUT++ predict the growth rates of peeling–ballooning modes to be lower across the whole toroidal mode number spectra, and the nonlinear BOUT++ simulation exhibits lower saturated fluctuation intensity as well, consistent with the experimentally observed lower ELM size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of Accelerated High Temperature Mechanical Testing Techniques

The Advanced Materials and Manufacturing Technologies (AMMT) Program focuses on advancing materials and manufacturing techniques for nuclear energy applications, particularly in the qualification of materials for high-temperature structural use. This report presents work on refining the creep testing of small specimen geometries. Efforts include the development of a new specimen geometry for sub-sized specimens, which were subjected to uniaxial creep tests. The results contribute to the understanding of material behavior under stress at elevated temperatures and offer potential improvements in creep data collection methods. These findings support ongoing advancements in material qualification processes essential for nuclear reactor applications.

36 MATERIALS SCIENCE↗

Simulation Accelerator

Under a NASA SBIR (Small Business Innovative Research) contract, (NAS5-30905), EAI Simulation Associates, Inc., developed a new digital simulation computer, Starlight(tm). With an architecture based on the analog model of computation, Starlight(tm) outperforms all other computers on a wide range of continuous system simulation. This system is used in a variety of applications, including aerospace, automotive, electric power and chemical reactors.

Source record↗

Modeling Continuous Online Refueling with SCALE 6.3.1 and Serpent-2 in the EIRENE Novel Molten Salt Reactor Design

The Molten Salt Reactor (MSR) is a Generation IV advanced fission reactor design in which the coolant, and in some cases the fuel itself, is in the form of liquid molten alkali-halide salts with a fluoride or chloride ionic base. In liquid-fueled MSRs, reactor refueling may be performed online, where fresh fuel salt is added to the core during operation without the need for shutdown periods. Refueling for these reactor designs is typically modeled using either a batch or continuous refueling approach, both of which can be simulated with the SCALE 6.3.1 and Serpent-2 code systems. However, the specific manner in which they are implemented can vary depending upon the desired rate of refueling, whether the refueling rate is constant or variable, and consideration of salt drainage for systems where the in-core salt volume is kept constant. In a novel MSR fuel cycle concept termed the “Sourdough” fuel cycle, fuel salt is allowed to “grow” within the core, with excess salt either being transferred to an external holding tank to maintain a constant core volume or diverted to an upper plenum within the core to allow for volume growth. In this work, continuous refueling was modeled in a thermal-spectrum, LEU-fueled, small MSR design operating with the Sourdough fuel cycle using the SCALE 6.3.1 and Serpent-2 codes. Two different continuous refueling approaches were simulated, with the performance of each being compared to determine which is most suitable for use with the Sourdough fuel cycle concept.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Implementation of an Experimental Setup Based on the Three-Omega Method for Thermal Conductivity Measurements of Molten Salts

Molten salt reactors (MSRs) are one of the most promising next-generation nuclear reactor designs which use molten salts as the fuel and coolant. To design, model, predict, license, and operate MSRs, the behavior of the fuel and coolant must be understood as a function of temperature and irradiation, i.e., burnup. Of particular importance and complexity is the reliable and reproducible measurement of thermal conductivity of molten salts at high temperatures (> 500 °C). Thermal conductivity provides information on the rate at which heat passes through a specified material and is dependent on chemical composition. Currently, in the national laboratory complex, several methods for measuring thermal conductivity on molten salt systems are under development, however, to date, consistent data have not been produced. Typically, thermal conductivity measurements of high-temperature melts have been conducted by using the Transient Hot Wire (THW) method where thermal conductivity is determined by the change in temperature over time at a certain distance from a heating wire. Another time-dependent technique is the laser flash method, commonly used for solid samples. Laser flash requires prior measurements of density and heat capacity to calculate the thermal conductivity of the sample, introducing additional measurement errors. Steady state and time-dependent measurements such as THW and laser flash have errors from convection and ambient temperature fluctuations (heat radiation). In order to produce high fidelity thermal conductivity data, a measurement technique needs to be developed that will reduce the effect of convection at high temperatures by measuring in a very fast time window. In addition, to measure on actinide and irradiated salts the sample size should be small. Therefore, this project aims to develop a frequency-domain technique based on the three-omega (3?) method, previously proven to work on solids, nanofluids and low temperature solar salts (25 to 250 °C). If the 3?-device can be designed with materials that can withstand high temperatures (>500 °C) and can be configured to reduce errors due to convection, then it would prove to be a valuable technique for ensuring thermal conductivity on molten salt systems needed for MSR development. For this proposal, the 3? instrument will be designed, assembled, and tested on non-radioactive salts before being used on actinide bearing salts. In addition, a thermal model will be developed to determine the thermal conductivity of the sample from the measured thermal oscillation of the sensor. Because 3? is a direct method that requires a small volume of sample and works in a very fast time window, it will aid in suppressing the interference due to convection which is ideal when working with high temperature liquid samples.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hygroscopic growth of UO 2 F 2 nanoparticles

Hygroscopicity is an important physicochemical property of aerosol that describes the ability of a particle to uptake water. The hygroscopic properties of uranyl fluoride (UO 2 F 2 ) aerosol generated from a UF 6 hydrolysis reactor was investigated for the first time using a custom-built Humidified Tandem Differential Mobility Analyzer (HTDMA). The HTDMA is capable of measuring UO 2 F 2 nanoparticle growth determined by mobility size over a wide range of atmospheric humidity from dry conditions at <10% relative humidity (RH) to 85% RH. The hygroscopic properties were determined for nanoparticles as small as 3.5 nm in this study. Although the largest size of UO 2 F 2 nanoparticles was 80 nm, monodisperse aerosol with a mobility diameter of up to approximately 500 nm can be investigated using the HTDMA. Anhydrous UO 2 F 2 nanoparticles with a mobility diameter of 3.5 nm were shown to be highly hygroscopic with a deliquescence relative humidity (DRH) of 10%. Hydrates with a larger mobility diameter from 10 to 80 nm were non-hygroscopic with no observable DRH and limited water uptake up to 85% RH. Here, these results demonstrate the hygroscopic properties of UO 2 F 2 nanoparticles are highly variable and based on both the mobility size and hydration state. Hygroscopicity affects the physicochemical properties of UO 2 F 2 nanoparticles, including the aerosol phase state and viscosity, with impacts on aerosol growth, coagulation, and deposition that is critical for understanding the fate and transport of UO 2 F 2 nanoparticles in the atmosphere.

Hygroscopicity↗

Characteristics of oxide-dispersion strengthened alloys produced by high-temperature severe deformation

This study is to explore an economically attractive and technically feasible processing method for oxide-nanoparticle strengthened alloys for fusion reactor application. Despite many scientific merits of the advanced oxide-dispersion strengthened (ODS) alloys, such as the nanostructured ferritic alloy (NFA) 14YWT, the only viable production path for a high-quality NFA is the high-power mechanical alloying process. This process is often a multi-day high-speed ball milling of alloy powder with a small quantity of yttria (Y 2 O 3 ) powder, followed by the milled-powder consolidation using extrusion or other methods and additional thermomechanical processing (TMP) for property control. This complex production path, including the low-temperature mechanical alloying in particular, has limited technical advancement toward the cost-effective and scale-up production of ODS alloy components. To overcome such a practical limitation, we proposed to explore alternative low-cost processing routes using traditional thermomechanical processing (TMP) method only. Further, a series of continuous TMP cycles, which were designed to impose high-temperature severe plastic deformation (HT-SPD) conditions to the consolidated powder mixtures, were applied to achieve the effective distribution of oxide particles in nanograin structure and thus desirable mechanical properties. Since the reduced-activation ferritic-martensitic (RAFM) alloy powders (Fe-10Cr and Fe-14Cr alloys with various Y contents) are available in our inventory, we focused to utilize the new solid-state synthesis approach for controlling oxide (oxygen source) dissolution and nanoscale clustering in nanograin structure in those alloys. A combination of powder consolidation at 900 °C and continuous thermomechanical activation at 600 °C yielded two essential ODS alloy microstructure contents–nanograin structure and nanoparticle distribution–and thus demonstrated a good combination of strength and ductility.

36 MATERIALS SCIENCE↗

Hygroscopic Growth of UO 2 F 2 Particles

Hygroscopicity is an important physicochemical property of aerosol that describes the ability of a particle to uptake water. The hygroscopic properties of uranyl fluoride (UO 2 F 2 ) aerosol generated from a UF 6 hydrolysis reactor was investigated for the first time using a custom-built Humidified Tandem Differential Mobility Analyzer (HTDMA). The HTDMA is capable of measuring UO 2 F 2 particle growth determined by mobility size over a wide range of atmospheric humidity from dry conditions at <10% relative humidity (RH) to 85% RH. The hygroscopic properties were determined for nanoparticles as small as 3.5 nm in this study. Anhydrous UO 2 F 2 particles with a mobility diameter of 3.5 nm were shown to be highly hygroscopic with a deliquescence relative humidity (DRH) of 10%. Hydrates with a larger mobility diameter of 10 to 80 nm were non-hygroscopic with no observable DRH and limited water uptake up to 85% RH. These results demonstrate the hygroscopic properties of UO 2 F 2 particles are highly variable and based on both the mobility size and hydration state. Hygroscopicity affects the physicochemical properties of UO 2 F 2 particles, including the aerosol phase state and viscosity, with impacts on aerosol growth, coagulation, and deposition that is critical for understanding the fate and transport of UO 2 F 2 particles in the atmosphere.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Plasma assisted distributed chemical production

The present development is a process to produce commodity chemicals such as methanol and syngas using an integrated plasma catalysis technology. The method comprises providing a fixed or fluidized bed reactor having a microwave plasma flame and a catalyst bed with a catalyst, wherein the catalyst is an alloyed bimetallic nanowire. In the process, the plasma flame fluidizes the catalyst thereby producing a more effective catalyst than the non-fluidized catalyst. It is anticipated that the reactor can have a throughput capacity of up to 30 Lpm/kW and can be effective for the conversion of CO 2 , CH 4 , air, water, and combinations thereof, through reactions such as pure CO 2 splitting, reverse water gas shift (RWGS) for CO production, methanol synthesis, and plasma reforming of methane, thereby making a system that would be attractive for small GTL units.

Sunkara, Mahendra↗

Fluid structural interaction simulation of the MP-1 plate performance irradiated in the advanced test reactor

Due to the high neutron fluxes they generate and increased heat transfer performance, plate type fuels are used in the U.S. high-performance research reactors. During irradiation, a significant amount of fission energy (i.e., ~200 MeV per atom) is released by the U-235 chain reaction that is carried away by the coolant. Through thermal hydraulic analysis, the coolant’s heat transfer capability is investigated to ensure that the desired fuel temperature can be maintained. At high temperatures, the plate undergoes elastic/plastic deformation, creep, and swelling as a result of both the temperature gradients and the fission gas production within the fuel. These effects are studied via fuel performance analysis. When the plate deformation is small enough that the coolant flow remains relatively unchanged, conducting these two types of analyses independently will suffice. But at high fuel burnups, the swollen plates may encroach into the coolant channels that separate the fuel plates from each other and cause these channels to narrow. Large reductions in channel gap size imperil cooling performance, causing fuel temperatures to rise. Further, if the plate deformation is asymmetric, the fuel centerline will shift toward one side of the channel, causing an uneven reduction in coolability. In addition, a boehmite (oxide) layer will, over time, grow on the plate surface, further obstructing the heat transfer process. To precisely predict fuel temperatures/deformation, a complete coupled analysis that considers coolant flow, heat transfer, oxide growth, elastic/plastic deformation, creep, and swelling is needed; however, this type of analysis method is not available in the literature. To fill the gap, this research developed a fluid structure interaction (FSI) approach to the fuel plate analysis, then successfully applied it to the Mini-Plate (MP)-1 experiment, which was irradiated in the Advanced Test Reactor (ATR) for both one and two cycles. The complete analysis coupled STAR-CCM+, a computational fluid dynamic (CFD) software for calculating flow, with Abaqus, a finite element analysis code for calculating plate deformation. Improvements in the results were found when comparing the fully coupled analysis to the independently conducted analyses but they were not significant due to the miniature size of the plates and the relatively short irradiation time. In the future, the fully coupled approach presented herein will be applied to full-size fuel plates with longer irradiation cycles once additional experiments become available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Representative-cell-based Modeling of HTTF in SAM

Although the current nuclear power market is primarily occupied by the light water reactors (LWRs), the TRISO-fueled helium-cooled graphite-moderated high temperature gas reactors (HTGRs) are drawing growing attention as the nuclear power industry marches towards more advanced systems. Compared to that of conventional LWRs, the thermal hydraulics of HTGR cores show extra complexity from their multi-scale heat transfer mechanisms with varying importance during different operation modes or transient stages. The modeling of HTGR cores for the system analysis purpose therefore faces the challenge of reaching reasonable fidelity and accuracy while maintaining sufficient simplicity. Previous development and validation efforts have demonstrated that a “2-D ring model” with reasonable performance can be implemented in the System Analysis Module (SAM) for prismatic HTGR cores. In the current project, an alternative modeling methodology based on “representative cells” is proposed for the typical prismatic core of HTGRs. Different from the previous ring model approach, the proposed methodology first separately models and then combines the small- and large-scale thermal conductions, by connecting representative local heat transfer units (cells) with effective core- wise thermal resistance. A modeling practice for an integral high-temperature test facility (HTTF) elaborates the modeling details. Steady-state validation of the resultant candidate model is performed against a higher-resolution benchmark from a 3D-1D coupled simulation, which shows satisfactory prediction with reasonably captured global parameters and well-represented temperature fields. A postulated pressurized conduction cooldown (PCC) is also simulated and analyzed, demonstrating the model’s capability of transient prediction with physically captured phenomena resolved in both small and large temporal and spatial scales. In general, the current work achieves a preliminary success in proposing a methodology using representative cells to model prismatic HTGR cores in SAM. Future efforts are envisioned with more validation activities and with potential modeling improvements to eventually achieve the high confidence on a high-fidelity robust modeling methodology with reasonable accuracy.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗

Energy-Efficient Microbial Electrochemical Lignin and Alkaline Hydroxide Recovery from DMR Black Liquor

Black liquor is a complex waste stream generated from the fast-expanding biorefinery industry, and the high management cost has been known as a major development barrier. In this work, we demonstrated a rationally designed 4-chamber microbial electrolysis process that was capable of recovering the value-added lignin, NaOH, and H2 products while removing waste organics from the Deacetylation and Mechanical Refining (DMR) black liquor, which could potentially reduce the disposal cost by up to 40%. The lab-scale MEC achieved 52% organic removal from a highly concentrated feedstock (28,960 ± 212 mg/L COD) and recovered 83.2 mM NaOH under a small voltage application (1.2 V). Lignin recovery was observed to be as high as 13.98 +/- 0.28 g/L, and 0.35 L/L/day H2 was produced with a faradaic efficiency of 93.5%. This work provided a good example on how complex liquid waste stream can be valorized via multi-functional microbial electrochemical technology, and how the traditionally unfavored pH gradient in an electrochemical reactor can be manipulated for lignin precipitation and alkaline recovery at the same time. The versatility and good performance in concurrent contaminants removal, lignin precipitation, H2 production and NaOH recovery demonstrate a good potential for microbial electrolysis in black liquor treatment.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗