Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “frequency ratio”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Python Instrument Control and Data Acquisition Suite for Reproducible Research

Tools that standardize and automate experimental data collection are needed for greater confidence in research results. The National Synchrotron Light Source-II (NSLS-II) has generated an open-source Python data acquisition, management, and analysis software suite that automates X-ray experiments and collects an experimental record that facilitates complete reproducibility. Here, we show that the NSLS-II tools are not only useful for X-ray science at large-scale facilities by presenting an add-on package that adapts these tools for use in a small laboratory with common physics and electrical engineering instruments. The composite software suite eases and automates the execution of experiments, records extensive metadata, stores data in portable containers, and speeds up the analysis through tools for comprehensive searches. In total, this software suite increases the reproducibility of laboratory experiments. We demonstrate the software via the evaluation of two lock-in amplifiers-the miniature ADA2200 and the ubiquitous Stanford Research Systems (SRS) SR810. The frequency resolution, signal-to-noise ratio, and dynamic reserve of the lock-in amplifiers are measured and presented. The usage of the software suite is described throughout these measurements so that the reader can implement the tools in their lab.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Pre-Sinkhole Seismicity at the Napoleonville Salt Dome: Implications for Local Seismic Monitoring of Underground Caverns

The formation of a sinkhole at the Napoleonville salt dome located near Bayou Corne in Assumption Parish, Louisiana, was preceded by an increase in seismicity, as observed 10 km away by a continuously recording broadband seismic sensor. The seismicity rate in the days prior to the sinkhole formation contained hundreds of long-period events, some of which could be detected using a power detector (short-term to long-term energy ratio) tuned to the duration and frequency content of the seismicity. A characteristic of some of these events is a resonance near 0.4 Hz, which could be due to either an acoustic resonance within the brine-filled cavity or possibly due to unsteady fluid flow through a crack. In this work, a more exhaustive catalog was obtained using some of the events detected by the power detector as templates in a correlation detector. Furthermore, the full pre-sinkhole sequence of detections could offer lessons in remote monitoring of underground caverns.

58 GEOSCIENCES↗
DOE OSTI2026

Final Report on the Measurement of the Positive Muon Anomalous Magnetic Moment at Fermilab to 127 ppb

This report details the final measurement of the muon magnetic anomaly, $a_μ=(g_μ-2)/2$, by the Muon $g-2$ experiment at Fermi National Accelerator Laboratory (FNAL), using positive muons collected from 2021 to 2023. The value of $a_μ$ is determined from the ratio of the anomalous spin precession frequency to the shielded proton precession frequency in the muon storage ring magnetic field, combined with external constants known at the 22 ppb level. The new dataset, containing over $2.5$ times the statistics of our previous results, yields $a_μ=116\,592\,0710(162)\times 10^{-12}$ (139 ppb), or $a_μ=116\,592\,0705(148)\times 10^{-12}$ (127 ppb) when combined with our previous results. The new experimental world average, dominated by the measurements at FNAL, is $a_μ^{\text{Exp}}=116\,592\,0715(145)\times 10^{-12}$ (124 ppb).

Aguillard, D. P. [U. Michigan, Ann Arbor] (ORCID:0↗
DOE OSTI2021

Characterizing patterns and variability of building electric load profiles in time and frequency domains

The rapid development of advanced metering infrastructure provides a new data source—building electrical load profiles with high temporal resolution. Electric load profile characterization can generate useful information to enhance building energy modeling and provide metrics to represent patterns and variability of load profiles. Such characterizations can be used to identify changes to building electricity demand due to operations or faulty equipment and controls. In this study, we proposed a two-path approach to analyze high temporal resolution building electrical load profiles: (1) time-domain analysis and (2) frequency-domain analysis. Furthermore, the commonly adopted time-domain analysis can extract and quantify the distribution of key parameters characterizing load shape such as peak-base load ratio and morning rise time, while a frequency-domain analysis can identify major periodic fluctuations and quantify load variability. We implemented and evaluated both paths using whole-year 15-minute interval smart meter data of 188 commercial office building in Northern California. The results from these two paths are consistent with each other and complementary to represent full dynamics of load profiles. The time- and frequency-domain analyses can be used to enhance building energy modeling by: (1) providing more realistic assumptions about building operation schedules, and (2) validating the simulated electric load profiles using the developed variability metrics against the real building load data.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Generalization of Deep-Learning Models for Classification of Local Distance Earthquakes and Explosions across Various Geologic Settings

Although accurately classifying signals from earthquakes and explosions at local distance (<250 km) remains an important task for seismic network operations, the growing volume of available seismic data presents a challenge for analysts using traditional source discrimination techniques. In recent years, deep-learning models have proven effective at discriminating between low-magnitude earthquakes and explosions measured at local distances, but it is not clear how well these models are capable of generalizing across different geological settings. To address the issue of generalization between regions, we train deep-learning models (convolutional neural networks [CNNs]) on time–frequency representations (scalograms) of three-component earthquake and explosion signals from eight different regions in the continental United States. We explore scenarios where models are trained on data from all regions, individual regions, or all but one region. We find that although CNN models trained on individual regions do not necessarily generalize well across different settings, models trained on multiple regions that include diverse path coverage generalize to new regions, with station-level accuracy of up to 90% or more for data sets from unseen regions. In general, CNN-based discrimination models significantly outperform models based on uncorrected P/S ratio (measured in the 10–18 Hz frequency band), even when CNN models are tested on data from entirely unseen regions.

58 GEOSCIENCES↗
DOE OSTI2022

Phonon spectra of a two-dimensional solid dusty plasma modified by two-dimensional periodic substrates

Here, phonon spectra of a two-dimensional (2D) solid dusty plasma modulated by 2D square and triangular periodic substrates are investigated using Langevin dynamical simulations. The commensurability ratio, i.e., the ratio of the number of particles to the number of potential well minima, is set to 1 or 2. The resulting phonon spectra show that propagation of waves is always suppressed due to the confinement of particles by the applied 2D periodic substrates. For a commensurability ratio of 1, the spectra indicate that all particles mainly oscillate at one specific frequency, corresponding to the harmonic oscillation frequency of one single particle inside one potential well. At a commensurability ratio of 2, the substrate allows two particles to sit inside the bottom of each potential well, and the resulting longitudinal and transverse spectra exhibit four branches in total. We find that the two moderate branches come from the harmonic oscillations of one single particle and two combined particles in the potential well. The other two branches correspond to the relative motion of the two-body structure in each potential well in the radial and azimuthal directions. The difference in the spectra between the square and triangular substrates is attributed to the anisotropy of the substrates and the resulting alignment directions of the two-body structure in each potential well.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Measurement of the muon anomalous precession frequency at the Muon g − 2 Experiment at Fermilab

The anomalous magnetic moment of the muon, $a_\mu = \frac{g-2}{2}$ is the fractional deviation of the muon $g$-factor from the value of 2. It emerges as the cumulative effect of the virtual particles participating in the muon interaction with a magnetic field via quantum loop corrections. Its value encodes all the possible interactions between the virtual particles and, for this reason, represents an important test of the Standard Model (SM). In particular, any deviation from the SM theoretical evaluation could be due to new physics contributions. The new Muon $g-2$ (E989) Experiment at Fermilab is currently operating to repeat and improve the previous E821 experiment at Brookhaven National Laboratory (BNL), aiming to reduce the experimental error by a factor of 4 to the final accuracy of 140 parts per billion (ppb). On April 7th, 2021, the E989 collaboration published the first result based on the first year of data taking (Run-1), measuring $a_\mu = 0.001~165~920~40(54)$ with a precision of 460 ppb. The measured value is consistent with the BNL measurement and strengthens the long-standing tension with the data-driven SM prediction to a combined discrepancy of 4.2$\sigma$. On the theory side, however, new efforts involving lattice-QCD techniques are starting to question the current consensus on the theoretical prediction, demanding new improvements on both the experimental and theoretical sides. The E989 collaboration is now finalizing the analysis of Run-2 and Run-3 data and a new publication is expected in the first half of 2023 with a combined statistical uncertainty of 200 ppb. The anomalous magnetic moment $a_\mu$ is measured as the ratio between the muon spin anomalous precession frequency, $\omega_a$, and the average magnetic field experienced by the muons as they circulate in the storage ring. This thesis presents a precession frequency analysis of the Run-1 data and an evaluation of the related systematic uncertainties. A new positron reconstruction developed for the analysis of the subsequent data-taking periods, aiming to reduce some of the major systematic uncertainties of the $\omega_a$ measurement, is presented. The author's involvement in the production of the Run-2/5 data and in the precise calibration of the detectors is discussed. Finally, the complete Run-1 $a_\mu$ result is presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Water vapor isotopic composition during TRACER IOP

The TRACER-Iso project took place from June 1 to September 30, 2022 as part of the intensive observation period (IOP) of the Tracking Aerosol Convection Interactions (TRACER) project at the main project site at the La Porte, Texas, municipal airport. The project sought to use measurements of the stable isotopic composition of water vapor to better understand cloud-aerosol interactions in deep convection in a polluted urban setting. Water vapor isotopologue and humidity measurements were determined using a Picarro L2130 Water Vapor Isotope Analyzer. The instrument was deployed in the guest van and consists of three main components: a water vapor isotope analyzer, and a Los Gatos Research Water Vapor Isotope Standard Source (WVISS), and a Los Gatos Research Dry Air Source. Ambient air samples were collected through an inlet located about one meter above roof of the shipping container. Samples were then delivered to the analyzer through teflon tubing by use of an external pump to minimize the transport time between the inlet and analyzer. The inlet tubing was surrounded by a Briskheat heat trace and insulating material to ensure ambient air samples did not fall below the dew point and result in condensation. The analyzer uses cavity ringdown spectroscopy to report isotopic ratios of ambient air samples at a frequency of 1 Hz. The calibration unit of the instrument was used in conjunction with the dry air source to periodically measure the δ values of standard waters. It uses a nebulizer to push small water droplets into a hot chamber that vaporizes the water without fractionation. This vapor was then transported to the analyzer using a built-in compressor and the dry air source, which allows each standard with known δ values to be measured at a wide range of humidity values for post-measurement calibration of ambient air samples. There are multiple sources of uncertainty introduced at different stages during the data collection and processing. These include (1) instrument precision, (2) uncertainty in the secondary standards, (3) humidity-correction uncertainty, and (4) VSMOW-SLAP calibration uncertainty. Uncertainty from each step is propagated in quadrature to calculate a total uncertainty of each isotopologue. Uncertainty was determined to be 2.3‰ for δD and 0.4‰ for δ18O.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hierarchical Power Flow Control in Smart Grids: Enhancing Rotor Angle and Frequency Stability with Demand-Side Flexibility

Large-scale integration of renewables in power systems gives rise to new challenges for keeping synchronization and frequency stability in volatile and uncertain power flow states. To ensure the safety of operation, the system must maintain adequate disturbance rejection capability at the time scales of both rotor angle and system frequency dynamics. This calls for flexibility to be exploited on both the generation and demand sides, compensating volatility and ensuring stability at the two separate time scales. This article proposes a hierarchical power flow control architecture that involves both transmission and distribution networks as well as individual buildings to enhance both small-signal rotor angle stability and frequency stability of the transmission network. The proposed architecture consists of a transmission-level optimizer enhancing system damping ratios, a distribution-level controller following transmission commands and providing frequency support, and a building-level scheduler accounting for quality of service and following the distribution-level targets. Furthermore, we validate the feasibility and performance of the whole control architecture through real-time hardware-in-loop tests involving real-world transmission and distribution network models along with real devices at the Stone Edge Farm Microgrid.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Multinuclear PFGSTE NMR description of 39 K, 23 Na, 7 Li, and 1 H specific activation energies governing diffusion in alkali nitrite solutions

While pulsed field gradient stimulated echo nuclear magnetic resonance (PFGSTE NMR) spectroscopy has found widespread use in the quantification of self-diffusivity for many NMR-active nuclei, extending this technique to uncommon nuclei with unfavorable NMR properties remains an active area of research. Potassium-39 ( 39 K) is an archetypical NMR nucleus exhibiting an unfavorable gyromagnetic ratio combined with a very low Larmor frequency. Despite these unfavorable properties, this work demonstrates that 39 K PFGSTE NMR experiments are possible in aqueous solutions of concentrated potassium nitrite. Further, analysis of the results indicates that 39 K NMR diffusometry is feasible when the nuclei exhibit spin–lattice and spin–spin relaxation coefficients on the order of 60–100 ms and 50–100 ms, respectively. The diffusivity of 39 K followed Arrhenius behavior, and comparative 23 Na, 7 Li, and 1 H PFGSTE NMR studies of equimolal sodium nitrite and lithium nitrite solutions led to correlations between the enthalpy of hydration with the activation energy governing self-diffusion of the cations and also of water. Realizing the feasibility of 39 K PFGSTE NMR spectroscopy has a widespread impact across energy sciences because potassium is a common alkali element in energy storage materials and other applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Gyrokinetic analysis of inter-edge localized mode transport mechanisms in a DIII-D pedestal

In this study, gyrokinetic simulations are used to study pedestal fluctuations for DIII-D discharge 174082 using the GENE code. Nonlinear local simulations indicate that electron heat flux has contributions from electron temperature gradient-driven transport but at levels insufficient to satisfy power balance. We show that microtearing modes (MTM) and neoclassical transport are likely to account for the remaining observed energy losses in the electron and ion channels, respectively. The MTM instabilities found in the simulations are consistent with the high-frequency fluctuations identified in the magnetic fluctuation data from Mirnov coils. The fluctuation data in this discharge also exhibit a low-frequency band of fluctuations. By modifying the equilibrium profiles and plasma β, simulations produce MHD modes, which may be responsible for these observed low-frequency fluctuations. We compare several metrics involving ratios of fluctuation amplitudes and transport quantities for both MTMs and MHD modes. This analysis suggests that the available data are consistent with the simultaneous activity of both MHD modes and MTMs provided that the former is limited largely to the particle transport channel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Real-space atomic dynamics in metallic liquids investigated by inelastic neutron scattering

Understanding the dynamics of liquids at the atomic level remains a major challenge. Even though viscosity is one of the most fundamental properties of liquids, its atomistic origin is not fully elucidated. Through inelastic neutron scattering experiment on levitated metallic liquid droplets, the time-dependent pair correlation function, the Van Hove function, was determined for Z⁢r 50 ⁢C⁢u 50 and Z⁢r 80 ⁢P⁢t 20 liquids at various temperatures. The time for change in local atomic connectivity, τ LC , which is the timescale of atomic bond cutting and forming, is estimated based on the exponential decay of the nearest neighbor peak of the Van Hove function. At high temperatures above the crossover temperature T A , τ LC is equal to the Maxwell relaxation time, τ M =η/G ∞ , where η is the macroscopic shear viscosity and G ∞ is the high-frequency shear modulus. Below T A the ratio of τ M /τ LC increases with decreasing temperature, indicating increased atomic cooperativity as predicted by molecular dynamics simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

SAIL-Isotopologue (SAIL-iso) Field Campaign Report

The SAIL-iso project took place from November 19, 2022 to June 3, 2023 as part of the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Surface Atmosphere Integrated Field Laboratory (SAIL) project at the main project site in Gothic, Colorado. The project sought to use measurements of the stable isotopic composition of water vapor to better understand cloud-aerosol interactions in a remote mountainous watershed. Water vapor isotopologue and humidity measurements were determined using a Los Gatos Research Triple Water Vapor Isotope Analyzer (TWVIA). The instrument was deployed in the guest van and consists of three main components: a water vapor isotope analyzer, a Los Gatos Research Water Vapor Isotope Standard Source (WVISS), and a Los Gatos Research Dry Air Source. Ambient air samples were collected through an inlet located about one meter above roof of the shipping container. Samples were then delivered to the analyzer through Teflon tubing by use of an external pump to minimize the transport time between the inlet and analyzer. The inlet tubing was surrounded by a Briskheat heat trace and insulating material to ensure ambient air samples did not fall below the dew point and result in condensation. The analyzer uses cavity ringdown spectroscopy to report isotopic ratios of ambient air samples at a frequency of 1 Hz. The calibration unit of the instrument was used in conjunction with the dry air source to periodically measure the δ values of standard waters. It uses a nebulizer to push small water droplets into a hot chamber that vaporizes the water without fractionation. This vapor was then transported to the analyzer using a built-in compressor and the dry air source, which allows each standard with known δ values to be measured at a wide range of humidity values for post-measurement calibration of ambient air samples. Multiple sources of uncertainty are introduced at different stages during the data collection and processing. These include (1) instrument precision, (2) uncertainty in the secondary standards, (3) humidity-correction uncertainty, and (4) Vienna Standard Mean Ocean Water – Standard Light Antarctic Precipitation (VSMOW-SLAP) calibration uncertainty. Uncertainty from each step is propagated in quadrature to calculate a total uncertainty of each isotopologue. Uncertainty was determined to be 2.3‰ for δD and 0.4‰ for δ 18 O.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The High Fraction of Thin Disk Galaxies Continues to Challenge ΛCDM Cosmology

Any viable cosmological framework has to match the observed proportion of early- and late-type galaxies. In this contribution, we focus on the distribution of galaxy morphological types in the standard model of cosmology (Lambda cold dark matter, ΛCDM). Using the latest state-of-the-art cosmological ΛCDM simulations known as Illustris, IllustrisTNG, and EAGLE, we calculate the intrinsic and sky-projected aspect ratio distribution of the stars in subhalos with stellar mass M * > 10 10 M ⊙ at redshift z = 0. There is a significant deficit of intrinsically thin disk galaxies, which however comprise most of the locally observed galaxy population. Consequently, the sky-projected aspect ratio distribution produced by these ΛCDM simulations disagrees with the Galaxy And Mass Assembly (GAMA) survey and Sloan Digital Sky Survey at ≥12.52σ (TNG50-1) and ≥14.82σ (EAGLE50) confidence. The deficit of intrinsically thin galaxies could be due to a much less hierarchical merger-driven build-up of observed galaxies than is given by the ΛCDM framework. It might also arise from the implemented sub-grid models, or from the limited resolution of the above-mentioned hydrodynamical simulations. We estimate that an 8 5 times better mass resolution realization than TNG50-1 would reduce the tension with GAMA to the 5.58σ level. Finally, we show that galaxies with fewer major mergers have a somewhat thinner aspect ratio distribution. Given also the high expected frequency of minor mergers in ΛCDM, the problem may be due to minor mergers. In this case, the angular momentum problem could be alleviated in Milgromian dynamics because of a reduced merger frequency arising from the absence of dynamical friction between extended dark matter halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Unfolder-based single-stage AC-AC conversion system

An power converter includes an unfolder connected to a three-phase source and has an output connection with a positive terminal, a negative terminal and a neutral terminal. The unfolder creates two unipolar piece-wise sinusoidal DC voltage waveforms offset by a half of a period. A three-input converter connected to the unfolder produces a quasi-sinusoidal output voltage across output terminals. Switches of the converter selectively connect the positive, negative and neutral inputs across the output terminals. A PWM controller controls a first duty ratio and a second duty ratio for the converter based on a phase angle of the source and a modulation index generated from an error signal related to a control variable. The duty ratios are time varying with a fundamental frequency of the source. The modulation index relates to output voltage of the converter, peak voltage or current of the source and/or peak current at the output terminals.

Teeneti, Chakridhar Reddy↗
DOE OSTI2025

Adaptation of Wallace's approach to the specific heat of elemental solids with significant intrinsic anharmonicity, particularly the light actinide metals

The quasiharmonic approximation is the most common method for modeling the specific heat of solids; however, it fails to capture the effects of intrinsic anharmonicity. In this study, we introduce the “elastic softening approximation,” an alternative approach to modeling intrinsic anharmonic effects on thermodynamic quantities, which is grounded in Wallace's thermodynamic framework that tracks entropy changes resulting from the continuous change (e.g., softening) of phonons as a function of temperature. A key finding of our study is a correlation between Poisson's ratio and the differential rate of phonon softening at finite frequencies, compared to lower frequencies relevant to elastic moduli measurements. We observe that elemental solids such as α -Be, diamond, Al, Cu, In, W, Au, and Pb, which span a wide range of Poisson's ratios and exhibit varying degrees of intrinsic anharmonicity, consistently follow this trend. When applied to α -U, α -Pu, and δ -Pu, our method reveals unusually large anharmonic phonon contributions at elevated temperatures across all three light actinide metals. These findings are attributed to the unique combination of enhanced covalency and softer elastic moduli inherent in the actinides, potentially influenced by their 5 f -electron bonding. Published by the American Physical Society 2025

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploration of Potential Superconducting Multi-Mode Cavity Architectures for Quantum Computing

This thesis describes the investigation of superconducting multi-mode cavity architecture for superconducting transmon-based quantum computing. The dissertation highlights useful features of radio-frequency cavities used for quantum computing, with a brief discussion on the advantages of superconducting cavities. In the subsequent section, the concept of transmon is introduced and the mechanism of coupling with a cavity is analyzed. Once the foundations of the topic are laid down, the thesis focuses on optimizing a multi-mode superconducting rf cavity design originally developed for high-energy physics applications. Such section articulates in two main parts, the first part concerning the "bare" cavity remake via finite-elements eigenmode simulations using a computer-aided design software called CST Studio Suite®. In the second part, a transmon qubit is physically inserted into the modified cavity to assess the qubit-cavity coupling of the new design. In evaluating said coupling, two distinct analysis methods are used, namely the black-box quantization method and the energy participation ratio method, both implemented using Ansys® High-Frequency Electromagnetic-Field Simulator, or HFSS™. Results from the two evaluations, compared together, show that the optimized design meets the requisites to be used for quantum computing purposes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Nonlinear modeling of the scaling law for the m/n = 3/2 error field penetration threshold

The scaling law for the m/n=3/2 error field (EF) penetration threshold is predicted numerically based on nonlinear single-fluid and two-fluid modeling using the TM1 code. The simulated penetration threshold of radial magnetic field b r at the plasma edge is scaled to the electron density n e , temperature T e , viscous time τ μ , toroidal field B t and the natural frequency ω in the form of b r /B t ∝n e αn T e αT τ μ αμ B t αB ω αω by scanning these parameters separately. Here, α n , α T , α μ , α B and α ω are the scaling coefficients on n e , T e , τ μ , B t and ω, respectively. Single-fluid modeling shows that the 3/2 EF threshold scales as b r /B t ∝n e 0.56 T e 0.6 τ μ -0.59 B t -1.15 ω, which is similar with the analytical scaling law in both the Rutherford and visco-resistive regimes. Yet, two-fluid modeling shows that the scaling law differs significantly in particular regarding the dependence on plasma rotation. In detail, the scaling coefficient α n on density decreases from 0.67 to 0.56 and α T on temperature decreases from 0.67 to 0.32, while α μ on viscous time is around -0.45 and α B on toroidal field decreases slightly from -1.15 to -1, when the ratio |ω E /ω *e | between plasma rotation frequency ω E and diamagnetic drift frequency ω *e varies from 0 to 10. Scans of the plasma rotation reveals that the penetration threshold linearly depends on the perpendicular rotation frequency (or natural frequency) ω ⊥e =ω E +ω *e , and there is a minimum in the required field amplitude when ω ⊥e 0. In addition, the enduring mystery of non-zero penetration threshold at zero plasma natural frequency in EF experiments is resolved by two-fluid simulations. We report that the very small island and smooth bifurcation in EF penetration near zero frequency is hard to detect in the experiment, leading to a finite penetration threshold within the capability of the experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗