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.
Point 7 of a list of works conducted by the USSR in accordance with the Joint Soviet-American Research Program on improving methods of heat probes from satellites is discussed. Numerous calculations of transmission functions, along with tables are enclosed.
The primary goal of this collaborative project, lead by Nicholas Gruber at UCLA, is to better quantify the air-sea fluxes of pre-industrial and anthropogenic carbon dioxide CO2 for the global ocean with a newly developed inverse method. The method involves using ocean general circulation models (OGCMs) to characterize the subsurface distribution of tracers injected for different regions of the ocean surface. By linearly combining a series of these basis (or Greens) functions, one should be able to recreate any arbitrary ocean tracer field. In our case, we are using the global database of ocean dissolved inorganic carbon (DIC) from the WOCE/JGOFS global CO2 survey from the 1990s. The observed DIC field has been partitioned into pre-industrial and anthropogenic components using the empirical C* technique. We then solve for an optimal combination of the regional basis functions that best match the two sets observations in a least-squares sense with singular value decomposition. This in turn provides solutions for the regional air-sea CO2 fluxes (pre-industrial and anthropogenic) that are independent of more traditional estimates derived from surface water disequilibria in the partial pressure of CO2 and wind speed gas exchange parameterizations. We are apply the method to a suite of nine different OGCMs to test the robustness of the derived surface fluxes fields to differences in ocean circulation. Together with a collegue Keith Lindsay at NCAR, S. Doney's role in this project was to create the appropriate steady-state and temporally evolving basis functions from the NCAR CSM-1 ocean circulation model. Because of a slow start across the project in defining a common set of experimental protocols, we requested and were given a one year, no cost extension. Both the steady-state and temporally evolving basis functions were delivered to the UCLA group for further analysis. Additional information is included in the original extended abstract.
Developing the technological response to realizing an efficient atmosphere revitalization system for future crewed spacecraft and space habitats requires identifying and describing functional trade spaces. Mission concepts and requirements dictate the necessary functions; however, the combination and sequence of those functions possess significant flexibility. Us-ing a closed loop environmental control and life support (ECLS) system architecture as a starting basis, a functional unit operations approach is developed to identify trade spaces. Generalized technological responses to each trade space are discussed. Key performance parameters that apply to functional areas are described.
Two meshless local Petrov-Galerkin (MLPG) methods based on two different trial functions but that use a simple linear test function were developed for beam and column problems. These methods used generalized moving least squares (GMLS) and radial basis (RB) interpolation functions as trial functions. These two methods were tested on various patch test problems. Both methods passed the patch tests successfully. Then the methods were applied to various beam vibration problems and problems involving Euler and Beck's columns. Both methods yielded accurate solutions for all problems studied. The simple linear test function offers considerable savings in computing efforts as the domain integrals involved in the weak form are avoided. The two methods based on this simple linear test function method produced accurate results for frequencies and buckling loads. Of the two methods studied, the method with radial basis trial functions is very attractive as the method is simple, accurate, and robust.
A fundamental input to models of light scattering from planetary regoliths is the mean phase function of the regolith particles. Using the known size distribution for typical lunar soils, the mean phase function and mean linear polarization for a regolith volume element of spherical particles of any composition were calculated from Mie theory. The two contour plots given here summarize the changes in the mean phase function and linear polarization with changes in the real part of the complex index of refraction, n - ik, for k equals 0.01, the visible wavelength 0.55 micrometers, and the particle size distribution of the typical mature lunar soil 72141. A second figure is a similar index-phase surface, except with k equals 0.1. The index-phase surfaces from this survey are a first order description of scattering by lunar-like regoliths of spherical particles of arbitrary composition. They form the basis of functions that span a large range of parameter-space.
The dynamics of vibrations in flexible structures can be conventiently modeled in terms of frequency response models. For structural control such models capture the distributed parameter dynamics of the elastic structural response as an irrational transfer function. For most flexible structures arising in aerospace applications the irrational transfer functions which arise are of a special class of pseudo-meromorphic functions which have only a finite number of right half place poles. Computational algorithms are demonstrated for design of multiloop control laws for such models based on optimal Wiener-Hopf control of the frequency responses. The algorithms employ a sampled-data representation of irrational transfer functions which is particularly attractive for numerical computation. One key algorithm for the solution of the optimal control problem is the spectral factorization of an irrational transfer function. The basis for the spectral factorization algorithm is highlighted together with associated computational issues arising in optimal regulator design. Options for implementation of wide band vibration control for flexible structures based on the sampled-data frequency response models is also highlighted. A simple flexible structure control example is considered to demonstrate the combined frequency response modeling and control algorithms.
Ab initio self-consistent-field plus configuration-interaction calculations are reported for the X(sup 1)Sigma(sup +), A(sup 1)Sigma(sup +), and B(sup 1)Pi states of LiH using a 22(sugma)12(pi)7(delta) function Slater basis set. The resulting dissociation energies, with the experimental values in parentheses, are D(sub e)(X(sup 1)Sigma(sup +)) = 19 972 (20288)/ cm, and D(sub e)(A(sup 1)Sigma(sup +)) = 9042 (8682)/ cm, and D(sup e)(B(sup 1)Pi) = 284 (288) /cm. This is the first ab initio treatment to quantitatively account for the binding in the B(sup 1)Pi state. Calculated dipole moments and electronic transition moments for the X(sup 1)Sigma(sup+)-A(sup 1)Sigma(sup +), X(sup 1)Sigma(sup +)-B((sup 1)Pi, and A(sup 1)Sigma(sup +)- B(sup 1)PI band systems are in excellent agreement with existing theoretical and experimental data. Radiative transition probabilities and lifetimes, including both the bound-bound and bound-free contributions, are computed for all vibrational levels of the A(sup 1)Sigma(sup+) and B(sup 1)Pi states. Comparison with previous results using experimentally based potentials provides insight into the sensitivity of the radiative lifetimes to the detailed nature of the uppermost region of the potentials. Our calculated lifetimes for the lower vibrational levels of the A(sup 1)Sigma(sup +) state are within the experimental uncertainty. Our calculated lifetimes for the three vibrational levels of the B(sup 1)Pi state are in excellent agreement with those of Zemke and Stwalley (values in parentheses), increasing with (upsilon)' from 11.3 (11.3) nsec at (upsilon)' = 0, to 17.0 (17.0) nsec at (upsilon)' = 1, and then to 23.5 (24.0) nsec at (upsilon)' = 2.
We report on our effort to develop algorithms to detect, size, and classify geological features in planetary and asteroidal data sets.
The requirements for very accurate ab initio quantum chemical prediction of dissociation energies are examined using a detailed investigation of the nitrogen molecule. Although agreement with experiment to within 1 kcal/mol is not achieved even with the most elaborate multireference CI (configuration interaction) wave functions and largest basis sets currently feasible, it is possible to obtain agreement to within about 2 kcal/mol, or 1 percent of the dissociation energy. At this level it is necessary to account for core-valence correlation effects and to include up to h-type functions in the basis. The effect of i-type functions, the use of different reference configuration spaces, and basis set superposition error were also investigated. After discussing these results, the remaining sources of error in our best calculations are examined.
The requirements for very accurate ab initio quantum chemical prediction of dissociation energies are examined using a detailed investigation of the nitrogen molecule. Although agreement with experiment to within 1 kcal/mol is not achieved even with the most elaborate multireference CI (configuration interaction) wave functions and largest basis sets currently feasible, it is possible to obtain agreement to within about 2 kcal/mol, or 1 percent of the dissociation energy. At this level it is necessary to account for core-valence correlation effects and to include up to h-type functions in the basis. The effect of i-type functions, the use of different reference configuration spaces, and basis set superposition error were also investigated. After discussing these results, the remaining sources of error in our best calculations are examined.
On the basis of 12-month moving averages (12-mma) of monthly mean sunspot number (R), sunspot cycle 24 had its minimum amplitude (Rm = 1.7) in December 2008. At 12 mo past minimum, R measured 8.3, and at 18 mo past minimum, it measured 16.4. Thus far, the maximum month-to-month rate of rise in 12-mma values of monthly mean sunspot number (AR(t) max) has been 1.7, having occurred at elapsed times past minimum amplitude (t) of 14 and 15 mo. Compared to other sunspot cycles of the modern era, cycle 24?s Rm and AR(t) max (as observed so far) are the smallest on record, suggesting that it likely will be a slow-rising, long-period sunspot cycle of below average maximum amplitude (RM). Supporting this view is the now observed relative strength of cycle 24?s geomagnetic minimum amplitude as measured using the 12-mma value of the aa-geomagnetic index (aam = 8.4), which also is the smallest on record, having occurred at t equals 8 and 9 mo. From the method of Ohl (the inferred preferential association between RM and aam), one predicts RM = 55 +/- 17 (the ?1 se prediction interval) for cycle 24. Furthermore, from the Waldmeier effect (the inferred preferential association between the ascent duration (ASC) and RM) one predicts an ASC longer than 48 mo for cycle 24; hence, maximum amplitude occurrence should be after December 2012. Application of the Hathaway-Wilson-Reichmann shape-fitting function, using an RM = 70 and ASC = 56 mo, is found to adequately fit the early sunspot number growth of cycle 24.
A universal algorithm is presented for calculating electrostatic fields in an infinite periodic structure consisting of electrodes of arbitrary shape which are located in mirror-symmetrical manner along the axis of electron-beam propagation. The method is based on the theory of R-functions, and the differential operators which are derived on the basis of the functions. Numerical results are presented and the accuracy of the results is examined.
The equilibrium bond distance in methane was optimized using coupled-pair functional and contracted CI wave functions, and a Gaussian basis that includes g-type functions on carbon and d-type functions on hydrogen. The resulting bond distance, when corrected for core-valence correlation effects, agrees with the experimental value of 2.052 a(0) to within the experimental uncertainty of 0.002 a(0). The main source of error in the best previous studies, which showed discrepancies with experiment of 0.007 a(0) is shown to be basis set incompleteness. In particular, it is important that the basis set be close to saturation, at least for the lower angular quantum numbers.
The equilibrium bond distance in methane has been optimized using coupled-pair functional and contracted CI wave functions, and a Gaussian basis that includes g-type functions on carbon and d-type functions on hydrogen. The resulting bond distance, when corrected for core-valance correlation effects, agrees with the experimental value of 2.052 a(0) to within the experimental uncertainty of 0.002 a(0). The main source of error in the best previous studies, which showed discrepancies with experiment of 0.007 a(0) is shown to be basis set incompleteness. In particular, it is important that the basis set be close to saturation, at least for the lower angular quantum numbers.
A recently proposed scheme for using natural orbitals from atomic configuration interaction wave functions as a basis set for linear combination of atomic orbitals (LCAO) calculations is extended for the calculation of molecular properties. For one-electron properties like multipole moments, which are determined largely by the outermost regions of the molecular wave function, it is necessary to increase the flexibility of the basis in these regions. This is most easily done by uncontracting the outermost Gaussian primitives, and/or by adding diffuse primitives. A similar approach can be employed for the calculation of polarizabilities. Properties which are not dominated by the long-range part of the wave function, such as spectroscopic constants or electric field gradients at the nucleus, can generally be treated satisfactorily with the original atomic natural orbital sets.
A recently proposed scheme for using natural orbitals from atomic configuration interaction (CI) wave functions as a basis set for linear combination of atomic orbitals (LCAO) calculations is extended for the calculation of molecular properties. For one-electron properties like multipole moments, which are determined largely by the outermost regions of the molecular wave function, it is necessary to increase the flexibility of the basis in these regions. This is most easily done by uncontracting the outmost Gaussian primitives, and/or by adding diffuse primitives. A similar approach can be employed for the calculation of polarizabilities. Properties which are not dominated by the long-range part of the wave function, such as spectroscopic constants or electric field gradients at the nucleus, can generally be treated satisfactorily with the original atomic natural orbital (ANO) sets.
Implementing specific image processing algorithms via very large scale integrated systems offers a potent solution to the problem of handling high data rates. Two algorithms stand out as being particularly critical -- geometric map transformation and filtering or correlation. These two functions form the basis for data calibration, registration and mosaicking. VLSI presents itself as an inexpensive ancillary function to be added to almost any general purpose computer and if the geometry and filter algorithms are implemented in VLSI, the processing rate bottleneck would be significantly relieved. A set of image processing functions that limit present systems to deal with future throughput needs, translates these functions to algorithms, implements via VLSI technology and interfaces the hardware to a general purpose digital computer is developed.
The geometric structure of Cr(CO)6 is optimized at the modified coupled-pair functional (MCPF), single and double excitation coupled-cluster (CCSD) and CCSD(T) levels of theory (including a perturbational estimate for connected triple excitations), and the force constants for the totally symmetric representation are determined. The geometry of Cr(CO)5 is partially optimized at the MCPF, CCSD and CCSD(T) levels of theory. Comparison with experimental data shows that the CCSD(T) method gives the best results for the structures and force constants, and that remaining errors are probably due to deficiencies in the one-particle basis sets used for CO. A detailed comparison of the properties of free CO is therefore given, at both the MCPF and CCSD/CCSD(T) levels of treatment, using a variety of basis sets. With very large one-particle basis sets, the SSCD(T) method gives excellent results for the bond distance, dipole moment and harmonic frequency of free CO. The total binding energies of Cr(CO)6 and Cr(CO)5 are also determined at the MCPF, CCSD and CCSD(T) levels of theory. The CCSD(T) method gives a much larger total binding energy than either the MCPF or CCSD methods. An analysis of the basis set superposition error (BSSE) at the MCPF level of treatment points out limitations in the one-particle basis used here and in a previous study. Calculations using larger basis sets reduced the BSSE, but the total binding energy of Cr(CO)6 is still significantly smaller than the experimental value, although the first CO bond dissociation energy of Cr(CO)6 is well described. An investigation of 3s3p correlation reveals only a small effect. The remaining discrepancy between the experimental and theoretical total binding energy of Cr(CO)6 is probably due to limitations in the one-particle basis, rather than limitations in the correlation treatment. In particular an additional d function and an f function on each C and O are needed to obtain quantitative results. This is underscored by the fact that even using a very large primitive se (1042 primitive functions contracted to 300 basis functions), the superposition error for the total binding energy of Cr(CO)6 is 22 kcal/mol at the MCPF level of treatment.