Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditional neutrality”

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

Interfacial electroneutrality controls transport of asymmetric salts through charge-patterned mosaic membranes

Membranes that selectively enhance target solute permeation while rejecting competing species are essential for precision separations. This study introduces charge-patterned mosaic membranes (CMMs) that selectively transport divalent asymmetric salts by leveraging a net-neutral membrane–solution interface. This mechanism, dictated by the charge ratio of positive and negative domains on the membrane surface and the balance of cations and anions in the salt, is supported by analytical, numerical, and experimental results. Analytical solutions identified cationic domain coverages (f + ) of 33%, 50%, and 66% as optimal for the selective transport of +2:−1 salts, +1:−1 salts, and +1:−2 salts, respectively, under conditions where the pattern size (L) is significantly larger than the Debye length. Numerical simulations and experiments using CMMs with alternating charged-stripes inkjet-printed onto nanostructure copolymer substrates confirmed these findings. By varying stripe widths to control f + , pressure-driven filtration experiments demonstrated selective enrichment of MgCl 2 and K 2 SO 4 at the predicted f + values, with deviations from these values leading to salt rejection. These results highlight the pivotal role of a net-neutral interface in enabling asymmetric salt enrichment. This study positions CMMs as a versatile platform for tuning ion selectivity, addressing challenges in resource recovery, water treatment, and precision separations.

additive manufacturing↗
NASA NTRS2012

The Limit of Free Magnetic Energy in Active Regions

By measuring from active-region magnetograms a proxy of the free energy in the active region fs magnetic field, it has been found previously that (1) there is an abrupt upper limit to the free energy the field can hold that increases with the amount of magnetic field in the active region, the active region fs magnetic flux content, and (2) the free energy is usually near its limit when the field explodes in a CME/flare eruption. That is, explosive active regions are concentrated in a main-sequence path bordering the free-energy ]limit line in (flux content, free-energy proxy) phase space. Here, from measurement of Marshall Space Flight Center vector magnetograms, we find the magnetic condition that underlies the free ]energy limit and the accompanying main sequence of explosive active regions. Using a suitable free ]energy proxy measured from vector magnetograms of 44 active regions, we find that (1) in active regions at and near their free ]energy limit, the ratio of magnetic-shear free energy to the non ]free magnetic energy the potential field would have is approximately 1 in the core field, the field rooted along the neutral line, and (2) this ratio is progressively less in active regions progressively farther below their free ]energy limit. This shows that most active regions in which this core-field energy ratio is much less than 1 cannot be triggered to explode; as this ratio approaches 1, most active regions become capable of exploding; and when this ratio is 1 or greater, most active regions are compelled to explode. From these results we surmise the magnetic condition that determines the free ]energy limit is the ratio of the free magnetic energy to the non-free energy the active region fs field would have were it completely relaxed to its potential ]field configuration, and that this ratio is approximately 1 at the free-energy limit and in the main sequence of explosive active regions.

Moore, Ron↗
DOE OSTI2021

Concentrating Rare Earth Elements in Coal mine drainage Using Coal Combustion Products through Abandoned Mine Land Reclamation

Rare earth elements (REEs) (including scandium, yttrium and a group of 15 lanthanides) are often considered to be critical components in the productions of renewable energy hardware, electric vehicles, health care and military equipment, and consumer electronic products. The demand of REEs has been projected to be growing at an annual rate of 5-9% in the next 25 years. In 2011, the global demand of total rare earth oxides (REOs) was estimated to be approximately 105,000 tons, which is expected to grow to 210,000 tons by 2025. China overwhelmingly dominates the current worldwide rare earth productions but has strategically restricted its exports, causing significant instability for the global market. In response to the increasing demand for REEs and the supply dominance of China, identifying alternative sources of REEs has become a critical issue for the United States and other countries. Coal, coal ash, and coal mine drainage (CMD) are considered to be the alternative sources of REEs. In the U.S., high REE concentrations have been reported to be closely associated with coal deposits, including the Appalachian Basins. When surface and/or groundwater come in contact with geologic strata containing sulfide minerals exposed by coal mining, the accelerated oxidation of sulfide minerals in the presence of ferric iron and/or oxygen can produce sulfuric acid. The process promotes the weathering of REE-bearing rocks and minerals in the host geologic strata. Compared to average river water and seawater, the concentrations of REEs can be orders of magnitude higher in CMD. In this study, we demonstrated a trap-extract-precipitate (TEP) process that can effectively recover REEs from CMD. The three-stage TEP process uses alkaline industrial by-products to capture REEs from CMD and then applies an extraction/precipitation procedure to produce a feedstock that can be economically processed to produce marketable rare earth oxides. The alkaline industrial by-products tested in this study include the residual from a water softening process (DRWP sludge) and two types of stabilized flue gas desulfurization materials (sFGDs). sFGD material is a mixture of lime (CaO) and two coal combustion by-products, calcium sulfite FGD by-product and fly ash. The objectives of this study are to (1) validate the effectiveness and feasibility; (2) determine mechanisms controlling the rare earth recovery, (3) quantify the associated economic and environmental benefits, and (4) evaluate the full-scale application. To achieve these objectives, tasks to be carried out in this proposed project are organized into three phases. In the first phase, the research team collaborated with Ohio Department of Natural Resources, American Electric Power, The Wilds (a nonprofit wildlife conservation organization), and a private landowner to carry out field investigations aimed to screen and evaluate the seasonal changes of rare earths in the CMD discharges that have high recovery potentials. Next, the recovery of REEs from CMD was tested using a series of lab-scale column and batch tests under, respectively, percolation and completely mixed conditions. Results obtained from these lab-scale studies show that all three tested solids are very effective in retaining REEs. Over 98% of the CMD REEs that contacted the solids were captured before the solids exhausted their neutralization capacities. We also determined an extraction process using a non-acid, organic ligand extraction solution that can effectively remobilize the retained REEs from the spent solids (over 90%). The REE concentrate (>7.5 wt. % of total REEs) is then formed in an aeration process. The TEP process uses environmentally benign industrial by-products and a naturally-occurring organic ligand to mitigate CMD and recover REEs. Techno-economic analysis (TEA) and life-cycle assessment (LCA) was carried out in the third phase. The engineering-economic costs and net energy, net CO2 emissions, and water and other requirements were investigated to understand the economic and environmental implications of this process. This work uses mass and energy balances from laboratory-scale experiments to estimate the economic costs and environmental impacts. The results suggest that passive treatment systems that use DRWP sludge are preferred over those that use sFGD material, because of lower economic costs ($89,300/yr with a unit cost of $86/gT-REE vs. $89,800/yr, or $278/gT-REE) and improved environmental performance across all indicators from two different impact assessment methods. These differences are largely attributable to the larger capacity of DRWP sludge in the passive treatment application. We envision this TEP process can be integrated with abandoned mine land (AML) reclamation to create an approach that can add economic incentives for AML reclamation, remediate CMD discharge, and eliminate public safety hazards and threats to local environment and ecological systems posed by AMLs. It can restore lands and communities that are adversely impacted by legacy mining.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1981

Additional extensions to the NASCAP computer code, volume 3

The ION computer code is designed to calculate charge exchange ion densities, electric potentials, plasma temperatures, and current densities external to a neutralized ion engine in R-Z geometry. The present version assumes the beam ion current and density to be known and specified, and the neutralizing electrons to originate from a hot-wire ring surrounding the beam orifice. The plasma is treated as being resistive, with an electron relaxation time comparable to the plasma frequency. Together with the thermal and electrical boundary conditions described below and other straightforward engine parameters, these assumptions suffice to determine the required quantities. The ION code, written in ASCII FORTRAN for UNIVAC 1100 series computers, is designed to be run interactively, although it can also be run in batch mode. The input is free-format, and the output is mainly graphical, using the machine-independent graphics developed for the NASCAP code. The executive routine calls the code's major subroutines in user-specified order, and the code allows great latitude for restart and parameter change.

Mandell, M. J.↗
NASA NTRS1984

The effects of the diffuse radiation fields due to multiple scattering and thermal reradiation by dust on the dynamics and thermodynamics of a dusty cometary atmosphere

A self-consistent model of a dirty, clathrate cometary nucleus is extended to account for diffuse radiation fields caused by multiple scattering and thermal reradiation of the solar continuum by dust. The model is configured to fit conditions expected for the various spacecraft which will encounter Halley's comet at 0.89 AU. The atmosphere is assumed a chemically reactive dust-gas mixture in quasi-steady spherically symmetric expansion. The effect of electron-neutral ion collisions is accounted for, along with rate constants of the various species of clathrate ice particles and radiative transfer of the solar input in the UV, visible and near-IR intervals. The opacity of the circumnuclear dust is projected to be an order of magnitude greater than previous estimates, which severely impacts the potential visibility of the nucleus to the cameras of the flyby satellites.

Marconi, M. L.↗
NASA NTRS1986

Simulation of the thermospheric tides by use of the NCAR thermospheric general circulation model

Numerical calculations of the thermospheric tidal winds and temperatures at equinox are presented. The calculations were made using the National Center for Atmospheric Research Thermospheric General Circulation Model (TGCM), which includes the effects of viscosity, conductivity, diffusion, ion drag, winds, and temperature gradients. The thermospheric diurnal and semidiurnal tides are excited in situ by solar heating and by ion-neutral momentum coupling. The semidiurnal tidal calculations also include the effects of upward propagating waves generated by heating in the lower atmosphere. This semidiurnal propagating component is modeled by use of the classical tidal perturbations as lower boundary conditions. The model is tuned by adjusting the propagating tidal forcing term until calculated semidiurnal wind and temperature fields best approximate incoherent scatter observations. The tidal TGCM results are consistent with previous theoretical work and successfully reproduce high altitude temperature and meridional velocity data, but they give significantly lower magnitudes for velocities and temperatures near 160 km than are seen by observations.

Fensen, C. G.↗
NASA NTRS1989

Theoretical study of the seasonal behavior of the global ionosphere at solar maximum

The seasonal behavior of the global ionosphere was studied using a time-dependent three-dimensional physical model (developed by Shunk and his coworkers) of the ionosphere at altitudes between 120 and 800 km. This model accounts for field-aligned diffusion, cross-field electrodynamic drifts both the equatorial region and at high latitudes, interhemispheric flow, thermospheric winds, polar wind escape, energy-dependent chemical reactions, neutral composition changes, ion production due to solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow, and local heating and cooling processes. The model studies were carried out for both June and December solstice conditions at solar maximum and for low geomagnetic activity. The ionospheric features predicted by the model agreed qualitatively with the available measurements.

Sojka, J. J.↗
NASA NTRS1990

Properties of large scale plasma flow during the early stage of the plasmaspheric refilling

The objective is to better characterize the macroscopic properties of the interhemisphere plasma flow by solving a more complete set of hydrodynamic equations than that solved previously. Specifically, the ion continuity, momentum and energy equations were solved for the plasma flow along the closed magnetic field lines. During the initial stage of the supersonic outflow in the equatorial region, the ions cool substantially. Using the hydrodynamic model for the large-scale plasma flow, the dynamics of shocks was examined which form in the geomagnetic flux tubes during the early stages of refilling. These shocks are more like those forming in neutral gases than the electrostatic shocks driven by microinstabilities involving ion-ion interaction. Therefore, the shocks seen in the hydrodynamic model are termed as hydrodynamic shocks. Such shocks are generally unsteady and therefore the usual shock jump conditions given by Rankine-Hugoniot relations are not strictly applicable to them. The density, flow velocity and temperature structures associated with the shocks are examined for both asymmetrical and symmetrical flows. In the asymmetrical flow the outflow from one of two conjugate ionospheres is dominant. On the other hand, in the symmetrical case outflows from the two ionospheric sources are identical. Both cases are treated by a two-stream model. In the late type of flow, the early-time refilling shows a relaxation type of oscillation, which is driven by the large-scale interactions between the two identical streams. After this early stage, the resulting temperature structure shows some interesting features. In the equatorial region the streams are isothermal, but in the off-equatorial regions the streams have quite different temperatures, and also densities and flow velocities. The dense and slow stream is found to be warmer than the low-density fast stream. In the late stage of refilling, the temperature is found to steadily increase from the conjugate ionospheres towards the equator; the equatorial temperature is found to be as high as about 8000 K compared to the ionospheric temperature of 3600 K.

Singh, Nagendra↗
NASA NTRS1958

Stability of Cylindrical and Conical Shells of Circular Cross Section, with Simultaneous Action of Axial Compression and External Normal Pressure

We consider in this report the determination of the upper limit of critical loads in the case of simultaneous action of a compressive force, uniformly distributed over plane cross sections, and of isotropic external normal pressure on cylindrical or conical shells of circular cross section. As a starting point we use the differential equations for neutral equilibrium of conical shells which have been used for the solution of the problem of stability of conical shells under torsion and under axial compression; upon solution of the problem it is possible to satisfy all boundary conditions, in contrast to the report where no attention is paid to the fulfillment of the boundary conditions, and to the report where only part of the boundary conditions are satisfied by solution of the problem according to Galerkin's method. Approximate formulas are used for the determination of the critical external normal pressure with simultaneous action of longituninal compression. Let us note that the formulas suggested in reference 5 are not well founded and may lead, in a number of cases, to a substantial mistake in the magnitude of the critical load.

Mushtari, K. M.↗
NASA NTRS2005

NASA GRC and MSFC Space-Plasma Arc Testing Procedures

Tests of arcing and current collection in simulated space plasma conditions have been performed at the NASA Glenn Research Center (GRC) in Cleveland, Ohio, for over 30 years and at the Marshall Space Flight Center (MSFC) in Huntsville, Alabama, for almost as long. During this period, proper test conditions for accurate and meaningful space simulation have been worked out, comparisons with actual space performance in spaceflight tests and with real operational satellites have been made, and NASA has achieved our own internal standards for test protocols. It is the purpose of this paper to communicate the test conditions, test procedures, and types of analysis used at NASA GRC and MSFC to the space environmental testing community at large, to help with international space-plasma arcing-testing standardization. To be discussed are: 1.Neutral pressures, neutral gases, and vacuum chamber sizes. 2. Electron and ion densities, plasma uniformity, sample sizes, and Debuy lengths. 3. Biasing samples versus self-generated voltages. Floating samples versus grounded. 4. Power supplies and current limits. Isolation of samples from power supplies during arcs. 5. Arc circuits. Capacitance during biased arc-threshold tests. Capacitance during sustained arcing and damage tests. Arc detection. Prevention sustained discharges during testing. 6. Real array or structure samples versus idealized samples. 7. Validity of LEO tests for GEO samples. 8. Extracting arc threshold information from arc rate versus voltage tests. 9. Snapover and current collection at positive sample bias. Glows at positive bias. Kapon (R) pyrolisis. 10. Trigger arc thresholds. Sustained arc thresholds. Paschen discharge during sustained arcing. 11. Testing for Paschen discharge threshold. Testing for dielectric breakdown thresholds. Testing for tether arcing. 12. Testing in very dense plasmas (ie thruster plumes). 13. Arc mitigation strategies. Charging mitigation strategies. Models. 14. Analysis of test results. Finally, the necessity of testing will be emphasized, not to the exclusion of modeling, but as part of a complete strategy for determining when and if arcs will occur, and preventing them from occurring in space.

Ferguson, Dale C.↗
DOE OSTI2023

Analysing the effects of heating and gas puffing in Proto-MPEX helicon and auxiliary heated plasmas

The Material Plasma Exposure Experiment (MPEX) is being constructed at Oak Ridge National Laboratory to investigate critical fusion reactor issues, such as plasma–material interactions (PMIs) under reactor-relevant conditions and time scales. The linear device Proto-MPEX was used as a test bed to address anticipated research and development issues associated with heating scenarios and establish the physics basis for MPEX. The SOLPS-ITER code suite has been applied to understand plasma and neutral transport in Proto-MPEX and to increase confidence in predictive simulations for MPEX. Coupling between COMSOL and SOLPS is performed to implement a 2D electron heating profile of the helicon source. The simulations show reasonable agreement with the experimental data for plasma with helicon and auxiliary electron cyclotron heating (ECH). Both Bohm and constant diffusion (D$_\perp$: 0.5 m 2 s –1 and $\chi_\perp$: 1 m 2 s –1 ) simulations show similar levels of agreement with respect to the sparse experimental data available, assuming a few per cent impurity concentration. ECH significantly increases the target electron temperature, however, the target electron density is reduced compared to helicon-only heated plasmas due to an increase in flow velocity and radial losses. The simulations show that further increasing the ECH power results in an increase in the target electron density due to increased recycling flux and ionization. The results indicate that ECH significantly enhances the target heat flux, with ECH power of 50 kW increasing the target heat flux from 0.4 to 17 MW m –2 . It is found that a small amount of gas puffing (GP) near the target plate can further increase the target heat fluxes at the higher ECH power cases, but the target heat flux is reduced at higher GP conditions due to a significant reduction in electron temperature via radiation. ECH and GP scenarios can generate a higher target flux, facilitating improved PMI studies with more reactor-relevant plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1978

Specific spacecraft evaluation: Special report

Charged and neutral particle transport from an 8 cm mercury ion thruster to the surfaces of the P 80-1 spacecraft and to the Teal Ruby sensor and the ECOM-501 sensor of that spacecraft were investigated. Laboratory measurements and analyses were used to examine line-of-sight and nonline-of sight particle transport modes. The recirculation of Hg(+) ions in the magnetic field of the earth was analyzed for spacecraft velocity and Earth magnetic field vector configurations which are expected to occur in near Earth, circular, high inclination orbits. For these magnetic field and orbit conditions and for expected ion release distribution functions, in both angles and energies, the recirculation/re-interception of ions on spacecraft surfaces was evaluated. The refraction of weakly energetic ions in the electric fields of the thruster plasma plume and in the electric fields between this plasma plume and the material boundaries of the thruster, the thruster sputter shield, and the various spacecraft surfaces were examined. The neutral particle transport modes of interest were identified as sputtered metal atoms from the thruster beam shield. Results, conclusions, and future considerations are presented.

Sellen, J. M., Jr.↗
DOE OSTI2021

Low coordination number copper catalysts for electrochemical CO 2 methanation in a membrane electrode assembly

The electrochemical conversion of CO 2 to methane provides a means to store intermittent renewable electricity in the form of a carbon-neutral hydrocarbon fuel that benefits from an established global distribution network. The stability and selectivity of reported approaches reside below technoeconomic-related requirements. Membrane electrode assembly-based reactors offer a known path to stability; however, highly alkaline conditions on the cathode favour C-C coupling and multi-carbon products. In computational studies herein, we find that copper in a low coordination number favours methane even under highly alkaline conditions. Experimentally, we develop a carbon nanoparticle moderator strategy that confines a copper-complex catalyst when employed in a membrane electrode assembly. In-situ XAS measurements confirm that increased carbon nanoparticle loadings can reduce the metallic copper coordination number. At a copper coordination number of 4.2 we demonstrate a CO 2 -to-methane selectivity of 62%, a methane partial current density of 136 mA cm –2 , and > 110 hours of stable operation.

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

Self-consistent modeling of tokamak edge plasma transport with lithium sources

Magnetic confinement fusion devices require effective heat and particle exhaust solutions on the divertor plates to operate sustainably, especially under reactor-relevant conditions. Liquid lithium divertors have been proposed to address two major challenges: control of excessive heat flux to plasma-facing components through vapor shielding and minimization of core plasma contamination from impurities. The National Spherical Torus Experiment-Upgrade (NSTX-U) will explore lithium as a divertor material due to its potential to meet both objectives. We present a self-consistent coupling framework between the plasma boundary transport code UEDGE and the lithium wall transport code Wall–Li to evaluate the feasibility and operational limits of lithium-based divertors. The model aims to optimize lithium sourcing levels to prevent core plasma contamination via fuel dilution while ensuring divertor protection through vapor shielding. This integrated framework, applicable to any tokamak with lithium sources, dynamically adjusts lithium sourcing based on local plasma conditions and surface temperature. The coupled model is tested using NSTX-like geometry and plasma conditions to assess its performance and reliability. Wall–Li calculates lithium fluxes from plasma-facing components, incorporating physical sputtering, thermally enhanced sputtering, and evaporation driven by surface temperature and ion flux. These fluxes are reintroduced into UEDGE as neutral lithium atoms, enabling simulation of their transport and distribution within the plasma. UEDGE computes plasma and neutral transport, surface heat flux, and iteratively feeds this information back to Wall–Li. A small time step is employed to ensure numerical stability and convergence, enabling accurate simulations over typical tokamak discharge durations. This integrated modeling approach provides a robust tool for identifying operational regimes that balance effective lithium sourcing with minimal core plasma contamination, offering critical insights for optimizing lithium-based divertor systems in current and future fusion devices.

Magnetic confinement fusion↗
DOE OSTI2020

Enhanced energy confinement after series of pellets in Wendelstein 7-X

A series of ice pellets was injected into the advanced stellarator Wendelstein 7-X (W7-X). Although the pellets were small and slow, deep and efficient particle fueling could be observed experimentally. The most striking feature appearing after the injection of the pellets, however, was a transient increase in the energy confinement time. This transient phase resembled in several aspects modes of enhanced confinement after gas-puff or pellet injection, as observed in other fusion experiments. All experimental attempts, to prolong this phase, failed. In this paper, discharges are described that show the enhanced energy confinement, and some conditions are summarized which seem to be essential in order to generate it. The focus here is on deep particle fueling by pellets, and shaping of the density profiles during and after the series of pellets. During this time, neutral gas particle re-fueling at the plasma edge is reduced, while density profile peaking and low impurity radiation losses are present.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The diurnal susceptibility of subtropical clouds to aerosols

The diurnal susceptibility of clouds and their radiative properties to aerosols is examined during their Lagrangian transition from subtropical stratocumulus to shallow cumulus regimes. Using large-eddy simulations, we analyze the six-day evolution of an air mass along a 3800 km observed trajectory from the coast of Peru toward the equator. Pristine and polluted scenarios are simulated with forcing imposed from weather reanalysis. The polluted scenario exhibits stronger diurnal variations in cloud water, cloud fraction, and albedo, with enhanced nighttime entrainment and suppressed precipitation. The overall response of cloud properties and outgoing shortwave radiation to droplet number concentration follows a distinct diurnal pattern: strong positive cloud adjustments dominate at night and in the morning, while weak negative adjustments prevail in the afternoon. This cycle is driven by the competition between precipitation suppression, which enhances cloud water and coverage, and entrainment drying, which depletes them. In polluted conditions, enhanced entrainment leads to a deeper and more decoupled boundary layer that cannot be sustained by surface fluxes in the afternoon, resulting in negative cloud adjustments. The enhanced entrainment rate under polluted conditions is caused by the reduced sedimentation of cloud and precipitation water from the entrainment zone. While the Twomey effect dominates the diurnal average albedo response, the diurnal variation in the competing cloud adjustments lead to a near-neutral net susceptibility in the afternoon, highlighting the critical role of diurnally varying processes in aerosol-cloud interactions.

Geosciences↗
NASA NTRS1979

Nighttime ionization by energetic particles at Wallops Island in the altitude region 120 to 200 km

Five Nike Apache rockets, each including an energetic particle spectrometer and an electron density-electron temperature experiment, have been launched from Wallops Island (L = 2.6) near midnight under varying geomagnetic conditions. On the most recent of these (5 January 1978) an additional spectrometer with a broom magnet, and a 391.4 nm photometer were flown. The data from this flight indicate that the energetic particle flux consists predominantly of protons, neutral hydrogen and possibly other energetic nuclei. The energy spectrum becomes much softer and the flux more intense with increasing Kp for 10-100 keV. The pitch angle distribution at 180 km is asymmetrical with a peak at 90 deg indicating that the majority of particles are near their mirroring altitude. Ionization rates are calculated based on the measured energy spectrum and mirror height distribution. The resulting ionization rate profile is found to be nearly constant with altitude in the region 120 to 200 km. The measured energetic particle flux and calculated ionization rate from the five flights are found to vary with magnetic activity (based on the Kp and Dst indexes) in the same way as the independently derived ionization rates deduced from the electron density profile.

Voss, H. D.↗
NASA NTRS1979

Harmonics of 60 Hz in power systems caused by geomagnetic disturbances

Simultaneous VLF/ULF observations carried out near Winnipeg, Manitoba show that geomagnetic disturbances control the behavior of harmonics of 60 Hz man-made electric power. The harmonics of 60 Hz detected by the VLF receiver are at multiples of 180 Hz, indicating that they originated from a 3 phase ac power system. Under geomagnetically quiet conditions, only odd harmonics of 70 Hz were detected. In disturbed conditions, both odd and even harmonics were excited. The strength of each harmonic changed concurrently with geomagnetic pulsation (ULF) activity. These findings seem to indicate that a portion of telluric currents shunted into the power line system through the neutrals of the Y-connected transformers give rise to a dc bias to the transformer core materials and that it distorts their hysteresis loops, activating harmonics of 60 Hz power. A mathematical proof is given that a hysteresis loop having a point of symmetry generates odd harmonics only, whereas loops lacking in point-symmetry generally give rise to both odd and even harmonics. A general formula was obtained to calculate the strength of each harmonic based on the shape of the hysteresis loop.

Hayashi, K.↗