Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Forming Limit Curve”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Direct Observations of Solute Dispersion in Rocks With Distinct Degree of Sub‐Micron Porosity

Abstract The transport of chemical species in rocks is affected by their structural heterogeneity to yield a wide spectrum of local solute concentrations. To quantify such imperfect mixing, advanced methodologies are needed that augment the traditional breakthrough curve analysis by probing solute concentration within the fluids locally. Here, we demonstrate the application of asynchronous, multimodality imaging by X‐ray computed tomography (XCT) and positron emission tomography (PET) to the study of passive tracer experiments in laboratory rock cores. The four‐dimensional concentration maps measured by PET reveal specific signatures of the transport process, which we have quantified using fundamental measures of mixing and spreading. We observe that the extent of solute spreading correlate strongly with the strength of subcore‐scale porosity heterogeneity measured by XCT, while dilution is enhanced in rocks containing substantial sub‐micron porosity. We observe that the analysis of different metrics is necessary, as they can differ in their sensitivity to the strength and forms of heterogeneity. The multimodality imaging approach is uniquely suited to probe the fundamental difference between spreading and mixing in heterogeneous media. We propose that when multi‐dimensional data is available, mixing and spreading can be independently quantified using the same metric. We also demonstrate that one‐dimensional transport models have limited predictive ability toward the internal evolution of the solute concentration, when the model is solely calibrated against the effluent breakthrough curves. The data set generated in this study can be used to build realistic digital rock models and to benchmark transport simulations that account deterministically for rock property heterogeneity.

Kurotori, Takeshi [Department of Chemical Engineer↗
DOE OSTI2024

Modeling Planar Electrodes and Zero–Gap Membrane Electrode Assemblies for CO 2 Electrolysis

Multiphysics modeling enables probing of conditions inside a CO 2 electrolyzer that are difficult to measure, such as local concentrations and pH, as well as rapid testing of possible design changes. A one-dimensional model for a zero-gap membrane electrode assembly (MEA) CO 2 electrolyzer was developed with the assumption that catalyst layers interact with the membrane ionomer such that the ionomer affects the underlying kinetics. The kinetics for bicarbonate reacting to form hydrogen are fit using a planar electrode model for silver with an ionomer coating. The MEA model results are validated against experimental studies for current density and product selectivity. Flooding of the cathode is modeled using saturation curves, and results show that blocked pores in the microporous layer play a significant role in limiting the mass transport at high potentials (>2.8 V). Sensitivity studies showed that CO Faradaic efficiency can be increased by decreasing catalyst layer thickness and porosity, and decreasing KHCO 3 concentration.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Binary pseudo-random array standards for calibration of 3D optical surface profilers used for metrology with aspheric x-ray optics

High-accuracy surface metrology is vitally important in manufacturing ultra-high-quality free-form mirrors designed to manipulate x-ray light with nanometer-scale wavelengths. The current and potential capabilities of x˗ray mirror manufacturing are limited by inherent imperfections of the integrated metrology tools. Metrology tools are currently calibrated with super-polished flat test-standard/reference mirrors. This is acceptable for fabrication of slightly curved x-ray optics. However, for even moderately curved aspherical x-ray mirrors the flat-reference calibration is not sufficiently accurate. For micro-stitching interferometry developed for surface measurements with curved x-ray mirrors, the tool aberration errors are known to be transferred into the optical surface topography of x-ray mirrors. Our approach to improving metrology is to thoroughly calibrate the measuring tool and apply the results of the calibration to deconvolution of the measured data. Here we explore the application of a recently developed technique for calibrating the instrument transfer function (ITF) of 3D optical surface profilers to metrology with significantly curved x-ray optics. The technique, based on test standards patterned with two-dimensional (2D) binary pseudo-random arrays (BPRAs), employs the unique properties of the BPRA patterns in the spatial frequency domain. The inherent 2D power spectral density of the pattern has a deterministic white-noise-like character that allows direct determination of the ITF with uniform sensitivity over the entire spatial frequency range and field of view of an instrument. The high efficacy of the technique has been previously demonstrated in application to metrology with flat and slightly curved optics. Here, we concentrate on development of an efficient fabrication process for production of highly randomized (HR) BPRA test standards on flat and 500-mm spherical optical substrates. We also compare and discuss the results of the ITF calibration of an interferometric microscope when using the HR BPRA standards on flat and curved substrates.

Munechika, Keiko↗
DOE OSTI2021

The Encounter of the Parker Solar Probe and a Comet-like Object Near the Sun: Model Predictions and Measurements

The Parker Solar Probe (PSP) aims to explore the nascent solar wind close to the Sun. Meanwhile, PSP is also expected to encounter small objects like comets and asteroids. Here, we survey the ephemerides to find the chance of a recent encounter and then model the interaction between released dusty plasmas and solar wind plasmas. On 2019 September 2, a comet-like object, the 322P/Solar and Heliosphere Observatory, just passed its perihelion flying to a heliocentric distance of 0.12 au and swept by PSP at a relative distance as close as 0.025 au. We present the dynamics of the dust particles released from 322P, forming a curved dust tail. Along the path of PSP in the simulated inner heliosphere, the states of plasma and magnetic field are sampled and illustrated, with the magnetic field sequences from simulation results being compared directly with the in situ measurements from PSP. Through the comparison, we suggest that 322P might be at a deficient activity level releasing limited dusty plasmas on its way to becoming a "rock comet." We also present images of solar wind streamers as recorded by the Wide-field Imager for Solar Probe Plus, showing an indication of dust bombardment for the images superposed with messy trails. We observe from the Large Angle and Spectrometric Coronagraph that 322P was transiting from a dimming region to a relatively bright streamer during its perihelion passage, and perform a simulation to confirm that 322P was flying from relatively faster to slower solar wind streams, modifying the local plasma states of the streams.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Assessment of irradiation embrittlement effect on fatigue life of a pressurized-water reactor pressure vessel using the fracture toughness master curve approach

The reactor pressure vessel (RPV) is a critical structural component in pressurized-water reactors, and it is designed to withstand extreme conditions, such as high pressures, elevated temperatures, and prolonged radiation exposure. Ensuring RPV integrity is essential for the safe and reliable long-term operation of nuclear power plants, especially as aging mechanisms such as fatigue and irradiation embrittlement pose increased risks. Fatigue, caused by cyclic thermal and mechanical loading, can lead to crack initiation in localized high-stress regions. Simultaneously, neutron irradiation, particularly in the beltline region, progressively reduces fracture toughness, increasing susceptibility to brittle fracture. These combined effects of fatigue and irradiation embrittlement potentially impact the RPV structural integrity, necessitating fitness-for-service assessments. This study applies the fracture toughness master curve approach to evaluate the impact of irradiation embrittlement on RPV fatigue life. A coupled thermo-mechanical stress analysis identifies critical stress locations under normal service transients, pinpointing regions most vulnerable to fatigue crack initiation and growth. Stress intensity factors for postulated flaws at these locations are calculated, enabling an assessment of fatigue life under irradiated and unirradiated conditions. The results indicate that neutron irradiation embrittlement accelerates the conditions in which a critical crack can form and lead to failure, particularly at lower temperatures. The failure occurs where reduced fracture toughness limits the material’s resistance to crack growth. Axial cracks at nozzle corners are the most life limiting without irradiation, while circumferential cracks demonstrate longer fatigue lives. The findings highlight the importance of incorporating irradiation effects into fatigue life predictions to ensure the long-term structural integrity of RPVs.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

U-splines: Splines over unstructured meshes

U-splines are a novel approach to the construction of a spline basis for representing smooth objects in Computer-Aided Design (CAD) and Computer-Aided Engineering (CAE). A spline is a piecewise-defined function that satisfies continuity constraints between adjacent cells in a mesh. U-splines differ from existing spline constructions, such as Non-Uniform Rational B-splines (NURBS), subdivision surfaces, T-splines, and hierarchical B-splines, in that they can accommodate local variation in cell size, polynomial degree, and smoothness simultaneously over more varied mesh configurations. Mixed cell types (e.g., triangle and quadrilateral cells in the same mesh) and T-junctions are also supported, although the continuity of interfaces with triangle and tetrahedral cells is limited in the present work. The U-spline algorithm introduces a new technique for using local null space solutions to construct basis functions for the global spline null space problem. The U-spline construction is presented for curves, surfaces, and volumes with higher dimensional generalizations possible. Lastly, a set of requirements are given to ensure that the U-spline basis is positive, forms a partition of unity, is complete, and is locally linearly independent.

42 ENGINEERING↗
DOE OSTI2023

Effective Landau-type model of a Hf x Zr 1 - x O 2 -graphene nanostructure

Here, to describe the charge-polarization coupling in the nanostructure formed by a thin Hf x Zr 1-x O 2 film with a single-layer graphene as a top electrode, we develop the “effective” Landau-Ginzburg-Devonshire model. This approach is based on the parametrization of the Landau expansion coefficients for the polar (FE) and antipolar (AFE) orderings in thin Hf 1-x Zr x O 2 films from a limited number of polarization-field curves and hysteresis loops. The Landau expansion coefficients are nonlinearly dependent on the film thickness h and Zr/[Hf+Zr] ratio x, in contrast to h-independent and linearly x-dependent expansion coefficients of a classical Landau energy. We explain the dependence of the Landau expansion coefficients by the strong nonmonotonic dependence of the Hf 1-x Zr x O 2 film polar properties on the film thickness, grain size and surface energy. The proposed Landau free energy with five “effective” expansion coefficients, which are interpolation functions of x and h, describes the continuous transformation of polarization dependences on applied electric field and hysteresis loop shapes induced by the changes of x and h in the range 0 < x < 1 and 5 nm < h < 35 nm. Using the effective free energy, we demonstrated that the polarization of Hf 1-x Zr x O 2 film influences strongly on the graphene conductivity, and the full correlation between the distribution of polarization and charge carriers in graphene is revealed. In accordance with our modeling, the polarization of the (5 – 25) nm thick Hf 1-x Zr x O 2 films, which are in the ferroelectric-like or antiferroelectric-like states for the chemical compositions 0.35 ≤ x ≤ 0.95, determine the concentration of carriers in graphene and can control its field dependence. The result can be promising for creation of next generation Si-compatible nonvolatile memories and graphene-ferroelectric FETs, because the working voltages applied to the Hf 1-x Zr x O 2 film (which acts as a gate) can be relatively low (less than 2 V). These low voltages are sufficient to induce the pronounced hysteresis of ferroelectric polarization in the Hf 1-x Zr x O 2 gate, which, due to the strong electric coupling, induces the hysteresis of the graphene charge.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ultra-stable and poison tolerant oxygen evolution activity enabled by surface In 2 O 3- x (OH) y of Co 3 In 2 S 2 large single crystals

Water is an earth-abundant source for clean hydrogen production via electrochemical water electrolysis (WE). However, the surface poisoning that occurs in aqueous electrolytes drastically deactivates the electrocatalytic performance of electrodes. Here, in this study, we report electrochemically formed In 2 O 3-x (OH) y on the surface of a large (1–1.5 mm long, 0.5–0.6 mm wide and 0.3–0.5 mm thick) single crystal of Weyl semimetal Co 3 In 2 S 2 (Co 3 In 2 S 2 /In 2 O 3-x (OH) y ) as an ultra-stable and poison tolerant electrode for the oxygen evolution reaction (OER) in 1 M KOH, addressing a bottleneck in WE. The OER activity of the powder form of Co 3 In 2 S 2 is limited by its aerophilic nature. Remarkably, the single-crystal electrodes maintained their high activity for a continuous operational period of 5 h in 1 M KOH electrolyte with/without 10 mM strong surface-poisoning ligands i.e., potassium cyanide, bipyridine, and ethylenediaminetetraacetate disodium salt. The electrodes exhibited stable OER activity for 1000 h at 100 mA cm -2 (1.73 V vs. RHE). The temperature-dependent OER polarization curves (10–70 °C) unambiguously revealed surface poisoning through the suppression of precatalytic Co-redox peaks on the bipyridine poisoned electrode, which led to the stabilization of surface Co-sites. The X-ray photoelectron spectroscopy analyses of pristine, poisoned and post-electrocatalytic single-crystal Co 3 In 2 S 2 electrodes revealed the existence of an In 2 O 3-x (OH) y surface phase, which could be the potential heterostructure for the origin of ultra-stable and poison tolerant OER activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Pattern formation and transport for externally driven active matter on periodic substrates (a)

Abstract We investigate the transport of interacting active run-and-tumble particles moving under an external drift force through a periodic array of obstacles for increasing drive amplitudes. For high activity where the system forms a motility-induced phase-separated state, there are several distinct dynamic phases including a low drive pinned cluster phase, an intermediate uniform fluid, and a higher drive stripe crystal state. The transitions between the phases are correlated with signatures in the transport curves, differential mobility, and power spectra of the velocity fluctuations. In contrast, in the low activity regime the transport curves and power spectra undergo little change as a function of drive. We argue that in the high activity limit, the behavior is similar to that of driven solids on periodic substrates, while in the low activity limit the system behaves like a driven fluid.

36 MATERIALS SCIENCE↗
DOE OSTI2023

SN 2022ann: a Type Icn supernova from a dwarf galaxy that reveals helium in its circumstellar environment

Here, we present optical and near-infrared (NIR) observations of the Type Icn supernova (SN Icn) 2022ann, the fifth member of its newly identified class of SNe. Its early optical spectra are dominated by narrow carbon and oxygen P-Cygni features with absorption velocities of ∼800 km s −1 ; slower than other SNe Icn and indicative of interaction with a dense, H/He-poor circumstellar medium (CSM) that is outflowing slower than typical Wolf–Rayet wind velocities of >1000 km s −1 . We identify helium in NIR spectra 2 weeks after maximum and in optical spectra at 3 weeks, demonstrating that the CSM is not fully devoid of helium. Unlike other SNe Icn, the spectra of SN 2022ann never develop broad features from SN ejecta, including in the nebular phase. Compared to other SNe Icn, SN 2022ann has a low luminosity (o-band absolute magnitude of ∼−17.7), and evolves slowly. The bolometric light curve is well-modelled by 4.8 M ⊙ of SN ejecta interacting with 1.3 M ⊙ of CSM. We place an upper limit of 0.04 M ⊙ of 56 Ni synthesized in the explosion. The host galaxy is a dwarf galaxy with a stellar mass of 10 7.34 M ⊙ (implied metallicity of log(Z/Z ⊙ ) ≈ 0.10) and integrated star-formation rate of log (SFR) = −2.20 M ⊙ yr −1 ; both lower than 97 per cent of galaxies observed to produce core-collapse supernovae, although consistent with star-forming galaxies on the galaxy Main Sequence. The low CSM velocity, nickel and ejecta masses, and likely low-metallicity environment disfavour a single Wolf–Rayet progenitor star. Instead, a binary companion is likely required to adequately strip the progenitor and produce a low-velocity outflow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Potential bill impacts of dynamic electricity pricing on California utility customers

The rapid growth of renewable generation is creating challenges for the California grid in the form of the “duck curve,” with increasingly steep ramping required for conventional generation resources in the morning and evening, and growing curtailment of solar resources in midday periods. Time-varying electricity tariffs have received considerable attention as a tool to address these challenges, with a renewed recent focus on the potential for dynamic tariffs that vary to reflect conditions on the grid in near-real time. Consideration of dynamic tariffs may raise concerns about the financial impact on utility customers, especially for those who have limited flexibility to modify their electricity consumption in response. Specific areas of concern include electricity bills, bill volatility, and equity implications related to cost shifting among customer groups. In this paper we leverage smart meter data for more than 400,000 California utility customers, spanning residential, commercial, industrial, and agricultural customers, to assess potential customer bill impacts arising from a multi-component dynamic tariff . Specifically, we compute impacts on customer bills and bill volatility under the assumption of fully inelastic demand, i.e., where customers do not change their consumption patterns in response to the tariff. We also assess various approaches designing subscription load shapes that customers can pre-purchase as a hedge that may provide a measure of protection against large negative impacts, while still incentivizing the modification of loads on the margin. We compare and contrast the relative impacts on different customer classes and discuss benefits and pitfalls of different dynamic tariff structures and subscription load shapes.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

High pressure melt line of nickel using a generalized embedded atomic method potential

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. Further, the potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: $T$ m = $1969.23+19.15P-0.012P$ 2 . In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dislocation‐Driven Formation of Oriented Macroperiodic Metastructures of Curved Single Crystal Lattices in Glass

Abstract Single crystals fabricated in glass by localized heating can develop uniquely deformed lattices stabilized by the surrounding amorphous medium. The development of lattice curvature appears to be intrinsic to the crystal growth process in some systems, while the result of the locally changing crystallography in others. In this work, a model laser‐fabricated rotating lattice Sb 2 S 3 crystal grown in stoichiometric glass is used to demonstrate fabrication of novel macroperiodic metastructures that utilize intrinsic lattice curvature superimposed with subtle crystallographic influences. The limited availability of slip systems drives the lattice curvature magnitude to vary with crystal growth direction, maximizing for lattices aligned with the predominant Burgers vector along with corresponding increases in dislocation density. Misaligned lattice orientations form smaller secondary lattice curvatures arising from misaligned Burgers vectors with further elastic contributions. Over extended crystal growth, these secondary components align the lattice to rotate about either the <001> or <010> crystal axes forming repeating metastructures of lattice orientation with periodicity 20–160 microns in length. The mechanistic approach used in this work may be expanded to other systems with known slip systems to better understand and design macroperiodic metastructures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On plasticity-enhanced interfacial toughness in bonded joints

Here, the performance and reliability of many structures and components depend on the integrity of interfaces between dissimilar materials. Interfacial toughness Γ is the key material parameter that characterizes resistance to interfacial crack growth, and Γ is known to depend on many factors including temperature. For example, previous work showed that the toughness of an epoxy/aluminum interface decreased 40 % as the test temperature was increased from -60 °C to room temperature (RT). Interfacial integrity at elevated temperatures is of considerable practical importance. Recent measurements show that instead of continuing to decrease with increasing temperature, Γ increases when test temperature is above RT. Cohesive zone finite element calculations of an adhesively bonded, asymmetric double cantilever beam specimen of the type used to measure Γ suggest that this increase in toughness may be a result of R-curve behavior generated by plasticity-enhanced toughening during stable subcritical crack growth with interfacial toughness defined as the critical steady-state limit value. In these calculations, which used an elastic-perfectly plastic epoxy model with a temperature-dependent yield strength, the plasticity-enhanced increase in Γ above its intrinsic value Γ o depended on the ratio of interfacial strength σ* to the yield strength σyb of the bond material. There is a nonlinear relationship between Γ/Γ o and σ*/σ yb with the value Γ/Γ o increasing rapidly above a threshold value of σ*/σ yb . The predicted increase in toughness can be significant. For example, there is nearly a factor of two predicted increase in Γ/Γ o during micrometer-scale crack-growth when σ*/σ yb = 2 (a reasonable choice for σ*/σ yb ). Furthermore, contrary to other reported results, plasticity-enhanced toughening can occur prior to crack advance as the cohesive zone forms and the peak stress at the tip of the original crack tip translates to the tip of the fully formed cohesive zone. These results suggest that plasticity-enhanced toughening should be considered when modeling interfaces at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Advanced Functional Materials ZnO@ZIF-8 Heterostructures for Neuromorphic Synaptic Devices

Metal-Organic Frameworks have shown great promise for memristive synaptic devices because of their adaptable electrical properties and inherent porosity, which help in efficient ion transport and storage of functional molecules. However, the limited conductivity and sensitivity to moisture associated with most MOF materials hinder their use in electronic applications. In this study, MOFs derived from metabolites have been synthesized, using ZnO as a self-sacrificial metallic source to grow Zeolitic Imidazolate Framework-8 (ZIF-8), forming ZnO@ZIF-8 heteronanostructures for neuromorphic applications. These hybrids ZnO@ZIF-8 act as promising candidates for advanced electronic applications. Through systematic investigations, including a comparison of ZnO@ZIF-8 bilayer devices with ZnO single-layer counterparts, a significant discovery emerged. The bilayer synaptic devices outperformed their single layer counterparts in electrical characteristics, as evidenced by Current-Voltage curve analysis and Long-Term Potentiation/Long Term Depression (LTP/LTD) measurements. Additionally, the ZnO layer thickness (10, 30, 50 nm) has been emphasized and the effects of annealing at 300°C in a vacuum were carefully studied for Neuromorphic Synaptic Devices. Remarkably, the 10 nm ZnO layer exhibited the lowest set voltage and a consistent On/Off ratio. This work marks a significant advancement in understanding bilayer synaptic devices, highlighting the crucial role of ZIF-8 in shaping the future of neuromorphic systems

Kolleboyina, Jayaramulu↗
DOE OSTI2026

Prediction of Creep-Induced Strain Using a Symbolic Regression-Based Model

Material creep under high-temperature conditions limits the lifetime and safety of structural systems such as advanced nuclear reactors. Conventional creep testing is slow and often produces inconsistent results across nominally identical experiments, making lifetime prediction uncertain. Here, to address these challenges, this work develops a data-driven symbolic regression (SR) model that consolidates results from duplicate creep tests and predicts the remaining strain-time curve of an ongoing experiment. The method uses piece-wise multi-objective SR with physical constraints to generate analytic, interpretable functions describing transient creep strain. Applied to Inconel Alloy 617 data, the approach achieved relative mean absolute errors of 1.0–9.5%, providing closed-form predictions of strain evolution. These results demonstrate a first step toward reducing the duration and cost of long-term creep testing while retaining physically interpretable model forms.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Enhanced Carbon Capture Behavior of Carbon Fibers via Ionic Liquid Modification

Carbon-based materials are widely deployed in carbon capture but are only limited to physisorption procedures. Further extending functionalized carbon materials to CO 2 chemisorption under low CO 2 concentration is highly desirable yet challenging. In this study, a carbon fiber composed of solely ultra-micropores (~0.45 nm) was deployed as the precursor to avoid the pore blocking effect, which was modified by superbase-derived ionic liquids (ILs) containing strong interaction sites with CO 2 . By forming a thin coating layer on the surface, the as-afforded surface-functionalized fiber materials demonstrated enhanced CO 2 uptake capacity and improved CO 2 sorption kinetics, as evaluated by both the volumetric method and thermogravimetric analysis, as well as the calculated energy distribution curves. The achievements made in this work provide guidance on the functionalization of carbon-based materials towards enhanced CO 2 chemisorption by forming a thin layer of selected IL coating.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗