Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY FUNCTIONAL METHOD”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1961

ANALYSIS OF CRYOGENIC PROPELLANT FEED SYSTEMS FOR ELECTROTHERMAL ENGINES.

An analytical investigation was made of the problems associated with the long-term storage and control of flow of cryogenic propellants for electrothermal engines. Storage of 35 to 174 pounds of propellant for three years is made difficult by the high ratio of surface area to storage tank volume, and the resulting high rate of heat leak. Vented and non-vented storage, storage as high pressure gas, and the use of refrigeration were considered. Recommended storage methods are: for ethane and ammonia vented storage - by venting propellant stored in excess of mission requirements; for hydrogen and helium - by non-vented storage of cold supercritical vapor; and for methane by storage as a non-vented two-phase mixture. The optimum propellant density in the storage tank and optimum insulation thickness were also determined. Design charts were developed for hydrogen and helium, for estimating the maximum pressure expected in non-vented storage tanks, as function of mission requirements. Separation of the constituent phases of liquid-vapor mixtures is achieved by throttling the two-phase mixture to a lower pressure (and temperature) and vaporizing the liquid in the mixture by heat exchange with the contents of the storage tank. The time limitations for continuous use of this method are discussed. Linde super insulations SI-44 or SI -91 are recommended for insulation. Precooling the insulation to temperatures below the steady state level with cold helium vapor before or during fueling of the storage tank eliminates the flow of heat to the tank contents during countdown and take-off. The masses, sizes, and power requirements of the recommended systems are low. The controls are simple and reliable. From the equations and charts in this report the components of the system for a particular mission can be sized. Recommendations are given for experimental work required for the development and testing of such components.

ELECTROTHERMAL ENGINE↗
NASA NTRS1988

Electron density-sensitive line ratios for Fe XII in the solar corona

The collision strengths for electron-impact excitation of fine-structure transitions from the ground 3s23p3 configuration to the excited 3s3p4 configuration in Fe XII are presented. These are calculated by transforming the K-matrices obtained in a nonrelativistic R-matrix method to K-matrices in intermediate coupling using term-coupling coefficients which take into account the effect of breakdown of target LS coupling. Extensive configuration-interaction wave functions are used to represent the target states. The effective collision strengths are determined assuming a Maxwellian distribution of electron energies. Rydberg series of resonances converging to the excited state thresholds are explicitly included in the calculation. The relative level populations and emission-line strengths for Fe XII are calculated using the new atomic data presented in this paper and that published by Tayal, Henry, and Pradhan (1987) and Tayal and Henry (1986) over a wide range of electron densities and temperatures (1-2 x 10 to the 6th K). Present results for collision strengths and emission-line strengths are compared with the previous calculation of Flower (1977). Some significant differences are noted.

Tayal, S. S.↗
NASA NTRS1993

An X-ray method for detecting substructure in galaxy clusters - Application to Perseus, A2256, Centaurus, Coma, and Sersic 40/6

We use the moments of the X-ray surface brightness distribution to constrain the dynamical state of a galaxy cluster. Using X-ray observations from the Einstein Observatory IPC, we measure the first moment FM, the ellipsoidal orientation angle, and the axial ratio at a sequence of radii in the cluster. We argue that a significant variation in the image centroid FM as a function of radius is evidence for a nonequilibrium feature in the intracluster medium (ICM) density distribution. In simple terms, centroid shifts indicate that the center of mass of the ICM varies with radius. This variation is a tracer of continuing dynamical evolution. For each cluster, we evaluate the significance of variations in the centroid of the IPC image by computing the same statistics on an ensemble of simulated cluster images. In producing these simulated images we include X-ray point source emission, telescope vignetting, Poisson noise, and characteristics of the IPC. Application of this new method to five Abell clusters reveals that the core of each one has significant substructure. In addition, we find significant variations in the orientation angle and the axial ratio for several of the clusters.

Mohr, Joseph J.↗
NASA NTRS2004

Using Redundancy To Reduce Errors in Magnetometer Readings

A method of reducing errors in noisy magnetic-field measurements involves exploitation of redundancy in the readings of multiple magnetometers in a cluster. By "redundancy"is meant that the readings are not entirely independent of each other because the relationships among the magnetic-field components that one seeks to measure are governed by the fundamental laws of electromagnetism as expressed by Maxwell's equations. Assuming that the magnetometers are located outside a magnetic material, that the magnetic field is steady or quasi-steady, and that there are no electric currents flowing in or near the magnetometers, the applicable Maxwell 's equations are delta x B = 0 and delta(raised dot) B = 0, where B is the magnetic-flux-density vector. By suitable algebraic manipulation, these equations can be shown to impose three independent constraints on the values of the components of B at the various magnetometer positions. In general, the problem of reducing the errors in noisy measurements is one of finding a set of corrected values that minimize an error function. In the present method, the error function is formulated as (1) the sum of squares of the differences between the corrected and noisy measurement values plus (2) a sum of three terms, each comprising the product of a Lagrange multiplier and one of the three constraints. The partial derivatives of the error function with respect to the corrected magnetic-field component values and the Lagrange multipliers are set equal to zero, leading to a set of equations that can be put into matrix.vector form. The matrix can be inverted to solve for a vector that comprises the corrected magnetic-field component values and the Lagrange multipliers.

Kulikov, Igor↗
NASA NTRS1972

A comparative study of nonparametric methods for pattern recognition

The applied research discussed in this report determines and compares the correct classification percentage of the nonparametric sign test, Wilcoxon's signed rank test, and K-class classifier with the performance of the Bayes classifier. The performance is determined for data which have Gaussian, Laplacian and Rayleigh probability density functions. The correct classification percentage is shown graphically for differences in modes and/or means of the probability density functions for four, eight and sixteen samples. The K-class classifier performed very well with respect to the other classifiers used. Since the K-class classifier is a nonparametric technique, it usually performed better than the Bayes classifier which assumes the data to be Gaussian even though it may not be. The K-class classifier has the advantage over the Bayes in that it works well with non-Gaussian data without having to determine the probability density function of the data. It should be noted that the data in this experiment was always unimodal.

Hahn, S. F.↗
NASA NTRS1975

Two dimensional thermal and charge mapping of power thyristors

The two dimensional static and dynamic current density distributions within the junction of semiconductor power switching devices and in particular the thyristors were obtained. A method for mapping the thermal profile of the device junctions with fine resolution using an infrared beam and measuring the attenuation through the device as a function of temperature were developed. The results obtained are useful in the design and quality control of high power semiconductor switching devices.

Hu, S. P.↗
NASA NTRS1975

Satellite disintegration dynamics

The subject of satellite disintegration is examined in detail. Elements of the orbits of individual fragments, determined by DOD space surveillance systems, are used to accurately predict the time and place of fragmentation. Dual time independent and time dependent analyses are performed for simulated and real breakups. Methods of statistical mechanics are used to study the evolution of the fragment clouds. The fragments are treated as an ensemble of non-interacting particles. A solution of Liouville's equation is obtained which enables the spatial density to be calculated as a function of position, time and initial velocity distribution.

Dasenbrock, R. R.↗
NASA NTRS1983

Modeling of steady, rotational, transonic winds from rotating stars and galaxies

The theory of steady transonic winds from condensed bodies is extended to general, two-dimensional, axisymmetric systems. A stream function is used to reduce the gasdynamics equations to a single, second-order differential equation plus an algebraic equation for the density. The approach extends Parker's (1958) quasi-one-dimensional theory, which uses Bernoulli's theorem, to a complete two-dimensional calculation which includes vorticity and rotation. The conceptual basis of the stream function approach is described in detail, and a numerical method for solving the resulting equations is presented. The applications illustrate solutions for stellar (spherical source surface) and galactic (oblate spheroidal source surface) models. Among other things, it is found that for rapidly rotating stellar models the Coriolis force dominates the centrifugal terms and the streamlines bend toward the rotation axis as a consequence of the axisymmetry. For both the rotating and nonrotating galaxy models, the streamlines bend toward the equator.

Kopriva, D. A.↗
NASA NTRS1994

Nonlinear calculations of the time evolution of black hole accretion disks

Based on previous works on black hole accretion disks, I continue to explore the disk dynamics using the finite difference method to solve the highly nonlinear problem of time-dependent alpha disk equations. Here a radially zoned model is used to develop a computational scheme in order to accommodate functional dependence of the viscosity parameter alpha on the disk scale height and/or surface density. This work is based on the author's previous work on the steady disk structure and the linear analysis of disk dynamics to try to apply to x-ray emissions from black candidates (i.e., multiple-state spectra, instabilities, QPO's, etc.).

Luo, C.↗
NASA NTRS1994

Probabilistic Simulation of Stress Concentration in Composite Laminates

A computational methodology is described to probabilistically simulate the stress concentration factors (SCF's) in composite laminates. This new approach consists of coupling probabilistic composite mechanics with probabilistic finite element structural analysis. The composite mechanics is used to probabilistically describe all the uncertainties inherent in composite material properties, whereas the finite element is used to probabilistically describe the uncertainties associated with methods to experimentally evaluate SCF's, such as loads, geometry, and supports. The effectiveness of the methodology is demonstrated by using is to simulate the SCF's in three different composite laminates. Simulated results match experimental data for probability density and for cumulative distribution functions. The sensitivity factors indicate that the SCF's are influenced by local stiffness variables, by load eccentricities, and by initial stress fields.

Chamis, C. C.↗
NASA NTRS2003

Electromagnetic Scattering From a Polygonal Thin Metallic Plate Using Quadrilateral Meshing

The problem of electromagnetic (EM) scattering from irregularly shaped, thin, metallic flat plates in free space is solved using the electric field integral equation (EFIE) approach in conjunction with the method of moments (MoM) with quadrilateral meshing. An irregularly shaped thin plate is discretized into quadrilateral patches and the unknown electric surface current over the plate is expressed in terms of proper basis functions over these patches. The basis functions for the electric surface current density that satisfy the proper boundary conditions on these quadrilateral patches are derived. The unknown surface current density on these quadrilateral patches is determined by setting up and solving the electric field integral equation by the application of the MoM. From the knowledge of the surface current density, the EM scattering from various irregularly shaped plates is determined and compared with the earlier published results. The novelty in the present approach is the use of quadrilateral patches instead of well known and often used triangular patches. The numerical results obtained using the quadrilateral patches compare favorably with measured results.

Deshpande, Manohar D.↗
NASA NTRS2017

Planetary Probe Entry Atmosphere Estimation Using Synthetic Air Data System

This paper develops an atmospheric state estimator based on inertial acceleration and angular rate measurements combined with an assumed vehicle aerodynamic model. The approach utilizes the full navigation state of the vehicle (position, velocity, and attitude) to recast the vehicle aerodynamic model to be a function solely of the atmospheric state (density, pressure, and winds). Force and moment measurements are based on vehicle sensed accelerations and angular rates. These measurements are combined with an aerodynamic model and a Kalman-Schmidt filter to estimate the atmospheric conditions. The new method is applied to data from the Mars Science Laboratory mission, which landed the Curiosity rover on the surface of Mars in August 2012. The results of the new estimation algorithm are compared with results from a Flush Air Data Sensing algorithm based on onboard pressure measurements on the vehicle forebody. The comparison indicates that the new proposed estimation method provides estimates consistent with the air data measurements, without the use of pressure measurements. Implications for future missions such as the Mars 2020 entry capsule are described.

Karlgaard, Chris↗
NASA NTRS2018

Behavior of Langmuir Probes in Non-Equilibrium Plasmas

Langmuir probes are diagnostic tools used to determine electron temperature, number density, and plasma potential. Irving Langmuir first used an electrostatic probe in the 1920s to find these characteristics in ionized gases. Single, double, and triple Langmuir probes are commonly used in plasma diagnostics because of their relative simplicity. In the single probe, a swept voltage is applied between the probe tip and circuit common to acquire a waveform showing the collected current as a function of applied voltage. A double Langmuir probe consists of two tips, both inserted into the plasma, with a voltage applied between them. As this voltage is swept, a current-voltage characteristic is measured. In a triple probe three probe tips are electrically coupled to each other with constant non-swept voltages applied between each of the tips. The voltages are selected to represent three points on the single Langmuir probe I-V curve. Elimination of the voltage sweep makes it possible to measure time-varying plasma properties in transient plasmas. Triple Langmuir probe measurements have been widely employed for various types of plasmas, including pulsed and time-varying plasmas such as those seen in pulsed plasma thrusters (PPTs), dense plasma focus devices, plasma flows, and fusion experiments. The typical Langmuir probe analysis for determining electron temperature and number density of the plasma (for a single, double, or triple Langmuir probe) includes an assumption that the plasma is in thermal equilibrium. While the this assumption may be justified for some applications, it is unlikely that it is fully justifiable for pulsed and time-varying plasmas or for the entire time a plasma device is in use. In the present work, we model the responses of Langmuir probes as they are inserted into a range of simple equilibrium and non-equilibrium plasmas. We return to basic governing equations of probe current collection and compute the current to the probes for a distribution function consisting of two Maxwellian distributions with different temperatures (the two-temperature Maxwellian). A variation of this method is also employed, where one of the Maxwellians is offset from zero (in velocity space) to add a suprathermal beam of electrons to the tail of the main Maxwellian distribution (the bump-on-the-tail distribution function). For a range of parameters in these non-Maxwellian distributions, we compute the current collection to the probes. Comparing the distribution function that was assumed a priori with the plasma density and temperature one would infer when applying standard probe theory to analyze the collected currents serves to illustrate the effect a non- Maxwellian plasma would have on results interpreted using the equilibrium probe current collection theory, allowing us to state the magnitudes of these deviations as a function of the assumed distribution function properties.

Polzin, Kurt A.↗
NASA NTRS2023

Estimating the Subsolar Magnetopause Position from Soft X-Ray Images Using a Low-Pass Image Filter

The Lunar Environment heliospheric X-ray Imager (LEXI) and Solar wind Magnetosphere Ionosphere Link Explorer (SMILE) missions will image the Earth’s dayside magnetopause and cusps in soft X-rays after their respective launches in the near future, to specify global magnetic reconnection modes for varying solar wind conditions. To support the success of these scientific missions, it is critical to develop techniques that extract the magnetopause locations from the observed soft X-ray images. In this research, we introduce a new geometric equation that calculates the subsolar magnetopause position ( R s ) from a satellite position, the look direction of the instrument, and the angle at which the X-ray emission is maximized. Two assumptions are used in this method: (1) The look direction where soft X-ray emissions are maximized lies tangent to the magnetopause, and (2) the magnetopause surface near the subsolar point is almost spherical and thus R s is nearly equal to the radius of the magnetopause curvature. We create synthetic soft X-ray images by using the Open Geospace General Circulation Model (OpenGGCM) global magnetohydrodynamic model, the galactic background, the instrument point spread function, and Poisson noise. We then apply the fast Fourier transform and Gaussian low-pass filters to the synthetic images to remove noise and obtain accurate look angles for the soft X-ray peaks. From the filtered images, we calculate R 2 and its accuracy for different LEXI locations, look directions, and solar wind densities by using the OpenGGCM subsolar magnetopause location as ground truth. Our method estimates R s with an accuracy of <0.3 R E when the solar wind density exceeds >10 cm -3 . The accuracy improves for greater solar wind densities and during southward interplanetary magnetic fields. The method captures the magnetopause motion during southward interplanetary magnetic field turnings. Consequently, the technique will enable quantitative analysis of the magnetopause motion and help reveal the dayside reconnection modes for dynamic solar wind conditions. This technique will support the LEXI and SMILE missions in achieving their scientific objectives.

Hyangpyo Kim↗
NASA NTRS1982

Solution to the backward-Kolmogorov equation for a nonstationary oscillation problem

The transition probability density function of a Markovian approximation of the response amplitude of an oscillator under nonstationary excitation is determined in an analytical form. A solution is presented for the associated backward-Kolmogorov equation by transforming the equation into a form amenable to solution by the method of the separation of variables. This procedure results in a boundary value problem which is then solved by using an infinite series of Laguerre polynomials. It is found that the infinite series solution is equivalent to a closed-form solution involving a Bessel function.

Solomos, G. P.↗
NASA NTRS1987

Effect of an electron scattering cloud on X-ray oscillations produced by beaming

The effect of a scattering cloud on the amplitude of oscillations produced by a rotating beam of X-rays is investigated using analytical and Monte Carlo methods. The scattering cloud was modeled as a uniform density sphere, and the source was represented as an anistropic distribution of radiation emerging from a point at the center of the scattering cloud. The intensity distribution produced by the source beam is examined as a function of optical depth. The relation between electron scattering optical depth and the forward-backward ratio is studied. It is observed that the scattering in a central corona of various optical depths reduces the amplitude of the oscillation. The data suggest that the quasi-periodic oscillations observed in the X-ray intensities of some luminous low-mass X-ray binaries are caused by oscillations in the luminosity of the X-ray star.

Brainerd, J.↗
NASA NTRS1998

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

VanOverbeke, Thomas J.↗
NASA NTRS2023

Optimal Estimation Inversion of Ionospheric Electron Density from GNSS-POD Limb Measurements: Part I-Algorithm and Morphology

GNSS-LEO radio links from Precise Orbital Determination (POD) and Radio Occultation (RO) antennas have been used increasingly in characterizing the global 3D distribution and variability of ionospheric electron density (N e ). In this study, we developed an optimal estimation (OE) method to retrieve N e profiles from the slant total electron content (hTEC) measurements acquired by the GNSS-POD links at negative elevation angles (ε < 0°). Although both OE and onion-peeling (OP) methods use the Abel weighting function in the N e inversion, they are significantly different in terms of performance in the lower ionosphere. The new OE results can overcome the large N e oscillations, sometimes negative values, seen in the OP retrievals in the E-region ionosphere. In the companion paper in this Special Issue, the HmF2 and NmF2 from the OE retrieval are validated against ground-based ionosondes and radar observations, showing generally good agreements in NmF2 from all sites. Nighttime hmF2 measurements tend to agree better than the daytime when the ionosonde heights tend to be slightly lower. The OE algorithm has been applied to all GNSS-POD data acquired from the COSMIC-1 (2006–2019), COSMIC-2 (2019–present), and Spire (2019–present) constellations, showing a consistent ionospheric N e morphology. The unprecedented spatiotemporal sampling of the ionosphere from these constellations now allows a detailed analysis of the frequency–wavenumber spectra for the N e variability at different heights. In the lower ionosphere (~150 km), we found significant spectral power in DE1, DW6, DW4, SW5, and SE4 wave components, in addition to well-known DW1, SW2, and DE3 waves. In the upper ionosphere (~450 km), additional wave components are still present, including DE4, DW4, DW6, SE4, and SW4. The co-existence of eastward- and westward-propagating wave4 components implies the presence of a stationary wave4 (SPW4), as suggested by other earlier studies. Further improvements to the OE method are proposed, including a tomographic inversion technique that leverages the asymmetric sampling about the tangent point associated with GNSS-LEO links.

GNSS constellation↗