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 883 records · Page 49

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Construction of explicit and implicit symmetric TVD schemes and their applications

A one-parameter family of explicit and implicit total variation diminishing (TVD) schemes is developed which permits incorporation of an expanded group of slope and flux limiters. The numerical technique is intended for use in calculations which include a time-differencing scheme and an optional Lax-Wendroff scheme. Methods of extending the TVD models to nonlinear scalar equations and systems of hyperbolic conservation equations are described. Sample results are presented from calculations of shocked flows around NACA 0012 and NACA 0018 airfoils.

Yee, H. C.↗

NASA NTRS1986

On the search for trajectories for a cometary nucleus sample return mission

The steps in the procedure for the generation of single trajectories for a comet nucleus sample return mission described by Stuhlinger et al. (1985) are discussed. The following assumptions are made for such a mission: (1) the spacecraft is equipped with an ion propulsion system supplied by solar cells, and (2) the comet has a low inclination and a period of about six years. The calculation of a specific trajectory incorporates the following main parameters: (1) the spacecraft mass at the beginning of the mission and the hyperbolic excess velocity, (2) the mass as a function of the injection energy, and (3) the available thrust level. The limitation of a supply of solar electric power to the region within 3 AU of the sun and the placement of permanent in-situ experiments at the comet are discussed. It is concluded that single trajectories for a nucleus sample return mission can be generated. A further step would be an integration of single solutions into classes of trajectories.

Schobert, D.↗

NASA NTRS1986

Entropy conservative finite element schemes

The question of entropy stability for discrete approximations to hyperbolic systems of conservation laws is studied. The amount of numerical viscosity present in such schemes is quantified and related to their entropy stability by means of comparison. To this end, two main ingredients are used: entropy variables and the construction of certain entropy conservative schemes in terms of piecewise-linear finite element approximations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only if, they contain more numerical viscosity than the abovementioned entropy conservation ones.

Tadmor, E.↗

NASA NTRS1987

The line integral approach to radarclinometry

Advances over previous work in radarclinometry (Wildey, 1984 and 1986) are proposed, including incorporation of a more exact treatment of the dependence on terrain orientation of the solid angle of integration of the microwave specific energy comprising individual pixel power, and the integration of the theory in the cross-ground-range direction. The local curvature hypothesis, making possible formulation of the theory in terms of a line integral, is expanded to include the possibilities of local cylindricity, local biaxial ellipsoidal hyperbolicity, and least-squares local sphericity. It is suggested that the banding noted in topographic maps which are produced on the basis of variation in the proportioning between curvature assumptions may be due to the lack of radiometric calibration.

Wildey, Robert L.↗

NASA NTRS1987

Extension of a semi-implicit shock-capturing algorithm for 3-D fully coupled, chemically reacting flows in generalized coordinates

A semiimplicit high-resolution shock-capturing method for multidimensional systems of hyperbolic conservation laws with stiff source terms has been developed by Yee and Shinn (1987). The goal of this work is to extend this method to solve the three-dimensional fully coupled Navier-Stokes equations for a hypersonic chemically reacting flow in generalized coordinates. In this formulation, the global continuity equation was replaced by all the species continuity equations. The shock-capturing technique is a second-order-accurate, symmetric total-variation-diminishing method which accounts fully and directly for the coupling among the fluid and all the species. To verify the current approach, it was implemented into an existing computer code which contained the MacCormack method. Test results for a five-species reacting flow are shown to be oscillation-free around the shock, and the time spent per iteration only doubles when compared to the result using classical way of supplying numerical dissipation. The extra computation is more than justified by the elimination of spurious oscillation and nonlinear instability associated with the classical shock-capturing schemes in computing hypersonic reacting flows.

Shinn, J. L.↗

NASA NTRS1987

Uniformly high order accurate essentially non-oscillatory schemes. III

In the present continuation of development and analysis efforts towards essentially nonoscillatory shock-capturing approximations of hyperbolic conservation laws, the Godunov (1959) scheme and its MUSCL second-order-accurate extension are generalized to an arbitrary order of accuracy by means of a hierarchy of uniformly high-order-accurate schemes. The design method, which involves the solution's essentially nonoscillatory piecewise polynomial reconstruction on the basis of its cell averages, followed by time-evolution through an approximate solution of the resulting initial value problem and its averaging over each cell, employs an adaptive stencil of grid points to yield schemes that are highly nonlinear.

Harten, Ami↗

NASA NTRS1987

Steady state computations for wave propagation problems

The behavior of difference approximations of hyperbolic partial differential equations as time t goes to infinity is studied. The rate of convergence to steady state is analyzed theoretically and experimentally for the advection equation and the linearized Euler equations. The choice of difference formulas and boundary conditions strongly influences the rate of convergence in practical steady-state calculations. In particular it is shown that upwind difference methods and characteristic boundary conditions have very attractive convergence properties.

Engquist, Bjorn↗

NASA NTRS1987

The numerical viscosity of entropy stable schemes for systems of conservation laws. I

Discrete approximations to hyperbolic systems of conservation laws are studied. The amount of numerical viscosity present in such schemes is quantified and related to their entropy stability by means of comparison.To this end, conservative schemes which are also entropy-conservative are constructed. These entropy-conservative schemes enjoy second-order accuracy; moreover, they can be interpreted as piecewise-linear finite-element methods, and hence can be formulated on various mesh configurations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only they contain more viscosity than that present in the above-mentioned entropy-conservative ones.

Tadmor, Eitan↗

NASA NTRS1987

Fronts propagating with curvature dependent speed: Algorithms based on Hamilton-Jacobi formulations

New numerical algorithms are devised (PSC algorithms) for following fronts propagating with curvature-dependent speed. The speed may be an arbitrary function of curvature, and the front can also be passively advected by an underlying flow. These algorithms approximate the equations of motion, which resemble Hamilton-Jacobi equations with parabolic right-hand-sides, by using techniques from the hyperbolic conservation laws. Non-oscillatory schemes of various orders of accuracy are used to solve the equations, providing methods that accurately capture the formation of sharp gradients and cusps in the moving fronts. The algorithms handle topological merging and breaking naturally, work in any number of space dimensions, and do not require that the moving surface be written as a function. The methods can be used also for more general Hamilton-Jacobi-type problems. The algorithms are demonstrated by computing the solution to a variety of surface motion problems.

Osher, Stanley↗

NASA NTRS1987

Analytical investigation of the dynamics of tethered constellations in Earth orbit (phase 2)

Simulation of two short distance crawling maneuvers of the elevator both with and without environmental perturbations acting upon the system is discussed. These simulation runs were performed in order to provide results useful for the interpretation of the data from the tests, on the ground, of a scaled down engineering model of the elevator. In these simulation runs the elevator crawls along the tether in accordance to the developed mirror image motion control law (MIMCL). Results from the simulation of the 4 km long maneuver run were compared to those obtained by adopting the modified hyperbolic tangent control law (MHTCL). A preprocessor was developed for setting up the initial conditions of a tethered system with L platforms, M longitudinal dampers, and N lumped masses (platforms plus tether beads). A short test run of the 4-platform system with 3 longitudinal dampers and 10 lumper mass without any perturbation acting upon the system is illustrated. In support of the Tether Applications Working Group, SAO is preparing a catalog of tether simulations, has prepared a set of simulator test cases, obtained results from SKYHOOK, and solicited results from simulators at other institutions, and prepared a paper on a specific analytic solution.

Lorenzini, Enrico C.↗

NASA NTRS1987

SHARP simulation of discontinuities in highly convective steady flow

For steady multidimesional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

NASA NTRS1988

Total Variation Diminishing (TVD) schemes of uniform accuracy

Explicit second-order accurate finite-difference schemes for the approximation of hyperbolic conservation laws are presented. These schemes are nonlinear even for the constant coefficient case. They are based on first-order upwind schemes. Their accuracy is enhanced by locally replacing the first-order one-sided differences with either second-order one-sided differences or central differences or a blend thereof. The appropriate local difference stencils are selected such that they give TVD schemes of uniform second-order accuracy in the scalar, or linear systems, case. Like conventional TVD schemes, the new schemes avoid a Gibbs phenomenon at discontinuities of the solution, but they do not switch back to first-order accuracy, in the sense of truncation error, at extrema of the solution. The performance of the new schemes is demonstrated in several numerical tests.

Hartwich, PETER-M.↗

NASA NTRS1988

Accurate boundary conditions for exterior problems in gas dynamics

The numerical solution of exterior problems is typically accomplished by introducing an artificial, far field boundary and solving the equations on a truncated domain. For hyperbolic systems, boundary conditions at this boundary are often derived by imposing a principle of no reflection. However, waves with spherical symmetry in gas dynamics satisfy equations where incoming and outgoing Riemann variables are coupled. This suggests that natural reflections may be important. A reflecting boundary condition is proposed based on an asymptotic solution of the far field equations. Nonlinear energy estimates are obtained for the truncated problem and numerical experiments presented to validate the theory.

Hagstrom, Thomas↗

NASA NTRS1988

Nozzle fabrication technique

This invention relates to techniques for fabricating hour glass throat or convergent divergent nozzle shapes, and more particularly to new and improved techniques for forming rocket nozzles from electrically conductive material and forming cooling channels in the wall thereof. The concept of positioning a block of electrically conductive material so that its axis is set at a predetermined skew angle with relation to a travelling electron discharge machine electrode and thereafter revolving the body about its own axis to generate a hyperbolic surface of revolution, either internal or external is novel. The method will generate a rocket nozzle which may be provided with cooling channels using the same control and positioning system. The configuration of the cooling channels so produced are unique and novel. Also the method is adaptable to nonmetallic material using analogous cutting tools, such as, water jet, laser, abrasive wire and hot wire.

Wells, Dennis L.↗

NASA NTRS1988

The solution of the Elrod algorithm for a dynamically loaded journal bearing using multigrid techniques

A numerical solution to a theoretical model of vapor cavitation in a dynamically loaded journal bearing is developed utilizing a multigrid iteration technique. The method is compared with a noniterative approach in terms of computational time and accuracy. The computational model is based on the Elrod algorithm, a control volume approach to the Reynolds equation which mimics the Jakobsson-Floberg and Olsson cavitation theory. Besides accounting for a moving cavitation boundary and conservation of mass at the boundary, it also conserves mass within the cavitated region via a smeared mass or striated flow extending to both surfaces in the film gap. The mixed nature of the equations (parabolic in the full film zone and hyperbolic in the cavitated zone) coupled with the dynamic aspects of the problem create interesting difficulties for the present solution approach. Emphasis is placed on the methods found to eliminate solution instabilities. Excellent results are obtained for both accuracy and reduction of computational time.

Woods, Claudia M.↗

NASA NTRS1988

Non-oscillatory spectral Fourier methods for shock wave calculations

A non-oscillatory spectral Fourier method is presented for the solution of hyperbolic partial differential equations. The method is based on adding a nonsmooth function to the trigonometric polynomials which are the usual basis functions for the Fourier method. The high accuracy away from the shock is enhanced by using filters. Numerical results confirm that no oscillations develop in the solution. Also, the accuracy of the spectral solution of the inviscid Burgers equation is shown to be higher than a fixed order.

Cai, Wei↗

NASA NTRS1988

Calculation of non-stationary aerodynamic forces in the near sonic range

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.

Teipel, I.↗

NASA NTRS1987

Vectorizable implicit algorithms for the flux-difference split, three-dimensional Navier-Stokes equations

The computational efficiency of four vectorizable implicit algorithms is assessed when applied to calculate steady-state solutions to the three-dimensional, incompressible Navier-Stokes equations in general coordinates. Two of these algorithms are characterized as hybrid schemes; that is, they combine some approximate factorization in two coordinate directions with relaxation in the remaining spatial direction. The other two algorithms utilize an approximate factorization approach which yields two-factor algorithms for three-dimensional systems. All four algorithms are implemented in identical high-resolution upwind schemes for the flux-difference split Navier-Stokes equations. These highly nonlinear schemes are obtained by extending an implicit Total Variation Diminishing (TVD) scheme recently developed for linear one-dimensional systems of hyperbolic conservation laws to the three-dimensional Navier-Stokes equations. The computation of vortical flow over a sharp-edged, thin delta wing has been chosen as a common numerical test case. The convergence of the algorithms is discussed and the accuracy of the computed flow-field results is assessed. The validity of the present results are demonstrated by a comparison with experimental data.

Hartwich, P. M.↗