"OPTIMAL SIMULATION TECHNIQUES FOR DISTRIBUTED ENERGY STORE RAILGUNS WITH SOLID STATE SWITCHES"

Similar documents
Electromagnetic Railgun

Research on High Power Railguns at the Naval Research Laboratory

DEVELOPMENT OF AN ULTRA-COMPACT EXPLOSIVELY DRIVEN MAGNETIC FLUX COMPRESSION GENERATOR SYSTEM

PULSED POWER SWITCHING OF 4H-SIC VERTICAL D-MOSFET AND DEVICE CHARACTERIZATION

ANALYSIS OF A PULSED CORONA CIRCUIT

COM DEV AIS Initiative. TEXAS II Meeting September 03, 2008 Ian D Souza

David L. Lockwood. Ralph I. McNall Jr., Richard F. Whitbeck Thermal Technology Laboratory, Inc., Buffalo, N.Y.

An experimental system was constructed in which

ANALYSIS OF SWITCH PERFORMANCE ON THE MERCURY PULSED- POWER GENERATOR *

FLASH X-RAY (FXR) ACCELERATOR OPTIMIZATION BEAM-INDUCED VOLTAGE SIMULATION AND TDR MEASUREMENTS *

INVESTIGATION OF A HIGH VOLTAGE, HIGH FREQUENCY POWER CONDITIONING SYSTEM FOR USE WITH FLUX COMPRESSION GENERATORS

Strategic Technical Baselines for UK Nuclear Clean-up Programmes. Presented by Brian Ensor Strategy and Engineering Manager NDA

Investigation of a Forward Looking Conformal Broadband Antenna for Airborne Wide Area Surveillance

A NEW BROADBAND PULSED HIGH VOLTAGE MONITOR *

Design, Construction, and Testing of an Inductive Pulsed-Power Supply for a Small Railgun

IB2-1 HIGH AVERAGE POWER TESTS OF A CROSSED-FIELD CLOSING SWITCH>:< Robin J. Harvey and Robert W. Holly

PULSED BREAKDOWN CHARACTERISTICS OF HELIUM IN PARTIAL VACUUM IN KHZ RANGE

EFFECT OF TRANSFORMER LEAKAGE INDUCTANCE ON THE THREE PHASE CAPACITIVE INPUT RECTIFIER

TRANSMISSION LINE AND ELECTROMAGNETIC MODELS OF THE MYKONOS-2 ACCELERATOR*

DEVELOPMENT OF STITCH SUPER-GTOS FOR PULSED POWER

Modeling of Ionospheric Refraction of UHF Radar Signals at High Latitudes

Technology Maturation Planning for the Autonomous Approach and Landing Capability (AALC) Program

A RENEWED SPIRIT OF DISCOVERY

Robotics and Artificial Intelligence. Rodney Brooks Director, MIT Computer Science and Artificial Intelligence Laboratory CTO, irobot Corp

Signal Processing Architectures for Ultra-Wideband Wide-Angle Synthetic Aperture Radar Applications

DESIGN OF A 16 kv, 100 ka, 2Hz POWER SUPPLY FOR HIGH-FIELD, REPETITIVELY PULSED, SPLIT-PAIR MAGNETS

Solar Radar Experiments

Simulating the Difference between a DES and a Simple Railgun using SPICE

0.18 μm CMOS Fully Differential CTIA for a 32x16 ROIC for 3D Ladar Imaging Systems

Student Independent Research Project : Evaluation of Thermal Voltage Converters Low-Frequency Errors

FLASH X-RAY (FXR) ACCELERATOR OPTIMIZATION INJECTOR VOLTAGE-VARIATION COMPENSATION VIA BEAM-INDUCED GAP VOLTAGE *

Innovative 3D Visualization of Electro-optic Data for MCM

PHASING CAPABILITY. Abstract ARRAY. level. up to. to 12 GW. device s outpu antenna array. Electric Mode. same physical dimensions.

U.S. Army Training and Doctrine Command (TRADOC) Virtual World Project

9.4 A HIGH CURRENT PULSER FOR EXPERIMENT 11225, "NEUTRINO ELECTRON ELASTIC SCATTERING" C. Dalton, G. Krausse, and J. Sarjeant

Durable Aircraft. February 7, 2011

SILICON CARBIDE FOR NEXT GENERATION VEHICULAR POWER CONVERTERS. John Kajs SAIC August UNCLASSIFIED: Dist A. Approved for public release

Development of a charged-particle accumulator using an RF confinement method FA

Parametric Analyses Using a Computational System Model of an Electromagnetic Railgun

OPTICAL EMISSION CHARACTERISTICS OF HELIUM BREAKDOWN AT PARTIAL VACUUM FOR POINT TO PLANE GEOMETRY

A Comparison of Two Computational Technologies for Digital Pulse Compression

8.2. Washington, D. C delivered 65 kj into a matched load with 63 ns FWHM. Peak power was about 1 TW.

Report Documentation Page

Synthetic Behavior for Small Unit Infantry: Basic Situational Awareness Infrastructure

Rump Session: Advanced Silicon Technology Foundry Access Options for DoD Research. Prof. Ken Shepard. Columbia University

August 9, Attached please find the progress report for ONR Contract N C-0230 for the period of January 20, 2015 to April 19, 2015.

NEURAL NETWORKS IN ANTENNA ENGINEERING BEYOND BLACK-BOX MODELING

REPORT DOCUMENTATION PAGE. A peer-to-peer non-line-of-sight localization system scheme in GPS-denied scenarios. Dr.

Adaptive CFAR Performance Prediction in an Uncertain Environment

RAVEN, A 5 kj, 1.5 MV REPETITIVE PULSER* G. J. Rohwein Sandia National Laboratories Albuquerque, New Mexico 87185

HIGH TEMPERATURE (250 C) SIC POWER MODULE FOR MILITARY HYBRID ELECTRICAL VEHICLE APPLICATIONS

Department of Energy Technology Readiness Assessments Process Guide and Training Plan

REGULATED CAPACITOR CHARGING CIRCUIT USING A HIGH REACTANCE TRANSFORMER 1

Presentation to TEXAS II

FAA Research and Development Efforts in SHM

UNCLASSIFIED UNCLASSIFIED 1

EFFECTS OF ELECTROMAGNETIC PULSES ON A MULTILAYERED SYSTEM

The Energy Spectrum of Accelerated Electrons from Waveplasma Interactions in the Ionosphere

MATLAB Algorithms for Rapid Detection and Embedding of Palindrome and Emordnilap Electronic Watermarks in Simulated Chemical and Biological Image Data

Report Documentation Page

Ocean Acoustics and Signal Processing for Robust Detection and Estimation

IREAP. MURI 2001 Review. John Rodgers, T. M. Firestone,V. L. Granatstein, M. Walter

Acoustic Monitoring of Flow Through the Strait of Gibraltar: Data Analysis and Interpretation

CFDTD Solution For Large Waveguide Slot Arrays

Wavelet Shrinkage and Denoising. Brian Dadson & Lynette Obiero Summer 2009 Undergraduate Research Supported by NSF through MAA

Coherent distributed radar for highresolution

INTEGRATIVE MIGRATORY BIRD MANAGEMENT ON MILITARY BASES: THE ROLE OF RADAR ORNITHOLOGY

ACCELERATOR FAST KICKER R&D WITH ULTRA COMPACT 50MVA NANO-SECOND FID PULSE GENERATOR

SA Joint USN/USMC Spectrum Conference. Gerry Fitzgerald. Organization: G036 Project: 0710V250-A1

Future Trends of Software Technology and Applications: Software Architecture

Radar Detection of Marine Mammals

Fuse and Load Testing With Mid-Sized, High Energy Density Flux Compression Generators

Mathematics, Information, and Life Sciences

Improving the Detection of Near Earth Objects for Ground Based Telescopes

Hybrid QR Factorization Algorithm for High Performance Computing Architectures. Peter Vouras Naval Research Laboratory Radar Division

N C-0002 P13003-BBN. $475,359 (Base) $440,469 $277,858

Advances in SiC Power Technology

Underwater Intelligent Sensor Protection System

AFRL-VA-WP-TP

GLOBAL POSITIONING SYSTEM SHIPBORNE REFERENCE SYSTEM

Thermal Simulation of a Silicon Carbide (SiC) Insulated-Gate Bipolar Transistor (IGBT) in Continuous Switching Mode

Design of Synchronization Sequences in a MIMO Demonstration System 1

REPORT DOCUMENTATION PAGE

Army Acoustics Needs

Active Denial Array. Directed Energy. Technology, Modeling, and Assessment

3. Faster, Better, Cheaper The Fallacy of MBSE?

14. Model Based Systems Engineering: Issues of application to Soft Systems

David Siegel Masters Student University of Cincinnati. IAB 17, May 5 7, 2009 Ford & UM

Measurement of Ocean Spatial Coherence by Spaceborne Synthetic Aperture Radar

Loop-Dipole Antenna Modeling using the FEKO code

Operational Domain Systems Engineering

Bistatic Underwater Optical Imaging Using AUVs

ULTRA FAST, HIGH REP RATE, HIGH VOLTAGE SPARK GAP PULSER

Key Issues in Modulating Retroreflector Technology

EVALUATION OF RESISTORS FOR TRANSIENT HIGH-VOLTAGE APPLICATIONS

W. L. Bird University of Texas, Austin, Center for Electro-Mechanics, Taylor Hall 167, Austin, TX 78712

Modeling Antennas on Automobiles in the VHF and UHF Frequency Bands, Comparisons of Predictions and Measurements

L. B. Gordon Space Power Institute 231 Leach Center Auburn University, Alabama 36849

Two-Way Time Transfer Modem

REPORT DOCUMENTATION PAGE

Transcription:

"OPTIMAL SIMULATION TECHNIQUES FOR DISTRIBUTED ENERGY STORE RAILGUNS WITH SOLID STATE SWITCHES" James B. Cornette USAF Wright Laboratory WL/MNMW c/o Institute for Advanced Technology The University of Texas at Austin 4030-2 West Braker Lane Austin, Texas 78759-5329 Dr. Richard A. Marshall Institute for Advanced Technology The University of Texas at Austin 4030-2 West Braker Lane Austin, Texas 78759-5329 ABSTRACT The objective of this paper is to present an optimal design methodology to determine the best firing strategy, energy store sizing, energy store spacing and maximum system efficiency for a Distributed Energy Store (DES) railgun. System simulations/designs will be based on the assumption that switching of the energy storage units will be accomplished using solid state devices. Candidate semiconductor technologies are promising in relation to solving the high energy, low weight requirements of a railgun system and other pulsed power systems requiring high energy, compact switching. A simulation code has been developed and used to produce non-dimensional data files that are then scaled to physical railgun values based on input parameters. Energy stores are assumed to be capacitive in nature with diodes to prevent negative currents and crowbar diodes to prevent voltage reversal of the capacitors. The main thrust of this simulation effort is to produce a DES design that optimizes the efficiency of the conversion of stored electrical energy to projectile kinetic energy, while also considering the abilities of near term solid state switching devices. INTRODUCTION Impetus for this work is derived from the Focused Technology Program (FfP) goals for launching of a projectile with a railgun[ll.' Under the FTP program, a projectile is to be launched under the following parameter conditions: Projectile Mass: Exit Velocity: Barrel Length: 4.44 kg 3.0 km/s 7.0m The FTP goals also state a barrel lifetime of 1000 shots and a multi-shot mission profile is given. For this work we will only consider the single shot mission. Many other parameters are required to simulate or design a system that will meet the FTP goals. The selection of those parameters, the type of power supply, the methodology for simulating the system, and the results of the simulation and conversion program that outputs physical railgun parameters are discussed. DES SYSTEM 1285

Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number. 1. REPORT DATE JUL 1995 2. REPORT TYPE N/A 3. DATES COVERED - 4. TITLE AND SUBTITLE Optimal Simulation Techniques For Distributed Energy Store Railguns With Solid State Switches 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR(S) 5d. PROJECT NUMBER 5e. TASK NUMBER 5f. WORK UNIT NUMBER 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) USAF Wright Laboratory WL/MNMW c/o Institute for Advanced Technology The University of Texas at Austin 4030-2 West Braker Lane Austin, Texas 78759-5329 8. PERFORMING ORGANIZATION REPORT NUMBER 9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR S ACRONYM(S) 12. DISTRIBUTION/AVAILABILITY STATEMENT Approved for public release, distribution unlimited 11. SPONSOR/MONITOR S REPORT NUMBER(S) 13. SUPPLEMENTARY NOTES See also ADM002371. 2013 IEEE Pulsed Power Conference, Digest of Technical Papers 1976-2013, and Abstracts of the 2013 IEEE International Conference on Plasma Science. Held in San Francisco, CA on 16-21 June 2013. U.S. Government or Federal Purpose Rights License. 14. ABSTRACT The objective of this paper is to present an optimal design methodology to determine the best firing strategy, energy store sizing, energy store spacing and maximum system efficiency for a Distributed Energy Store (DES) railgun. System simulations/designs will be based on the assumption that switching of the energy storage units will be accomplished using solid state devices. Candidate semiconductor technologies are promising in relation to solving the high energy, low weight requirements of a railgun system and other pulsed power systems requiring high energy, compact switching. A simulation code has been developed and used to produce non-dimensional data files that are then scaled to physical railgun values based on input parameters. Energy stores are assumed to be capacitive in nature with diodes to prevent negative currents and crowbar diodes to prevent voltage reversal of the capacitors. The main thrust of this simulation effort is to produce a DES design that optimizes the efficiency of the conversion of stored electrical energy to projectile kinetic energy, while also considering the abilities of near term solid state switching devices. 15. SUBJECT TERMS 16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF ABSTRACT SAR a. REPORT unclassified b. ABSTRACT unclassified c. THIS PAGE unclassified 18. NUMBER OF PAGES 6 19a. NAME OF RESPONSIBLE PERSON

Standard Form 298 (Rev. 8-98) Prescribed by ANSI Std Z39-18

A DES system for a railgun is one that incorporates the use of many power supplies along the barrel instead of a single power supply placed at the breech, as in most conventional systems. One benefit of the DES configuration is that the barrel current is not flowing from the breech to the projectile at all times, as in a breech fed railgun. An introduction to the circuit equations that are pertinent to railgun operation are presented by 1. P. Barber[2) and J. B. Comette[3l, and the equations describing DES operation are presented well by J. V. Parker and R. A. Marshall[4,5]. Power Supplies Power supplies for this work are assumed to be capacitor banks. Capacitors have proven high current capability, availability, and operating voltages that match well to railgun systems. However, the simulation can be modified to model other types of energy storage devices. The capacitor banks are assumed to have crowbar diodes to prevent negative voltages from being impressed onto the capacitors and a diode in the forward current leg (see Figure 1) of each power supply to prevent negative currents from flowing into the capacitors. This forward leg diode precludes individual banks from discharging into the other capacitor banks. A single capacitor bank connected to the breech of a railgun is depicted in Figure 1. LBank R Bank L' + v CAP crowbar diode -1 diode.. I BARREL + v A Capacitors Figure 1. Single Capacitor Bank, Breech Connected to an Railgun. The breech voltage of the barrel is represented by VB, and the armature voltage of the projectile is VA. The breech voltage is governed by the following equation and the armature voltage is assumed constant for this work. The type of armature that each user wishes to simulate will cause the value of VA to vary tremendously and is up to the user to determine. The ability to interactively alter VA prior to each run is an important capability of this simulation. V =IR' V [ alai 1 al ax] B X+ A + X ax at + ax at Assuming L'=aLtax, and L' & R' are independent of projectile position: a I V 8 =IR'x+VA +L'x Jt +IL'v 1286

In the DES case, each bank will see a different VB depending on their connection point along the length of the barrel. For this simulation each capacitor stage has a series inductance and resistance that can be set at the beginning of each run and banks are spaced evenly down the length of the barrel as discussed in the next section. In this case, the FfP goals set the barrel length to 7.0 m. From the equations of motion and Maxwell's equations, the force exerted on the projectile due to the interaction of the current density vector and the magnetic flux density vector result in the commonly known Lorentz force. After determining the magnetic energy in terms if current (I) and position (x) and differentiating with respect to x we have the electric force on the projectile: and assuming constant I for small ax: F = _!_[/ 2 al + 2/L ai J e 2 ax ax Which is iri terms of the barrel inductance gradient (L'=aL;ax in J.l.H/m) and the current (1). This force is the basis for determining the optimum arrangement of power supplies along the barrel and constant force on the projectile, L' and I must be constant. The L' is assumed constant for this study and equal to 0.5 J.l.H/m. The power supplies are assumed to store equal energy and are placed at equal distances down the barrel. In order to transfer the energy stored in each capacitor bank most efficiently to the projectile as the projectile's velocity increases, the open circuit voltage and capacitance of the banks are allowed to vary. By changing the time constant of the banks as the projectile moves down the barrel, the energy stages toward the muzzle discharge faster than those closer to the breech. Allowing the capacitance and charge voltage of the energy stages to change is the most efficient method for the use of capacitor banks in a DES system (see figure 2), however, it is not as practical to build in reality. Therefore, the code has been altered to allow the user to choose whether or not they would like equal capacitance and voltage per stage as in a more realistic system. One design procedure would be to determine the optimum system using varying voltage and capacitance, then run the same system with equal voltage and capacitance. After several iterations, one can arrive at an approximation of the optimum system with a practical capacitor bank. The loss parameters of armature voltage and rail resistance are interactively set by the user at the beginning of each run. The armature voltage is constant and the rail resistance is the gradient in 0/m for one rail. SIMULATIONS The DES circuit solver used here is based on the work of J. V. Parker[Sl. The circuit equations are solved in dimensionless form. Once an optimized solution is determined, as related to the users input parameters, a physical railgun design must be translated from the dimensionless data. The original code has been modified as stated above to allow equal voltage and capacitance per stage, increasing the energy storage of any of the stages, and increasing the number of energy stages to 100. 1287

A second code has been written that takes the non-dimensional data from the FORTRAN circuit solving simulation and then produces a physical railgun data set, given certain input parameters. This conversion routine is a MATLAB based code that is capable of up to 100 individual energy stages. Output from the MA TLAB code shows the user how efficiently the energy was transferred to the projectile and various other parameters. A summary of those parameters are presented below: Number of stages (N) Circuit input parameters Energy stored Energy transferred to projectile System energy losses Energy losses per stage Initial stage volts (I to N) Stage R, L, C's, etc. Maximum gun current Maximum gun di/dt Maximum projectile acceleration Maximum bore pressure Maximum back voltage (L'Iv) Maximum stage currents (1 to N) Maximum stage di/dt's (1 to N) Many other parameters can be displayed and plotted for the user. EXAMPLE SIMULATION RUN As an example of the input, output and data plotting capabilities of the DES simulation and conversion program, a single run utilizing the ability to alter the energy storage of any stage and using the more efficient variable capacitance and voltage for each stage is summarized with the following results. below: Example input parameters, required before running the simulation, for the FfP goals are given Number of Stages = 28 Bore size (square) IOOmm Barrel Length = 7.0m Mass of projectile 4.44 kg Injection Velocity 0.098 km/s Exit Velocity 3.0 km/s L' = 0.5 f.lh/m Armature Voltage = 9.577 v R' = 122.61 ~/m Stage #1 Energy = 2.97 * Ei,i=1,28 An output summary for the FTP run is given below: Total Stored Energy = 28.576 MJ Max. Stage Current = 3.998 MA Exit Energy = 19.98 MJ Max. Stage di/dt = 37.981 ka/jls Exit Time = 4.759 ms Max. back Voltage = 5.413 kv Gun losses = 0.291 MJ Stage #1 C = 246.98 mf Stage Losses = 2.116 MJ Stage#l L = 0.330 J.t.H Energy Efficiency = 69.92% Stage#l R = 49.95 ~ Max. Bore Pressure = 57.95 kpsi Other Stage C's 2.9 to 41 mf Max. Acceleration = 91.80 kgee Other Stage L's = 0.980 JlH Max. Gun Current = 3.998 MA Other Stage R's = 148.35~ Exit Current = 3.441 MA When increasing the energy storage of any bank, it is assumed that equal energy storage banks are connected in parallel (i.e., Increase the energy storage of stage#l by 3.0, Cl =3*Cold, Ll=Lold/3, Rl=Rold/3 ). A plot of the total gun current and the individual stage currents versus real time is presented in Figure 2. Notice that the first energy stage current is much higher than the others. This simulation run is an example of increasing the energy storage of the first stage in order to reach maximum gun current 1288

faster. In this example the energy storage of the first stage is approximately three times that of the other stages. Gun current (MA) 4 3.5 3,..., ~ 2.5 2... 1.5 a 1 5 u 0.5 0 Total GJJ11 Current (MA) -- Stage Gurrents (MA) -0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 Time (ms) Figure 2. Total Gun Current and Stage Currents versus Time Projectile acceleration in kgee is shown if Figure 3. If considering solid state devices for the closing switches for such a DES system, the individual stage maximum di/dt's are an important factor (among many others) that will determine if they can successfully perform the switching. Therefore, an example of individual stage di/dt's in ka/j.ls are shown in Figure 4. The maximum di/dt capability of potential solid state devices today is approximately 20 ka/j.l.s. 100r-----r-----r-----r-----r-----r-----r-----r-----r-----~----, j. c:.q I 30 20 10 0 0~---0~.5~--~1~--~1~.~5----2~--~2~.5~--~3~--~3~.~5----~4----~4~.5~--~5 Time (ms) Figure 3. Projectile Acceleration 1289

40-35~ 30 ~~ 25!I <n : i g 20 l e: l ~15n ~ i! 10 ; I! i 5!\ Jl ' i\~.' 1... Stage di/dts (ka/micro-s) ' ~ ~ h ~1\ f\ 1\ \ ; ' I \ ~1: \ ~~ \\ l l ; \ \ : ;! \ l ; I \: \ \ \ : i \ : \i : i \ i \ i 1 l ~ ~ 1\ i\ ~\ \ H :: n I 'i q \ q :! \. i ~...-v.._.. VUwe/ ;,t ~VV\:J'[;T:I<~.fl -5~--~ ~--~ ~--~~--~----~----~ ----~ ----~ -------~~ v_~v--~ 0 0.5 1.5 2 2.5 3 3.5 4 4.5 5 Time (ms) Figure 4. Stage dildt's SUMMARY A practical DES railgun simulation code has been developed and used to complete first order system design for the FfP requirements. The methodology assumes equal energy storage in each stage, equal distance between stages, equal resistance and inductance per stage, and varying capacitance and charge voltage in each stage. Provisions have been made for the user to compare a system based on the above stated assumptions to one that has equal capacitance and voltage in each stage. The energy storage of any stage can also be adjusted to reach the maximum driving current faster than would normally occur in a DES system of equal energy stages. Future work will be directed towards adding solid state device characteristics to the code to enable design of switching schemes for each stage. Once complete, the long term goal is to construct an easy to use interactive interface application for users to design a complete DES system according to the type of energy storage units, switching devices, and barrel parameters desired. REFERENCES [I] H. D. Fair, "Applications of Electric Launch Systems", IEEE Transactions on Magnetics, Vol. 29, No.1, January 1993. [2] J. P. Barber, "The Acceleration of Macroparticles and Hypervelocity Electromagnetic Accelerator", Ph. D. dissertation, The Australian National University, Canberra, 1972. [3) J. B. Cornette, "Demonstration of a Battery Charged Capacitor Pulsed Power System for Rapid Fire Electromagnetic Launcher Experiments", Master's Thesis, August 1990, Wright Laboratory Armament Directorate Technical Report, WL-TR-92-7033. [4) R. A. Marshall, J.P. Barber, "The 10 km/s, 10 kg Railgun", IEEE Transactions on Magnetics, Vol. 27, No. 1, January 1991. [5] J. V. Parker, "Electromagnetic Projectile Acceleration Using Distributed Energy Stores", Journal Applied Physics, 53(10), October 1982. 1290