Search NASA⌕ Search

SEARCH · Search NASA

Results for “hyperbolicity”

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 721 records · Page 40

TIGGERC: Turbomachinery interactive grid generator energy distributor and restart code

A two dimensional multi-block grid generator was developed for a new design and analysis system for studying multi-blade row turbomachinery problems with an axisymmetric viscous/inviscid 'average passage' through flow code. TIGGERC is a mouse driven, fully interactive grid generation program which can be used to modify boundary coordinates and grid packing. TIGGERC generates grids using a hyperbolic tangent or algebraic distribution of grid points on the block boundaries and the interior points of each block grid are distributed using a transfinite interpolation approach. TIGGERC generates a blocked axisymmetric H grid, C grid, I grid, or O grid for studying turbomachinery flow problems. TIGGERC was developed for operation on small high speed graphic workstations.

Miller, David P.↗

Grid generation research at OSU

In the last two years, effort was concentrated on: (1) surface modeling; (2) surface grid generation; and (3) 3-D flow space grid generation. The surface modeling shares the same objectives as the surface modeling in computer aided design (CAD), so software available in CAD can in principle be used for solid modeling. Unfortunately, however, the CAD software cannot be easily used in practice for grid generation purposes, because they are not designed to provide appropriate data base for grid generation. Therefore, we started developing a generalized surface modeling software from scratch, that provides the data base for the surface grid generation. Generating surface grid is an important step in generating a 3-D space for flow space. To generate a surface grid on a given surface representation, we developed a unique algorithm that works on any non-smooth surfaces. Once the surface grid is generated, a 3-D space can be generated. For this purpose, we also developed a new algorithm, which is a hybrid of the hyperbolic and the elliptic grid generation methods. With this hybrid method, orthogonality of the grid near the solid boundary can be easily achieved without introducing empirical fudge factors. Work to develop 2-D and 3-D grids for turbomachinery blade geometries was performed, and as an extension of this research we are planning to develop an adaptive grid procedure with an interactive grid environment.

Nakamura, S.↗

Numerical methods for the simulation of complex multi-body flows with applications for the integrated Space Shuttle vehicle

The following papers are presented: (1) numerical methods for the simulation of complex multi-body flows with applications for the Integrated Space Shuttle vehicle; (2) a generalized scheme for 3-D hyperbolic grid generation; (3) collar grids for intersecting geometric components within the Chimera overlapped grid scheme; and (4) application of the Chimera overlapped grid scheme to simulation of Space Shuttle ascent flows.

Chan, William M.↗

Direct simulation of high-speed mixing layers

A computational study of a nonreacting high-speed mixing layer is performed. A higher order algorithm with sufficient grid points is used to resolve all relevant scales. In all cases, a temporal free-stream disturbance is introduced. The resulting flow is time-sampled to generate a statistical cross section of the flow properties. The studies are conducted at two convective Mach numbers, three free-stream turbulence intensities, three Reynolds numbers, and two types of initial profiles-hyperbolic tangent (tanh) and boundary layer. The boundary-layer profile leads to more realistic predictions of the transition processes. The predicted transition Reynolds number of 0.18 x 10(exp 6) compares well with experimental data. Normalized vortex spacings for the boundary-layer case are about 3.5 and compare favorably with the 1.5 to 2.5 found in experimental measurements. The tanh profile produces spacings of about 10. The growth rate of the layer is shown to be moderately affected by the initial disturbance field, but comparison with experimental data shows moderate agreement. For the boundary-layer case, it is shown that noise at the Strouhal number of 0.007 is selectively amplified and shows little Reynolds number dependence.

Mukunda, H. S.↗

A numerical resolution study of high order essentially non-oscillatory schemes applied to incompressible flow

High order essentially non-oscillatory (ENO) schemes, originally designed for compressible flow and in general for hyperbolic conservation laws, are applied to incompressible Euler and Navier-Stokes equations with periodic boundary conditions. The projection to divergence-free velocity fields is achieved by fourth order central differences through Fast Fourier Transforms (FFT) and a mild high-order filtering. The objective of this work is to assess the resolution of ENO schemes for large scale features of the flow when a coarse grid is used and small scale features of the flow, such as shears and roll-ups, are not fully resolved. It is found that high-order ENO schemes remain stable under such situations and quantities related to large-scale features, such as the total circulation around the roll-up region, are adequately resolved.

Weinan, E.↗

Arrival and departure impulsive Delta V determination for precessing Mars parking orbits

An attempt is made to develop a method for realistic estimation of the initial LEO mass. The method takes into account the actual geometry between the inbound and outbound hyperbolic asymptotes and the parking orbit, along with precession effects caused by the oblateness of Mars, in calculating the arrival and departure Delta V values. Three mission scenarios alternative to the arbitrarily assumed two tangential periapsis burns are described: a tangential periapsis arrival and an in-plane departure; an in-plane arrival and in-plane departure; and a tangential periapsis arrival and a 3D departure. Results obtained by the method under consideration compared well with a trajectory integration code, where the differences in the initial LEO orbit mass were within one percent, for all three cases. The method is considered to be an ideal tool for preliminary mission design, since it reduces the analysis computation time with minimal loss in accuracy.

Desai, Prasun N.↗

A staggered mesh finite difference scheme for the computation of compressible flows

A simple high resolution finite difference technique is presented to approximate weak solutions to hyperbolic systems of conservation laws. The method does not rely on Riemann problem solvers and is therefore easy to extend to a wide variety of problems. The overall performance (resolution and CPU requirements) is competitive, with other state-of-the-art techniques offering sharp nonoscillatory shocks and contacts. Theoretical results confirm the reliability of the approach for linear systems and nonlinear scalar equations.

Sanders, Richard↗

On steady field-aligned double-adiabatic flow

The paper summarizes the four general double-adiabatic invariants governing ideal MHD flow in a narrow tube and renders them in a form suitable for the calculation of flow in tubes of nonconstant cross-sectional area. Results for converging-diverging flux tubes are developed and compared to single-adiabatic results which are identical to those for ordinary gasdynamic flow. It is found that in certain parameter regimes the governing equation for the vector potential is elliptic and can be reduced to Laplace's equation by simple stretching of one coordinate. In other regimes the equation is hyperbolic, indicating the presence of wave patterns. The relevance of these results to flow over magnetopause undulations, including flux transfer events, is discussed.

Sonnerup, B. U. O.↗

Determining Mars parking orbits which ensure tangential periapsis burns at arrival and departure

A method is presented which finds Mars parking orbits which allow tangential periapsis burns at both arrival and departure. This method accounts for the actual geometry at both arrival and departure between the hyperbolic asymptotes and the orbital plane, along with the precession effects caused by the oblateness of Mars. Thus, realistic Delta-V values (and hence initial low-earth orbit masses) are obtained for these orbits. The results obtained from the present method compare very well with a trajectory integration program while only requiring CPU times of about one minute. Therefore, due to the computational efficiency and accuracy, the present method would be an ideal tool to use in preliminary mission design, since it provides the opportunity to incorporate realistic Mars parking orbits effects.

Desai, Prasun N.↗

The Space Transportation System summer environment on launch pad

This paper describes a 2D flow and thermalanalysis to determine the solar effect on the Space Shuttle launch components subsequent to the external tank (ET) loading operation in extremely hot conditions. An existing CFD code Parabolic Hyperbolic or Elliptical Numerical Integration Code Series was used in the study. The analysis was done for a 2D slice between planes perpendicular to the longitudinal axis of the STS and passing through the lower portions of the Redesigned Solid Rocket Motors (RSRMs), the ET, and the wing of the Orbiter. The results are presented as local and average values of the heat transfer coefficient, and the Nusselt number, and the surface temperature around the RSRMs and the ET. Solar heating effects increased the surface temperatures of the RSRMs by 9-11 F. Higher prelaunch surface temperatures measured on the east and west RSRMs (in the inboard region between the RSRMs and the ET) during 19 most recent launches of the STS are correlated as a function of the ambient temperature.

Ahmad, R. A.↗

Ulysses dust measurements near Jupiter

Submicrometer- to micrometer-sized particles were recorded by the Ulysses dust detector within 40 days of the Jupiter flyby. Nine impacts were recorded within 50 Jupiter radii with most of them recorded after closest approach. Three of these impacts are consistent with particles on prograde orbits around Jupiter and the rest are believed to have resulted from gravitationally focused interplanetary dust. From the ratio of the impact rate before the Jupiter flyby to the impact rate after the Jupiter flyby it is concluded that interplanetary dust particles at the distance of Jupiter move on mostly retrograde orbits. On 10 March 1992, Ulysses passed through an intense dust stream. The dust detector recorded 126 impacts within 26 hours. The stream particles were moving on highly inclined and apparently hyperbolic orbits with perihelion distances of greater than 5 astronomical units. Interplanetary dust is lost rather quickly from the solar system through collisions and other mechanisms and must be almost continuously replenished to maintain observed abundances. Dust flux measurements, therefore, give evidence of the recent rates of production from sources such as comets, asteroids, and moons, as well as the possible presence of interstellar grains.

Gruen, Eberhard↗

Development and verification of methods for predicting flow rates through leaks in valves and couplings

This is the final report of a research effort which addresses the title problem. The report discusses two broad models of flows, which represent the following extreme cases: (1) inertia-dominated flow, where friction is relatively insignificant; and (2) friction-dominated flow where inertia is insignificant. In class (2), the leak channel might consist of the gap between a scratch in a plastic seal and a polished metal plate against which the seal is pressed. Here, the cross section of the leak channel is modeled as a flat bottomed crescent. A publication generated under the present grant period presents an exact solution of the equations of fully-developed laminar pipe flow of a liquid in the case of a crescent beneath a hyperbolic arc. A Master's thesis project supported by the present grant presents the corresponding solution beneath a circular arc. A second publication reviews the flow of a gas through the same channel, which may be analyzed by a standard one-dimensional model (Fanno flow) for an engineering approximation. Finally, the report discusses the design and progress in the fabrication of a leak-test cell, in which one may measure the flow of fluid through a controlled flaw in a seal. The aim of such measurements is to furnish data for comparison with the predictions of the theory.

Russell, John M.↗

Comet nongravitational forces and meteoritic impacts

We have considered those comets whose original orbits have been determined to be hyperbolic when only planetary perturbations are accounted for. It is found that formally unbound incident trajectories correlate most confidently with orbits that have small perihelion distances and move in a retrograde sense relative to planetary motion. Arguments are presented that these results are not due to measurement error or to selection effects. We conclude that the phenomenon is attributable to enhanced volatility leading to abnormally large nongravitational forces. Since the effect is absent in the prograde small-perihelia population, increased insolation is not the sole explanation. It is suggested that the significance of the retrograde correlation is connected with a larger energy of relative motion between retrograde comets and a population of prograde ecliptic meteoroids which impact the comet mantle exposing the underlying volatiles. The subsequent enhanced outgassing is the cause of the larger nongravitational forces.

Matese, John J.↗

The geocentric particulate distribution: Cometary, asteroidal, or space debris?

Definition of the Low Earth Orbit (LEO) particulate environment has been refined considerably with the analysis of data from NASA's Long Duration Exposure Facility (LDEF). Measurements of the impact rates from particulates ranging from sub-micron to millimetres in dimension and, especially, information on their directionality has permitted new scrunity of the sources of the particulates. Modelling of the dynamics of both bound (Earth orbital) and unbound (hyperbolic interplanetary) particulates intercepting LDEF's faces leads to the conclusion that the source is dominantly interplanetary for particle dimensions of greater than some 5 microns diameter; however the anisotropy below this dimension demands lower velocities and is compatible with an orbital component. Characteristics of the LDEF interplanetary component are compatible with familiar meteoroid sources and deep space measurements. Understanding of the orbital component which exceeds the interplanetary flux by a factor of 4 is less clear; although the very small particulates in orbit have been associated with space debris (Lawrance and Brownlee, 1986) this data conflicts with other measurements (McDonnell, Carey and Dixon, 1984) at the same epoch. By analysis of trajectories approaching the Earth and its atmosphere, we have shown that a significant contribution could be captured by aerocapture, i.e., atmospheric drag, from either asteroidal or cometary sources; such enhancement is unlikely however to provide the temporal and spatial fluctuations observed by the LDEF Interplanetary Dust Experiment (Mullholland et al. 1992). A further new mechanism is also examined, that of aerofragmentation capture, where an atmospheric grazing trajectory, which would not normally lead to capture, leads to fragmentation by thermal or mechanical shock; the microparticulates thus created can be injected in large numbers, but only into short-lifetime orbits. The concentration in one particular orbit plane, could explain the temporal fluctuations seen on LDEF; space debris could also explain the phenomenon.

Mcdonnell, J. A. M.↗

Research in computational fluid dynamics and analysis of algorithms

Recently, higher-order compact schemes have seen increasing use in the DNS (Direct Numerical Simulations) of the Navier-Stokes equations. Although they do not have the spatial resolution of spectral methods, they offer significant increases in accuracy over conventional second order methods. They can be used on any smooth grid, and do not have an overly restrictive CFL dependence as compared with the O(N(exp -2)) CFL dependence observed in Chebyshev spectral methods on finite domains. In addition, they are generally more robust and less costly than spectral methods. The issue of the relative cost of higher-order schemes (accuracy weighted against physical and numerical cost) is a far more complex issue, depending ultimately on what features of the solution are sought and how accurately they must be resolved. In any event, the further development of the underlying stability theory of these schemes is important. The approach of devising suitable boundary clusters and then testing them with various stability techniques (such as finding the norm) is entirely the wrong approach when dealing with high-order methods. Very seldom are high-order boundary closures stable, making them difficult to isolate. An alternative approach is to begin with a norm which satisfies all the stability criteria for the hyperbolic system, and look for the boundary closure forms which will match the norm exactly. This method was used recently by Strand to isolate stable boundary closure schemes for the explicit central fourth- and sixth-order schemes. The norm used was an energy norm mimicking the norm for the differential equations. Further research should be devoted to BC for high order schemes in order to make sure that the results obtained are reliable. The compact fourth order and sixth order finite difference scheme had been incorporated into a code to simulate flow past circular cylinders. This code will serve as a verification of the full spectral codes. A detailed stability analysis by Carpenter (from the fluid Mechanics Division) and Gottlieb gave analytic conditions for stability as well as asymptotic stability. This had been incorporated in the code in form of stable boundary conditions. Effects of the cylinder rotations had been studied. The results differ from the known theoretical results. We are in the middle of analyzing the results. A detailed analysis of the effects of the heating of the cylinder on the shedding frequency had been studied using the above schemes. It has been found that the shedding frequency decreases when the wire was heated. Experimental work is being carried out to affirm this result.

Gottlieb, David↗

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.↗

On essentially non-oscillatory schemes on unstructured meshes: Analysis and implementation

A few years ago, the class of Essentially Non-Oscillatory Schemes for the numerical simulation of hyperbolic equations and systems was constructed. Since then, some extensions have been made to multidimensional simulations of compressible flows, mainly in the context of very regular structured meshes. In this paper, we first recall and improve the results of an earlier paper about non-oscillatory reconstruction on unstructured meshes, emphasizing the effective calculation of the reconstruction. Then we describe a class of numerical schemes on unstructured meshes and give some applications for its third order version. This demonstrates that a higher order of accuracy is indeed obtained, even on very irregular meshes.

Abgrall, R.↗

High-order cyclo-difference techniques: An alternative to finite differences

The summation-by-parts energy norm is used to establish a new class of high-order finite-difference techniques referred to here as 'cyclo-difference' techniques. These techniques are constructed cyclically from stable subelements, and require no numerical boundary conditions; when coupled with the simultaneous approximation term (SAT) boundary treatment, they are time asymptotically stable for an arbitrary hyperbolic system. These techniques are similar to spectral element techniques and are ideally suited for parallel implementation, but do not require special collocation points or orthogonal basis functions. The principal focus is on methods of sixth-order formal accuracy or less; however, these methods could be extended in principle to any arbitrary order of accuracy.

Carpenter, Mark H.↗