Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical 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 91 records · Page 5

Numerical computation of space shuttle orbiter flow field

The complete inviscid viscous real gas flow around the Space Shuttle Orbiter was computed. Real gas effects are important in predicting the reentry environment around the Orbiter because the high temperatures within the shock layer cause the air to dissociate and ionize, thus invalidating the perfect gas assumption. It is shown that real gas effects has a significant influence on the Orbiter aerodynamics. The approach utilizes a time dependent Navier-Stokes code to compute the subsonic nose portion of the flow field. This nose solution provides the initial conditions for the parabolized Navier-Stokes (PNS) code.

Tannehill, J. C.↗

Numerical computation of viscous flow around bodies and wings moving at supersonic speeds

Research in aerodynamics is discussed. The development of equilibrium air curve fits; computation of hypersonic rarefield leading edge flows; computation of 2-D and 3-D blunt body laminar flows with an impinging shock; development of a two-dimensional or axisymmetric real gas blunt body code; a study of an over-relaxation procedure forthe MacCormack finite-difference scheme; computation of 2-D blunt body turbulent flows with an impinging shock; computation of supersonic viscous flow over delta wings at high angles of attack; and computation of the Space Shuttle Orbiter flowfield are discussed.

Tannehill, J. C.↗

Numerical computation of space shuttle orbiter flow field

A new parabolized Navier-Stokes (PNS) code has been developed to compute the hypersonic, viscous chemically reacting flow fields around 3-D bodies. The flow medium is assumed to be a multicomponent mixture of thermally perfect but calorically imperfect gases. The new PNS code solves the gas dynamic and species conservation equations in a coupled manner using a noniterative, implicit, approximately factored, finite difference algorithm. The space-marching method is made well-posed by special treatment of the streamwise pressure gradient term. The code has been used to compute hypersonic laminar flow of chemically reacting air over cones at angle of attack. The results of the computations are compared with the results of reacting boundary-layer computations and show excellent agreement.

Tannehill, John C.↗

Downward continuation of the free-air gravity anomalies to the ellipsoid using the gradient solution and terrain correction: An attempt of global numerical computations

The formulas for the determination of the coefficients of the spherical harmonic expansion of the disturbing potential of the earth are defined for data given on a sphere. In order to determine the spherical harmonic coefficients, the gravity anomalies have to be analytically downward continued from the earth's surface to a sphere-at least to the ellipsoid. The goal is to continue the gravity anomalies from the earth's surface downward to the ellipsoid using recent elevation models. The basic method for the downward continuation is the gradient solution (the g sub 1 term). The terrain correction was also computed because of the role it can play as a correction term when calculating harmonic coefficients from surface gravity data. The fast Fourier transformation was applied to the computations.

Wang, Y. M.↗

The validation of numerical computations

Information is given in viewgraph form on the validation of computational fluid dynamics codes. Topics covered include the types of validation required, aerodynamic heating to a slender 5 degree cone, heat transfer on cones with isentropic compression surfaces, the validation of turbulence data, modeling of thermodynamic and transport properties, real gas code validation, the Flight Dynamics Laboratory validation efforts, and reacting gas experimental data in low density flow.

Buck, Melvin L.↗

Numerical computation of complex multi-body Navier-Stokes flows with applications for the integrated Space Shuttle launch vehicle

An enhanced grid system for the Space Shuttle Orbiter was built by integrating CAD definitions from several sources and then generating the surface and volume grids. The new grid system contains geometric components not modeled previously plus significant enhancements on geometry that has been modeled in the old grid system. The new orbiter grids were then integrated with new grids for the rest of the launch vehicle. Enhancements were made to the hyperbolic grid generator HYPGEN and new tools for grid projection, manipulation, and modification, Cartesian box grid and far field grid generation and post-processing of flow solver data were developed.

Chan, William M.↗

Determination of forces in a magnetic bearing actuator - Numerical computation with comparison to experiment

Calculations of the forces exerted on a journal by a magnetic bearing actuator are presented, along with comparisons to experimentally measured forces. The calculations are based on two-dimensional solutions for the flux distribution in the metal parts and free space, using finite but constant permeability in the metals. Above a relative permeability of 10,000 the effects of changes in permeability are negligible, but below 10,000 decreases in permeability cause significant decreases in the force. The calculated forces are shown to depend on the metal permeability more strongly when the journal is displaced from its centered position. The predicted forces in the principal attractive direction are in good agreement with experiment when a relatively low value of permeability is chosen. The forces measured normal to the axis of symmetry when the journal is displaced from that axis, however, are significantly higher than predicted by theory, even with a value of relative permeability larger than 5000. These results indicate a need for further work including nonlinear permeability distributions.

Knight, J. D.↗

Numerical Computation of Low-Speed Concurrent Flow Flame Spread in Mixed Buoyant and Forced Flow

The effect of low-speed mixed convection (forced plus buoyant) on concurrent flow flame spread over a thin solid is examined. Computations are carried out using an existing model. Results indicate that seemingly small levels of gravity can significantly alter flame spread rates. Starting with a purely forced flow condition, as gravity is added, the entrainment due to buoyancy lengthens the flame and increases the spread rate significantly. Buoyancy has an influence on the extinction limits. At low speed, the presence of a small gravity level widens the flammability limit.

Ferkul, Paul V.↗

Numerical computation of aerodynamics and heat transfer in a turbine cascade and a turn-around duct using advanced turbulence models

The objective of this research is to develop turbulence models to predict the flow and heat transfer fields dominated by the curvature effect such as those encountered in turbine cascades and turn-around ducts. A Navier-Stokes code has been developed using an explicit Runge-Kutta method with a two layer k-epsilon/ARSM (Algebraic Reynolds Stress Model), Chien's Low Reynolds Number (LRN) k-epsilon model and Coakley's LRN q-omega model. The near wall pressure strain correlation term was included in the ARSM. The formulation is applied to Favre-averaged N-S equations and no thin-layer approximations are made in either the mean flow or turbulence transport equations. Anisotropic scaling of artificial dissipation terms was used. Locally variable timestep was also used to improve convergence. Detailed comparisons were made between computations and data measured in a turbine cascade by Arts et al. at Von Karman Institute. The surface pressure distributions and wake profiles were predicted well by all the models. The blade heat transfer is predicted well by k-epsilon/ARSM model, as well as the k-epsilon model. It's found that the onset of boundary layer transition on both surfaces is highly dependent upon the level of local freestream turbulence intensity, which is strongly influenced by the streamline curvature. Detailed computation of the flow in the turn around duct has been carried out and validated against the data by Monson as well as Sandborn. The computed results at various streamwise locations both on the concave and convex sides are compared with flow and turbulence data including the separation zone on the inner well. The k-epsilon/ARSM model yielded relatively better results than the two-equation turbulence models. A detailed assessment of the turbulence models has been made with regard to their applicability to curved flows.

Lakshminarayana, B.↗

Numerical Computation of the Chemically Reacting Flow around the National Aero-Space Plane

This final report summarizes the research accomplished. The research performed during the grant period can be divided into the following major areas: (1) Computation of chemically reacting Supersonic combustion ramjet (scramjet) flowfields. (2) Application of a two-equation turbulence model to supersonic combustion flowfields. (3) Computation of the integrated aerodynamic and propulsive flowfields of a generic hypersonic space plane. (4) Computation of hypersonic flows with finite-catalytic walls. (5) Development of an upwind Parabolized Navier-Stokes (PNS) code for thermo-chemical nonequilibrium flows.

Tannehill, J. C.↗

Real time numerical computation of corrosion rates from corrosion sensors

A method for determining a corrosion rate includes receiving an input signal including corrosion data from a sensor. A signal-to-noise ratio (SNR) and a type of noise is determined from the input signal. A corrosion rate computation method is selected based on the SNR and the type of noise. The corrosion rate is determined, using a processor, based on a selected corrosion rate computation method.

36 MATERIALS SCIENCE↗

Multiple burn fuel-optimal orbit transfers: Numerical trajectory computation and neighboring optimal feedback guidance

This report describes current work in the numerical computation of multiple burn, fuel-optimal orbit transfers and presents an analysis of the second variation for extremal multiple burn orbital transfers as well as a discussion of a guidance scheme which may be implemented for such transfers. The discussion of numerical computation focuses on the use of multivariate interpolation to aid the computation in the numerical optimization. The second variation analysis includes the development of the conditions for the examination of both fixed and free final time transfers. Evaluations for fixed final time are presented for extremal one, two, and three burn solutions of the first variation. The free final time problem is considered for an extremal two burn solution. In addition, corresponding changes of the second variation formulation over thrust arcs and coast arcs are included. The guidance scheme discussed is an implicit scheme which implements a neighboring optimal feedback guidance strategy to calculate both thrust direction and thrust on-off times.

Chuang, C.-H.↗

Computational and numerical analysis of AC optimal power flow formulations on large-scale power grids

We report that alternating current optimal power flow (AC-OPF) is a fundamental tool in electric utilities to determine optimal operation of the various resources. Typically, the AC-OPF problem uses power balance formulation containing voltages and power equations. Yet, there is no comprehensive comparison of the different AC-OPF formulations, especially for large-scale networks. This paper presents a detailed comparative evaluation of different formulations of the AC-OPF problem on networks ranging from 9-bus to 25,000 buses. Three different formulations: 1) power balance with polar voltages, 2) power balance with Cartesian voltages, and 3) current balance with Cartesian voltages are discussed in detail by comparing their characteristics, and numerical and computational performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Radiative transfer effects on reflected shock waves. I - Transparent gas.

Analytical and numerical calculation of the effects of radiative cooling on the flow field and heat transfer behind a reflected shock wave. The analytical solutions are developed by an expansion procedure about the Newtonian limit - i.e., the flow variables are expanded in the small parameter epsilon representing the initial density ratio across the shock front. Solutions carried out through the zeroth- and first-orders in epsilon show an accuracy to the order of epsilon squared - namely, for the specific conditions considered the analytically calculated enthalpy changes agree to within 2% with the numerically computed changes. The numerical procedure consists of a finite-difference calculation with pressure gradients neglected. The analytical calculations show that pressure-gradient effects may be neglected, except at very long times when the shocked gas has cooled by a large amount. The present calculations show much greater shock-wave attenuation than does a previous numerical computation. This difference is attributed to a better representation of the thermodynamic properties in the present calculations.

Su, F. Y.↗

External audio for IBM-compatible computers

Numerous applications benefit from the presentation of computer-generated auditory stimuli at points discontiguous with the computer itself. Modification of an IBM-compatible computer for use of an external speaker is relatively easy but not intuitive. This modification is briefly described.

Washburn, David A.↗

Parity-odd and even trispectrum from axion inflation

The four-point correlation function of primordial scalar perturbations has parity-even and parity-odd contributions and the parity-odd signal in cosmological observations is opening a novel window to look for new physics in the inflationary epoch. We study the distinct parity-odd and even prediction from the axion inflation model, in which the inflaton couples to a vector field via a Chern-Simons interaction, and the vector field is considered to be either approximately massless (m A << Hubble scale H) or very massive (m A ~ H). The parity-odd signal arises due to one transverse mode of the vector field being predominantly produced during inflation. We adopt the in-in formalism to evaluate the correlation functions. Considering the vector field mode function to be dominated by its real part up to a constant phase, we simplify the formulas for numerical computations. The numerical studies show that the massive and massless vector fields give significant parity-even signals, while the parity-odd contribution is about one to two orders of magnitude smaller.

79 ASTRONOMY AND ASTROPHYSICS↗