SEARCH · Engineering Papers
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.
A Poisson equation on a rectangular domain is solved by coupling two methods: the domain is divided in two squares, a finite element approximation is used on the first square and a spectral discretization is used on the second one. Two kinds of matching conditions on the interface are presented and compared. In both cases, error estimates are proved.
A Poisson equation on a rectangular domain is solved by coupling two methods: the domain is divided in two squares, a finite element approximation is used on the first square and a spectral discretization is used on the second one. Two kinds of matching conditions on the interface are presented and compared. In both cases, error estimates are proved.
A new method is proposed for selecting steady state shapes of disk galaxies as 'most probable states' of a large number of stars, given only the total energy and total angular momentum. A partial differential equation is derived for the mean gravitational potential; it is closely related to the 'sinh-Poisson' equation for the mean-field description of a line vortex system or electrostatic guiding-center plasma. A 'water bag' approximation to the distribution function for bound stars renders the equation analytically tractable, but accurate solution of it may require numerical integration, in view of its general nonlinearity.
Reflection from a one-dimensional random medium of discrete scatterers is considered. The discrete scattering medium is modeled by a Poisson impulse process with concentration lambda. By employing the Markov property of the Poisson impulse process, an exact functional integro-differential equation of the Kolmogorov-Feller type is found for the average reflected power. Approximate solutions to this equation are obtained by regular perturbation and two variable expansion techniques in the limit of small lambda. The regular perturbation results is valid for small slab thicknesses, while the two-variable result is uniformly valid for any thickness. The two-variable result shows that as the slab size becomes infinite all of the incident power is reflected on the average.
This paper handles a problem encountered in the solution by way of Chebyshev spectral collocation discretization of the Poisson equation with homogeneous Neumann boundary conditions. When pseudospectral Chebyshev methodology is utilized to discretize the partial differential equation these integrals have to be approximated by an appropriate quadrature formula. The correct quadrature formula is developed for these problems, taking into account the degree of the polynomials involved. It is shown that this formula leads to a well-conditioned Chebyshev approximation to the differential equations and that the compatibility condition is automatically satisfied.
Since 1750, the number of cataclysmic volcanic eruptions (volcanic explosivity index (VEI)>=4) per decade spans 2-11, with 96 percent located in the tropics and extra-tropical Northern Hemisphere. A two-point moving average of the volcanic time series has higher values since the 1860's than before, being 8.00 in the 1910's (the highest value) and 6.50 in the 1980's, the highest since the 1910's peak. Because of the usual behavior of the first difference of the two-point moving averages, one infers that its value for the 1990's will measure approximately 6.50 +/- 1, implying that approximately 7 +/- 4 cataclysmic volcanic eruptions should be expected during the present decade (2000-2009). Because cataclysmic volcanic eruptions (especially those having VEI>=5) nearly always have been associated with short-term episodes of global cooling, the occurrence of even one might confuse our ability to assess the effects of global warming. Poisson probability distributions reveal that the probability of one or more events with a VEI>=4 within the next ten years is >99 percent. It is approximately 49 percent for an event with a VEI>=5, and 18 percent for an event with a VEI>=6. Hence, the likelihood that a climatically significant volcanic eruption will occur within the next ten years appears reasonably high.
Many optical imaging, ranging, and communications systems rely on the estimation of the arrival time of an optical pulse. Recently, such systems have been increasingly employing photon-counting photodetector technology, which changes the statistics of the observed photocurrent. This requires time-of-arrival estimators to be developed and their performances characterized. The statistics of the output of an ideal photodetector, which are well modeled as a Poisson point process, were considered. An analytical model was developed for the mean-square error of the maximum likelihood (ML) estimator, demonstrating two phenomena that cause deviations from the minimum achievable error at low signal power. An approximation was derived to the threshold at which the ML estimator essentially fails to provide better than a random guess of the pulse arrival time. Comparing the analytic model performance predictions to those obtained via simulations, it was verified that the model accurately predicts the ML performance over all regimes considered. There is little prior art that attempts to understand the fundamental limitations to time-of-arrival estimation from Poisson statistics. This work establishes both a simple mathematical description of the error behavior, and the associated physical processes that yield this behavior. Previous work on mean-square error characterization for ML estimators has predominantly focused on additive Gaussian noise. This work demonstrates that the discrete nature of the Poisson noise process leads to a distinctly different error behavior.
Electron density measurements have been made in steady-state plasmas in a spherical inertial electrostatic confinement (IEC) discharge using microwave interferometry. Plasma cores interior to two cathodes, having diameters of 15 and 23 cm, respectively, were probed over a transverse range of 10 cm with a spatial resolution of about 1.4 cm for buffer gas pressures from 0.2 to 6 Pa in argon and deuterium. The transverse profiles are generally flat, in some cases with eccentric symmetric minima, and give mean densities of from approx. = 0.4 to 7x 10(exp 10)/cu cm, the density generally increasing with the neutral gas pressure. Numerical solutions of the 1-D Poisson equation for EC plasmas are reviewed and energy distribution functions are identified which give flat transverse profiles. These functions are used with the plasma approximation to obtain solutions which also give densities consistent with the measurements, and a double potential well solution is obtained which has minima qualitatively similar to those observed. Explicit consideration is given to the compatibility of the solutions interior and exterior to the cathode, and to grid transparency. Deuterium fusion neutron emission rates were also measured and found to be isotropic, to within the measurement error, over two simultaneous directions. Anisotropy was observed in residual emissions during operation with non-fusing hydrogen-1. The deuterium rates are consistent with predictions from the model.
Electron density measurements have been made in steady-state plasmas in a spherical inertial electrostatic confinement (IEC) discharge using microwave interferometry. Plasma cores interior to two cathodes, having diameters of 15 and 23 cm, respectively, were probed over a transverse range of 10 cm with a spatial resolution of about 1.4 cm for buffer gas pressures from 0.2 to 6 Pa in argon and deuterium. The transverse profiles are generally flat, in some cases with eccentric symmetric minima, and give mean densities of from approx. = 0.4 to 7 x 10(exp 10)/cu cm, the density generally increasing with the neutral gas pressure. Numerical solutions of the one-dimensional Poisson equation for IEC plasmas are reviewed and energy distribution functions are identified which give flat transverse profiles. These functions are used with the plasma approximation to obtain solutions which also give densities consistent with the measurements, and a double potential well solution is obtained which has minima qualitatively similar to those observed. Explicit consideration is given to the compatibility of the solutions interior and exterior to the cathode, and to grid transparency. Deuterium fusion neutron emission rates were also measured and found to be isotropic, to within the measurement error, over two simultaneous directions. Anisotropy was observed in residual emissions during operation with nonfusing hydrogen- 1. The deuterium rates are consistent with predictions from the model.
A pressurized cylindrical shell containing a longitudinal crack is considered. The shear modulus of the sheet may be evaluated from the measured Young's moduli and the Poisson's ratios rather than being an independent material constant. Two examples, one for a mildly orthotropic (titanium) and the other for a strongly orthotropic (graphite) material approximately satisfying the condition of special orthotropy are given. The results show that the stress intensity factors are rather strongly dependent on the degree of orthotropy.
A pressurized cylindrical shell containing a longitudinal crack is considered. It is assumed that the material has a special orthotropy - namely, that the shear modulus of the sheet may be evaluated from the measured Young's moduli and the Poisson's ratios and is not an independent material constant. Two examples, one for a mildly orthotropic (titanium) and the other for a strongly orthotropic (graphite) material approximately satisfying the condition of special orthotropy are given. The results show that the stress intensity factors are rather strongly dependent on the degree of orthotropy.
This paper reports on a theoretical study of the limitations on silicon solar-cell efficiency for both n(+)-p and n(+)-p-p(+) type cells. Detailed calculations have been made of solar-cell operation using a general computer analysis program for semiconductor devices. The computer program, which simultaneously solves Poisson's equation and the electron and hole quasi-Fermi level equations, provides an accurate numerical solution of solar-cell operation without limiting assumptions or approximations. It is found that minority-carrier lifetime and heavy doping effects in the n(+) surface region present serious limitations to efficiency in low-resistivity silicon solar cells.
We present a hierarchical approach to the 'atomistic' simulation of aggressively scaled sub-0.1 micron MOSFET's. These devices are so small that their characteristics depend on the precise location of dopant atoms within them, not just on their average density. A full-scale three-dimensional drift-diffusion atomistic simulation approach is first described and used to verify more economical, but restricted, options. To reduce processor time and memory requirements at high drain voltage, we have developed a self-consistent option based on a solution of the current continuity equation restricted to a thin slab of the channel. This is coupled to the solution of the Poisson equation in the whole simulation domain in the Gummel iteration cycles. The accuracy of this approach is investigated in comparison to the full self-consistent solution. At low drain voltage, a single solution of the nonlinear Poisson equation is sufficient to extract the current with satisfactory accuracy. In this case, the current is calculated by solving the current continuity equation in a drift approximation only, also in a thin slab containing the MOSFET channel. The regions of applicability for the different components of this hierarchical approach are illustrated in example simulations covering the random dopant-induced threshold voltage fluctuations, threshold voltage lowering, threshold voltage asymmetry, and drain current fluctuations.
Experiments were performed to determine the effect of cyclic loading on bearing cage strength. A long term working tensile load of approximately 1300 N (300 lbs) was found to be the likely maximum. Higher loads caused a decrease in cage tensile strength after the 125,000 cycle testing period. Poisson's ratio in compression was found to be highly dependent upon the direction of the fiberglass plies. At room temperature the value was 0.15 with the plies and 0.68 across the plies. At -196 C (-321 F), the value with the plies was 0.20. The results of the analyses conducted have again demonstrated the critical need for improved lubrication in the high pressure oxygen turbopump bearings. Lubricant films with low shear strength and low friction coefficients promote cage stability and decrease ball/cage forces during marginal operating conditions. The analysis of the effect of combined bearing loads on ball/cage loads has identified a radial load of 3600 N (800 lbs) as the maximum for the current clearance of the balls and cage pockets. Liquid oxygen impinging on the cage in the direction of rotation was found to enhance cage stability.
I have performed a reliability & maintainability analysis for the Amine Swingbed payload system. The Amine Swingbed is a carbon dioxide removal technology that has gone through 2,400 hours of International Space Station on-orbit use between 2013 and 2016. While the Amine Swingbed is currently an experimental payload system, the Amine Swingbed may be converted to system hardware. If the Amine Swingbed becomes system hardware, it will supplement the Carbon Dioxide Removal Assembly (CDRA) as the primary CO2 removal technology on the International Space Station. NASA is also considering using the Amine Swingbed as the primary carbon dioxide removal technology for future extravehicular mobility units and for the Orion, which will be used for the Asteroid Redirect and Journey to Mars missions. The qualitative component of the reliability and maintainability analysis is a Failure Modes and Effects Analysis (FMEA). In the FMEA, I have investigated how individual components in the Amine Swingbed may fail, and what the worst case scenario is should a failure occur. The significant failure effects are the loss of ability to remove carbon dioxide, the formation of ammonia due to chemical degradation of the amine, and loss of atmosphere because the Amine Swingbed uses the vacuum of space to regenerate the Amine Swingbed. In the quantitative component of the reliability and maintainability analysis, I have assumed a constant failure rate for both electronic and nonelectronic parts. Using this data, I have created a Poisson distribution to predict the failure rate of the Amine Swingbed as a whole. I have determined a mean time to failure for the Amine Swingbed to be approximately 1,400 hours. The observed mean time to failure for the system is between 600 and 1,200 hours. This range includes initial testing of the Amine Swingbed, as well as software faults that are understood to be non-critical. If many of the commercial parts were switched to military-grade parts, the expected mean time to failure would be 2,300 hours. Both calculated mean times to failure for the Amine Swingbed use conservative failure rate models. The observed mean time to failure for CDRA is 2,500 hours. Working on this project and for NASA in general has helped me gain insight into current aeronautics missions, reliability engineering, circuit analysis, and different cultures. Prior my internship, I did not have a lot knowledge about the work being performed at NASA. As a chemical engineer, I had not really considered working for NASA as a career path. By engaging in interactions with civil servants, contractors, and other interns, I have learned a great deal about modern challenges that NASA is addressing. My work has helped me develop a knowledge base in safety and reliability that would be difficult to find elsewhere. Prior to this internship, I had not thought about reliability engineering. Now, I have gained a skillset in performing reliability analyses, and understanding the inner workings of a large mechanical system. I have also gained experience in understanding how electrical systems work while I was analyzing the electrical components of the Amine Swingbed. I did not expect to be exposed to as many different cultures as I have while working at NASA. I am referring to both within NASA and the Houston area. NASA employs individuals with a broad range of backgrounds. It has been great to learn from individuals who have highly diverse experiences and outlooks on the world. In the Houston area, I have come across individuals from different parts of the world. Interacting with such a high number of individuals with significantly different backgrounds has helped me to grow as a person in ways that I did not expect. My time at NASA has opened a window into the field of aeronautics. After earning a bachelor's degree in chemical engineering, I plan to go to graduate school for a PhD in engineering. Prior to coming to NASA, I was not aware of the graduate Pathways program. I intend to apply for the graduate Pathways program as positions are opened up. I would like to pursue future opportunities with NASA, especially as my engineering career progresses.
It is shown that a rigorous mathematical treatment of the boundary conditions for rigid no-slip surfaces requires, in general, that two passes be made in solving Poisson's equation. A Fourier solution procedure was the key to this realization. It is also shown that a corresponding method exists for finite and, for the same local approximation, the usual one step method gives the same result for the geometry considered here. The more fundamental two step method gives a new insight into the physical significance of boundary vorticity.
Power-limited communication links may be susceptible to significant degradation if intentional or unintentional pulsed high level radio frequency interference (RFI) is present. Pulsed RFI is, in fact, of current interest to NASA in studies relating to its Tracking and Data Relay Satellite System (TDRSS). The present paper examines the impact of pulsed RFI on the error probability performance of a power-limited satellite communication link: the assumed modulation scheme is PN coded binary PSK. The composite effects of thermal noise, pulsed CW and pulsed Gaussian noise are analyzed, where RFI arrivals are assumed to follow Poisson statistics. Under the assumption that the satellite repeater is ideal and that integrate and dump filtering is employed at the ground receiver, an exact error probability expression and associated approximations are derived. Computed results are generated using an arbitrarily specified RFI model.
NASA Space Shuttle Challenger SIR-B ocean scenes are used to derive directional wave spectra for which speckle noise is modeled as a function of Rayleigh random phase coherence downrange and Poisson random amplitude errors inherent in the Doppler measurement of along-track position. A Fourier filter that preserves SIR-B image phase relations is used to correct the stationary and dynamic response characteristics of the remote sensor and scene correlator, as well as to subtract an estimate of the speckle noise component. A two-dimensional map of sea surface elevation is obtained after the filtered image is corrected for both random and deterministic motions.