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.
A new numerical method used to drastically reduce the computation time required to solve the Navier-Stokes equations at flight Reynolds numbers is described. The new method makes it possible and practical to calculate many important three-dimensional, high Reynolds number flow fields on computers.
A reliable method is presented for calculating the flowfield about a cascade of arbitrary 2-D airfoils. The method approximates the three-dimensional flow in a turbomachinery blade row by correcting for streamtube convergence and radius change in the throughflow direction. The method is a fully conservative solution of the full potential equation incorporating the finite volume technique on a body-fitted periodic mesh, with an artificial density imposed in the transonic region to ensure stability and the capture of shock waves. Comparison of results for several supercritical blades shows good agreement with their hodograph solutions. Other calculations for these profiles as well as standard NACA blade sections indicate that this is a useful scheme for analyzing both the design and off-design performance of turbomachinery blading.
Interferometry methods were applied to the investigation of steady and unsteady flows in large scale transonic wind tunnels. Holographic interferometry was demonstrated to provide reliable flow visualization and quantitative results for a number of two-dimensional flows. These conclusions were based on extensive comparisons with results obtained by other means. Data obtained on a NACA 64A010 airfoil with an oscillating flap installed in the Ames 11-foot transonic tunnel are presented. Interferograms were recorded at a free stream Mach number of 0.8, flap frequency of 30 Hertz and chord Reynolds numbers of 6.6 x 10 to the 6th and 12.3 x 10 to the 6th. The interferometric results were reduced to dynamic surface pressures, Mach contours and wake flow profiles. A new interferometry method that is capable of providing real-time interferometry data is also discussed.
Compressible Euler equations are solved for two-dimensional problems by a preconditioned conjugate gradient-like technique. An approximate Riemann solver is used to compute the numerical fluxes to second order accuracy in space. Two ways to achieve parallelism are tested, one which makes use of parallelism inherent in triangular solves and the other which employs domain decomposition techniques. The vectorization/parallelism in triangular solves is realized by the use of a recording technique called wavefront ordering. This process involves the interpretation of the triangular matrix as a directed graph and the analysis of the data dependencies. It is noted that the factorization can also be done in parallel with the wave front ordering. The performances of two ways of partitioning the domain, strips and slabs, are compared. Results on Cray YMP are reported for an inviscid transonic test case. The performances of linear algebra kernels are also reported.
An analysis method for predicting the inplane stress states in anisotropic finite plates with an elliptical notch is presented. This method can be used to analyze plates with arbitrary notch orientations with respect to the plate material axes. The analysis results have been validated using finite element analysis results for unnotched composite plates and experimental and finite element analysis results for stiffened composite panels with a skin that has orthotropic properties. The good agreement between these results, until the panel exhibits nonlinear response either due to bending or initiation of damage, indicates that the present analysis method can be used to determine accurately the inplane stress states and stress concentrations at and around an elliptical notch.
An optimization method has been developed to design axisymmetric body shapes such as fuselages, nacelles, and external fuel tanks with increased transition Reynolds numbers in subsonic compressible flow. The new design method involves a constraint minimization procedure coupled with analysis of the inviscid and viscous flow regions, and linear stability analysis of the compressible boundary layer. Boundary-layer transition is predicted by a 'hybrid' transition criterion based on Granville's transition criterion and a criterion using linear stability theory coupled with the e(n)-method. A tiptank of a business-jet is used as an example to illustrate that the method can be utilized to design an axisymmetric body shape with extensive natural laminar flow. Boundary layer transition is predicted to occur at a transition Reynolds number of 6.04 x 10 to the 6th on the original tiptank. On the designed body shape, a transition Reynolds number of 7.22 x 10 to the 6th is predicted using compressible linear stability theory coupled with e(n)-method.
This paper describes matrix methods that have been developed for calculating compressible flow on a blade-to-blade surface of revolution. The methods have been fully tested to date only for the design of plane cascades to prescribed blade surface distributions; the methods will be illustrated here for that problem only. Similar methods are presently being applied to both the direct and indirect problems and for flow on arbitrary surfaces of revolution in annular cascades with stream sheet thickness variations. It is believed that by such methods, both the direct and indirect calculations can be reduced to about 60 to 90 seconds of computing.
Development of the plasma sprayed graded, layered ZRO2/CoCrAlY seal system for gas turbine engine blade tip seal applications up to 1589 K (2400 F) surface temperature was continued. The effect of changing ZRO2/CoCrAlY ratios in the intermediate layers on thermal stresses was evaluated analytically with the goal of identifying the materials combinations which would minimize thermal stresses in the seal system. Three methods of inducing compressive residual stresses in the sprayed seal materials to offset tensile thermal stresses were analyzed. The most promising method, thermal prestraining, was selected based upon potential, feasibility and complexity considerations. The plasma spray equipment was modified to heat, control and monitor the substrate temperature during spraying. Specimens were fabricated and experimentally evaluated to: (1) substantiate the capability of the thermal prestrain method to develop compressive residual stresses in the sprayed structure and (2) define the effect of spraying on a heated substate on abradability, erosion and thermal shock characteristics of the seal system. Thermal stress analysis, including residual stresses and material properties variations, was performed and correlated with thermal shock test results. Seal system performance was assessed and recommendations for further development were made.
An optimization method has been developed to design axisymmetric body shapes such as fuselages, nacelles, and external fuel tanks with increased transition Reynolds numbers in subsonic compressible flow. The new design method involves a constraint minimization procedure coupled with analysis of the inviscid and viscous flow regions and linear stability analysis of the compressible boundary-layer. In order to reduce the computer time, Granville's transition criterion is used to predict boundary-layer transition and to calculate the gradients of the objective function, and linear stability theory coupled with the e(exp n)-method is used to calculate the objective function at the end of each design iteration. Use of a method to design an axisymmetric body with extensive natural laminar flow is illustrated through the design of a tiptank of a business jet. For the original tiptank, boundary layer transition is predicted to occur at a transition Reynolds number of 6.04 x 10(exp 6). For the designed body shape, a transition Reynolds number of 7.22 x 10(exp 6) is predicted using compressible linear stability theory coupled with the e(exp n)-method.
Currently, the dominating constraint in many high performance computing applications is data capacity and bandwidth, in both internode communications and even moreso in intranode data motion. A new approach to address this limitation is to make use of data compression in the form of a compressed data array. Storing data in a compressed data array and converting to standard IEEE-754 types as needed during a computation can reduce the pressure on bandwidth and storage. However, repeated conversions (lossy compression and decompression) introduce additional approximation errors, which need to be shown to not significantly affect the simulation results. Here, we extend recent work that analyzed the error of a single use of compression and decompression of the ZFP compressed data array representation to the case of time-stepping and iterative schemes, where an advancement operator is repeatedly applied in addition to the conversions. We show that the accumulated error for iterative methods involving fixed-point and time evolving iterations is bounded under standard constraints. An upper bound is established on the number of additional iterations required for the convergence of stationary fixed-point iterations. An additional analysis of traditional forward and backward error of stationary iterative methods using ZFP compressed arrays is also presented. The results of several 1D, 2D, and 3D test problems are provided to demonstrate the correctness of the theoretical bounds.
A fine-mesh method incorporating two new operators, which drastically reduces the computation time, has been developed for solving the time-dependent Navier-Stokes equations at flight Reynolds numbers. The approach time-splits the equations into a hyperbolic part and a parabolic part, solves the hyperbolic part by a new explicit numerical method based on characteristics theory, and solves the parabolic part by a new efficient implicit parabolic method. The method has reduced the computation time by one and two orders of magnitude from that required previously to solve for the interaction of a shock wave with a boundary layer on a flat plate.
A method for predicting the wall-normal velocity in compressible boundary layer flow is presented which is the extension of the fully implicit, spectral collocation method developed by Pruett and Streett (1991). The spectrally accurate method under consideration is capable of extracting the wall-normal velocity from the transformation from the physical plane to a computational plane, without interpolation, leaving the continuity equation free as a check on the quality of the solution. When used in conjunction with a highly accurate spectral collocation method for the compressible boundary-layer equations, the method gives a discrete solution which satisfies the continuity equation nearly to machine precision.
We present a general framework for compressing unstructured scientific data with known local connectivity. A common application is simulation data defined on arbitrary finite element meshes. The framework employs a greedy topology preserving reordering of original nodes which allows for seamless integration into existing data processing pipelines. This reordering process depends solely on mesh connectivity and can be performed offline for optimal efficiency. However, the algorithm’s greedy nature also supports on-the-fly implementation. The proposed method is compatible with any compression algorithm that leverages spatial correlations within the data. The effectiveness of this approach is demonstrated on a large-scale real dataset using several compression methods, including MGARD, SZ, and ZFP.
Three test schemes were examined for testing graphite/epoxy (Narmco T300/5208) composite material specimens to failure in compression, including an adaptation of the IITRI "wedge grip" compression fixture, a face-supported-compression fixture, and an end-loaded-coupon fixture. The effects of specimen size, specimen support arrangement and method of load transfer on compressive behavior of graphite/epoxy were investigated. Compressive stress strain, strength, and modulus data obtained with the three fixtures are presented with evaluations showing the effects of all test parameters, including fiber orientation. The IITRI fixture has the potential to provide good stress/strain data to failure for unidirectional and quasi-isotropic laminates. The face supported fixture was found to be the most desirable for testing + or - 45 s laminates.
Systems, apparatuses, and methods for implementing a multi-tiered approach to cache compression are disclosed. A cache includes a cache controller, light compressor, and heavy compressor. The decision on which compressor to use for compressing cache lines is made based on certain resource availability such as cache capacity or memory bandwidth. This allows the cache to opportunistically use complex algorithms for compression while limiting the adverse effects of high decompression latency on system performance. To address the above issue, the proposed design takes advantage of the heavy compressors for effectively reducing memory bandwidth in high bandwidth memory (HBM) interfaces as long as they do not sacrifice system performance. Accordingly, the cache combines light and heavy compressors with a decision-making unit to achieve reduced off-chip memory traffic without sacrificing system performance.
Compressible, inviscid secondary injection and primary flows calculated using integral relation method for blunt bodies - application to thrust vector control
The objective of this work is to compare a high-order solver with a low-order solver for performing Large-Eddy Simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the highorder method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.
The objective of this work is to compare a high-order solver with a low-order solver for performing large-eddy simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the high-order method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.