A Fortran computer code for the calculation of unsteady supersonic flow field by the method of perturbation
Fortran computer code for analysis of unsteady supersonic flow field around oscillating body of revolution
Hsu, P.-T.↗
SEARCH · Search NASA
Results for “Unsteady Flow”
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.
Spatial adaption procedures on unstructured meshes for accurate unsteady aerodynamic flow computation
Spatial adaption procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaption procedures were developed and implemented within a two-dimensional unstructured-grid upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in a high gradient region or the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational costs. A detailed description is given of the enrichment and coarsening procedures and comparisons with alternative results and experimental data are presented to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady transonic results, obtained using spatial adaption for the NACA 0012 airfoil, are shown to be of high spatial accuracy, primarily in that the shock waves are very sharply captured. The results were obtained with a computational savings of a factor of approximately fifty-three for a steady case and as much as twenty-five for the unsteady cases.
Rausch, Russ D.↗
Spatial adaption procedures on unstructured meshes for accurate unsteady aerodynamic flow computation
Spatial adaption procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaption procedures were developed and implemented within a two-dimensional unstructured-grid upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in high gradient regions of the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational cost. The paper gives a detailed description of the enrichment and coarsening procedures and presents comparisons with alternative results and experimental data to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady transonic results, obtained using spatial adaption for the NACA 0012 airfoil, are shown to be of high spatial accuracy, primarily in that the shock waves are very sharply captured. The results were obtained with a computational savings of a factor of approximately fifty-three for a steady case and as much as twenty-five for the unsteady cases.
Rausch, Russ D.↗
Spatial adaptation procedures on tetrahedral meshes for unsteady aerodynamic flow calculations
Spatial adaptation procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaptation procedures were developed and implemented within a three-dimensional, unstructured-grid, upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in high gradient regions of the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational cost. A detailed description of the enrichment and coarsening procedures are presented and comparisons with experimental data for an ONERA M6 wing and an exact solution for a shock-tube problem are presented to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady results, obtained using spatial adaptation procedures, are shown to be of high spatial accuracy, primarily in that discontinuities such as shock waves are captured very sharply.
Rausch, Russ D.↗
Spatial adaptation procedures on tetrahedral meshes for unsteady aerodynamic flow calculations
Spatial adaptation procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaptation procedures were developed and implemented within a three-dimensional, unstructured-grid, upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in high gradient regions of the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational cost. The paper gives a detailed description of the enrichment and coarsening procedures and presents comparisons with experimental data for an ONERA M6 wing and an exact solution for a shock-tube problem to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady results, obtained using spatial adaptation procedures, are shown to be of high spatial accuracy, primarily in that discontinuities such as shock waves are captured very sharply.
Rausch, Russ D.↗
An alternating-direction implicit algorithm for unsteady potential flow
An implicit finite-difference scheme is presented for the efficient computation of unsteady potential flow about airfoils. The formulation uses density and the velocity potential as dependent variables and is cast in conservation form to assure the theoretically correct determination of shockwave location and speed. To enable boundary conditions to be imposed directly on the airfoil surface, a time-varying sheared-rectilinear coordinate transformation is employed. Calculated time-history solutions on a pulsating airfoil are compared with the results of other unsteady transonic codes, including a previous method of the authors. The present method is demonstrated to be unconditionally stable and to give accurate solutions with sharply resolved shocks.
Chipman, R.↗
Efficient solution of unsteady transonic flows about airfoils
An implicit finite difference procedure was developed for the efficient solution of unsteady transonic flow fields. Sample computations illustrate applications of procedures to aerodynamic problems. Solutions are presented that illustrate three types of shock wave motion that can result from airfoil control surface oscillations. The significant effect of wind tunnel wall conditions on these shock wave motions is demonstrated. Solutions are also presented for a simple aeroelastic problem in which the flow field equations and the structural motion equations are integrated simultaneously in time. Both stable and unstable aeroelastic interactions are considered. The procedure is adapted to compute unsteady aerodynamic force coefficients by the indicial method.
Ballhaus, W. F.↗
A harmonic analysis method for unsteady transonic flow and its application to the flutter of airfoils
A finite difference method for solving the unsteady transonic flow about harmonically oscillating wings is investigated. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. The differential equation for the unsteady velocity potential is linear with spatially varying coefficients and with the time variable eliminated by assuming harmonic motion. A study is presented of the shock motion associated with an oscillating airfoil and its representation by the harmonic procedure. The effects of the shock motion and the resulting pressure pulse are shown to be included in the harmonic pressure distributions and the corresponding generalized forces. Analytical and experimental pressure distributions for the NACA 64A010 airfoil are compared for Mach numbers of 0.75, 0.80 and 0.842. A typical section, two-degree-of-freedom flutter analysis of a NACA 64A010 airfoil is performed. The results show a sharp transonic bucket in one case and abrupt changes in instability modes.
Ehlers, F. E.↗
Discrete Adjoint-Based Design for Unsteady Turbulent Flows On Dynamic Overset Unstructured Grids
A discrete adjoint-based design methodology for unsteady turbulent flows on three-dimensional dynamic overset unstructured grids is formulated, implemented, and verified. The methodology supports both compressible and incompressible flows and is amenable to massively parallel computing environments. The approach provides a general framework for performing highly efficient and discretely consistent sensitivity analysis for problems involving arbitrary combinations of overset unstructured grids which may be static, undergoing rigid or deforming motions, or any combination thereof. General parent-child motions are also accommodated, and the accuracy of the implementation is established using an independent verification based on a complex-variable approach. The methodology is used to demonstrate aerodynamic optimizations of a wind turbine geometry, a biologically-inspired flapping wing, and a complex helicopter configuration subject to trimming constraints. The objective function for each problem is successfully reduced and all specified constraints are satisfied.
Nielsen, Eric J.↗
Prediction of unsteady transonic flow around missile configurations
This paper describes the preliminary development of a method for predicting the unsteady transonic flow around missiles at transonic and supersonic speeds, with the final goal of developing a computer code for use in aeroelastic calculations or during maneuvers. The basic equations derived for this method are an extension of those derived by Klopfer and Nixon (1989) for steady flow and are a subset of the Euler equations. In this approach, the five Euler equations are reduced to an equation similar to the three-dimensional unsteady potential equation, and a two-dimensional Poisson equation. In addition, one of the equations in this method is almost identical to the potential equation for which there are well tested computer codes, allowing the development of a prediction method based in part on proved technology.
Nixon, D.↗
The numerical simulation of multistage turbomachinery flows
The effect of the unsteady flow field in a multistage compressor on the time-averaged performance was assessed. The energy transport by the unsteady deterministic flow field was taken into account. The magnitude of the body force resulting from the aerodynamic loading on a blade row was compared to the gradient of the stress tensor associated with the unsteady time-resolved flow field generated by the blade row. The magnitude of the work performed by these forces was compared to the divergence of the energy correlations produced by the unsteady time-resolved flow field. The stress tensor and the energy correlations are non-negligible in the end wall regions. The results suggest that the turbulence is the primary source of flow mixing in the midspan region.
Adamczyk, John J.↗
Unsteady MHD flow in arbitrary ducts by point-matching method.
Unsteady MHD conducting fluid flow in arbitrary ducts under transverse magnetic field, studying effects of pressure gradient impulsive change
Jeng, D. R.↗
Control of unsteady separated flow associated with the dynamic stall of airfoils
An effort to understand and control the unsteady separated flow associated with the dynamic stall of airfoils was funded for three years through the NASA cooperative agreement program. As part of this effort a substantial data base was compiled detailing the effects various parameters have on the development of the dynamic stall flow field. Parameters studied include Mach number, pitch rate, and pitch history, as well as Reynolds number (through two different model chord lengths) and the condition of the boundary layer at the leading edge of the airfoil (through application of surface roughness). It was found for free stream Mach numbers as low as 0.4 that a region of supersonic flow forms on the leading edge of the suction surface of the airfoil at moderate angles of attack. The shocks which form in this supersonic region induce boundary-layer separation and advance the dynamic stall process. Under such conditions a supercritical airfoil profile is called for to produce a flow field having a weaker leading-edge pressure gradient and no leading-edge shocks. An airfoil having an adaptive-geometry, or dynamically deformable leading edge (DDLE), is under development as a unique active flow-control device. The DDLE, formed of carbon-fiber composite and fiberglass, can be flexed between a NACA 0012 profile and a supercritical profile in a controllable fashion while the airfoil is executing an angle-of-attack pitch-up maneuver. The dynamic stall data were recorded using point diffraction interferometry (PDI), a noninvasive measurement technique. A new high-speed cinematography system was developed for recording interferometric images. The system is capable of phase-locking with the pitching airfoil motion for real-time documentation of the development of the dynamic stall flow field. Computer-aided image analysis algorithms were developed for fast and accurate reduction of the images, improving interpretation of the results.
Wilder, M. C.↗
Aeroacoustic sensitivity analysis and optimal aeroacoustic design of turbomachinery blades
During the first year of the project, we have developed a theoretical analysis - and wrote a computer code based on this analysis - to compute the sensitivity of unsteady aerodynamic loads acting on airfoils in cascades due to small changes in airfoil geometry. The steady and unsteady flow though a cascade of airfoils is computed using the full potential equation. Once the nominal solutions have been computed, one computes the sensitivity. The analysis takes advantage of the fact that LU decomposition is used to compute the nominal steady and unsteady flow fields. If the LU factors are saved, then the computer time required to compute the sensitivity of both the steady and unsteady flows to changes in airfoil geometry is quite small. The results to date are quite encouraging, and may be summarized as follows: (1) The sensitivity procedure has been validated by comparing the results obtained by 'finite difference' techniques, that is, computing the flow using the nominal flow solver for two slightly different airfoils and differencing the results. The 'analytic' solution computed using the method developed under this grant and the finite difference results are found to be in almost perfect agreement. (2) The present sensitivity analysis is computationally much more efficient than finite difference techniques. We found that using a 129 by 33 node computational grid, the present sensitivity analysis can compute the steady flow sensitivity about ten times more efficiently that the finite difference approach. For the unsteady flow problem, the present sensitivity analysis is about two and one-half times as fast as the finite difference approach. We expect that the relative efficiencies will be even larger for the finer grids which will be used to compute high frequency aeroacoustic solutions. Computational results show that the sensitivity analysis is valid for small to moderate sized design perturbations. (3) We found that the sensitivity analysis provided important insight into how airfoils should be modified to improve their aeroelastic stability. Using this insight, we redesigned an aeroelastically unstable cascade to produce a stable cascade.
Hall, Kenneth C.↗
Unsteady Turbopump Flow Simulations
The objective of the current effort is two-fold: 1) to provide a computational framework for design and analysis of the entire fuel supply system of a liquid rocket engine; and 2) to provide high-fidelity unsteady turbopump flow analysis capability to support the design of pump sub-systems for advanced space transportation vehicle. Since the space launch systems in the near future are likely to involve liquid propulsion system, increasing the efficiency and reliability of the turbopump components is an important task. To date, computational tools for design/analysis of turbopump flow are based on relatively lower fidelity methods. Unsteady, three-dimensional viscous flow analysis tool involving stationary and rotational components for the entire turbopump assembly has not been available, at least, for real-world engineering applications. Present effort is an attempt to provide this capability so that developers of the vehicle will be able to extract such information as transient flow phenomena for start up, impact of non-uniform inflow, system vibration and impact on the structure. Those quantities are not readily available from simplified design tools. In this presentation, the progress being made toward complete turbo-pump simulation capability for a liquid rocket engine is reported. Space Shuttle Main Engine (SSME) turbo-pump is used as a test case for the performance evaluation of the hybrid MPI/Open-MP and MLP versions of the INS3D code. Relative motion of the grid system for rotor-stator interaction was obtained by employing overset grid techniques. Time-accuracy of the scheme has been evaluated by using simple test cases. Unsteady computations for SSME turbopump, which contains 106 zones with 34.5 Million grid points, are currently underway on Origin 2000 systems at NASA Ames Research Center. Results from these time-accurate simulations with moving boundary capability and the performance of the parallel versions of the code will be presented.
Centin, Kiris C.↗
A comparison of two incompressible Navier-Stokes algorithms for unsteady internal flow
A comparative study of two different incompressible Navier-Stokes algorithms for solving an unsteady, incompressible, internal flow problem is performed. The first algorithm uses an artificial compressibility method coupled with upwind differencing and a line relaxation scheme. The second algorithm uses a fractional step method with a staggered grid, finite volume approach. Unsteady, viscous, incompressible, internal flow through a channel with a constriction is computed using the first algorithm. A grid resolution study and parameter studies on the artificial compressibility coefficient and the maximum allowable residual of the continuity equation are performed. The periodicity of the solution is examined and several periodic data sets are generated using the first algorithm. These computational results are compared with previously published results computed using the second algorithm and experimental data.
Wiltberger, N. Lyn↗
Axial flow rotor unsteady response to circumferential inflow distortions
The unsteady response of an axial flow fan rotor to steady, circumferential inflow velocity and stagnation pressure distortions is assessed by two different methods. These are: (1) investigation of the unsteady normal force and pitching moment on a chordwise element of a rotor blade; and (2) investigation of the variation of the stagnation pressure distortion between the inlet and exit of the rotor. Experimental measurements of these unsteady characteristics are presented as a function of the geometry of the rotor -- stagger angle, solidity and steady angle of incidence -- for sinusoidally varying circumferential distortions with different numbers of distortion cycles. These measurements are compared with several theoretical analyses. While these comparisons indicate some of the deficiencies which exist in the theories, the existence of an unsteady cascade effect and the ability of the theories to adequately predict the trend of the unsteady response due to variations in reduced frequency, rotor stagger angle, solidity and mean incidence angle is clearly demonstrated.
Bruce, E. P.↗
Applications of the Method of Space-Time Conservation Element and the Solution Element to Unsteady Chemically Reactive Flows
This document reports the conclusion and findings of our research activities for this grant. The goal of the project is the development and application of the method of Space-Time Conservation Element and Solution Element, or the CE/SE method, to simulate chemically reacting flows. The product of this project will be a high-fidelity, time-accurate flow solver analyzing unsteady flow fields advanced propulsion concepts, including the low-emission turbojet engine combustion and flow fields of the Pulse Detonation Engines (PDE). Based on the documents and computer software of the CE/SE method that we have received from the CE/SE working group at NASA Lewis, we have focused our research effort on addressing outstanding technical issues related to the extension of the CE/SE method for unsteady, chemically reactive flows. In particular, we have made progresses in the following three aspects: (1) Derivation of the governing equations for reacting flows; (2) Numerical treatments of stiff source terms; and (3) Detailed simulations of ZND detonation waves.
Yu, Sheng-Tao↗