SEARCH · Search NASA
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.
Indirect method considers effect of nonlinearity and obtains approximate probability densities of certain characteristics, from these densities various statistical properties are calculated. Direct method determines response moments directly. Neither method is restricted by system motion with regard to whether or not it is stationary.
The equations of motion necessary to compute the transient response of multi-mass flexible rotors are formulated to include unbalance, rotor acceleration, and flexible damped nonlinear bearing stations. A method of calculating the unbalance response of flexible rotors from a modified Myklestad-Prohl technique is discussed in connection with the method of solution for the transient response. Several special cases of simplified rotor-bearing systems are presented and analyzed for steady-state response, stability, and transient behavior. These simplified rotor models produce extensive design information necessary to insure stable performance to elastic mounted rotor-bearing systems under varying levels and forms of excitation. The nonlinear journal bearing force expressions derived from the short bearing approximation are utilized in the study of the stability and transient response of the floating bush squeeze damper support system. Both rigid and flexible rotor models are studied, and results indicate that the stability of flexible rotors supported by journal bearings can be greatly improved by the use of squeeze damper supports. Results from linearized stability studies of flexible rotors indicate that a tuned support system can greatly improve the performance of the units from the standpoint of unbalanced response and impact loading. Extensive stability and design charts may be readily produced for given rotor specifications by the computer codes presented in this analysis.
A fourth order box method for calculating high accuracy numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations is presented. The method is the natural extension of the second order Keller Box scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary layer equations. Numerical results for high accuracy test cases show the method to be significantly faster than other higher order and second order methods.
An approximate method for calculating the flow properties of gas-particle mixture flowing over blades in a cascade is studied. Using an analytical method, the solid particle trajectory and location of collisions between the solid particles and the blade surfaces are determined. In addition an experimental investigation of the trajectories and velocities of solid particles suspended in a fluid passing through an axial flow compressor cascade was performed. The cascade blades were made of 2024 aluminum alloy and the solid particles used were quartz sand with average diameter of 165 microns. The following parameters were investigated: blade pressure distribution, and total pressure loss coefficient, both depending on the degree of erosion. In addition, a theoretical estimation of the blade erosion is compared with the experimental data.
When the en route noise of a representative aircraft powered by an eight-blade SR-7 propeller was previously calculated, the noise level was cited as a possible concern associated with the acceptance of advanced turboprop aircraft. Some potential methods for reducing the en route noise were then investigated and are reported. Source noise reductions from increasing the blade number and from operating at higher rotative speed to reach a local minimum noise point were investigated. Greater atmospheric attenuations for higher blade passing frequencies were also indicated. Potential en route noise reductions from these methods were calculated as 9.5 dB (6.5 dB(A)) for a 10-blade redesigned propeller and 15.5 dB (11 dB(A)) for a 12-blade redesigned propeller.
A recent calculation of the antisymmetric stretch frequency for the rectangular structure of quartet O4(+) using the singles and doubles quadratic configuration interaction method with a perturbational estimate of connected triple excitations (QCISD(T)) method gave a value of 3710 cm(exp -1). This anomalous frequency is shown to be a consequence of symmetry breaking effects, which occur even though the QCISD(T) solution derived from a delocalized SCF reference function lies energetically well below the two localized (symmetry-broken) solutions at the equilibrium geometry. The symmetry breaking is almost eliminated at the CCSD level of theory, but the small remaining symmetry breaking effects are magnified at the CCSD(T) level of theory so that the antisymmetric stretch frequency is still significantly in error. The use of the Brueckner coupled cluster method, however, leads to a symmetrical solution which is free of symmetry breaking effects, with an antisymmetric stretch frequency of 1322 cm(exp -1), in good agreement with our earlier calculations using the complete active space self consistent field/complete active space state interaction (CASSCF/CASSI) method.
Herein we present numerical results of ground-state energies of 9 molecules in the well-established G2 molecule set given by the Gutzwiller conjugate gradient minimization (GCGM) method. The method, beyond the commonly used Gutzwiller approximation, was recently developed based on Gutzwiller variational wave functions. We find that compared to benchmark data given by full configuration interaction, GCGM total energies are reasonably well reproduced with the minimum basis set. To include the dynamical correlation beyond the minimal basis calculations, we adopt the local density approximation for the dynamical correlation energy $E_c$. By comparing the results with benchmark data given by experiments and large-basis configuration interaction, the GCGM total energies with $E_c$ are in general better reproduced, but discrepancies are still observed for some dimers.
Current efforts to reduce the sonic boom of a future High Speed Civil Transport (HSCT) by careful shaping have led to the need for more accurate predictions of the near-field flow conditions of the configuration. A fully three-dimensional Euler finite volume code is used to predict sonic boom pressure signatures for two low boom concepts - one designed to cruise at Mach 2 and the other at Mach 3. Calculations were carried out using a grid topology that has been modified to reduce the inaccuracies caused by grid spreading often suffered with CFD methods when calculations several body lengths downstream become necessary. Comparisons of CFD results and experimental wind tunnel signatures are shown. Ground signatures are predicted by extrapolating the pressures predicted by the Euler code with an extrapolation method based on the Whitham theory.
Jurisdictions are increasingly adopting compensation structures for distributed PV and PV-battery systems that price exported energy lower than energy consumed onsite (net billing rates). Standard methods for calculating the bill savings from PV and PV-battery systems were developed for net metering structures, and applying these same methods to net billing structures (such as using Typical Meteorological Year weather with actual year load) introduces bias errors that underestimate PV-battery system bill savings by between 1.5\% and 9\%, depending on the utility rate. We assess the magnitude of these errors and compare them to other sources of uncertainty when estimating the bill savings from PV and PV-battery systems under more complex utility rates.
We present an evaluation of a new nongeometric technique for calculating right and left ventricular volumes. This method calculates ventricular chamber volumes from multiple cross-sectional echocardiographic views taken from a single point as the echo beam is tilted progressively through the ventricle. Right and left ventricular volumes are calculated from both the approximate short axis and approximate apical position on 20 in vitro human hearts and compared with the actual chamber volumes. The results for both ventricles from both positions are excellent. Correlation coefficients are > 0.95 for all positions; the standard errors are in the range of 5 to 7 mL and the slopes and intercepts for the regression lines are not significantly different from 1 and 0, respectively (except for the left ventricular short-axis intercept). For all positions, approximately 6 to 8 views are needed for peak accuracy (7.5 degrees to 10 degrees separation). This approach offers several advantages. No geometric assumptions about ventricular shape are made. All images are acquired from a single point (or window), and the digitized points can be used to make a three-dimensional reconstruction of the ventricle. Also, during the calculations a volume distribution curve for the ventricle is produced. The shape of this curve can be characteristic for certain situations (ie, right ventricle, short axis) and can be used to make new simple equations for calculating volume. We conclude that this is an accurate nongeometric method for determining both right and left ventricular volumes in vitro.
Accurate estimation of the partial atomic charges on metal centers is useful for understanding electronic and catalytic properties of materials. However, different methods of calculating these charges may give quite different results; this issue has been more widely studied for molecules than for solids. Here we study the charges on the metal centers of a test set of 18 solids containing transition metals by using density functional theory with several density functionals (PBE, PBE+U, TPSS, revTPSS, HLE17, revM06-L, B3LYP, B3LYP*, and other exchange-modified B3LYP functionals) and four charge models (Bader, Hirshfeld, CM5, and DDEC6). The test set contains 12 systems with nonmagnetic metal centers (eight metal oxides (MO 2 ), two metal sulfides (MS 2 ), and two metal selenides (MSe 2 )) and six ferromagnetic transition metal complexes. Our study shows that, among the four types of charges, Bader charges are the highest and Hirshfeld charges are the lowest for all the systems, regardless of the functional being used. The CM5 charges are bigger than DDEC6 charges for MX 2 with M = Ti or Mo and X = S or Se, but for the other 14 cases they are lower. We found that the most of the systems are sensitive to the Hubbard U parameters in PBE+U and to the percentage X of Hartree–Fock exchange in exchange-modified B3LYP; as we increase U or X, the charges on the metal atoms in MX 2 increase steadily. Testing different density functionals shows charges calculated with higher Hubbard U parameters in PBE+U are comparable to B3LYP (with 20% Hartree–Fock exchange). Among four meta-GGA functionals studied, the charges with HLE17 have the closest agreement with B3LYP. In conclusion, the variation of charges with choice of charge model is greater than the variation with choice of density functional.
Calculating flutter boundaries and estimating stresses in infinite spanwise array of panels using exact method of Zeydel
The method for calculating the pitch angle diffusion coefficient is investigated for cosmic rays in a static random magnetic field, using the resonance integral method. The pitch angle diffusion coefficient may be derived from the Vlasov equation via ensemble averaging.
A rapid and accurate method of calculating optimal impulsive transfers in the restricted problem of three bodies has been developed. The technique combines a multi-conic method of trajectory integration with primer vector theory and an accelerated gradient method of trajectory optimization. A unique feature is that the state transition matrix and the primer vector are found analytical without additional integrations or differentiations. The method has been applied to the determination of optimal two and three impulse transfers between the L2 libration point and circular orbits about both the earth and the moon.
The results of numerous investigators have shown the importance of labyrinth seal inlet swirl on the calculated dynamic stiffness of labyrinth seals. These results have not included any calculation of inlet leakage swirl as a function of geometry and sealing conditions of the given seal. This paper outlines a method of calculating the inlet swirl at a given seal by introducing a radial chamber to predict the gas swirl as it goes from the stage tip down to the seal location. For a centrifugal compressor, this amounts to including the flow path from the impeller discharge, down the back of the disk or front of the cover, then into the shaft seal or eye packing, respectively. The solution includes the friction factors of both the disk and stationary wall with account for mass flow rate and calculation of radial pressure gradients by a free vortex solution. The results of various configurations are discussed and comparisons made to other published results of disk swirl.
The effects of steady and impulse-type residual accelerations on dopant distributions during directional solidification in 2D and 3D 'generic' models of the Bridgman-Stockbarger technique are investigated using numerical methods. The calculations are based on the thermophysical properties of molten germanium doped with a low concentration of gallium. A Chebyshev collocation pseudospectral method is used for the solution of the governing momentum-, mass-, species-, and heat-transfer equations. Only convection caused by temperature gradients is considered. It is found that lateral nonuniformity in composition is very sensitive to the orientation of the steady component of the residual gravity vector and to the particular operating conditions under consideration. It is also found that laterally or radially averaged composition profiles are alone insufficient to describe the extent of residual convection in a spacecraft environment. The effects of impulse-type disturbances can be severe and can extend for times on the order of 1000 sec after the termination of the impulse.
Three methods of calculating correction factors developed to match steady lifting distributions of pressure, forces and moments on airfoil sections, and total forces and moments. Data for rectangular supercritical wing tested in NASA Langley Research Center Transonic Dynamics Tunnel used to determine correction factors to match distributions of surface pressure for mach numbers ranging from 0.266 to 0.8. These correction factors also applied to calculations of unsteady aerodynamic flow, and comparisons made with oscillatory experimental data for range of reduced frequencies at several mach numbers. Methods of combining analytical and experimental data should lead to more-accurate predictions of performance of airfoils. Also applied to such airfoils as those in gas-turbine engines, fans, and wind turbines.
A simple and easily applied method of calculating the dimensions and horsepower necessary for an airship to have any given performance is described and illustrated by examples. The method includes means for estimating the changes in performance or in size when modifications or new features are introduced into a design, involving increase or saving in weights, or changes in resistance or propulsive efficiency.