SEARCH · Engineering Papers

Results for “effective diffusivity coefficient”

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

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Comprehensive Dynamics in a Polyelectrolyte Complex Coacervate

The linear viscoelastic response, LVR, of a hydrated polyelectrolyte complex coacervate, PEC, was evaluated over a range of frequencies, temperatures, and salt concentrations. The PEC was a nearly stoichiometric blend of a quaternary ammonium poly([3-(methacrylamido)propyl]trimethylammonium chloride), PMAPTAC, and poly(2-acrylamido-2-methyl-1-propanesulfonic acid sodium salt), PAMPS, an aliphatic sulfonate, selected because they remain fully charged over the conditions of use. Narrow molecular weight distribution polyelectrolytes were prepared using fractionation techniques. A partially deuterated version of PMAPTAC was incorporated to determine its coil radius of gyration, R g , within PECs using small-angle neutron scattering. Chain dimensions were determined to be Gaussian with a Kuhn length of 2.37 nm, which remained constant from 25 to 65 °C. The LVR for a series of matched molecular weight PECs, mostly above the entanglement threshold, exhibited crossovers of modulus versus frequency classically attributed to the reptation time, relaxation between entanglements, and the relaxation of a Kuhn length of units (the “monomer” time). The scaling for zero shear viscosity, η 0 , versus chain length, N, was η 0 ~ N 3.1 , in agreement with “sticky reptation” theory. The lifetime and activation energy, E p , of a pair between polyanion and polycation repeat units, Pol + Pol – , were determined from diffusion coefficients of salt ions within the PEC. The activation energy for LVR of salt-free PECs was 2E p , showing that the key mechanism limiting the dynamics of undoped PECs is pair exchange. An FTIR technique was used to distinguish whether SCN– acts as a counterion or a co-ion within PECs. Doping of PECs with NaSCN breaks Pol + Pol – pairing efficiently, which decreases effective cross-linking and decreases viscosity. Finally, an equation was derived that quantitatively predicts this effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Modeling the Size Distribution and Chemical Composition of Secondary Organic Aerosols during the Reactive Uptake of Isoprene-Derived Epoxydiols under Low-Humidity Condition

Reactive uptake of isoprene epoxydiols (IEPOX), which are isoprene oxidation products, onto acidic sulfate aerosols is recognized to be an important mechanism for the formation of isoprene-derived secondary organic aerosol (SOA). While a mechanistic understanding of IEPOX-SOA formation exists, several processes affecting their formation remain uncertain. Evaluating mechanistic IEPOX-SOA models with controlled laboratory experiments under longer atmospherically relevant timescales is critical. Here, we implement our latest understanding of IEPOX-SOA formation within a box model to simulate the measured reactive uptake of IEPOX on polydisperse ammonium bisulfate seed aerosols within an environmental Teflon chamber. The model is evaluated with single-particle measurements of size distribution, volume, density, and composition of aerosols due to IEPOX-SOA formation at timescales of hours. Additionally, we find that the model can simulate the growth of particles due to IEPOX multiphase chemistry, as reflected in increases of the mean particle size and volume concentrations, and a shift of number size distribution to larger sizes. The model also predicts the observed evolution of particle number mean diameter and total volume concentrations at the end of the experiment. We show that in addition to self-limiting effects of IEPOX-SOA coatings, accounting for the molar balance between inorganic and organic sulfate and mass accommodation coefficient of IEPOX are important parameters governing modeling of IEPOX-SOA formation. Thus, models which do not account for the molar sulfate balance and/or diffusion limitations within IEPOX-SOA coatings are likely to predict too high IEPOX-SOA formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Turbulent drifts of impurity ions as an explanation for anomalous radial transport in the far-SOL of DIII-D

Abstract Successful fusion reactor operation relies on minimal core contamination by impurities, otherwise too much power may be radiated and harm performance. This requires reliable predictions of impurity transport from the scrape-off layer (SOL) into the core, beyond the traditional ‘anomalous’ diffusion approach. We report a set of far-SOL tungsten transport simulations that demonstrate the role of turbulent drifts on radial impurity transport. A turbulent plasma background is simulated using the gyrokinetic SOL code Gkeyll. Tungsten ions are followed within the plasma background using only their drifts. We find that tungsten tends to travel radially outwards with velocities between v r = 300–1200 m s −1 primarily due to polarization drift. We also extract an anomalous radial diffusion coefficient that varies from D r anom = 5–20 m 2 s −1 . These results are compared to and agree with previous interpretive modeling results. We also show how the turbulent polarization drift can transport some tungsten ions from the wall inwards with effective pinch velocities up to 10 000 m s −1 . We conclude that turbulent drifts are a likely explanation for historically anomalous radial impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Avalanche gain and its effect on energy resolution in GEM-based detectors

Here, we present avalanche gain and associated resolution measurements recorded with a 4 He:CO 2 (70:30) gas mixture and pure SF 6 , a Negative Ion (NI) gas. SF 6 is of particular interest to the directional dark matter detection community, as its low thermal diffusion helps to retain recoil ionization track features over long drift lengths. With the aid of a general form of the reduced first Townsend coefficient (RFTC), multiple GEM-based detector data sets are used to study the high-gain behavior of the 4 He:CO 2 gas mixture. The high-gain data is well described purely in terms of the reduced electric field strength and the number of GEMs, and the robust relationship between the RFTC and the average, reduced, electric field strength across the GEMs is emphasized. The associated (pulse-height) resolution measurements are used to discuss the variance of the avalanche distribution and to describe and estimate the lower limits of energy resolution one should expect to measure using a simple relationship with the RFTC. In the end, a description of avalanche gain, its effect on energy resolution, and the contributing experimental parameters in GEM-based detectors is developed over a broad parameter space for further use.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

The physics basis to integrate an MHD stable, high-power hybrid scenario to a cool divertor for steady-state reactor operation

Abstract Coupling a high-performance core to a low heat flux divertor is a crucial step for ITER and a Fusion Pilot Plant or DEMO. Experiments in DIII-D recently expanded the steady-state hybrid scenario to high density and divertor impurity injection to study the feasibility of a radiating mantle solution. This work presents the physics basis for trade-offs between density, current drive and stability to tearing modes (TMs) at high β. EC power is crucial to tailor the plasma profiles into a passively stable state, and to eject impurities from the core. Off-axis EC depositions decrease the heating efficiency, but calculated electron heat transport coefficients show that this effect is partially mitigated by improved confinement inside the EC deposition. Additionally, the reduction in pressure is recovered by increasing the density. This favourable scaling of confinement with density was observed in high power plasmas for years, and this work provides a comprehensive explanation. ELITE predictions indicate that a path in peeling-ballooning stability opens up for certain conditions of density, power, q 95 and shaping, allowing the edge pressure to continue increasing without encountering a limit. In the core, calculated anomalous fast-ion diffusion coefficients are consistent with density fluctuation measurements in the toroidicity-induced Alfvén eigenmode range, showing that smaller fast-ion losses contribute to the enhanced confinement at high density. The edge integration study shows that divertor heat loads can be reduced with Ne and Ar injection, but this eventually triggers a cascade of n = 1, 2, 3 core TMs. We can now show that impurity radiation in the core is small and it is not the cause for the drop in confinement at high Ar and Ne injection rates. The overlap between the core TMs is consistent with the loss of pressure as estimated by the Belt model for the coupled rational surfaces. Optimization of these trade-offs has achieved plasmas with sustained H 98y2 = 1.7, f GW = 0.7 and ∼85% mantle radiation. The scenario and its variations at higher density and on- vs off-axis EC heating has been studied as a candidate for an integrated solution for several reactor designs, such as ITER, ARC, and the ARIES-ACT1 case, showing promising results in terms of fusion power and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Structuring, stochastic behavior, and charge storage capacity of redox-active microemulsions formulated with mixtures of toluene and ionic liquid as oil phase

Oil/water microemulsions (µEs) are promising electrolytes for redox flow batteries (RFBs) as they simultaneously improve charge capacity and ionic conductivity. Here, we report the successful formulation of bicontinuous µEs where the oil phase is a solution of trihexyltetradecylphosphonium bis(trifluoromethylsulfonyl)amide ionic liquid in toluene with redox-active ferrocene. We examined the effect of supporting electrolyte anion (NO 3 – , Cl – , ClO 4 – ) on the structure, reactivity, transport, and electrolytic performance of redox-active µEs. Neutron scattering and nuclear magnetic resonance showed that the domain size increased as Cl – > NO 3 – > ClO 4 – while the ferrocene diffusion coefficient increased as NO 3 – > Cl – > ClO 4 – . Scanning electrochemical microscopy indicated anion-dependent current fluctuations during electrolysis, with ClO 4 – exhibiting the least high-frequency oscillations, which correlate to the highest charge and discharge capacity and reversibility. All ionic liquid containing systems improved the performance of toluene-based µEs, highlighting new design principles for these electrolytes.

diffusion↗

DOE OSTI2026

Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas

Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.

Plasma kinetic theory↗

DOE OSTI2025

Optimizing energy yield of monolithic perovskite/silicon tandem solar cells in real-world Conditions: The impact of luminescent coupling

Efficient light management is key to maximizing power conversion efficiency (PCE) in monolithic perovskite/silicon tandem solar cells. Achieving peak efficiency requires closely matched current generation in all junctions, especially in integrated configurations. However, real-world conditions vary significantly due to factors such as sunlight spectrum, diffuse-to-direct sunlight ratio, angular distribution of light, subcell temperature coefficients, and ground reflection. This study introduces a comprehensive optical and device simulation to optimize perovskite/silicon tandem cells, considering experimental luminescent coupling (LC) efficiency and its dependence on working conditions, alongside variations in radiative recombination, effect of temperature on absorptivity spectra, and cloud cover. Our results show potential energy yield improvements of up to 1.4 % with LC, based on current perovskite radiative recombination records, and up to 4 % with direct bandgap materials. Although radiative recombination's dependence on excitation intensity reduces output power and requires thicker absorbers, LC compensates for these losses. LC also lowers the optimized bandgap for the perovskite top cell from 1.72 eV to 1.64–1.68 eV, or even lower in regions with redshifted irradiance. Additionally, optimization revealed that thinner silicon bottom cells require a higher perovskite top cell bandgap, impacting the balance between fabrication cost and cell stability.

14 SOLAR ENERGY↗

DOE OSTI2024

Unraveling Interdiffusion Phenomena and the Role of Nanoscale Diffusion Barriers in the Copper–Gold System

Diffusion is one of the most fundamental concepts in materials science, playing a pivotal role in materials synthesis, forming, and degradation. Of particular importance is solid state interdiffusion of metals which defines the usable parameter space for material combinations in the form of alloys. This parameter space can be explored on the macroscopic scale by using diffusion couples. However, this method reaches its limit when going to low temperatures, small scales, and when testing ultrathin diffusion barriers. Therefore, this work transfers the principle of the diffusion couples to small scales by using core–shell nanowires and in situ heating. This allows us to delve into the interdiffusion dynamics of copper and gold, revealing the interplay between diffusion and the disorder–order phase transition. Our in situ TEM experiments in combination with chemical mapping reveal the interdiffusion coefficients of Cu and Au at low temperatures and highlight the impact of ordering processes on the diffusion behavior. The formation of ordered domains within the solid-solution is examined using high-resolution imaging and nanodiffraction including strain mapping. In addition, we examine the effectiveness of ultrathin Al 2 O 3 barrier layers to control interdiffusion of the diffusion couple. Our findings indicate that a 5 nm thick layer serves as an efficient diffusion barrier. Furthermore, this research provides valuable insights into the interdiffusion behavior of Cu and Au on the nanoscale, offering potential applications in the development of miniaturized integrated circuits and nanodevices.

alloys↗

DOE OSTI2024

Sorption kinetics and stability of conventional adsorbents for mercury remediation

In-situ remediation of mercury at numerous contaminated sites worldwide is a challenging and costly endeavor due to the persistency of this contaminant. In this study, we evaluated eight commercially available sorbent media ranging from carbon-, clays- and silica-based materials (PBC– Biochar, eSorb – Sorbster, nsPAC – Powdered Activated Carbon, fsPAC – Powdered Activated Carbon with Mackinawite, F300 – Filtrasorb 300, Si-SH – Silica Thiol, eBind – RemBind, Q-Clay – Organoclay PM-199), for their effectiveness in sorbing mercury (Hg 2+ ) and mercury complexed with dissolved organic matter (Hg-DOM). Under the chosen experimental conditions of this study, results showed that in the absence of DOM, the kinetic rates of Hg 2+ sorption onto the evaluated sorbents were in the order of 0.31 min-1 (Si-SH) to 2.98 min -1 (nsPAC), whereas in the presence of DOM, the rates varied from 0.16 min-1 (F300) to 0.95 min -1 (nsPAC). The measured sorption capacity for Hg 2+ in the absence of DOM varied from 3.02 mg/g (Q-Clay) to 35.15 mg/g (Si-SH), whereas in the presence of DOM, calculated partition coefficient (KD) ranged from 69.7 mL/g (Q-Clay) to 41,510 mL/g (Si-SH). Furthermore, kinetic data suggest liquid film diffusion was the rate-limiting steps governing mercury sorption onto the studied media. Overall, the obtained study parameters (kinetics/isotherm) are particularly important in informing robust engineering designs for deployment of the vast majority of evaluated sorbents. Thus, sorbent-based strategies offer viable solutions for cost-effective cleanup of mercury at industrially contaminated sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Quantitative nanoscale tracking of oxygen vacancy diffusion inside single ceria grains by in situ transmission electron microscopy

Oxygen vacancy formation and migration in ceria is critical to its electrochemical and catalytic properties in systems for chemical and energy transformation, but its quantification is rather challenging especially at atomic-scale because of disordered distribution. Here we report a rational approach to track oxygen vacancy diffusion in single grains of pure and Sm-doped ceria at -20 °C to 160 °C using in situ (scanning) transmission electron microscopy ((S)TEM). To create a gradient in oxygen vacancy concentration, a small region (~30 nm in diameter) inside a ceria grain is reduced to the C-type CeO 1.68 phase by the ionization or radiolysis effect of a high-energy electron beam. The evolution in oxygen vacancy concentration is then mapped through lattice expansion measurement using scanning nano-beam diffraction or 4D STEM at a spatial resolution better than 2 nm; this allows direct determination of local oxygen vacancy diffusion coefficients in a very small domain inside pure and Sm-doped ceria at different temperatures. Further, the activation energies for oxygen transport are determined to be 0.59, 0.66, 1.12, and 1.27 eV for pure CeO 2 , Ce 0.94 Sm 0.06 O 1.97 , Ce 0.89 Sm 0.11 O 1.945 , and Ce 0.8 Sm 0.2 O 1.9 , respectively, implying that activation energy increases due to impurity scattering. The results are qualitatively supported by density functional theory (DFT) calculations. In addition, our in situ TEM investigation reveals that dislocations impede oxygen vacancy diffusion by absorbing oxygen vacancies from the surrounding areas and pinning them locally. With more oxygen vacancies absorbed, dislocations show extended strain fields with local tensile zone sandwiched between the compressed ones. Therefore, dislocation density should be reduced in order to minimize the resistance to oxygen vacancy diffusion at low temperatures.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Tritium diffusion and formation in the bulk and defective surface of γ-LiAlO2 pellets: First-principles investigation

In tritium-producing burnable absorber rods (TPBAR), γ-LiAlO2 is used in the form of an annular ceramic pellet enriched with the 6Li isotope. When irradiated in a pressurized water reactor (PWR), the 6Li pellets absorb neutrons and produce tritium (3H) through 6Li + n 3H + α. The 3H chemically reacts with the metal getter where it is captured and leads to formation of a metal hydride. For TPBARs to enable effective tritium production in PWRs, we investigated the 3H diffusion pathways in the bulk and surface of γ-LiAlO2 with different concentrations of lithium defects. The calculated results for bulk and low-index surfaces, thermal conductivity, 3H activation energy barriers, and the 3H diffusion coefficients in γ-LiAlO2 are in good agreement with the available experimental values. In the bulk, our results show that the smallest activation energy barrier is 0.63 eV for substitutional 3H diffusion with a diffusion coefficient of 3.25x10-12 m2/s. After 3H diffused from bulk to the surface, it could form different species (such as 3H2, 3H2O, C3H4) depending on the surface structure, vacancy types and the impurity carbon. Our results indicate that the 3H2 is the main product from γ-LiAlO2 pellets. As the number of VLi vacancies and 3H atoms increases under irradiation (3H-rich condition), 3H2O release could increase from the surface of γ-LiAlO2 pellets.

Jia, Ting↗

DOE OSTI2025

Charm and bottom hadrons in hot hadronic matter

Heavy quarks, and the hadrons containing them, are excellent probes of the QCD medium formed in high-energy heavy-ion collisions, as they provide essential information on the transport properties of the medium and how quarks color-neutralize into hadrons. Large theoretical and phenomenological efforts have been dedicated thus far to assess the diffusion of charm and bottom quarks in the quark–gluon plasma and their subsequent hadronization into heavy-flavor (HF) hadrons. However, the fireball formed in heavy-ion collisions also features an extended hadronic phase, and therefore any quantitative analysis of experimental observables needs to account for the rescattering of charm and bottom hadrons. This is further reinforced by the presence of a QCD cross-over transition and the notion that the interaction strength is maximal in the vicinity of the pseudo-critical temperature. We review existing approaches for evaluating the interactions of open HF hadrons in a hadronic heat bath and the pertinent results for scattering amplitudes, spectral functions and transport coefficients. While most of the work to date has focused on D -mesons, we also discuss excited states as well as HF baryons and the bottom sector. Both the HF hadro-chemistry and bottom observables will play a key role in future experimental measurements. We also conduct a survey of transport calculations in heavy-ion collisions that have included effects of hadronic HF diffusion and assess its impact on various observables.

effective field theories↗

DOE OSTI2025

Multiphase Processing of the Water-Soluble and Insoluble Phases of Biomass Burning Organic Aerosol

Biomass burning is one of the most significant sources of organic aerosol in the atmosphere. Biomass burning organic aerosol (BBOA) has been observed to undergo liquid– liquid phase separation (LLPS) to give core–shell morphology with the hydrophobic phase encapsulating the hydrophilic phase, potentially impacting the evolution of light-absorbing components, i.e., brown carbon (BrC), through multiphase processes. Here, we demonstrate how multiphase processing differs between the watersoluble (i.e., hydrophilic) and insoluble (i.e., hydrophobic) phases of BBOA in terms of reactive uptake of ozone in a coated-wall flow tube. Effects of relative humidity (RH) and ultraviolet (UV) irradiation were investigated. Experimental timeseries were used to inform simulations using multilayer kinetic modeling. Among non-irradiated thin films, the uptake coefficient was greatest for the water-soluble phase at 75% RH (3 × 10 –5 , corresponding to a diffusion coefficient of BrC, D BrC , of 3 × 10 –9 cm 2 s –1 ) and least for the same phase at 0% RH (1 × 10 –5 , corresponding to D BrC of 1 × 10 –10 cm 2 s –1 ). The uptake coefficient for the water-insoluble phase fell between these two (about 1.5 × 10 –5 ), regardless of RH, and the corresponding D BrC increased only slightly (8 × 10 –10 cm 2 s –1 at 0% RH to 9 × 10 –10 cm 2 s –1 at 75% RH). The uptake coefficients of both phases at 0% RH decreased significantly after UV irradiation, consistent with a transition from viscous liquid to solid and supported by qualitative microscopy observations. Modeling multiphase ozone oxidation of primary BrC components in the atmosphere demonstrated, first, that LLPS may extend the lifetime of water-soluble BBOA encapsulated by water-insoluble species by a factor of 1.5 at moderate to high RH and, also, that UV irradiation may extend the lifetime of both phases by more than a factor of 2.5.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

The design of a photoionization front experiment using the Z-Machine as a driving source and estimated measurements

Radiation-driven heat fronts are present in the early universe during reionization, the circumstellar medium of supernovae, and in high-energy-density physics experiments. Dedicated experiments to observe and diagnose the behavior of these types of heat fronts can improve our understanding of these phenomena. A simulation study of photoionization fronts using the HELIOS-CR radiation hydrodynamics code provides an experimental design for the Z-Machine at Sandia National Laboratory using a measurement-calibrated input radiation flux to drive the photoionization front. The simulations use detailed atomic physics and non-diffusive radiation transport in 1D to determine an optimal gas pressure of 0.75 atm for an experiment in N gas as well as the effects of increasing the thickness of the window that seals the gas cell. Post-processing of these simulations demonstrates that ratios of atomic rate coefficients place the heat front in a physics regime where photoionization dominates the energy deposition. Here, to see the sensitivity of the simulations to changes in the model and spatial grid, this analysis performed resolution, atomic model detail, and radiation transport angular grid studies showing less than 10% deviation from the nominal model for increased complexity, when possible. An effort to emulate 3D geometric effects on the radiation flux using an artificial attenuation scheme has shown that, even for conservative estimates of the flux, simulations still produce a photoionization front. Estimations of a streaked, visible spectroscopy measurement using SPECT3D showed that line emission measurements are present early in time and that later in time thermal emission should become dominant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Space-Time Finite Element Tensor Network Approach for the Time-Dependent Convection–Diffusion–Reaction Equation with Variable Coefficients

In this paper, we present a new space-time Galerkin-like method, where we treat the discretization of spatial and temporal domains simultaneously. This method utilizes a mixed formulation of the tensor-train (TT) and quantized tensor-train (QTT) (please see Section Tensor-Train Decomposition), designed for the finite element discretization (Q1-FEM) of the time-dependent convection–diffusion–reaction (CDR) equation. We reformulate the assembly process of the finite element discretized CDR to enhance its compatibility with tensor operations and introduce a low-rank tensor structure for the finite element operators. Recognizing the banded structure inherent in the finite element framework’s discrete operators, we further exploit the QTT format of the CDR to achieve greater speed and compression. Additionally, we present a comprehensive approach for integrating variable coefficients of CDR into the global discrete operators within the TT/QTT framework. The effectiveness of the proposed method, in terms of memory efficiency and computational complexity, is demonstrated through a series of numerical experiments, including a semi-linear example.

convection–diffusion–reaction equation↗

DOE OSTI2021

On Stochastic Heating and Its Phase-space Signatures in Low-beta Kinetic Turbulence

We revisit the theory of stochastic heating of ions and investigate its phase-space signatures in kinetic turbulence of relevance to low- β portions of the solar wind. In particular, we retain a full scale-dependent approach in our treatment, and we explicitly consider the case in which electric-field fluctuations can be described by a generalized Ohm's law that includes Hall and thermoelectric effects. These two electric-field terms provide the dominant contributions to stochastic ion heating when the ion-Larmor scale is much smaller than the ion skin depth, ρ i << d i , which is the case at β << 1. Employing well-known spectral scaling laws for Alfvén-wave and kinetic-Alfvén-wave turbulent fluctuations, we obtain scaling relations characterizing the field-perpendicular particle-energization rate and energy diffusion coefficient associated with stochastic heating in these two regimes. Phase-space signatures of ion heating are then investigated using three-dimensional hybrid-kinetic simulations of continuously driven Alfvénic turbulence at low β (namely, β i = β e = 0.3 and β i = β e = 1/9). In these simulations, energization of ions parallel to the magnetic field is subdominant compared to its perpendicular counterpart ( Q ∥,i << Q ⊥,i ), and the fraction of turbulent energy that goes into ion heating is ≈75% at β i = 0.3 and ≈40% at β i ≃ 0.1. The phase-space signatures of ion energization are consistent with Landau-resonant collisionless damping and a ( β -dependent) combination of ion-cyclotron and stochastic heating. We demonstrate good agreement between our scale-dependent theory and various signatures associated with the stochastic portion of the heating. We discuss briefly the effect of intermittency on stochastic heating and the implications of our work for the interpretation of stochastic heating in solar-wind spacecraft data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Physically Based Modeling of PEMFC Cathode Catalyst Layers: Effective Microstructure and Ionomer Structure–Property Relationship Impacts

Abstract This work presents a pseudo-two-dimensional proton exchange membrane fuel cell (PEMFC) model incorporating Nafion ionomer structure–property relationships in the cathode catalyst layer (CL) to capture and explain losses at low Pt loading. Structural data from neutron reflectometry and thin film Nafion conductivity measurements predict variations in the oxygen diffusion coefficient and ionic conductivity with changing CL ionomer thickness and Pt loading. By including these structure–property relationships, predicted polarization curves agree closely with previously published experimental data from cells with Pt loadings between 0.025 and 0.2 mg/cm2. Results demonstrate that structure–property relationships based on physically measurable ionomer and CL properties provide a feasible interpretation of PEMFC CL phenomena for a range of Pt loadings and help explain previously unaccounted-for losses at low Pt. Results also show that simulations must account for surface species coverage variations in order to properly capture the kinetic losses. Finally, results suggest that an increase in ionomer thickness surrounding the C/Pt surfaces may lead to improved cell performance due to improved ionic conductivity.

Electrochemistry↗