Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mean free path”

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 199 records · Page 11

Extrapolation of thermal conductivity in non-equilibrium molecular dynamics simulations to bulk scale

Predictions of the bulk scale thermal conductivity of solids using non-equilibrium molecular dynamics (MD) simulations have relied on the linear extrapolation of the thermal resistivity versus the reciprocal of the system length in the simulations. Several studies have reported deviation of the extrapolation from linearity near the micro-scale, raising a concern of its applicability to large systems. To investigate this issue, present work conducted extensive MD simulations of silicon with two different potentials (EDIP and Tersoff-II) for unprecedented length scales up to 10.3 μm and simulation times up to 530 ns. For large systems ≥0.35 μm in size the non-linearity of the extrapolation of the reciprocal of the thermal conductivity is mostly due to ignoring the dependence of the thermal conductivity on temperature. To account for such dependence, the present analysis fixes the temperature range for determining the gradient for calculating the thermal conductivity values. However, short systems ≤0.23 μm in size show significant non-linearity in the calculated thermal conductivity values using a temperature window of 500 ± 10 K from the simulations results with the EDIP potential. Since these system sizes are shorter than the mean phonon free path in EDIP (~0.22 μm), the nonlinearity may be attributed to phonon transport. For the MD simulations with the Tersoff-II potential there is no significant non-linearity in the calculated thermal conductivity values for systems ranging in size from 0.05 to 5.4 μm.

42 ENGINEERING↗

Grain Boundary-Limited Thermal Transport in Suspended Thin Graphite across an Unexplored Thickness Regime

Here, we present systematic thermal conductivity (κ) measurements of suspended thin graphite ribbons, 234–527 nm thick, using a four-probe 3ω method. Unlike recent reports of phonon hydrodynamics and exceptionally high κ in micrometer-thick graphite ( Science, 2020), we observe significantly lower κ and no signatures of collective phonon flow in this intermediate thickness regime. Instead, our measured κ lies between few-layer graphene and bulk graphite. These results agree with a first-principles-informed Peierls–Boltzmann transport model with spatially resolved Monte Carlo sampling. Additionally, the temperature for the peak κ shifts lower with increasing thickness, due to the interplay of phonon-boundary and phonon-isotope scattering. Incorporating grain boundary scattering into simulations is necessary to replicate the experimental trends. These findings delineate the boundary between ballistic, hydrodynamic, and diffusive transport regimes in graphite and underscore the dominant role of disorder and geometry in phonon transport in quasi-two-dimensional materials, offering insights for nanoscale thermal management.

Boltzmann transport model↗

Ghost imaging second harmonic generation microscopy

Second harmonic generation (SHG) microscopy is useful for visualizing interfaces and sub-structures within a wide range of materials due to the propensity for SHG to occur in non-centrosymmetric environments. However, since SHG is a nonlinear process generally necessitating small focal sizes for higher peak powers, a raster scanning approach is usually needed to build an SHG image over a significant sample size. While raster scanning is effective, there is a cost in terms of the time needed to acquire the image and, also, some materials cannot withstand the higher optical intensities within the small focal volume. As such, we describe a SHG microscopy approach based on ghost imaging (GI), which enables imaging data to be collected in parallel rather than sequentially as in raster scanning techniques. We experimentally demonstrate the approach and combine GI-SHG with compressive sensing to make further substantial gains in reducing the amount of sampling required for image reconstruction. Furthermore, GI-SHG is shown to have significant advantages for imaging in highly scattering environments, partly because GI is a background-free approach requiring spatial correlations between photons that travel two paths, with one path entirely devoid of sample interaction. This basic property of GI means that only the photons that travel unimpeded through the sample preserve the spatial correlations needed to reconstruct the image, while more scattered photons do not contribute to the overall GI signal. Finally, we compare the image quality and sampling properties of three different reconstruction algorithms used for compressive sensing.

47 OTHER INSTRUMENTATION↗

Field-driven ion pairing dynamics in concentrated electrolytes

Here, we investigate ion pairing dynamics in electrolytes driven far from equilibrium using molecular simulations and nonequilibrium rate theory. Focusing on 0.5M LiPF 6 in water and acetonitrile under uniform electric fields, we compute transition path theory observables, including reactive fluxes and mean first-passage times of ion pairing. Moreover, we introduce a dynamical proxy of free-ion population, where its field-induced change is strongly correlated with the nonlinear enhancement of conductivity, yielding an increase of 40% at 50 mV/Å in acetonitrile, compared to that of less than 10% in aqueous electrolytes. Further kinetic analysis elucidates that Onsager’s classical theory substantially overestimates field-induced enhancement of ion pair dissociation in molecular electrolytes. This discrepancy arises from solvent-mediated dynamical pathways and field-induced dielectric decrement that suppress ion pair dissociation within explicit solvents, highlighting that a faithful description of molecular details is essential. Our results provide a molecular interpretation of nonlinear electrolyte transport beyond continuum theories and establish a general framework for quantifying nonequilibrium reaction kinetics in condensed phase systems.

Moon, Seokjin [University of California, Berkeley,↗

Actinium–DOTA coordination in water from hybrid ML/MM: Structure, free energies, and water-exchange pathways

Quantitative simulation of trivalent ƒ-block chelates in water remains challenging because bonded and non-bonded force-field models make different approximations for coordination structure, exchange dynamics, and ion–ligand interactions in highly charged systems. Here, we develop a hybrid machine-learning/molecular-mechanics (ML/MM) framework for Ac 3+ –DOTA in explicit solvent by training an E(3)-equivariant neural network potential (MACELES) on mechanically embedded QM/MM data for Ac aquo and Ac–DOTA species and coupling it to NAMD 2.14 with particle-mesh Ewald electrostatics. Nanosecond ML/MM trajectories remain numerically stable and preserve chelate integrity, yielding a compact DOTA inner shell with an inner-sphere water coordination number of CN Ac,O w ≈ 1.7 arising from a dynamic equilibrium between one- and two-water states (37.5% and 59.9% of frames; three waters 2.5%). A 5 ns potential of mean force shows two low-lying basins at CN Ac,O w ≈ 1 and CN Ac,O w ≈ 2. DFT end-state free energies are consistent with the ML/MM profile, and DFT minimum-energy paths provide a qualitative electronic-structure reference for the observed basin connectivity. State-resolved kinetics reveal picosecond water-exchange pathways that couple hydration changes to transient DOTA arm fluctuations, and training-set comparisons show that temperature-matched Ac–DOTA data optimize energy/force accuracy while more diverse solvated data improve charge prediction. Overall, the present hybrid ML/MM model provides a practical description of Ac 3+ –DOTA hydration thermodynamics and short-time exchange behavior in explicit water at MD-like cost.

Actinium↗

Nitrogen derived from Combined Algal Processing supports algae cultivation for biofuels

Algae biomass, a biofuel feedstock with potential to reduce global CO 2 emissions, is currently too expensive Nitrogen, a significant input and life cycle cost, requires research on recycling and additional sources to reduce costs. Integral to nitrogen recycling is the knowledge of actual nitrogen requirements in outdoor ponds, thus data from the Unified Field Studies of the Algae Testbed Public-Private Partnership was analyzed. The data showed that Nannochloropsis oceanica used 14–24 mg/L NH 4 + per harvest cycle for robust growth. Nitrogen recycled directly from the algae conversion process will likely be the primary source of nitrogen for algae cultivation. Harvested algae can be fractionated via Combined Algal Processing (CAP), where algae are pretreated, fermented, and the fuel feedstocks (e. g. ethanol, lipids) are removed leaving a nitrogen rich slurry. Chlorella vulgaris grew robustly in nitrogen-free media supplemented with CAP residues through 4 harvest cycles with no inhibition. The amount of nitrogen as free ammonium available for recycling directly without further processing was insufficient (~4% of total needed), thus the nitrogen fixing bacteria, Azotobacter vinelandii, was investigated as a means for increasing ammonium concentration in CAP residues. This organism was shown to be proficient at growth on sugars liberated from algae biomass during pretreatment and sustained algae growth using A. vinelandii supernatants was also demonstrated. Using Azotobacter to increase ammonium concentrations on-site using an inexpensive biological process to generate the necessary ammonium for outdoor algae cultivation improves sustainability of algae biofuels by avoiding the energy intensive steps of typical ammonium production and transport. Furthermore, these results show that byproducts of algae conversion begin to close nitrogen demand and, with further supplementation by diazotrophic bacteria or protein hydrolysis, may completely replace current energyintensive and costly industrial nitrogen sources offering a path to increased economic feasibility and improved life cycle analysis of algae-based biofuels.

09 BIOMASS FUELS↗

AmeriFlux FLUXNET-1F CA-SCC Scotty Creek Landscape

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site CA-SCC Scotty Creek Landscape. This is the FLUXNET version of the carbon flux data for the site CA-SCC Scotty Creek Landscape produced by applying the standard ONEFlux (1F) software. Site Description - The Scotty Creek flux tower is located in an organic-rich boreal forest-wetland landscape about 50 km south of Fort Simpson in the Taiga Plains of the Mackenzie watershed. The tower was installed in 2013 and operates an open-path EC system year-round running on solar power only. Flux footprints contain about 50 % forested peat plateaus and 50 % wetlands (i.e., collapse-scar bogs). The forests are underlain by permafrost, while the treeless wetlands are permafrost-free. The tower itself is located on a forested peat plateau. Black spruce tree density on plateaus is sparse and the mean canopy height is ca. 5 m.

Sonnentag, Oliver↗

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Determination of Stress-Free Lattice Spacing (D₀) for Residual Stress Relaxation Measurement in Ni-based Superalloys by Neutron Diffraction

The residual stress and its relaxation in critical engineering components are key to structural materials reliability. The determination of residual stress in bulk engineering structures like Ni-based superalloys can be conducted by measuring the lattice strain with penetrating neutron diffraction. The state-of-the-art engineering diffraction beamline VULCAN at the Oak Ridge National Laboratory (ORNL) Spallation Neutron Source (SNS) provides the critical capability to evaluate the residual strain/stress relaxation in engineering component thanks to the high flux and event data features. By the means of residual stress/strain calculation, it is critical to accurately determine the stress-free lattice spacing d₀, which can be altered by the change of chemistry, unrelieved stress, and measurement scheme. Here we reported d₀ determination for determining residual stress relaxation and distribution in Alloy 718 superalloys after different quenching treatments. Selection of locations for d₀ measurement was discussed by considering the neutron path, attenuation and sample alignment. A dynamic d₀ that resulted from atom diffusion and chemistry change was estimated as function of time for the in-situ relaxation characterization. It demonstrated how the dynamic d₀ values may influence on the strain calculation in different thermally treated samples. The results highlight the importance of dynamic d₀ for in-situ relaxation involving chemical changes and provide guidance to the dynamic d₀ measurement and simulation.

36 MATERIALS SCIENCE↗

Dotted-line FLEET for two-component velocimetry

Femtosecond laser electronic excitation tagging (FLEET) is a powerful unseeded velocimetry technique typically used to measure one component of velocity along a line, or two or three components from a dot. In this Letter, we demonstrate a dotted-line FLEET technique which combines the dense profile capability of a line with the ability to perform two-component velocimetry with a single camera on a dot. Our set-up uses a single beam path to create multiple simultaneous spots, more than previously achieved in other FLEET spot configurations. We perform dotted-line FLEET measurements downstream of a highly turbulent, supersonic nitrogen free jet. Dotted-line FLEET is created by focusing light transmitted by a periodic mask with rectangular slits of 1.6 × 40 mm 2 and an edge-to-edge spacing of 0.5 mm, then focusing the imaged light at the measurement region. Up to seven symmetric dots spaced approximately 0.9 mm apart, with mean full-width at half maximum diameters between 150 and 350 µm, are simultaneously imaged. Both streamwise and radial velocities are computed and presented in this Letter.

47 OTHER INSTRUMENTATION↗

Coupling High-Throughput Experiments and Regression Algorithms to Optimize PGM-Free ORR Electrocatalyst Synthesis

Over the past decades, significant improvement has been achieved in the performance of platinum group metal-free (PGM-free) materials as an alternative to Pt-based electrocatalysts for oxygen reduction reaction (ORR). However, further progress in ORR activity requires evaluation of precursors and synthesis approaches. In response to this challenge, we generated a first of its kind experimental data set of 36 samples using high-throughput synthesis and activity measurements. Several control parameters (e.g., Fe precursor identity, the precursor content, and pyrolysis temperature) were varied. We then developed several state-of-the-art machine learning (ML) based regression models to predict ORR activity, dependent on selected synthesis variables. Through an iterative algorithm, higher prediction accuracy (smaller root-mean-square error) was achieved. We identified that gradient boosting regression (GBR) and support vector regression (SVR), among several methods, work best for this data set. Aided by our ML-based surrogate models, we decided to alter catalyst synthesis conditions, which resulted in a 36% increase in measured ORR activity in comparison to the maximum ORR mass activity value of 21.9 A/g catalyst in the original data set. Overall, this combined experiment and machine learning approach represents a promising path forward toward developing highly efficient next-generation ORR electrocatalysts and, more generally, functional materials.

25 ENERGY STORAGE↗

Observed and Simulated Variability of Droplet Spectral Dispersion in Convective Clouds Over the Amazon

Abstract In this study, the variability of the spectral dispersion of droplet size distributions (DSDs) in convective clouds is investigated. Analyses are based on aircraft measurements of growing cumuli near the Amazon basin, and on numerical simulations of an idealized ice‐free cumulus. In cleaner clouds, the relative dispersion , defined as the ratio of the standard deviation to the mean value of the droplet diameter, is negatively correlated with the ratio of the cloud water content ( ) to the adiabatic liquid water content ( ), while no strong correlation between and is seen in polluted clouds. Bin microphysics numerical simulations suggest that these contrasting behaviors are associated with the effect of collision‐coalescence in cleaner clouds, and secondary droplet activation in polluted clouds, in addition to the turbulent mixing of parcels that experienced different paths within the cloud. Collision‐coalescence simultaneously broadens the DSDs and decreases , explaining the inverse relationship between and in cleaner clouds. Secondary droplet activation broadens the DSDs but has little direct impact on . The combination of a rather modest DSD broadening due to weak collision‐coalescence with enhanced droplet activation in both diluted and highly undiluted cloud regions may contribute to maintain a relatively uniform within polluted clouds. These findings can be useful for parameterizing the shape parameter ( ) of gamma DSDs in bulk microphysics cloud‐resolving models. It is shown that emulating the observed relationship improves the estimation of the collision‐coalescence rate in bulk microphysics simulations compared to the bin simulations.

54 ENVIRONMENTAL SCIENCES↗

Correlations between conduction electrons in dense plasmas

Most treatments of electron-electron correlations in dense plasmas either ignore them entirely (random phase approximation) or neglect the role of ions (jellium approximation). In this work, we go beyond both these approximations to derive a formula for the electron-electron static structure factor which properly accounts for the contributions of both ionic structure and quantum-mechanical dynamic response in the electrons. The result can be viewed as a natural extension of the quantum Ornstein-Zernike theory of ionic and electronic correlations, and it is suitable for dense plasmas in which the ions are classical and the conduction electrons are quantum-mechanical. Here, the corresponding electron-electron pair distribution functions are compared with the results of path integral Monte Carlo simulations, showing good agreement whenever no strong electron resonance states are present. We construct approximate potentials of mean force which describe the effective screened interaction between electrons. Significant deviations from Debye-Hückel screening are present at temperatures and densities relevant to high-energy density experiments involving warm and hot dense plasmas. The presence of correlations between conduction electrons is likely to influence the electron-electron contribution to the electrical and thermal conductivity. It is expected that excitation processes involving the conduction electrons (e.g., free-free absorption) will also be affected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Revisiting the pseudo-supercritical path method: An improved formulation for the alchemical calculation of solid–liquid coexistence

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid–liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. Furthermore, a detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology.

25 ENERGY STORAGE↗

Emission Inventories from Natural Gas Storage Facilities using Regional Frequency Comb Laser Monitoring and Aircraft Flyovers. Final Report

The project centers around a ground-based regional methane sensor developed by the CU/NIST team members under the DOE ARPA-E MONITOR program, and sophisticated aircraft measurement techniques developed by the UCD/Scientific Aviation team (responsible for a recent Science publication quantifying the Aliso Canyon storage field release). The dual frequency-comb spectrometer is an invisible, eye-safe laser capable of measuring atmospheric methane concentrations along beam paths 1+ miles in length with high precision and stability (<1 ppb methane over 1 mile). The CU/NIST team has developed a unique approach to using a single, central spectrometer to locate and size methane emissions as small as 6 scfh from specific gas field structures within the 1+ mile range of the laser. Deployment involves the stationing of small cubic mirrors at strategic locations across the gas storage site. The mirrors direct light from the spectrometer back to a detector, where methane concentrations are recorded based on light absorption at specific wavelengths. The measurements are coupled with high resolution gas transport models based on meteorological measurements (wind, etc.) to determine the precise location and emission rate of methane sources. The calibration-free and continuous nature of dual frequency-comb spectrometer measurements means that a large area can be monitored continuously for temporal variability of emissions over long periods of time. The aircraft measurements involve spiral flights around storage facilities with a sensitive ethane/methane measurement package. Upwind/downwind concentration comparisons and wind data provide total facility emissions data. During this project, a new micrometeorological package was integrated onto the Scientific Aviation aircraft to provide even greater emissions quantification capability. In this project, the ground-based measurement system was deployed for 11 months (to capture seasonal variability) at a West Coast storage facility, together with a campaign of twice monthly aircraft total facility emissions flights. The data from these two modalities was combined to quantify the methane emissions from the facility with first-of-its-kind temporal and spatial detail. Concurrent with this deployment, the aircraft team performed flights to quantify emissions at a wider array of previously un-surveyed storage sites. The ground-based measurement system was then moved to an additional site of differing total storage and delivery capacity. Aircraft mass balance flights continued, repeatedly surveying (twice monthly) this site, as well as surveying a large number of previously un-assessed sites spanning a variety of size, reservoir type and age. Assessment of ground-based, temporally and spatially detailed methane emissions information alongside aircraft-derived methane total facility emissions information led to an improved understanding of emissions, relevant for inventories of the natural gas storage sector. In particular, this work provided temporal detail (seasonal of emissions and frequency of leaks) to stakeholders and personnel responsible for updates to the EPA’s GHGI. The impact of the proposed work has been a dramatic improvement in the temporal, spatial, and site- and type- specific detail and quantification of emissions inventories for the natural gas storage sector. This work directly benefits the EPA’s GHGI, directly ameliorates the environmental impacts and public health and safety of the natural gas storage sector, and aids policymakers and industry to make sound choices with respect to management and regulation.

03 NATURAL GAS↗

Experimental trials of species-specific bat flight responses to an ultrasonic deterrent

Unintended consequences of increasing wind energy production include bat mortalities from wind turbine blade strikes. Ultrasonic deterrents (UDs) have been developed to reduce bat mortalities at wind turbines. Our goal was to experimentally assess the species-specific effectiveness of three emission treatments from the UD developed by NRG Systems. We conducted trials in a flight cage measuring approximately 60 m × 10 m × 4.4 m (length × width × height) from July 2020 to May 2021 in San Marcos, Texas, USA. A single UD was placed at either end of the flight cage, and we randomly selected one for each night of field trials. Trials focused on a red bat species group (Lasiurus borealis and Lasiurus blossevillii; n = 46) and four species: cave myotis (Myotis velifer; n = 57), Brazilian free-tailed bats (Tadarida brasiliensis; n = 73), evening bats (Nycteceius humeralis; n = 53), and tricolored bats (Perimyotis subflavus; n = 17). The trials occurred during three treatment emissions: low (emissions from subarrays at 20, 26, and 32 kHz), high (emissions from subarrays at 38, 44, and 50 kHz), and combined (all six emission frequencies). We placed one wild-captured bat into the flight cage for each trial, which consisted of an acclimation period, a control period with the UD powered off, and the three emission treatments (order randomly selected), each interspersed with a control period. We tracked bat flight using four thermal cameras placed outside the flight cage. We quantified the effectiveness of each treatment by comparing the distances each bat flew from the UD during each treatment vs. the control period using quantile regression. Additionally, we conducted an exploratory analysis of differences between sex and season and sex within season using analysis of variance. Broadly, UDs were effective at altering the bats’ flight paths as they flew farther from the UD during treatments than during controls; however, results varied by species, sex, season, and sex within season. For the red bat group, bats flew farther from the UD during all treatments than during the control period at all percentiles (p < 0.001), and treatments were comparable in effectiveness. For cave myotis, all percentile distances were farther from the UD during each of the treatments than during the control, except the 90th percentile distance during high, and low was most effective. For evening bats and Brazilian free-tailed bats, results were inconsistent, but high and low were most effective, respectively. For tricolored bats, combined and low were significant at the 10th–75th percentiles, high was significant at all percentiles, and combined was most effective. Results suggest UDs may be an effective means of reducing bat mortalities due to wind turbine blade strikes. We recommend that continued research on UDs focus on low emission treatments, which have decreased sound attenuation and demonstrated effectiveness across the bat species evaluated in this study.

17 WIND ENERGY↗

A spall and diffraction study of nanosecond pressure release across the iron ε-α phase boundary

The extreme response of polycrystalline iron at high pressures and high strain rates is revealed by means of high-power laser pulses. The compression portion of the pulse coupled with x-ray diffraction identifies the expected body-centered cubic (α) to hexagonal close packed (ε) displacive transformation. Upon release, observation shows that the complete reverse transformation takes approximately 8 ns and that the structure returns to its initial microstructural configuration, in a reversible transformation path. This is in good agreement with molecular dynamics (MD) simulations which predict an inverse dependence between transformation time and strain rate. The grain size is reduced from μm to nm range during compression and begins increasing back to the original grain size on decompression. The kinetics of the transition is dictated by heterogenous nucleation as it follows the Johnson-Mehl-Avrami-Kolmogorov equation with the appropriate time exponent of ~1. This is confirmed by MD simulations which also identify profuse twinning and dislocation generation. Here, the tensile pulse generated upon reflection at the free surface is captured by time-resolved free surface velocity measurements from which a peak tensile stress of 7 GPa is obtained, in stark contrast with its quasi-static value of ~200 MPa. At these strain rates, the strength of grain interiors, which is determined by twinning and slip exceeds the strength of the boundaries, and failure initiates preferentially in the latter.

36 MATERIALS SCIENCE↗

Conceptual Basis and Techno-Economic Modeling for Integrated Algal Biorefinery Conversion of Microalgae to Fuels and Products (2019 NREL TEA Update: Highlighting Paths to Future Cost Goals via a New Pathway for Combined Algal Processing)

The report documents the conceptual basis for a new potential Combined Algal Processing design strategy which may allow more flexibility in accommodating different algal biomass feedstock compositions, by enabling upgrading of both protein and carbohydrates in a single step, without a strict requirement for either component to be in soluble or monomeric form, while maintaining effective wet lipid extraction techniques to enable high lipid recoveries. In light of previously-established constraints around algal biomass costs (which are significantly higher than lignocellulosic terrestrial biomass), the present CAP processing strategy reflects an integrated biorefinery concept producing both fuels and value-added chemical coproducts as a means to improve profitability and generate coproduct revenues to help drive down the minimum fuel selling price (MFSP) towards economically viable levels. Namely, this report highlights an integrated CAP biorefinery process and associated technical targets that would be required to achieve U.S. Department of Energy target MFSP goals of $2.5/gallon gasoline equivalent by 2030. This is accomplished by a process involving low-cost seasonal storage of algal biomass during high-growth seasons, rapid flash hydrolysis pretreatment of the biomass, solvent extraction of pretreated biomass, cleanup and fractionation of lipids into triglyceride and free fatty acid fractions, and a series of thermochemical conversion steps to upgrade carbohydrates and protein to hydrocarbon fuels. These steps include mild oxidative treatment (MOT), a process originally investigated at NREL for upgrading lignin, followed by catalytic ketonization and hydrotreating of MOT products to fuels. Isolated triglycerides are sent to a coproduct train, with the base case focused on upgrading to polyurethane foams as a high-value, high-market-volume coproduct.

09 BIOMASS FUELS↗