Computational methods for calculation of rates of evaporation in the Ta-Hf-C system
Computational methods for calculation of evaporation rates in tantalum-hafnium-carbon system
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.
Computational methods for calculation of evaporation rates in tantalum-hafnium-carbon system
Matrix method for calculating source function, mean intensity, and flux in model stellar atmosphere
A simple method for calculating the characteristics of the Dutch roll motion of an airplane is obtained by arranging the lateral equations of motion in such form and order that an iterative process is quickly convergent.
This report describes a method that can calculate transient aerodynamic heating and transient surface temperatures at supersonic and hypersonic speeds. This method can rapidly calculate temperature and heating rate time-histories for complete flight trajectories. Semi-empirical theories are used to calculate laminar and turbulent heat transfer coefficients and a procedure for estimating boundary-layer transition is included. Results from this method are compared with flight data from the X-15 research vehicle, YF-12 airplane, and the Space Shuttle Orbiter. These comparisons show that the calculated values are in good agreement with the measured flight data.
Two approaches to accelerating the method of complete linearization for calculating NLTE model stellar atmospheres are suggested. The first one, the so-called Kantorovich variant of the Newton-Raphson method, consists of keeping the Jacobi matrix of the system fixed, which allows us to calculate the costly matrix inversions only a few times and then keep them fixed during the subsequent computations. The second method is an application of the Ng acceleration. Both methods are extremely easy to implement with any model atmosphere code based on complete linearization. It is demonstrated that both methods, and especially their combination, yield a rapidly and globally convergent algorithm, which takes 2 to 5 times less computer time, depending on the model at hand and the required accuracy, than the ordinary complete linearization. Generally, the time gain is more significant for more complicated models. The methods were tested for a broad range of atmospheric parameters, and in all cases they exhibited similar behavior. Ng acceleration applied on the Kantorovich variant thus offers a significant improvement of the standard complete-linearization method, and may now be used for calculating relatively involved NLTE model stellar atmospheres.
Analytical and numerical methods of calculation of the energy and three dimensional EPS characteristics are reported. The angular and lateral functions of electrons in EPS have been obtained by the Landau and small angle approximations A and B and compared with earlier data. A numerical method of solution of cascade equations for the EPS distribution function moments has been constructed. Considering the equilibrium rms angle as an example, errors appearing when approximating the elementary process cross sections by their asymptotic expressions are analyzed.
A method is described for calculation of the radiation pattern of large aperture antennas. A piece-wise linear approximation of the aperture field using overlapping pyramidal basis functions allows the radiation pattern of an aperture antenna to be calculated as though it were a two-dimensional array. The calculation of radiation pattern data versus theta and phi, suitable for 3-D or contour plot algorithms, is achieved by locating the array in the yz-plane and performing a summation over the aperture field data sampled on a square grid. A FORTRAN subroutine is provided for performing radiation pattern calculations. Numerical results are included to demonstrate the accuracy and convergence of the method. These numerical results indicate that typical accuracies of + or - 0.1 dB for Directivity, + or - dB for the 1st Sidelobe Level, and + - 2dB for the 2nd Sidelobe Level can be obtained with an aperture grid of 45x45 points and requires approximately 0.02 seconds CPU time per far-field data point on a VAX 11/750 with a floating point accelerator.
Direct numerical method for calculation of supersonic inviscid flow about axisymmetric blunt body at large angle of attack
The use of Murman and Krupp's (1971) method for calculating steady, inviscid, transonic flows with imbedded shock waves is discussed in application to two-dimensional problems for transonic flows past sharp and moderately blunt nosed geometries. A relaxation algorithm is used to solve the equations.
A method for calculating the induced velocities at the blades of an inclined propeller is presented. The calculations are based upon an assumed wake consisting of one skewed helical vortex per blade. The blade circulation is assumed constant with respect to both radial and azimuth positions. The restriction of uniform radial circulation can be removed by superposition; however, the assumption of constant azimuth-wise circulation restricts the analysis to small propeller inclinations. Numerical values have been obtained for one wake skew angle. These values indicate large effects due to wake vortex spacing and to the number of blades.
An approximate method of calculating the deformations of wings of uniform thickness having swept, M or W, Delta, and swept-tip plan forms is presented. The method employs an adjustment to the elementary beam theory to account for the effect of the triangular root portion of a swept wing on the deformation of the outboard section of the wing. To demonstrate the general applicability of the method, the modified elementary theory is applied to the more complex M or W, Delta, and swept-tip plan forms as well as to swept plan forms. For the purpose of calculating angles of attack, it is shown that the unmodified elementary beam theory applied to that part of the wing outboard of the root triangle produces satisfactory results. However, for calculating deflections it is necessary to include the effects of the root-triangle deformation.
The development of a potential-flow/boundary-layer method for calculating subsonic and transonic turbulent flow past airfoils with trailing-edge separation is reported. A moment-of-momentum integral boundary-layer method is used which employs the law-of-the-wall/law-of-the-wake velocity profile and a two-layer eddy-viscosity model and ignores the laminar sublayer. All integrals across the boundary layer are obtained in closed form. Separation is assumed to occur when the shearing-stress velocity vanishes. A closed-form solution is derived for separated-flow regions where the shearing stress is negligible. In the potential-flow method, the exact form of the airfoil boundary condition is used, but it is applied at the chord line rather than the airfoil surface. This allows the accurate computation of flow about airfoils at large angles of attack but permits the use of body-oriented Cartesian computational grids. The governing equation for the perturbation velocity potential contains several terms in addition to the classical small-disturbance terms.
Approximative method for calculating laminar and turbulent heat transfer over bodies of revolution at angle of attack
Time-dependant method for calculating supersonic blunt body flow fields with sharp corners and embedded shock waves
A method is presented for calculating the aerodynamic heating and shear stresses at the wall for tangent ogive noses that are slender enough to maintain an attached nose shock through that portion of flight during which heat transfer from the boundary layer to the wall is significant. The lower entropy of the attached nose shock combined with the inclusion of the streamwise pressure gradient yields a reasonable estimate of the actual flow conditions. Both laminar and turbulent boundary layers are examined and an approximation of the effects of (up to) moderate angles-of-attack is included in the analysis. The analytical method has been programmed in FORTRAN 4 for an IBM 360/91 computer.
A method is presented for calculating the aerodynamic heating and shear stresses at the wall for tangent ogive noses that are slender enough to maintain an attached nose shock through that portion of flight during which heat transfer from the boundary layer to the wall is significant. The lower entropy of the attached nose shock combined with the inclusion of the streamwise pressure gradient yields a reasonable estimate of the actual flow conditions. Both laminar and turbulent boundary layers are examined and an approximation of the effects of (up to) moderate angles-of-attack is included in the analysis. The analytical method has been programmed in FORTRAN IV for an IBM 360/91 computer.
The Brinkley, Huff, and White methods for chemical-equilibrium calculations were modified and extended in order to permit an analytical comparison. The extended forms of these methods permit condensed species as reaction products, include temperature as a variable in the iteration, and permit arbitrary estimates for the variables. It is analytically shown that the three extended methods can be placed in a form that is independent of components. In this form the Brinkley iteration is identical computationally to the White method, while the modified Huff method differs only'slightly from these two. The convergence rates of the modified Brinkley and White methods are identical; and, further, all three methods are guaranteed to converge and will ultimately converge quadratically. It is concluded that no one of the three methods offers any significant computational advantages over the other two.
A discrete singularity method has been developed for calculating the potential flow around two-dimensional airfoils. The objective was to calculate velocities at any arbitrary point in the flow field, including points that approach the airfoil surface. That objective was achieved and is demonstrated here on a Joukowski airfoil. The method used combined vortices and sources ''submerged'' a small distance below the airfoil surface and incorporated a near-field subvortex technique developed earlier. When a velocity calculation point approached the airfoil surface, the number of discrete singularities effectively increased (but only locally) to keep the point just outside the error region of the submerged singularity discretization. The method could be extended to three dimensions, and should improve nonlinear methods, which calculate interference effects between multiple wings, and which include the effects of force-free trailing vortex sheets. The capability demonstrated here would extend the scope of such calculations to allow the close approach of wings and vortex sheets (or vortices).