SEARCH · Engineering Papers

Results for “hyperbolic”

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 955 records · Page 53

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Faster Algorithm For Computation Of Incompressible Flow

Improved algorithm yields faster numerical solutions of Navier-Stokes equations of steady or unsteady three-dimensional flow of incompressible fluid. In artificial-compressibility method, unsteady flow treated as incompressible in advancing from one time step to next, but at each time step (or in steady state), fluid treated as having variable compressibility enabling propagation of flow field, and subiterations performed in increments of pseudotime until effects of compressibility subside. Directly couples pressure and velocity fields at same time step and converts elliptic incompressible Navier-Stokes equations to hyperbolic form more amenable to numerical integration.

Rogers, S. E.↗

NASA NTRS1991

Overview Of Methods For Computation Of Shocks

NASA technical memorandum provides systematic overview of class of conservative finite-difference shock-capturing numerical-integration methods for solution of hyperbolic conservation laws. Unified and generalized formulation for one class of methods presented.

Yee, H. C.↗

NASA NTRS1990

Spectral methods for time dependent problems

Spectral approximations are reviewed for time dependent problems. Some basic ingredients from the spectral Fourier and Chebyshev approximations theory are discussed. A brief survey was made of hyperbolic and parabolic time dependent problems which are dealt with by both the energy method and the related Fourier analysis. The ideas presented above are combined in the study of accuracy stability and convergence of the spectral Fourier approximation to time dependent problems.

Tadmor, Eitan↗

NASA NTRS1990

Spurious frequencies as a result of numerical boundary treatments

The stability theory for finite difference Initial Boundary-Value approximations to systems of hyperbolic partial differential equations states that the exclusion of eigenvalues and generalized eigenvalues is a sufficient condition for stability. The theory, however, does not discuss the nature of numerical approximations in the presence of such eigenvalues. In fact, as was shown previously, for the problem of vortex shedding by a 2-D cylinder in subsonic flow, stating boundary conditions in terms of the primitive (non-characteristic) variables may lead to such eigenvalues, causing perturbations that decay slowly in space and remain periodic time. Characteristic formulation of the boundary conditions avoided this problem. A more systematic study of the behavior of the (linearized) one-dimensional gas dynamic equations under various sets of oscillation-inducing legal boundary conditions is reported.

Abarbanel, Saul↗

NASA NTRS1990

An investigation of design optimization using a 2-D viscous flow code with multigrid

Computational fluid dynamics (CFD) codes have advanced to the point where they are effective analytical tools for solving flow fields around complex geometries. There is also a need for their use as a design tool to find optimum aerodynamic shapes. In the area of design, however, a difficulty arises due to the large amount of computer resources required by these codes. It is desired to streamline the design process so that a large number of design options and constraints can be investigated without overloading the system. There are several techniques which have been proposed to help streamline the design process. The feasibility of one of these techniques is investigated. The technique under consideration is the interaction of the geometry change with the flow calculation. The problem of finding the value of camber which maximizes the ratio of lift over drag for a particular airfoil is considered. In order to test out this technique, a particular optimization problem was tried. A NACA 0012 airfoil was considered at free stream Mach number of 0.5 with a zero angle of attack. Camber was added to the mean line of the airfoil. The goal was to find the value of camber for which the ratio of lift over drag is a maximum. The flow code used was FLOMGE which is a two dimensional viscous flow solver which uses multigrid to speed up convergence. A hyperbolic grid generation program was used to construct the grid for each value of camber.

Doria, Michael L.↗

NASA NTRS1990

Wide acceptance angle, high concentration ratio, optical collector

The invention is directed to an optical collector requiring a wide acceptance angle, and a high concentration ratio. The invention is particularly adapted for use in solar collectors of cassegrain design. The optical collector system includes a parabolic circular concave primary mirror and a hyperbolic circular convex secondary mirror. The primary mirror includes a circular hole located at its center wherein a solar collector is located. The mirrored surface of the secondary mirror has three distinct zones: a center circle, an on-axis annulus, and an off-axis section. The parabolic shape of the primary mirror is chosen so that the primary mirror reflects light entering the system on-axis onto the on-axis annulus. A substantial amount of light entering the system off-axis is reflected by the primary mirror onto either the off-axis section or onto the center circle. Subsequently, the off-axis sections reflect the off-axis light toward the solar collector. Thus, off-axis light is captured which would otherwise be lost to the system. The novelty of the system appears to lie in the configuration of the primary mirror which focuses off-axis light onto an annular portion of the secondary mirror to enable capture thereof. This feature results in wide acceptance angle and a high concentration ratio, and also compensates for the effects of non-specular reflection, and enables a cassegrain configuration to be used where such characteristics are required.

Kruer, Mark Arthur↗

NASA NTRS1990

Asymptotic-induced numerical methods for conservation laws

Asymptotic-induced methods are presented for the numerical solution of hyperbolic conservation laws with or without viscosity. The methods consist of multiple stages. The first stage is to obtain a first approximation by using a first-order method, such as the Godunov scheme. Subsequent stages of the method involve solving internal-layer problems identified by using techniques derived via asymptotics. Finally, a residual correction increases the accuracy of the scheme. The method is derived and justified with singular perturbation techniques.

Garbey, Marc↗

NASA NTRS1990

Dynamics of groups around interacting double ellipticals: Measuring dark matter haloes

Binary galaxies, as binary stars, are important to measure masses, as suggested by Page (1952). Because three orbit parameters are measurable for galaxies at one instant of time, severe uncertainties remain in the orbit and mass determinations. These uncertainties can partly be overcome by statistical studies of selected samples and/or n-body simulations. Close double galaxies (and isolated galaxies) could also be useful to estimate dynamical masses if we can find test particles around them. Interacting elliptical pairs or dumb-bell galaxies are found with a large range, between 0-1200 km s(exp -1), of relative radial velocities. Standard 2-body orbit calculations, highly uncertain due to projection factors, suggest for the largest velocity differences very large galaxy masses, if the systems are bound and stationary. However, recent n-body simulations model these binaries as galaxies captured from hyperbolic orbits, requiring masses of order a few times 10(exp 11) solar maximum (Borne et al. 1988), but producing systems that are short lived. A different picture appears when we study observationally the dynamical mass of interacting double ellipticals using faint satellite galaxies. These satellites contribute little luminosity and, presumably, little mass to the system. The authors present results of two such groups, basically forming systems of test particles, around the dumb-bells NGC 4782/3 and IC 5049. They also briefly discuss the satellite group around the central dumb-bell in the cluster Sersic 40/6. Apparently, they detect large quantities of dark matter in the vicinity of these dumb-bell galaxies, because the system masses of approx. 4.5 times 10(exp 13) solar mass and 8 times 10(exp 13) solar mass for NGC 4782/3 and IC 5049, respectively, are quite high. Likewise, the mass of the Sersic 40/6 inner core is 7 times 10(exp 13) solar mass. The possibility that a common massive dark matter halo increases the merging times of these types of galaxies is suggested. Here, the authors assume H sub o = 100 km s (exp -1) Mpc(exp -1).

Quintana, H.↗

NASA NTRS1989

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗

NASA NTRS1991

Computation of supersonic and low subsonic cascade flows using an explicit Navier-Stokes technique and the kappa-epsilon turbulence model

A fully explicit two-dimensional flow solver, based on a four-stage Runge-Kutta scheme, was developed and utilized to predict two-dimensional viscous flow through turbomachinery cascades for which experimental data is available. The formulation is applied to the density averaged Navier-Stokes equations. Several features of the technique improve the ability of the code to predict high Reynolds number flows on highly stretched grids. These include a low Reynolds number compressible form of the k-epsilon turbulence model, anisotropic scaling of artificial dissipation terms and locally varying timestep evaluation based on hyperbolic and parabolic stability considerations. Comparisons between computation and experiment are presented for both a supersonic and a low-subsonic compressor cascade. These results indicate that the code is capable of predicting steady two-dimensional viscous cascade flows over a wide range of Mach numbers in reasonable computational times.

Kunz, R. F.↗

NASA NTRS1991

The direct simulation of high-speed mixing-layers without and with chemical heat release

A direct numerical simulation of high speed reacting and non-reacting flows for H2-air systems is presented. The calculations are made for a convective Mach number of 0.38 with hyperbolic tangent initial profile and finite rate chemical reactions. A higher-order numerical method is used in time accurate mode to time advance the solution to a statistical steady state. About 600 time slices of all the variables are then stored for statistical analysis. It is shown that most of the problems of high-speed combustion with air are characterized by relatively weak heat release. The present study shows that: (1) the convective speed is reduced by heat release by about 10 percent at this convective Mach number M(sub c) = 0.38; (2) the variation of the mean and rms fluctuation of temperature can be explained on the basis of temperature fluctuation between the flame temperature and the ambient; (3) the growth rate with heat release is reduced by 7 percent; and (4) the entrainment is reduced by 25 percent with heat release. These differences are small in comparison with incompressible flow dynamics, and are argued to be due to the reduced importance of heat release in comparison with the large enthalpy gradients resulting from the large-scale vortex dynamics. It is finally suggested that the problems of reduced mixing in high-speed flows are not severely complicated by heat release.

Sekar, B.↗

NASA NTRS1991

Nozzle fabrication technique

A block of electrically conductive material which is to be formed into a body with internal and/or external surfaces that approximate hyperboloids of one sheet is placed so that its axis is set at a predetermined skew angle with relation to a traveling EDM electrode wire. The electrode wire is then moved into cutting proximity of the body wire. Thereafter, by revolving the body about its own axis, the external and/or internal surfaces of the body will be cut into an approximate hyperbolic surface of revolution depending upon whether the body is positioned with the cutting wire outside of the body or in a previously formed longitudinal passage in the body. As an alternative technique, elongated channels can also be cut into the wall of the body by successively orienting the body to a selected number of angular positions, with the electrode wire being either outside of the body or in a previously formed passage in the body. At each of these angular positions, the electrode wire is moved orthogonally with respect to the axis of the wire, while both the body axis skew angle and the rotational position about that axis is controlled by cutting a channel or groove in the body to relieve stresses in the body material or to convey a coolant fluid.

Wells, Dennis L.↗

NASA NTRS1990

Particle orbits in two-dimensional equilibrium models for the magnetotail

Assuming that there exist an equilibrium state for the magnetotail, particle orbits are investigated in two-dimensional kinetic equilibrium models for the magnetotail. Particle orbits in the equilibrium field are compared with those calculated earlier with one-dimensional models, where the main component of the magnetic field (Bx) was approximated as either a hyperbolic tangent or a linear function of z with the normal field (Bz) assumed to be a constant. It was found that the particle orbits calculated with the two types of models are significantly different, mainly due to the neglect of the variation of Bx with x in the one-dimensional fields.

Karimabadi, H.↗

NASA NTRS1990

Long-time behaviour of absorbing boundary conditions

A new class of computational far-field boundary conditions for hyperbolic partial differential equations was recently introduced by the authors. These boundary conditions combine properties of absorbing conditions for transient solutions and properties of far-field conditions for steady states. This paper analyses the properties of the wave equation coupled with these new boundary conditions: well-posedness, dissipativity and convergence in time.

Engquist, B.↗

NASA NTRS1990

A two-step iterative method for evolving nonlinear acoustic systems to a steady-state

A new approach for evolving two-dimensional nonlinear acoustic systems with flow to a steady state is presented. The approach is a two-step iterative method which is tested on a benchmark acoustic problem for which an exact analytical solution is available. Results are also calculated for a nonlinear acoustic problem for which an exact analytical solution is not known. Results indicate that the two-step method represents a powerful, efficient, and stable method for evolving two-dimensional acoustic systems to a steady state, and that the method is applicable to any number of spatial dimensions and to other hyperbolic systems. It is noted that for the benchmark problem only a single iteration on the method is required when the transient and steady-state field are of the same order of magnitude; however, four iterations are required when the steady-state field is several orders of magnitude smaller than the transient field. This method requires six iterations before achieving a steady state for the nonlinear test problem.

Watson, Willie R.↗

NASA NTRS1991

Numerical prediction of the unsteady flowfield around the F-18 aircraft at large incidence

This paper describes a numerical method capable of solving the steady and unsteady viscous flow around complete aircraft configurations at high angles of attack. This method is used to simulate the external flow around the F-18 aircraft, including deflected control surfaces. The current technique employs a generalized overset zonal grid scheme to decompose the computational space around the aircraft. The grid around various components of the aircraft are created numerically using a three-dimensional hyperbolic grid generation procedure. The Reynolds-averaged Navier-Stokes equations are integrated using a time-accurate, implicit procedure. Results for the turbulent flow around the F-18 aircraft at 30 degrees angle of attack show the details of the flowfield structure, including the unsteadiness created by the vortex burst and the resulting fluctuating airloads exerted on the vertical tail. The computed results agree fairly well with flight data for surface pressure, surface flow pattern, vortex burst location, and the dominant frequency for tail load fluctuations.

Rizk, Yehia M.↗

NASA NTRS1991

Finite-volume application of high-order ENO schemes to two-dimensional boundary-value problems

Finite-volume applications of high-order accurate ENO schemes to two-dimensional boundary-value problems are studied. These schemes achieve high-order spatial accuracy, in smooth regions, by a piecewise polynomial approximation of the solution from cell averages. In addition, this spatial operation involves an adaptive stencil algorithm in order to avoid the oscillatory behavior that is associated with interpolation across steep gradients. High-order TVD Runge-Kutta methods are employed for time integration, thus making these schemes best suited for unsteady problems. Fifth- and sixth-order accurate applications are validated through a grid refinement study involving the solutions of scalar hyperbolic equations. A previously proposed extension for the Euler equations of gas dynamics is tested, including its application to solutions of boundary-value problems involving solid walls and curvilinear coordinates.

Casper, Jay↗

NASA NTRS1991

Fractal dimension in nonhyperbolic chaotic scattering

In chaotic scattering there is a Cantor set of input-variable values of zero Lebesgue measure (i.e., zero total length) on which the scattering function is singular. For cases where the dynamics leading to chaotic scattering is nonhyperbolic (e.g., there are Kolmogorov-Arnol'd-Moser tori), the nature of this singular set is fundamentally different from that in the hyperbolic case. In particular, for the nonhyperbolic case, although the singular set has zero total length, strong evidence is presented to show that its fractal dimension is 1.

Lau, Yun-Tung↗