Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer performance”

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

Ascent Heating Thermal Analysis on Spacecraft Adaptor Fairings

When the Crew Exploration Vehicle (CEV) is launched, the spacecraft adaptor (SA) fairings that cover the CEV service module (SM) are exposed to aero heating. Thermal analysis is performed to compute the fairing temperatures and to investigate whether the temperatures are within the material limits for nominal ascent aeroheating case. The ascent heating is analyzed by using computational fluid dynamics (CFD) and engineering codes at Marshall Space Flight Center. The aeroheating environment data used for this work is known as Thermal Environment 3 (TE3) heating data. One of the major concerns is with the SA fairings covering the CEV SM and the SM/crew launch vehicle (CLV) flange interface. The TE3 heating rate is a function of time, wall temperature, and the spatial locations. The implementation of the TE3 heating rate as boundary conditions in the thermal analysis becomes challenging. The ascent heating thermal analysis on SA fairings and SM/CLV flange interface are performed using two commercial software packages: Cullimore & Ring (C&R) Thermal Desktop (TD) 5.1 and MSC Patran 2007r1 b. TD is the pre-and post-processor for SINDA, which is a finite-difference-based solver. In TD, the geometry is built and meshed, the boundary conditions are defined, and then SINDA is used to compute temperatures. MSC Pthermal is a finite-element- based thermal solver. MSC Patran is the pre- and post-processor for Pthermal. Regarding the boundary conditions, the convection, contact resistance, and heat load can be imposed in different ways in both programs. These two software packages are used to build the thermal model for the same analysis to validate each other and show the differences in the modeling details.

Wang, Xiao Yen↗

Task III: Development of an Effective Computational Methodology for Body Force Representation of High-speed Rotor 37

A framework for an effective computational methodology for characterizing the stability and the impact of distortion in high-speed multi-stage compressor is being developed. The methodology consists of using a few isolated-blade row Navier-Stokes solutions for each blade row to construct a body force database. The purpose of the body force database is to replace each blade row in a multi-stage compressor by a body force distribution to produce same pressure rise and flow turning. To do this, each body force database is generated in such a way that it can respond to the changes in local flow conditions. Once the database is generated, no hrther Navier-Stokes computations are necessary. The process is repeated for every blade row in the multi-stage compressor. The body forces are then embedded as source terms in an Euler solver. The method is developed to have the capability to compute the performance in a flow that has radial as well as circumferential non-uniformity with a length scale larger than a blade pitch; thus it can potentially be used to characterize the stability of a compressor under design. It is these two latter features as well as the accompanying procedure to obtain the body force representation that distinguish the present methodology from the streamline curvature method. The overall computational procedures have been developed. A dimensional analysis was carried out to determine the local flow conditions for parameterizing the magnitudes of the local body force representation of blade rows. An Euler solver was modified to embed the body forces as source terms. The results from the dimensional analysis show that the body forces can be parameterized in terms of the two relative flow angles, the relative Mach number, and the Reynolds number. For flow in a high-speed transonic blade row, they can be parameterized in terms of the local relative Mach number alone.

Tan, Choon-Sooi↗

Efficient spares matrix multiplication scheme for the CYBER 203

This work has been directed toward the development of an efficient algorithm for performing this computation on the CYBER-203. The desire to provide software which gives the user the choice between the often conflicting goals of minimizing central processing (CPU) time or storage requirements has led to a diagonal-based algorithm in which one of three types of storage is selected for each diagonal. For each storage type, an initialization sub-routine estimates the CPU and storage requirements based upon results from previously performed numerical experimentation. These requirements are adjusted by weights provided by the user which reflect the relative importance the user places on the resources. The three storage types employed were chosen to be efficient on the CYBER-203 for diagonals which are sparse, moderately sparse, or dense; however, for many densities, no diagonal type is most efficient with respect to both resource requirements. The user-supplied weights dictate the choice.

Lambiotte, J. J., Jr.↗

Vortex ring bubbles

Toroidal bubbles with circulation are studied numerically and by means of a physically motivated model equation. Two series of computations are performed by a boundary-integral method. One set shows the starting motion of an initially spherical bubble as a gravitationally driven liquid jet penetrates through the bubble from below causing a toroidal geometry to develop. The jet becomes broader as surface tension increases and fails to penetrate if surface tension is too large. The dimensionless circulation that develops is not very dependent on the surface tension. The second series of computations starts from a toroidal geometry, with circulation determined from the earlier series, and follows the motion of the rising and spreading vortex ring. Some modifications to the boundary-integral formulation were devised to handle the multiply connected geometry. The computations uncovered some unexpected rapid oscillations of the ring radius. These oscillations and the spreading of the ring are explained by the model equation which provides a more general description of vortex ring bubbles than previously available.

Lundgren, T. S.↗

Laminar-Turbulent Transition Upstream of the Entropy-Layer Swallowing Location in Hypersonic Boundary Layers

Numerical and experimental studies have demonstrated that modal growth of planar Mack modes is responsible for laminar-turbulent transition on sharp cones at hypersonic speeds. However, the physical mechanisms that lead to transition onset upstream of the entropy-layer swallowing location over sufficiently blunt geometries are not well understood as yet. Modal amplification is too weak or nonexistent to initiate transition at moderate-to-large bluntness values. Nonmodal analysis shows that, with increasing nose bluntness, both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification. However, because of the relatively weak signature of the nonmodal traveling disturbances within the boundary-layer region, the route to transition onset subsequent to the nonmodal growth remains unclear. Thus, nonlinear parabolized stability equations (NPSE) and direct numerical simulations (DNS) have been used to investigate the potential transition mechanisms over a 7-degree blunt cone that was tested in the AFRL Mach-6 high-Reynoldsnumber facility. Computations are performed to separately follow the nonlinear development of two classes of inflow disturbances, namely, a pair of oblique traveling waves with equal but opposite angles with respect to the mean flow direction and a planar traveling wave. Results in both cases show an excellent agreement between the NPSE and DNS predictions, establishing that the NPSE is an accurate and efficient technique for predicting the nonlinear development for these particular nonmodal traveling disturbances. Computations reveal that the oblique mode interactions lead to the generation of stationary streaks inside the boundary layer that, in turn, facilitate the growth of a subharmonic sinuous disturbance. For relatively modest amplitudes of the inflow disturbance, the oblique-mode breakdown can lead to transition at the measured location of transition onset during the experiment. On the other hand, the nonlinear development of a planar traveling wave leads to the formation of inclined structures just above the boundary-layer edge and these structures are strongly reminiscent of the transitional events observed during blunt cone experiments by using schlieren flow visualizations.

Entropy↗

Time-accurate Navier-Stokes computations of self-excited two-dimensional unsteady cavity flows

Time-accurate Navier Stokes computations were performed to study self-excited oscillatory shear layer behavior over a two-dimensional rectangular cavity and to gain insight into the vortex-acoustic coupling that sustains these oscillations. Results were obtained solving the full, unsteady compressible Navier-Stokes equations for Reynolds-averaged turbulent air flow past a cavity with L/d = 2. Conditions duplicated test conditions in experiments by Rossiter at Mach numbers varied from 0.6 to 1.2. The computational scheme was the USA 2D series code developed by Rockwell's Science Center with an algebraic turbulence model for the shear layer and recirculating region. Two physical conditions were simulated: an isolated cavity situated (1) in an infinite plane wall and (2) in an infinite channel with an opposing parallel wall. Computed results for (1) showed excellent agreement with Rossiter's measured amplitudes and frequencies. A tuned acoustic duct resonance effect was observed in (2), demonstrating a controlled change in the acoustic feedback mechanism and resultant amplification in computed oscillations.

Dougherty, N. Sam↗

Effect of Distributed Patch of Smooth Roughness Elements on Transition in a High-Speed Boundary Layer

Surface roughness is known to have a substantial impact on the aerothermodynamic loading of hypersonic vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, smooth-shaped roughness elements on the laminar boundary layer over a 7-degree half-angle, circular cone for flow conditions corresponding to a selected trajectory point from the ascent phase of the HIFiRE-1 flight experiment. For peak-to-valley roughness heights of 50 percent or less in comparison with the thickness of the unperturbed boundary layer, the computations converge to a stationary flow, suggesting that the flow is globally stable. Analysis of convective instabilities in the wake of the roughness patch indicates two dominant families of unstable disturbances, namely, a high frequency mode that corresponds to Mack mode waves modified by the wake and a lower frequency mode that corresponds to shear layer instabilities associated with the streaks in the roughness wake. Even though the peak growth rate of the later mode is more than 35 percent greater than the peak growth rates of the Mack modes, the latter modes achieve higher amplification ratios, and hence, are likely to dominate the onset of transition, which is estimated to occur slightly later than that in the unperturbed, i.e., smooth surface boundary layer. Additional computations are performed to investigate the effects of various roughness patch configurations on a Mach 3.5 flat plate boundary layer, to help guide an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at NASA Langley Research Center. In this case, the cumulative reinforcement of basic state distortion over the length of the roughness patch is predicted to yield a significantly earlier transition than that over a smooth plate or a plate with a shorter length roughness patch.

Meelan Choudhari↗

Analyzing a 35-Year Hourly Data Record: Why So Difficult?

At the Goddard Distributed Active Archive Center, we have recently added a 35-Year record of output data from the North American Land Assimilation System (NLDAS) to the Giovanni web-based analysis and visualization tool. Giovanni (Geospatial Interactive Online Visualization ANd aNalysis Infrastructure) offers a variety of data summarization and visualization to users that operate at the data center, obviating the need for users to download and read the data themselves for exploratory data analysis. However, the NLDAS data has proven surprisingly resistant to application of the summarization algorithms. Algorithms that were perfectly happy analyzing 15 years of daily satellite data encountered limitations both at the algorithm and system level for 35 years of hourly data. Failures arose, sometimes unexpectedly, from command line overflows, memory overflows, internal buffer overflows, and time-outs, among others. These serve as an early warning sign for the problems likely to be encountered by the general user community as they try to scale up to Big Data analytics. Indeed, it is likely that more users will seek to perform remote web-based analysis precisely to avoid the issues, or the need to reprogram around them. We will discuss approaches to mitigating the limitations and the implications for data systems serving the user communities that try to scale up their current techniques to analyze Big Data.

computational performance↗

SSME LOX post flow analysis/fluid structure interaction. Volume 1: Flow analysis

The realization of measures to improve the performance of the Space Shuttle is, to a large extent, dependent on an improved understanding of the fluid flow phenomena occurring in the main engine. The overall arrangement of the primary components of the Space Shuttle Main Engine (SSME) are presented. The impingement of the hot gases from the transfer ducts onto the LOX posts causes them to be subjected to severe thermal and gas dynamic loads, which in the past have resulted in the occasional breakage of some elements of the outer row of posts during test firings of the engine, particularly at higher power levels. Large velocities in the gap between the LOX posts may also be a contributing factor in causing dynamic stability problems. The deforming structural response of the posts to the pressure loading may likely affect the gas flowfield by producing a moving flowfield boundary, thereby creating a dynamically coupled unsteady fluid-structure system. The objective was to investigate the three-dimensional, turbulent flow around a simplified SSME LOX post array using an existing Reynolds averaged Navier-Stokes flow solver and a suitable turbulence model to parameterize the turbulent shear stresses. Numerical computations were performed to analyze the effect on the flowfield of varying the spacing between the LOX posts, which were modeled as rigid, three-dimensional circular cylinders. The methodology used in the computations is described. Results are presented.

Burke, Roger W.↗

PURDU-WINCOF: A computer code for establishing the performance of a fan-compressor unit with water ingestion

A model for predicting the performance of a multi-spool axial-flow compressor with a fan during operation with water ingestion was developed incorporating several two-phase fluid flow effects as follows: (1) ingestion of water, (2) droplet interaction with blades and resulting changes in blade characteristics, (3) redistribution of water and water vapor due to centrifugal action, (4) heat and mass transfer processes, and (5) droplet size adjustment due to mass transfer and mechanical stability considerations. A computer program, called the PURDU-WINCOF code, was generated based on the model utilizing a one-dimensional formulation. An illustrative case serves to show the manner in which the code can be utilized and the nature of the results obtained.

Leonardo, M.↗

Aerodynamic shape optimization via sensitivity analysis on decomposed computational domains

Direct and iterative method considered to be most applicable to large systems of linear equations arising in discrete sensitivity analysis are assessed. Based on a single-domain grid, computations are performed using a banded matrix solver and an iterative solver, the generalized minimum residual (GMRES) method. The banded matrix solver is found to be generally the most economical method for those applications where the number of right-hand sides is large (i.e., a large number of design variables or a large number of adjoint vectors). For systems of equations that are too large to be solved by direct methods, an approach is proposed whereby the computational domain is divided into small subdomains, and each subdomain is solved separately.

Eleshaky, Mohamed E.↗

Prediction of STS-107 Hypervelocity Flow Fields about the Shuttle Orbiter with Various Wing Leading Edge Damage

Computations were performed for damaged configurations of the Shuttle Orbiter in support of the STS-107 Columbia accident investigation. Two configurations with missing wing leading-edge reinforced carbon-carbon (RCC) panels were evaluated at conditions just prior to the peak heating trajectory point. The initial configuration modeled the Orbiter with an approximate missing RCC panel 6 to determine whether this damage could result in anomalous temperatures measured during the STS-107 reentry. This missing RCC panel 6 computation was found to produce heating augmentation factors of 5 times the nominal heating rates on the side fuselage with lesser heat increases on the front of the OMS pod. This is consistent with the thermocouple and resistance temperature detector sensors from the STS-107 re-entry which observed off nominal high early in the re-entry trajectory. A second damaged configuration modeled the Orbiter with missing RCC panel 9 and included ingestion of the flow into the outboard RCC channel. This computation lowered the level (only 2 times nominal) and moved the location of the heating augmentation on the leeside fuselage relative to the missing RCC panel 6 configuration. The lesser heating augmentation for missing RCC panel 9 was confined near the wing fuselage juncture. Near nominal heating was predicted on the remainder of the side fuselage with some lower than nominal heating on the front surface of the OMS pod. These results for missing RCC panel 9 are consistent with data from the STS-107 re-entry where the heating augmentation was observed to move off the side fuselage and OMS pod sensors at later times in the trajectory. As this solution requires supersonic mass ingestion into the RCC channel, it is probably not an appropriate model prior to penetration of the flow through the spar into the wing structure. It may, however, be representative of the conditions at later times and could account for the movement of the heating signature on the side fuselage.

Pulsonetti, Maria V.↗

EOS MLS Level 1B Data Processing, Version 2.2

A computer program performs level- 1B processing (the term 1B is explained below) of data from observations of the limb of the Earth by the Earth Observing System (EOS) Microwave Limb Sounder (MLS), which is an instrument aboard the Aura spacecraft. This software accepts, as input, the raw EOS MLS scientific and engineering data and the Aura spacecraft ephemeris and attitude data. Its output consists of calibrated instrument radiances and associated engineering and diagnostic data. [This software is one of several computer programs, denoted product generation executives (PGEs), for processing EOS MLS data. Starting from level 0 (representing the aforementioned raw data, the PGEs and their data products are denoted by alphanumeric labels (e.g., 1B and 2) that signify the successive stages of processing.] At the time of this reporting, this software is at version 2.2 and incorporates improvements over a prior version that make the code more robust, improve calibration, provide more diagnostic outputs, improve the interface with the Level 2 PGE, and effect a 15-percent reduction in file sizes by use of data compression.

Perun, Vincent↗

FDNS code to predict wall heat fluxes or wall temperatures in rocket nozzles

This report summarizes the findings on the NASA contract NAG8-212, Task No. 3. The overall project consists of three tasks, all of which have been successfully completed. In addition, some supporting supplemental work, not required by the contract, has been performed and is documented herein. Task 1 involved the modification of the wall functions in the code FDNS to use a Reynolds Analogy-based method. Task 2 involved the verification of the code against experimentally available data. The data chosen for comparison was from an experiment involving the injection of helium from a wall jet. Results obtained in completing this task also show the sensitivity of the FDNS code to unknown conditions at the injection slot. Task 3 required computation of the flow of hot exhaust gases through the P&W 40K subscale nozzle. Computations were performed both with and without film coolant injection. The FDNS program tends to overpredict heat fluxes, but, with suitable modeling of backside cooling, may give reasonable wall temperature predictions. For film cooling in the P&W 40K calorimeter subscale nozzle, the average wall temperature is reduced from 1750 R to about 1050 R by the film cooling. The average wall heat flux is reduced by a factor of three.

Karr, Gerald R.↗

Fiber pushout test: A three-dimensional finite element computational simulation

A fiber pushthrough process was computationally simulated using three-dimensional finite element method. The interface material is replaced by an anisotropic material with greatly reduced shear modulus in order to simulate the fiber pushthrough process using a linear analysis. Such a procedure is easily implemented and is computationally very effective. It can be used to predict fiber pushthrough load for a composite system at any temperature. The average interface shear strength obtained from pushthrough load can easily be separated into its two components: one that comes from frictional stresses and the other that comes from chemical adhesion between fiber and the matrix and mechanical interlocking that develops due to shrinkage of the composite because of phase change during the processing. Step-by-step procedures are described to perform the computational simulation, to establish bounds on interfacial bond strength and to interpret interfacial bond quality.

Mital, Subodh K.↗

Fiber pushout test - A three-dimensional finite element computational simulation

A fiber pushthrough process was computationally simulated using three-dimensional finite element method. The interface material is replaced by an anisotropic material with greatly reduced shear modulus in order to simulate the fiber pushthrough process using a linear analysis. Such a procedure is easily implemented and is computational very effective. It can be used to predict fiber pushthrough load for a composite system at any temperature. The average interface shear strength obtained from pushthrough load can easily be separated into its two components: one that comes from frictioal stresses and the other that comes from chemical adhesion between fiber and the matrix and mechanical interlocking that develops due to shrinkage of the composite because of phase change during the processing. Step-by-step procedures are described to perform the computational simulation, to establish bounds on interfacial bond strength and to interpret interfacial bond quality.

Mital, Subodh K.↗

The nondeterministic divide

The nondeterministic divide partitions a vector into two non-empty slices by allowing the point of division to be chosen nondeterministically. Support for high-level divide-and-conquer programming provided by the nondeterministic divide is investigated. A diva algorithm is a recursive divide-and-conquer sequential algorithm on one or more vectors of the same range, whose division point for a new pair of recursive calls is chosen nondeterministically before any computation is performed and whose recursive calls are made immediately after the choice of division point; also, access to vector components is only permitted during activations in which the vector parameters have unit length. The notion of diva algorithm is formulated precisely as a diva call, a restricted call on a sequential procedure. Diva calls are proven to be intimately related to associativity. Numerous applications of diva calls are given and strategies are described for translating a diva call into code for a variety of parallel computers. Thus diva algorithms separate logical correctness concerns from implementation concerns.

Charlesworth, Arthur↗

High performance architecture for robot control

Practical aspects of the design and implementation of a modular, high performance, parallel computer control system for telerobots are discussed. Topics of consideration include system architecture, operator interface, and control execution. In a laboratory environment, a telerobotics test control configuration is used to obtain measurements on communications and control loop timing for use in an effective full scale operational system design. The feasibility of the selected architectural approach has been successfully demonstrated. The modularity of the software and hardware enables ease of transport for use in the operational system. The distributed partioning of the control algorithms and the performance measurements acquired during control system implementation are discussed.

Byler, E.↗