Search NASA⌕ Search

SEARCH · Search NASA

Results for “First principles calculations”

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 145 records · Page 8

Constraints on the Recent Rate of Lunar Regolith Accumulation from Diviner Observations

Many large craters on the lunar nearside show radar CPR signatures consistent with the presence of blocky ejecta blankets, to distances pre dicted to be covered by continuous ejecta. However, most of these sur faces show limited enhancements in both derived rock abundance and rock-free regolith temperatures calculated from Diviner nighttime infrar ed observations. This indicates that the surface rocks are covered by a layer of thermally insulating regolith material. By matching the results of one-dimensional thermal models to Diviner nighttime temperat ures, we have constrained the thermophysical properties of the upper regolith, and the thickness of regolith overlying proximal ejecta. We find that for all of the regions surveyed (all in the nearside highla nds), the nighttime cooling curves are best fit by a density profile that varies exponentially with depth, consistent with a linear mixture of rocks and regolith fines, with increasing rock content with depth . Our results show significant spatial variations in the density e-folding depth, H, among young crater ejecta regions, indicating differen ces in the thickness of accumulated regolith. However, away from youn g craters, the average regional "equilibrium" value of H (Heq) is remarkably consistent, and is on the order of 5 cm. As expected, near-rim ejecta associated with young craters show lower values of H, indicating a high rock content in the shallow subsurface; for older craters, the average value of H approaches the regional value of Heq. Calculat ed H values for young craters show a clear correlation with published ages, providing the first observational constraint on the recent rate of lunar regolith accumulation. In addition, this result may help to resolve the apparent discrepancy between ages calculated from small crater counts on melt ponds versus counts on continuous ejecta (e.g., King crater; Ashley et al., 2011, LPSC 42, abstract 2437). This method could, in principle, be extended to other airless bodies (e.g., aste roids), which would in turn constrain the recent impactor flux.

Ghent, R. R.↗

Principles of equivalence, Eotvos experiments, and gravitational red-shift experiments - The free fall of electromagnetic systems to post-post-Coulombian order

Free fall in an external static gravitational field of a composite test body consisting of electromagnetically interacting charged particles is analyzed using the TH-epsilon-mu formalism developed by Lightman and Lee (1973). The Lightman and Lee calculation that was carried out to 'post-Coulombian' order is extended to 'post-post-Coulombian' order, and violations of the weak equivalence principle are found which result from the magnetostatic self-energy of the test body. The current generation of Eotvos experiments is shown to test the cited formalism at first order in the external gravitational potential to an accuracy of 4 by 10 to the -10th power for electrostatic energies and 6 by 10 to the -6th power for magnetostatic energies. It is concluded that although the formalism can be tested to an accuracy of 0.04 for electrostatic energies at second order in the gravitational potential, a significant second-order test for magnetostatic energies would require Eotvos experiments accurate to better than one part in 10 to the 15th power.

Haugan, M. P.↗

Recent studies of magnetic canopies

Two current studies are described which stem from Giovanelli's seminal studies of the spreading of chromospheric fields near active regions and active-region network. First, improved observational techniques are described for obtaining magnetograms in the Ca II 8542 A, Fe I 8688 A, and C I 9111 A lines which at least in principle allow for more accurate treatment of instrumental noise and allow better inference of field orientation. Second, a generalized response function is developed for calculating theoretical magnetograph signals from arbitrary line-of-sight variations of magnetic field, and initial applications to two-dimensional, potential-field models of network fields are described. Preliminary indications are that potential-field models can better explain the presence of low-lying, diffuse horizontal fields than can thin flux-tube models, but fail to predict a differential response between the different lines.

Jones, H. P.↗

Solar flux forecasting using mutual information with an optimal delay

Solar flux F(sub 10.7) directly affects the atmospheric density, thereby changing the lifetime and prediction of satellite orbits. For this reason, accurate forecasting of F(sub 10.7) is crucial for orbit determination of spacecraft. Our attempts to model and forecast F(sub 10.7) uncovered highly entangled dynamics. We concluded that the general lack of predictability in solar activity arises from its nonlinear nature. Nonlinear dynamics allow us to predict F(sub 10.7) more accurately than is possible using stochastic methods for time scales shorter than a characteristic horizon, and with about the same accuracy as using stochastic techniques when the forecasted data exceed this horizon. The forecast horizon is a function of two dynamical invariants: the attractor dimension and the Lyapunov exponent. In recent years, estimation of the attractor dimension reconstructed from a time series has become an important tool in data analysis. In calculating the invariants of the system, the first necessary step is the reconstruction of the attractor for the system from the time-delayed values of the time series. The choice of the time delay is critical for this reconstruction. For an infinite amount of noise-free data, the time delay can, in principle, be chosen almost arbitrarily. However, the quality of the phase portraits produced using the time-delay technique is determined by the value chosen for the delay time. Fraser and Swinney have shown that a good choice for this time delay is the one suggested by Shaw, which uses the first local minimum of the mutual information rather than the autocorrelation function to determine the time delay. This paper presents a refinement of this criterion and applies the refined technique to solar flux data to produce a forecast of the solar activity.

Ashrafi, S.↗

The effect of ionization on the infrared absorption spectra of PAHs: A preliminary report

The emission lines observed in many interstellar IR sources at 3.28, 6.2, 7.7, 8.7, and 11.3 microns are theorized to originate from polycyclic aromatic hydrocarbons (PAHs). These assignments are based on analyses of lab IR spectra of neutral PAHs. However, it is likely that in the interstellar medium that PAHs are ionized, i.e., are positively charged. Besides, as pointed out by Allamandola et al., although the IR emission band spectrum resembles what one might expect from a mixture of PAHs, it does not match in details such as frequency, band profile, or relative intensities predicted from the absorption spectra of any known PAH molecule. One source of more information to test the PAH theory is ab initio molecular orbital theory. It can be used to compute, from first principles, the geometries, vibrational frequencies, and vibrational intensities for model PAH compounds which are difficult to study in the lab. The Gaussian 86 computer program was used to determine the effect of ionization on the infrared absorption spectra of several small PAHs: naphthalene and anthracene. A preliminary report is presented of the results of these calculations.

Defrees, Doug J.↗

Can Vertical Profiles of Tropospheric Methane on Titan Be Derived from Radio-Occultation Soundings?

The intensity of the received signal at Earth in the radio occultations of Titan is attenuated both by refractive defocusing and pressure-induced absorption from N2-N2 and CH4-N2 pairs. Because the absorption strength is different for the two sets of pairs, matching the retrieved absorptivity profile can in principle yield the vertical variation in gaseous methane in the troposphere. There are two factors that make this difficult. The first is the propagation of noise in the phase and amplitude of the received signal in the absorption retrieval. The phase data is first inverted to retrieve vertical profiles of refractivity, from which the refractive defocusing is calculated. This is then subtracted from the observed. intensity attenuation of the received signal to generate a profile of atmospheric absorption. The second problem is the uncertainty in the pressure-induced absorption coefficients. Laboratory data at radio wavelengths is only available near room temperature (see, e.g., [1] for N2-N2), and the extrapolation to the low temperatures in Titan's troposphere is not well established. Ab initio calculations by Borysow et al. [2, 3] provide absorption coefficients at low temperatures and long wavelengths, but their accuracy has come into question. We present examples from Cassini radio occultations of Titan to illustrate the difficulties. For methane mole fractions in the lower troposphere comparable to that inferred from the Huygens probe (approximately 0.05), it will be difficult to separate the contributions of N2-N2 collisions from those of N2-CH4, collisions to the retrieved absorption. However, higher concentrations of CH4 and/or a higher signal-to-noise ratio from a future uplink experiment could result in a successful separation of the two components. However, key to this are highly accurate estimates of the absorption from a combination of laboratory measurements at love temperatures and long wavelengths, and possibly improved theoretical calculations.

Flasar, F. Michael↗

Mirroring within the Fokker-Planck formulation of cosmic ray pitch angle scattering in homogeneous magnetic turbulence

The Fokker-Planck coefficient for pitch angle scattering, appropriate for cosmic rays in homogeneous, stationary, magnetic turbulence, is computed from first principles. No assumptions are made concerning any special statistical symmetries the random field may have. This result can be used to compute the parallel diffusion coefficient for high energy cosmic rays moving in strong turbulence, or low energy cosmic rays moving in weak turbulence. Becuase of the generality of the magnetic turbulence which is allowed in this calculation, special interplanetary magnetic field features such as discontinuities, or particular wave modes, can be included rigorously. The reduction of this results to previously available expressions for the pitch angle scattering coefficient in random field models with special symmetries is discussed. The general existance of a Dirac delta function in the pitch angle scattering coefficient is demonstrated. It is proved that this delta function is the Fokker-Planck prediction for pitch angle scattering due to mirroring in the magnetic field.

Goldstein, M. L.↗

Application of theoretical models to active and passive remote sensing of saline ice

The random medium model is used to interpret the polarimetric active and passive measurements of saline ice. The ice layer is described as a host ice medium embedded with randomly distributed inhomogeneities, and the underlying sea water is considered as a homogeneous half-space. The scatterers in the ice layer are modeled with an ellipsoidal correlation function. The orientation of the scatterers is vertically aligned and azimuthally random. The strong permittivity fluctuation theory is employed to calculate the effective permittivity and the distorted Born approximation is used to obtain the polarimetric scattering coefficients. We also calculate the thermal emissions based on the reciprocity and energy conservation principles. The effects of the random roughness at the air-ice, and ice-water interfaces are accounted for by adding the surface scattering to the volume scattering return incoherently. The above theoretical model, which has been successfully applied to analyze the radar backscatter data of the first-year sea ice near Point Barrow, AK, is used to interpret the measurements performed in the CRRELEX program.

Han, H. C.↗

Towards an MHD Theory for the Standoff Distance of Earth's Bow Shock

A magnetohydrodynamic (MHD) theory is developed for the standoff distance a(s) of the bow shock and the thickness Delta(ms) of the magnetosheath, using the empirical Spreiter et al. relation Delta(ms) = kX and the MHD density ratio X across the shock. The theory includes as special cases the well-known gasdynamic theory and associated phenomenological MHD-like models for Delta(ms) and As. In general, however, MHD effects produce major differences from previous models, especially at low Alfev (Ma) and Sonic (Ms) Mach numbers. The magnetic field orientation Ma, Ms and the ratio of specific heats gamma are all important variables of the theory. In contrast, the fast mode Mach number need play no direct role. Three principle conclusions are reached. First the gasdynamic and phenomenological models miss important dependences of field orientation and Ms generally provide poor approximations to the MHD results. Second, changes in field orientation and Ms are predicted to cause factor of approximately 4 changes in Delta(ms) at low Ma. These effects should be important when predicting the shock's location or calculating gramma from observations. Third, using Spreiter et al.'s value for k in the MHD theory leads to maxima a(s) values at low Ma and nominal Ms that are much smaller than observations and MHD simulations require. Resolving this problem requires either the modified Spreiter-like relation and larger k found in recent MHD simulations and/or a breakdown in the Spreiter-like relation at very low Ma.

Carins, Iver H.↗

SN 1991T - Gamma-Ray Observatory's first supernova?

Consideration is given to the explosion of the Type Ia supernova SN 1991T in the nearby galaxy NGC 4527 detected in gamma-ray lines by the recently launched GRO. The dominant gamma-line and continuum features of the new 'delayed detonation' model FDEFA1 are calculated and compared to those for standard deflagration models W7 and cdtg7. It is shown that there are many useful hard photon discriminants of the Type Ia explosion mechanism that can, in principle, be detected by the OSSE and COMPTEL instruments on the GRO. Either SN 1991T, if bright enough, or one of the several Type Ia supernovae expected to be within the GRO's range during its active life, may make it possible to settle the detonation/deflagration debate, verify the generic thermonuclear white dwarf model of Type Ia explosions, and calibrate the Type Ia B(max)/847 keV line flux ratio.

Burrows, Adam↗

Investigation of radar backscattering from second-year sea ice

The scattering properties of second-year ice were studied in an experiment at Mould Bay in April 1983. Radar backscattering measurements were made at frequencies of 5.2, 9.6, 13.6, and 16.6 GHz for vertical polarization, horizontal polarization and cross polarizations, with incidence angles ranging from 15 to 70 deg. The results indicate that the second-year ice scattering characteristics were different from first-year ice and also different from multiyear ice. The fading properties of radar signals were studied and compared with experimental data. The influence of snow cover on sea ice can be evaluated by accounting for the increase in the number of independent samples from snow volume with respect to that for bare ice surface. A technique for calculating the snow depth was established by this principle and a reasonable agreement has been observed. It appears that this is a usable way to measure depth in snow or other snow-like media using radar.

Lei, Guang-Tsai↗

Analysis of the feasibility of an experiment to measure carbon monoxide in the atmosphere

The feasibility of measuring atmospheric carbon monoxide from a remote platform using the correlation interferometry technique was considered. It has been determined that CO data can be obtained with an accuracy of 10 percent using this technique on the first overtone band of CO at 2.3 mu. That band has been found to be much more suitable than the stronger fundamental band at 4.6 mu. Calculations for both wavelengths are presented which illustrate the effects of atmospheric temperature profiles, inversion layers, ground temperature and emissivity, CO profile, reflectivity, and atmospheric pressure. The applicable radiative transfer theory on which these calculations are based is described together with the principles of the technique.

Bortner, M. H.↗

Stirling engine design manual

This manual is intended to serve both as an introduction to Stirling engine analysis methods and as a key to the open literature on Stirling engines. Over 800 references are listed and these are cross referenced by date of publication, author and subject. Engine analysis is treated starting from elementary principles and working through cycles analysis. Analysis methodologies are classified as first, second or third order depending upon degree of complexity and probable application; first order for preliminary engine studies, second order for performance prediction and engine optimization, and third order for detailed hardware evaluation and engine research. A few comparisons between theory and experiment are made. A second order design procedure is documented step by step with calculation sheets and a worked out example to follow. Current high power engines are briefly described and a directory of companies and individuals who are active in Stirling engine development is included. Much remains to be done. Some of the more complicated and potentially very useful design procedures are now only referred to. Future support will enable a more thorough job of comparing all available design procedures against experimental data which should soon be available.

Martini, W. R.↗

Architecture Synthesis and Reduced-Cost Architectures for Human Exploration Missions

Development of architectures for human exploration missions has been pursued in the international aerospace community for a long time. This paper attempts a different approach and way of looking at architectures. Most of the emphasis is on lunar architectures with a brief look at Mars. The first step is to set forth overarching gods in order to understand origins of requirements. Then, principles and guidelines are developed for architecture formulation. It is argued that safety and cost are the primary factors. Alternative mission profiles are examined for adherence to the principles, and specific architectures formulated according to the guidelines. The guidelines themselves indicate preferred evolution paths from lunar to Mars architectures. Results of example calculations are given to illustrate the process, and an evolution path is recommended. Safety and cost criteria tend to conflict, but it is shown that cost-efficient architectures can be enhanced for good safety ratings at modest cost.

Woodcock, Gordon↗

Uncertainty Quantification of Classical Theories of Dendritic Growth Kinetics Applied to Nickel-Based Alloys

The solidification velocity in a model nickel-alloy single crystal during laser spot melting was recently characterized using synchrotron X-ray imaging. The measured solidification velocity was found to exceed the absolute stability threshold predicted by the Kurz-Giovanola-Trivedi (KGT) model. The discrepancies between the model and experiments motivate the further assessment of accurate material properties. This work quantifies the impact material property uncertainty has on model predictions of the absolute stability threshold velocities. Properties from the literature are reviewed and compared to those calculated using computational thermodynamics to provide uncertainty estimates on input properties to the KGT model. Global sensitivity analysis is used to quantify the influence of each uncertain input property on the predicted threshold velocity. This work supports the understanding of the nickel-alloy solidification during powder bed fusion additive manufacturing and identifies the solidification material properties that are the most important to assess from first-principles computations and experiments.

Computational thermodynamics↗

Structure of supersonic jet flow and its radiated sound

The present paper explores the use of large-eddy simulations as a tool for predicting noise from first principles. A high-order numerical scheme is used to perform large-eddy simulations of a supersonic jet flow with emphasis on capturing the time-dependent flow structure representating the sound source. The wavelike nature of this structure under random inflow disturbances is demonstrated. This wavelike structure is then enhanced by taking the inflow disturbances to be purely harmonic. Application of Lighthill's theory to calculate the far-field noise, with the sound source obtained from the calculated time-dependent near field, is demonstrated. Alternative approaches to coupling the near-field sound source to the far-field sound are discussed.

Mankbadi, Reda R.↗

A Comprehensive Reliability Methodology for Assessing Risk of Reusing Failed Hardware Without Corrective Actions with and Without Redundancy

This paper deals with the development of a reliability methodology to assess the consequences of using hardware, without failure analysis or corrective action, that has previously demonstrated that it did not perform per specification. The subject of this paper arose from the need to provide a detailed probabilistic analysis to calculate the change in probability of failures with respect to the base or non-failed hardware. The methodology used for the analysis is primarily based on principles of Monte Carlo simulation. The random variables in the analysis are: Maximum Time of Operation (MTO) and operation Time of each Unit (OTU) The failure of a unit is considered to happen if (OTU) is less than MTO for the Normal Operational Period (NOP) in which this unit is used. NOP as a whole uses a total of 4 units. Two cases are considered. in the first specialized scenario, the failure of any operation or system failure is considered to happen if any of the units used during the NOP fail. in the second specialized scenario, the failure of any operation or system failure is considered to happen only if any two of the units used during the MOP fail together. The probability of failure of the units and the system as a whole is determined for 3 kinds of systems - Perfect System, Imperfect System 1 and Imperfect System 2. in a Perfect System, the operation time of the failed unit is the same as that of the MTO. In an Imperfect System 1, the operation time of the failed unit is assumed as 1 percent of the MTO. In an Imperfect System 2, the operation time of the failed unit is assumed as zero. in addition, simulated operation time of failed units is assumed as 10 percent of the corresponding units before zero value. Monte Carlo simulation analysis is used for this study. Necessary software has been developed as part of this study to perform the reliability calculations. The results of the analysis showed that the predicted change in failure probability (P(sub F)) for the previously failed units is as high as 49 percent above the baseline (perfect system) for the worst case. The predicted change in system P(sub F) for the previously failed units is as high as 36% for single unit failure without any redundancy. For redundant systems, with dual unit failure, the predicted change in P(sub F) for the previously failed units is as high as 16%. These results will help management to make decisions regarding the consequences of using previously failed units without adequate failure analysis or corrective action.

Putcha, Chandra S.↗

Perturbation theory for a resistivity shunted Josephson element

We present a systematic perturbation theory, extendable in principle to all orders of magnitude, for the solution of the equations of motion of an ideal Josephson element shunted by a resistance and driven by a dc current source and a small time-dependent source. We present second-order results for the case in which the time dependence is that of a single sinusoid, and these results are compared with other numerical and analytical calculations. Near, but not on, the first constant voltage step where the perturbation theory appears divergent, the current-voltage characteristic is calculated by means of a nonperturbative adiabatic procedure. The impedance and responsivity agree with earlier results.

Thompson, E. D.↗