SEARCH · Engineering Papers

Results for “Second-order”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2026

One-dimensional neutron diffraction from layered graphite: Reciprocal space structure and grating behavior

In this work we report observations of one-dimensional neutron diffraction from highly oriented pyrolytic graphite (HOPG), where the scattering angle varies continuously with incident angle following classical grating-like behavior. The 2D polycrystalline structure of HOPG—with highly aligned layers along the 𝑐 axis but random in-plane rotations—creates planes of scattering intensity in reciprocal space at 𝑄 𝑐 =𝑛⁢(2⁢𝜋/𝑑) where 𝑑=3.35Å is the interlayer spacing. As the Ewald sphere sweeps through reciprocal space during sample rotation, it continuously intersects these planes, producing the observed angular dispersion. We observe both first-order (𝑛=1) and second-order (𝑛=2) diffraction at conventional scattering angles (25°–70°), with peak positions that remain temperature-independent between 10 K and 294 K and follow quantitative agreement with momentum conservation 𝑄 𝑐 =𝑘⁢[sin⁡𝜓−sin⁡𝜓 𝑓 ]=𝑛⁢(2⁢𝜋/𝑑). X-ray diffraction under similar conditions shows no comparable behavior, confirming that sharp nuclear-vacuum contrast is essential. While diffraction intensities are weak (∼10 −6 of Bragg peaks), the observations demonstrate how the interplay of atomic-scale periodicity, nuclear contrast, and structural disorder enables observation of continuous diffraction curves at thermal neutron wavelengths, illustrating how HOPG's unique microstructure determines its scattering properties beyond conventional Bragg diffraction.

2-dimensional systems↗

DOE OSTI2025

Interacting Dirac magnons in the van der Waals ferromagnet CrBr 3

We study the effects of magnon-magnon interactions in the two-dimensional van der Waals ferromagnet CrBr 3 focusing on its honeycomb lattice structure. Motivated by earlier theoretical predictions of temperature-induced spectral shifts and van Hove singularities in the magnon dispersion [S. S. Pershoguba et al., Phys. Rev. X 8, 011010 (2018)], we go beyond the commonly used thermal magnon approximation by applying second-order perturbation theory in a fully numerical framework. Our analysis uncovers significant deviations from previous analysis: in particular, the predicted singularities are absent, consistent with recent inelastic neutron scattering measurements [S. E. Nikitin et al ., Phys. Rev. Lett. 129, 127201 (2022)]. Moreover, we find that the temperature dependence of the renormalized magnon spectrum exhibits a distinct 𝑇 3 behavior for the optical magnon branch, while retaining 𝑇 2 behavior for the acoustic or down magnon band. This feature sheds light on the collective dynamics of Dirac magnons and their interactions. We further compare the honeycomb case with a triangular Bravais lattice, relevant for ferromagnetic monolayer MnBi 2 ⁢Te 4 , and show that both systems lack singular features while displaying quite distinct thermal trends.

Holstein-Primakoff method↗

DOE OSTI2025

Dissipative Phase Transition in the Two-Photon Dicke Model

We explore the dissipative phase transition of the two-photon Dicke model, a topic that has garnered significant attention recently. Our analysis reveals that while single-photon loss does not stabilize the intrinsic instability in the model, the inclusion of two-photon loss restores stability, leading to the emergence of superradiant states, which coexist with the normal vacuum states. Using a second-order cumulant expansion for the photons, we derive an analytical description of the system in the thermodynamic limit, which agrees well with the exact calculation results. Additionally, we present the Wigner function for the system, shedding light on the breaking of the 𝑍4 symmetry inherent in the model. These findings offer valuable insights into stabilization mechanisms in open quantum systems and pave the way for exploring complex nonlinear dynamics in two-photon Dicke models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Tuning of altermagnetism by strain

For all collinear altermagnets, we sort out piezomagnetic free-energy invariants allowed in the nonrelativistic limit and relativistic piezomagnetic invariants bilinear in the Néel vector $\mathbf{L}$ and magnetization $\mathbf{M}$, which include strain-induced Dzyaloshinskii-Moriya interaction. The symmetry-allowed responses are fully determined by the nonrelativistic spin Laue group. In the nonrelativistic limit, two distinct mechanisms are discussed: the band-filling mechanism, which exists in metals and is illustrated using the simple two-dimensional Lieb lattice model, and the temperature-dependent exchange-driven mechanism, which is illustrated using first-principles calculations for transition-metal fluorides. The leading second-order nonrelativistic term in the strain-induced magnetization is also obtained for CrSb. Piezomagnetism due to the strain-induced Dzyaloshinskii-Moriya interaction is calculated from first principles for transition-metal fluorides, MnTe, and CrSb. Finally, we discuss triplet superconducting correlations supported by altermagnets and protected by inversion rather than time-reversal symmetry. We apply the nonrelativistic classification of Cooper pairs to describe the interplay between strain and superconductivity in the two-dimensional Lieb lattice and in bulk rutile structures. Here, we show that triplet superconductivity is, on average, unitary in an unstrained altermagnet, but becomes non-unitary under piezomagnetically active strain.

FOS: Physical sciences↗

DOE OSTI2025

Ferroelectric phase transition in group-IV monochalcogenides from an equivariant machine learned force field

Group-IV monochalcogenides are a class of layered ferroelectric semiconductors that have demonstrated spontaneous intrinsic polarization above room temperature. Here, in this study, we use the multi-atomic cluster expansion (MACE) machine learning architecture to train and test a force field capable of modeling the structural properties and second-order ferroelectric-to-paraelectric phase transition in a Group-IV monochalcogenide, GeSe. The model captures the double-well potential energy surface associated with the onset of macroscopic polarization in bulk GeSe within 12.5 meV/atom, as well as near-equilibrium properties like the phonon dispersion. The development of this quantitatively accurate force field enables long-time molecular dynamics simulations, which predict the critical temperature of the ferroelectric-to-paraelectric phase transition in bulk GeSe to be T c = 600 K. This study demonstrates the capabilities of equivariant force-fields to accurately describe phenomena associated with structural symmetry breaking.

ferroelectricity↗

DOE OSTI2026

High-order cumulants and correlation functions near the critical point from molecular dynamics

We present a systematic investigation of particle-number fluctuations in the crossover region near the critical end point of a first-order phase transition using molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our prior studies to third- and fourth-order cumulants in both coordinate- and momentum-space acceptances and integrated correlation functions (factorial cumulants). We find that, even near the critical point, non-Gaussian cumulants equilibrate on timescales comparable to those of the second-order cumulants, but show stronger finite-size effects. The presence of interactions and of the critical point leads to strong deviations of the cumulants from the ideal-gas baseline in coordinate space; these deviations are expected to persist in momentum space in the presence of collective expansion. In particular, the kurtosis becomes strongly negative, κσ 2 ≪ − 1 , on the crossover side of the critical point. However, this signal is significantly diluted once an efficiency cut used to distinguish protons from baryons is applied, leading to |κσ 2 | ≲ 1 even in the presence of the critical point. We discuss our results in the context of ongoing measurements of proton-number cumulants in heavy-ion collisions in RHIC-BES-II.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Correlation-aware binning for small-angle neutron scattering via Gaussian-process inference

Binning in small-angle neutron scattering (SANS) is typically performed empirically, with fixed parameters chosen for convenience rather than statistical optimality. Such practices often fail to balance statistical precision and spatial resolution, leading to inconsistencies across instruments and datasets. Here we establish a correlation-aware framework that determines the optimal bin width from first principles by extending the classical Freedman–Diaconis (FD) rule to account for inter-bin correlations with a Gaussian process. In this formulation, the scattering intensity is treated as a smooth stochastic field whose statistical coherence is described by a covariance matrix. Analytical expressions of errors derived from this model yield closed-form criteria that separate the total deviation into contributions from counting noise, aliasing distortion and curvature-dependent correlation effects. Expressed in reduced variables, the resulting dimensionless error surface reveals a continuous transition from the uncorrelated FD regime to the correlation-dominated limit, providing a unified description of noise suppression and resolution control. Because the formulation depends only on the profile characteristics of scattering intensity I(Q), specifically its average intensity and first- and second-order derivatives, it applies generally to any SANS measurement regardless of sample, instrument or geometry. Experimental validation using small- and ultra-small-angle neutron scattering data confirms the predicted scaling behavior, demonstrating that correlation-aware inference systematically reduces mean-squared error and enables information-efficient reproducible data reduction across materials and instruments.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗

DOE OSTI2025

Temperature and Strain Sensing Characteristics of a 128° YX-Cut LiNbO 3 Rayleigh-Mode SAW Sensor From Room to Cryogenic Temperatures

Accurate, passive, and wireless monitoring of cryogenic hardware is essential for high-energy physics, space propulsion, and biomedical instrumentation. Here, this study quantifies the coupled temperature-strain behavior of Rayleigh-mode surface acoustic-wave (SAW) delay-line sensors fabricated on 128° YX-cut LiNbO 3 . A nonlinear finite element (FE) model incorporating Varshni-based elastic constants, higher-order thermal expansion, and temperature-dependent piezo- and dielectric coefficients was developed and validated experimentally between 280K and 80K. Free-standing (first test condition) and bonded/wired (second test condition) devices exhibited indistinguishable thermal responses; the average temperature coefficient of delay (TCD) in the critical cryogenic range from 130K down to 80K differed by only 0.15 ppm/K (0.32%), confirming that bonding-induced stress is negligible. Over 280-80K the measured TCD was 61.77 ppm/K, while the FE model predicted an equivalent temperature coefficient of frequency (TCF) of −62.74 ppm/K with an overall coefficient of determination R2 = 0.998. In the critical cryogenic interval 130-80K the TCD fell to 47.66 ppm/K, indicating improved thermal stability at low temperature. Controlled loading (0-300 με) revealed a strain coefficient of delay (SCD) that rises from 0.53 ± 0.02 ppm/με at 300K to 1.05 ± 0.02 ppm/με at 80K. This modest sensitivity confirms that, for temperature sensing, strain is a second-order perturbation above 135K but must be compensated at deeper cryogenic levels. Overall, this work establishes a predictive multiphysics model together with repeatable wired measurements that confirm the suitability of SAW sensors for temperature and strain monitoring in extreme cryogenic environments, while also providing a baseline for future wireless implementations.

25 ENERGY STORAGE↗

DOE OSTI2026

Deep learning-assisted modeling for χ (2) nonlinear optics

Modeling second-order (χ(2)) nonlinear optical processes remains computationally expensive due to the need to resolve fast field oscillations and simulate wave propagation using methods such as the split-step Fourier method (SSFM). This can become a bottleneck in real-time applications, such as high-repetition-rate laser systems requiring rapid feedback and control. We present a long short-term memory-based surrogate model trained on SSFM simulations generated from a start-to-end model of the photocathode drive laser at SLAC National Accelerator Laboratory’s Linac Coherent Light Source II. The model achieves over 250× speedup while maintaining high fidelity, enabling future real-time optimization and laying the foundation for data-integrated modeling frameworks and digital twins of laser systems.

Accelerator Physics (physics.acc-ph)↗

DOE OSTI2026

Explicit Monotone Stable Super-Time-stepping Methods for Finite Time Singularities

We explore a novel way to numerically resolve the scaling behavior of finite-time singularities in solutions of nonlinear parabolic PDEs. The Runge–Kutta–Legendre (RKL) and Runge–Kutta–Gegenbauer (RKG) super-time-stepping methods were originally developed for nonlinear complex physics problems with diffusion. These are multistage single step second-order, forward-in-time methods with no implicit solves. The advantage is that the time-step size for stability scales with stage number 𝑠 as $\mathcal{O}$⁡(𝑠 2 ). Many interesting nonlinear PDEs have finite-time singularities, and the presence of diffusion often limits one to using implicit or semi-implicit time-step methods for stability constraints. Finite-time singularities are particularly challenging due to the large range of scales that one desires to resolve, often with adaptive spatial grids and adaptive time steps. Here, in this study, we show two examples of nonlinear PDEs for which the self-similar singularity structure has time and space scales that are resolvable using the RKL and RKG methods, without forcing even smaller time steps. Compared to commonly used implicit numerical methods, we achieve a significantly smaller run time while maintaining comparable accuracy. We also prove numerical monotonicity for both the RKL and RKG methods under their linear stability conditions for the constant coefficient heat equation, in the case of infinite domain and periodic boundary condition, leading to a theoretical guarantee of the superiority of the RKL and RKG methods over traditional super-time-stepping methods, such as the Runge-Kutta-Chebyshev and the orthogonal Runge-Kutta-Chebyshev methods. Code can be found at https://github.com/ZT220501/SRK-Singularity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

A Unidirectional Two-Compartment Neuron Circuit with On-chip STDP learning

Most neuromorphic chips implement the single-compartment point neuron model where synapse circuits connect directly to a leaky integrate and fire (LIF) soma circuit. However, when using a biologically plausible soma circuit (e.g., Hodgkin-Huxley neuron model), an interface circuitry, such as a current conveyor circuit, is needed to transmit synaptic current to the soma circuit. This is especially true for ultra-low power neuron circuits, where membrane capacitance is on the order of 20 fF. This need for an interface circuit arises because the parasitic capacitance and leakage current caused by fabrication mismatch and second-order effects of the output transistors in the synapse circuits can disturb the spiking dynamics of the soma circuit if connected without an interface. Using an interface circuit to isolate the soma’s membrane capacitor from synapses resolves this issue. We propose to use a unidirectional resistor (a transconductance circuit) to connect the synapse and soma circuits instead of conventional current conveyor circuits. Using a biologically plausible spike pattern detection model, we show that the on-chip spike-timing-dependent plasticity (STDP) learning performance of the proposed unidirectional two-compartment neuron circuit is similar to a single-compartment circuit (with a current conveyor as an interface) and additionally, it is more power-efficient and biologically plausible. The chip is fabricated in the Taiwan Semiconductor Manufacturing Company (TSMC) 250 nm technology node and comprises a single neuron circuit.

Gautam, Ashish [ORNL]↗

DOE OSTI2026

Did You Win the GPU Cloud Lottery? Benchmarking from TFLOPS to Tokens/$

Cloud GPUs are commonly assumed to deliver consistent performance for a given GPU model. This assumption does not always hold: cloud providers employ diverse system configurations and virtualization mechanisms, and GPUs themselves exhibit non-negligible manufacturing variability (the silicon lottery). In this work, we present a large-scale measurement study of GPU performance variability across 11 cloud providers, covering over 3,500 physical GPUs and 6,800 benchmark runs. Our hierarchical analysis shows that while execution-level variation stays below 9%, performance varies by up to 38% across devices and providers for the same GPU model. Regression analysis indicates that driver- and OS-related software factors contribute less than 1% of the variance; instead, silicon lottery effects dominate observed performance variation, and cloud providers further amplify them through persistent, systematic second-order effects.

Slynko, Platon [Silicon Data, New York, USA] (ORCI↗

DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗

DOE OSTI2026

Emergent pair density wave order across a Lifshitz transition

We numerically investigate the telltale signs of pair-density-wave order (PDW) in the Kondo-Heisenberg chain by focusing on the momentum resolved spectrum in different parameter regimes. Density matrix renormalization group calculations reveal that this phase is characterized by a dispersion with two minima and four Fermi points, indicating the emergence of an effective next-nearest-neighbor hopping that arises as a second-order effect to avoid magnetic frustration. The pairs appear in the spectrum as in-gap bound states with weight concentrated in the hole pockets. The low-energy physics can be understood by means of a generalized t−J model with next-nearest-neighbor hopping. Our results offer a guide for searching for experimental signatures, and for other models that can realize PDW physics.

Yang, Luhang [University of Tennessee Knoxville]↗

DOE OSTI2025

Hybrid HEFA-HDCJ Process for the Production of Jet Fuel Blendstocks

09 BIOMASS FUELS↗

DOE OSTI2025

Using quantum noise correlation analysis to measure ion temperature

Quantum noise correlation analysis was tested at the proof-of-concept level as a technique to measure ion temperature in a plasma. If eventually successful, this technique could enable a compact, inexpensive, and robust ion temperature diagnostic suitable for a burning plasma environment. Ion temperature is a key parameter determining the fusion performance of a burning plasma, as the fusion cross-section has a strong dependence on ion temperature. This ion temperature diagnostic would require only a small optical view of the plasma through a port to passively record impurity line emission. The instrumentation would be remote from the reactor behind the neutron and bio-shielding. The technique relies solely on quantum correlations of the photons emitted by a plasma impurity to measure ion temperature; there is no grating dispersion of an emission line-width or pulse-height analysis of photon energy. This measurement innovation was tested with instrumentation consisting of two single-photon detectors with high timing resolution, a time-tagging unit, and simple light collection optics. This instrumentation measures the second-order correlation between the light intensity falling on the two detectors. The next steps beyond the proof-of-concept level will be development of diagnostic designs for application of this technique to high-temperature and burning plasmas. Arrays of single-photon avalanche detectors to multiplex measurements of photon correlation will be required to reduce signal integration time to an acceptable duration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Using Separation-Enhanced Isotope Ratio Mass Spectrometry to Enable Increased Renewable Carbon Content in Transportation Fuels (CRADA 525)

Stable isotope ratio measurements of carbon atoms using isotope ratio mass spectrometry (IRMS) can be an effective tool for quantifying biogenic carbon in co-processed fuels, with results approaching the precision and accuracy of accelerator mass spectrometry (AMS). The lower cost of an IRMS may enable deployment to refineries, improving access and analysis turnaround times (≤2 hours), and, by extension, provide data that can allow process optimization to maximize renewable carbon in desired refinery products. This project explored the integration of chemical separation with IRMS analyses to enable highly detailed tracking of biogenic carbon into fuel product streams separated by boiling point range, chemical class, or specific compound. Forty-nine fuels and fuel components of fossil and biogenic origin, spanning gasoline and diesel boiling point ranges, were received from three refiners and were analyzed for their δ 13 C values via IRMS. Results spanned a 13 C range from ca. 10‰ to 44‰ and reflect materials derived from sustainable sources (e.g., C4 or C3 plants, animal-based pathways, syngas) or from fossil-derived fuels. Common ranges are approximately 18‰ to 9‰ and approximately 30‰ to 20‰ for C4 and C3 plants, respectively, and approximately 34‰ to 24‰ and approximately 70‰ to 33‰ for petroleum-derived fuels and methane, respectively. Fuel-like standards were developed and tested using direct-injection elemental analyzer (EA) IRMS for liquid fuels. This method was compared with the published methods, yielding statistically similar results. Four blend curve sets were produced ranging from 0% to 100% of a fuel containing biogenic carbon, focusing on 0% to 10% biogenic carbon. Linear fits were the most applicable for two of the four blend curve sets; however, two sets were found to exhibit slightly quadratic behavior, which was more pronounced in low biogenic blend samples, necessitating second-order fits. The origin of the slight quadratic behavior remains unclear; however, the discussion points to possible interpretations. CanmetENERGY thoroughly characterized a majority of the samples using one- and two-dimensional gas chromatography (GC and GC×GC, respectively) and other analyses. Selected samples were subjected to solid phase extraction (SPE) for saturate, olefin, aromatic, and polar (SOAP) analysis, and the resulting solvent-diluted fractions containing saturates and aromatics were returned to Pacific Northwest National Laboratory (PNNL), where the solvent was removed via evaporation or physical separation using GC techniques. Characterization and separations provided an understanding of saturate and aromatic content, as well as boiling point ranges for each sample and sample fraction. Samples resulting from SPE were examined using EA-IRMS and gas chromatography combustion IRMS (GC-C-IRMS) analyses. Both approaches suggest that the range in values between end-members can be increased by selecting the paraffinic or aromatic fraction of the end-member or by selecting among individual compounds resulting from GC separation of the paraffinic fractions. Considerable work remains to put these approaches into practice and statistically validate the benefit for using a fraction or individual compound over bulk analysis of a sample. However, initial results suggest that separations provide advantages for samples having blend ratios of less than 10% biogenic blendstocks. 13 C results showed statistically similar biofuel blend results to those obtained at PNNL, although additional work is needed to obtain better reproducibility. Select samples were sent to Los Alamos National Laboratory (LANL) for IRMS measurements and Beta Analytics for AMS measurements. This work suggests that IRMS and AMS yield closely comparable results and in some circumstances, IRMS could serve as a surrogate for AMS. While additional work is needed to better resolve statistical advantages for separations and better show the comparable nature of IRMS and AMS in both the biogenic carbon analysis of bulk chemical classes, initial results from this study suggest that these should be pursued in order to proliferate this approach for quantifying biogenic carbon in transportation fuels to the refinery level, thereby potentially enabling process optimization in co-processing scenarios.

09 BIOMASS FUELS↗

DOE OSTI2025

Hybrid HEFA-HDCJ Process for the Production of Jet Fuel Blendstocks

The hydrotreatment of bio-oil derived from the pyrolysis and biocrude from hydrothermal liquefaction of lignocellulosic materials to produce hydrocarbons faces significant technological challenges, mainly due to the high reactivity and poor thermal stability of bio-oil, resulting in the formation of large quantities of coke. This problem has been addressed by existing PNNL patents with a two-step hydrotreatment technology in which the bio-oil is first stabilized with a noble hydrogenation metal (often Pt or Ru). Then, in the second step, the bio-oil is deoxygenated with a Ni-Mo or Co-Mo sulfide catalyst. The main problem with this approach is that the Pt/Ru catalysts deactivate easily in the presence of S or other impurities, which are commonly present in pyrolysis oils. In this project, we explored technological solutions to mitigate coke formation, avoiding the use of Pt/Ru catalysts. Our strategy is based on three actions: (1) Bio-oil stabilization in the presence of alcohols. In this project, we studied the stabilization with butanol. (2) the use of a cosolvent to solubilize the bio-oil. Because coke formation reactions are second-order reactions a reduction in the concentration of reactive bio-oil molecules. In this case, we used yellow greases as a co-solvent. (3) Separation of bio-oil reactive fractions. In this project, we studied the removal of water-soluble fractions. Our batch co-hydrotreatment studies confirmed that the addition of butanol and methanol and the blend with lipids effectively contributed to mitigating coke formation (reducing coke yield to about 1 wt.% %). The removal of sugars did not have a noticeable effect on the overall coke yield, suggesting that coke precursors are present in all bio-oil fractions. Our analytical work suggests that they may be concentrated in the water-insoluble/CH 2 Cl 2 insoluble fractions of pyrolysis oils. Although the technological strategies tested resulted in significant coke reductions, the levels achieved were not sufficiently low to ensure a reliable operation in continuous, fixed-bed trickle-bed reactors. Long runs of more than 100 hours (maximum: 255 h) of co-processing time on stream were achieved in a continuous 40 mL reactor. When the same test was conducted in a larger 400 mL reactor, pressure drop increases associated with coke formation were observed. This increase in coke formation could be due to larger temperature gradients in the bed. Hydrodeoxygenation tests in moving bed reactors and using more active hydrogenation catalysts (for example Ni) could lead to more reliable operations. Unfortunately, our team did not have access to such experimental setups. The technoeconomic analysis suggests that although alcohol use is an effective means to reduce coke formation, the use of alcohol increases production cost. Thus, its use needs to be minimized. A delicate balance needs to be found between the use of technological solutions that allow the reliable operation of the system (stabilization with Ni catalysts, use of small quantities of solvents, processing in moving bed reactors) with a tolerable level of coke formation for the hydrodeoxygenation reactor used and that result in minimum production costs.

09 BIOMASS FUELS↗