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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Enthalpy-Based Thermal Evolution of Loops: II. Improvements to the Model

This paper further develops the zero-dimensional (0D) hydrodynamic coronal loop model "Enthalpy-based Thermal Evolution of Loops" (EBTEL) originally proposed by Klimchuk et al (2008), which studies the plasma response to evolving coronal heating. It has typically been applied to impulsive heating events. The basis of EBTEL is the modelling of mass exchange between the corona and transition region and chromosphere in response to heating variations, with the key parameter being the ratio of transition region to coronal radiation. We develop new models for this parameter that now include gravitational stratification and a physically motivated approach to radiative cooling. A number of examples are presented, including nanoflares in short and long loops, and a small flare. It is found that while the evolution of the loop temperature is rather insensitive to the details of the model, accurate tracking of the density requires the inclusion of our new features. In particular, we are able to now obtain highly over-dense loops in the late cooling phase and decreases to the coronal density arising due to stratification. The 0D results are compared to a 1D hydro code (Hydrad). The agreement is acceptable, with the exception of the flare case where some versions of Hydrad can give significantly lower densities. This is attributed to the method used to model the chromosphere in a flare. EBTEL is suitable for general use as a tool for (a) quick-look results of loop evolution in response to a given heating function and (b) situations where the modelling of hundreds or thousands of elemental loops is needed. A single run takes a few seconds on a contemporary laptop.

Cargill, P. J.↗
NASA NTRS2003

Reduction of Defects in Germanium-Silicon

Crystals grown without contact with a container have far superior quality to otherwise similar crystals grown in direct contact with a container. In addition to float-zone processing, detached- Bridgman growth is a promising tool to improve crystal quality, without the limitations of float zoning or the defects introduced by normal Bridgman growth. Goals of this project include the development of the detached Bridgman process to be reproducible and well understood and to quantitatively compare the defect and impurity levels in crystals grown by these three methods. Germanium (Ge) and germanium-silicon (Ge-Si) alloys are being used. At MSFC, we are responsible for the detached Bridgman experiments intended to differentiate among proposed mechanisms of detachment, and to confirm or refine our understanding of detachment. Because the contact angle is critical to determining the conditions for detachment, the sessile drop method was used to measure the contact angles as a function of temperature and composition for a large number of substrates made of potential ampoule materials. Growth experiments have used pyrolytic boron nitride (pBN) and fused silica ampoules with the majority of the detached results occurring predictably in the pBN. Etch pit density (EPD) measurements of normal and detached Bridgman-grown Ge samples show a two order of magnitude improvement in the detached-grown samples. The nature and extent of detachment is determined by using profilometry in conjunction with optical and electron microscopy. The stability of detachment has been analyzed, and an empirical model for the conditions necessary to achieve sufficient stability to maintain detached growth for extended periods has been developed. We have investigated the effects on detachment of ampoule material, pressure difference above and below the melt, and Si concentration; samples that are nearly completely detached can be grown repeatedly in pBN. Current work is concentrated on developing a method to make in situ pressure measurements in the growth ampoules.

Source record↗
NASA NTRS1974

A solution to the problem of optimizing the fuel bias for a liquid propellant rocket by an application of the central limit theorem

A method of determining the fuel bias for a bipropellant liquid rocket that minimizes outage associated penalties on payload potential is presented. A fuel bias so derived is normally called the optimum fuel bias. The subjects discussed are: (1) probability density function of outage, (2) computer program listing, and (3) choosing the optimum fuel bias.

Viera, W. J.↗
NASA NTRS1981

Interpretation of phase velocity measurements of wind-generated surface waves

Differences in experimental results for surface gravity wave phase speeds obtained by cross spectra and slope-height spectra are compared. It is shown that, for a simple two-dimensional model containing both dispersive and non-dispersive waves, the observed differences can be explained readily. The slope-height technique measures a weighted average of the wavenumber squared; for frequencies sufficiently large compared with the dominant frequency, the computed phase-speed variation with frequency is representative of the wave having the larger wavenumber - the dispersive wave. For the cross-spectral method, it is shown that the small-probe-separation assumption usually employed is not valid for frequencies much larger than the dominant frequency, except at the singular point where both waves have exactly equal spectral densities and the phase function represents an average of the two modes. For all other cases, the phase function approaches that of the wave having the larger spectral density, and essentially ignores the presence of the other wave, even if both modes are relatively close in their contributions to the total spectral density.

Dudis, J. J.↗
NASA NTRS1990

High speed turbulent reacting flows: DNS and LES

Work on understanding the mechanisms of mixing and reaction in high speed turbulent reacting flows was continued. Efforts, in particular, were concentrated on taking advantage of modern computational methods to simulate high speed turbulent flows. In doing so, two methodologies were used: large eddy simulations (LES) and direct numerical simulations (DNS). In the work related with LES the objective is to study the behavior of the probability density functions (pdfs) of scalar properties within the subgrid in reacting turbulent flows. The data base obtained by DNS for a detailed study of the pdf characteristics within the subgrid was used. Simulations are performed for flows under various initializations to include the effects of compressibility on mixing and chemical reactions. In the work related with DNS, a two-dimensional temporally developing high speed mixing layer under the influence of a second-order non-equilibrium chemical reaction of the type A + B yields products + heat was considered. Simulations were performed with different magnitudes of the convective Mach numbers and with different chemical kinetic parameters for the purpose of examining the isolated effects of the compressibility and the heat released by the chemical reactions on the structure of the layer. A full compressible code was developed and utilized, so that the coupling between mixing and chemical reactions is captured in a realistic manner.

Givi, Peyman↗
NASA NTRS1991

The analysis of apparent optical depth profiles for interstellar absorption lines

Attention is given to the apparent optical depth method, a procedure for analyzing interstellar absorption lines. Observed absorption-line profiles are converted into profiles of apparent optical depth, and apparent column density per unit velocity. By comparing the latter for a given interstellar species having two or more absorption lines which differ in the product, the presence or absence of unresolved saturated structure in the profiles can be directly inferred. The method is illustrated using absorption-line data from the IUE satellite for the highly ionized gas toward HD 64760. Additional illustrations and a study of the method's accuracy are provided through a series of numerical simulations of multicomponent interstellar absorption situations. The method is compared to the standard growth curve method for deriving interstellar column densities. The principal value of the apparent optical depth method is that the absorption-line data are directly converted into a form that provides for direct scientific interpretations of the physical conditions in the interstellar absorbing medium as a function of velocity.

Savage, Blair D.↗
NASA NTRS1992

FDDO and DSMC analyses of rarefied gas flow through 2D nozzles

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas expanding through a two-dimensional nozzle and into a surrounding low-density environment. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, the variable hard sphere model is used as a molecular model and the no time counter method is employed as a collision sampling technique. The results of both the FDDO and the DSMC methods show good agreement. The FDDO method requires less computational effort than the DSMC method by factors of 10 to 40 in CPU time, depending on the degree of rarefaction.

Chung, Chan-Hong↗
NASA NTRS1998

Pulse Shaped 8-PSK Bandwidth Efficiency and Spectral Spike Elimination

The most bandwidth-efficient communication methods are imperative to cope with the congested frequency bands. Pulse shaping methods have excellent effects on narrowing bandwidth and increasing band utilization. The position of the baseband filters for the pulse shaping is crucial. Post-modulation pulse shaping (a low pass filter is located after the modulator) can change signals from constant envelope to non-constant envelope, and non-constant envelope signals through non-linear device (a SSPA or TWT) can further spread the power spectra. Pre-modulation pulse shaping (a filter is located before the modulator) will have constant envelope. These two pulse shaping methods have different effects on narrowing the bandwidth and producing bit errors. This report studied the effect of various pre-modulation pulse shaping filters with respect to bandwidth, spectral spikes and bit error rate. A pre-modulation pulse shaped 8-ary Phase Shift Keying (8PSK) modulation was used throughout the simulations. In addition to traditional pulse shaping filters, such as Bessel, Butterworth and Square Root Raised Cosine (SRRC), other kinds of filters or pulse waveforms were also studied in the pre-modulation pulse shaping method. Simulations were conducted by using the Signal Processing Worksystem (SPW) software package on HP workstations which simulated the power spectral density of pulse shaped 8-PSK signals, end to end system performance and bit error rates (BERS) as a function of Eb/No using pulse shaping in an AWGN channel. These results are compared with the post-modulation pulse shaped 8-PSK results. The simulations indicate traditional pulse shaping filters used in pre-modulation pulse shaping may produce narrower bandwidth, but with worse BER than those in post-modulation pulse shaping. Theory and simulations show pre- modulation pulse shaping could also produce discrete line power spectra (spikes) at regular frequency intervals. These spikes may cause interference with adjacent channel and reduce power efficiency. Some particular pulses (filters), such as trapezoid and pulses with different transits (such as weighted raised cosine transit) were found to reduce bandwidth and not generate spectral spikes. Although a solid state power amplifier (SSPA) was simulated in the non-linear (saturation) region, output power spectra did not spread due to the constant envelope 8-PSK signals.

Tao, Jian-Ping↗
NASA NTRS1999

3D-Stereoscopic Analysis of Solar Active Region Loops: SoHo/EIT Observations at Temperatures of 1.5-2.5 MK - 2

In this paper we study the three-dimensional (3D) structure of hot (T(sub e) approximately equals 1.5 - 2.5 MK) loops in solar active region NOAA 7986, observed on 1996 August 30 with the Extreme-ultraviolet Imaging Telescope (EIT) onboard the Solar and Heliospheric Observatory (SoHO). This complements a first study on cooler (T(sub e) approximately equals 1.0 - 1.5 MK) loops of the same active region, using the same method of Dynamic Stereoscopy to reconstruct the 3D geometry. We reconstruct the 3D-coordinates x(s), y(s), z(s), the density n(sub e)(s), and temperature profile T(sub e)(s) of 35 individual loop segments (as function of the loop coordinate s) using EIT 195 A and 284 A images. The major findings are: (1) All loops are found to be in hydrostatic equilibrium, in the entire temperature regime of T(sub e) = 1.0 - 2.5 MK; (2) The analyzed loops have a height of 2-3 scale heights, and thus only segments extending over about one vertical scale height have sufficient emission measure contrast for detection; (3) The temperature gradient over the lowest scale height is of order dT/ds is approximately 1 - 4 K/km; (4) The radiative loss rate is found to exceed the conductive loss rate by about two orders or magnitude, making thermal conduction negligible to explain the temperature structure of the loops; (5) A steady-state can only be achieved when the heating rate E(sub H) matches the radiative loss rate in hydrostatic equilibrium, requiring a heat deposition length lambda(sub H) of the half density scale height lambda, predicting a scaling law with the loop base pressure, EH varies as p(sub 0 exp 2). This favors coronal heating mechanisms that operate near the loop footpoints; (6) We find a reciprocal correlation between the loop pressure p(sub 0) and loop length L, i.e. p(sub 0) varies as 1/L, implying a scaling law of the steady-state requirement with loop length, i.e. E(sub H ) varies as 1/L(exp 2). The heating rate shows no correlation with the loop-aligned magnetic field component B(sub z) at the footpoints, but is correlated with the azimuthal field B(sub phi) = Bz(RDelta Phi/L) of a twisted loop, and is thus consistent with heating mechanisms based on field-aligned currents.

Aschwanden, Markus J.↗
NASA NTRS1975

Infrared forbidden lines in H II regions and planetary nebulae

Emissivity coefficients are calculated as functions of electron density and temperature for the IR forbidden lines between 2 and 300 microns of the elements and ionization stages most abundant in H II regions and planetary nebulae. The effect of self-absorption is investigated, and a method is presented for estimating when self-absorption may be important. The IR line spectrum of the Orion Nebula is predicted using the emissivity coefficients, and ionic abundances are derived for the Ar III, S IV, and Ne II forbidden lines in a number of planetary nebulae. The sulfur abundance and sulfur/oxygen ratio are discussed for H II regions and planetary nebulae. The average sulfur/oxygen ratio is found to be about 0.05 for gaseous nebulae.

Simpson, J. P.↗
NASA NTRS1986

Practical method for highly accurate large-scale surface calculations

An accurate and efficient film linearized muffin-tin orbital (FLMTO) technique for surface electronic-structure calculations is presented which uses only 60-70 basis functions, as opposed to the 300 functions used in the linear augmented plane-wave method. Calculations for three different (3d and 4d) transition-metal films resulted in high quality results for five-layer slabs of Cu(001), Fe(001), and Ru(001), in addition to good results for the work functions and projected density of states. By retaining the LMTO small basis size, computer time and memory are reduced, making practical the study of systems with a larger number of atoms in the two-dimensional unit cell.

Fernando, G. W.↗
NASA NTRS2017

Association Between Vascular Density and Loss of Protective RAS During Early NPDR by Fractal Dimension

Purpose: Our hypothesis predicts that blood vessels within the retina increase in density during early-stage nonproliferative diabetic retinopathy (NPDR), based on previous results of a small retrospective study. For the current prospective study, the remodeling of arteries and veins during progression of early NPDR is assessed by a repertoire of parameters that includes the fractal dimension (D(sub f) ). In complex structures such as branching vascular trees, D(sub f) is a sensitive measure of space-filing capacity. The renin-angiotensin system (RAS) is implicated in DR pathogenesis and the function of circulating angiogenic cells (CACs), a critical bone marrow-derived population instrumental in vascular repair. Methods: Arterial and venous branching patterns were extracted from images of 6 normal controls and 3 early NPDR subjects (2 moderate, 1 mild) acquired by Heidelberg Spectralis (Registered Trademark) OCT following fluorescein angiography (FA). The vascular branching patterns were analyzed by NASAs VESsel GENeration Analysis (VESGEN) software, in which skeletonized representations were generated automatically to yield D(sub f) by the box-counting method. For binary 2D images, D(sub f) varies between limiting Euclidean dimensions of 1 and 2. Peripheral blood of diabetics and controls was collected for CD34+ CAC isolation. The gene expression of RAS in CACs was assessed by qPCR for Mas receptor to Ang-(1-7). The vasoreparative function of the CACs was measured by migration ability toward CXCL12 (SDF-1). Results: By D(sub f), venous and arterial densities were 1.370 +/- 0.006 and 1.329 +/- 0.016 for early NPDR, compared to 1.318 +/- 0.012 and 1.320 +/- 0.036 for control. The space filling capacity in early NPDR measured by D(sub f), a sensitive parameter, therefore demonstrated a pronounced increase for veins, but not for arteries. Mas receptor mRNA in CACs was increased in diabetics without DR but reduced with onset of NPDR, indicating possible loss of compensation of protective RAS during early DR. Migratory dysfunction of CD34+ cells was further associated with DR. Conclusions: As assessed by the fractal dimension in our preliminary study, the space-filling capacity of veins, but not arteries, was greater in early NPDR than in control. Larger patient populations will be examined as we complete our ongoing longitudinal study. Results further suggest the protective RAS axis within diabetic CACs is lost early in DR and is associated with increased vascular remodeling as evidenced by VESGEN analysis.

vascular↗
NASA NTRS2004

On Bayesian Inductive Inference & Predictive Estimation

We investigate Bayesian inference and the Principle of Maximum Entropy (PME) as methods for doing inference under uncertainty. This investigation is primarily through concrete examples that have been previously investigated in the literature. We find that it is possible to do Bayesian inference and PME inference using the same information, despite claims to the contrary, but that the results are not directly comparable. This is because Bayesian inference yields a probability density function (pdf) over the unknown model parameters, whereas PME yields point estimates. If mean estimates are extracted from the Bayesian pdfs, the resulting parameter estimates can differ radically from the PME values and also from the Maximum Likelihood values. We conclude that these differences are due to the Bayesian inference not assuming anything beyond the given prior probabilities and the data, whereas PME implicitly assumes that the given constraints are the only constraints that are operating. Since this assumption can be wrong, PME values may have to be revised when subsequent data shows evidence for more constraints. The entropy concentration previously "proved" by E. T. Jaynes is shown to be in error. Further, we show that PME is a generalized form of independence assumption, and so can be a very powerful method of inference when the variables being investigated are largely independent of each other.

Cheeseman, Peter↗
NASA NTRS1981

Statistics of some atmospheric turbulence records relevant to aircraft response calculations

Methods for characterizing atmospheric turbulence are described. The methods illustrated include maximum likelihood estimation of the integral scale and intensity of records obeying the von Karman transverse power spectral form, constrained least-squares estimation of the parameters of a parametric representation of autocorrelation functions, estimation of the power spectra density of the instantaneous variance of a record with temporally fluctuating variance, and estimation of the probability density functions of various turbulence components. Descriptions of the computer programs used in the computations are given, and a full listing of these programs is included.

Mark, W. D.↗
NASA NTRS1990

Linear Calibration of Radiographic Mineral Density Using Video-Digitizing Methods

Radiographic images can provide quantitative as well as qualitative information if they are subjected to densitometric analysis. Using modern video-digitizing techniques, such densitometry can be readily accomplished using relatively inexpensive computer systems. However, such analyses are made more difficult by the fact that the density values read from the radiograph have a complex, nonlinear relationship to bone mineral content. This article derives the relationship between these variables from the nature of the intermediate physical processes, and presents a simple mathematical method for obtaining a linear calibration function using a step wedge or other standard.

Martin, R. Bruce↗
NASA NTRS1990

Linear Calibration of Radiographic Mineral Density Using Video-Digitizing Methods

Radiographic images can provide quantitative as well as qualitative information if they are subjected to densitometric analysis. Using modem video-digitizing techniques, such densitometry can be readily accomplished using relatively inexpensive computer systems. However, such analyses are made more difficult by the fact that the density values read from the radiograph have a complex, nonlinear relationship to bone mineral content. This article derives the relationship between these variables from the nature of the intermediate physical processes, and presents a simple mathematical method for obtaining a linear calibration function using a step wedge or other standard.

Martin, R. Bruce↗
NASA NTRS2014

Reconstructing the Sky Location of Gravitational-Wave Detected Compact Binary Systems: Methodology for Testing and Comparison

The problem of reconstructing the sky position of compact binary coalescences detected via gravitational waves is a central one for future observations with the ground-based network of gravitational-wave laser interferometers, such as Advanced LIGO and Advanced Virgo. Different techniques for sky localization have been independently developed. They can be divided in two broad categories: fully coherent Bayesian techniques, which are high latency and aimed at in-depth studies of all the parameters of a source, including sky position, and "triangulation-based" techniques, which exploit the data products from the search stage of the analysis to provide an almost real-time approximation of the posterior probability density function of the sky location of a detection candidate. These techniques have previously been applied to data collected during the last science runs of gravitational-wave detectors operating in the so-called initial configuration. Here, we develop and analyze methods for assessing the self consistency of parameter estimation methods and carrying out fair comparisons between different algorithms, addressing issues of efficiency and optimality. These methods are general, and can be applied to parameter estimation problems other than sky localization. We apply these methods to two existing sky localization techniques representing the two above-mentioned categories, using a set of simulated inspiralonly signals from compact binary systems with a total mass of equal to or less than 20M solar mass and nonspinning components. We compare the relative advantages and costs of the two techniques and show that sky location uncertainties are on average a factor approx. equals 20 smaller for fully coherent techniques than for the specific variant of the triangulation-based technique used during the last science runs, at the expense of a factor approx. equals 1000 longer processing time.

technique used↗