Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tailoring”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Analysis of Laminar Convective Heat Transfer Over Structured Non-Wetting Surfaces

Structured non-wetting surfaces provide alternating no-slip and partial slip boundary conditions to the fluid flow which, in turn, affects the convective heat transfer performance over the surfaces. In this work, an analytical model is developed for the interfacial Nusselt number, the overall Nusselt number and a thermal hydraulic performance factor for fluid flow in a cylinder patterned with structured non-wetting surfaces, for the two cases of uniform wall heat flux and uniform wall temperature. In addition, by considering the stability of the Cassie state of wettability and its transition to the Wenzel state for flow over superhydrophobic surfaces, the present model overcomes certain limitations of the previously reported studies in the literature. Based on the analytical formulations and the stability constraints, the present paper provides optimum design maps for tailoring structured non-wetting surfaces for maximizing convective heat transfer and the combined thermal-hydraulic performance in applications. Use of the design maps on example cases is also discussed. It is shown that the use of structured non-wetting surfaces is most effective for low Reynolds numbers and/or small cylinder radius.

42 ENGINEERING↗
DOE OSTI2021

Unified descriptor for enhanced critical heat flux during pool boiling of hemi-wicking surfaces

Boiling heat transfer is dictated by interfacial phenomena at the three-phase contact line where vapor bubbles form on the surface. Structured surfaces have shown significant enhancement in critical heat flux (CHF) during pool boiling by tailoring interfacial phenomena. This CHF enhancement has been primarily explained by two structural effects: roughness, which extends the contact line length at the bubble base, and wickability, the ability to imbibe liquid through surface structures by capillary pumping. Here, in this work, we show that CHF enhancement on structured surfaces cannot be described by roughness or wickability alone. This result was confirmed using systematically designed micropillar surfaces with controlled roughness and wickability. Further, we performed a scaling analysis and derived a unified descriptor, which represents the combined effects of thin film density and volumetric wicking rate. This unified descriptor shows a reasonable correlation with CHF values with our experiments and literature data. This work provides important insights in understanding the role of surface structures on CHF enhancement, thereby providing guidelines for the systematic design of surface structures for enhanced pool boiling heat transfer.

42 ENGINEERING↗
DOE OSTI2024

Thermal conductivity of 3D-printed block-copolymer-inspired structures

This study primarily focuses on examining the impact that geometric structure has on thermal conductivity of multi-phase constructs in different 3D-printed poly(lactic acid), PLA, samples. The investigated structures are inspired by morphologies formed by diblock copolymers: lamellae, hexagonally packed cylinders, and gyroid. This research also investigates how volume percentage and material combination influence the thermal conductivity of these structures. Further, the samples can be tailored to simulate various thermal management structures observed in practical applications, such as thermal interface materials in electronic devices. Thermal conductivity ratio is controlled using air, the least conductive material at 0.026 W/(m K), PLA at 0.136 W/(m K), and thermal paste at 5.11 W/(m K). Different models were tested against thermal conductivity measurements in order to capture the effect of material type (PLA-Air versus PLA-Thermal Paste), volume percentage, structure, and orientation. Simple, effective medium models were good predictions of thermal conductivity in lamellar structures, but it was necessary to develop models for conduction through cylindrical and gyroid structures. Finally, all results were normalized to find a universal model that is independent of structure and material. This approach provides a simple method to predict how to reduce or enhance transport properties and heat management capabilities of 3D printed objects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Numerical simulation of compressible fluid-dynamics in the chamber of inertial fusion energy systems

Here, this paper aims to establish new and innovative modeling capabilities for analyzing chambers in Inertial Fusion Energy (IFE) systems. IFE is emerging as a promising method to achieve fusion power production, but several challenges must be overcome to develop an IFE pilot plant or deploy commercial IFE systems. These challenges are both theoretical and technical, encompassing a deeper understanding of the underlying physical phenomena and the development of new technologies and materials. One of the needs is to develop mathematical models to describe IFE systems and numerical tools to simulate them. This paper contributes to this endeavor by presenting a new OpenFOAM solver for IFE systems, focusing on gas dynamics in their chambers. The analysis and development of chamber designs will play a significant role in the transition from single-shot experiments to high-repetition rates, as there is a need to protect the chamber walls from the intense radiation fields produced by fusion reactions. A promising design option, normally referred to as thick wall chamber design, consists in using lithium or molten salt jet arrays within the chamber. A critical phenomenon is the venting of high-pressure gases from the center to the external part of the chamber, passing through the blanket array. This process involves the propagation and attenuation of strong pressure waves, requiring suitable modeling approaches for compressible fluid-dynamics. The solver proposed in this work implements a multi-material hydrodynamics model tailored to accurately describe the non-linear propagation of pressure waves while avoiding numerical oscillation issues typical of high-velocity compressible simulation. This solver is verified against numerical test cases, validated against experimental data, and applied to the analysis of the High-Yield Lithium-Injection Fusion-Energy (HYLIFE-I) concept. The relevance of this paper is threefold. Firstly, it contributes to developing and testing modeling approaches for compressible fluid-dynamics phenomena, with specific focus on the new and unexplored topic of IFE thick-liquid-wall blanket modeling. Secondly, it marks one of the first applications of the OpenFOAM library in the research field of IFE systems. Finally, the investigated problem is of practical interest for IFE developers, as it provides useful indications about relevant phenomena in pressure wave propagation in the chamber of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experimental characterization of a triply periodic minimal surface PCM-to-air thermal energy storage device

Here, the expansion of energy production has intensified research into thermal energy storage (TES) to manage variability and improve system efficiency. Heat exchanger design is critical to achieving high power densities in TES systems. This study presents a high-surface-area phase change material (PCM)-to-air heat exchanger fabricated via resin-based stereolithography with a novel gyroid-based geometry tailored for enhanced performance. Material properties of the commercial PCM and resin were characterized using analytical techniques. A controlled air loop was used to evaluate heat transfer and pressure drop at various flow rates and inlet temperatures. Results show a strong dependence on the inlet temperature difference (ΔT) relative to the PCM melting point. During charging at the highest flow rate, increasing the inlet air temperature from ΔT = 5 °C to 20 °C above the melting point increased the average heat transfer rate by 187%. During discharging, decreasing the inlet temperature by the same ΔT below the melting point led to a 232% increase. Notably, the high-surface-area design enabled nearly symmetric charging and discharging behavior, a novel result for PCM-based TES systems which are often restricted by natural convection and other effects. The overall heat transfer coefficient was calculated and compared to values from standard design correlations. The maximum thermal effectiveness reached 95% at ΔT = 20 °C and a moderate flow rate of 34 m 3 /h. Peak coefficients of performance (COP) of 7 during charging and 6.4 during discharging were observed at ΔT = 20 °C and a low flow rate of 20 m 3 /h. These results demonstrate the viability of additively manufactured geometries for advanced TES applications.

25 ENERGY STORAGE↗
DOE OSTI2026

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Modeling The Effects of Loading Scenario and Thermal Expansion Coefficient on Potential Failure of Cryo-compressed Hydrogen Vessels

A multiscale thermomechanical model for a simplified Type-3 cryogenic compressed hydrogen (H2) storage vessel is developed in this paper. The model accounts for the temperature-dependent elastic-plastic behavior of the vessel carbon/epoxy composite overwrap and aluminum alloy liner. The homogenized thermo-elastic-plastic behavior for the individual laminas of the vessel layup is obtained by an incremental Eshelby-Mori-Tanka approach associated with a micromechanical failure criterion to predict lamina failure while a standard elastic-plastic constitutive model is used to describe the behavior of a typical aluminum alloy assumed for the liner. The vessel response to external loadings is achieved by a finite element method. Four loading scenarios representing four thermomechanical cycles applied to the vessel are analyzed to evaluate constituent and lamina stresses as well as the associate failure criterion during the cycle according to these scenarios. The model can provide helpful guidance to mitigate thermal stresses by an adequate selection of loading scenario, optimizing the layup and by tailoring thermomechanical properties of the resin matrix.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermodynamic parametric analysis of refueling heavy-duty hydrogen fuel-cell electric vehicles

Reliable design and safe operation of heavy-duty hydrogen refueling stations are essential for the successful deployment of heavy-duty fuel cell electric vehicles (FCEVs). Fueling heavy-duty FCEVs is different from light-duty vehicles in terms of the dispensed hydrogen quantities and fueling rates, requiring tailored fueling station design for each vehicle class. In particular, the selection and design of the onboard hydrogen storage tank system and the fueling performance requirements influence the safe design of hydrogen fueling stations. A thermodynamic modeling and analysis are performed to evaluate the impact of various fueling parameters and boundary conditions on the fueling performance of heavy-duty FCEVs. Here, we studied the effect of dispenser pressure ramp rate and precooling temperature, initial tank temperature and pressure, ambient temperature, and onboard storage design parameters, such as onboard storage pipe diameter and length, on the fueling rate and final vehicle state-of-charge, while observing prescribed tank pressure and temperature safety limits. An important finding was the sensitivity of the temporal fueling rate profile and the final tank state of charge to the design factors impacting pressure drop between the dispenser and vehicle tank, including onboard storage pipe diameter selection, and flow coefficients of nozzle, valves, and fittings. The fueling rate profile impacts the design and cost of the hydrogen precooling unit upstream of the dispenser.

08 HYDROGEN↗
DOE OSTI2024

Surface exsolving perovskite ceramics as catalyst for microwave methane pyrolysis to co-generate hydrogen and carbon nanotube

In this study, we have successfully developed a novel surface exsolving perovskite ceramic, which demonstrates the ability to simultaneously generate CO x -free hydrogen and carbon nanotubes (CNTs) through methane pyrolysis under the influence of microwave irradiation. We conducted an extensive survey and optimization of various perovskite materials specifically tailored for microwave applications. Among the materials investigated, nickel-doped strontium titanium oxide (STON) emerged as the most promising candidate, exhibiting both a satisfactory methane conversion rate and excellent responsiveness to microwave irradiation. Refinement of STON was carried out by fine-tuning the Ni content, optimizing the reduction dwell time, and adjusting the reduction temperature. Notably, SrTiNi 0.08 O 3 (STON8), demonstrated an impressive initial methane conversion rate of up to 40%. Transmission Electron Microscopy (TEM) provided visual evidence of the correlation between Ni content and reduction temperature, with respect to the exsolved Ni metal particle size. This finding highlights the immense potential of surface exsolving perovskite ceramics as a highly effective catalyst for the simultaneous production of CNTs and COx-free hydrogen via methane pyrolysis under microwave irradiation. Furthermore, it represents a step forward in the field of catalytic materials and microwave-driven processes.

08 HYDROGEN↗
DOE OSTI2025

Rapid synthesis of phase-engineered tungsten carbide electrocatalysts via flash joule heating for high-current-density hydrogen evolution

Fabricating durable and high-performance electrocatalysts operating at high current densities for industrial acidic hydrogen evolution remains a daunting challenge. Tailoring the phase composition of electrocatalysts is a promising strategy to harness synergistic effects and improve charge transfer, thereby optimizing their performance. Here, this work presents a fast, green method based on flash joule heating (FJH) to synthesize phase-engineered tungsten carbide electrocatalysts for the acidic hydrogen evolution reaction (HER) at high current densities. Tungsten carbide electrodes with varying FJH treatment durations (3, 10, 30, and 60 s) are fabricated to fine-tune the mixture of tungsten monocarbide (WC) and tungsten semicarbide (W 2 C) phases. Results show that samples with a 30-s treatment (TC-3) exhibit an optimal balance between these phases, leading to a low overpotential of 180.97 and 387 mV at current densities of 10 mA/cm 2 and 4 A/cm 2 , respectively. TC-3 exhibits significantly lower charge transfer resistance compared to the other electrocatalysts, which can be attributed to its optimal phase ratio. Notably, the TC-3 electrocatalyst remains stable for over 9 days at 4 A/cm 2 due to their controlled phases and excellent corrosion-resistant properties. This work highlights a new method to fabricate cost-effective, high-performance tungsten carbide electrocatalysts with well-controlled phase compositions.

Acidic hydrogen evolution reaction↗
DOE OSTI2026

Toward the application of the Risk-Based Maintenance approach to a hydrogen-fueled manufacturing plant

Safety, reliability, and maintenance of hydrogen-based equipment, including Risk-Based Maintenance (RBM), have recently gained increasing attention, as hydrogen, while key to decarbonizing hard-to-abate sectors, raises safety concerns. However, two literature gaps limit the development of risk-based strategies for hydrogen technologies. First, existing RBM methodologies do not account for hydrogen-specific components, such as electrolyzers, which differ from conventional equipment. Second, available studies on electrolyzer reliability remain largely qualitative or laboratory-scale. This study addresses these gaps by proposing an adapted RBM methodology tailored to hydrogen-fueled manufacturing facilities. Integrating qualitative tools (e.g., FMEA) with probabilistic models (e.g., Bayesian Networks) enables comprehensive hazardous scenarios identification and likelihood estimation. A conceptual layout of a glass furnace supplied by a 3 MW PEM electrolyzer is considered as a case study to demonstrate the feasibility of the adapted RBM framework in identifying high-risk components and enabling maintenance prioritization, potentially improving plant safety.

08 HYDROGEN↗
DOE OSTI2026

Facile synthesis of Co(OH)2 nanoneedle arrays grown on stainless steel for industrial electrochemical oxygen evolution reaction

The electrochemical oxygen evolution reaction (OER) is a critical half-reaction in a variety of energy conversion and storage applications, however, OER suffers from sluggish kinetics due to the four proton-electron transfer processes. To remedy this, interfacial engineering proves an effective strategy to design high active OER catalysts. Herein, we report a facile synthesis of Co(OH)2 nanoneedle (Co-NN) arrays grown on stainless steel (SS) via a one-step hydrothermal reaction, and the resultant impressive OER performance. Particularly, the Co-NN/SS needs 267, 307, and 576 mV overpotentials to reach 10, 100, and 1000 mA cm−2 in 1 M KOH & room temperature, and the overpotentials drop to 199, 257, and 283 mV to achieve the same current densities in 30 wt% KOH & 80 °C. Additionally, an alkaline water electrolysis (AWE) cell coupled with Co-NN/SS anode and PtRu/NiMo cathode only requires 1.65 and 1.73 V iR-free cell voltage to reach 1.0 and 2.0 A cm−2, respectively. The sterling OER performance could be attributed to four possible reasons, the Fe in SS surface tailoring the electronic properties of Co(OH)2, the mixed metal oxides on SS surface accelerating surface reconstruction of Co(OH)2 nanoneedles into CoOOH active species, a Ni-rich surface layer formation cooperatively boosting the OER activity, and a local electric field effect concentrating reactants on the tip surface and accelerating the mass transfer. This work provides a facile approach for synthesizing highly efficient OER catalysts for industrial applications by interfacial engineering.

Lyu, Xiang [ORNL] (ORCID:0000000208673248)↗
DOE OSTI2023

An operando synchrotron study on the effect of wire melting state on solidification microstructures of Inconel 718 in wire-laser directed energy deposition

Directed energy deposition (DED) with a coaxial wire-laser configuration has gained significant attention in recent years for the production of large-scale metallic components because of its low directional dependence, fast deposition rate, high feedstock efficiency, and low manufacturing costs. This work studies the coaxial wire-laser DED process of Inconel 718 alloy under a stable deposition condition with a relatively low input volumetric energy density (55.5J/mm 3 ). Post characterization reveals a cluster of refined grains at the center-bottom region of the as-printed track. Operando high-energy synchrotron X-ray experiments and multi-physics modeling are applied innovatively to study the fundamental mechanism responsible for the formation of this microstructure. The X-ray diffraction experiment provides direct evidence, which is supported by the simulation, that the feeding wire can reach the melt pool bottom and release solid particles (primarily carbides) near the mushy zone owing to insufficient melting. Consequently, these sub-micron sized particles suppress the growth of large columnar grains and cause the formation of unique microstructural heterogeneity. In conclusion, this discovery offers new opportunities for tailoring the solidification microstructure by controlling the melting state of the feedstock wire in DED process, in addition to commonly known factors such as the thermal gradient and solidification velocity.

42 ENGINEERING↗
DOE OSTI2026

Influence of rolling reduction and annealing on recrystallization and grain structure in Ta-2.5W alloys

The microstructural evolution of wrought Tantalum - 2.5 wt% Tungsten (Ta-2.5W) alloys during thermomechanical processing is critical for optimizing their mechanical reliability in demanding applications such as aerospace, chemical processing, and nuclear technology. Despite the widespread use of Ta-W alloys, a comprehensive understanding of how rolling reduction, annealing temperature, and elemental inhomogeneity interact to determine recrystallization behavior and grain refinement remains incomplete. Here, in this study, we systematically investigate the effects of cold rolling and subsequent annealing on the microstructure of Ta-2.5W, with particular attention to grain orientation, stored energy, and elemental banding. Our results demonstrate that higher rolling reduction rates lower the onset and completion temperatures for recrystallization, resulting in finer and more homogeneous grain structures. Electron backscatter diffraction (EBSD) analysis reveals that grains with a 〈111〉 parallel to the plate normal possess higher stored energy and nucleate recrystallization more readily than grains with a 〈001〉 parallel to the plate normal. Elemental mapping shows that tungsten inhomogeneity leads to localized bands of accelerated recrystallization and hardness variation. These findings provide new insights into the mechanisms of microstructure refinement in Ta-2.5W alloys, offering guidance for tailoring processing routes to achieve superior performance in demanding engineering environments.

Annealing↗
DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗
DOE OSTI2020

Lignin extraction and upgrading using deep eutectic solvents

Deep eutectic solvents (DESs) are a class of green solvents exhibiting low vapor pressure, high thermal stability,low toxicity, and biodegradability. These features make them attractive solvents for lignocellulosic biomassprocessing, especially for lignin valorization. DESs not only solubilize lignin from lignocellulose biomass, butalso confer extracted lignin with properties favorable to its further upgrading. This review summarizes recentadvances in biomass fractionation using DESs for lignin extraction. It also covers lignin upgrading in DESs interms of depolymerization and functionalization. The challenges of tailoring lignin properties by DESs for differentvalorization routes are also discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

From acetone fractionation to lignin-based phenolic and polyurethane resins

In this paper, acetone fractionation was used to increase the homogeneity of softwood kraft lignin for product development. The acetone soluble and insoluble fractions were analyzed in terms of elemental composition, chemical structure, and thermal properties. The phenolic hydroxyl content of the acetone soluble fraction increased by 24%. The acetone soluble fraction of kraft lignin was used to fully substitute phenol in phenolic resin formulation. The resulting formulated lignin-formaldehyde adhesive showed higher adhesion strength (4.3 MPa) compared to a lab formulated phenol-formaldehyde resin (3.4 MPa) with 100% wood failure after lap shear test. On the other hand, the acetone insoluble lignin fraction, which had 23% higher aliphatic hydroxyl content compared to the starting lignin, was used to entirely replace petroleum-based polyol in polyurethane resin formulation. Then a bio-derived solvent (dihydrolevoglucosenone, Cyrene) was used to dissolve the acetone insoluble lignin fraction to synthesize lignin-based polyurethane resin. Tailoring the structure of kraft lignin through solvent fractionation resulted in two fractions, which were used to replace 100% of two major petrochemicals (phenol and polyol) in the formulations of phenolic and PU resins, that would otherwise be challenging due to the heterogeneity of softwood kraft lignin.

36 MATERIALS SCIENCE↗