Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Analysis and Computation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Stability Analysis of Streaks Induced by Optimized Vortex Generators

Numerical computations are performed to investigate the potential for transition control in an axisymmetric boundary layer via fully realizable, streamwise stationary streaks induced by an azimuthally periodic array of surface mounted vortex generators (VGs). Previous work has shown that suitable streaks of this type can significantly reduce the growth of Mack’s second mode instabilities, but large streak amplitudes can make the flow susceptible to previously absent streak instabilities that can become the leading cause of transition. Here, we use the adjoint capabilities of the SU2 flow solver to optimize the VG shape to maximize the reduction in the growth of second-mode disturbances while also preventing the streak amplitudes from reaching large enough values to precipitate an earlier onset of transition via streak instabilities. The geometry and the freestream flow conditions are selected to match a relevant trajectory lo-cation from the HIFiRE-1 flight experiment. Results show that the optimized VGs can increase the mean streak amplitude by 117% with respect to a manually developed baseline design. The stability of this optimized basic state is analyzed via the plane-marching parabolized stability equations, predicting a fully laminar flow over the entire cone, or equivalently, yielding transition delay of 130% versus the 17% for the baseline VGs.

Boundary layer transition↗
NASA NTRS2022

Stability Analysis of Streaks Induced By Optimized Vortex Generators

Numerical computations are performed to investigate the potential for transition control in an axisymmetric boundary layer via fully realizable, streamwise stationary streaks induced by an azimuthally periodic array of surface mounted vortex generators (VGs). Previous work has shown that suitable streaks of this type can significantly reduce the growth of Mack’s second mode instabilities, but large streak amplitudes can make the flow susceptible to previously absent streak instabilities that can become the leading cause of transition. Here, we use the adjoint capabilities of the SU2 flow solver to optimize the VG shape to maximize the reduction in the growth of second-mode disturbances while also preventing the streak amplitudes from reaching large enough values to precipitate an earlier onset of transition via streak instabilities. The geometry and the freestream flow conditions are selected to match a relevant trajectory lo-cation from the HIFiRE-1 flight experiment. Results show that the optimized VGs can increase the mean streak amplitude by 117% with respect to a manually developed baseline design. The stability of this optimized basic state is analyzed via the plane-marching parabolized stability equations, predicting a fully laminar flow over the entire cone, or equivalently, yielding transition delay of 130% versus the 17% for the baseline VGs.

Boundary layer transition↗
NASA NTRS1972

Flow field predictions for a slab delta wing at incidence

Theoretical results are presented for the structure of the hypersonic flow field of a blunt slab delta wing at moderately high angle of attack. Special attention is devoted to the interaction between the boundary layer and the inviscid entropy layer. The results are compared with experimental data. The three-dimensional inviscid flow is computed numerically by a marching finite difference method. Attention is concentrated on the windward side of the delta wing, where detailed comparisons are made with the data for shock shape and surface pressure distributions. Surface streamlines are generated, and used in the boundary layer analysis. The three-dimensional laminar boundary layer is computed numerically using a specially-developed technique based on small cross-flow in streamline coordinates. In the rear sections of the wing the boundary layer decreases drastically in the spanwise direction, so that it is still submerged in the entropy layer at the centerline, but surpasses it near the leading edge. Predicted heat transfer distributions are compared with experimental data.

Conti, R. J.↗
NASA NTRS1972

Testing relativistic gravity theories using radio tracking data from planetary orbiting spacecraft.

We present a thorough analysis of a computational method for determining the numerical values of the relativity and other related dynamical parameters using two-way Doppler and ranging data from planetary orbiting spacecraft. The computational method consists of two parts. From Doppler data we first determine the earth-planet components of the position of the orbiting spacecraft relative to the center of gravity of the planet to high accuracy; adding the observed spacecraft range yields a range value to the center of the planet. These constructed earth-planet range data, referred to as normal points, are then treated as raw data in a regression analysis combined with planetary radar delay and meridian circle measurements to solve for the significant solar system dynamical parameters. The major errors sources in the planetary orbiter process are enumerated and their individual effects on the overall accuracy of the normal point accuracies are presented. The accuracies of the parameter estimates as a function of time, data sampling, and a priori assumptions are illustrated.

Jordan, J. F.↗
NASA NTRS1978

Computational flow development for unsteady viscous flows: Foundation of the numerical method

A procedure is presented for effective consideration of viscous effects in computational development of high Reynolds number flows. The procedure is based on the interpretation of the Navier-Stokes equations as vorticity transport equations. The physics of the flow was represented in a form suitable for numerical analysis. Lighthill's concept for flow development for computational purposes was adapted. The vorticity transport equations were cast in a form convenient for computation. A statement for these equations was written using the method of weighted residuals and applying the Galerkin criterion. An integral representation of the induced velocity was applied on the basis of the Biot-Savart law. Distribution of new vorticity, produced at wing surfaces over small computational time intervals, was assumed to be confined to a thin region around the wing surfaces.

Bratanow, T.↗
NASA NTRS2016

Direct Numerical Simulation of an Airfoil with Sand Grain Roughness on the Leading Edge

As part of a computational study of acoustic radiation due to the passage of turbulent boundary layer eddies over the trailing edge of an airfoil, the Lattice-Boltzmann method is used to perform direct numerical simulations of compressible, low Mach number flow past an NACA 0012 airfoil at zero degrees angle of attack. The chord Reynolds number of approximately 0.657 million models one of the test conditions from a previous experiment by Brooks, Pope, and Marcolini at NASA Langley Research Center. A unique feature of these simulations involves direct modeling of the sand grain roughness on the leading edge, which was used in the abovementioned experiment to trip the boundary layer to fully turbulent flow. This report documents the findings of preliminary, proof-of-concept simulations based on a narrow spanwise domain and a limited time interval. The inclusion of fully-resolved leading edge roughness in this simulation leads to significantly earlier transition than that in the absence of any roughness. The simulation data is used in conjunction with both the Ffowcs Williams-Hawkings acoustic analogy and a semi-analytical model by Roger and Moreau to predict the farfield noise. The encouraging agreement between the computed noise spectrum and that measured in the experiment indicates the potential payoff from a full-fledged numerical investigation based on the current approach. Analysis of the computed data is used to identify the required improvements to the preliminary simulations described herein.

Ribeiro, Andre F. P.↗
NASA NTRS1980

Analysis and testing of numerical formulas for the initial value problem

Three computer programs for evaluating and testing numerical integration formulas used with fixed stepsize programs to solve initial value systems of ordinary differential equations are described. A program written in PASCAL SERIES, takes as input the differential equations and produces a FORTRAN subroutine for the derivatives of the system and for computing the actual solution through recursive power series techniques. Both of these are used by STAN, a FORTRAN program that interactively displays a discrete analog of the Liapunov stability region of any two dimensional subspace of the system. The derivatives may be used by CLMP, a FORTRAN program, to test the fixed stepsize formula against a good numerical result and interactively display the solutions.

Brown, R. L.↗
NASA NTRS1986

A general numerical analysis program for the superconducting quasiparticle mixer

A user-oriented computer program SISCAP (SIS Computer Analysis Program) for analyzing SIS mixers is described. The program allows arbitrary impedance terminations to be specified at all LO harmonics and sideband frequencies. It is therefore able to treat a much more general class of SIS mixers than the widely used three-frequency analysis, for which the harmonics are assumed to be short-circuited. An additional program, GETCHI, provides the necessary input data to program SISCAP. The SISCAP program performs a nonlinear analysis to determine the SIS junction voltage waveform produced by the local oscillator. The quantum theory of mixing is used in its most general form, treating the large signal properties of the mixer in the time domain. A small signal linear analysis is then used to find the conversion loss and port impedances. The noise analysis includes thermal noise from the termination resistances and shot noise from the periodic LO current. Quantum noise is not considered. Many aspects of the program have been adequately verified and found accurate.

Hicks, R. G.↗
NASA NTRS1997

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

Applications are described of high-performance computing methods to the numerical simulation of complete jet engines. The methodology focuses on the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field elements. New partitioned analysis procedures to treat this coupled three-component problem were developed. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. The NASA-sponsored ENG10 program was used for the global steady state analysis of the whole engine. This program uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 was developed as well as the capability for the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames.

Felippa, C. A.↗
NASA NTRS1982

Numerical analysis of confined turbulent flow

The considered investigation is concerned with the development of an efficient computational method for obtaining a physical understanding of an internal turbulent field. The employed approach makes use of a 'two equation' type model for the turbulence to obtain the numerical solution of a two-dimensional confined turbulent flow. The mean flow governing equations are considered along with the governing equation of the mean temperature and concentrations, and the boundary conditions. The numerical procedure for solving the turbulent flow is discussed, taking into account an approximation to the nonlinear terms, and the inner and outer coupling. Attention is given to a stability convergence analysis, the stability characteristics, and computational examples.

Lin, A.↗
NASA NTRS1989

Structural zooming research and development of an interactive computer graphical interface for stress analysis of cracks

Engineering problems sometimes involve the numerical solution of boundary value problems over domains containing geometric feature with widely varying scales. Often, a detailed solution is required at one or more of these features. Small details in large structures may have profound effects upon global performance. Conversely, large-scale conditions may effect local performance. Many man-hours and CPU-hours are currently spent in modeling such problems. With the structural zooming technique, it is now possible to design an integrated program which allows the analyst to interactively focus upon a small region of interest, to modify the local geometry, and then to obtain highly accurate responses in that region which reflect both the properties of the overall structure and the local detail. A boundary integral equation analysis program, called BOAST, was recently developed for the stress analysis of cracks. This program can accurately analyze two-dimensional linear elastic fracture mechanics problems with far less computational effort than existing finite element codes. An interactive computer graphical interface to BOAST was written. The graphical interface would have several requirements: it would be menu-driven, with mouse input; all aspects of input would be entered graphically; the results of a BOAST analysis would be displayed pictorially but also the user would be able to probe interactively to get numerical values of displacement and stress at desired locations within the analysis domain; the entire procedure would be integrated into a single, easy to use package; and it would be written using calls to the graphic package called HOOPS. The program is nearing completion. All of the preprocessing features are working satisfactorily and were debugged. The postprocessing features are under development, and rudimentary postprocessing should be available by the end of the summer. The program was developed and run on a VAX workstation, and must be ported to the SUN workstation. This activity is currently underway.

Gerstle, Walter↗
DOE OSTI2025

Boosting efficiency and reducing graph reliance: Basis adaptation integration in Bayesian multi-fidelity networks

The computational cost of high-fidelity numerical models makes outer-loop analysis, which requires repeated interrogation of the model such as uncertainty quantification, computationally demanding. Multi-fidelity methods, which construct a surrogate model using data from an ensemble of models of varying cost and accuracy, can substantially reduce the cost of outer-loop analysis. However, these methods can be difficult to apply when the model ensemble does not admit a clear hierarchy a priori and the correlations between models are low. Consequently, in this paper, we present a multi-fidelity method that leverages dimension reduction to enhance the correlation between models, thereby reducing the amount of data needed to train a surrogate from an unordered ensemble of models. Our method utilizes basis adaptation to build low-dimensional polynomial chaos expansions of each model and employs Multi-fidelity Networks to encode the relationships among models. We show that the resulting method exhibit two notable advantages over its counterpart: (1) enhanced accuracy (both reduced bias and variance); and (2) reduced dependency on the graph structure encoding relationships among models. We demonstrate the approach on an analytical test problem and a challenging finite element model for a spent nuclear fuel. Our method produces a surrogate model that is significantly more accurate than either a single-fidelity surrogate or a multi-fidelity surrogate constructed without basis adaptation.

42 ENGINEERING↗
NASA NTRS1977

Synchrotron emissivity near the electron cyclotron and upper hybrid frequencies

The spontaneous synchrotron emissivity from a high density, weakly relativistic plasma in thermal equilibrium is discussed. Thermal effects on the dielectric polarization of the plasma are included in the analysis, and the result is compared with computations of the emissivity based upon the cold plasma approximation. A numerical analysis is performed for frequencies in the vicinity of both the electron cyclotron frequency and the upper hybrid resonance. Significant modifications to the synchrotron emissivity in the cold plasma approximation are found in both of these frequency regimes.

Freund, H. P.↗
NASA NTRS1981

Numerical analysis of the scramjet inlet flow field using two-dimensional Navier-Stokes equations

A computer code has been developed to solve the full two-dimensional Navier-Stokes equations in a scramjet inlet. The analysis uses a numerical coordinate transformation which generates a set of boundary fitted curvilinear coordinates. The explicit finite-difference algorithm of MacCormack is used to solve the governing equations. A two-layer eddy viscosity model is used for the turbulent flow. The code can analyze both inviscid and viscous flows with multiple struts in the flow field. Detailed results are presented for two model problems and two scramjet inlets with one and two struts. The application of the two-dimensional analysis in the preliminary design of the actual scramjet inlet is briefly discussed.

Kumar, A.↗
NASA NTRS1973

Numerical analysis of stiffened shells of revolution. Volume 1: Theory manual for STARS-2S, 2B, 2V digital computer programs

The theoretical analysis background for the STARS-2 (shell theory automated for rotational structures) program is presented. The theory involved in the axisymmetric nonlinear and unsymmetric linear static analyses, and the stability and vibrations (including critical rotation speed) analyses involving axisymmetric prestress are discussed. The theory for nonlinear static, stability, and vibrations analyses, involving shells with unsymmetric loadings are included.

Svalbonas, V.↗
NASA NTRS1975

Numerical nonlinear inelastic analysis of stiffened shells of revolution. Volume 1: Theory manual for STARS-2P digital computer program

The theoretical analysis background for the STARS-2P nonlinear inelastic program is discussed. The theory involved is amenable for the analysis of large deflection inelastic behavior in axisymmetric shells of revolution subjected to axisymmetric loadings. The analysis is capable of considering such effects as those involved in nonproportional and cyclic loading conditions. The following are also discussed: orthotropic nonlinear kinematic hardening theory; shell wall cross sections and discrete ring stiffeners; the coupled axisymmetric large deflection elasto-plastic torsion problem; and the provision for the inelastic treatment of smeared stiffeners, isogrid, and waffle wall constructions.

Svalbonas, V.↗
NASA NTRS1992

Numerical analysis of the hot-gas-side and coolant-side heat transfer in liquid rocket engine combustors

The objectives of this paper are to develop computational methods to predict the hot-gas-side and coolant-side heat transfer, and to use these methods in parametric studies to recommend optimized design of the coolant channels for regeneratively cooled liquid rocket engine combustors. An integrated numerical model which incorporates computational fluid dynamics (CFD) for the hot-gas thermal environment, and thermal analysis for the coolant channels, was developed. The mode was validated by comparing predicted heat fluxes with those of hot-firing test and industrial design methods. Parametric studies were performed to find a strategy for optimized combustion chamber coolant channel design.

Wang, T. S.↗