SEARCH · Engineering Papers

Results for “conservative”

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 NTRS1986

Divergence study of a high-aspect ratio, forward-swept wing

An experimental wind-tunnel study to determine the divergence characteristics of a high-aspect ratio, forward-swept wing has been conducted in the NASA Langley Research Center (LaRC) Transonic Dynamics Tunnel (TDT). The rectangular wing used for this study had a panel aspect ratio of 9.16 (lambda = 0 deg.) and the sweep angle could be set at lambda = 0 deg., -15 deg., -30 deg., -45 deg., or -60 deg. A rectangular wing tip shape was tested at each of these sweep angles. In addition, a tip shape parallel to the freestream flow was tested for a wing sweep angle of lambda = -45 deg. The root of the wing was cantilever mounted to the wall of the wind tunnel. Divergence conditions were measured at M = 0.4 for each sweep angle and tip configuration tested. Subcritical response techniques were used to extrapolate to the divergence conditions during the wind-tunnel test. The primary objective of this test was to obtain data which could be used to verify for this configuration the divergence prediction capability of an aeroelastic analysis code. Subsonic lifting surface theory (kernel function) aerodynamics are utilized by this particular code. The analytical predictions of divergence were found to be significantly conservative at all forward sweep angles. At lambda = -45 deg., the analysis was 14 percent conservative. The effect of the two tip shapes on the divergence dynamic pressure was predicted accurately by the analysis. The divergence condition for the tip shape parallel to the flow occurred at a dynamic pressure 14 percent higher than the divergence condition with a rectangular tip shape.

Cole, S. R.↗

NASA NTRS1986

Multigrid methods in structural mechanics

Although the application of multigrid methods to the equations of elasticity has been suggested, few such applications have been reported in the literature. In the present work, multigrid techniques are applied to the finite element analysis of a simply supported Bernoulli-Euler beam, and various aspects of the multigrid algorithm are studied and explained in detail. In this study, six grid levels were used to model half the beam. With linear prolongation and sequential ordering, the multigrid algorithm yielded results which were of machine accuracy with work equivalent to 200 standard Gauss-Seidel iterations on the fine grid. Also with linear prolongation and sequential ordering, the V(1,n) cycle with n greater than 2 yielded better convergence rates than the V(n,1) cycle. The restriction and prolongation operators were derived based on energy principles. Conserving energy during the inter-grid transfers required that the prolongation operator be the transpose of the restriction operator, and led to improved convergence rates. With energy-conserving prolongation and sequential ordering, the multigrid algorithm yielded results of machine accuracy with a work equivalent to 45 Gauss-Seidel iterations on the fine grid. The red-black ordering of relaxations yielded solutions of machine accuracy in a single V(1,1) cycle, which required work equivalent to about 4 iterations on the finest grid level.

Raju, I. S.↗

NASA NTRS1986

A minimum entropy principle in the gas dynamics equations

Let u(x bar,t) be a weak solution of the Euler equations, governing the inviscid polytropic gas dynamics; in addition, u(x bar, t) is assumed to respect the usual entropy conditions connected with the conservative Euler equations. We show that such entropy solutions of the gas dynamics equations satisfy a minimum entropy principle, namely, that the spatial minimum of their specific entropy, (Ess inf s(u(x,t)))/x, is an increasing function of time. This principle equally applies to discrete approximations of the Euler equations such as the Godunov-type and Lax-Friedrichs schemes. Our derivation of this minimum principle makes use of the fact that there is a family of generalized entrophy functions connected with the conservative Euler equations.

Tadmor, E.↗

NASA NTRS1985

The effects of strain heating in lithospheric stretching models

The deformation by stretching of a continental type lithosphere has been formulated so that the problem can be solved by a continuum mechanical approach. The deformation, stress state, and temperature distribution are constrained to satisfy the physical laws of conservation of mass, energy, momentum, and an experimentally defined rheological response. The conservation of energy equation including a term of strain energy dissipation is given. The continental lithosphere is assumed to have the rheology of an isotropic, incompressible, nonlinear viscous, two layered solid.

Stanton, M.↗

NASA NTRS1986

Fracture analysis of radial scientific instrument module registration fittings of the space telescope

Various pieces of the registration fittings for the Radial Scientific instrument (SI) module of the Space Telescope were examined from a fracture mechanics point of view and deemed to be fail-safe or else have had maximum allowable flaw sizes specified for them. The results of these analyses are summarized in tabular form. In many instances the applied stress levels were so low that the threshold stress intensity factor range was never reached. In most of the others the allowable flaw sizes were large enough to be detected by visual inspection. However, for some parts, such as the flexures connecting the aluminum cover to the ball retainer in the fitting at point A, the flaw sizes were rather small. Eddy current tests are capable of detecting flaws of this size (0.022 inches x 0.1 inches), so for those which have been so tested these small flaws should represent no danger of going undetected. In every instance approximations were made to err on the conservative side. These were pointed out in the discussions of the analyses for each fitting. One conservative approximation that was not mentioned, however, is the fact that retardation was not included in the crack propagation computations.

Springfield, C. W., Jr.↗

NASA NTRS1986

Embedded mesh solution of the 2-D Euler equations - Evaluation of interface formulations

Solution of the steady 2-D Euler equations using mesh embedding, or local grid refinement, with a cell-centered finite volume scheme is investigated. Embedded regions which are topologically similar to the global grid are considered. An isoenergetic model for the governing equations is used in Jameson's finite volume multistage scheme with modifications to the boundary conditions and smoothing. A detailed study of the embedding interface flux and smoothing formulations is conducted. Taylor expansion analysis reveals that local second order spatial accuracy is not possible if a conservative interface flux formulation is used. The analysis also gives constraints for local first order accuracy. An energy stability analysis indicates that downwind weighting of interface fluxes causes local instabilities. Analysis shows that conservative interface smoothing formulations must have a locally convective component, but that correct interface formulations allow globally dissipative smoothing. Embedded mesh solutions obtained with this scheme are presented for a transonic airfoil. They show that if embedding interfaces are close to the shocks, then small modifications in the interface location can have large effects on converge and solution accuracy.

Allmaras, S. R.↗

NASA NTRS1985

Foundations for the numerical solution of the Euler equations

The Navier-Stokes equations represent an extremely good model of the physical phenomena encountered in most aeronautical problems. However, the computational resource needed to solve the Navier-Stokes equations are so large that even with today's supercomputers, it is necessary to make use of simpler models. A large number of external aerodynamic problems can be accurately described by a simpler model. This model consists of an outer inviscid flow plus a boundary-layer thickness correction for the vehicle shape. The outer inviscid model may be represented by the potential equation or by the Euler equation. The present paper provides the foundations for the numerical solution of the Euler equations. The governing equations are considered, taking into account conservation laws, the medium, the differential form of the conservation laws, generalized solutions, shock-fitting, and characteristics. Attention is also given to initial and boundary conditions, existence and uniqueness, and rotational phenomena.

Salas, M. D.↗

NASA NTRS1985

AOTV bluff body flow - Relaxation algorithm

Preliminary investigations have led to the consideration of two classes of energy efficient aeroassisted orbital transfer vehicle (AOTV) concepts, with high and low L/D. The AOTV represent reentry space vehicles which are to fly hypersonically in the continuum upper reaches of the atmosphere. Problems arise in connection with convective and nonequilibrium radiative heating. As the possibilities for obtaining needed information by experimentation are limited, reliance must be placed on computational techniques to simulate testing of design concepts. A new relaxation algorithm for the compressible Navier-Stokes equations has been developed in the context of recent work with a conservative flux difference eigenvector split implicit upwind scheme CSCM (conservative supracharacteristics method). Attention is given to the CSCM implicit method, and early results for axisymmetric bluff body flow including complete forebody and base flow for a free flight range model drag brake AOTV.

Lombard, C. K.↗

NASA NTRS1986

Transonic potential flow in hyperbolic nozzles

The full potential equation for the classical problem of transonic flow through a hyperbolic nozzle (with or without a shock wave) is solved in conservation form using the finite volume method of Jameson and Caughey (1977). Either a firstor a second-order numerical viscosity is added in the direction of the flow, explicitly, in conservation form. A multigrid alternating direction implicit method is used to solve the difference equations, and the results obtained are compared with analytical and numerical results from previous researches.

Park, M.↗

NASA NTRS1987

Service-Life Extension of Explosive Escape Devices

Chemical and functional tests yield conservative service-life estimates. Approach to extension of service lives of explosive devices in aircraft escape system developed, supported by testing of representative candidate devices to evaluate quantitatively effects of service, age, and degradation, and to enable responsible, conservative service-life determinations. Five types of explosive components evaluated: rigid and flexible explosive transfer lines; one-way transfers; flexible, linear-shaped charges; and initiation-handles. Extension of service in realistic manner provides both cost savings and increased system reliability.

Bement, L. J.↗

NASA NTRS1987

Isotope Exchange in Oxide Catalyst

Replacement technique maintains level of CO2/18 in closed-cycle CO2 lasers. High-energy, pulsed CO2 lasers using rare chemical isotopes must be operated in closed cycles to conserve gas. Rare isotopes operated in closed cycles to conserve gas. Rare isotopes as CO2/18 used for improved transmission of laser beam in atmosphere. To maintain laser power, CO2 must be regenerated, and O2 concentration kept below few tenths of percent. Conditions achieved by recombining CO and O2.

Hess, Robert V.↗

NASA NTRS1986

Transformation of two and three-dimensional regions by elliptic systems

Efforts in transferring computational work from the LRC computer to the IRIS Graphics Workstation at MSU are reported and the computation of a conservative solution of a simple hyperbolic equation on an overlapping grid is discussed. Several conclusions concerning computations on overlapping grids are apparent. Problems only occur when there is a major difference in grid spacing on the individual component grids. In the case of hyperbolic equations, it is necessary that both interpolation and extrapolation be applied at the grid boundaries. When interpolated values are used at outflow boundary points, excessive oscillations in the numerical solution may be the result. The same conclusions would be valid for more complicated systems of hyperbolic equations such as the Euler equations for inviscid flow. Some of the solution values would be extrapolated at the overlap boundary, the exact number depending on the number of characteristics pointing out of the overlap region. It is also possible that similar boundary conditions may be needed for some parabolic equations such as high Reynolds number viscous flow equations. Efforts were also expended on the development of three-dimensional conservative interpolation procedures. Finally, the investigation of grid smoothing procedures were initiated during this reporting period. It was decided that the first grid smoothing algorithms will be based on the concepts of variational grid generation.

Mastin, C. W.↗

NASA NTRS1986

Transport of absolute angular momentum in quasi-axisymmetric equatorial jet streams

It is well known that prograde equatorial jet stresses cannot occur in an axisymmetric inviscid fluid, owing to the constraints of local angular momentum conservation. For a viscous fluid, the constraints of mass conservation prevent the formation of any local maximum of absolute angular momentum (m) without a means of transferring m against its gradient (delta m) in the meridional plane. The circumstances under which m can be diffused up-gradient by normal molecular viscosity are derived, and illustrated with reference to numerical simulations of axisymmetric flows in a cylindrical annulus. Viscosity is shown to act so as to tend to expel m from the interior outwards from the rotation axis. Such an effect can produce local super-rotation even in a mechanically isolated fluid. The tendency of viscosity to result in the expulsion of m is shown to be analogous in certain respects to a vorticity-mixing hypothesis for the effects of non-axisymmetric eddies of the zonally-averaged flow. It is shown how the advective and diffusive transport of m by non-axisymmetric eddies can be represented by the Transformed Eulerian Mean meridional circulation and the Eliassen-Palm (EP) flux of Andrews and McIntyre respectively, in the zonal mean. Constraints on the form and direction of the EP flux in an advective/diffusive flow for such eddies are derived, by analogy with similar constraints on the diffusive flux of m due to viscosity.

Read, P. L.↗

NASA NTRS1985

Nonlinear potential analysis techniques for supersonic aerodynamic design

A numerical method based on the conservation form of the full potential equation has been applied to the problem of three-dimensional supersonic flows with embedded subsonic regions. The governing equation is cast in a nonorthogonal coordinate system, and the theory of characteristics is used to accurately monitor the type-dependent flow field. A conservative switching scheme is employed to transition from the supersonic marching procedure to a subsonic relaxation algorithm and vice versa. The newly developed computer program can handle arbitrary geometries with fuselage, canard, wing, flow through nacelle, vertical tail and wake components at combined angles of attack and sideslip. Results are obtained for a variety of configurations that include a Langley advanced fighter concept with fuselage centerline nacelle, Rockwell's Advanced Tactical Fighter (ATF) with wing mounted nacelles, and the Shuttle Orbiter configuration. Comparisons with available experiments were good.

Shankar, V.↗

NASA NTRS1986

On orbital allotments for geostationary satellites

The following satellite synthesis problem is addressed: communication satellites are to be allotted positions on the geostationary arc so that interference does not exceed a given acceptable level by enforcing conservative pairwise satellite separation. A desired location is specified for each satellite, and the objective is to minimize the sum of the deviations between the satellites' prescribed and desired locations. Two mixed integer programming models for the satellite synthesis problem are presented. Four solution strategies, branch-and-bound, Benders' decomposition, linear programming with restricted basis entry, and a switching heuristic, are used to find solutions to example synthesis problems. Computational results indicate the switching algorithm yields solutions of good quality in reasonable execution times when compared to the other solution methods. It is demonstrated that the switching algorithm can be applied to synthesis problems with the objective of minimizing the largest deviation between a prescribed location and the corresponding desired location. Furthermore, it is shown that the switching heuristic can use no conservative, location-dependent satellite separations in order to satisfy interference criteria.

Gonsalvez, David J. A.↗

NASA NTRS1987

Research in computational fluid dynamics

The numerical integration of quasi-one-dimensional unsteady flow problems which involve finite rate chemistry are discussed, and are expressed in terms of conservative form Euler and species conservation equations. Hypersonic viscous calculations for delta wing geometries is also examined. The conical Navier-Stokes equations model was selected in order to investigate the effects of viscous-inviscid interations. The more complete three-dimensional model is beyond the available computing resources. The flux vector splitting method with van Leer's MUSCL differencing is being used. Preliminary results were computed for several conditions.

Murman, Earll M.↗

NASA NTRS1986

On symmetric and upwind TVD schemes

The performance of the upwind and symmetric total variation diminishing (TVD) schemes in viscous and inviscid airfoil steady-state calculations is considered, and the extension of the implicit second-order-accurate TVD scheme for hyperbolic systems of conservative laws in curvilinear coordinates is discussed. For two-dimensional steady-state applications, schemes are implemented in a conservative noniterative alternating direction implicit form, and results illustrate that the algorithm produces a fairly good solution for an RAE2822 airfoil calculation. The study demonstrates that the symmetric TVD scheme is as accurate as the upwind TVD scheme, while requiring less computational effort than it.

Yee, H. C.↗

NASA NTRS1987

A coupled mode formulation by reciprocity and a variational principle

A coupled mode formulation for parallel dielectric waveguides is presented via two methods: a reciprocity theorem and a variational principle. In the first method, a generalized reciprocity relation for two sets of field solutions satisfying Maxwell's equations and the boundary conditions in two different media, respectively, is derived. Based on the generalized reciprocity theorem, the coupled mode equations can then be formulated. The second method using a variational principle is also presented for a general waveguide system which can be lossy. The results of the variational principle can also be shown to be identical to those from the reciprocity theorem. The exact relations governing the 'conventional' and the new coupling coefficients are derived. It is shown analytically that the present formulation satisfies the reciprocity theorem and power conservation exactly, while the conventional theory violates the power conservation and reciprocity theorem by as much as 55 percent and the Hardy-Streifer (1985, 1986) theory by 0.033 percent, for example.

Chuang, Shun-Lien↗