Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Comparison”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

NaCl aggregation in water at elevated temperatures and pressures: Comparison of classical force fields

The properties of water vary dramatically with temperature and density. This can be exploited to control its effectiveness as a solvent. Thus, supercritical water is of keen interest as solvent in many extraction processes. The low solubility of salts in lower density supercritical water has even been suggested as a means of desalination. The high temperatures and pressures required to reach supercritical conditions can present experimental challenges during collection of required physical property and phase equilibria data, especially in salt-containing systems. Molecular simulations have the potential to be a valuable tool for examining the behavior of solvated ions at these high temperatures and pressures. However, the accuracy of classical force fields under these conditions is unclear. We have, therefore, undertaken a parametric study of NaCl in water, comparing several salt and water models at 200 bar–600 bar and 450 K–750 K for a range of salt concentrations. We report a comparison of structural properties including ion aggregation, hydrogen bonding, density, and static dielectric constants. All of the force fields qualitatively reproduce the trends in the liquid phase density. An increase in ion aggregation with decreasing density holds true for all of the force fields. The propensity to aggregate is primarily determined by the salt force field rather than the water force field. This coincides with a decrease in the water static dielectric constant and reduced charge screening. While a decrease in the static dielectric constant with increasing NaCl concentration is consistent across all model combinations, the salt force fields that exhibit more ionic aggregation yield a slightly smaller dielectric decrement.

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

Direct comparison of wire, foil, and hybrid X-pinches on a 200 kA, 150 ns current driver

Wire X-pinches (WXPs) have been studied comprehensively as fast (~1 ns pulse width), small (~1 μm) x-ray sources, created by twisting two or more fine wires into an “X” to produce a localized region of extreme magnetic pressure at the cross-point. Recently, two alternatives to the traditional WXP have arisen: the hybrid X-pinch (HXP), composed of two conical electrodes bridged by a thin wire or capillary, and the laser-cut foil X-pinch (LCXP), cut from a thin foil using a laser. We present a comparison of copper wire, hybrid, and laser-cut foil X-pinches on a single experimental platform: UC San Diego’s ~200 kA, 150 ns rise time GenASIS driver. All configurations produced 1–2 ns pulse width, ≤5 μm soft x-ray (Cu L-shell, ~1 keV) sources (resolutions diagnostically limited) with comparable fluxes. WXP results varied with linear mass and wire count, but consistently showed separate pinch and electron-beam-driven sources. LCXPs produced the brightest (~1 MW), smallest (≤5 μm) Cu K-shell sources, and spectroscopic data showed both H-like Cu K α lines indicative of source temperatures ≥2 keV, and cold K α (~8050 eV) characteristic of electron beam generated sources, which were not separately resolved on other diagnostics (within 1–2 ns and ≤200 μm). HXPs produced minimal K-shell emission and reliably single, bright, and small L-shell sources after modifications to shape the early current pulse through them. Benefits and drawbacks for each configuration are discussed to provide potential X-pinch users with the information required to choose the configuration best suited to their needs.

42 ENGINEERING↗
DOE OSTI2021

Comparison of the dielectric and magnetocaloric properties of bulk and film of GdFe 0.5 Cr 0.5 O 3

Here, we report a comparison of the magnetic, magnetocaloric, and dielectric properties of 50% iron substituted GdCrO 3 (GdFe 0.5 Cr 0.5 O 3 ) bulk pellet and 960 nm thick film of GdFe 0.5 Cr 0.5 O 3 (GFCO). The 960 nm film was synthesized on a platinized-silicon substrate by chem. soln. deposition and spin-coating methods. The X-ray diffraction scans of the bulk sample and the film as well as the morphology of the film as examined by the field-emission scanning electron microscope indicate phase-pure and polycrystalline nature of these samples. XPS was used to deermine the valence states of Gd, Fe, and Cr. The temp. dependence of the dielectric const. from 225 to 700 K shows peaks at T C = 525 K for the bulk and ≈450 K for the film due to ferroelectric to paraelectric transitions, since elec. polarization vs elec. field hysteresis loops are observed. at room temp. The dielectric studies in the bulk GFCO for T > T C indicate a relaxor-like behavior. The measurements of the magnetization (M) of the samples as a function of temp. (5-350 K) and magnetic field (H) up to 7 T (=70 kOe) depict hysteresis behavior at low temps. due to the canted antiferromagnetic order of Fe 3+ /Cr 3+ below the Néel temp. of ≈275 K. The M vs H isotherms at various temps. are used to det. and compare the magnetic entropy change (-ΔS) and relative cooling power (RCP) of the two samples, yielding (-ΔS) = 30.7 J/kg K (18.8 J/kg K) and RCP = 566.5 J/kg (375 J/kg) for the bulk (960 nm film) samples of GFCO at 7 K and 7 T, resp. The plot of RCP vs T shows that magnetic cooling for this system is most effective for T < 30 K. Comparatively smaller magnitudes of (-ΔS) and RCP for the film vis-a-vis the bulk sample of GFCO scale with its reduced magnetization. This suggests that further improvements in the quality of the films are needed to improve their magnetization and hence their magnetocaloric properties, possibly making them useful for on-chip cooling in miniaturized devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Code-to-code comparison and validation of the radiation-hydrodynamics capabilities of the FLASH code using a laboratory astrophysical jet

The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and behavior of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities in the publicly available FLASH radiation-hydrodynamics code, we present detailed code-to-code comparisons between FLASH and the HYDRA code developed at Lawrence Livermore National Laboratory using previously published HYDRA simulations from Grava et al. [Phys. Rev. E 78, 016403 (2008)]. That study describes a laser experiment that produced a jet-like feature that the authors compare to astrophysical jets. Importantly, the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment included detailed x-ray interferometric measurements of electron number densities and a time-integrated measurement of the soft x-ray spectrum. Despite markedly different methods for treating the computational mesh, and different equations of state and opacity models, the FLASH results resemble the results from HYDRA and, most importantly, the experimental measurements of electron density. Having validated the FLASH code in this way, we use the code to further investigate and understand the formation of the jet seen in the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment and discuss its relation to the Wan et al. [Phys. Rev. E 55, 6293 (1997)] experiment at the NOVA laser.

Orban, Chris (ORCID:0000000312004538)↗
DOE OSTI2022

System-Level Comparison of Sodium and Salt Systems in Support of the Gen3 Liquids Pathway

This paper presents the comparison of the annual performance and the Levelised Cost of Energy (LCOE) of two high- temperature solar power tower configurations using a 565MWth chloride salt receiver and a 540 MWth liquid sodium receiver. Each power tower consists of a sodium or salt receiver, a 175 m high concrete tower, and a two-tank chloride salt system for thermal energy storage (TES). The power plants were simulated using SolarTherm to estimate their energy per year, capacity factor, and LCOE. The simulation was verified against the System Advisor Model (SAM), showing that the LCOE calculation differences are within ±1.0%. Results showed that the sodium-receiver system yields 69.60 USD/MWh real LCOE, versus 78.55 USD/MWh for the chloride salt system. Results also showed that sodium wins here due to a smaller and more efficient receiver, enabled by higher flux limits and reduced field size, leading to less mirror reflective area and site improvement cost.

capacity factor↗
DOE OSTI2024

Large-database cross-verification and validation of tokamak transport models using baselines for comparison

State-of-the-art 1D transport solvers ASTRA and TRANSP are verified, then validated across a large database of semi-randomly selected, time-dependent DIII-D discharges. Various empirical models are provided as baselines to contextualize the validation figures of merit using statistical hypothesis tests. For predicting plasma temperature profiles, no statistically significant advantage is found for the ASTRA and TRANSP simulators over a baseline empirical (two-parameter) model. For predicting stored energy, a significant advantage is found for the simulators over a baseline empirical model based on confinement time scaling. Uncertainty in the results due to diagnostic and profile fitting uncertainties is approximated and determined to be insignificant due in part to the large quantity of discharges employed in the study. Advantages are discussed for validation methodologies like this one that employ (1) large databases and (2) baselines for comparison that are specific to the intended use-case of the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Comparison of tungsten versus molybdenum for double shell capsules using machine learning design optimization

Double shell targets are an alternative ignition platform for inertial confinement fusion. One design consideration for double shell targets is the choice of inner shell material to help trap radiation emitted by the hot fuel to aid ignition. Materials such as molybdenum and tungsten are of interest for the inner shell layer of the targets. While molybdenum has a lower density that could inhibit instability growth and allow for radiography and code benchmarking, tungsten has a higher density that could provide more compression and confinement. These tradeoffs have been explored using optimized designs for each material. Our previous work [Vazirani et al., “Coupling 1D xRAGE simulations with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 28, 122709 (2021); Vazirani et al., “Coupling multi-fidelity xRAGE with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 30, 062704 (2023); and Vazirani et al., “Bayesian batch optimization for molybdenum versus tungsten inertial confinement fusion double shell target design,” Stat. Anal. Data Min. 17, e11698 (2024)] resulted in a multi-fidelity Bayesian optimization framework to find yield-optimized double shell target geometries. By leveraging simulations of varying fidelities (one-dimensional and two-dimensional) to inform one another, the multi-fidelity optimization was able to optimize a design in the highest fidelity with significantly fewer simulations than would be used in a systematic parameter scan. In this work, we apply the multi-fidelity Bayesian optimization to explore the optimized designs of double shell targets with molybdenum and tungsten inner shells as well as the physics producing the high performing implosions. A physics exploration of all the simulations used in this study shows trends in designs that contribute to high yields, ion temperatures, and fuel areal densities. Comparison of molybdenum and tungsten simulations shows that they can produce similar implosion conditions with different geometries, which would be important to study in experiments. Graded density layers produce varying performances with the two materials but continue to be of interest for future studies along with studies of doped inner shell materials and applied surface roughness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Direct comparison of gyrokinetic and fluid scrape-off layer simulations

Typically, fluid simulations are used for tokamak divertor design. However, fluid models are only valid if the scrape-off layer (SOL) is highly collisional. This assumption is valid in many present-day experiments but is questionable in the upstream SOL of some high-power scenarios envisioned for burning plasmas and fusion pilot plants. This paper reports on comparisons between fluid and kinetic simulations of the SOL for upstream parameters and geometry representative of the Spherical Tokamak for Energy Production fusion pilot plant. The SOLPS-ITER (fluid) and Gkeyll (gyrokinetic) codes are operated in a two-dimensional axisymmetric mode, which replaces turbulence with ad-hoc diffusivities. In kinetic simulations, we observe that the ions in the upstream SOL experience significant mirror trapping. This substantially increases the upstream temperature and has important implications for impurity dynamics. We show that the mirror force, which is excluded in SOLPS’s fluid equations, enhances the electrostatic potential drop along the field line in the SOL. We also show that the assumption of equal main ion and impurity temperatures, which is made in commonly used fluid codes, is invalid for the regimes explored here. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions [Kotschenreuther et al., in 29th IAEA 29 Fusion Energy Conference (IAEA, London, UK, 2023)]. This effect can be dramatic, reducing the midplane impurity density by orders of magnitude. These results indicate that in lower collisionality SOL’s the tolerable downstream impurity densities may be higher than would be predicted by fluid simulations, allowing for higher radiated power while avoiding unacceptable core contamination. Our results highlight the importance of kinetic simulations for divertor design and optimization for fusion pilot plants.

Computational fluid dynamics↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

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

A comparison of the dielectric function for p-point water models

We compute the longitudinal dielectric function ɛl(k) for p-point water models, p = 3 (SPC/E, TIP3P), p = 4 (TIP4P-EW, OPC), and p = 5 (TIP5P-E) from the charge–charge fluctuation function SZZ(k), obtained from two methods: average of the charge in k-space (method:sum) and from the pair distribution function in Fourier space (method:integral). The latter requires a continuation to small k-values, which we thoroughly discuss. We conclude with a detailed comparison of the longitudinal dielectric function for the different p-point models and its main characteristics: small k and large k limits, position of the poles, and its (complex) zeros.

Shevy, Noah [Ames National Laboratory, Department ↗
DOE OSTI2021

Comparison of Irradiated TRISO Fuel Radioactivity from Multiple Advanced Reactor Designs

The SCALE code system was used to model, deplete, and compare several different TRISO-fueled reactor designs: a helium-cooled prismatic reactor, a helium-cooled pebble-bed reactor (PBR), and a Fluoride-Lithium-Beryllium (FLIBE) molten-salt-cooled PBR. The purpose of this comparison was to understand how differences in the reactor designs affect the radioactivity of the fuel after discharge and whether those differences are significant. First, the various reactor designs were build and depleted in the TRITON module for each design and fuel enrichment. Then, the TRITON outputs were used to create burn-up dependent reactor libraries. These libraries were then used by ORIGEN to determine the activities of discharged fuel for each reactor, which were compared to generic Westinghouse 17x17 fuel. Overall, the results showed that short term decays are dominated by reactors with higher operating powers, and the reactor type, initial fuel enrichment, and maximum burn-up are of only secondary importance. Although this analysis only focuses on activities in becquerels, these dependencies are consistent with the expected behavior of decay heat. However, analysis of long- term time periods post-irradiation shows that the reactor type and maximum burn-up have strong impacts on the activities; initial fuel enrichment has a secondary impact while operating power is inconsequential. These results would be useful for analyses, such as dose assessment and modeling in post- release scenarios; normal fuel handling operations; and spent fuel transport, storage and disposal. Of particular interest, the results in this report show that analyses that focus on spent nuclear fuel of advanced reactors need to consider each parameter carefully. Unsurprisingly, if the correct operating power is not used in short term analyses, the results will not be correct. Perhaps unexpectedly, however, if the correct reactor type is not used, then the long term results will also be incorrect, especially for areas such as permanent disposal. Even though this report focuses on the total activity of nuclear fuel, it provides initial results on the effects of various input parameters and also provides a framework to extend the work into other analyses of spent fuel from advanced reactors, especially those employing TRISO fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Design and Analysis of the Open-Surface Slow Li Flow Divertor and Comparison to the Fast Li Flow Divertor

In the ongoing U.S. project, “Liquid Metal Plasma Facing Components,” sponsored by the U.S. Department of Energy, efforts have been taken to develop two open-surface divertor designs for the Fusion Nuclear Science Facility using liquid lithium (Li) as a heat and particle flux removal media. The main focus of this study is the design and analysis of a slow (~1 mm/s) and thin (<1 mm) open-surface Li flow divertor with a Li-cooled substrate, which is then compared with an earlier design of a fast (up to 10 m/s) and thick (~0.5 cm) Li flow divertor with the substrate cooled with helium. The slow Li flow divertor design is based on the original LiWall concept developed at the Princeton Plasma Physics Laboratory. Such a thin and slow Li layer can remove the particle flux by reducing the recycling flux, while the heat flux is removed mainly through the heat sink located beneath. In the present study, the heat sink is provided through a Li cooling flow inside the substrate of reduced activation ferritic/martensitic steel. By performing a multiphysics analysis with COMSOL that included liquid-metal magnetohydrodynamics (MHD), heat transfer, and structural mechanics, the impact of various factors on the divertor heat removal capability, such as Li flow velocity, MHD effects, and inlet velocity boundary condition, were examined. Based on comparisons of the two divertor designs, it was shown that the fast-flow divertor significantly outperformed the slow-flow design, whose heat removal capability was limited to ~1 to 2 MW/m2.

Jiang, Yuchen [ORNL] (ORCID:0000000283280284)↗
DOE OSTI2023

A comparison of Eulerian and semi-Lagrangian approaches for modeling stream water quality

This paper describes and compares some of the advantages and limitations of Eulerian and Lagrangian approaches to water quality modeling and introduces a mixed Eulerian-Lagrangian (or semi-Lagrangian) methodology that captures the strengths of both approaches. The semi-Lagrangian modeling approach is applied to advection-dominated rivers, and flexibly ensures unconditional stability for all time step durations and grid segmentations. The semi-Lagrangian modeling approach is demonstrated by applying it to estimate the dissolved oxygen concentrations in the Sava River in Slovenia, focusing on aspects of the methodology and findings that would be of broad interest to managers of water quality in fluvial water bodies. Results of comparisons of the semi-Lagrangian model with the Eulerian-based QUAL2K model in steady and non-steady scenarios demonstrate that while both models are fully capable of producing satisfactory results when optimally configured, the semi-Lagrangian approach offers accuracy and stability without sensitivity to the interaction of time step size and computational grid segmentation scheme.

Sava River↗
DOE OSTI2026

Comparison of juvenile white sturgeon blade strike survival through a conventional and a novel hydropower turbine

Hydropower is a globally important renewable energy source with profound impacts on aquatic organisms, including direct injury and mortality to fish that move through hydropower facilities. While it is expected that the characteristics of turbine design and operation influence the type, frequency, and magnitude of injuries suffered by entrained fish, few direct comparisons exist to illustrate differences in fish survival outcomes for conventional and novel turbine designs. This study evaluated blade strike-associated injury and mortality rates for juvenile white sturgeon (Acipenser transmontanus) passed through a model-scale turbine equipped with a runner having conventional blade profiles (thin, straight leading edges) and the same turbine equipped with a runner designed for improved fish survival (thick, slanted leading edges). In both trials, head and turbine runner rotational speeds were matched to produce five blade peripheral speeds between 15.0 and 27.6 m/s. High-speed video was captured for all turbine passage events. Fish were assessed for injuries following passage, and mortalities were assessed immediately after passage and after 48 h. Conventional runner passage resulted in 42% to 78% survival after 48 h, with approximately one-third of all tested fish killed by severing. Under the same test conditions, immediate and 48 h survival rates through the novel runner were 100% except at the highest speed condition (95.6% survival at 48 h). These results clearly indicate that turbine design has a profound effect on fish survival, and suggest that thoughtful redesign of hydropower equipment could significantly ameliorate the risks of blade strike mortality.

13 HYDRO ENERGY↗
DOE OSTI2023

Ground and excited states of even-numbered Hubbard ring at half-filling: comparison of the extended Gutzwiller approach with exact diagonalization

It remains a great challenge in condensed matter physics to develop a method to treat strongly correlated many-body systems with balanced accuracy and efficiency. We introduce an extended Gutzwiller (EG) method incorporating a manifold technique, which builds an effective manifold of the many-body Hilbert space, to describe the ground- and excited-state properties of strongly correlated electrons. Here we systematically apply an EG projector onto the ground and excited states of a non-interacting system. Diagonalization of the true Hamiltonian within the manifold formed by the resulting EG wavefunctions gives the approximate ground and excited states of the correlated system. To validate this technique, we implement it on even-numbered fermionic Hubbard rings at half-filling with periodic boundary conditions, and compare the results with the exact diagonalization (ED) method. The EG method is capable of generating high-quality ground and low-lying excited state wavefunctions, as evidenced by the high overlaps of wavefunctions between the EG and ED methods. Favorable comparisons are also achieved for other quantities including the total energy, the double occupancy, the total spin and the staggered magnetization. With the capability of accessing the excited states, the EG method can capture the essential features of the one-electron removal spectral function that contains contributions from states deep in the excited spectrum. Finally, we provide an outlook on the application of this method on large extended systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Understanding JET-C quiescent phases with edge harmonic magnetohydrodynamic activity and comparison with behaviour under ITER-like wall conditioning

An analysis of edge localised mode-free (quiescent) H-mode discharges exhibiting edge harmonic magnetoydrodynamic activity in the JET-carbon wall machine is presented. It is observed that the otherwise quiescent pulses with multiple-n harmonic oscillations are sustained until a threshold in pedestal electron density and collisionality is crossed. The macroscopic pedestal parameters associated with the quiescent phase are compared with those of a database of JET-ELMy discharges with both carbon and ITER-like wall (ILW). This comparison provides the identification of the existence regions in the relevant pedestal and global plasma parameters for edge harmonic oscillations (EHOs) in JET plasmas. Although the ELMy database scans pedestal collisionality and β values typical of ET-carbon quiescent operation, shaping and current are not simultaneously compatible with EHO existence. Nevertheless, ILW operation with JET-carbon quiescent-like parameters could in principle be achieved, and improved pedestal performance could be observed in more recent JET-ILW pulses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Spectro-temporal comparisons of optical emission, absorption, and laser-induced fluorescence for characterizing ns and fs laser-produced plasmas

We performed simultaneous measurement of absorption, emission, and laser-induced fluorescence spectroscopic signatures for determining nanosecond and femtosecond laser-produced plasma’s physical properties throughout its lifecycle. Plasmas are produced by focusing either ~6 ns, 1064 nm pulses from an Nd:YAG or ~35 fs, ~800 nm pulses from a Ti:Sapphire laser on an Inconel target that contains Al as a minor alloying addition. A continuous-wave narrowband tunable laser was used for performing absorption and fluorescence spectroscopy while a fast-gated detection system was used for emission spectroscopy. The temporal evolution of emission, fluorescence, and absorbance of Al transitions are compared for both ns and fs laser-produced plasms. Time-resolved absorbance was also used for evaluating linewidth, lineshape, temperature, and column-averaged atomic number density at late times of ns and fs plasma evolution. Our results demonstrate that lower and excited-state populations of fs laser produced plasmas are short-lived in comparison to those in ns plasmas. The lower state population is observed to reach a maximum value earlier in time for the fs plasma versus the ns plasma, while the kinetic temperature for the ns plasma was higher than for the fs plasma at most times of the plasma evolution.

(280.4788) Optical Sensing, (300.2140) Emission Sp↗
DOE OSTI2020

Effect of net direct current on the properties of radio frequency sheaths: simulation and cross-code comparison

In order to understand, predict and control ion cyclotron range of frequency (ICRF) interactions with tokamak scrape-off layer plasmas, computational tools which can model radio frequency (RF) sheaths are needed. In particular, models for the effective surface impedance and DC rectified sheath potentials may be coupled with full wave RF simulation codes to predict self-consistent wave fields near surfaces and the resulting power dissipation and plasma-material interactions from ion sputtering. In this study, previous work assuming zero net DC current flow through the sheath is generalized to allow the surface to collect net positive or negative current, as is often observed in experiments. The waveforms, DC potential and RF admittance are investigated by means of analytical theory, nonlinear fluid and particle-in-cell (PIC) codes. Cross-code comparisons provide detailed model verification and elucidate the roles of ion and electron kinetics. When the sheath draws negative (positive) DC current, the voltage rectification is reduced (increased) compared with the zero-current case, and the magnitude of both the real and imaginary parts of the admittance are increased (reduced). In conclusion, a previous four-input parametrization of the sheath rectification and admittance properties is generalized to include a fifth parameter describing the DC sheath current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗