Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “compressible flows”

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 793 records · Page 44

Mars Helicopter Ingenuity Rotor Geometry

The NASA Mars Helicopter Ingenuity performed its maiden flight on April 19, 2021, and opened new possibilities for planetary exploration. Ingenuity has ended its mission after experiencing an anomaly during its 72nd flight on January 18, 2024. Due to the low density and low temperature of the surface atmosphere on Mars, significant challenges must be overcome to achieve aerodynamic performance. The low density of the Martian atmosphere and the relatively small-scale rotor result in flows with very low Reynolds number, reducing the lifting force and lifting efficiency, respectively. Until recently, compressible flow at Reynolds numbers around Re = 10,000 was of moderate practical interest but with the success of Ingenuity, the aerodynamic regime is appreciating considerable attention. The purpose of this report is to document the geometry (planform, twist, and airfoils) of the Mars Helicopter rotor, for use in future research on rotor performance in the Mars atmosphere.

Mars↗

Numerical simulation of flows in curved diffusers with cross-sectional transitioning using a three-dimensional viscous analysis

A three dimensional analysis for fully viscous, subsonic, compressible flow is evaluated. An approximate form of the Navier Stokes equations is solved by an implicit spatial marching technique. Calculations were made for flow in a circular S duct and in the F 16 inlet duct. The computed total pressure contours and secondary flow velocity vectors are presented. Qualitative comparisons with experiment are shown for both ducts. The analysis is used to show how the cross section transitioning in the F 16 inlet suppresses the development of a secondary flow vortex.

Towne, C. E.↗

General design method for 3-dimensional, potential flow fields. Part 2: Computer program DIN3D1 for simple, unbranched ducts

The general design method for three-dimensional, potential, incompressible or subsonic-compressible flow developed in part 1 of this report is applied to the design of simple, unbranched ducts. A computer program, DIN3D1, is developed and five numerical examples are presented: a nozzle, two elbows, an S-duct, and the preliminary design of a side inlet for turbomachines. The two major inputs to the program are the upstream boundary shape and the lateral velocity distribution on the duct wall. As a result of these inputs, boundary conditions are overprescribed and the problem is ill posed. However, it appears that there are degrees of compatibility between these two major inputs and that, for reasonably compatible inputs, satisfactory solutions can be obtained. By not prescribing the shape of the upstream boundary, the problem presumably becomes well posed, but it is not clear how to formulate a practical design method under this circumstance. Nor does it appear desirable, because the designer usually needs to retain control over the upstream (or downstream) boundary shape. The problem is further complicated by the fact that, unlike the two-dimensional case, and irrespective of the upstream boundary shape, some prescribed lateral velocity distributions do not have proper solutions.

Stanitz, J. D.↗

Accounting For Compressibility In Viscous Flow In Pipes

Method developed to account for effects of compressibility in viscous flows through long, circular pipes of uniform diameter. Based on approximation of variations in density and velocity across pipe cross section by profile equations developed for boundary-layer flow between flat plates.

Steinle, Frank W.↗

Wave-interactions in supersonic and hypersonic flows

Work completed under the current grant comprises the start of a theoretical and computational attack on the subharmonic route to secondary instabilities in compressible flows. The total flow field in this problem is made up of the following components: (1) a steady streamwise mean boundary layer flow which depends only on the normal space component y; (2) a two-dimensional time dependent T-S wave which moves with wavespeed c and has no spanwise dependence; and (3) a fully three-dimensional, time dependent T-S wave whose streamwise wavenumber is half of the streamwise wavenumber associated with the two-dimensional T-S wave in b. If a frame of reference is adopted which moves with the wavespeed c of the 2-D T-S wave, the time dependence of this portion of the flow can be eliminated. The effective steady mean flow in this problem is now the sum of the original parallel steady mean flow and the initial 2-D T-S instability. Dependence on the streamwise coordinate x in this mean flow can be extracted by assuming normal mode expansions involving complex exponentials and the streamwise wavenumber a. However, it is important to note that, because this is a wave-wave interaction problem, unlike the usual linear instability case, both the complex exponential, and its complex conjugate, must be retained in describing the 2-D T-S wave. The role of the perturbation to the steady mean flow is now played by the 3-D time dependent T-S wave. In treating this wave, normal modes in the streamwise and spanwise directions and time may be used. Consistent with the subharmonic nature of this transition route, the streamwise wavenumber is a/2, and complex conjugates of the complex exponential must be employed. This is not the case with the modes giving z and t dependence with wavespeed o and spanwise wavenumber B as the effective mean flow quantities are independent of z and their time dependence is accounted for by the moving frame of reference. Consequently, the wave-wave interaction which will produce mean flow modification occurs only through the streamwise exponentials.

Lakin, William D.↗

The liquid-conduction, vapor-flow heat pipe model in Sockeye

A single-phase heat pipe flow model implemented in the heat pipe application Sockeye is described. This model solves one-dimensional, compressible flow equations for the vapor phase in the center of a heat pipe, which are coupled to the two-dimensional heat conduction equation for the wick, liquid, and cladding, as well as to an ordinary differential equation tracking the working fluid inventory in the evaporator section of the heat pipe. This model is demonstrated with several test problems, including comparisons to analytic solutions for the vapor flow fields, analytic curves for sonic and capillary limitations of heat pipe operation, and some experimental data. The numerical solution gives excellent agreement for verification problems and good agreement with experimental results. Also, demonstrations show that the model is very robust, allowing for full simulations of heat pipe transients, including frozen startup, sonic-limited (supersonic) flow, and heat pipe shutdown.

42 ENGINEERING↗

A computational design tool for scramjet combustor film cooling and fuel mixing predictions

A methodology for the analysis of a film cooled combustor of a hydrogen-fueled, supersonic combustion ramjet (scramjet) is presented. The SHIP3D parabolic flow code, based on an extension of the pressure-correction method of Patankar and Spalding for compressible flow is used to compute the flow field in the combustor. The effect of film cooling on the overall vehicle performance is then calculated using the one-dimensional, multistep cycle code, CPIPE. Results are presented for a conceptual combustor at Mach 16 flight conditions. It is seen that film cooling can have a significant effect on the engine performance.

Kamath, Pradeep S.↗

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗

The effects of compressibility on dynamic stall

In this paper, typical computational predictions and experimental measurements of compressible flow past an airfoil at dynamic stall conditions are studied and compared to develop an insight into the effect of compressibility on dynamic stall. The dependency of the critical Mach number on airfoil leading edge curvature, camber, and angle of attack is investigated. Evidence is presented to show that a local region of supersonic flow occurs on an oscillating airfoil, even for a freestream Mach number as low as 0.2, if the boundary layer remains attached and the angle of attack is sufficiently high; that a shock terminates this local supersonic bubble; and that the vorticity that this shock generates grows rapidly and becomes very unstable as the angle of attack increases beyond the value at which the maximum local flow speed first exceeds the speed of sound. It is suggested that these shock-induced effects compete with the dynamic viscous effects occurring in the boundary layer in determining the onset of separation, which can lead to premature dynamic stall and can significantly reduce the maximum dynamic lift that can otherwise be obtained.

Fung, K.-Y.↗

Generalization and extension of the law of acoustic energy conservation in a nonuniform flow

An exact conservation equation is derived which generalizes the familiar acoustic energy equations. The new relation is valid for arbitrary disturbances to a viscous, compressible flow. It is suggested by a development of the acoustic energy equation by means of a regular perturbation expansion of the general energy equation of fluid mechanics. A perturbation energy density and flux are defined and identified as the exact physical quantities whose leading order perturbation representations are the usual acoustic energy density and flux. The conservation equation governing the perturbation energy quantities is shown to yield previously known results for several special cases.

Myers, M. K.↗

Efficient Parallel Algorithm For Direct Numerical Simulation of Turbulent Flows

A distributed algorithm for a high-order-accurate finite-difference approach to the direct numerical simulation (DNS) of transition and turbulence in compressible flows is described. This work has two major objectives. The first objective is to demonstrate that parallel and distributed-memory machines can be successfully and efficiently used to solve computationally intensive and input/output intensive algorithms of the DNS class. The second objective is to show that the computational complexity involved in solving the tridiagonal systems inherent in the DNS algorithm can be reduced by algorithm innovations that obviate the need to use a parallelized tridiagonal solver.

Moitra, Stuti↗

Control law synthesis and sensor design for active flutter suppression.

Methods are presented for representing unsteady aerodynamic loadings, valid for arbitrary motion, for the two-dimensional typical section with a trailing-edge control surface in incompressible flow and a three-dimensional lifting surface with leading- or trailing-edge control surfaces in subsonic compressible flow. Loads for the two-dimensional incompressible case are obtained via analytic continuation of the Theodorsen function into the complex plane, with guidance from a time-domain approximation of the Wagner indicial function. For the three-dimensional case, it is proposed that oscillatory generalized aerodynamic forces be approximated by Pade fractions, thereby permitting an obvious continuation into the complex plane. A theoretical justification for this procedure is briefly outlined, and some practical aspects of its implementation are discussed.

Lyons, M. G.↗

A method for calculating strut and splitter plate noise in exit ducts: Theory and verification

Portions of a four-year analytical and experimental investigation relative to noise radiation from engine internal components in turbulent flow are summarized. Spectra measured for such airfoils over a range of chord, thickness ratio, flow velocity, and turbulence level were compared with predictions made by an available rigorous thin-airfoil analytical method. This analysis included the effects of flow compressibility and source noncompactness. Generally good agreement was obtained. This noise calculation method for isolated airfoils in turbulent flow was combined with a method for calculating transmission of sound through a subsonic exit duct and with an empirical far-field directivity shape. These three elements were checked separately and were individually shown to give close agreement with data. This combination provides a method for predicting engine internally generated aft-radiated noise from radial struts and stators, and annular splitter rings. Calculated sound power spectra, directivity, and acoustic pressure spectra were compared with the best available data. These data were for noise caused by a fan exit duct annular splitter ring, larger-chord stator blades, and turbine exit struts.

Fink, M. R.↗

A finite element approach for solution of the 3D Euler equations

Prediction of thermal deformations and stresses has prime importance in the design of the next generation of high speed flight vehicles. Aerothermal load computations for complex three-dimensional shapes necessitate development of procedures to solve the full Navier-Stokes equations. This paper details the development of a three-dimensional inviscid flow approach which can be extended for three-dimensional viscous flows. A finite element formulation, based on a Taylor series expansion in time, is employed to solve the compressible Euler equations. Model generation and results display are done using a commercially available program, PATRAN, and vectorizing strategies are incorporated to ensure computational efficiency. Sample problems are presented to demonstrate the validity of the approach for analyzing high speed compressible flows.

Thornton, E. A.↗

Coupled flow, thermal and structural analysis of aerodynamically heated panels

A finite element approach to coupling flow, thermal and structural analyses of aerodynamically heated panels is presented. The Navier-Stokes equations for laminar compressible flow are solved together with the energy equation and quasi-static structural equations of the panel. Interactions between the flow, panel heat transfer and deformations are studied for thin stainless steel panels aerodynamically heated by Mach 6.6 flow.

Thornton, Earl A.↗

Coupled flow, thermal and structural analysis of aerodynamically heated panels

A finite element approach for coupling flow, thermal and structural analyses of aerodynamically heated panels is presented. The Navier-Stokes equations for laminar compressible flow are solved together with the energy equation and quasi-static structural equations of the panel. Interactions between the flow, panel heat transfer and deformations are studied for thin stainless steel panels aerodynamically heated by Mach 6.6 flow.

Thornton, Earl A.↗

Recent computation of viscous effects in transonic flow

The paper describes some improvements to a method of solving the time-dependent Reynolds-averaged Navier-Stokes equations for two-dimensional compressible flows for describing the coupling between viscous and inviscid regions in transonic flows. Mesh patching is implemented at the sublayer-boundary layer interfaces and the boundary layer-inviscid flow interface. Mesh stretching is geometric and is small in regions of large gradients and large in regions of small gradients, thus preserving second-order accuracy. Mesh lines are arranged so as to accommodate MacCormack's rapid solver for hyperbolic systems and a wider class of turbulent transport models. Results of computations of transonic flows over lifting airfoils are compared with experimental data and other solutions.

Deiwert, G. S.↗