Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Abel Inversion”

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 37 records · Page 2

Flame Chemiluminescence Rate Constants for Quantitative Microgravity Combustion Diagnostics

Absolute excited state concentrations of OH(A), CH(A), and C2(d) were determined in three low pressure premixed methane-air flames. Two dimensional images of chemiluminescence from these states were recorded by a filtered CCD camera, processed by Abel inversion, and calibrated against Rayleigh scattering, Using a previously validated 1-D flame model with known chemistry and excited state quenching rate constants, rate constants are extracted for the reactions CH + O2 (goes to) OH(A) + CO and C2H + O (goes to) CH(A) + CO at flame temperatures. Variations of flame emission intensities with stoichiometry agree well with model predictions.

Luque, Jorge↗

Concentration Measurements in Self-Excited, Momentum-Dominated Helium Jets

Flow structure of momentum-dominated pure helium jets discharged vertically into ambient air was investigated using high-speed rainbow schlieren deflectometry (RSD) technique. Effects of the operating parameters, i.e., Reynolds number (Re) and Richardson number (Ri), on the oscillatory behavior of the flow were examined over a range of experimental conditions. To seek the individual effect of these parameters, one of them was fixed and the other was varied with certain constraints. Measurements revealed highly periodic oscillations in the laminar region as well as high regularity in transition and turbulent regions. Maximum spectral power profiles at different axial locations indicated the oscillation amplitude increasing until the breakdown of the jet in the turbulent regime. The transition from the laminar to turbulent flow was also investigated. Fast Fourier transform analysis performed in the transition regime showed that the flow oscillates at a unique frequency, which was the same in the upstream laminar flow region. Measured deflection angle data were used in Abel inversion algorithm to construct the helium concentration fields. Instantaneous helium concentration contours revealed changes in the flow structure and evolution of vortical structures during an oscillation cycle. Temporal evolution plots of helium concentration at different axial location showed repeatable oscillations at all axial and radial locations up to the turbulent regime. A cross-correlation technique, applied to find the spatial displacements of the vortical structures, provided correlation coefficient peaks between consecutive schlieren images. Results show that the vortical structure convected and accelerated only in the axial direction.

Yildirim, Bekir Sedat↗

Concentration Measurements in Self-Excited Momentum Dominated Low-Density Gas Jets

Flow structure of self-excited, laminar, axisymmetric, momentum-dominated helium jets discharged vertically into ambient air was investigated using high-speed rainbow schlieren deflectometry technique. Measurements were obtained at temporal resolution of 1 ms and spatial resolution of 0.19 mm for two test cases with Richardson number of 0.034 and 0.018. Power spectra revealed that the oscillation frequency was independent of spatial coordinates, suggesting global oscillations in the flow. Abel inversion algorithm was used to reconstruct the concentration field of helium. Instantaneous concentration contours revealed changes in the flow field and evolution of vortical structures during an oscillation cycle. Temporal evolution plots of helium concentration at different axial locations provided detailed information about the instability in the flow field.

Yildirim, B. S.↗

Full Field Scalar Measurements in a Pulsating Helium Jet using Rainbow Schlieren Deflectometry

The flow structure of a pulsating helium jet was investigated using quantitative rainbow schlieren deflectometry. The operating parameters included the tube inside diameter, the jet Reynolds number and the jet Richardson number. The jet structure was characterized by the frequency spectrum, temporal evolution of the oscillations and mean and root-mean-square profiles of the species mole fraction, which in this case, was helium. Experiments were conducted using a variable nozzle facility. Angular deflection data were obtained using rainbow schlieren deflectometry across full field of color images taken at a temporal resolution of 60Hz. The flicker cycle was analyzed by instantaneous rainbow schlieren images, contour plots of deflection angle and helium mole fraction placed sequentially in time. Two observe the flow downstream, the tube was lowered with respect to the optical setup and images were taken at different axial planes. In this way, the laminar, transition and turbulent regions of the jet flow were observed. Abel inversion algorithm was utilized to reconstruct the refractive index field from the measurements of the beam deflection angle. The concentration field was then generated from the refractive index field. The phenomenon of vortex initiation and propagation in the flicker cycle was described by correlating the ray deflection angle and concentration contour plots. Experiments in the flickering jet reveal global oscillation in the flow field. The effect of jet exit Reynolds number and jet Richardson number on the flickering frequency was analyzed. The effect of jet Richardson number was more pronounced than that of the jet exit Reynolds number. The flow field was studies quantitatively in terms of temporal evolution and statistical description of helium mole fraction.

Pasumarthi, Kasyap S.↗

Local Neutral Density and Plasma Parameter Measurements in a Hollow Cathode Plume

In order to understand the cathode and keeper wear observed during the Extended Life Test (ELT) of the DS1 flight spare NSTAR thruster and provide benchmarking data for a 2D cathode/cathode-plume model, a basic understanding of the plasma and neutral gas parameters in the cathode orifice and keeper region of the cathode plume must be obtained. The JPL cathode facility is instrumented with an array of Langmuir probe diagnostics along with an optical diagnostic to measure line intensity of xenon neutrals. In order to make direct comparisons with the present model, a flat plate anode arrangement was installed for these tests. Neutral density is deduced from the scanning probe data of the plasma parameters and the measured xenon line intensity in the optical regime. The Langmuir probes are scanned both axially, out to 7.0 cm downstream of the keeper, and radially to obtain 2D profile of the plasma parameters. The optical fiber is housed in a collimating stainless steel tube, and is scanned to view across the cathode plume along cuts in front of the keeper with a resolution of 1.5 mm. The radial intensities are unfolded using the Abel inversion technique that produces radial profiles of local neutral density. In this paper, detailed measurements of the plasma parameters and the local neutral densities will be presented in the cathode/keeper plume region for a 1.5 cm diameter NEXIS cathode at 25A of discharge current at several different strengths of applied magnetic field.

Plasma density↗

Storm/Quiet Ratio Comparisons Between TIMED/SABER NO (sup +)(v) Volume Emission Rates and Incoherent Scatter Radar Electron Densities at E-Region Altitudes

Broadband infrared limb emission at 4.3 microns is measured by the TIMED/SABER instrument. At night, these emission observations at E-region altitudes are used to derive the so called NO+(v) Volume Emission Rate (VER). NO+(v) VER can be derived by removing the background CO2(v3) 4.3 microns radiance contribution using SABER-based non-LTE radiation transfer models, and by performing a standard Abel inversion on the residual radiance. SABER observations show that NO+(v) VER is significantly enhanced during magnetic storms in accordance with increased ionization of the neutral atmosphere by auroral electron precipitation, followed by vibrational excitation of NO+ (i.e., NO+(v)) from fast exothermic ion-neutral reactions, and prompt infrared emission at 4.3 m. Due to charge neutrality, the NO+(v) VER enhancements are highly correlated with electron density enhancements, as observed for example by Incoherent Scatter Radar (ISR). In order to characterize the response of the storm-time E-region from both SABER and ISR measurements, a Storm/Quiet ratio (SQR) quantity is defined as a function of altitude. For SABER, the SQR is the ratio of the storm-to-quiet NO+(v) VER. SQR is the storm-to-quiet ratio of electron densities for ISR. In this work, we compare SABER and ISR SQR values between 100 to 120 km. Results indicate good agreement between these measurements. SQR values are intended to be used as a correction factor to be included in an empirical storm-time correction to the International Reference Ionosphere model at E-region altitudes.

Fernandez, J. R.↗

Ionospheric Mapping with the GPS/MET

The GPS/MET experiment, which placed a GPS receiver in a low-Earth orbit tracking the GPS in an occultation geometry, has collected thousands of occultations since its launch in April of 1995. Each occultation can be inverted to give an electron density profile in the ionosphere, temperature and pressure profiles in the lower mesosphere, stratosphere and upper troposphere, and water vapor density profiles in the lower troposphere. This paper gives a summary of the ionospheric effects on the GPS/MET signal and examines some of the retrieved electron density profiles. We examine the bending induced by the ionosphere on the occulting signal and the resulting separation of the two GPS links corresponding to the L1 and L2 phase signals. We also examine the amplitude scintillation caused by sharp layers at the bottom of the ionosphere. We briefly describe the Abel inversion method for obtaining the index of refraction and show several examples of electron density profiles and comparisons derived from the Parameterized Ionospheric Model (PIM) and to incoherent scatter radar measurements.

GPS↗

Collisional-Radiative Models of Neutral and Singly-Ionized Xenon in Hall Thrusters: Experimental Validation and Model Investigations

The recently developed collisional radiative (CR) models for neutral xenon (Xe I) and singly charged xenon ion (Xe II) described by Chaplin et al. in the paper "Insights from Collisional-Radiative Models of Neutral and Singly-Ionized Xenon in Hall Thrusters" also submitted to the 2021 AIAA Propulsion and Energy Forum are more thoroughly investigated. Comparisons between the existing and recently developed JPL Xe I CR models are done and reasons for discrepancies between the two models are explored. An experiment to validate the Xe II CR model is carried out in the anode region of the UCLA Plasma interactions facility. The Xe II CR model is applied to measured xenon ion emission line intensities to determine electron temperature and density. Spatial resolution of these optical measurements is obtained through an Abel inversion for direct comparison to Langmuir probe results. Although the measured electron densities were on the order of $10^{17}$ m$^{-3}$ and the Xe II CR model is most suitable for use with densities above $10^{18}$ m$^{-3}$ due to the strong dependence of the line ratios on electron temperature in that regime, initial results from the validation experiment show reasonable agreements between the optical and Langmuir probe results. The 441 nm to 605 nm line ratio appears to have the most promise for non-intrusive electron temperature measurements with the Xe II CR model. Shortcomings of the Xe II CR model revealed in the validation experiment are to be explored further.

Lobbia, Robert B.↗

A Theoretical and Computational Revisit of Conversions Between Whitham’s F-Function and Equivalent Area

This paper proves mathematically that the integral transforms between Whitham’s F-function and equivalent area are the inverse transforms of each other if and only if the slope of the equivalent area at the origin is zero. This mathematical fact contradicts the accepted unconditional inverse relation between Whitham’s F-function and equivalent area in the sonic boom research literature. Piecewise linear approximations of an F-function and of the second derivative of an equivalent area are used to derive numerical formulas for conversions between Whitham’s F-function and equivalent area. Numerical results are included to show convergence of the numerical conversions as the maximum length of the segments for piecewise linear approximations goes to zero. These numerical conversions are approximately the inverse transforms of each other when the second derivative of an equivalent area is continuous and the slope of the equivalent area at the origin is zero.

Whitham's F-function↗

Current methods of radio occultation data inversion

The methods of Abel integral transform and ray-tracing inversion have been applied to data received from radio occultation experiments as a means of obtaining refractive index profiles of the ionospheres and atmospheres of Mars and Venus. In the case of Mars, certain simplifications are introduced by the assumption of small refractive bending in the atmosphere. General inversion methods, independent of the thin atmosphere approximation, have been used to invert the data obtained from the radio occultation of Mariner 5 by Venus; similar methods will be used to analyze data obtained from Jupiter with Pioneers F and G, as well as from the other outer planets in the Outer Planet Grand Tour Missions.

Kliore, A. J.↗

Some procedures for computerized electronic data processing of absorption measurements from artificial earth satellites

The processing of data obtained from solar absorption radiation measurements is discussed. The position of the satellite was obtained by numerical integration of the differential equations of motion using initial conditions. The position of the sun was calculated as a function of time, and the tangential elevation was determined approximately from the positions of the satellite and the sun. The coefficients of an approximation formula and of a data smoothing process were determined, and the inversion of an Abel integral equation is solved analytically.

Trinkkeller, B.↗

Structure of the atmosphere of Jupiter from Pioneer 10 and 11 radio occultation measurements

A technique for incorporating effects of the oblateness of Jupiter's atmosphere into the data analysis procedure for radio-occultation measurements is described which makes use of a spherical harmonic representation of the gravity field to compute the shape of the planet. With this technique, the center of refraction is located by the radius of curvature and the normal direction at the closest approach point of the ray. The present technique, a subsequent approximation involving the use of a constant center of refraction, and the Abel integral transform inversion method for spherical planets are applied to Pioneer 10 and 11 data. All the intermediate results obtained show a temperature inversion between 10 and 100 mb, with temperatures from 130 to 170 K at 10 mb and from 80 to 120 K at 100 mb. Comparison of the radio-occultation profiles with radiative-convective equilibrium models and the temperature profile based on Pioneer 10 IR radiometer data indicates good agreement between the models and the occultation results.

Kliore, A. J.↗

Electron beam dispersion measurements in nitrogen using two-dimensional imaging of N2(+) fluorescence

Experimental results are presented related to the radial spread of fluorescence excited by 10 and 20 KeV electron beams passing through nonflowing rarefied nitrogen at 293 K. An imaging technique for obtaining species distributions from measured beam-excited fluorescence is described, based on a signal inversion scheme mathematically equivalent to the inversion of the Abel integral equation. From fluorescence image data, measurements of beam radius, integrated signal intensity, and spatially resolved distributions of N2(+) first-negative-band fluorescence-emitting species have been made. Data are compared with earlier measurements and with an heuristic beam spread model.

Clapp, L. H.↗

Ion beam probing of electrostatic fields

The determination of a cylindrically symmetric, time-independent electrostatic potential V in a magnetic field B with the same symmetry by measurements of the deflection of a primary beam of ions is analyzed and substantiated by examples. Special attention is given to the requirements on canonical angular momentum and total energy set by an arbitrary, nonmonotone V, to scaling laws obtained by normalization, and to the analogy with ionospheric sounding. The inversion procedure with the Abel analysis of an equivalent problem with a one-dimensional fictitious potential is used in a numerical experiment with application to the NASA Lewis Modified Penning Discharge. The determination of V from a study of secondary beams of ions with increased charge produced by hot plasma electrons is also analyzed, both from a general point of view and with application to the NASA Lewis SUMMA experiment. Simple formulas and geometrical constructions are given for the minimum energy necessary to reach the axis, the whole plasma, and any point in the magnetic field. The common, simplifying assumption that V is a small perturbation is critically and constructively analyzed; an iteration scheme for successively correcting the orbits and points of ionization for the electrostatic potential is suggested.

Persson, H.↗

Optimal Estimation Inversion of Ionospheric Electron Density from GNSS-POD Limb Measurements: Part II-Validation and Comparison Using NmF2 and hmF2

A growing number of SmallSat/CubeSat constellations with high-rate (50–100 Hz) global navigation satellite system radio occultations (GNSS-RO) as well as low-rate (1 Hz) precise orbit determination (GNSS-POD) limb-viewing capabilities provide unprecedented spatial and temporal sampling rates for ionospheric studies. In the F-region electron density (N 𝑒 ) retrieval process, instead of the conventional onion-peeling (OP) inversion, an optimal estimation (OE) inversion technique was recently developed using total electron content measurements acquired by GNSS-POD link. The new technique is applied to data acquired from the COSMIC-1, COSMIC-2, and Spire constellations. Although both OE and OP techniques use the Abel weighting function in N 𝑒 inversion, OE significantly differs in its performance, especially in the lower F- and E-regions. In this work, we evaluate and compare newly derived data sets using F2 peak properties with other space-based and ground-based observations. We determine the F2 peak N 𝑒 (NmF2) and its altitude (hmF2), and compare them with the OP-retrieved values. Good agreement is observed between the two techniques for both NmF2 and hmF2. In addition, we also utilize autoscaled F2 peak measurements from a number of worldwide Digisonde stations (∼30). The diurnal sensitivity and latitudinal variability of the F2 peak between the two techniques are carefully studied at these locations. Good agreement is observed between OE-retrieved NmF2 and Digisonde-measured NmF2. However, significant differences appear between OE-retrieved hmF2 and Digisonde-measured hmF2. During the daytime, Digisonde-measured hmF2 remains ∼25–45 km below the OE-retrieved hmF2, especially at mid and high latitudes. We also incorporate F-region N 𝑒 measurements from two incoherent scatter radar observations at high latitudes, located in the North American (Millstone Hill) and European (EISCAT at Tromso) sectors. The radar measurements show good agreement with OE-retrieved values. Although there are several possible sources of error in the ionogram-derived N 𝑒 profiles, our further analysis on F1 and F2 layers indicates that the low Digisonde hmF2 is caused by the autoscaled method, which tends to detect a height systematically below the F2 peak when the F1 layer is present.

GNSS radio occultation↗

Hierarchically Parallelized Constrained Nonlinear Solvers with Automated Substructuring

This paper develops a parallelizable multilevel multiple constrained nonlinear equation solver. The substructuring process is automated to yield appropriately balanced partitioning of each succeeding level. Due to the generality of the procedure,_sequential, as well as partially and fully parallel environments can be handled. This includes both single and multiprocessor assignment per individual partition. Several benchmark examples are presented. These illustrate the robustness of the procedure as well as its capability to yield significant reductions in memory utilization and calculational effort due both to updating and inversion.

Padovan, Joe↗