SEARCH · Search NASA
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.
Experiments done with the aim of improving the accuracy and efficiency of explicit boundary algorithms for quasilinear hyperbolic systems of partial differential equations in matrix form are described. Existing methods are criticized and a new approach to the problem is presented.
The development of a mathematical method for calculating nonstationary supersonic flow in the near-sonic range is described. A perturbation formula is derived based on the exact stationary values; it is applicable to the equation for potential. The problem can thus be divided into stationary and nonstationary fields. The pressure distribution in an oscillating profile is determined, based on hyperbolic differential equations. It is shown that there are important corollaries concerning the application of linear theory. With suitable extrapolations, linear theory can be used up to about Mach 0.8. Linear theory is not applicable, however, when determining the moment coefficients; for this case, a special technique is described.
In these lecture notes we describe the construction, analysis, and application of ENO (Essentially Non-Oscillatory) and WENO (Weighted Essentially Non-Oscillatory) schemes for hyperbolic conservation laws and related Hamilton- Jacobi equations. ENO and WENO schemes are high order accurate finite difference schemes designed for problems with piecewise smooth solutions containing discontinuities. The key idea lies at the approximation level, where a nonlinear adaptive procedure is used to automatically choose the locally smoothest stencil, hence avoiding crossing discontinuities in the interpolation procedure as much as possible. ENO and WENO schemes have been quite successful in applications, especially for problems containing both shocks and complicated smooth solution structures, such as compressible turbulence simulations and aeroacoustics. These lecture notes are basically self-contained. It is our hope that with these notes and with the help of the quoted references, the reader can understand the algorithms and code them up for applications.
The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts to give a unified and generalized formulation of a class of high-resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic nonlinear, spatially second and third-order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.
The paper is concerned with the application of the Murman and Cole (1971) relaxation scheme to steady, inviscid transonic flow problems in two and three dimensions. This scheme, which automatically accounts for weak shock waves, uses separate difference operators in elliptic and hyperbolic regions. The details of the scheme are described in terms of the original small disturbance formulation of Murman and Cole. In particular, Murman's recent (1973) introduction of fully conservative difference operators to obtain the correct shock jumps is examined. The extension to treating the exact isentropic equation is then covered with special attention given to Jameson's (to appear) rotated difference scheme for supersonic flow regions. The bulk of the discussion is related to two-dimensional procedures, and some comparisons with experiment are made, with emphasis on the effects of viscosity and wind-tunnel walls. Application of the Murman-Cole scheme is then discussed for small disturbances in three dimensions.
A new method is described for the determination of optimal spacecraft trajectories in an inverse-square field using finite, fixed thrust. The method employs a recently developed optimization technique which uses a piecewise polynomial representation for the state and controls, and collocation, thus converting the optimal control problem into a nonlinear programming problem, which is solved numerically. This technique has been modified to provide efficient handling of those portions of the trajectory which can be determined analytically, i.e., the coast arcs. Among the problems that have been solved using this method are optimal rendezvous and transfer (including multirevolution cases) and optimal multiburn orbit insertion from hyperbolic approach.
A novel approach for the solution of inviscid flow problems for subsonic compressible flows is described. The approach is based on canonical forms of the equations, in which subsystems governed by hyperbolic operators are separated from those governed by elliptic ones. The discretizations used as well as the iterative techniques for the different subsystems are inherently different. Hyperbolic parts, which describe, in general, propagation phenomena, are discretized using upwind schemes and are solved by marching techniques. Elliptic parts, which are directionally unbiased, are discretized using h-elliptic central discretizations, and are solved by pointwise relaxations together with coarse grid acceleration. The resulting discretization schemes introduce artificial viscosity only for the hyperbolic parts of the system; thus a smaller total artificial viscosity is used, while the multigrid solvers used are much more efficient. Solutions of the subsonic compressible Euler equations are achieved at the same efficiency as the full potential equation.
Lambert’s problem subject to a continuous acceleration is solved using the family of generalized logarithmic spirals. Thanks to the existence of two first integrals related to the energy and angular momentum surprising analogies with the Keplerian case are found. A minimum-energy spiral transfer exists. Increasing the value of the constant of the generalized energy yields pairs of conjugate spiral trajectories. The properties of such spirals are strongly connected with the properties of conjugate Keplerian orbits. When the generalized constant of the energy reaches a critical value the two solutions degenerate into a pair of parabolic spirals, one of which connects the two vectors through infinity. From that point the spiral transfers become hyperbolic. Generalized logarithmic spirals admit closed-form solutions to all the required magnitudes including the time of flight, providing a deep insight into the dynamics of the problem. In addition, the maximum acceleration along the transfer is found analytically so the solutions that violate the design constraints on the maximum thrust acceleration can be rejected without any further computations. When the time of flight is fixed there is still a degree of freedom in the solution, related to a control parameter. Resonant transfers appear naturally thanks to the symmetry properties of the generalized logarithmic spirals. The problem of designing a low-thrust transfer between two bodies can be reduced to solving the corresponding spiral Lambert’s problem. In order to show the versatility of the method it is applied to the design of an asteroid tour and to explore launch opportunities to Mars.\
New nonlinear partial differential equations containing terrain curvature and its rate of change are derived that describe the flow of an atmospheric density current. Unlike the classical hydraulic-type equations for density currents, the new equations are valid for two-dimensional, gradually varied flow over highly curved terrain, hence suitable for computing unsteady (or steady) flows over arbitrary mountain/valley profiles. The model assumes the atmosphere above the density current exerts a known arbitrary variable pressure upon the unknown interface. Later this is specialized to the varying hydrostatic pressure of the atmosphere above. The new equations yield the variable velocity distribution, the interface position, and the pressure distribution that contains a centrifugal component, often significantly larger than its hydrostatic component. These partial differential equations are hyperbolic, and the characteristic equations and characteristic directions are derived. Using these to form a characteristic mesh, a hypothetical unsteady curved-flow problem is calculated, not based upon observed data, merely as an example to illustrate the simplicity of their application to unsteady flows over mountains.
The objective of this paper is to investigate the effects of the numerical dissipation and the resolution of the solution on coarser grids for multigrid with the Euler equation approximations. The convergence is accomplished by multi-stage explicit time-stepping to steady state accelerated by FAS multigrid. A theoretical investigation is carried out for linear hyperbolic equations in one and two dimensions. The spectra reveals that for stability and hence robustness of spatial discretizations with a small amount of numerical dissipation the grid transfer operators have to be accurate enough and the smoother of low temporal accuracy. Numerical results give grid independent convergence in one dimension. For two-dimensional problems with a small amount of numerical dissipation, however, only a few grid levels contribute to an increased speed of convergence. This is explained by the small numerical dissipation leading to dispersion. Increasing the mesh density and hence making the problem over resolved increases the number of mesh levels contributing to an increased speed of convergence. If the steady state equations are elliptic, all grid levels contribute to the convergence regardless of the mesh density.
An analysis is presented of supercircular entry into a planet's atmosphere giving particular attention to the corridor through which spacecraft must be guided in order to accomplish various maneuvers. A dimensionless parameter based on conditions at the conic perigee altitude is introduced for characterizing supercircular entries and conveniently pre-scribing corridor widths associated with elliptic, parabolic, or hyperbolic approach trajectories. The analysis applies to vehicles of arbitrary weight, shape, and size. Illustrative calculations are made for Venus, Earth, Mars, Jupiter, and Titan. For nonlifting vehicles having fixed aerodynamic coefficients, curves are presented of dimensionless parameters from which can be calculated the maximum deceleration, maximum rate of laminar convective heating, and total laminar heat absorbed during single-pass entry at velocities up to twice circular velocity. For lifting vehicles, curves are presented of the maximum deceleration and overshoot boundary of an entry corridor; equations are presented for estimating laminar aerodynamic heating from the maximum deceleration. It is shown that the corridor width is independent of vehicle weight, dimensions, and drag coefficient, provided these are the same at the overshoot boundary as at undershoot. The corridors of certain planets can be broadened markedly by the application of aerodynamic lift; for example, the 10-earth-g corridor width for single- pass, nonlifting, parabolic entry is increased from 0 miles for Jupiter, 7 for Earth, and 8 for Venus, to 52, 51, and 52 miles, respectively, by employing a lift-drag ratio of 1. The use of aerodynamic lift does not increase appreciably the corridors of Mars and Titan. All corridor widths decrease rapidly as the entry velocity is increased. Terminal guidance requirements on accuracy of velocity and flight path angle for successfully entering various corridors are compared with analogous requirements for putting a satellite into orbit, for hitting the moon from the earth, and for achieving ICBM accuracy. Consideration is given to the terminal guidance problem involved in using a planet's atmosphere-rather than rocket fuel-to effect orbital transfers from heliocentric to planetocentric motion, thereby converting a hyperbolic approach trajectory to an elliptic orbit about the target planet. This fuel saving maneuver appears technologically feasible for certain planetary voyages, and implies the possibility of achieving a large reduction in required Earth lift-off weight of chemical propulsion systems.
We propose a one-parameter family of low-dissipation centred numerical schemes for solving hyperbolic equations in conservative or non-conservative form, using finite volume or discontinuous Galerkin finite element methods. The new schemes spring out from the multi-dimensional FORCE method and are determined by a single parameter α ≥ 1. Given an increasing sequence of real numbers 1 ≤ α 1 < α 2 < . . . < α K , there corresponds a sequence of numerical schemes with stability restriction associated to a decreasing sequence of Courant numbers 1 > c 1 > c 2 > . . . > c K > 0 and a decreasing sequence of corresponding numerical viscosity functions d 1 > d 2 > . . . > d K . For a given Courant number c k ≤ 1 there is a real number α k > 0 and a corresponding stable scheme with minimal numerical viscosity d k . The proposed schemes suit very well the family of high-order discontinuous Galerkin finite element methods and the recently proposed class of ADER-TR schemes, whose orders of accuracy define decreasing sequences of Courant numbers, as the order of accuracy increases. The centred methods of this paper are stable in 2D and 3D in the frame of simultaneous updating formulae, unlike other centred methods, such as 1D FORCE, which are not. Furthermore, the schemes are highly accurate for slowly-moving waves, which is precisely the kind of waves that traditional centred methods smear disastrously. Additional features of the proposed schemes include ease of implementation and applicability to any hyperbolic system either in conservative or non-conservative form. Here, the proposed schemes are analysed and computationally assessed through a suite of test problems for a linear model system, for the Euler equations in one and two space dimensions, and for the Baer-Nunziato equations for compressible two-phase flow.
Pyranda Is fortran 90 code which solves compact finite difference equations on a three dimensional fortran arrays, domain decomposed using the Message Passing Interface (MPI). The fortran source solves the linear system of equations which forms a pentadtagona1 matrix. Derivatives are used to construct hyperbolic partlal differential equation (PDE) operators, namely gradient, divergence, Laplacian, and curl operators. 5. Solution Method: The library provides a python Interface which facilitates problem setup and description. However, the core solver Is written in the fortran 90 language. The solver couples MPI communication and linear solves together so that a direct (non-iterative) solution to the linear system of equations can be computed with minimal communication. Domain decomposition is limited to a uniform data size but provides perfect load balancing. 6. Computer for which the code Is written:
Analytical procedures for extracting piecewise smooth solutions of hyperbolic systems from raw oscillatory data obtained by pseudospectral methods are developed. The validity of the approach is demonstrated for the case of linear problems with constant coefficients, and plausibility arguments are presented which indicate its applicability to nonlinear operators when the steady state has been achieved. Numerical results for the development of an oblique shock when a wedge is inserted at zero angle of attack into a uniform supersonic flow of an ideal gas (the time-dependent two-dimensional Euler equations discretized in space by the pseudospectral Chebyshev method) are presented in tables and graphs.
The Galerkin criterion within a finite element Weighted Residuals formulation is employed to establish an implicit solution algorithm for an initial-value partial differential equation. Numerical solutions of a transient parabolic and a hyperbolic equation, obtained using linear, quadratic and two cubic finite element basis functions, are employed to quantize accuracy and confirm and refine theoretical convergence rate estimates. The linear basis algorithm for the hyperbolic equation displays excellent accuracy on a coarse computational grid and a high-order convergence rate with discretization refinement. Good accuracy and a strong convergence rate in surface flux are determined for a nonhomogeneous Neumann boundary constraint applied to a parabolic equation. The results amply demonstrate the impact of the nondiagonal finite element initial-value matrix structure on solution accuracy and/or convergence rate.
Extensions of a second order high resolution explicit method for the numerical computation of weak solutions of one dimensonal hyperbolic conservation laws are discussed. The main objectives were (1) to examine the shock resoluton of Harten's method for a two dimensional shock reflection problem, (2) to study the use of a high resolution scheme as a post-processor to an approximate steady state solution, and (3) to construct an implicit in the delta-form using Harten's scheme for the explicit operator and a simplified iteration matrix for the implicit operator.
A constant coefficient hyperbolic system in one space variable, with zero initial data is discussed. Dissipative boundary conditions are imposed at the two points x = + or - 1. This problem is discretized by a spectral approximation in space. Sufficient conditions under which the spectral numerical solution is stable are demonstrated - moreover, these conditions have to be checked only for scalar equations. The stability theorems take the form of explicit bounds for the norm of the solution in terms of the boundary data. The dependence of these bounds on N, the number of points in the domain (or equivalently the degree of the polynomials involved), is investigated for a class of standard spectral methods, including Chebyshev and Legendre collocations.
Smoothed particle hydrodynamics (SPH) discretization techniques are generalized to develop a method, smoothed particle interpolation (SPI), for solving initial value problems of systems of a nonhydrodynamical nature. Under this approach, SPH is viewed as strickly as interpolation scheme and, as such, suitable for solving general hyperbolic and parabolic equations. The SPI method is tested on (1) the wave equation with inhomogenous sound speed and (2) Burger's equation. The efficiency of SPI is studied by comparing SPI solutions to those obtained with standard finite-difference methods. It is shown that the power of SPI arises when the smoothing particles are free to move.