Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffuse”

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

Combining Distributed Electrochemical Measurements with a Semi-Empirical Model to Identify Local Variations in Overpotentials for PEM Water Electrolyzers

Abstract A segmented cell and PCB approach was used to demonstrate an advanced technique for measuring and understanding the electrochemical behavior of PEM water electrolyzers. By employing a semi-empirical model, discrete modes of overpotential and their spatial distributions were identified in an operating cell. Distributed area-specific resistance (DASR measurements indicate that local dehydration of the polymer electrolyte membrane is the primary driver of increased polarization, behaving similarly to mass transport-limited currents at low flow rates due to ineffective gas-phase bubble removal. To improve quantification, the overall ohmic overpotential is divided into background (η Ohm,b ) and mass transport-coupled (η Ohm,mt ) components. Two types of anode diffusion media, a baseline PTL and a thin foil liquid/gas diffusion layer (LGDL), were examined. The results show significant variations in η Ohm,b and highlight that ineffective bubble removal leads to pronounced local dehydration. The effect of cell temperature on various overpotentials was also examined. Negligible differences in diffusion overpotential were observed when comparing PEWE performance using the PTL and LGDL, despite major differences in two-phase fluid transport within each diffusion media. This underscores the critical role of the catalyst layer's micro-porous structure in diffusion processes. The study highlights the complex nature and intricacies of mass transport limitations within PEWEs.

Roy, Anirban (ORCID:0000000163803386)↗
DOE OSTI2025

Analysis of Silver Release from Furnace-Tested TRISO Particles

High temperature testing of intact TRISO particles previously irradiated in the AGR-5/6/7 experiment was performed in the Furnace for Irradiated TRISO Testing (FITT) to directly confirm diffusive release of silver and europium from intact TRISO particles. Testing was conducted from 1,100–1,600°C for exposure times up to 100 h to directly confirm silver through layer release below safety testing temperatures. The FITT analysis showed highest levels of silver release at 1,300–1,400°C which confirms athermal release behaviors previous observed in step-wise and transient safety tests. Additionally, release was non-uniform with some particles releasing a majority of their inventory while others did not appear to release silver under identical testing conditions, which was consistent with historic observations. An assessment of the effective silver diffusion coefficient in the SiC layer was conducted and indicated maximum values in the 1,300–1,400°C range. The magnitude of the calculated diffusion coefficients also exceeded currently accepted diffusion coefficients. The release behavior of europium from intact particles was also analyzed in FITT for at 1,450–1,550°C for 500 h to 750 h. Direct confirmation of europium release below safety testing temperatures was confirmed absent contributions from matrix release. The analysis indicated europium release follows a general Arrhenius behavior and suggests general uniform release and indicates irradiation conditions influence observed release response. Calculated diffusion kinetics agree well with historic experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Self-Consistent Relativistic Electron Scattering using the Sherlock Scattering Model for X-ray Diagnostics

We present on a new, self-consistent, arbitrary-temperature Romberg integration scheme for modeling electron scattering in materials in a LANL Lagrangian Shock Hydro (LSH) code. Electron beam-target interactions are fundamental to a wide range of scientific and technological applications. When high-energy electron beams hit their target, they may scatter, deposit energy, or ionize the source. These processes govern the behavior and outcomes in nanotechnology manufacturing, electron microscopy, and modern X-ray diagnostics. Simulating these interactions is essential for interpreting experimental results, predicting material responses, and designing efficient tools and experiments. At Los Alamos, this is done using a LSH code, which is a multi-dimension, multi-material, massively parallel, multi-physics code used to simulate applications from asteroid impacts to electron beam interactions. By effectively and efficiently modeling the way that electrons scatter from the beam we can bolster these simulations and more accurately predict experimental outcomes. The model currently implemented in the LSH of interest is based on work by Papp and does not self-consistently preserve momentum in the slightly relativistic regime; here we adopt a model proposed by Braams and Karney and implement a Romberg integration scheme to compute the diffusion tensor. In this paper we will provide background on the Braams-Karney diffusion tensor as well as the Romberg integration scheme we employed to numerically solve for it. We will show that our integration scheme is accurate in solving for the set of scalar potentials used to re-express the diffusion tensor in differential form, and in solving for the diffusion coefficients in the larger LSH code. By using this diffusion tensor rather than the existing Papp one, and numerically integrating it with a Romberg method, we produce much more accurate, self-consistent results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Mercury’s Chaotic Secular Evolution as a Subdiffusive Process

Abstract Mercury’s orbit can destabilize, generally resulting in a collision with either Venus or the Sun. Chaotic evolution can causeg 1 to decrease to the approximately constant value ofg 5 and create a resonance. Previous work has approximated the variation ing 1 as stochastic diffusion, which leads to a phenomological model that can reproduce the Mercury instability statistics of secular andN-body models on timescales longer than 10 Gyr. Here we show that the diffusive model significantly underpredicts the Mercury instability probability on timescales less than 5 Gyr, the remaining lifespan of the solar system. This is becauseg 1 exhibits larger variations on short timescales than the diffusive model would suggest. To better model the variations on short timescales, we build a new subdiffusive phenomological model forg 1 . Subdiffusion is similar to diffusion but exhibits larger displacements on short timescales and smaller displacements on long timescales. We choose model parameters based on the behavior of theg 1 trajectories in theN-body simulations, leading to a tuned model that can reproduce Mercury instability statistics from 1–40 Gyr. This work motivates fundamental questions in solar system dynamics: why does subdiffusion better approximate the variation ing 1 than standard diffusion? Why is there an upper bound ong 1 , but not a lower bound that would prevent it from reachingg 5 ?

Astronomy & Astrophysics↗
DOE OSTI2025

Generative learning of densities on manifolds

A generative modeling framework is proposed that combines diffusion models and manifold learning to efficiently sample data densities on manifolds. The approach utilizes Diffusion Maps to uncover possible low-dimensional underlying (latent) spaces in the high-dimensional data (ambient) space. Two approaches for sampling from the latent data density are described. The first is a score-based diffusion model, which is trained to map a standard normal distribution to the latent data distribution using a neural network. The second one involves solving an Itô stochastic differential equation in the latent space. Additional realizations of the data are generated by lifting the samples back to the ambient space using Double Diffusion Maps , a recently introduced technique typically employed in studying dynamical system reduction; here the focus lies in sampling densities rather than system dynamics. The proposed approaches enable sampling high dimensional data densities restricted to low-dimensional, a priori unknown manifolds. The efficacy of the proposed framework is demonstrated through a benchmark problem and a material with multiscale structure.

Double diffusion maps↗
DOE OSTI2025

Water Dynamics of Superacid Aromatic Proton Exchange Membranes for Fuel Cell Applications

Proton exchange membranes (PEMs) with high conductivity are of critical importance for the development of fuel cells, electrolyzers, and other electrochemical technologies. In this research, poly(1,1,2,2-tetrafluoro-2-phenoxyethane-1-sulfonic acid) (PTPS) with an aromatic polymer main chain and a perfluorinated superacidic polymer side chain was synthesized. The water dynamics of PTPS were characterized across various length scales using a combination of Fourier-transform infrared spectroscopy (FTIR) and nuclear magnetic resonance (NMR) and compared with Nafion, a standard perfluorinated PEM, and sulfonated poly(ether sulfone) (SPES 40), an aromatic PEM without perfluorinated superacid side chains. The T 1 and T 2 relaxation times of water in the samples probed by NMR increase from SPES 40 to PTPS to Nafion, indicating that the local motion of the water molecules becomes faster. This trend corresponds well with the relative fraction of bulk-like water determined using FTIR. At larger length scales, the diffusion coefficient of water was characterized using pulsed-field gradient NMR (PFG-NMR). At a longer diffusion time (Δ = 100 ms), PTPS has a smaller diffusion coefficient compared with both Nafion and SPES 40, due to restricted diffusion, and this effect is also evident in the proton conductivity of the hydrated membranes. From this comparison, it is apparent that the aromatic backbone and side chain type greatly influence the water dynamics in PEMs at various length scales and the water dynamics significantly impact the bulk proton conductivity. These insights will lead to new designs for aromatic PEMs and help to identify bottlenecks in current materials.

25 ENERGY STORAGE↗
DOE OSTI2026

Integrated Strategies for Overcoming Resolution Limits in Electron Beam Lithography of Chemically Amplified Resists

Electron beam lithography (EBL) of chemically amplified resists (CARs) faces fundamental challenges, including stochastic electron scattering and acid diffusion, that limit resolution and reproducibility. Using SU-8 as a model CAR, this study systematically investigated complementary strategies to address these challenges, combining multipass exposure, proximity effect correction (PEC) with midrange correction factors, base quencher incorporation, and post-exposure bake (PEB) suppression. Monte Carlo simulations and calibrated PEC modeling revealed that extending the point spread function to include a midrange scattering component significantly improved critical dimension (CD) control across varying pattern densities, correcting deviations that conventional two-term PEC failed to capture. Multipass exposure, particularly 4-pass writing with a 25% offset, redistributed the dose to average stochastic beam and scattering fluctuations, reducing line-width roughness by more than 50% and yielding more uniform nanoscale features. Photoacid confinement was investigated by adding urea as a base quencher, which successfully reduced acid diffusion but introduced substantial sensitivity penalties without improving ultimate resolution or Z-factor performance, underscoring the trade-offs of chemical versus physical confinement. Suppressing PEB most directly minimized acid diffusion, resulting in improved Z-factors and reproducible 30 nm half-pitch dense line/space patterns. Overall, these results demonstrated that PEC with midrange correction, multipass strategies, quencher additives, and PEB-free processing addresses different aspects of the EBL process window and that their integration provides a comprehensive framework for managing stochastic scattering, diffusion, and chemical amplification effects. This framework advances dense nanoscale patterning in CARs and establishes guiding principles for optimizing resist design and process strategies in high-resolution EBL and potentially other advanced lithographies, such as extreme ultraviolet (EUV) lithography.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Large eddy simulation of a supersonic lifted hydrogen flame: Impacts of Lewis, turbulent Schmidt and Prandtl numbers

Parametric large eddy simulations (LES) of a supersonic lifted hydrogen flame are reported. The emphases are on two aspects: impacts of (1) Lewis number (Lei of the ith species) and (2) turbulent Schmidt and Prandtl numbers (⁠Sc t and Pr t ⁠) on supersonic turbulent flame and flow structures. Five cases are considered: species-specific Le i , Sc t = Sr t = 1.0 (C 0 ); unity Le i , Sc t = Pr t = 1.0 (C 1 ); species-specific Sc t = 0.5, Pr t = 1.0 (C 2 ); species-specific Le i , Sc t = 1.0, Pr t = 0.5 (C 3 ); and species-specific Le i , Sc t = Pr t = 0.5 (C 4 ). Numerical results of instantaneous and/or time-averaged species mole fractions, mixture fraction, heat release rate, flame base location, and mixed modes of premixed and diffusion combustion are compared between cases C 0 and C 1 . Differences in auto-ignition locations and strengths and flame structures and stabilization specify the impacts of Lewis number. They are triggered by different predictions of species mass and thermal diffusions at fuel-coflow and/or coflow-ambient air mixing layers. These differences are rationalized by a scale analysis of mass/thermal diffusion and convection for case C 0 , which suggests the relatively low but non-negligible former against the latter. Cases C 0 and C 2 –C 4 barely see differences in terms of instantaneous and/or time-averaged temperature, velocity, and mixed combustion modes except for further downstream areas where combustion occurs. Both Sc t and Pr t impose less significant influences than Lewis number, as sub-grid scale mass/thermal diffusion is subordinate to its resolved counterpart according to their scale analysis for case C 4 .

42 ENGINEERING↗
DOE OSTI2025

CO 2 in Ionene–Ionic Liquid Composite Membranes

Abstract Ionene – ionic liquid (IL) composites are promising materials for CO 2 separation, yet a molecular‐level understanding of their structure and its impact on CO 2 speciation, solubility, rotation, and diffusivity remains unclear. Herein, using multimodal nuclear magnetic resonance (NMR), time‐of‐flight secondary ion mass spectrometry (ToF‐SIMS), atomic force microscopy (AFM), and molecular dynamics (MD) simulations, we reveal that the composites contain IL‐rich domains extending across hundreds of nanometres within the ionene matrix, and these bicontinuous domains span the entire membrane depth. CO 2 also absorbs into the ionene matrix, with the distribution between two CO 2 species varying with temperature and time. The rotational correlation times of these two species are on the timescale of 0.1 and 1 ns, respectively. As IL content increases, the ionic domains expand, resulting in higher CO 2 solubility due to enhanced molecular dynamics and increased free volume in both ionene backbones and IL‐rich regions. Although CO 2 diffusion in the membranes is an order of magnitude slower than in bulk IL, the activation energy for CO 2 diffusion remains comparable. Ionene‐IL composites represent a promising platform for designing CO 2 separation membranes, offering enhanced CO 2 diffusion and selectivity through IL‐rich domains, and increased CO 2 solubility and mechanical integrity from the ionene matrix.

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

Statistical distributions for transient transport

Here, this paper introduces the use of statistical distributions based on transport differential equations for clear distinction of transport modes within transient kinetic experiments. More specifically, novel techniques are developed for the transient data obtained through the Temporal Analysis of Products (TAP) reactor and are applicable to experiments where pulse response into a gas flow is used. The methodology allows distinguishing between two domains of diffusion transport in heterogeneous catalytic systems, i.e., Knudsen and non-Knudsen diffusion, using statistical fingerprints, and finding the transition domain. Two distribution parameters were obtained that directly result in coefficients that correspond to the concentration and the rate of transport. Using a linear relationship between the rate and concentration coefficients, Knudsen diffusion is revealed when the rate of transport is constant and non-Knudsen diffusion is confirmed when the rate of transport coefficient is a function of the concentration coefficient. As a result, accurate transport information can be extracted from experimental data even in the presence of comprising instrument drift or noise particularly when analyzing higher pressure pulse responses with complex transport. This enables more direct investigation of experiments influenced by gas-phase reactions.

TAP reactor↗
DOE OSTI2025

Mobility assessment of the BCC and carbide phases in the C-Nb, C-U and Nb-U systems

Uranium carbides with refractory metal additions are considered for Gen IV nuclear reactors and nuclear thermal propulsion as fuels for their high-temperature and corrosion resistant properties. Understanding kinetic effects that dictate microstructural evolution during fabrication and operating conditions is essential to advance technological development of these fuels. This work presents the development of an atomic mobility database for C-Nb-U systems based off available experimental data supported with ab-initio methods. The mobility assessments and uncertainty quantification (using Markov chain Monte Carlo) were conducted in the Kawin software. Carbon diffusion is considered dominant, as metal diffusion is much slower, with niobium diffusion being even slower and rate limiting than uranium metal. We provide a comprehensive and self-consistent thermo-kinetic database that is validated by diffusion couple simulations through Kawin. In conclusion, this enables prediction of microstructural and phase evolution critical for the development and lifetime assessment of next generation nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

First in operando transient grating spectroscopy measurements on tungsten during high flux plasma operation in PISCES-RF

This work details the construction and first in operando transient grating spectroscopy measurements conducted on PISCES-RF during plasma operation. A preliminary study on a tungsten sample with varied plasma species (Ar and D 2 ), ion flux (4 × 10 21 to 6 × 10 22 ion/m 2 /s), ion energy (5 to 80 eV), and surface temperature (20 to 200 °C) is presented. Prior to plasma exposure, the thermal diffusivity and surface acoustic wave (SAW) speed, (6.9 ± 0.2) × 10 −5 m 2 /s and (2.66 ± 0.02) km/s, were measured in situ with a thermal grating period of 12.5 μm. During Ar plasma exposure, these two quantities vary according to the sample surface temperature. At low flux and varied ion energy, D 2 plasma exposure also results in thermal diffusivity and SAW speed varying with surface temperature alone. At high flux, D 2 plasma exposure results in no convincing change to the thermal diffusivity, but the SAW speed is reduced. After returning to 20 °C, post plasma exposures, the thermal diffusivity and SAW speed almost recover to the initial value prior to exposure.

Surface acoustic wave SAW↗
DOE OSTI2025

Ultrafast Oxygen Conduction in Sillén Oxychlorides

Oxygen ion conductors are crucial for enhancing the efficiency of various clean energy technologies, including fuel cells, solid oxide air batteries, electrolyzers, membranes, sensors, and more. In this study, a structure-similarity analysis of ≈62k oxygen-containing compounds identified the MBi 2 O 4 X (M = rare-earth element, X = halogen element) family as promising candidates for fast oxygen transport. Among these, LaBi 2 O 4 Cl is found as an ultrafast oxygen conductor with an ultralow migration barrier of 0.1 eV based on ab initio studies. Its 2D layered structure, featuring a “triple fluorite” layer, supports diffusion of both oxygen vacancies and interstitials. In addition to vacancy diffusion with a 0.1 eV barrier, ab initio studies show interstitial diffusion exhibits a modest barrier of 0.6–0.8 eV. Frenkel pairs are found to be the dominant defects in intrinsic LaBi 2 O 4 Cl, facilitating significant vacancy-mediated oxygen diffusion at elevated temperatures. With 2.8% oxygen vacancies, LaBi 2 O 4 Cl is predicted to achieve a conductivity of 0.3 S/cm at 25 °C in a single crystal. Experimental synthesis and characterization of polycrystalline LaBi 2 O 4 Cl and Sr-doped LaBi 2 O 4 Cl revealed conductivity exceeding that of YSZ and LSGM below 400 °C, with lower activation energies, achieving a total conductivity of 0.1−0.2 mS/cm at 300 °C. Here, while these results confirm its potential of fast oxygen transport, we suggest further experimental optimization of LaBi 2 O 4 Cl, including aliovalent doping and microstructure refinement, could significantly enhance its performance, facilitating fast oxygen conduction approaching room temperature.

Defects↗
DOE OSTI2026

Insights into Rotational and Translational Dynamics in Mixtures of Ethylene Glycol and Choline Chloride Using Nuclear Magnetic Resonance Techniques

This work examines molecular dynamics and interactions in ethylene glycol–choline chloride (EG–ChCl) mixtures across 0–33 mol % ChCl, spanning the true eutectic region near 17–20 mol % and the commonly used 1:2 formulation. We combine pulsed-field-gradient (PFG) diffusion, fast-field-cycling (FFC) relaxometry, temperature-dependent 13 C T 1 , and nuclear Overhauser effect spectroscopy (NOESY) to disentangle local from macroscopic dynamics. PFG and FFC show that both translational and average rotational motions largely track the strong increase in viscosity with ChCl content, with ethylene glycol consistently diffusing faster than the choline cation and no global dynamical anomaly at the eutectic composition. More subtle, site-specific composition effects nevertheless emerge. The ratio of the diffusion coefficient of the hydroxyl group of choline to the diffusion coefficient of the methyl group of choline displays a shallow minimum in the 17–25 mol % region, indicating a modest change in how the hydroxyl-bearing end of choline samples the underlying translational motion relative to the methyl groups. 13 C T 1 analysis shows that rotational correlation times at 25 °C generally increase with ChCl, reflecting viscosity-coupled slowing, while the CH 2 –N α site exhibits a small but reproducible deviation from this monotonic trend near the eutectic. NOESY spectra at similar compositions reveal enhanced cross-relaxation between EG and choline protons, consistent with increased headgroup–solvent contact density rather than a wholesale structural rearrangement. Overall, our multitechnique study demonstrates that EG–ChCl dynamics are predominantly viscosity-dominated, with the eutectic region acting as a subtle dynamical crossover where specific choline segments become maximally coupled to the hydrogen-bond network. These insights refine the structure–dynamics picture of choline-chloride DESs and provide practical guidance for tuning composition in electrochemical, separation, and catalytic applications.

diffusion↗
DOE OSTI2025

Atomic-Scale Insights into Carbon Dissolution in α-, γ-, and θ-Al 2 O 3 : Phase-dependent Transport Dynamics from First-Principles Calculations

α-Al 2 O 3 exhibits superior carburizing corrosion resistance compared to metastable γ-Al 2 O 3 and θ-Al 2 O 3 phases in high-temperature CO 2 environments, yet its atomic-scale origins remain unclear. Using first-principles density functional theory, we systematically investigate carbon dissolution and diffusion in α-Al 2 O 3 , γ-Al 2 O 3 , and θ-Al 2 O 3 , including the effects of oxygen (O) and aluminum (Al) vacancies. Our results show that α-Al 2 O 3 consistently exhibits higher carbon solution enthalpies than γ-Al 2 O 3 and θ-Al 2 O 3 in both pristine and defective structures, indicating lower intrinsic carbon solubility in α-Al 2 O 3 . Vacancies significantly enhance carbon incorporation: O vacancies reduce solution enthalpy, while Al vacancies further amplify this effect, with a strong preference for carbon at Al vacancy sites. Carbon diffusion barriers are also highest in α-Al 2 O 3 , reflecting slower carbon mobility. Al vacancies increase diffusion barriers across all phases, while O vacancies raise barriers in α- and γ-Al 2 O 3 but slightly lower them in θ-Al 2 O 3 . These results reveal a dual mechanism behind the carburizing resistance of α-Al 2 O 3 : reduced carbon solubility and elevated diffusion barriers. Furthermore, this work provides atomic-scale insights to guide the design of alumina-based materials with improved carburizing resistance through phase selection and defect engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2025

DNA-PAINT Imaging with Hydrogel Imprinting and Clearing

Hydrogel-embedding is a versatile technique in fluorescence microscopy, offering stabilization, optical clearing, and the physical expansion of biological specimens. DNA-PAINT is a super-resolution microscopy approach based on the diffusion and transient binding of fluorescently labeled oligos, but its feasibility in hydrogels has not yet been explored. In this study, we demonstrate that polyacrylamide hydrogels support sufficient diffusion for effective DNA-PAINT imaging. Using acrydite-anchored oligonucleotides imprinted from patterned DNA origami nanostructures and microtubule filaments in fixed cells, we find that hydrogel embedding preserves docking strand positioning at the nanoscale. Sample clearing via protease treatment had minor structural effects on the microtubule structure and enhanced diffusion and accessibility to hydrogel-imprinted docking strands. Our work demonstrates promising potential for diffusion and binding-based fluorescence imaging applications in hydrogel-embedded samples.

DNA origami↗
DOE OSTI2025

In-Plane Anisotropy of Charge Density Wave Fluctuations in 1⁢𝑇−TiSe 2

We report measurements of anisotropic triple-𝑞 charge density wave (CDW) fluctuations in the transition metal dichalcogenide 1⁢𝑇−TiSe 2 over a large volume of reciprocal space with x-ray diffuse scattering. Above the transition temperature, 𝑇 CDW , the in-plane diffuse scattering is marked by ellipses which reveal that the in-plane fluctuations are anisotropic. In addition, the out-of-plane diffuse scattering is characterized by rodlike structures which indicate that the CDW fluctuations in neighboring layers are largely decoupled. Furthermore, our analysis of the diffuse scattering line shapes and orientations suggests that the three charge density wave components contain independent phase fluctuations with a hierarchy of length scales, leading to intricate fluctuation patterns that go beyond the conventional 2D-to-3D crossover picture.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

An Asymptotic Preserving Discontinuous Galerkin Method for a Linear Boltzmann Semiconductor Model

A key property of the linear Boltzmann semiconductor model is that as the collision frequency tends to infinity, the phase space density $f$ = $f$ ($x, v, t$) converges to an isotropic function $M (v)$$ρ$$(x, t)$, called the drift-diffusion limit, where $M$ is a Maxwellian and the physical density $ρ$ satisfies a second-order parabolic PDE known as the drift-diffusion equation. Numerical approximations that mirror this property are said to be asymptotic preserving. In this paper we build a discontinuous Galerkin method to the semiconductor model, and we show this scheme is both uniformly stable in $ε$, where 1/$ε$ is the scale of the collision frequency, and asymptotic preserving. Here in particular, we discuss what properties the discrete Maxwellian must satisfy in order for the schemes to converge in $ε$ to an accurate $h$-approximation of the drift-diffusion limit. Discrete versions of the drift-diffusion equation and error estimates in several norms with respect to $ε$ and the spacial resolution are also included.

97 MATHEMATICS AND COMPUTING↗