Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer experiments”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Software for Acoustic Rendering

SLAB is a software system that can be run on a personal computer to simulate an acoustic environment in real time. SLAB was developed to enable computational experimentation in which one can exert low-level control over a variety of signal-processing parameters, related to spatialization, for conducting psychoacoustic studies. Among the parameters that can be manipulated are the number and position of reflections, the fidelity (that is, the number of taps in finite-impulse-response filters), the system latency, and the update rate of the filters. Another goal in the development of SLAB was to provide an inexpensive means of dynamic synthesis of virtual audio over headphones, without need for special-purpose signal-processing hardware. SLAB has a modular, object-oriented design that affords the flexibility and extensibility needed to accommodate a variety of computational experiments and signal-flow structures. SLAB s spatial renderer has a fixed signal-flow architecture corresponding to a set of parallel signal paths from each source to a listener. This fixed architecture can be regarded as a compromise that optimizes efficiency at the expense of complete flexibility. Such a compromise is necessary, given the design goal of enabling computational psychoacoustic experimentation on inexpensive personal computers.

Miller, Joel D.↗
NASA NTRS2002

Tactical Defenses Against Systematic Variation in Wind Tunnel Testing

This paper examines the role of unexplained systematic variation on the reproducibility of wind tunnel test results. Sample means and variances estimated in the presence of systematic variations are shown to be susceptible to bias errors that are generally non-reproducible functions of those variations. Unless certain precautions are taken to defend against the effects of systematic variation, it is shown that experimental results can be difficult to duplicate and of dubious value for predicting system response with the highest precision or accuracy that could otherwise be achieved. Results are reported from an experiment designed to estimate how frequently systematic variations are in play in a representative wind tunnel experiment. These results suggest that significant systematic variation occurs frequently enough to cast doubts on the common assumption that sample observations can be reliably assumed to be independent. The consequences of ignoring correlation among observations induced by systematic variation are considered in some detail. Experimental tactics are described that defend against systematic variation. The effectiveness of these tactics is illustrated through computational experiments and real wind tunnel experimental results. Some tutorial information describes how to analyze experimental results that have been obtained using such quality assurance tactics.

DeLoach, Richard↗
NASA NTRS1993

The effect of grid aspect ratio on convergence

Vector stability theory and computational experiments are used to examine the effects of high aspect ratio grids on algorithm convergence. Emphasis is placed on identifying the degree of convergence deterioration that ADI schemes experience as grid aspect ratio increases and on deriving methods for circumventing this problem. Approximately factored implicit schemes are found to experience convergence deterioration because of nonoptimum local time stepping procedures and the increased need for viscous preconditioning. An enhanced algorithm is devised on the basis of improved selection of the local time step, appropriate definition of viscous preconditioning matrix, and proper implementation of boundary conditions. The novel algorithm provides uniformly efficient convergence at all aspect ratios for both Euler and Navier-Stokes computations for a variety of test problems.

Buelow, P. E. O.↗
NASA NTRS1992

Aeroelastic analysis of advanced propellers using an efficient Euler solver

A 3D Euler solver is coupled with a 3D structural dynamics model to investigate flutter of propfans. A hybrid scheme is used to reduce computational time for the Euler equations and a normal mode analysis is used for flutter calculations. Experimental and calculated flutter results are compared for an advanced propeller propfan which experienced flutter at transonic tip relative velocities. The predicted flutter calculations are in close agreement with the experimental data. A structural damping value of 0.5 percent was required to predict the behavior observed in the experiment. Computations show that the flutter behavior is dominated by the second mode, but coupling with the first mode is required. The addition of other modes to the calculations did not affect the flutter behavior.

Srivastava, R.↗
NASA NTRS1983

A versatile ground data display system for Spacelab experiments

A microprocessor-based data capture and display system suited to the needs of many Spacelab experimenters is described. Although developed specifically for the needs of the FAUST UV astronomy telescope, several features of the system may make it attractive to other groups who need to interact with their Spacelab data in real or near-real time. The display system captures data from the serial data stream originating in the Spacelab experiment computer input/output (ECIO) data channel, and can reformat and display several hundred variables per second in real time. The system offers a combination of CRT display and analog paper strip-chart output. A facility is also provided for simultaneously creating floppy-disk records of short segments of data according to prearranged criteria, to allow detailed examination of critical events in the operation of the experiment. The simplicity, economy, reliability, and flexibility of the display system suggest its use for monitoring a variety of low- and medium-data-rate Spacelab experiments.

Werthimer, D.↗
NASA NTRS1994

Effect of grid aspect ratio on convergence

The effects of high-aspect-ratio grids on algorithm convergence are considered by means of vector stability theory and computational experiments. The results indicate that approximately factored implicit schemes experience convergence deterioration because of nonoptimum local time-stepping procedures and increased need for viscous preconditioning. Based on this insight, an enhanced algorithm is devised using improved selection of the local time step, appropriate definition of the viscous preconditioning matrix, and proper implementation of the boundary conditions. The new algorithm provides uniformly efficient convergence at all aspect ratios for both Euler and Navier-Stokes computations for a variety of test problems.

Buelow, P. E. O.↗
NASA NTRS2019

Numerical Investigation and Optimization of a Flushwall Injector for Scramjet Applications at Hypervelocity Flow Conditions

An investigation utilizing Reynolds-averaged simulations (RAS) was performed in order to demonstrate the use of design and analysis of computer experiments (DACE) methods in Sandia’s DAKOTA software package for surrogate modeling and optimization. These methods were applied to a flow- path fueled with an interdigitated flushwall injector suitable for scramjet applications at hyper- velocity conditions and ascending along a constant dynamic pressure flight trajectory. The flight Mach number, duct height, spanwise width, and injection angle were the design variables selected to maximize two objective functions: the thrust potential and combustion efficiency. Because the RAS of this case are computationally expensive, surrogate models are used for optimization. To build a surrogate model a RAS database is created. The sequence of the design variables comprising the database were generated using a Latin hypercube sampling (LHS) method. A methodology was also developed to automatically build geometries and generate structured grids for each design point. The ensuing RAS analysis generated the simulation database from which the two objective functions were computed using a one-dimensionalization (1D) of the three-dimensional simulation data. The data were fitted using four surrogate models: an artificial neural network (ANN), a cubic polynomial, a quadratic polynomial, and a Kriging model. Variance-based decomposition showed that both objective functions were primarily driven by changes in the duct height. Multiobjective design optimization was performed for all four surrogate models via a genetic algorithm method. Optimal solutions were obtained at the upper and lower bounds of the flight Mach number range. The Kriging model predicted an optimal solution set that exhibited high values for both objective functions. Additionally, three challenge points were selected to assess the designs on the Pareto fronts. Further sampling among the designs of the Pareto fronts may be required to lower the surrogate model errors and perform more accurate surrogate-model-based optimization.

Shenoy, Rajiv R.↗
NASA NTRS1997

An Application of the Difference Potentials Method to Solving External Problems in CFD

Numerical solution of infinite-domain boundary-value problems requires some special techniques that would make the problem available for treatment on the computer. Indeed, the problem must be discretized in a way that the computer operates with only finite amount of information. Therefore, the original infinite-domain formulation must be altered and/or augmented so that on one hand the solution is not changed (or changed slightly) and on the other hand the finite discrete formulation becomes available. One widely used approach to constructing such discretizations consists of truncating the unbounded original domain and then setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The role of the ABC's is to close the truncated problem and at the same time to ensure that the solution found inside the finite computational domain would be maximally close to (in the ideal case, exactly the same as) the corresponding fragment of the original infinite-domain solution. Let us emphasize that the proper treatment of artificial boundaries may have a profound impact on the overall quality and performance of numerical algorithms. The latter statement is corroborated by the numerous computational experiments and especially concerns the area of CFD, in which external problems present a wide class of practically important formulations. In this paper, we review some work that has been done over the recent years on constructing highly accurate nonlocal ABC's for calculation of compressible external flows. The approach is based on implementation of the generalized potentials and pseudodifferential boundary projection operators analogous to those proposed first by Calderon. The difference potentials method (DPM) by Ryaben'kii is used for the effective computation of the generalized potentials and projections. The resulting ABC's clearly outperform the existing methods from the standpoints of accuracy and robustness, in many cases noticeably speed up the multigrid convergence, and at the same time are quite comparable to other methods from the standpoints of geometric universality and simplicity of implementation.

Ryaben 'Kii, Victor S.↗
NASA NTRS1992

The effect of secular resonances in the asteroid region between 2.1 and 2.4 AU

The asteroid region between 2.1 and 2.4 AU appears to be depopulated at inclinations i greater than 12 deg. This region is surrounded by the three main secular resonances nu(sub 5), nu(sub 6), and nu(sub 16) and is crossed by higher order secular resonances. Secular resonances appear to overlap in this region. Numerical integrations of the orbits of seventeen fictituous asteroids with initial inclinations 12 deg less than or equal to i less than or equal to 20 deg show the following: (1) this particular asteroid region is not depopulated in our computer experiment on timescales of 2.7 Myrs; (2) inclinations are pumped up by successive crossings through higher order secular resonances while eccentricities are not increased sufficiently to produce planet-crossers; (3) bodies located in the bordering nu(sub 6) resonance with semi-major axes a less than or equal to 2.4 AU become Earth-crossers on a time scale of 1 Myr; and (4) we confirm the result that modes due to higher order secular resonances must be eliminated when proper elements are computed.

Froeschle, Christiane↗
NASA NTRS1987

Goertler instability on an airfoil

An effective computational scheme was developed to study the growth/damping of Goertler vortices along walls of variable curvature. Computational experiments indicate that when the amplification rates for the u-, v-, and w-perturbations are the same, the finite difference approach to solve the initial value problem and the normal mode approach give identical results for the Blasius boundary layer on constant curvature concave walls. The growth of Goertler vortices was rapid in the concave regions and was followed by sharp damping in the convex region. However, multiple sets of counter-rotating vortices were formed and remained far downstream in the convex region. The current computational scheme can be easily extended to more realistic problems including variable pressure gradients and suction effects.

Kalburgi, Vijay↗
NASA NTRS1972

Comparison of some approximations for isotropic turbulence.

Study of several related turbulence approximations with regard to dynamical properties and agreement of numerical predictions with laboratory and computer experiments. The approximations considered include the direct-interaction equations (Kraichnan, 1964), Herring's (1966) self-consistent-field theory, a generalization of Edwards' (1964) theory, the abridged Lagrangian-history, direct-interaction approximation (Kraichnan, 1966), the test-field model (Kraichnan, 1971), and an approximation, not previously described, in which one velocity field passively suffers convection by another. Most of the cited approximations are representable by stochastic model equations for the velocity amplitude. Explicit constructions are given for the stochastic models, in a form that can be approximated on a digital computer. These constructions are used to discuss the physical and mathematical differences between the model dynamics and actual Navier-Stokes dynamics.-

Herring, J. R.↗
NASA NTRS1992

Optimization and analysis of large chemical kinetic mechanisms using the solution mapping method - Combustion of methane

A method of systematic optimization, solution mapping, as applied to a large-scale dynamic model is presented. The basis of the technique is parameterization of model responses in terms of model parameters by simple algebraic expressions. These expressions are obtained by computer experiments arranged in a factorial design. The developed parameterized responses are then used in a joint multiparameter multidata-set optimization. A brief review of the mathematical background of the technique is given. The concept of active parameters is discussed. The technique is applied to determine an optimum set of parameters for a methane combustion mechanism. Five independent responses - comprising ignition delay times, pre-ignition methyl radical concentration profiles, and laminar premixed flame velocities - were optimized with respect to thirteen reaction rate parameters. The numerical predictions of the optimized model are compared to those computed with several recent literature mechanisms. The utility of the solution mapping technique in situations where the optimum is not unique is also demonstrated.

Frenklach, Michael↗
NASA NTRS1989

Computational results for the effects of external disturbances on transition location of bodies of revolution from subsonic to supersonic speeds and comparisons with experimental data

Computational experiments have been performed for a few configurations in order to investigate the effects of external flow disturbances on the extent of laminar flow and wake drag. Theoretical results have been compared with experimental data for the AEDC cone, for Mach numbers from subsonic to supersonic, and for both free flight and wind tunnel environments. The comparisons have been found to be very satisfactory, thus establishing the utility of the present method for the design and development of laminar flow configurations and for the assessment of wind tunnel data. In addition, results of calculations concerning the effects of unit Reynolds numbers on transition are presented. In addition to the AEDC cone, computations have been performed for an ogive body of revolution at zero angle of attack and supersonic Mach numbers. Results are presented for transition Reynolds number and wake drag for external disturbances corresponding to free air and the test section of the AEDC-VKF tunnel. These results have been found to compare quite well with wind tunnel data for cases when surface suction is applied as well as when suction is absent.

Goradia, S. H.↗
NASA NTRS2003

X-33 Rev-F Turbulent Aeroheating Results From Test 6817 in NASA Langley 20-Inch Mach 6 Air Tunnel and Comparisons With Computations

Measurements and predictions of the X-33 turbulent aeroheating environment have been performed at Mach 6, perfect-gas air conditions. The purpose of this investigation was to compare measured turbulent aeroheating levels on smooth models, models with discrete trips, and models with arrays of bowed panels (which simulate bowed thermal protections system tiles) with each other and with predictions from two Navier-Stokes codes, LAURA and GASP. The wind tunnel testing was conducted at free stream Reynolds numbers based on length of 1.8 x 10(exp 6) to 6.1 x 10(exp 6) on 0.0132 scale X-33 models at a = 40-deg. Turbulent flow was produced by the discrete trips and by the bowed panels at ill but the lowest Reynolds number, but turbulent flow on the smooth model was produced only at the highest Reynolds number. Turbulent aeroheating levels on each of the three model types were measured using global phosphor thermography and were found to agree to within .he estimated uncertainty (plus or minus 15%) of the experiment. Computations were performed at the wind tunnel free stream conditions using both codes. Turbulent aeroheating levels predicted using the LAURA code were generally 5%-10% lower than those from GASP, although both sets of predictions fell within the experimental accuracy of the wind tunnel data.

Hollis, Brian R.↗
NASA NTRS2019

Surrogate Modeling and Optimization of a Combustor with an Interdigitated Flushwall Injector

Design and Analysis of Computer Experiments (DACE) methods are applied and used to perform surrogate modeling and optimization of a simplified combustor flowpath with an interdigitated flushwall injector. The objectives of the optimization are the thrust potential and combustion efficiency, which are evaluated across a range of flight Mach numbers, duct heights, spanwise spacings, and injection angles. The focus of this work is to highlight the application of a sequential learning approach, in order to learn about the responses of the objective functions over the design space and to identify local regions of interest for further analysis. This approach is contrasted to a previous effort where only a single sampling set was used to fit surrogate models and perform optimization. The optimal solutions resulting from the previous and present approaches are different, due to surrogate model-guided local refinement of the design space allowed by the sequential learning method. The values of the global error estimates between the previous and present approaches are comparable, but the sequential method proved more computationally cost-effective. Further refinement in the optimal regions might be needed to improve predictive capability of surrogate models and to obtain the optimal solutions sets.

Rajiv R Shenoy↗
NASA NTRS1995

Unsteady numerical simulations of the stability and dynamics of flames

In this report we describe the research performed at the Naval Research Laboratory in support of the NASA Microgravity Science and Applications Program over the past three years (from Feb. 1992) with emphasis on the work performed since the last microgravity combustion workshop. The primary objective of our research is to develop an understanding of the differences in the structure, stability, dynamics and extinction of flames in earth gravity and in microgravity environments. Numerical simulations, in which the various physical and chemical processes can be independently controlled, can significantly advance our understanding of these differences. Therefore, our approach is to use detailed time-dependent, multi-dimensional, multispecies numerical models to perform carefully designed computational experiments. The basic issues we have addressed, a general description of the numerical approach, and a summary of the results are described in this report. More detailed discussions are available in the papers published which are referenced herein. Some of the basic issues we have addressed recently are (1) the relative importance of wall losses and gravity on the extinguishment of downward-propagating flames; (2) the role of hydrodynamic instabilities in the formation of cellular flames; (3) effects of gravity on burner-stabilized flames, and (4) effects of radiative losses and chemical-kinetics on flames near flammability limits. We have also expanded our efforts to include hydrocarbon flames in addition to hydrogen flames and to perform simulations in support of other on-going efforts in the microgravity combustion sciences program. Modeling hydrocarbon flames typically involves a larger number of species and a much larger number of reactions when compared to hydrogen. In addition, more complex radiation models may also be needed. In order to efficiently compute such complex flames recent developments in parallel computing have been utilized to develop a state-of-the-art parallel flame code. This is discussed below in some detail after a brief discussion of the numerical models.

Kailasanath, K.↗