Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “neutral model”

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 523 records · Page 29

Neutralizer Hollow Cathode Simulations and Validation with Experiments

The fidelity of electric propulsion physics-based models depends largely on the validity of their predictions over a range of operating conditions and geometries. In general, increased complexity of the physics requires more extensive comparisons with laboratory data to identify the region(s) that lie outside the validity of the model assumptions and to quantify the uncertainties within its range of application. This paper presents numerical simulations of neutralizer hollow cathodes at various operating conditions and orifice sizes. The simulations were performed using a two-dimensional axisymmetric model that solves numerically a relatively extensive system of conservation laws for the partially-ionized gas in these devices. The results for the plasma are compared directly with Langmuir probe measurements. The computed keeper voltages are also compared with the observed values. Wherever model inputs and/or specific physics of the cathode discharge are uncertain, additional sensitivity calculations have been performed to quantify the uncertainties.

partially ionized gases↗

DSMC simulation of ionized rarefied flows

Recently a renewed interest has been exhibited in understanding contamination and space environmental effect in the LEO range of 200 to 600 km as well as in low density microelectronics manufacturing technologies. Realistic simulations must model the physics of the highly coupled effects of neutral and charged particle flows, thermodynamic nonequilibrium, surface charging, and electromagnetic field effects. The computational requirements are enormous for this level of modeling and are almost impossible on current single processor supercomputers. An effort was initiated to develop this capability on a massively parallel MIMD supercomputers. The Wake Shield Facility experiment will be used for test case calculations. Preliminary results indicate that parallel supercomputers may be the only computers capable of simulating the required level of physics and spatial resolution. Calculations which contain from 1 to 8 million simulation particles are presented.

Bartel, Timothy J.↗

Radiation Transport of Heliospheric Lyman-alpha from Combined Cassini and Voyager Data Sets

Heliospheric neutral hydrogen scatters solar Lyman-alpha radiation from the Sun with '27-day' intensity modulations observed near Earth due to the Sun's rotation combined with Earth's orbital motion. These modulations are increasingly damped in amplitude at larger distances from the Sun due to multiple scattering in the heliosphere, providing a diagnostic of the interplanetary neutral hydrogen density independent of instrument calibration. This paper presents Cassini data from 2003-2004 obtained downwind near Saturn at approximately 10 AU that at times show undamped '27-day' waves in good agreement with the single-scattering models of Pryor et al., 1992. Simultaneous Voyager 1 data from 2003- 2004 obtained upwind at a distance of 88.8-92.6 AU from the Sun show waves damped by a factor of -0.21. The observed degree of damping is interpreted in terms of Monte Carlo multiple-scattering calculations (e.g., Keller et al., 1981) applied to two heliospheric hydrogen two-shock density distributions (discussed in Gangopadhyay et al., 2006) calculated in the frame of the Baranov-Malama model of the solar wind interaction with the two-component (neutral hydrogen and plasma) interstellar wind (Baranov and Malama 1993, Izmodenov et al., 2001, Baranov and Izmodenov, 2006). We conclude that multiple scattering is definitely occurring in the outer heliosphere. Both models compare favorably to the data, using heliospheric neutral H densities at the termination shock of 0.085 cm(exp -3) and 0.095 cm(exp -3). This work generally agrees with earlier discussions of Voyager data in Quemerais et al., 1996 showing the importance of multiple scattering but is based on Voyager data obtained at larger distances from the Sun (with larger damping) simultaneously with Cassini data obtained closer to the Sun.

Pryor, W.↗

Pickup Protons: Comparisons using the Three-Dimensional MHD HHMS-PI model and Ulysses SWICS Measurements

We report the first comparisons of pickup proton simulation results with in situ measurements of pickup protons obtained by the SWICS instrument on Ulysses. Simulations were run using the three dimensional (3D) time-dependent Hybrid Heliospheric Modeling System with Pickup Protons (HHMS-PI). HHMS-PI is an MHD solar wind model, expanded to include the basic physics of pickup protons from neutral hydrogen that drifts into the heliosphere from the local interstellar medium. We use the same model and input data developed by Detman et al. (2011) to now investigate the pickup protons. The simulated interval of 82 days in 2003 2004, includes both quiet solar wind (SW) and also the October November 2003 solar events (the Halloween 2003 solar storms). The HHMS-PI pickup proton simulations generally agree with the SWICS measurements and the HHMS-PI simulated solar wind generally agrees with SWOOPS (also on Ulysses) measurements. Many specific features in the observations are well represented by the model. We simulated twenty specific solar events associated with the Halloween 2003 storm. We give the specific values of the solar input parameters for the HHMS-PI simulations that provide the best combined agreement in the times of arrival of the solar-generated shocks at both ACE and Ulysses. We show graphical comparisons of simulated and observed parameters, and we give quantitative measures of the agreement of simulated with observed parameters. We suggest that some of the variations in the pickup proton density during the Halloween 2003 solar events may be attributed to depletion of the inflowing local interstellar medium (LISM) neutral hydrogen (H) caused by its increased conversion to pickup protons in the immediately preceding shock.

Intriligator, Devrie S.↗

Gas-Grain Models for Interstellar Anion Chemistry

Long-chain hydrocarbon anions C(sub n) H(-) (n = 4, 6, 8) have recently been found to be abundant in a variety of interstellar clouds. In order to explain their large abundances in the denser (prestellar/protostellar) environments, new chemical models are constructed that include gas-grain interactions. Models including accretion of gas-phase species onto dust grains and cosmic-ray-induced desorption of atoms are able to reproduce the observed anion-to-neutral ratios, as well as the absolute abundances of anionic and neutral carbon chains, with a reasonable degree of accuracy. Due to their destructive effects, the depletion of oxygen atoms onto dust results in substantially greater polyyne and anion abundances in high-density gas (with n(sub H2) approx > / cubic cm). The large abundances of carbon-chain-bearing species observed in the envelopes of protostars such as L1527 can thus be explained without the need for warm carbon-chain chemistry. The C6H(-) anion-to-neutral ratio is found to be most sensitive to the atomic O and H abundances and the electron density. Therefore, as a core evolves, falling atomic abundances and rising electron densities are found to result in increasing anion-to-neutral ratios. Inclusion of cosmic-ray desorption of atoms in high-density models delays freeze-out, which results in a more temporally stable anion-to-neutral ratio, in better agreement with observations. Our models include reactions between oxygen atoms and carbon-chain anions to produce carbon-chain-oxide species C6O, C7O, HC6O, and HC7O, the abundances of which depend on the assumed branching ratios for associative electron detachment

Cordiner, M. A.↗

Post sunset behavior of the 6300 A atomic oxygen airglow emission

A theoretical model of the 6300 A OI airglow emission was developed based on the assumptions that both the charged and neutral portions of the Earth's upper atmosphere are in steady state conditions of diffusive equilibrium. Intensities of 6300 A OI emission line were calculated using electron density true height profiles from a standard C-4 ionosonde and exospheric temperatures derived from Fabry-Perot interferometer measurements of the Doppler broadened 6300 A emission line shape as inputs to the model. Reaction rate coefficient values, production mechanism efficiencies, solar radiation fluxes, absorption cross sections, and models of the neutral atmosphere were varied parametrically to establish a set of acceptable inputs which will consistently predict 6300 A emission intensities that closely agree with intensities observed during the post-sunset twilight period by an airglow observatory consisting of a Fabry-Perot interferometer and a turret photometer. Emission intensities that can only result from the dissociative recombination of molecular oxygen ions were observed during the latter portion of the observational period. Theoretical calculations indicate that contamination of the 6300 A OI emission should be on the order of or less than 3 percent; however, these results are very sensitive to the wavelengths of the individual lines and their intensities relative to the 6300 A OI intensity. This combination of a model atmosphere, production mechanism efficiencies, and quenching coefficient values was used when the dissociative photoexcitation and direct impact excitation processes were contributing to the intensity to establish best estimates of solar radiation fluxes in the Schumann--Runge continuum and associated absorption cross sections. Results show that the Jacchia 1971 model of the upper atmosphere combined with the Ackerman recommended solar radiation fluxes and associated absorption cross sections produces theoretically calculated intensities that more closely agree with the observed intensities than all the other combinations.

Smith, R. E.↗

Advancing Solid Earth Science through Improved Atmosphere Modeling

We proposed to investigate and develop better models for the effect of the hydrostatic and water vapor components of the neutral atmosphere on delay for VLBI and GPS by using a Numerical Weather Model to better simulate realistic atmosphere conditions. By using a raytrace calculation through the model atmosphere at the times of actual VLBI observations, the potential improvement in geodetic results can be evaluated. Also, by calculating the actual variation of delays with elevation and azimuth, the errors in current mapping function models can be assessed. The VLBI data to be initially analyzed are the fifteen days of the CONT02 sessions of 2002 October which included eight stations. There are three segments to the research. 1) The PSU/NCAR fifth generation mesoscale numerical weather model (MM5) will be used to provide the state of the atmosphere with highest horizontal resolution of 3 km. 2) A three-dimensional raytrace program will be developed to determine the delays through the model atmosphere at the times and in the directions of the VLBI observations for each of the sites. 3) The VLBI data will be analyzed using both standard models for the atmosphere mapping functions and the mapping functions derived from the NWM raytracing.

Niell, A. E.↗

The role of O-type neutral lines in magnetic merging during substorms and solar flares

An analytical model of O-type configurations was derived, for visualizing the geometry and for numerical treatment of plasma flows. It was found that at a certain distance from the neutral line the mean particle motion became decoupled from that of magnetic field lines (which obey the MHD condition). The decoupling distance depended on initial conditions in momentum space, suggesting that the MHD approximation which averages out such conditions may not suffice for describing plasma dynamics near the neutral line. It was also found that after inflowing particles are decoupled from the field line motion, they go over to a mode of runaway acceleration along the neutral line. It is concluded that if merging occurs at an X-O pair, two particle populations may be expected - low energy particles accelerated adiabatically by earthward convection past the X-type line and high-energy particles convected towards the O-type line. The second acceleration process depends critically on the rapidity of merging and is therefore expected to vary considerably from event to event.

Stern, D. P.↗

MPD Thruster Performance Analytic Models

Magnetoplasmadynamic (MPD) thrusters are capable of accelerating quasi-neutral plasmas to high exhaust velocities using Megawatts (MW) of electric power. These characteristics make such devices worthy of consideration for demanding, far-term missions such as the human exploration of Mars or beyond. Assessment of MPD thrusters at the system and mission level is often difficult due to their status as ongoing experimental research topics rather than developed thrusters. However, in order to assess MPD thrusters utility in later missions, some adequate characterization of performance, or more exactly, projected performance, and system level definition are required for use in analyses. The most recent physical models of self-field MPD thrusters have been examined, assessed, and reconfigured for use by systems and mission analysts. The physical models allow for rational projections of thruster performance based on physical parameters that can be measured in the laboratory. The models and their implications for the design of future MPD thrusters are presented.

Gilland, James↗

MPD Thruster Performance Analytic Models

Magnetoplasmadynamic (MPD) thrusters are capable of accelerating quasi-neutral plasmas to high exhaust velocities using Megawatts (MW) of electric power. These characteristics make such devices worthy of consideration for demanding, far-term missions such as the human exploration of Mars or beyond. Assessment of MPD thrusters at the system and mission level is often difficult due to their status as ongoing experimental research topics rather than developed thrusters. However, in order to assess MPD thrusters utility in later missions, some adequate characterization of performance, or more exactly, projected performance, and system level definition are required for use in analyses. The most recent physical models of self-field MPD thrusters have been examined, assessed, and reconfigured for use by systems and mission analysts. The physical models allow for rational projections of thruster performance based on physical parameters that can be measured in the laboratory. The models and their implications for the design of future MPD thrusters are presented.

Gilland, James↗

HMI Data Driven Magnetohydrodynamic Model Predicted Active Region Photospheric Heating Rates: Their Scale Invariant, Flare Like Power Law Distributions, and Their Possible Association With Flares

A data driven, near photospheric, 3 D, non-force free magnetohydrodynamic model predicts time series of the complete current density, and the resistive heating rate Q at the photosphere in neutral line regions (NLRs) of 14 active regions (ARs). The model is driven by time series of the magnetic field B observed by the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory (SDO) satellite. Spurious Doppler periods due to SDO orbital motion are filtered out of the time series for B in every AR pixel. Errors in B due to these periods can be significant. The number of occurrences N(q) of values of Q > or = q for each AR time series is found to be a scale invariant power law distribution, N(Q) / Q−s, above an AR dependent threshold value of Q, where 0.3952 < or = s < or = 0.5298 with mean and standard deviation of 0.4678 and 0.0454, indicating little variation between ARs. Observations show that the number of occurrences N(E) of coronal flares with a total energy released > or = E obeys the same type of distribution, N(E) / E−S, above an AR dependent threshold value of E, with 0.38 < or approx. S < or approx. 0.60, also with little variation among ARs. Within error margins the ranges of s and S are nearly identical. This strong similarity between N(Q) and N(E) suggests a fundamental connection between the process that drives coronal flares and the process that drives photospheric NLR heating rates in ARs. In addition, results suggest it is plausible that spikes in Q, several orders of magnitude above background values, are correlated with times of the subsequent occurrence of M or X flares.

Goodman, Michael L.↗

Magnetosphere of Saturn

Models of the Saturnian magnetosphere based on the application of magnetospheric scaling relations to a spin-aligned planetary magnetic dipole, that produces a surface equatorial field strength in the range 0.5 to 2 gauss, exhibit the following properties: (1) The orbit of Titan lies inside of the magnetosphere essentially all of the time, even when variations in the size of the magnetosphere resulting from solar wind pressure changes are taken into account; (2) the Brice-type planetary plasmasphere reaches a peak density of about 10 protons cm/3 at L approximately 7 (L = planetocentric distance in units of planetary radii); (3) Saturn's rings have a profound effect on the energetic particle population and the plasmaspheres derived from interstellar neutrals and Titan's torus; (4) the model calculation suggests that the Titan-derived plasmasphere may be self-amplifying with a feed-back factor greater than unity, which implies the possibility of a non-linearly saturated, highly inflated Saturnian magnetosphere; and (5) this same source can have important eroding effects on the outer edge of the rings determined by Brown-Lauzerotti sputtering rates.

Siscoe, G. L.↗

Helium escape from the Earth's atmosphere - The charge exchange mechanism revisited

We have studied the escape of neutral helium from the terrestrial atmosphere through exothermic charge exchange reactions between He(+) ions and the major atmospheric constituents N2, O2 and O. Elastic collisions with the neutral background particles were treated quantitatively using a recently developed kinetic theory approach. An interhemispheric plasma transport model was employed to provide a global distribution of He(+) ions as a function of altitude, latitude and local solar time and for different levels of solar ionization. Combining these ion densities with neutral densities from an MSIS model and best estimates for the reaction rate coefficients of the charge exchange reactions, we computed the global distribution of the neutral He escape flux. The escape rates show large diurnal and latitudinal variations, while the global average does not vary by more than a factor of three over a solar cycle. We find that this escape mechanism is potentially important for the overall balance of helium in the Earth's atmosphere. However, more accurate values for the reaction rate coefficients of the charge exchange reactions are required to make a definitive assessment of its importance.

Lie-Svendsen, O.↗

Etching in Chlorine Discharges Using an Integrated Feature Evolution-Plasma Model

Etching of semiconductor materials is reliant on plasma properties. Quantities such as ion and neutral fluxes, both in magnitude and in direction, are often determined by reactor geometry (height, radius, position of the coils, etc.) In order to obtain accurate etching profiles, one must also model the plasma as a whole to obtain local fluxes and distributions. We have developed a set of three models that simulates C12 plasmas for etching of silicon, ion and neutral trajectories in the plasma, and feature profile evolution. We have found that the location of the peak in the ion densities in the reactor plays a major role in determining etching uniformity across the wafer. For a stove top coil inductively coupled plasma (ICP), the ion density is peaked at the top of the reactor. This leads to nearly uniform neutral and ion fluxes across the wafer. A side coil configuration causes the ion density to peak near the sidewalls. Ion fluxes are thus greater toward the wall's and decrease toward the center. In addition, the ions bombard the wafer at a slight angle. This angle is sufficient to cause slanted profiles, which is highly undesirable.

Hwang, Helen H.↗

The Effects of Neutral Inertia on Ionospheric Currents in the High-Latitude Thermosphere Following a Geomagnetic Storm

Results of an experimental and theoretical investigation into the effects of the time dependent neutral wind flywheel on high-latitude ionospheric electrodynamics are presented. The results extend our previous work which used the National Center for Atmospheric Research Thermosphere/Ionosphere General Circulation Model (NCAR TIGCM) to theoretically simulate flywheel effects in the aftermath of a geomagnetic storm. The previous results indicated that the neutral circulation, set up by ion-neutral momentum coupling in the main phase of a geomagnetic storm, is maintained for several hours after the main phase has ended and may dominate height-integrated Hall currents and field-aligned currents for up to 4-5 hours. We extend the work of Deng et al. to include comparisons between the calculated time-dependent ionospheric Hall current system in the storm-time recovery period and that measured by instruments on board the Dynamics Explorer 2 (DE 2) satellite. Also, comparisons are made between calculated field-aligned currents and those derived from DE 2 magnetometer measurements. These calculations also allow us to calculate the power transfer rate (sometimes called the Poynting flux) between the magnetosphere and ionosphere. The following conclusions have been drawn: (1) Neutral winds can contribute significantly to the horizontal ionospheric current system in the period immediately following the main phase of a geomagnetic storm, especially over the magnetic polar cap and in regions of ion drift shear. (2) Neutral winds drive Hall currents that flow in the opposite direction to those driven by ion drifts. (3) The overall morphology of the calculated field-aligned current system agrees with previously published observations for the interplanetary magnetic field (IMF) B(sub Z) southward conditions, although the region I and region 2 currents are smeared by the TI(ICM model grid resolution. (4) Neutral winds can make significant contributions to the field-aligned current system when B(sub Z) northward conditions prevail following the main phase of a storm, but can account for only a fraction of the observed currents. (5) DE 2 measurements provide a demonstration of "local" (satellite-altitude) flywheel effects. (6) On the assumption that the magnetosphere acts as an insulator, we calculate neutral-wind-induced polarization electric fields of approx. 20-30 kV in the period immediately following the geomagnetic storm.

Deng, W.↗

Thermospheric and ionospheric structure of the Southern Hemisphere polar cap on October 21, 1981, as determined from Dynamics Explorer 2 satellite data

For a number of years, satellites have been employed to measure auroral particles and fields within the high-latitude thermosphere. In the present paper, data from orbit 1174 of the Dynamics Explorer 2 satellite on October 21, 1981, are utilized to calculate the structure of the ionosphere and thermosphere below the satellite altitude down to about 80 km. Attention is given to details regarding the DE 2 measurements, a satellite track model, the calculated ionospheric structure, and the calculated neutral gas heating rates. The investigation demonstrates that the technique of deriving characteristics of the ionosphere and atmosphere below a satellite track promises to be very fruitful for defining the characteristics of the lower thermosphere for use in large numerical models of the thermosphere and ionosphere.

Emery, B. A.↗

Planetary Boundary Layer Simulation Using TASS

Boundary conditions to an existing large-eddy simulation model have been changed in order to simulate turbulence in the atmospheric boundary layer. Several options are now available, including the use of a surface energy balance. In addition, we compare convective boundary layer simulations with the Wangara and Minnesota field experiments as well as with other model results. We find excellent agreement of modelled mean profiles of wind and temperature with observations and good agreement for velocity variances. Neutral boundary simulation results are compared with theory and with previously used models. Agreement with theory is reasonable, while agreement with previous models is excellent.

Schowalter, David G.↗