Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALCULATION METHODS”

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 1,081 records · Page 60

Ab initio multireference study of the BN molecule

The lowest 1Sigma(+) and 3Pi states of the BN molecule are studied using multireference configuration interaction (MRCI) and averaged coupled-pair functional (ACPF) methods and large atomic natural orbital (ANO) basis sets, as well as several coupled cluster methods. Our calculations strongly support a 3Pi ground state, but the a1Sigma(+) state lies only 381 +/- 100/cm higher. The a1Sigma(+) state wave function exhibits strong multireference character and, consequently, the predictions of the perturbationally-based single-reference CCSD(T) coupled cluster method are not as reliable in this case as the multireference results. The theoretical predictions for the spectroscopic constants of BN are in good agreement with experiment for the Chi3Pi state, but strongly suggest a misassignment of the fundamental vibrational frequency for the a1Sigma(+) state.

Martin, J. M. L.↗

Design Criteria for Low Profile Flange Calculations

An analytical method and a design procedure to develop flanged separable pipe connectors are discussed. A previously established algorithm is the basis for calculating low profile flanges. The characteristics and advantages of the low profile flange are analyzed. The use of aluminum, titanium, and plastics for flange materials is described. Mathematical models are developed to show the mechanical properties of various flange configurations. A computer program for determining the structural stability of the flanges is described.

Leimbach, K. R.↗

On the accuracy of shape sensitivity derivatives

Six methods for calculating shape sensitivity derivatives are compared for a two-material beam problem with a moving interface. It is found that as the finite-element mesh is refined, displacement sensitivity derivatives converge more slowly than the displacement themselves. Five of the methods agree fairly well, but the adjoint variational surface method provides substantially different results. However, the difference is found to reflect convergence from another direction to the same answer, rather than reduced accuracy.

Haftka, Raphael T.↗

surf3d: A 3-D finite-element program for the analysis of surface and corner cracks in solids subjected to mode-1 loadings

A computer program, surf3d, that uses the 3D finite-element method to calculate the stress-intensity factors for surface, corner, and embedded cracks in finite-thickness plates with and without circular holes, was developed. The cracks are assumed to be either elliptic or part eliptic in shape. The computer program uses eight-noded hexahedral elements to model the solid. The program uses a skyline storage and solver. The stress-intensity factors are evaluated using the force method, the crack-opening displacement method, and the 3-D virtual crack closure methods. In the manual the input to and the output of the surf3d program are described. This manual also demonstrates the use of the program and describes the calculation of the stress-intensity factors. Several examples with sample data files are included with the manual. To facilitate modeling of the user's crack configuration and loading, a companion program (a preprocessor program) that generates the data for the surf3d called gensurf was also developed. The gensurf program is a three dimensional mesh generator program that requires minimal input and that builds a complete data file for surf3d. The program surf3d is operational on Unix machines such as CRAY Y-MP, CRAY-2, and Convex C-220.

Raju, I. S.↗

A method of computing the transient temperature of thick walls from arbitrary variation of adiabatic-wall temperature and heat-transfer coefficient

A method of calculating the temperature of thick walls has been developed in which the time series and the response to a unit triangle variation of surface temperature concepts are used, together with essentially standard formulas for transient temperature and heat flow into thick walls. The method can be used without knowledge of the mathematical tools of its development. The method is particularly suitable for determining the wall temperature in one-dimensional thermal problems in aeronautics where there is a continuous variation of the heat-transfer coefficient and adiabatic-wall temperature. The method also offers a convenient means for solving the inverse problem of determining the heat-flow history when temperature history is known.

Hill, P R↗

Multipath analysis diffraction calculations

This report describes extensions of the Kirchhoff diffraction equation to higher edge terms and discusses their suitability to model diffraction multipath effects of a small satellite structure. When receiving signals, at a satellite, from the Global Positioning System (GPS), reflected signals from the satellite structure result in multipath errors in the determination of the satellite position. Multipath error can be caused by diffraction of the reflected signals and a method of calculating this diffraction is required when using a facet model of the satellite. Several aspects of the Kirchhoff equation are discussed and numerical examples, in the near and far fields, are shown. The vector form of the extended Kirchhoff equation, by adding the Larmor-Tedone and Kottler edge terms, is given as a mathematical model in an appendix. The Kirchhoff equation was investigated as being easily implemented and of good accuracy in the basic form, especially in phase determination. The basic Kirchhoff can be extended for higher accuracy if desired. A brief discussion of the method of moments and the geometric theory of diffraction is included, but seems to offer no clear advantage in implementation over the Kirchhoff for facet models.

Statham, Richard B.↗

Coherent reference generator phase stability

Approximate phase stability estimates for the coherent reference generator (CRG) unit in the DSN's Frequency and Timing Subsystem (FTS) are calculated. The method used involves estimating the phase noise introduced by CRG components based upon measurements made in the past on similar components in other parts of the FTS and obtaining the CRG phase noise from the component phase noises. Three estimates of phase stability are calculated: the fractional frequency change for a 5 C step in temperature, the phase noise spectral density, and the Allan standard deviation. It is found from these estimates that the CRG phase stability is better than that of the H-maser physics unit + receiver. Thus, the first step in improving FTS phase stability would be to make improvements in the H-maser physics unit + receiver. These results are corroborated by indirect clock stability estimates calculated from Doppler data.

Korwar, V. N.↗

Navier-Stokes simulations of tip vortices for fixed and rotating helicopter blades

Flowfield and tip vortex results are presented for a hovering rotor blade at subcritical and supercritical flow conditions for both nonlifting and lifting configurations. These results are calculated numerically by solving, in a time-accurate fashion, the unsteady thin layer Navier-Stokes equations written in rotor coordinates. The lifting calculations use an induced downwash correction, estimated from a free wake analysis, to the geometric angle of attack of the blades to account for the wake effect. Comparison of numerical results with the experimental data shows very good agreement for all cases considered. Alternate methods of calculating hovering rotor flowfield as steady state flowfield on isolated fixed-blade that have the same circulation distribution as that of rotor in hover are also explored.

Srinivasan, G. R.↗

Effective Uncertainty Quantification for Multi-Angle Polarimetric Aerosol Remote Sensing Over Ocean

Multi-angle polarimetric (MAP) measurements can enable detailed characterization of aerosol microphysical and optical properties and improve atmospheric correction in ocean color remote sensing. Advanced retrieval algorithms have been developed to obtain multiple geophysical parameters in the atmosphere–ocean system. Theoretical pixel-wise retrieval uncertainties based on error propagation have been used to quantify retrieval performance and determine the quality of data products. However, standard error propagation techniques in high-dimensional retrievals may not always represent true retrieval errors well due to issues such as local minima and the nonlinear dependence of the forward model on the retrieved parameters near the solution. In this work, we analyze these theoretical uncertainty estimates and validate them using a flexible Monte Carlo approach. The Fast Multi-Angular Polarimetric Ocean coLor (FastMAPOL) retrieval algorithm, based on efficient neural network forward models, is used to conduct the retrievals and uncertainty quantification on both synthetic HARP2 (Hyper-Angular Rainbow Polarimeter 2) and AirHARP (airborne version of HARP2) datasets. In addition, for practical application of the uncertainty evaluation technique in operational data processing, we use the automatic differentiation method to calculate derivatives analytically based on the neural network models. Both the speed and accuracy associated with uncertainty quantification for MAP retrievals are addressed in this study. Pixel-wise retrieval uncertainties are further evaluated for the real AirHARP field campaign data. The uncertainty quantification methods and results can be used to evaluate the quality of data products, as well as guide MAP algorithm development for current and future satellite systems such as NASA’s Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission.

PACE↗

A lifting surface theory for thrust augmenting ejectors

The circulation theory of airfoil lift has been applied to calculate the performance of thrust augmenting ejectors. The ejector shroud is considered to be 'flying' in the secondary velocity field induced by the entrainment of the primary jet, so that the augmenting thrust is viewed as analogous to the lift on an airfoil. Vortex lattice methods are utilized to compute the thrust augmentation from the force on the flaps. The augmentation is shown to be a function of the length and shape of the flaps, as well as their position and orientation. Predictions of this new theory are compared with the results of classical methods of calculating the augmentation by integration of the stream thrust.

Bevilaqua, P. M.↗

The role of magnetic fields in the collapse of protostellar gas clouds

The paper presents the results of a numerical calculation of the collapse of an idealized protostellar gas cloud including the effects of a 'frozen-in' magnetic field. The 'traditional' picture of magnetic effects on gas clouds and recent observational and theoretical work on the subject are summarized. Attention is given to the method of calculation and the results are interpreted. It is found that the central magnetic field in the collapsing cloud model follows a rho to the 1/2 power relation, and the discussion implies that this is a general result which should hold true for some range of initial conditions around those chosen. In addition, it is found that the outer envelope of the cloud will be held up by tension in the field lines.

Scott, E. H.↗

Calculation of the binomial survivor function

A method is presented for calculating the binomial SF (cumulative binomial distribution), binfc(k;p,n), especially for a large n, beyond the range of existing tables, where conventional computer programs fail because of underflow and overflow, and Gaussian or Poisson approximations yield insufficient accuracy for the purpose at hand. This method is used to calculate and sum the individual binomial terms while using multiplication factors to avoid underflow; the factors are then divided out of the partial sum whenever it has the potential to overflow. A computer program uses this technique to calculate the binomial SF for arbitrary inputs of k, p, and n. Two other algorithms are presented to determine the value of p needed to yield a specified SF for given values of k and n and calculate the value where p = SF for a given k and n. Reliability applications of each algorithm/program are given, e.g., the value of p needed to achieve a stated k-out-of-n:G system reliability and the value of p for which k-out-of-n:G system reliability equals p.

Bowerman, Paul N.↗

A study of the interaction between a wake vortex and an encountering airplane

An analytical study of the mutual interaction between a wake vortex and an encountering airplane has been conducted. The interaction was calculated on the basis of conservation of angular momentum of the wake vortex and the encountering airplane in the roll axis. For the nominal conditions assumed in the study, the interaction reduced the roll response in the lust second by 20 percent compared to the classical method of calculating the response. However, depending on the conditions assumed, the reduction varied from 0 to 30 percent which could be highly significant in determining separation standards for airplanes. Controlled experiments are needed to determine the exact amount of interaction.

Stewart, E. C.↗

The mathematical-physical problem of P prime(f) to N(h) inversion for anisotropic media

The requirements for a numerical inversion method to calculate electron density profiles from ionograms are discussed. A systematic investigation about the independence of the two magnetoionic components shows that the extent to which the range of ambiguity can be reduced depends mainly on the magnetic latitude. Error estimates and first-order corrections for less elaborate methods are included.

Paul, A. K.↗

Production of H-3 at large momentum in alpha-(C-12) collisions at 2A GeV

Calculations of the longitudinal and transverse momentum distributions for H-3 production at large momentum in alpha-(C-12) collisions near 2A GeV are compared to experiment. Triton exchange and final-state interactions are shown to represent large corrections to the impulse approximation for proton knockout. A method for calculating interference effects for inelastic fragmentation is discussed. Good agreement with experiment using a phenomenological overlap function for (H-3)-p is found successful in describing H-3 production in pion-induced reactions. Comparisons to momentum distributions obtained through (p, 2p) and (e, ep) reactions on He-4 are made.

Cucinotta, Francis A.↗

Enhancements to Linear Stability-Based, CFD-integrated Transition Prediction for High-Speed Flows

Combining linear stability calculations with computational fluid dynamics (CFD) simulations has great potential for the automated modeling of high-speed flows, especially when adequate information about the configuration and the disturbance environment is available. However, a significant impediment to the applicability of this technique is the lack of an efficient method to calculate the crucial amplification ratio corresponding to the onset of transition in hypersonic flows. This ratio, also known as the "transition N-factor," is dependent upon the freestream disturbance environment as well as the surface properties of the test article. In response to the need for an engineering solution to predict the transition N-factor within conventional hypersonic wind tunnels, this paper presents a data-driven correlation that expands the existing correlations from straight circular cones with a narrow range of half angles to a broader array of axisymmetric configurations. Furthermore, when tested against a chosen dataset that was not used in its calibration, the suggested correlation shows good predictive accuracy with an RMS error of only 6.9%. Although similar accuracy may also be achieved via existing correlations based on similar datasets, predictions based on the proposed correlation have the advantage of not requiring an extensive amount of configuration-specific data. Practical applications often have access to the input parameters needed for this correlation, such as the freestream disturbance intensity, Mach number, and body-based slenderness Reynolds number. Additionally, this correlation outperforms the traditional assumption of a constant N-factor, particularly for configurations with blunted nose geometries. The development of this correlation is grounded in an extensive dataset encompassing conical models with body half-angles varying between 5 degrees and 16 degrees, Mach numbers ranging from 5 to 14, and nosetip-based Reynolds numbers approaching the transition reversal limit for blunt-nosed cones.

CFD↗

Enhancements to Linear Stability-Based, CFD-integrated Transition Prediction for High-Speed Flows

Combining linear stability calculations with computational fluid dynamics (CFD) simulations has great potential for the automated modeling of high-speed flows, especially when adequate information about the configuration and the disturbance environment is available. However, a significant impediment to the applicability of this technique is the lack of an efficient method to calculate the crucial amplification ratio corresponding to the onset of transition in hypersonic flows. This ratio, also known as the "transition N-factor," is dependent upon the freestream disturbance environment as well as the surface properties of the test article. In response to the need for an engineering solution to predict the transition N-factor within conventional hypersonic wind tunnels, this paper presents a data-driven correlation that expands the existing correlations from straight circular cones with a narrow range of half angles to a broader array of axisymmetric configurations. Furthermore, when tested against a chosen dataset that was not used in its calibration, the suggested correlation shows good predictive accuracy with an RMS error of only 6.9%. Although similar accuracy may also be achieved via existing correlations based on similar datasets, predictions based on the proposed correlation have the advantage of not requiring an extensive amount of configuration-specific data. Practical applications often have access to the input parameters needed for this correlation, such as the freestream disturbance intensity, Mach number, and body-based slenderness Reynolds number. Additionally, this correlation outperforms the traditional assumption of a constant N-factor, particularly for configurations with blunted nose geometries. The development of this correlation is grounded in an extensive dataset encompassing conical models with body half-angles varying between 5 degrees and 16 degrees, Mach numbers ranging from 5 to 14, and nosetip-based Reynolds numbers approaching the transition reversal limit for blunt-nosed cones.

CFD↗

Flutter calculations in three degrees of freedom

The present paper is a continuation of the general study of flutter published in NACA reports nos. 496 and 685. The paper is mainly devoted to flutter in three degrees of freedom (bending, torsion, and aileron) for which a number of selected cases have been calculated and presented in graphical form. The results are analyzed and discussed with regard to the effects of structural damping, of fractional-span ailerons, and of mass-balancing. The analysis shows that more emphasis should be put on the effect of structural damping and less on mass-balancing. The conclusion is drawn that a definite minimum amount of structural damping, which is usually found to be present, is essential in the calculations for an adequate description of the flutter case. Theoretical flutter predictions are thus brought into closer agreement with the facts of experience. A brief discussion is included of a particular biplane that had experienced flutter at about 200 miles per hour. Some simplifications have been achieved in the method of calculation. (author)

Theodorsen, Theodore↗