Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computation time”

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,315 records · Page 73

Noniterative cross-flow integration for the pressure-split analysis of subsonic mixing-layer problems

In utilizing pressure-split methodology for the two-dimensional analysis of curved wall jets and the subsonic regions of underexpanded free jets, Dash et al (1984) employed the cross flow procedure of Bradshaw et al (1979, 1984). In assessing this procedure, it was noted that the iterative steps required between the continuity and cross flow momentum equation solutions could be eliminated by combining these equations into a unified equation for the cross flow velocity, yielding a significant reduction in computer time and eliminating possible convergence problems which occur in iterative approaches. Attention is presently given to the application of this noniterative procedure to a two-dimensional curved wall jet problem.

Dash, S. M.↗

A three-dimensional Navier-Stokes/Euler code for blunt-body flow computations

The formulation computation method of an improved version of the three-dimensional Navier-Stokes/Euler computation algorithm of Li (1981) for the numerical simulation of blunt-body reentry flows are discussed, and results for five sample problems are presented graphically. The vector notation of the coordinate systems is defined; the governing equations are presented in full; the Jacobian matrices and damping terms of the factorization technique (based on the alternating-direction implicit procedure of Beam and Warming, 1978) are explained; and the capabilities, limitations and proper use of the code are summarized. Examples presented include the equilibrium-flow problem for Shuttle-orbiter reentry at Mach 22 and angle of attack 40.8 deg and the perfect-gas problem of an aerobraking orbital-transfer vehicle with an ellipsoidal/60-deg cone and a toroidal sonic shoulder. The advantages of the improved code in terms of accuracy and computation time are indicated.

Li, C. P.↗

Study of the kinetics and mechanism of the thermal nitridation of SiO2

X-ray photoelectron spectroscopy (XPS) has been used to study the nitridation time and temperature dependence of the nitrogen distribution in thermally nitrided SiO2 films. The XPS data show that the maximum nitrogen concentration near the (SiO(x)N(y)/Si interface is initially at the interface, but moves 20-25 A away from the interface with increasing nitridation time. Computer modeling of the kinetic processes involved is carried out and reveals a mechanism in which diffusing species, initially consisting primarily of nitrogen, react with the substrate, followed by formation of the oxygen-rich oxynitride due to reaction of the diffusing oxygen displaced by the slower nitridation of the SiO2. The data are consistent with this mechanism provided the influence of the interfacial strain on the nitridation and oxidation kinetics is explicitly accounted for.

Vasquez, R. P.↗

Computation of unsteady transonic aerodynamics with steady state fixed by truncation error injection

A novel technique is introduced for efficient computations of unsteady transonic aerodynamics. The steady flow corresponding to body shape is maintained by truncation error injection while the perturbed unsteady flows corresponding to unsteady body motions are being computed. This allows the use of different grids comparable to the characteristic length scales of the steady and unsteady flows and, hence, allows efficient computation of the unsteady perturbations. An example of typical unsteady computation of flow over a supercritical airfoil shows that substantial savings in computation time and storage without loss of solution accuracy can easily be achieved. This technique is easy to apply and requires very few changes to existing codes.

Fung, K.-Y.↗

Unsteady lifting-line theory as a singular-perturbation problem

Unsteady lifting-line theory is developed for a wing of large aspect ratio oscillating at low frequency in inviscid incompressible flow. The wing is assumed to have a rigid chord but a flexible span. Use of the method of matched asymptotic expansions reduces the problem from a singular integral equation to quadrature. The pressure field and airloads, for a prescribed wing shape and motion, are obtained in closed form as expansions in inverse aspect ratio. A rigorous definition of unsteady induced downwash is also obtained. Numerical calculations are presented for an elliptic wing in pitch and heave; compared with numerical lifting-surface theory, computation time is reduced significantly. The present work also identifies and resolves errors in the unsteady lifting line theory of James (1975), and points out a limitation in that of Van Holten (1975, 1976, 1977).

Ahmadi, A.↗

Computing Composition/Depth Profiles From X-Ray Diffraction

Diffraction-intensity bands deconvolved relatively quickly. TIBAC constructs composition/depth profiles from X-ray diffraction-intensity bands. Intensity band extremely sensitive to shape of composition/depth profile. TIBAC incorporates straightforward transformation of intensity band that retains accuracy of earlier simulation models, but is several orders of magnitude faster in total computational time. TIBAC written in FORTRAN 77 for batch execution.

Wiedemann, K. E.↗

Derivatives of the Arithmetic-Geometric Mean

Developed for theoretical studies of lightning requiring estimates of electric fields in clouds, the technique is expected to reduce computation time and improve accuracy.

Taton, F. B.↗

Strategy for reflector pattern calculation: Let the computer do the work

Using high frequency approximations, the secondary pattern of a reflector antenna can be calculated by numerically evaluating a radiation integral I(u,v). In recent years, tremendous effort has been expended to reducing I(u,v) to Fourier integrals. These reduction schemes are invariably reflector geometry dependent. Hence, different analyses/computer software development must be carried out for different reflector shapes/boundaries. it is pointed out, that, as the computer power improves, these reduction schemes are no longer necessary. Comparable accuracy and computation time can be achieved by evaluating I(u,v) by a brute force FFT described in this note. Furthermore, there is virtually no restriction on the reflector geometry by using the brute force FFT.

Lam, P. T.↗

Lubrication of nonconformal contacts

Minimum film thickness results for piezoviscous-rigid regime of lubrication are developed for a compressible Newtonian fluid with Roelands viscosity. The results provide a basis for the analysis and design of a wide range of machine elements operating in the piezoviscous-rigid regime of lubrication. A new numerical method of calculating elastic deformation in contact stresses is developed using a biquadratic polynomial to approximate the pressure distribution on the whole domain analyzed. The deformation of every node is expressed as a linear combination of the nodal pressures whose coefficients can be combined into an influence coefficient matrix. This approach has the advantages of improved numerical accuracy, less computing time and smaller storage size required for influence matrix. The ideal elastohydrodynamic lubrication is extended to real bearing systems in order to gain an understanding of failure mechanisms in machine elements. The improved elastic deformation calculation is successfully incorporated into the EHL numerical scheme. Using this revised numerical technique and the flow factor model developed by Patir and Cheng (1978) the surface roughness effects on the elastohydrodynamic lubrication of point contact is considered. Conditions typical of an EHL contact in the piezoviscous-elastic regime entrained in pure rolling are investigated. Results are compared with the smooth surface solutions. Experiments are conducted to study the transient EHL effects in instrument ball bearings.

Jeng, Y. R.↗

Equivalent plate analysis of aircraft wing box structures with general planform geometry

A new equilvalent plate analysis formulation is described which is capable of modeling aircraft wing structures with a general planform such as cranked wing boxes. Multiple trapezoidal segments are used to represent such planforms. A Ritz solution technique is used in conjunction with global displacement functions which encompass all the segments. This Ritz solution procedure is implemented efficiently into a computer program so that it can be used by rigorous optimization algorithms for application in early preliminary design. A direct method to interface this structural analysis procedure with aerodynamic programs for use in aeroelastic calculations is described. This equivalent plate analysis procedure is used to calculate the static deflections and stresses and vibration frequencies and modes of an example wing configuration. The numerical results are compared with results from a finite element model of the same configuration to illustrate typical levels of accuracy and computation times resulting from use of this procedure.

Giles, G. L.↗

An E-plane analysis of aperture-matched horn antennas using the moment method and the uniform geometrical theory of diffraction

The moment method and the uniform geometrical theory of diffraction are utilized to obtain two separate solutions for the E-plane field pattern of an aperture-matched horn antenna. This particular horn antenna consists of a standard pyramidal horn with the following modifications: a rolled edge section attached to the aperture edges and a curved throat section. The resulting geometry provides significantly better performance in terms of the pattern, impedance, and frequency characteristics than normally obtainable. The moment method is used to calculate the E-plane pattern and BSWR of the antenna. However, at higher frequencies, large amounts of computation time are required. The uniform geometrical theory of diffraction provides a quick and efficient high frequency solution for the E-plane field pattern. In fact, the uniform geometrical theory of diffraction may be used to initially design the antenna; then, the moment method may be applied to fine tune the design. This procedure has been successfully applied to a compact range feed design.

Heedy, D. J.↗

Structural analysis of turbine blades using unified constitutive models

The utility of advanced constitutive models and structural analysis methods in predicting the cyclic life of an air-cooled turbine blade is assessed. Five structural analysis methods were exercised in calculating the cyclic stress-strain response at the airfoil critical location. The methods studied were a cyclic elastic finite-element analysis, nonlinear finite-element analyses based on classical inelastic models and the unified models of Bodner and Walker, and a simplified inelastic procedure. These analyses were compared in terms of computing times and of predicted crack initiation lives using the Strainrange Partitioning method.

Kaufman, A.↗

Improved 2-dimensional-kinetics computer program

The analytical capability of the existing Two-Dimensional Kinetics (TDK)/Boundary Layer Module (BLM) computer program for performance of Orbit Transfer Vehicle (OTV) thrust chambers was examined. Areas which need further examination and improvement are the thick boundary layer inviscid core flow interaction and the related thrust loss calculation. To warrant highly accurate results the provision of the best available program input data is mandatory, as well as, the use of sophisticated modeling techniques to reduce computation time with advanced error control criteria. The simulation of wall shocks, Mach-shocks, and shocks induced by large concave wall curvature is essential since this interaction with each other produces flow field changes which in turn affect the nozzle performance.

Nickerson, G. R.↗

Comparative efficiency of finite, boundary and hybrid element methods in elastostatics

The comparative computational efficiencies of the finite element (FEM), boundary element (BEM), and hybrid boundary element-finite element (HBFEM) analysis techniques are evaluated for representative bounded domain interior and unbounded domain exterior problems in elastostatics. Computational efficiency is carefully defined in this study as the computer time required to attain a specified level of solution accuracy. The study found the FEM superior to the BEM for the interior problem, while the reverse was true for the exterior problem. The hybrid analysis technique was found to be comparable or superior to both the FEM and BEM for both the interior and exterior problems.

Schwartz, C. W.↗

Application of a low order panel method to complex three-dimensional internal flow problems

An evaluation of the ability of a low order panel method to predict complex three-dimensional internal flow fields was made. The computer code VSAERO was used as a basis for the evaluation. Guidelines for modeling internal flow geometries were determined and the effects of varying the boundary conditions and the use of numerical approximations on the solutions accuracy were studied. Several test cases were run and the results were compared with theoretical or experimental results. Modeling an internal flow geometry as a closed box with normal velocities specified on an inlet and exit face provided accurate results and gave the user control over the boundary conditions. The values of the boundary conditions greatly influenced the amount of leakage an internal flow geometry suffered and could be adjusted to eliminate leakage. The use of the far-field approximation to reduce computation time influenced the accuracy of a solution and was coupled with the values of the boundary conditions needed to eliminate leakage. The error induced in the influence coefficients by using the far-field approximation was found to be dependent on the type of influence coefficient, the far-field radius, and the aspect ratio of the panels.

Ashby, D. L.↗

Research related to improved computer aided design software package

The comparative computational efficiencies of the finite element (FEM), boundary element (BEM), and hybrid boundary element-finite element (HVFEM) analysis techniques are evaluated for representative bounded domain interior and unbounded domain exterior problems in elastostatics. Computational efficiency is carefully defined in this study as the computer time required to attain a specified level of solution accuracy. The study found the FEM superior to the BEM for the interior problem, while the reverse was true for the exterior problem. The hybrid analysis technique was found to be comparable or superior to both the FEM and BEM for both the interior and exterior problems.

Walston, W. H., Jr.↗

Comparison of uniform perturbation and numerical solutions for some potential flows past slender bodies

Approximate solutions for potential flow past an axisymmetric slender body and past a thin airfoil are calculated using a uniform perturbation method and then compared with either the exact analytical solution or the solution obtained using a purely numerical method. The perturbation method is based upon a representation of the disturbance flow as the superposition of singularities distributed entirely within the body, while the numerical (panel) method is based upon a distribution of singularities on the surface of the body. It is found that the perturbation method provides very good results for small values of the slenderness ratio and for small angles of attack. Moreover, for comparable accuracy, the perturbation method is simpler to implement, requires less computer memory, and generally uses less computation time than the panel method. In particular, the uniform perturbation method yields good resolution near the regions of the leading and trailing edges where other methods fail or require special attention.

Wong, T.-C.↗

Calculation of recoil implantation profiles using known range statistics

A method has been developed to calculate the depth distribution of recoil atoms that result from ion implantation onto a substrate covered with a thin surface layer. The calculation includes first order recoils considering projected range straggles, and lateral straggles of recoils but neglecting lateral straggles of projectiles. Projectile range distributions at intermediate energies in the surface layer are deduced from look-up tables of known range statistics. A great saving of computing time and human effort is thus attained in comparison with existing procedures. The method is used to calculate recoil profiles of oxygen from implantation of arsenic through SiO2 and of nitrogen from implantation of phosphorus through Si3N4 films on silicon. The calculated recoil profiles are in good agreement with results obtained by other investigators using the Boltzmann transport equation and they also compare very well with available experimental results in the literature. The deviation between calculated and experimental results is discussed in relation to lateral straggles. From this discussion, a range of surface layer thickness for which the method applies is recommended.

Fung, C. D.↗