SEARCH · Search NASA

Results for “Higher order calculations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1974

Improved solution for potential flow about arbitrary axisymmetric bodies by the use of a higher-order surface source method. Part 2. User's manual for computer program

The surface-source method of calculating potential flow is improved by refining the underlying numerical analysis. The present analysis uses parabolic elements and linearly-varying source density which results in a large increase in computing speed and accuracy. The computer program including all relevant input and output is described.

Friedman, D. M.↗

NASA NTRS1989

Numerical simulation of rotorcraft

The objective of the research is to develop and validate accurate, user-oriented viscous CFD codes (with inviscid options) for three-dimensional, unsteady aerodynamic flows about arbitrary rotorcraft configurations. Unsteady, three-dimensional Euler and Navier-Stokes codes are developed, adapted, and extended to rotor-body combinations. Flow solvers are coupled with zonal grid topologies, including rotating and nonrotating blocks. Special grid clustering and wave-fitting techniques were developed to capture low-level radiating acoustic waves. Significant progress was made in computing the propagation of acoustic waves due to the interaction of a concentrated vortex and a helicopter airfoil. The need for higher-order schemes was firmly established in relatively inexpensive two-dimensional calculations. In three dimensions, the number of grid points required to capture the low-level acoustic waves becomes very large, so that large supercomputer memory becomes essential. Good agreement was obtained between the numerical results obtained with a thin-layer Navier-Stokes code and experimental data from a model rotor. In addition, several nonrotating configurations that are sometimes proposed to simulate rotor blade tips in conventional wind tunnels were examined, and the complex flow around the radical tip shape of the world's fastest helicopter is under investigation. These studies demonstrate the flexibility and power of CFD to gain physical insight, study novel ideas, and examine various possibilities that might be difficult or impossible to set up in physical experiments. As a prelude to studies of rotor-body aerodynamic interactions, a preliminary grid topology and moving-interface strategy were developed. A new Euler/Navier-Stokes code using these techniques computes the vortical wake directly, rather than modeling it, as in most previous rotorcraft studies. Several hover cases were run for conventional and advanced-geometry blades. Numerical schemes using multi-zones and/or adaptive grids appear to be necessary to simulate the complex vortical flows in rotor wakes.

Mccroskey, William J.↗

NASA NTRS1974

Analysis of finite-size phased arrays of circular waveguide elements

A derivation is presented for the calculation of the interelement mutual coupling in a finite-size planar array of waveguide-fed apertures covered by a multilayered dielectric and/or plasma. The general mutual admittance expression is evaluated for circular apertures and the mutual coupling calculations are verified experimentally for two transverse electric circular waveguide mode excited apertures. A parametric study of higher order mode aperture fields indicates that the only significant change in the circular aperture mutual coupling is due to the transverse magnetic mode, which introduces an additional phase shift. Qualitative agreement between calculations for a 183-element array of circular apertures and an infinite array establishes the validity of the finite-array theoretical model.

Bailey, M. C.↗

NASA NTRS2004

Increasing Accuracy in Computed Inviscid Boundary Conditions

A technique has been devised to increase the accuracy of computational simulations of flows of inviscid fluids by increasing the accuracy with which surface boundary conditions are represented. This technique is expected to be especially beneficial for computational aeroacoustics, wherein it enables proper accounting, not only for acoustic waves, but also for vorticity and entropy waves, at surfaces. Heretofore, inviscid nonlinear surface boundary conditions have been limited to third-order accuracy in time for stationary surfaces and to first-order accuracy in time for moving surfaces. For steady-state calculations, it may be possible to achieve higher accuracy in space, but high accuracy in time is needed for efficient simulation of multiscale unsteady flow phenomena. The present technique is the first surface treatment that provides the needed high accuracy through proper accounting of higher-order time derivatives. The present technique is founded on a method known in art as the Hermitian modified solution approximation (MESA) scheme. This is because high time accuracy at a surface depends upon, among other things, correction of the spatial cross-derivatives of flow variables, and many of these cross-derivatives are included explicitly on the computational grid in the MESA scheme. (Alternatively, a related method other than the MESA scheme could be used, as long as the method involves consistent application of the effects of the cross-derivatives.) While the mathematical derivation of the present technique is too lengthy and complex to fit within the space available for this article, the technique itself can be characterized in relatively simple terms: The technique involves correction of surface-normal spatial pressure derivatives at a boundary surface to satisfy the governing equations and the boundary conditions and thereby achieve arbitrarily high orders of time accuracy in special cases. The boundary conditions can now include a potentially infinite number of time derivatives of surface-normal velocity (consistent with no flow through the boundary) up to arbitrarily high order. The corrections for the first-order spatial derivatives of pressure are calculated by use of the first-order time derivative velocity. The corrected first-order spatial derivatives are used to calculate the second- order time derivatives of velocity, which, in turn, are used to calculate the corrections for the second-order pressure derivatives. The process as described is repeated, progressing through increasing orders of derivatives, until the desired accuracy is attained.

Dyson, Roger↗

NASA NTRS1984

Calculation of transonic flow in a linear cascade

Turbomachinery blade designs are becoming more aggressive in order to achieve higher loading and greater range. New analysis tools are required to cope with these heavily loaded blades that may operate with a thin separated region near the trailing edge on the suction surface. An existing, viscous airfoil code was adapted to cascade conditions in an attempt to provide this capability. Comparisons with recently obtained data show that calculated and experimental surface Mach numbers were in good agreement but loss coefficients and outlet air angles were not.

Donovan, L. F.↗

NASA NTRS1984

Calculation of transonic flow in a linear cascade

Donovan, L. F.↗

NASA NTRS1969

Rotational relaxation in parahydrogen and its mixtures with helium, neon, and argon at 300 deg K.

Ultrasonic‐velocity dispersion measurements have been performed in parahydrogen and its mixtures with helium, neon, and argon, all at 300°K. In each case the experimental dispersion curves can be matched successfully to those calculated under the assumption that the 0–2 rotational transition relaxes separately from the 2–4 and higher‐order terms. For pure pH 2 we find a relaxation time 𝜏 20 of 1.30 × 10 −8 sec for the 2 → 0 transition and a 𝜏 42 of 3.90 × 10 −8 sec for the 4 → 2 transition. Comparison with the quantum‐mechanical theories of Roberts and of Davison for H 2 –H 2 collisions using Morse potentials shows good agreement for 𝜏 20 over the temperature range of 75°–300°K. The Morse‐potential asymmetry parameter yielding the best fit is β = 0.113 for Roberts' calculation and 0.108 for Davison's. It is found that He–pH 2 collisions are more effective than pH 2 –pH 2 in producing the J = 0 to J = 2 transition, but less effective for higher‐order transitions. Collisions of neon with pH 2 are found to be more effective at room temperature for inducing the 0 → 2 and 2 → 4 transitions than either helium or argon.

Leonard M Valley↗

NASA NTRS2017

Further Validation of a CFD Code for Calculating the Performance of Two-Stage Light Gas Guns

Earlier validations of a higher-order Godunov code for modeling the performance of two-stage light gas guns are reviewed. These validation comparisons were made between code predictions and experimental data from the NASA Ames 1.5" and 0.28" guns and covered muzzle velocities of 6.5 to 7.2 km/s. In the present report, five more series of code validation comparisons involving experimental data from the Ames 0.22" (1.28" pump tube diameter), 0.28", 0.50", 1.00" and 1.50" guns are presented. The total muzzle velocity range of the validation data presented herein is 3 to 11.3 km/s. The agreement between the experimental data and CFD results is judged to be very good. Muzzle velocities were predicted within 0.35 km/s for 74% of the cases studied with maximum differences being 0.5 km/s and for 4 out of 50 cases, 0.5 - 0.7 km/s.

Ballistic Range↗

NASA NTRS2008

Radiative Susceptibility of Cloudy Atmospheres to Droplet Number Perturbations: 1. Theoretical Analysis and Examples from MODIS

Theoretical and satellite-based assessments of the sensitivity of broadband shortwave radiative fluxes in cloudy atmospheres to small perturbations in the cloud droplet number concentration (N) of liquid water clouds under constant water conditions are performed. Two approaches to study this sensitivity are adopted: absolute increases in N, for which the radiative response is referred to as absolute cloud susceptibility, and relative increases in N or relative cloud susceptibility. Estimating the former is more challenging as it requires an assumed value for either cloud liquid water content or geometrical thickness; both susceptibilities require an assumed relationship between the droplet volume and effective radius. Expanding upon previous susceptibility studies, present radiative calculations include the effect of AN perturbations on droplet asymmetry parameter and single-scattering albedo, in addition to extinction. Absolute cloud susceptibility has a strong nonlinear dependence on the droplet effective radius as expected, while relative cloud susceptibility is primarily dependent on optical thickness. Molecular absorption and reflecting surfaces both reduce the relative contribution of the cloud to the top-of-atmosphere (TOA) flux and therefore also reduce the TOA albedo susceptibility. Transmittance susceptibilities are negative with absolute values similar to albedo susceptibility, while atmospheric absorptance susceptibilities are about an order of magnitude smaller than albedo susceptibilities and can be either positive or negative. Observation-based susceptibility calculations are derived from MODIS pixel-level retrievals of liquid water cloud optical thickness, effective radius, and cloud top temperature; two data granule examples are shown. Susceptibility quantifies the aerosol indirect effect sensitivity in a way that can be easily computed from model fields. As such, susceptibilities derived from MODIS observations provide a higher-order test of model cloud properties used for indirect effect studies. MODIS-derived global distributions of cloud susceptibility and radiative forcing calculations are presented in a companion paper.

Platnick, Steven↗

NASA NTRS2002

Towards Arbitrary Accuracy Inviscid Surface Boundary Conditions

Inviscid nonlinear surface boundary conditions are currently limited to third order accuracy in time for non-moving surfaces and actually reduce to first order in time when the surfaces move. For steady-state calculations it may be possible to achieve higher accuracy in space, but high accuracy in time is required for efficient simulation of multiscale unsteady phenomena. A surprisingly simple technique is shown here that can be used to correct the normal pressure derivatives of the flow at a surface on a Cartesian grid so that arbitrarily high order time accuracy is achieved in idealized cases. This work demonstrates that nonlinear high order time accuracy at a solid surface is possible and desirable, but it also shows that the current practice of only correcting the pressure is inadequate.

Dyson, Rodger W.↗

NASA NTRS1988

K-epsilon turbulence model assessment with reduced numerical diffusion for coaxial jets

A numerical study has been conducted to evaluate the performance of the k-epsilon turbulence model for axisymmetric, unconfined, swirling, and nonswirling coaxial jet flows. Two lower order schemes, hybrid and power-law, and two higher order schemes, flux-spline and bounded skew upwind differencing, were employed in this investigation. The predicted results indicate that the higher order numerical schemes have greater potential for future model improvements and complex flow calculations in terms of storage, accuracy, and execution time. For the nonswirling flow, computations using an algebraic stress model have also been presented.

Nikjooy, M.↗

NASA NTRS1991

Aerosol characterization in smoke plumes from a wetlands fire

Results are presented from airborne measurements of aerosol mass loading, size distribution, and elemental composition obtained in a smoke plume from a fire at a Florida wildlife refuge. It is found that there was a high concentration of carbon-containing aerosols and salt crystals in the 0.1-0.2 micron size range; the composition and morphology of the aerosols differed with size. Aerosol mass concentrations are used in conjunction with those obtained for CO2 in order to calculate ratios of aerosol and CO2. The ratios are noted to be higher for the smoldering phase of the fire than for its flaming phase.

Woods, David C.↗

NASA NTRS1998

Utilization of Low Gravity Environment for Measuring Liquid Viscosity

The method of drop coalescence is used for determining the viscosity of highly viscous undercooled liquids. Low gravity environment is necessary in order to allow for examining large volumes affording much higher accuracy for the viscosity calculations than possible for smaller volumes available under 1 - g conditions. The drop coalescence method is preferred over the drop oscillation technique since the latter method can only be applied for liquids with vanishingly small viscosities. The technique developed relies on both the highly accurate solution of the Navier-Stokes equations as well as on data from experiments conducted in near zero gravity environment. Results are presented for method validation experiments recently performed on board the NASA/KC-135 aircraft. While the numerical solution was produced using the Boundary Element Method. In these tests the viscosity of a highly viscous liquid, glycerine at room temperature, was determined using the liquid coalescence method. The results from these experiments will be discussed.

Antar, Basil N.↗

NASA NTRS1990

Turbulent supersonic convection in three dimensions

Previous numerical calculations of two-dimensional, compressible convection are extended to three dimensions, using a higher order Godunov scheme. The results show that the flow readily becomes supersonic in the upper boundary layer, where shock structures form intermittently in the vicinity of the strong downflow lanes. The convection as a whole is strongly time-dependent and evolves on a time scale comparable to the sound crossing time. The motions in the upper layers are characterized by the rapid expansion of the upward-moving fluid elements. In the interior, most of the heat is carried by a small fraction of the fluid residing in strong, highly coherent downflows. The remaining fluid is dominated by small-scale, disorganized turbulent motions.

Malagoli, Andrea↗

NASA NTRS2005

Computing Gravitational Fields of Finite-Sized Bodies

A computer program utilizes the classical theory of gravitation, implemented by means of the finite-element method, to calculate the near gravitational fields of bodies of arbitrary size, shape, and mass distribution. The program was developed for application to a spacecraft and to floating proof masses and associated equipment carried by the spacecraft for detecting gravitational waves. The program can calculate steady or time-dependent gravitational forces, moments, and gradients thereof. Bodies external to a proof mass can be moving around the proof mass and/or deformed under thermoelastic loads. An arbitrarily shaped proof mass is represented by a collection of parallelepiped elements. The gravitational force and moment acting on each parallelepiped element of a proof mass, including those attributable to the self-gravitational field of the proof mass, are computed exactly from the closed-form equation for the gravitational potential of a parallelepiped. The gravitational field of an arbitrary distribution of mass external to a proof mass can be calculated either by summing the fields of suitably many point masses or by higher-order Gauss-Legendre integration over all elements surrounding the proof mass that are part of a finite-element mesh. This computer program is compatible with more general finite-element codes, such as NASTRAN, because it is configured to read a generic input data file, containing the detailed description of the finiteelement mesh.

Quadrelli, Marco↗

NASA NTRS1959

A Numerical Method for Obtaining Monoenergetic Neutron Flux Distributions and Transmissions in Multiple-Region Slabs

This method is investigated for semi-infinite multiple-slab configurations of arbitrary width, composition, and source distribution. Isotropic scattering in the laboratory system is assumed. Isotropic scattering implies that the fraction of neutrons scattered in the i(sup th) volume element or subregion that will make their next collision in the j(sup th) volume element or subregion is the same for all collisions. These so-called "transfer probabilities" between subregions are calculated and used to obtain successive-collision densities from which the flux and transmission probabilities directly follow. For a thick slab with little or no absorption, a successive-collisions technique proves impractical because an unreasonably large number of collisions must be followed in order to obtain the flux. Here the appropriate integral equation is converted into a set of linear simultaneous algebraic equations that are solved for the average total flux in each subregion. When ordinary diffusion theory applies with satisfactory precision in a portion of the multiple-slab configuration, the problem is solved by ordinary diffusion theory, but the flux is plotted only in the region of validity. The angular distribution of neutrons entering the remaining portion is determined from the known diffusion flux and the remaining region is solved by higher order theory. Several procedures for applying the numerical method are presented and discussed. To illustrate the calculational procedure, a symmetrical slab ia vacuum is worked by the numerical, Monte Carlo, and P(sub 3) spherical harmonics methods. In addition, an unsymmetrical double-slab problem is solved by the numerical and Monte Carlo methods. The numerical approach proved faster and more accurate in these examples. Adaptation of the method to anisotropic scattering in slabs is indicated, although no example is included in this paper.

Schneider, Harold↗

NASA NTRS1995

Building a better leapfrog

In stellar dynamical computer simulations, as well as other types of simulations using particles, time step size is often held constant in order to guarantee a high degree of energy conservation. In many applications, allowing the time step size to change in time can offer a great saving in computational cost, but variable-size time steps usually imply a substantial degradation in energy conservation. We present a meta-algorithm' for choosing time steps in such a way as to guarantee time symmetry in any integration scheme, thus allowing vastly improved energy conservation for orbital calculations with variable time steps. We apply the algorithm to the familiar leapfrog scheme, and generalize to higher order integration schemes, showing how the stability properties of the fixed-step leapfrog scheme can be extended to higher order, variable-step integrators such as the Hermite method. We illustrate the remarkable properties of these time-symmetric integrators for the case of a highly eccentric elliptical Kepler orbit and discuss applications to more complex problems.

Hut, Piet↗

NASA NTRS1975

Calculation of arbitrary-order diffraction efficiencies of thick gratings with arbitrary grating shape

A method for calculating arbitrary-order diffraction efficiencies of thick, lossless transmission gratings with arbitrary periodic grating shapes has been developed. A Fourier-series representation of the grating is employed, along with a coupled-mode theory of diffraction. For illustration, numerical values of the diffraction efficiencies at the first three Bragg angles are calculated for sinusoidal, square-wave, triangular, and saw-tooth gratings. Numerical results for the same grating shapes with the same parameters are also calculated for comparison, by extending Burckhardt's numerical method for analyzing thick sinusoidal gratings. The comparison shows that the coupled-mode theory provides results with relative computational ease and results that are in agreement with calculations obtained by extending the more-rigorous Burckhardt theory to nonsinusoidal grating shapes and to higher-order Bragg angles.

Su, S. F.↗