Search NASA⌕ Search

SEARCH · Search NASA

Results for “parabolic problems”

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 307 records · Page 17

Comparison Between Navier-Stokes and Thin-Layer Computations for Separated Supersonic Flow

In the numerical simulation of high Reynolds-number flow, one can frequently supply only enough grid points to resolve the viscous terms in a thin layer. As a consequence, a body-or stream-aligned coordinate system is frequently used and viscous terms in this direction are discarded. It is argued that these terms cannot be resolved and computational efficiency is gained by their neglect. Dropping the streamwise viscous terms in this manner has been termed the thin-layer approximation. The thin-layer concept is an old one, and similar viscous terms are dropped, for example, in parabolized Navier-Stokes schemes. However, such schemes also make additional assumptions so that the equations can be marched in space, and such a restriction is not usually imposed on a thin-layer model. The thin-layer approximation can be justified in much the same way as the boundary-layer approximation; it requires, therefore, a body-or stream-aligned coordinate and a high Reynolds number. Unlike the boundary-layer approximation, the same equations are used throughout, so there is no matching problem. Furthermore, the normal momentum equation is not simplified and the convection terms are not one-sided differenced for marching. Consequently, the thin-layer equations are numerically well behaved at separation and require no special treatment there. Nevertheless, the thin-layer approximation receives criticism. It has been suggested that the approximation is invalid at separation and, more recently, that it is inadequate for unsteady transonic flow. Although previous comparisons between the thin-layer and Navier-Stokes equations have been made, these comparisons have not been adequately documented.

Degani, David↗

On the Theory of the Laval Nozzle

In the present paper, the motion of a gas in a plane-parallel Laval nozzle in the neighborhood of the transition from subsonic to supersonic velocities is studied. In a recently published paper, F. I. Frankl, applying the holograph method of Chaplygin, undertook a detailed investigation of the character of the flow near the line of transition from subsonic to supersonic velocities. From the results of Tricomi's investigation on the theory of differential equations of the mixed elliptic-hyperbolic type, Frankl introduced as one of the independent variables in place of the modulus of the velocity, a certain specially chosen function of this modulus. He thereby succeeded in explaining the character of the flow at the point of intersection of the transition line and the axis of symmetry (center of the nozzle) and in studying the behavior of the stream function in the neighborhood of this point by separating out the principal term having, together with its derivatives, the maximum value as compared with the corresponding corrections. This principal term is represented in Frankl's paper in the form of a linear combination of two hypergeometric functions. In order to find this linear combination, it is necessary to solve a number of boundary problems, which results in a complex analysis. In the investigation of the flow with which this paper is concerned, a second method is applied. This method is based on the transformation of the equations of motion to a form that may be called canonical for the system of differential equations of the mixed elliptic-hyperbolic type to which the system of equations of the motion of an ideal compressible fluid refers. By studying the behavior of the integrals of this system in the neighborhood of the parabolic line, the principal term of the solution is easily separated out in the form of a polynomial of the third degree. As a result, the computation of the transitional part of the nozzle is considerably simplified.

Falkovich, S. V.↗

Semi-discrete Galerkin solution of the compressible boundary-layer equations with viscous-inviscid interaction

A semi-discrete Galerkin (SDG) method is under development to model attached, turbulent, and compressible boundary layers for transonic airfoil analysis problems. For the boundary-layer formulation the method models the spatial variable normal to the surface with linear finite elements and the time-like variable with finite differences. A Dorodnitsyn transformed system of equations is used to bound the infinite spatial domain thereby providing high resolution near the wall and permitting the use of a uniform finite element grid which automatically follows boundary-layer growth. The second-order accurate Crank-Nicholson scheme is applied along with a linearization method to take advantage of the parabolic nature of the boundary-layer equations and generate a non-iterative marching routine. The SDG code can be applied to any smoothly-connected airfoil shape without modification and can be coupled to any inviscid flow solver. In this analysis, a direct viscous-inviscid interaction is accomplished between the Euler and boundary-layer codes through the application of a transpiration velocity boundary condition. Results are presented for compressible turbulent flow past RAE 2822 and NACA 0012 airfoils at various freestream Mach numbers, Reynolds numbers, and angles of attack.

Day, Brad A.↗

Assessment of a parabolic analysis for axisymmetric internal flows in rocket and turbomachinery ducts

The flow paths in gas turbine passages encompass a wide range of flow properties such as Reynolds number and Mach number as well as many other variable flow conditions such as swirl, free-stream turbulence, and laminar/turbulent transition. An existing computer program, the Axisymmetric Diffuser Duct (ADD) code, which calculates compressible, turbulent, swirling flow through axisymmetric ducts has been modified to include the effects of free-stream turbulence and laminar/turbulent transition. The program has been evaluated on a matrix of test cases to determine its accuracy, robustness, and limits of applicability. This improved version of the ADD code calculates solutions which compare well with available data and can now be applied to a wider range of problems than previously possible. In addition a concept called local enhancement was developed and tested on a simple two-dimensional geometry in order to demonstrate a method to reduce computer time. In this concept the pressure distribution is calculated on a coarse grid using the ADD code and the viscous layer is locally enhanced using a boundary layer analysis. By applying this concept, an order of magnitude reduction in computer time was possible without any loss in accuracy.

Power, G. D.↗

Leading-edge receptivity for blunt-nose bodies

This research program investigates boundary-layer receptivity in the leading-edge region for bodies with blunt leading edges. Receptivity theory provides the link between the unsteady distrubance environment in the free stream and the initial amplitudes of the instability waves in the boundary layer. This is a critical problem which must be addressed in order to develop more accurate prediction methods for boundary-layer transition. The first phase of this project examines the effects of leading-edge bluntness and aerodynamic loading for low Mach number flows. In the second phase of the project, the investigation is extended to supersonic Mach numbers. Singular perturbation techniques are utilized to develop an asymptotic theory for high Reynolds numbers. In the first year, the asymptotic theory was developed for leading-edge receptivity in low Mach number flows. The case of a parabolic nose is considered. Substantial progress was made on the Navier-Sotkes computations. Analytical solutions for the steady and unsteady potential flow fields were incorporated into the code, greatly expanding the types of free-stream disturbances that can be considered while also significantly reducing the the computational requirements. The time-stepping algorithm was modified so that the potential flow perturbations induced by the unsteady pressure field are directly introduced throughout the computational domain, avoiding an artificial 'numerical diffusion' of these from the outer boundary. In addition, the start-up process was modified by introducing the transient Stokes wave solution into the downstream boundary conditions.

Kerschen, Edward J.↗

Large space system control technology status and accomplishments

In the area of antenna modeling and control, preliminary structural models were defined for two representative parabolic reflectors. A control system design evolved for attitude control of the reflectors. The controller design was based on a lumped control concept where the control hardware (sensor and actuators) was mounted at the base of the antenna. In the area of distributed control, static shape control techniques were worked out to establish a prescribed vehicle shape. The corresponding estimation process for determination of vehicle shape from selected sensor measurements was also developed. A model order reduction study was conducted to find the best preflight dynamical models for on board controller design. The problems caused by truncation of the vehicle dynamics required to minimize on board computations were investigated. Estimator designs were developed for on board detection of large structure model errors.

Rodriguez, G.↗

Crashworthy Seats Would Afford Superior Protection

Seats to prevent or limit crash injuries to astronauts aboard the crew vehicle of the Orion spacecraft are undergoing development. The design of these seats incorporates and goes beyond crash-protection concepts embodied in prior spacecraft and racing-car seats to afford superior protection against impacts. Although the seats are designed to support astronauts in a recumbent, quasi-fetal posture that would likely not be suitable for non-spacecraft applications, parts of the design could be adapted to military and some civilian aircraft seats and to racing car seats to increase levels of protection. The main problem in designing any crashworthy seat is to provide full support of the occupant against anticipated crash and emergency-landing loads so as to safely limit motion, along any axis, of any part of the occupant s body relative to (1) any other part of the occupant s body, (2) the spacecraft or other vehicle, and (3) the seat itself. In the original Orion spacecraft application and in other applications that could easily be envisioned, the problem is complicated by severe limits on space available for the seat, a requirement to enable rapid egress by the occupant after a crash, and a requirement to provide for fitting of the seat to a wide range of sizes and shapes of a human body covered by a crash suit, space suit, or other protective garment. The problem is further complicated by other Orion-application-specific requirements that must be omitted here for the sake of brevity. To accommodate the wide range of crewmember body lengths within the limits on available space in the original Orion application, the design provides for taller crewmembers to pull their legs back closer toward their chests, while shorter crewmembers can allow their legs to stretch out further. The range of hip-support seat adjustments needed to effect this accommodation, as derived from NASA s Human Systems Integration Standard, was found to define a parabolic path along which the knees must be positioned. For a given occupant, the specific position along the path depends on the distance from the heel to the back of the knee. The application of the concept of parabolic adjustment of the hip-support structure caused the seat pan to also take on a parabolic shape, yielding the unanticipated additional benefit that the seat pan fits the occupant s buttocks and thighs more nearly conformally than do seat pans of prior design. This more nearly conformal fit effectively eliminates a void between the occupant s body and the seat pan, thereby helping to prevent what, in prior seat designs, was shifting of the occupant s body into that void during an impact.

Gohmert, Dustin↗

Numerical solutions of transonic flows by parametric differentiation and integral equation techniques

The paper discusses results of an exploratory study of the advantages obtained by combining integral equation methods with the method of parametric differentiation in the treatment of transonic flow problems. In the proposed method, the nonlinear unsteady transonic flow equation for small perturbations is transformed into a linear equation by the use of the method of parametric differentiation. The linear equation is split into a pair of weakly coupled partial differential equations by writing the transformed perturbation potential as the sum of a steady component and an unsteady component. The solution of the steady equation as an integral equation is based on Ogana's treatment (1978). As a test case, the formulation is applied to predict the steady transonic flow over a nonlifting parabolic-arc airfoil.

Sivaneri, N. T.↗

Application of CFD to sonic boom near and mid flow-field prediction

A 3-D parabolized Navier-Stokes (PNS) code was used to calculate the supersonic overpressures from three different geometries at near- and mid-flow fields. Wind tunnel data is used for code validation. Comparison of the computed results with different grid refinements is shown. It is observed that a large number of grid points is needed to resolve the tail shock/expansion fan interaction. Therefore, an adaptive grid approach is employed to calculate the flow field. The agreement between the numerical results and the wind tunnel data confirms that computational fluid dynamics can be applied to the problem of sonic boom prediction.

Cheung, Samson H.↗

The semi-discrete Galerkin finite element modelling of compressible viscous flow past an airfoil

A method is developed to solve the two-dimensional, steady, compressible, turbulent boundary-layer equations and is coupled to an existing Euler solver for attached transonic airfoil analysis problems. The boundary-layer formulation utilizes the semi-discrete Galerkin (SDG) method to model the spatial variable normal to the surface with linear finite elements and the time-like variable with finite differences. A Dorodnitsyn transformed system of equations is used to bound the infinite spatial domain thereby permitting the use of a uniform finite element grid which provides high resolution near the wall and automatically follows boundary-layer growth. The second-order accurate Crank-Nicholson scheme is applied along with a linearization method to take advantage of the parabolic nature of the boundary-layer equations and generate a non-iterative marching routine. The SDG code can be applied to any smoothly-connected airfoil shape without modification and can be coupled to any inviscid flow solver. In this analysis, a direct viscous-inviscid interaction is accomplished between the Euler and boundary-layer codes, through the application of a transpiration velocity boundary condition. Results are presented for compressible turbulent flow past NACA 0012 and RAE 2822 airfoils at various freestream Mach numbers, Reynolds numbers, and angles of attack. All results show good agreement with experiment, and the coupled code proved to be a computationally-efficient and accurate airfoil analysis tool.

Meade, Andrew J., Jr.↗

Status Of The UPS Space-Marching Flow Solver

The status of the three-dimensional parabolized Navier-Stokes solver UPS is described. The UPS code, initiated at NASA Ames Research Center in 1986, continues to develop and evolve through application to supersonic and hypersonic flow fields. Hypersonic applications have motivated enhancement of the physical modeling capabilities of the code, specifically real gas modeling, boundary conditions, and turbulence and transition modeling. The UPS code has also been modified to enhance robustness and efficiency in order to be practically used in concert with an optimization code for supersonic transport design. These developments are briefly described along with some relevant results for generic test problems obtained during verification of the enhancements. Included developments and results have previously been published and widely disseminated domestically.

Lawerence, Scott L.↗

Gas-Liquid Separation Strategies in Microgravity Environment

Bubble entrainment in liquids represents a serious problem in the microgravity environment. Whenever bubbles are entrained in a liquid,they tend to remain stationary in the liquid bulk in the absence of any external forcing. This is due to the reduction or complete absence of the buoyancy force in the microgravity environment, Thus the buoyancy force can not the be exploited to place the bubbles at the top of the liquid volume as in Ig(sub o) conditions. This situation represents a serious drawback in many space based engineering and scientific applications. We have demonstrated in a series of low gravity experiments conducted during parabolic flight on board aircraft that bubbles can be controlled in such a manner as to increase,the probability of their expulsion from a liquid bulk. In these tests the liquid'bulk was made either to be contained within, or to flow through specially designed containers using capillary force alone. Such containers appear to facilitate bubble removal, from the liquid bulk. Different successful liquid flow configurations will be discussed and the efficacy of the resulting bubble expulsion mechanisms will be demonstrated.

Antar, Basil N.↗

Combined Numerical/Analytical Perturbation Solutions of the Navier-Stokes Equations for Aerodynamic Ejector/Mixer Nozzle Flows

In spite of rapid advances in both scalar and parallel computational tools, the large number of variables involved in both design and inverse problems make the use of sophisticated fluid flow models impractical, With this restriction, it is concluded that an important family of methods for mathematical/computational development are reduced or approximate fluid flow models. In this study a combined perturbation/numerical modeling methodology is developed which provides a rigorously derived family of solutions. The mathematical model is computationally more efficient than classical boundary layer but provides important two-dimensional information not available using quasi-1-d approaches. An additional strength of the current methodology is its ability to locally predict static pressure fields in a manner analogous to more sophisticated parabolized Navier Stokes (PNS) formulations. To resolve singular behavior, the model utilizes classical analytical solution techniques. Hence, analytical methods have been combined with efficient numerical methods to yield an efficient hybrid fluid flow model. In particular, the main objective of this research has been to develop a system of analytical and numerical ejector/mixer nozzle models, which require minimal empirical input. A computer code, DREA Differential Reduced Ejector/mixer Analysis has been developed with the ability to run sufficiently fast so that it may be used either as a subroutine or called by an design optimization routine. Models are of direct use to the High Speed Civil Transport Program (a joint government/industry project seeking to develop an economically.viable U.S. commercial supersonic transport vehicle) and are currently being adopted by both NASA and industry. Experimental validation of these models is provided by comparison to results obtained from open literature and Limited Exclusive Right Distribution (LERD) sources, as well as dedicated experiments performed at Texas A&M. These experiments have been performed using a hydraulic/gas flow analog. Results of comparisons of DREA computations with experimental data, which include entrainment, thrust, and local profile information, are overall good. Computational time studies indicate that DREA provides considerably more information at a lower computational cost than contemporary ejector nozzle design models. Finally. physical limitations of the method, deviations from experimental data, potential improvements and alternative formulations are described. This report represents closure to the NASA Graduate Researchers Program. Versions of the DREA code and a user's guide may be obtained from the NASA Lewis Research Center.

DeChant, Lawrence Justin↗

Three-layer interactive method for computing supersonic laminar separated flows

An interactive model for numerical computation of complicated two-dimensional flowfields including regions of reversed flow is proposed. The present approach is one of dividing the flowfield into three regions, in each of which a simplified mathematical model is applied: (1) outer, supersonic flow for which the full potential equation (hyperbolic) is used; (2) viscous, laminar layer in which the compressible boundary-layer model (parabolic) is used; and (3) recirculating flow modeled by the incompressible Navier-Stokes equations (elliptic). For matching of the numerical solutions in the three layers, two interaction models are developed: one for pressure interaction, the other for interaction between the shear layer and the recirculating flow. The uniform solution for the whole flowfield is then obtained by iteration of the local solutions under the constraints imposed by matching. The three-layer interactive model is used for solution of the flowfield past an asymmetric cavity. The method is shown to be capable of dealing with backflow without encountering problems at separation, characteristic to the boundary-layer approach.

Brandeis, J.↗

Curved edge modification of compact range reflector

For many years, the compact range has successfully been used to design and evaluate many electromagnetic radiation and scattering systems. However, it has had limited use for large structures because of the large discrepancy between the reflector and target zone sizes. The limitation results from the large edge diffracted signal which emanates from the termination of the reflector. Previous attempts to solve this problem have involved using serrated edges to diffuse the diffracted field. While this reduces the edge diffracted field in the target zone, it is not eliminated. In addition, one has introduced many new corners which also diffract into the target zone. A curved edge modification is presented which reduces the edge diffraction by an order of magnitude or more and in the process does not create new mechanisms which perturb the plane wave in the target zone. Using this curved modification, one is able to design a compact range reflector whose target zone cross section approaches the size of the reflector parabolic surface. This makes the compact range ideal for scale model measurements since larger targets (antennas) can be measured. It also makes it more practical to consider systems which can measure full scale vehicles.

Burnside, Walter D.↗

Application of CFD to sonic boom near and mid flow-field prediction

A three-dimensional parabolized Navier-Stokes (PNS) code has been used to calculate the supersonic overpressures from three different geometries at near- and mid-flow fields. Wind-tunnel data is used for code validation. Comparison of the computed results with different grid refinements is shown in this paper. It is observed that a large number of grid points is needed to resolve the tail shock/expansion fan interaction. Therefore, an adaptive grid approach is employed to calculate the flow field. The agreement between the numerical results and the wind-tunnel data confirms that computational fluid dynamics can be applied to the problem of sonic boom prediction.

Cheung, Samson H.↗

A 3D-PNS computer code for the calculation of supersonic combusting flows

A computer code has been developed based on the three-dimensional parabolized Navier-Stokes (PNS) equations which govern the supersonic combusting flow of the hydrogen-air system. The finite difference algorithm employed was a hybrid of the Schiff-Steger algorithm and the Vigneron, et al., algorithm which is fully implicit and fully coupled. The combustion of hydrogen and air was modeled by the finite-rate two-step combustion model of Rogers-Chinitz. A new dependent variable vector was introduced to simplify the numerical algorithm. Robustness of the algorithm was considerably enhanced by introducing an adjustable parameter. The computer code was used to solve a premixed shock-induced combustion problem and the results were compared with those of a full Navier-Stokes code. Reasonably good agreement was obtained at a fraction of the cost of the full Navier-Stokes procedure.

Chitsomboon, Tawit↗

An Analysis of the Corridor and Guidance Requirements for Supercircular Entry into Planetary Atmospheres

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.

Chapman, Dean R.↗