Engineering Papers⌕ Search

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.

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

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

On the method of pseudo compressibility for numerically solving incompressible flows

Pseudo compressibility is used for numerically solving incompressible flows to achieve computational efficiency. The use of pseudo compressibility results in a system of hyperbolic-type equations of motion that introduce waves of finite speed. The interactions of the wave propagation and the vorticity spreading are analyzed. A criterion governing the dependence of the pseudo compressiblity on the Reynolds number and the characteristic length of the flow geometry is obtained that allows for a proper convergence. It is demonstrated that the solution does tend to the incompressible limit. External and internal viscous flow test problems are presented to verify the theory.

Chang, J. L. C.↗
NASA NTRS1983

Wrinkling phenomenon in structures. II - Wrinkling criteria

The wrinkling phenomenon in elasticity is mathematically defined and described as a special type of instability associated with the loss of hyperbolicity of the governing equations when some of the characteristic speeds of elastic wave propagation become imaginary. The developed criteria of wrinkling are applied to such structural elements as membranes, strings, soft shells, soft rods, and laminated material. The cause of wrinkling in these elements is due to compression in longitudinal directions that exceeds the through-the-thickness shear modulus. For elements with internal fluid flow the critical compression is reduced due to the destabilizing effect of the flow. The most effective mathematical analysis technique is a combination of tensor analysis and index notations.

Zak, M.↗
NASA NTRS1984

A comparative study of the parabolized Navier-Stokes code using various grid-generation techniques

The parabolized Navier-Stokes (PNS) equations are used to calculate the flow-field characteristics about the hypersonic research aircraft X-24C. A comparison of the results obtained using elliptic, hyperbolic, and algebraic grid generators is presented. The outer bow shock is treated as a sharp discontinuity, and the discontinuities within the shock layer are captured. Surface pressures and heat-transfer results at angles of attack of 6 deg and 20 deg, obtained using the three grid generators, are compared. The PNS equations are marched downstream over the body in both Cartesian and cylindrical base coordinate systems, and the results are compared. A robust marching procedure is demonstrated by successfully using large marching step sizes with the implicit shock fitting procedure. A correlation is found between the marching step size, Reynolds number, and the angle of attack at fixed values of smoothing and stability coefficients for the marching scheme.

Kaul, U. K.↗
NASA NTRS1983

Dynamical evolution of the Oort cometary cloud

The dynamical evolution of comets in the Oort cloud under the influence of stellar perturbations has been modeled using Monte Carlo techniques. It is shown that the cloud has been depleted over the history of the solar system. Comets are lost from the cloud by direct ejection due to close stellar encounters, diffusion of aphelia to distances beyond the sun's sphere of influence, or diffusion of perihelia into the planetary region where Jupiter and Saturn perturbations either eject them on hyperbolic trajectories or capture them to short-period orbits. The population of the cloud is estimated to be 1.0 - 1.5 x 10 to the 12th comets and the total mass is on the order of 1.9 earth masses. In addition to random passing stars, less frequent encounters with giant molecular clouds may play a significant role in randomizing the orbits of comets in the cloud and reducing the effective radius of the sun's sphere of influence.

Weissman, P. R.↗
NASA NTRS1982

A minimal residual method for transonic potential flows

For transonic flow calculations, a combination of the successive line over-relaxation (SLOR) and the preconditioned conjugate gradient (CG) method has been suggested by Wong and Hafez (1981). This paper studies the method of minimal residual (MR) which avoids a combined iteration. This method is closely related to the CG method, may be regarded as a first-order gradient method, and is applicable to symmetric and nonsymmetric matrices. The problem is formulated as a nonlinear mixed elliptic-hyperbolic partial differential equation which includes an artificial viscosity and a switching function which is zero in subsonic regions and nonzero in supersonic regions. Alternatives to the SLOR method which provide faster convergence rates are introduced. The preconditioned MR algorithm is developed, and transonic potential flows around NACA 0012 airfoil are calculated for different Mach numbers and angles of attack. Preliminary results are presented, demonstrating that the MR algorithm requires no parameter estimation and rapidly converges for subsonic flows.

Wong, Y. S.↗
NASA NTRS1982

A high-resolution numerical technique for inviscid gas-dynamic problems with weak solutions

The shock resolution of Harten's (1982) second-order explicit method for one-dimensional hyperbolic conservation laws is investigated for a two-dimensional gas-dynamic problem. The possible extension to a high resolution implicit method for both one- and two-dimensional problems is also investigated. Applications of Harten's method to the quasi-one-dimensional nozzle problem with two nozzle shapes (divergent and convergent-divergent) and the two-dimensional shock-reflection problem resulted in high shock resolution steady-state numerical solutions.

Yee, H. C.↗
NASA NTRS1984

A spectral element method for fluid dynamics - Laminar flow in a channel expansion

A spectral element method that combines the generality of the finite element method with the accuracy of spectral techniques is proposed for the numerical solution of the incompressible Navier-Stokes equations. In the spectral element discretization, the computational domain is broken into a series of elements, and the velocity in each element is represented as a high-order Lagrangian interpolant through Chebyshev collocation points. The hyperbolic piece of the governing equations is then treated with an explicit collocation scheme, while the pressure and viscous contributions are treated implicitly with a projection operator derived from a variational principle. The implementation of the technique is demonstrated on a one-dimensional inflow-outflow advection-diffusion equation, and the method is then applied to laminar two-dimensional (separated) flow in a channel expansion. Comparisons are made with experiment and previous numerical work.

Patera, A. T.↗
NASA NTRS1984

Postinstability models in elasticity

It is demonstrated that the instability caused by the failure of hyperbolicity in elasticity and associated with the problem of unpredictability in classical mechanics expresses the incompleteness of the original model of an elastic medium. The instability as well as the ill-posedness of the Cauchy problem are eliminated by reformulating the original model.

Zak, M.↗
NASA NTRS1984

Relaxation and approximate factorization methods for the unsteady full potential equation

The unsteady form of the full potential equation is solved in conservation form, using implicit methods based on approximate factorization and relaxation schemes. A local time linearization for density is introduced to enable solution to the equation in terms of phi, the velocity potential. A novel flux-biasing technique is applied to generate proper forms of the artificial viscosity, to treat hyperbolic regions with shocks and sonic lines present. The wake is properly modeled by accounting not only for jumps in phi, but also for jumps in higher derivatives of phi obtained from requirements of density continuity. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. Results are presented for flows over airfoils, cylinders, and spheres. Comparisons are made with available Euler and full potential results.

Shankar, V.↗
NASA NTRS1984

Transformation of two and three-dimensional regions by elliptic systems

The research during this period continued to expand the class of numerical algorithms that can be accurately and efficiently implemented on overlapping grids. Whereas previous calculations have been used to solve elliptic equations and to find the steady-state solution of parabolic equations, the present work is aimed towards developing time-accurate solution techniques for parabolic and hyperbolic equations. The primary difficulty here is in the correct treatment of the interior boundary nodes that must be updated at each iteration. The implementation of explicit methods is straightforward. However, the common practice of lagging these values when using an implicit methods leads to inconsistencies in the difference equation. One way to avoid this problem is to alternately calculate with an implicit and an explicit method on each subgrid. With this procedure, the explicit method generates boundary values at the next time level which are then used by the implicit step. It can be shown that when a backward implicit method is combined with a forward explicit method, the composite method is second order accurate and unconditionally stable for linear problems. A second area in which progress can be reported is in the distribution of grid points on curves and surfaces.

Mastin, E. W.↗
NASA NTRS1984

Tidal disruption and the origin of the Moon

The dynamic problem of the tidal disruption of a rocky planetismal was solved by a direct integration of the fully three-dimensional, nonlinear equations of motion. The hypothesis that any object that passes within the Roche limit is disrupted was disproven. A time dependent solution was performed numerically, treating the planetismal as a fluid with a Murnaghan equation of state in the solid regions and zero pressure otherwise. Calculations show that a rocky body which passes by the Earth on a parabolic orbit with a perigee within the Roche limit is not tidally disrupted. Objects on hyperbolic orbits would experience even less tidal disruption. The results herein do not apply to bodies with very low viscosity. It is shown, however, that tidal disruption can be ruled out as a mechanism for reducing planetismal masses. Mechanisms for forming the Moon which rely upon tidal disruption are unlikely to be correct.

Mizuno, H.↗
NASA NTRS1984

Mechanical models of close approaches and collisions of large protoplanets

In models of lunar origin by great impact, attention is usually paid to the hydrodynamic expansion resulting from the great amount of thermal energy. However, the source of this disruption is inevitably significant in a close approach between large bodies, it is to be expected that significant departures from simple hyperbolic orbits would occur even before impact. These departures could arise from mechanical effects, and hence purely mechanical models are worth pursuing. The most interesting results obtained for approach offsets are small multiples of the planet radius and approach velocities of a few kilometers/second. In an interaction between Mars and Earth sized protoplanets, most of the material ends in collision, but a few percent end in elliptic orbits and a few percent escape. Another model considered is an offset collision, arising from a wide range of approach velocities and offsets.

Kaula, W. M.↗
NASA NTRS1984

Lunar origin from impact on the Earth

All theories of lunar origin involve events or processes which seemingly have low efficiencies or low probabilities or both. An impact-triggered fission lunar origin is presented. If the impact ejecta (a mixture of target and projectile) leave the impact site ballistically and are subsequently acted upon only by the gravity field of a spherical Earth, then the ejecta either reimpacts the Earth or escapes on a hyperbolic trajectory. Hence the need for a second burn. Three possible resolutions are considered: pressure gradient acceleration, non-central gravity, and viscous spreading.

Stevenson, D. J.↗
NASA NTRS1985

The eigenvalues of the pseudospectral Fourier approximation to the operator sin (2x) d/dx

It is shown that the eigenvalues Z sub i of the pseudospectral Fourier approximation to the operator sin(2x) curly d/curly dx satisfy (R sub e) (Z sub i) = + or - 1 or (R sub e)(Z sub I) = 0. Whereas this does not prove stability for the Fourier method, applied to the hyperbolic equation U sub t = sin (2x)(U sub x) - pi x pi; it indicates that the growth in time of the numerical solution is essentially the same as that of the solution to the differential equation.

Tal-Ezer, H.↗
NASA NTRS1985

Recovering pointwise values of discontinuous data within spectral accuracy

The pointwise values of a function, f(x), can be accurately recovered either from its spectral or pseudospectral approximations, so that the accuracy solely depends on the local smoothness of f in the neighborhood of the point x. Most notably, given the equidistant function grid values, its intermediate point values are recovered within spectral accuracy, despite the possible presence of discontinuities scattered in the domain. (Recall that the usual spectral convergence rate decelerates otherwise to first order, throughout). To this end, a highly oscillatory smoothing kernel is employed in contrast to the more standard positive unit-mass mollifiers. In particular, post-processing of a stable Fourier method applied to hyperbolic equations with discontinuous data, recovers the exact solution modulo a spectrally small error. Numerical examples are presented.

Gottlieb, D.↗
NASA NTRS1985

Uniformly high-order accurate non-oscillatory schemes, 1

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws was begun. These schemes share many desirable properties with total variation diminishing schemes (TVD), but TVD schemes have at most first order accuracy, in the sense of truncation error, at extreme of the solution. A uniformly second order approximation was constucted, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell.

Harten, A.↗
NASA NTRS1985

Shock capturing

Recent developments which have improved the understanding of how finite difference methods resolve discontinuous solutions to hyperbolic partial differential equations are discussed. As a result of this understanding improved shock capturing methods are currently being developed and tested. Some of these methods are described and numerical results are presented showing their performance on problems containing shocks in one and two dimensions. A conservative difference scheme is defined. Conservation implies that, except in very special circumstances, shocks must be spread over at least two grid intervals. These two interval shocks are actually attained in one dimension if the shock is steady and an upwind scheme is used. By analyzing this case, the reason for this excellent shock resolution can be determined. This result is used to provide a mechanism for improving the resolution of two dimensional steady shocks. Unfortunately, this same analysis shows that these results cannot be extended to shocks which move relative to the computing grid. Total variation diminishing (TVD) finite difference schemes and flux limiters are introduced to deal with money shocks and contact discontinuities.

Davis, S. F.↗