Engineering Papers⌕ Search

Engineering topics

Koshak, William J.

Publications and source records attributed to Koshak, William J..

65 records · Page 4

A Multipole Expansion Method for Analyzing Lightning Field Changes

Changes in the surface electric field are frequently used to infer the locations and magnitudes of lightning-caused changes in thundercloud charge distributions. The traditional procedure is to assume that the charges that are effectively deposited by the flash can be modeled either as a single point charge (the Q-model) or a point dipole (the P-model). The Q-model has 4 unknown parameters and provides a good description of many cloud-to-ground (CG) flashes. The P-model has 6 unknown parameters and describes many intracloud (IC) discharges. In this paper, we introduce a new analysis method that assumes that the change in the cloud charge can be described by a truncated multipole expansion, i.e., there are both monopole and dipole terms in the unknown source distribution, and both terms are applied simultaneously. This method can be used to analyze CG flashes that are accompanied by large changes in the cloud dipole moment and complex IC discharges. If there is enough information content in the measurements, the model can also be generalized to include quadrupole and higher order terms. The parameters of the charge moments are determined using a 3-dimensional grid search in combination with a linear inversion, and because of this, local minima in the error function and the associated solution ambiguities are avoided. The multipole method has been tested on computer simulated sources and on natural lightning at the NASA Kennedy Space Center and USAF Eastern Range.

Koshak, William J.↗

Electro-Optic Lightning Detector

Electric field measurements are fundamental to the study of thunderstorm electrification, thundercloud charge structure, and the determination of the locations and magnitudes of charges deposited by lightning. Continuous field observations can also be used to warn of impending electrical hazards. For example, the USAF Eastern Range (ER) and NASA Kennedy Space Center (KSC) in Florida currently operate a ground-based network of electric field mill sensors to warn against lightning hazards to space vehicle operations/launches. The sensors provide continuous recordings of the ambient field. Others investigators have employed flat-plate electric field antennas to detect changes In the ambient field due to lightning. In each approach, electronic circuitry is used to directly detect and amplify the effects of the ambient field on an exposed metal conductor (antenna plate); in the case of continuous field recordings, the antenna plate is alternately shielded and unshielded by a grounded conductor. In this work effort, an alternate optical method for detecting lightning-caused electric field changes is Introduced. The primary component in the detector is an anisotropic electro-optic crystal of potassium di-hydrogen phosphate (chemically written as KH2PO4 (KDP)). When a voltage Is placed across the electro-optic crystal, the refractive Indices of the crystal change. This change alters the polarization state of a laser light beam that is passed down the crystal optic axis. With suitable application of vertical and horizontal polarizers, a light transmission measurement is related to the applied crystal voltage (which in turn Is related to the lightning caused electric field change). During the past two years, all critical optical components were procured, assembled, and aligned. An optical housing, calibration set-up, and data acquisition system was integrated for breadboard testing. The sensor was deployed at NASA Marshall Space Flight Center (MSFC) in the summer of 1998 to collect storm data. Because solid-state technology is used, future designs of the sensor will be significantly scaled down In physical dimension and weight compared to the present optical breadboard prototype. The use of fiber optics would also provide significant practical improvements.

Koshak, William J.↗

Calculation of the Bulk Electromagnetic Properties of Thunderclouds Using a Two-Space Scattering Formalism

A two-space scatterer formalism and equivalent medium approach of Twersky are used to obtain formulae for the bulk electromagnetic parameters of a thundercloud medium that is illuminated by lightning optical emissions. A modified WKB approximation is applied to derive the two-space scattering amplitude. Optical constants of dry air, moist atmosphere, and water,.are inserted into the formulae to generate numerical values via iteration. The two-space scatterer formalism results are close to those obtained from free- or sin-le-space formalisms for the dilute case. The numerical values of the bulk parameters are required to successfully transform the scattering problem to that of an equivalent obstacle excited by an incident wave traveling in K - space but radiating in k(sub 1) -space.

Phanord, Dieudonne D.↗

On the Retrieval of Lightning Radio Sources from Time-of-Arrival Data

We examine the problem of retrieving three-dimensional lightning locations from radio frequency Time-Of-Arrival (TOA) measurements. Arbitrary antenna locations are considered. By judiciously differencing measurements that are related to the location of the antennas and their excitation times, the problem is converted from the initial spherical nonlinear form to a system of linear equations. In the linear formalism, the source location and time-of-occurrence is viewed geometrically as an intersection of hyperplanes in the four-dimensional Minkowski space (x,y,z,t). The linear equations are solved to obtain explicit analytic expressions for the location and time variables. Retrieval errors are not interpreted with conventional Geometrical Dilution of Precision (GDOP) arguments as discussed by Holmes and Reedy (1951), but with more recent inversion analyses considered by Twomey (1977). Measurement errors are propagated analytically so that the specific effect of these errors on the solution is clarified. The sensitivity of the solution on the number of antennas used, antenna network geometry, source position, and measurement differencing schemes are discussed in terms of the eigenvalues of the linear system.

Koshak, William J.↗

Electro-optic Lightning Detector

The design, alignment, calibration, and field deployment of a solid-state lightning detector is described. The primary sensing component of the detector is a potassium dihydrogen phosphate (KDP) electro-optic crystal that is attached in series to a flat plate aluminum antenna; the antenna is exposed to the ambient thundercloud electric field. A semiconductor laser diode (lambda = 685 nm), polarizing optics, and the crystal are arranged in a Pockels cell configuration. Lightning-caused electric field changes are related to small changes in the transmission of laser light through the optical cell. Several hundred lightning electric field change excursions were recorded during five thunderstorms that occurred in the summer of 1998 at the NASA Marshall Space Flight Center (MSFC) in northern Alabama.

Koshak, William J.↗

Aircraft electric field measurements: Calibration and ambient field retrieval

An aircraft locally distorts the ambient thundercloud electric field. In order to determine the field in the absence of the aircraft, an aircraft calibration is required. In this work a matrix inversion method is introduced for calibrating an aircraft equipped with four or more electric field sensors and a high-voltage corona point that is capable of charging the aircraft. An analytic, closed form solution for the estimate of a (3 x 3) aircraft calibration matrix is derived, and an absolute calibration experiment is used to improve the relative magnitudes of the elements of this matrix. To demonstrate the calibration procedure, we analyze actual calibration date derived from a Lear jet 28/29 that was equipped with five shutter-type field mill sensors (each with sensitivities of better than 1 V/m) located on the top, bottom, port, starboard, and aft positions. As a test of the calibration method, we analyze computer-simulated calibration data (derived from known aircraft and ambient fields) and explicitly determine the errors involved in deriving the variety of calibration matrices. We extend our formalism to arrive at an analytic solution for the ambient field, and again carry all errors explicitly.

Koshak, William J.↗

Diffusion model for lightning radiative transfer

A one-speed Boltzmann transport theory, with diffusion approximations, is applied to study the radiative transfer properties of lightning in optically thick thunderclouds. Near-infrared (lambda = 0.7774 micrometers) photons associated with a prominent oxygen emission triplet in the lightning spectrum are considered. Transient and spatially complex lightning radiation sources are placed inside a rectangular parallelepiped thundercloud geometry and the effects of multiple scattering are studied. The cloud is assumed to be composed of a homogeneous collection of identical spherical water droplets, each droplet a nearly conservative, anisotropic scatterer. Conceptually, we treat the thundercloud like a nuclear reactor, with photons replaced by neutrons, and utilize standard one-speed neutron diffusion techniques common in nuclear reactor analyses. Valid analytic results for the intensity distribution (expanded in spherical harmonics) are obtained for regions sufficiently far from sources. Model estimates of the arrival-time delay and pulse width broadening of lightning signals radiated from within the cloud are determined and the results are in good agreement with both experimental data and previous Monte Carlo estimates. Additional model studies of this kind will be used to study the general information content of cloud top lightning radiation signatures.

Koshak, William J.↗

Time-averaged current analysis of a thunderstorm using ground-based measurements

The amount of upward current provided to the ionosphere by a thunderstorm that appeared over the Kennedy Space Center (KSC) on July 11, 1978, is reexamined using an analytic equation that describes a bipolar thunderstorm's current contribution to the global circuit in terms of its generator current, lightning currents, the altitudes of its charge centers, and the conductivity profile of the atmosphere. Ground-based measurements, which were obtained from a network of electric field mills positioned at various distances from the thunderstorm, were used to characterize the electrical activity inside the thundercloud. The location of the lightning discharges, the type of lightning, and the amount of charge neutralized during this thunderstorm were computed through a least squares inversion of the measured changes in the electric fields following each lightning discharge. These measurements provided the information necessary to implement the analytic equation, and consequently, a time-averaged estimate of this thunderstorm's current contribution to the global circuit was calculated. From these results the amount of conduction current supplied to the ionosphere by this small thunderstorm was computed to be less than 25% of the time-averaged generator current that flowed between the two vertically displaced charge centers.

Driscoll, Kevin T.↗

A linear method for analyzing lightning field changes

A constrained, least-squares method for analyzing multiple-station measurements of lightning field changes (delta Es) is introduced. Previous methods have attempted to fit the spatial pattern of lightning delta Es using nonlinear models, such as a point charge (Q) or a point dipole (P) model. With the linear method, the delta Es are described not by models but by a general volume charge distribution that is deposited on a large (40 x 40 x 20 cu km) Cartesian grid above the measuring network. A linear system of equations is used to relate the measured delta Es to the charges that are deposited at each grid point. With this approach, the information content of the measurements can be quantified by an eigenanalysis of the covariance matrix of the linear system. Constraints can be used to reduce the infinity of possible solutions to the linear system and also to reduce systematic biases that can be introduced by the method of solution. It is shown that a Landweber iterative method, derived from the general method of steepest descent, can be used to solve the linear system and that the resulting volume charge distributions are generally consistent with computer-simulated charge sources, when these sources are over the measuring network. The Landweber iteration has also provided solutions for natural lightning events that are consistent with Q- and P-model results.

Koshak, William J.↗

Analysis of lightning field changes during active Florida thunderstorms

A computer algorithm for deriving accurate values of lightning-caused changes in cloud electric fields under active storm conditions was developed and applied to data obtained during two thunderstorms from a network of ground-based electric field mills at the NASA Kennedy Space Center and the Cape Canaveral Air Force Station. The resulting field changes were analyzed using a least-squares optimization procedure and point-charge (Q) and point-dipole (P) models. The values and the time variations of the Q-model parameters under active storm conditions were found to be similar to those reported by Maier and Krider (1986) for small storms, when the computations were carried out with the same analysis criteria and comparable biases. The parameters of the P solutions were found to vary with time within the storm interval and from storm to storm.

Koshak, William J.↗