Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion current”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stochastic fluctuations in relativistic fluids: Causality, stability, and the information current

We develop a general formalism for introducing stochastic fluctuations around thermodynamic equilibrium which takes into account, for the first time, recent developments in the causality and stability properties of relativistic hydrodynamic theories. The method is valid for any covariantly stable theory of relativistic viscous fluid dynamics derived from a covariant maximum entropy principle. We illustrate the formalism with some applications, showing how it could be used to consistently introduce fluctuations in a model of relativistic heat diffusion and in conformally invariant Israel-Stewart theory in a general hydrodynamic frame. Furthermore, the latter example is used to study the hydrodynamic frame dependence of the symmetric two-point function of fluctuations of the energy-momentum tensor.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Selenium Migration and Local Structures in Cu‐Doped CdSeTe Solar Cells after Aging

Selenium grading plays a critical role in state-of-the-art Cadmium Telluride photovoltaic cells by enhancing long-wavelength absorption and extending minority carrier lifetimes —key to enabling the current performance record of 23.08%. However, very little is understood about selenium motion. In this study, a comprehensive, multimodal, and multiscale approach is employed to investigate Se migration and local structural changes in copper (Cu)-doped CdSeTe solar cells subjected to accelerated stress. X-ray fluorescence (XRF) microscopy shows unexpected levels of Se diffusion after 500 h under heat (75°C) and light (0.8 suns, 80 mW/cm 2 ), suggesting the coexistence of fast and slow diffusion channels even at low temperatures, with unexpectedly low activation energies (<0.85 eV). X-ray Absorption Near Edge Structure (XANES) analysis indicates a preferential migration of Se atoms to anionic lattice sites and a reduction in Se-Cl co-passivation at Te-terminated dislocation cores. Furthermore, these findings point to a reconfiguration of Se local environments and highlight the potential role of extended structural defects in enabling Se transport at low temperatures. Additionally, XANES results suggest that the presence of metallic Cu across the absorber layer may contribute to back-contact degradation and reduced hole density in both fresh and aged devices.

14 SOLAR ENERGY↗
DOE OSTI2025

Creating Favorable Pt/Co Interfaces via a Two‐Step Approach for Constructing Highly Durable PtCo Intermetallic Fuel Cell Catalysts

Structurally ordered PtCo intermetallics are one of the most promising oxygen-reduction catalysts in proton exchange membrane fuel cells (PEMFCs) due to their intrinsically improved catalytic activity and stability relative to PtCo solid-solution alloys. However, increasing the heating temperature to achieve a desirable high degree of ordering results in severe particle agglomeration and low mass activity and stability. Herein, a two-step synthesis approach is developed to create an L1 2 -Pt 3 Co intermetallic structure with an increased ordering degree and well-dispersed ultrafine particles. The first step of the synthesis yields ultrafine Pt nanoparticles that are well-dispersed on the ZIF-8-derived carbon support. The second adsorption step enables us to fine-tune the Pt and Co interfaces, assisted by optimal amino acids, to establish a favorable Co-rich environment around fine Pt nanoparticles, facilitating Co diffusion into the Pt crystalline under mild thermal conditions (<800 °C). In conclusion, this two-step ordered L1 2 -Pt 3 Co catalyst is systematically evaluated using membrane electrode assemblies under heavy-duty vehicle (HDV) conditions and demonstrated exceptional performance and durability, retaining 1.35 A cm -2 only a 7% loss in current density at 0.7 V after an extensive accelerated stress test of 150,000 voltage cycles.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Lithium‐Containing Hybrid SEI Layer Enabling High Mass Loading and Anode‐Less Sodium Metal Batteries

Abstract The continuous rupturing and rebuilding of unstable solid electrolyte interphase (SEI) layer during cycling would block Na + diffusion and induce Na dendrite formation, ultimately limiting the practical application of high‐energy‐density sodium metal batteries. Herein, a hybrid SEI layer containing Li‐species is dexterously constructed on the surface of sodium metal anode. Li‐containing inorganic components (Li 3 N, LiF, and Li 2 CO 3 ) are introduced to stabilize the Na/electrolyte interface and enhance the mechanical and diffusion kinetic properties of the SEI layer, which can reduce the side reactions and gas generation, regulate Na + flux during cycling and promote rapid Na + migration for uniform dendrite‐free Na deposition. As a result, the constructed Na symmetric cells achieve low overpotential and long cycle life of 5900, 1800, and 500 h at current densities of 3, 10, and 30 mA cm −2 , respectively. Furthermore, the full cells paired with the Na₃V₂(PO₄)₃ cathode demonstrate high specific capacity and excellent cycle stability, even at an ultra‐high cathode loading of 39.3 mg cm −2 and a low N/P ratio (negative/positive electrode capacity ratio of 1.21).

Xia, Li [State Key Laboratory of Physical Chemistr↗
DOE OSTI2021

Fictitious phase separation in Li layered oxides driven by electro-autocatalysis

Layered oxides widely used as lithium-ion battery electrodes are designed to be cycled under conditions that avoid phase transitions. Although the desired single-phase composition ranges are well established near equilibrium, operando diffraction studies on many-particle porous electrodes have suggested phase separation during delithiation. Notably, the separation is not always observed, and never during lithiation. These anomalies have been attributed to irreversible processes during the first delithiation or reversible concentration-dependent diffusion. However, these explanations are not consistent with all experimental observations such as rate and path dependencies and particle-by-particle lithium concentration changes. Here, we show that the apparent phase separation is a dynamical artefact occurring in a many-particle system driven by autocatalytic electrochemical reactions, that is, an interfacial exchange current that increases with the extent of delithiation. We experimentally validate this population-dynamics model using the single-phase material L ix (Ni 1/3 Mn 1/3 Co 1/3 )O 2 (0.5 < x < 1) and demonstrate generality with other transition-metal compositions. Operando diffraction and nanoscale oxidation-state mapping unambiguously prove that this fictitious phase separation is a repeatable non-equilibrium effect. We quantitatively confirm the theory with multiple-datastream-driven model extraction. More generally, our study experimentally demonstrates the control of ensemble stability by electro-autocatalysis, highlighting the importance of population dynamics in battery electrodes (even non-phase-separating ones).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Protective coatings with nanoscale multilayer architecture: current state and main trends

Binary systems (layers) and multicomponent layers with bilayer thicknesses varying from microns to a few nanometers are considered. A strong dependence of the microstructure formation on the deposition conditions (the magnitude of the bias potential, the flow of inert and reactive gases, the ion beam power (current density) supplied to the target) is found. General regularities in the change in the functional properties and structural state of the coatings are formulated. The conditions for the formation of a columnar structure and growth texture are noted; various mechanisms and their interpretation are proposed. Finally, the relationship among the size of nanograins, diffusion at the interphase boundary, and effects determining the physico-mechanical, tribological characteristics, as well as resistance to corrosion and oxidation, is shown.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Jet wake from linearized hydrodynamics

We explore how to improve the hybrid model description of the particles originating from the wake that a jet produced in a heavy ion collision leaves in the droplet of quark-gluon plasma (QGP) through which it propagates, using linearized hydrodynamics on a background Bjorken flow. Jet energy and momentum loss described by the hybrid model become currents sourcing linearized hydrodynamics. By solving the linearized hydrodynamic equations numerically, we investigate the development of the wake in the dynamically evolving droplet of QGP, study the effect of viscosity, scrutinize energy-momentum conservation, and check the validity of the linear approximation. We find that linearized hydrodynamics works better in the viscous case because diffusive modes damp the energy-momentum perturbation produced by the jet. We calculate the distribution of particles produced from the jet wake by using the Cooper-Frye prescription and find that both the transverse momentum spectrum and the distribution of particles in azimuthal angle are similar in shape in linearized hydrodynamics and in the hybrid model. Their normalizations are different because the momentum-rapidity distribution in the linearized hydrodynamics analysis is more spread out, due to sound modes. Since the Bjorken flow has no transverse expansion, we explore the effect of transverse flow by using local boosts to add it into the Cooper-Frye formula. After including the effects of transverse flow in this way, the transverse momentum spectrum becomes harder: more particles with transverse momenta bigger than 2 GeV are produced than in the hybrid model. Although we defer implementing this analysis in a jet Monte Carlo, as would be needed to make quantitative comparisons to data, we gain a qualitative sense of how the jet wake may modify jet observables by computing proxies for two example observables: the lost energy recovered in a cone of varying open angle, and the fragmentation function. We find that linearized hydrodynamics with transverse flow effects added improves the description of the jet wake in the hybrid model in just the way that comparison to data indicates is needed. Our study illuminates a path to improving the description of the wake in the hybrid model, highlighting the need to take into account the effects of both transverse flow and the broadening of the energy-momentum perturbation in spacetime rapidity on particle production.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Digitization modeling of a CGEM detector based on Garfield++ simulation

Abstract Background The Cylindrical Gas Electron Multiplier Inner Tracker (CGEM-IT) is one of the candidates for the replacement of the current BESIII inner drift chamber, which has developed serious aging effects due to beam-related backgrounds. Purpose To build a full digitization model, it is necessary to understand the electron drift behavior and signal induction. Methods A simulation study based on Garfield++ and ANSYS software has been performed. Results Parameters describing the Lorentz angle, diffusion effect, drift time, multiplication, and induction are obtained from the Garfield++ simulation. Based on these results, a preliminary digitization model is implemented in the BESIII offline software system.

Zhao, Jing-Yi (ORCID:0000000220287286)↗
DOE OSTI2025

Lowering the barrier to access information-rich transient kinetic data for machine learning methods

Transient kinetic data contain a wealth of information about intrinsic features of a catalyst as well as the reaction mechanism. Currently, high volume transient data is underutilized, and data science methods could both increase the value of information that can be extracted from this data, integrate experimental with theoretical data sources, and accelerate the pace of catalyst technology advancement. Transient kinetic characterizations with simple probe molecules exhibiting reversible adsorption, irreversible adsorption and bulk-surface diffusion are presented as training components for similar experiments with more complex surface reactions. In conclusion, by increasing the availability and accessibility of transient kinetic data through details of its structure and acquisition, we aim to decrease the barrier for data scientists to apply machine learning methods to this valuable data source.

Catalysis↗
DOE OSTI2026

Multiscale Characterization of Electrode-Induced Degradation in Perovskite Solar Cells

The stability of metal-halide-perovskite (MHP) solar cells must be understood and improved for the commercial viability of MHP technologies. Here, we apply multiscale characterization methods to study degradation modes, specifically electrode corrosion, for p-i-n MHP partial device stacks and full devices that are stored in the dark under an inert atmosphere. Our multiscale characterization approaches include full-device electro-optical performance using current-voltage (JV) curves and spatial imaging with electroluminescence (EL) and photoluminescence (PL). We further correlate interface properties using cross-sectional Kelvin probe force microscopy, which maps the nanoscale electric field properties, and electron microscopy, which demonstrates structural and chemical features. Devices stored as a full device stack degrade primarily by metal (Ag) electrode diffusion into the absorber, with formation of AgI byproducts and Ag accumulation near the indium tin oxide (ITO) contact. This causes decomposition of the perovskite absorber domains, loss of the potential drop at the electron transport layer (ETL)/perovskite interface near the metal contact, and increased equivalent resistance at the perovskite/hole transport layer (HTL) interface near the ITO contact. The devices stored without metal show a different degradation pathway dominated by corrosion of the ITO, creating voids at the ITO electrode surface with diffusion of In and Sn into the absorber. We conclude that metal electrode-induced degradation is the most severe degradation pathway under dark storage, but that ITO corrosion and absorber instability must also be mitigated. We further demonstrate mitigation of these degradation pathways by changes to the device stack, including a SnO x blocking layer at the ETL side and replacing ITO with FTO at the HTL side. These results provide a useful demonstration of specific dark degradation pathways at each electrode interface, as well as a unique multiscale example that links degradation of chemical, structural, and electrical interface properties to the full-device electro-optical characteristics.

14 SOLAR ENERGY↗
DOE OSTI2021

Computational general relativistic force-free electrodynamics: II. Characterization of numerical diffusivity

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe’s most extreme flows of energy. The discretization of Maxwell’s equations – needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling – introduces numerical diffusion. It acts as a source of dissipation independent of the system’s physical constituents. Understanding the numerical diffusion of scientific codes is the key to classifying their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field (j ∥ ), in combination with the breakdown of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g., the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of j ∥ .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Diffusion Synthetic Acceleration for Heterogeneous Domains, Compatible with Voids

A standard approach to solving the S N transport equations is to use source iteration with diffusion synthetic acceleration (DSA). Although this approach is widely used and effective on many problems, there remain some practical issues with DSA preconditioning, particularly on highly heterogeneous domains. For large-scale parallel simulation, it is critical that both (a) preconditioned source iteration converges rapidly and (b) the action of the DSA preconditioner can be applied using fast, scalable solvers, such as algebraic multigrid (AMG). For heterogeneous domains, these two interests can be at odds. In particular, there exist DSA diffusion discretizations that can be solved rapidly using AMG, but they do not always yield robust/fast convergence of the larger source iteration. Conversely, there exist robust DSA discretizations where source iteration converges rapidly on difficult heterogeneous problems, but fast parallel solvers like AMG tend to struggle applying the action of such operators. Moreover, very few current methods for the solution of deterministic transport are compatible with voids. This paper develops a new heterogeneous DSA preconditioner based on only preconditioning the optically thick subdomains. The resulting method proves robust on a variety of heterogeneous transport problems, including a linearized hohlraum mesh related to inertial confinement fusion. Moreover, the action of the preconditioner is easily computed using O(1) AMG iterations, convergence of the transport iteration typically requires 2 to 5× fewer iterations than current state-of-the-art “full” DSA, and the proposed method is trivially compatible with voids. On the hohlraum problem, rapid convergence is obtained by preconditioning less than 3% of the mesh elements with five to ten AMG iterations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Editors’ Choice—Flooded by Success: On the Role of Electrode Wettability in CO 2 Electrolyzers that Generate Liquid Products

Economic operation of carbon dioxide (CO 2 ) electrolyzers generating liquid products will likely require high reactant conversions and product concentrations, conditions anticipated to challenge existing gas diffusion electrodes (GDEs). Notably, electrode wettability will increase as lower surface tension products (e.g., formic acid, alcohols) are introduced into electrolyte streams, potentially leading to flooding. To understand the hydraulically stable operating envelopes in mixed aqueous-organic liquid domains, we connect intrinsic electrode wettability descriptors to operating parameters such as electrolyte flow rate and current. We first measure contact angles of water-organic product dilutions on polytetrafluoroethylene (PTFE) and graphite surfaces as planar analogues for GDE components. We then use material balances around the reactive gas-liquid-solid interface to calculate product mass fractions as functions of water sweep rate and current. Product composition maps visualize the extent to which changes in cell performance influence capillary pressure, a determinant of GDE saturation. Analyses suggest that formic acid mixtures pose little risk for GDE flooding across a wide range of conditions, but effluents containing <30% alcohol by mass may cause flooding. This study reveals opportunities to integrate microstructural features and oleophobic surface treatments into GDEs to repel aqueous-organic mixtures and expand the window of stable operating conditions.

25 ENERGY STORAGE↗
DOE OSTI2024

Understanding and Tailoring Diffusion and Co-Adsorption Inside the Confined Pores of Metal-Organic Frameworks (Final Scientific/Technical Report for Award DE-SC0019902)

The aim of this program was to gain a fundamental understanding of the behavior of various guest molecules in nano-confined environments, such as metal organic frameworks (MOFs), using a combination of novel synthesis, ab initio modeling, and in situ characterization. Through this project, we developed a concise understanding of the mechanisms that control adsorption/desorption of gaseous molecules and their mixtures, leading to design/synthesis guidelines for MOFs with desired functionality. We further developed methods to disentangle kinetic from thermodynamic effects during adsorption, as well as to characterize the interactions at play. In the first funding cycle, the focus was on the unambiguously characterization of co-adsorption and diffusion of gasses/vapors and their mixtures. In the second funding cycle, the focus was on characterizing the effects of the nano-confinement on the kinetics and thermodynamics of adsorption processes inside MOFs, again with an emphasis on mixtures of gasses and vapors. The nano-confinement can tip the thermodynamic vs. kinetic balance, and current understanding and theory based on single-component analysis can lead to incorrect predictions for mixtures. This is of particular interest in real-world applications, where gasses/vapors are typically mixed, contain impurities, or are often exposed to humid conditions. Our main findings were: (i) within confined environments the adsorption behavior of mixed gasses/vapors can be drastically different from the “sum” of the corresponding single phases; (ii) co-adsorption is often competitive and detrimental to performance, but it can also be cooperative and beneficial; (iii) in some co-adsorbed gasses/vapors, molecules that are strongly bound in the single-component phase can be replaced by molecules that are nominally weaker bound (molecular exchange) due to guest-guest interactions that lower the kinetic barriers and favor the final adsorption state; (iv) kinetic and thermodynamic effects can be precisely controlled through pore-size engineering and synthesis; and, (v) kinetic effects can be identified and disentangled from thermodynamic effects during adsorption through a series of sequential and simultaneous gas loading measurements. The short-term goal of this program was the controlling and understanding of common MOF systems in real-world situations where gasses/vapors are mixed, which will have an important impact on industrial processes and applications from gas storage and sequestration to catalysis and sensors. The long-term goals include the development of theoretical and experimental methods for gaining a fundamental understanding of adsorption/reaction processes within MOFs, as well as new guidelines for synthesizing MOFs with tailored physical and chemical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

I-V curve analysis of single-nanoparticle Al-doped SrTiO 3

Strontium Titanate (SrTiO 3 or STO) is an oxide-semiconductor material which shows promise for photocatalytic water splitting, but remains poorly understood. To better understand its carrier dynamics, STO nanoparticles deposited onto an Au substrate were placed inside of a Scanning Electron Microscope (SEM) and then contacted with an Au probe. Current-voltage (I-V) measurements were taken with and without electron beam illumination. Fitting to the curves, we estimate the doping density ($N_D \approx 10^{22.1 \pm 1.2}$ cm$^{-3}$) and built-in potential ($\phi_{bi} \approx 0.84 \pm 0.45$ eV). We then used these I-V measurements to estimate the diffusion time ($\tau$) of our samples and observed an exponential dependence of $\tau$ on sample size in the range of 350-550 nm.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantitative Prediction of Perovskite Stability Using Accelerated Testing and Machine Learning

PV technologies based on hybrid perovskites offer the potential for reducing solar cell costs, but they are particularly vulnerable to degradation by environmental factors such as moisture, oxygen, and illumination. Commercialization will require not only stable materials and device architectures but also accelerated testing protocols and models that can predict degradation from the accelerated testing data. Here, we report results from in situ photoluminescence (PL) and photoconductivity (PC) measurements during perovskite degradation with simultaneous optical transmittance (Tr) measurements or reflected dark field (DF) imaging. From PL, PC, Tr, and DF, we determine (respectively) the steady-state quasi-Fermi level splitting, the mean effective carrier diffusion length, the extent of conversion of perovskite to higher bandgap degradation products, and the extent of scattering from domains with different orientation or composition, all as a function of time during degradation. Simultaneous measurement of PL-PC-Tr or PL-PC-DF on perovskite absorbers in an environmental chamber over a wide range of humidity, oxygen, temperature, and illumination levels yields a rich data set. We use machine learning to develop a model that accurately predicts (within 10%) the time for the diffusion length to decrease to 85% of its initial value. The model takes the environmental conditions and the first few measurements from the PL-PC-Tr or PL-PC-DF experiment as input. Thus, the model provides a framework to interpret the results of accelerated testing of absorber materials. One of the dominant features in the model of degradation for CH 3 NH 3 PbI 3 is the initial rate at which transmittance increase. For devices with opaque contacts, the pixel-averaged rate of change of the intensity in dark-field images can be used in place of transmittance. Results are also presented on full PV devices with in-situ current-voltage (JV) measurements. Here, the data from simultaneous PL-DF-JV under environmental stresses reveal intimate connection between degradation and shunts and provide a framework to extend accelerated testing to devices. Further, the presentation reveals underlying universal behavior in degradation pathways over a broad range of environmental stresses and perovskite compositions.

14 SOLAR ENERGY↗
DOE OSTI2025

Ensemble‐Based, Large‐Eddy Reconstruction of Wind Turbine Inflow in a Near‐Stationary Atmospheric Boundary Layer Through Generative Artificial Intelligence

ABSTRACT To validate the second‐by‐second dynamics of turbines in field experiments, it is necessary to accurately reconstruct the winds going into the turbine. Current time‐resolved inflow reconstruction techniques estimate wind behavior in unobserved regions using relatively simple spectral‐based models of the atmosphere. Here, we develop a technique for time‐resolved inflow reconstruction that is rooted in a large‐eddy simulation model of the atmosphere. Our “large‐eddy reconstruction” technique blends observations and atmospheric model information through a diffusion model machine learning algorithm, allowing us to generate probabilistic ensembles of reconstructions for a single 10‐min observational period. Our generated inflows can be used directly by aeroelastic codes or as inflow boundary conditions in a large‐eddy simulation. We verify the second‐by‐second reconstruction capability of our technique in three synthetic field campaigns, finding positive Pearson correlation coefficient values () between ground‐truth and reconstructed streamwise velocity, as well as smaller positive correlation coefficient values for unobserved fields (spanwise velocity, vertical velocity, and temperature). We validate our technique in three real‐world case studies by driving large‐eddy simulations with reconstructed inflows and comparing to independent inflow measurements. The reconstructions are visually similar to measurements, follow desired power spectra properties, and track second‐by‐second behavior ().

17 WIND ENERGY↗
DOE OSTI2021

Collisional effects on resonant particles in quasilinear theory

A careful examination of the effects of collisions on resonant wave-particle interactions leads to an alternate interpretation and deeper understanding of the quasilinear operator originally formulated by Kennel and Engelmann (Phys. Fluids vol. 9, 1966, pp. 2377- 2388) for collisionless, magnetized plasmas, and widely used to model radio frequency heating and current drive. The resonant and nearly resonant particles are particularly sensitive to collisions that pitch angle scatter them out of and into resonance. As a result, the resonant particle-wave interactions occur in the center of a narrow collisional boundary when the collision frequency nu is very small compared to the wave frequency omega. The diffusive nature of the pitch angle scattering combined with the wave-particle resonance condition enhances the collision frequency by (omega/nu)2/3 >>1, resulting in an effective resonant particle collision time of tau_int ~ (nu /omega)2/3 nu <<1/ nu . A rigorous collisional boundary layer analysis generalizes the standard quasilinear operator to a form that is fully consistent with Kennel-Englemann, but allows replacing the delta function appearing in the diffusivity with a simple integral (having the appropriate delta function limit) retaining the new physics associated with the narrow boundary layer, while preserving the entropy production principle. The limitations of the collisional boundary layer treatment are also estimated, and indicate that substantial departures from Maxwellian are not permitted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗