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.
Computer program incorporates technique that reduces variance of forward Monte Carlo method for given amount of computer time in determining radiation environment in complex organic and inorganic systems exposed to significant amounts of radiation.
Flight data from five aircraft were used to assess the effects of sampling rate and record length reductions on estimates of stability and control derivatives produced by a maximum likelihood estimation method. Derivatives could be extracted from flight data with the maximum likelihood estimation method even if there were considerable reductions in sampling rate and/or record length. Small amplitude pulse maneuvers showed greater degradation of the derivative maneuvers than large amplitude pulse maneuvers when these reductions were made. Reducing the sampling rate was found to be more desirable than reducing the record length as a method of lessening the total computation time required without greatly degrading the quantity of the estimates.
An approach to the solution of matrix problems resulting from integral equations of mathematical physics is presented. Based on the inherent smoothness in such equations, the problem is reformulated using a set of orthogonal basis vectors, leading to an equivalent coefficient problem which can be of lower order without significantly impairing the accuracy of the solution. This approach was evaluated using a two-dimensional Neumann problem describing the inviscid, incompressible flow over an airfoil. Two different kinds of mode functions were investigated, namely eigenfunction series and Fourier series. The method using Fourier series was found preferable. It uses all of the coefficients from a Fast Fourier Transform algorithm in an approximate method which exploits the known structure of the transformed coefficient matrix and very promising results for the flow over a realistic airfoil are obtained. On the basis of the results presented here, an order of magnitude reduction in this computer time can be expected for such problems as compared with the time for a direct matrix solution.
The author has identified the following significant results. The minimum entropy model may provide several useful advantages over traditional techniques for processing LANDSAT data. Total computer time to conduct a complete pattern recognition process is reduced. Subjective (transformed image), as well as statistically derived information is made available to the analyst/user much earlier in the analysis process. A rapid feedback loop in which numerous training set combinations can be tested for difference and representativeness is available. Additional tests of LANDSAT data processing using the minimum entropy model are clearly justified.
A numerical minimization scheme is used in conjunction with two-dimensional and three-dimensional inviscid transonic flow analysis codes to provide procedures for wing leading edge aerodynamic design. The procedures are demonstrated in the design of a new leading edge to improve C-141 cruise performance. For the high aspect ratio moderately swept C-141 wing, the 2-D procedure is shown to yield results which are in close agreement with those obtained using the 3-D technique. Although the 2-D approach uses much less computation time than the 3-D technique, the latter requires fewer manhours than the former. Comparisons of predicted and wind tunnel measured performance improvements are presented which verify the design procedures.
Lift interference effects are discussed based on Bland's (1968) integral equation. A mathematical existence theory is utilized for which convergence of the numerical method has been proved for general (square-integrable) downwashes. Airloads are computed using orthogonal airfoil polynomial pairs in conjunction with a collocation method which is numerically equivalent to Galerkin's method and complex least squares. Convergence exhibits exponentially decreasing error with the number n of collocation points for smooth downwashes, whereas errors are proportional to 1/n for discontinuous downwashes. The latter can be reduced to 1/n to the m+1 power with mth-order Richardson extrapolation (by using m = 2, hundredfold error reductions were obtained with only a 13% increase of computer time). Numerical results are presented showing acoustic resonance, as well as the effect of Mach number, ventilation, height-to-chord ratio, and mode shape on wind-tunnel interference. Excellent agreement with experiment is obtained in steady flow, and good agreement is obtained for unsteady flow.
This paper describes a procedure for the automatic iteration of an inverse boundary-layer technique to a prescribed pressure distribution in a separated flow. The technique is demonstrated by computation of two transonic airfoil flows and two externally generated shock boundary-layer interaction flows. These results are compared to experimental data and to solutions of the full Navier-Stokes equations. These comparisons indicate that substantial economies can be obtained by applying methods like the present, in lieu of full Navier-Stokes methods, in zonal calculation schemes for design purposes. The optimization technique leading to convergence is described in detail and a table of typical computation time is presented.
A computer program using a variation of the Monte Carlo method was developed to predict radiative exchange factors in three-dimensional configurations with curved surfaces. The model allows for arbitrary emission, absorption, and reflection characteristics and for radiation exchange between surfaces as well as external irradiation through openings. It was shown that exchange factors can be determined efficiently by using the normal Monte Carlo method for closed configurations and/or high surface emissivities, while energy partitioning is vastly superior in the case of open configurations and/or low surface emissivities. Comparison with some experiments performed on the Space Shuttle heat rejection system demonstrates the model's accuracy as well as its superiority over a simplified diffuse-surface analysis, which requires similar amounts of computer time.
A numerical scheme is presented for the computation of transonic aerodynamic loads on flexible wings. The method consists of iteratively applying the loads computed by a 3D transonic aerodynamics code to a structural model to obtain elastic twist, and then recomputing the loads. Because this iteration is performed concurrently with the iterations performed in computing the aerodynamics, flexible loads are obtained in roughly the same amount of computing time as required to obtain rigid loads. Applications of this method to a flexible supercritical transonic transport wing are presented and compared with model test data.
A simple thermal analysis program for Space Shuttle payloads has been developed to accommodate the user who requires an easily understood but dependable analytical tool. The thermal analysis program includes several thermal subprograms traditionally employed in spacecraft thermal studies, a data management system for data generated by the subprograms, and a master program to coordinate the data files and thermal subprograms. The language and logic used to run the thermal analysis program are designed for the small user. In addition, analytical and storage techniques which conserve computer time and minimize core requirements are incorporated into the program.
Two- and three-dimensional turbulent wakes which develop under the influence of streamline curvature and rotation were calculated with three turbulence closure models. The first model comprised of transport equations for the turbulent kinetic energy and the rate of energy dissipation. The second model comprised of equations for the rate of turbulent kinetic energy dissipation and Reynolds stresses, but the effects of the convection and diffusion in the Reynolds stress transport equation were handled collectively. The third model utilizes equations of turbulent kinetic energy dissipation and Reynolds stresses in nearly exact form. All of the three models were modified for the effects of streamline curvature. The second and third models include the effects of rotation through the rotation-originated redistribution term in the transport equation of Reynolds stresses. The numerical results for the wakes of airfoils, cascades and turbomachinery rotor blades demonstrate that the second and third models provide accurate predictions, but computer time and storage can be considerably saved with the second model.
Numerical studies clarifying the advantages and disavantages of conventional versus upwind convective finite elements are presented along with lumped versus consistent formulations for practical conduction forced-convection analysis. A finite-element procedure for treatment of negligible capacitance fluid nodes is presented. The procedure is based on procedures used in finite-element structural dynamics to treat nodes with negligible structural mass. Two finite-element programs and a finite-difference lumped-parameter program used in the studies are discussed. Evaluation studies utilizing three convection and two combined conduction-convection problems are then presented and discussed. Additionally, the computational time saving offered by the finite element procedure is considered for a practical combined conduction-convection problem.
There is increased interest in, and concern about, deposition and corrosion phenomena in combustion systems containing inorganic condensible vapors and particles (salts, ash). To meet the need for a computationally tractable deposition rate theory general enough to embrace multielement/component situations of current and future gas turbine and magnetogasdynamic interest, a multicomponent chemically 'frozen' boundary layer (CFBL) deposition theory is presented and its applicability to the special case of Na2SO4 deposition from seeded laboratory burner combustion products is demonstrated. The coupled effects of Fick (concentration) diffusion and Soret (thermal) diffusion are included, along with explicit corrections for effects of variable properties and free stream turbulence. The present formulation is sufficiently general to include the transport of particles provided they are small enough to be formally treated as heavy molecules. Quantitative criteria developed to delineate the domain of validity of CFBL-rate theory suggest considerable practical promise for the present framework, which is characterized by relatively modest demands for new input information and computer time.
The paper presents the evaluation of algorithms having potential for reducing computer time to integrate large systems of equations governing nonlinear transient heat transfer in complex structures. A series of test problems is being developed to test two algorithms. The Shuttle Orbiter frame model having 190 nodes was the subject of an analytical/experimental study. Plans are to perform calculations for the complete wing model having about 4000 nodes as well as a regional (plug) model having about 200 nodes.
Cluster Compression Algorithm (CCA) preprocesses Landsat image data immediately following satellite data sensor (receiver). Data are reduced by extracting pertinent image features and compressing this result into concise format for transmission to ground station. This results in narrower transmission bandwidth, increased data-communication efficiency, and reduced computer time in reconstructing and analyzing image. Similar technique could be applied to other types of recorded data to cut costs of transmitting, storing, distributing, and interpreting complex information.
The aerodynamic characteristics of wings with leading-edge vortex separation were predicted using a method based on a flow model with free vortex elements which are allowed to merge into a concentrated core. The calculated pressure distribution is more accurate than that predicted by methods with discrete vortex filaments alone. In addition, the computer time is reduced approximately by half.
For a nearly Gaussian transmitted pulse shape scattered from a nearly Gaussian distributed sea surface, a small argument series expansion of one term lead to a several term power series expression for the mean return waveform. Specific expressions are given for the first four terms. These results, which require less computer time than numerical convolution, are useful for data analysis from current or past radar altimeters and for design studies of future systems. Several representative results are presented for an idealized Seasat-1 radar altimeter.
A computer program to calculate the radiation properties of the reflector antennas is presented. It can be used for paraboloidal, spherical, or ellipsoidal reflector surfaces and is easily modified to handle any surface that can be expressed analytically. The program is general enough to allow any arbitrary location and pointing angle for the feed antenna. The effect of blockage due to the feed horn is also included in the computations. The computer program is based upon the technique of tracing the rays from the feed antenna to the reflector to an aperture plane. The far field radiation properties are then calculated by performing a double integration over the field points in the aperture plane. To facilitate the computation of double intergral, the field points are first aligned along the equispaced straight lines in the aperture plane. The computation time is relatively insensitive to the absolute size of the aperture and even though no limits on the largest reflector size have been determined, the program was used for reflector diameters of 1000 wavelenghts.