Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power modeling and estimation”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Advanced Sensor Deployment for Distribution System State Estimation and Fault Identification: Preprint

Distribution systems are currently facing steep operational challenges as a result of the rapidly increasing integration of renewables and other distributed energy resources (DERs)at both the primary and secondary circuit levels. Distribution utilities and system operators have traditionally had some visibility of their primary circuits using low-frequency supervisory control and data acquisition systems, and they have had very poor if not zero visibility of the secondary circuits where the presence of DERs is constantly increasing. Therefore, this paper presents simulation studies to demonstrate the benefits of an advanced, high-fidelity sensor technology, called as the Meta-Alert System (MAS), developed by Electrical Grid Monitoring Ltd. (EGM), on the distribution grid. First, a reliable model of the EGM sensors is developed, and then two use cases - distribution system state estimation (DSSE) and fault identification - are simulated to evaluate the performance of the MAS technology. Simulation results on the Electric Power Research Institute J1feeder demonstrate that the MAS can effectively participate in system-level DSSE programs and can detect and locate faults faster than traditional distribution protection schemes.

distribution system monitoring↗
DOE OSTI2025

Shallow Geothermal Resources for Cooling Applications at the University of Hawai'i

Drilling activities account for 30% to 57% of the cost to develop and install a geothermal plant. Therefore, an accurate representation of the cost to drill a well is paramount in techno-economic analysis to determine the feasibility of a geothermal power project. In 2022, the National Renewable Energy Laboratory (NREL) endeavored to revise the U.S. Department of Energy (DOE) GeoVision baseline drilling cost curves due to extensive improvement in drilling rates at the Utah Frontier Observatory Research in Geothermal Energy (FORGE) demonstration site. That effort did not culminate in the recommendation of new curves because the actual project costs did not match the reported performance improvements and were at or above the GeoVision baseline. The need for another iteration of this analysis has arisen from industry record drilling performance reported by recent commercial field-scale and demonstration projects, including Fervo Energy's Cape Station, the Utah FORGE 16B(78)-32 demonstration and the Geysers Power Company's GDC-36 demonstration. Therefore, in this work, we have estimated the resulting industry average rate of penetration (ROP) and bit life and applied these parameters as inputs to the Well Cost Simplified model used in the GeoVision analysis. The resulting revised cost curves show a significant decline from the GeoVision baseline. For vertical wells, the magnitude of cost decline ranges between 12% and 24% while for deviated wells, cost reductions between 18% and 26% are estimated. The revised cost curves are in good agreement with actual cost data and therefore, quantify the economic impact of the utilization of (and advances in) polycrystalline diamond compact (PDC) bit technology and the application of physics-based methodologies that optimize mechanical specific energy.

building cooling↗
DOE OSTI2025

2025 Geothermal Drilling Cost Curves Update: Preprint

Drilling activities account for 30% to 57% of the cost to develop and install a geothermal plant. Therefore, an accurate representation of the cost to drill a well is paramount in techno-economic analysis to determine the feasibility of a geothermal power project. In 2022, the National Renewable Energy Laboratory (NREL) endeavored to revise the U.S. Department of Energy (DOE) GeoVision baseline drilling cost curves due to extensive improvement in drilling rates at the Utah Frontier Observatory Research in Geothermal Energy (FORGE) demonstration site. That effort did not culminate in the recommendation of new curves because the actual project costs did not match the reported performance improvements and were at or above the GeoVision baseline. The need for another iteration of this analysis has arisen from industry record drilling performance reported by recent commercial field-scale and demonstration projects, including Fervo Energy’s Cape Station, the Utah FORGE 16B(78)-32 demonstration and the Geysers Power Company’s GDC-36 demonstration. Therefore, in this work, we have estimated the resulting industry average rate of penetration (ROP) and bit life and applied these parameters as inputs to the Well Cost Simplified model used in the GeoVision analysis. The resulting revised cost curves show a significant decline from the GeoVision baseline. For vertical wells, the magnitude of cost decline ranges between 12% and 24% while for deviated wells, cost reductions between 18% and 26% are estimated. The revised cost curves are in good agreement with actual cost data and therefore, quantify the economic impact of the utilization of (and advances in) polycrystalline diamond compact (PDC) bit technology and the application of physics-based methodologies that optimize mechanical specific energy.

15 GEOTHERMAL ENERGY↗
NASA NTRS1994

The Uranian aurora and its relationship to the magnetosphere

About 32 h of Voyager Ultraviolet Spectrometer (UVS) observations of Uranus H2 band airglow emission (875 less than or equal to lambda less than or equal to 1115 A) have been analyzed using the singular value decomposition (SVD) approach to inversion, producing an intensity map showing aurora at both magnetic poles. An H Lyman alpha aurora may also be present but is difficult to separate from scattered solar and local interstellar medium components. SVD analysis of variance shows that the intensity estimate is significantly larger than the error estimate over both Uranographic poles and part of the equatorial region, fortuitously including both magnetic polar regions. The Goddard Space Flight Center Q(sub 3) magnetic field model correctly predicts that the aurora should be larger in area and emit more power at the weaker N magnetic pole than at the stronger S magnetic pole. However, the auroral emissions are quite localized in magnetic longitude and so do not form complete auroral ovals. The brightest auroral emission at each magnetic pole is confined to a range of approximately 90 deg of magnetic longitude centered on the magnetotail direction, at moderate magnetic L parameter (5 less than or equal to L less than or equal to 10), but some emission at each pole is distributed over a range of more than 180 deg of longitude. The magnetic longitudes of the aurora are completely inconsistent with the 'windshield wiper' effect for either ions or electrons, indicating that some other effect, such as rapid depletion of the population of precipitating particles of highly localized strong pitch-angle diffusion, may be acting to localize emission. The low apparent L of the precipitating particles indicates that their energies may be less than or equal to 10 keV. Hence magnetospheric convection is likely to be important, and thus particles exciting the aurora may not remain on constant L shells. The precipitating particles may be a relatively low-energy population at high L that is heated to aurora-exciting energy by adiabatic compression during convection to low L. We estimate that the total auroral power output at H Lyman alpha and shorter wavelengths is about 3 x 10(exp 9) to 7 x 10(exp 9) W, requiring about 10 times that much power for excitation.

Herbert, Floyd↗
DOE OSTI2024

Wake Effects in Lower Carbon Future Scenarios

In August 2022, the U.S. Congress passed the Inflation Reduction Act (IRA), which intended to accelerate U.S. decarbonization, clean energy manufacturing, and deployment of new power and end-use technologies. The National Renewable Energy Laboratory has examined possible scenarios for growth by 2050 resulting from the IRA and other emissions reduction drivers and defined several possible scenarios for large-scale wind deployment. These scenarios incorporate large clusters of turbines operating as wind farms grouped around existing or likely transmission lines which will result in wind farm wakes. Using a numerical weather prediction (NWP) model, we assess these wake effects in a domain in the U. S. Southern Great Plains for a representative year with four scenarios in order to validate the simulations, estimate the internal wake impact, and quantify the cluster wake effect. Herein, we present a validation of the ”no wind farm” scenario and quantify the internal waking effect for the ”ONE” wind farm scenario. Future work will use the “MID” scenario (more than 8000 turbines) and the “HI” scenario (more than 16,000 turbines) to quantify the effect of cluster wakes or inter-farm wakes on power production.

17 WIND ENERGY↗
DOE OSTI2021

Analyzing Gas Turbine-Generator Performance of the Hybrid Power System

We report kinetic energy available from the turbomachinery may boost short-term ramp rates during load shifting between generator and fuel cell that may be exploited to optimize the performance of turbine-fuel cell hybrid systems during load following. The paper starts with modeling a gas turbine, gearbox, and generator sitting on the rigid shafts as a two-mass model primarily for estimating the rotational inertia available from the turbogenerator. The gas turbine and its environment impose constraints of available space allowances, measurements from high-speed rotating parts, lack of accessibility to the shaft of measure, making replacements and/or required installation adjustments of the torque transducer difficult. It is demonstrated that in the absence of torque measurement, the data containing only a slow excursion of the rotor motion in response to a slight load change is sufficient to determine inertia and damping coefficients from the model. Furthermore, the generator’s reactive power capability was investigated because it serves as an asset in integrated hybrid systems to meet dynamic reactive power demands, reducing the burden on the fuel cell’s power electronics control in generating reactive power.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Cosmology From CMB Lensing and Delensed EE Power Spectra Using 2019-2020 SPT-3G Polarization Data

From CMB polarization data alone we reconstruct the CMB lensing power spectrum, comparable in overall constraining power to previous temperature-based reconstructions, and an unlensed E -mode power spectrum, with clear detections of the third through tenth acoustic peaks. The observations, taken in 2019 and 2020 with the South Pole Telescope (SPT) and the SPT-3G camera, cover 1500 deg 2 at 95, 150, and 220 GHz with arcminute resolution and roughly 4.9 µ K-arcmin coadded noise in polarization. The power spectrum estimates, together with systematic parameter estimates and a joint covariance matrix, follow from a Bayesian analysis using the Marginal Unbiased Score Expansion (MUSE) method. The E -mode spectrum at ℓ > 2000 and lensing spectrum at L > 350 are the most precise to date. Assuming the ΛCDM model, and using only these SPT data and priors on τ and absolute calibration from Planck, we find H 0 = 66.81 ± 0.81 km/s/Mpc, comparable in precision to the Planck determination and in 5.4 σ tension with the most precise H 0 inference derived via the distance ladder. We also find S 8 ≡ σ 8 (Ω m /0.3) 0.5 = 0.850 ± 0.017, providing further independent evidence of a slight tension with low-redshift structure probes. The ΛCDM model provides a good simultaneous fit to the combined Planck, ACT, and SPT data, and thus passes a powerful test. Combining these CMB datasets with BAO observations, we explore extensions to the ΛCDM model. We find that the effective number of neutrino species, spatial curvature, and primordial helium fraction are consistent with standard model values, and that the 95% confidence upper limit on the neutrino mass sum is 0.075 eV, close to the minimum sum expected from observations of solar and atmospheric neutrino oscillations. The SPT data are consistent with the somewhat weak (< 3 σ ) preference for excess lensing power seen in Planck and ACT data relative to predictions of the ΛCDM model given the combined Planck, ACT, and BAO data sets. Finally, we also detect at greater than 3 σ the influence of non-linear evolution in the CMB lensing power spectrum and discuss it in the context of the S 8 tension. Forthcoming SPT-3G analyses will feature deeper and wider observations in temperature and polarization, providing even tighter constraints and more powerful tests of the ΛCDM model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Sensor selection and tool wear prediction with data‐driven models for precision machining

Abstract Estimation of tool wear in precision machining is vital in the traditional subtractive machining industry to reduce processing cost, improve manufacturing efficiency and product quality. In this vein, fusion of time and frequency‐domain features of commonly sensed signals can provide an early indication of tool wear and improve its prediction accuracy for prognostics and health management. This paper presents a data‐driven methodology and a complete tool chain for the inference of precision machining tool wear from fused machine measurements, such as cutting force, power, audio and vibration signals, and quantify the usefulness of each measurement. Indicators of tool wear are extracted from time‐domain signal statistics, frequency‐domain analysis, and time‐frequency domain analysis. Correlation coefficients between the extracted features (indicators) and the tool wear are used to select the most informative features. Principal Component Analysis and Partial Least‐Squares are used to reduce the dimensionality of the feature space. Regression models, including linear regression, support vector regression, Decision tree regression, neural network regression and Gaussian process regression, are used to predict the tool wear using data from a Haas milling machine performing spiral boss face milling. The performance of the regression models based on subsets of sensors validates the preliminary estimates about the saliency of the sensors. The experimental results show that the proposed methods can predict the machine tool wear precisely, with readily available sensor measurements. Neural network and Gaussian process regression were able to achieve good estimates of tool wear at different machine operating conditions. The most informative signal in predicting tool wear was shown to be the vibration signal. Time‐frequency domain features were the most informative features among the combination of features of three domains. In addition, using partial least squares components extracted from the original features of signals led to higher prediction accuracy.

Han, Seulki↗
DOE OSTI2020

Mitigating Cascading Outages in Severe Weather Using Simulation-Based Optimization

Severe weather events can trigger cascading power outages and lead to significant losses. In this work, we investigate cascading outage mitigation under severe weather conditions. Given day-ahead weather forecasts and component failure models, we aim to identify a set of power lines that can be hardened to minimize the expected impact of potential cascading outages. Since the expected load shedding cannot be expressed as an explicit function of line hardening decisions and system states, we developed a cascading outage simulator to estimate the expected value of load shedding under various initial weather-related disruption scenarios generated using a weather forecast. To avoid massive enumeration of all possible combinations of line hardening decisions and reduce the simulation efforts, we employed an efficient simulation-based optimization approach that quickly identifies the (near) optimal line hardening decisions in the presence of both large simulation noises due to the highly variable initial disturbances and system states, and significant randomness in the subsequent cascades. Furthermore, the algorithm is also able to utilize parallel computing to dramatically reduce computation time to support decision making in preparation for severe weather conditions. We performed a case study on the Northeast Power Coordinating Council (NPCC) 140-bus system model to demonstrate that our approach can significantly improve power grid resilience to adverse weather events.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1994

Durability evaluation of ceramic components using CARES/LIFE

The computer program CARES/LIFE calculates the time-dependent reliability of monolithic ceramic components subjected to thermomechanical and/or proof test loading. This program is an extension of the CARES (Ceramics Analysis and Reliability Evaluation of Structures) computer program. CARES/LIFE accounts for the phenomenon of subcritical crack growth (SCG) by utilizing the power law, Paris law, or Walker equation. The two-parameter Weibull cumulative distribution function is used to characterize the variation in component strength. The effects of multiaxial stresses are modeled using either the principle of independent action (PIA), the Weibull normal stress averaging method (NSA), or the Batdorf theory. Inert strength and fatigue parameters are estimated from rupture strength data of naturally flawed specimens loaded in static, dynamic, or cyclic fatigue. Application of this design methodology is demonstrated using experimental data from alumina bar and disk flexure specimens which exhibit SCG when exposed to water.

Nemeth, Noel N.↗
DOE OSTI2023

Causal CCF Parameter Estimations 2020

This report documents the quantitative results of the causal common-cause failure (CCF) parameter estimations for the failure cause groups “component,” “design,” “environment,” “human,” and “other,” based on CCF data through 2020 in the U.S. Nuclear Regulatory Commission (NRC) CCF database: https://rads.inl.gov/Pages/CCF.aspx. This report utilizes the same data period (2006–2020) and CCF templates as INL/EXT-21-62940, Revision 1, CCF Parameter Estimations, 2020 Update. The 2015 causal CCF prior distributions for the specific failure cause groups (instead of the 2015 generic CCF prior distributions) were used in this report to estimate the associated causal CCF parameters. All the 2015 causal CCF prior distributions and generic CCF prior distributions were developed in INL/EXT-21-43723, Developing Generic Prior Distributions for Common Cause Failure Alpha Factors and Causal Alpha Factors, using CCF data from 1997 to 2015. These quantitative results were developed to support the causal alpha factor model and should be used as appropriate in probabilistic risk assessment (PRA) studies such as the NRC Significance Determination Process for commercial nuclear power plants in the United States.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Confronting interacting dark radiation scenarios with cosmological data

Dark radiation (DR) is generally predicted in new physics scenarios that address fundamental puzzles of the Standard Model or tensions in the cosmological data. Cosmological data have the sensitivity to constrain not only the energy density of DR but also whether it is interacting. In this paper, we present a systematic study of five types of interacting DR (free-streaming, fluid, decoupling, instantaneous decoupling, and recoupling DR) and their impact on cosmological observables. We modify the Boltzmann hierarchy to describe all these types of interacting DR under the relaxation time approximation. We, for the first time, robustly calculate the collision terms for recoupling scalar DR and provide a better estimation of the recoupling transition redshift. We demonstrate the distinct features of each type of DR on the cosmic microwave background and matter power spectra. We perform Markov-chain Monte Carlo scans using the Planck 2018 data and baryon acoustic oscillation data. Assuming no new physics in the standard model neutrino sector, we find no statistically significant constraints on the couplings of DR, although there is a slight preference for the fluidlike limit of all the cases. In the case of instantaneous decoupling DR, this limit corresponds to a late transition redshift around recombination. Furthermore, the ΔN eff constraint varies marginally depending on the type of DR.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Common-Cause Component Group Modeling Issues in Probabilistic Risk Assess

Common cause failures (CCFs) have been recognized as significant risk contributors in probabilistic risk assessments (PRAs) for commercial nuclear power plants (NPPs) since 1980s. A series of reports including those of United States Nuclear Regulatory Commission (NRC) regulations (NUREGs) have been published from late 1980s to early 2000s to provide guidelines for performing CCF event data analysis and modeling CCF in PRA. However, there are still some issues existed in CCF modeling and CCF parameter estimations. One of such issues is the modeling of multiple common-cause component group (CCCG) for the same component in a PRA model. This paper exams the CCCG definition and the requirements from the PRA standard, investigates the CCCG modeling practices and issues existing in PRAs, and presents systematic approaches and guidelines to model CCCG in PRA properly and consistently.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS2016

Modeling the Swift Bat Trigger Algorithm with Machine Learning

To draw inferences about gamma-ray burst (GRB) source populations based on Swift observations, it is essential to understand the detection efficiency of the Swift burst alert telescope (BAT). This study considers the problem of modeling the Swift / BAT triggering algorithm for long GRBs, a computationally expensive procedure, and models it using machine learning algorithms. A large sample of simulated GRBs from Lien et al. is used to train various models: random forests, boosted decision trees (with AdaBoost), support vector machines, and artificial neural networks. The best models have accuracies of greater than or equal to 97 percent (less than or equal to 3 percent error), which is a significant improvement on a cut in GRB flux, which has an accuracy of 89.6 percent (10.4 percent error). These models are then used to measure the detection efficiency of Swift as a function of redshift z, which is used to perform Bayesian parameter estimation on the GRB rate distribution. We find a local GRB rate density of n (sub 0) approaching 0.48 (sup plus 0.41) (sub minus 0.23) per cubic gigaparsecs per year with power-law indices of n (sub 1) approaching 1.7 (sup plus 0.6) (sub minus 0.5) and n (sub 2) approaching minus 5.9 (sup plus 5.7) (sub minus 0.1) for GRBs above and below a break point of z (redshift) (sub 1) approaching 6.8 (sup plus 2.8) (sub minus 3.2). This methodology is able to improve upon earlier studies by more accurately modeling Swift detection and using this for fully Bayesian model fitting.

gamma-ray burst: general – gamma-rays: general â↗
DOE OSTI2025

Manufacturing Cost Analysis for PEM Electrolyzers and Perspectives for Future Cost Reduction

Electrolyzer capital costs strongly influence the total levelized cost of hydrogen production and have implications for hydrogen deployment. Current electrolyzer costs are high, and large cost reductions may be needed to achieve competitive hydrogen costs and targets. Understanding pathways for cost reduction via R&D and deployment is a critical research area for informed energy planning and enabling hydrogen use. This work presents bottom-up cost estimates of polymer electrolyte membrane (PEM) electrolyzer systems tied to design specifications and discusses perspectives for cost reduction opportunities based on ongoing research. We use a detailed manufacturing and process model for a 1 MW PEM electrolyzer stack and balance of plant (BOP) for rigorous cost estimation. This allows for robust estimates of component and manufacturing costs and examination of key cost contributors. Stack costs are dominated by material costs such as iridium and platinum catalysts, especially at high manufacturing rates; power electronics and hydrogen purification equipment are the largest contributors to BOP cost. At higher manufacturing rates, better equipment utilization could reduce stack costs significantly, and we estimate that experience and bulk purchasing will allow for cost reductions to some BOP components. Still, many well-established BOP technologies and stack material costs are less likely to see significant cost reductions at high manufacturing rates. As such, manufacturing scale is limited in how much it can reduce electrolyzer costs, and additional advances for cost reduction may be needed to achieve cost targets. It will likely take many combined strategies to achieve significant cost reductions for electrolyzers and enable low-cost hydrogen production. We can use our manufacturing cost model to quantify potential cost reductions from the considerations described above and demonstrate pathways to lower cost electrolyzers. This allows for better understanding of cost reduction strategies and enables more informed research, development, and deployment for electrolyzers.

cost↗
DOE OSTI2021

Model-based Cyber-attack Detection for Voltage Source Converters in Island Microgrids

With the upgrading of communication and control in microgrids, cyber threats on the power converters are increasing. Here, in this paper, a novel model-based cyber-attack detection methodology is proposed for each voltage source converter in microgrids. The Harmonics State Space Matrix (H-Matrix) is used to build the closed-loop transfer function, providing model-based estimation with the grid voltage and control reference. Then, the space phase model (SPM) is used to calculate residual to detect cyber-attacks in voltage source converter (VSC). As controller and converter parameters are included in the H-matrix, the proposed method also could detect cyber-attacks when voltage fluctuation occurs in the grid. To verify the feasibility of the proposed method, several attacks targeting VSC in an islanded microgrid are simulated in MATLAB. Simulation and detection results are introduced to demonstrate the resilience and feasibility of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Dark Energy Survey Year 3 results: Simulation-based cosmological inference with wavelet harmonics, scattering transforms, and moments of weak lensing mass maps. II. cosmological results

Here, we present a simulation-based cosmological analysis using a combination of Gaussian and non-Gaussian statistics of the weak lensing mass (convergence) maps from the first three years of the Dark Energy Survey. We implement the following: (1) second and third moments; (2) wavelet phase harmonics; (3) the scattering transform. Our analysis is fully based on simulations, spans a space of seven 𝑤 Cold Dark Matter (𝑤⁢ CDM) cosmological parameters, and forward models the most relevant sources of systematics inherent in the data: masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration. We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a simulation-based inference approach. Including and combining different non-Gaussian statistics is a powerful tool that strongly improves constraints over Gaussian statistics (in our case, the second moments); in particular, the figure of merit (𝑆 8 , Ω m ) is improved by 70% (Λ ⁢CDM) and 90% (𝑤 ⁢CDM). When all the summary statistics are combined, we achieve a 2% constraint on the amplitude of fluctuations parameter 𝑆 8 ≡ 𝜎 8 ⁢(Ω m /0.3) 0.5 , obtaining 𝑆 8 = 0.794 ±0.017 (Λ⁢ CDM) and 𝑆 8 = 0.817 ±0.021 (𝑤 ⁢CDM), and a ∼10% constraint on Ω m , obtaining Ω m =0.259 ±0.025 (Λ ⁢CDM) and Ω m = 0.273 ±0.029 (𝑤⁢ CDM). In the context of the 𝑤⁢ CDM scenario, these statistics also strengthen the constraints on the parameter 𝑤, obtaining 𝑤 <−0.72. The constraints from different statistics are shown to be internally consistent (with a 𝑝-value>0.1 for all combinations of statistics examined). We compare our results to other weak lensing results from the first three years of the Dark Energy Survey data, finding good consistency; we also compare with results from external datasets, such as planck constraints from the cosmic microwave background, finding statistical agreement, with discrepancies no greater than <2.2⁢𝜎.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Modeling of a 43 GHz microplasma within a photonic crystal

The electron density (n e ) and argon metastable density (1 s 5 ) of a 43GHz microplasma are obtained using the zero-dimensional plasma kinetics solver (ZDPlasKin®) for an incident power of 160 and 1000mW and argon pressures of 10–600Torr (1.3 ×10 3 to 8.0 × 10 4 Pa). These simulations are compared with previously published experimental data. To determine the self-consistent electric field in the modeled plasma, the three-dimensional millimeter wave fields are computed as a function of electron density using ANSYS EM19.2, HFSS®. This electromagnetic field model is coupled to ZDPlasKin such that any increase in the simulated plasma density correctly attenuates the simulated electric field within the microplasma. The electron density is found to be sensitive to argon gas temperature, so a two-zone temperature model was needed to obtain agreement with experimental measurements. The temperature in the central core of the microplasma was determined by a previous experimental study. That temperature was used as an input to the model for the simulation of volume recombination losses. The outer regions of the microplasma are assumed to be in equilibrium with the walls (300K). This second temperature was used in the model to determine diffusion losses. The modeled electron and metastable densities are of the order 10 20 and 10 18 m –3 , respectively. This is in good agreement with those measured experimentally as long as the two-zone temperature model is used. In the absence of a hot gas core, the modeled three-body recombination rates are excessive and the simulation severely under-predicts the electron density and over-estimates the metastable density. We conclude that the millimeter wave microplasma has a hot core (2500K at 600 Torr) that rarifies the argon gas and effectively reduces three-body recombination. Furthermore, this allows one to achieve high electron density on the order of 10 20 m –3 with only 100mW of wave power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗