Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear potential 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.

At least 271 records · Page 15

PROGRAM VSAERO: A computer program for calculating the non-linear aerodynamic characteristics of arbitrary configurations: User's manual

VSAERO is a computer program used to predict the nonlinear aerodynamic characteristics of arbitrary three-dimensional configurations in subsonic flow. Nonlinear effects of vortex separation and vortex surface interaction are treated in an iterative wake-shape calculation procedure, while the effects of viscosity are treated in an iterative loop coupling potential-flow and integral boundary-layer calculations. The program employs a surface singularity panel method using quadrilateral panels on which doublet and source singularities are distributed in a piecewise constant form. This user's manual provides a brief overview of the mathematical model, instructions for configuration modeling and a description of the input and output data. A listing of a sample case is included.

Maskew, B.↗

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗

International Space Station Bipropellant Plume Contamination Model Update for Short Thruster Pulse Widths

The International Space Station (ISS) Bipropellant Plume Contamination Model has been a vital tool for characterizing the thruster plume-induced contamination environment at the ISS but was not intended to be used for very short thruster pulse widths. This paper presents an updated model that ensures full start-up and shut-down phases are modeled for each thruster firing, when the majority of liquid phase contaminant mass is released. The updated ISS Bipropellant Plume Contamination Model prevents potential under-prediction of thruster plume-induced contamination due to visiting vehicle proximity operations and provides a way to take advantage of thruster start-up and shut-down data performance data gathered during thruster test programs, if available. The International Space Station (ISS) Bipropellant Plume Contamination Model developed by Boeing Space Environments is a semi-empirical model anchored in flight experiment data and has been a vital tool for characterizing the thruster plume-induced contamination environment at the ISS.[1] The current model utilizes flight experiment data from the Plume Impingement Contamination (PIC) and Shuttle Plume Impingement Flight Experiment (SPIFEX) studies, which include Orbiter 3870 N Primary Reaction Control System (PRCS) and Russian 130 N thrusters operating in pulse mode with 80-100 ms.[2] As the next generation of crew and cargo visiting vehicles are developed and arrive at ISS, minimum pulse widths of vehicle thrusters used for ISS proximity operations have decreased significantly below 80 ms. An update to the ISS Bipropellant Contamination Model is needed to prevent potential under-prediction of thruster plume-induced contamination for these very short pulse widths. Contamination due to thruster plumes occurs in the liquid phase (i.e., unburned or partially burned propellant in the plume). Liquid phase releases primarily occur during thruster start-up and shut-down phases, with the steady state phase contributing a small amount to the total contaminant mass released. The current ISS Bipropellant Plume Contamination Model includes functional dependencies on thruster parameters (thrust, mass flow rate, specific impulse (Isp)), distance to receiver surface, and angle off plume centerline. For modeling purposes, contaminant mass released scales linearly with pulse width. This is a conservative approach for pulse widths greater than 80 ms but could under-predict contamination for very short pulse widths (i.e., with little or no steady state phase). This paper presents a model update to add a functional dependency on commanded pulse width to ensure the portion of the pulse spent in start-up and shut-down phases is appropriately modeled, when the majority of liquid phase contaminant mass is released. The updated model prevents potential under-prediction of thruster plume-induced contamination due to visiting vehicle proximity operations while providing a way to take advantage of thruster start-up and shut-down data gathered during test programs, if available. Developing the plume model for a specific thruster using the ISS Bipropellant Plume Contamination Model must be done in consideration of all available thruster performance data. Example thruster data will be used to illustrate this point and options for implementing the model update with existing ISS Bipropellant Plume Contamination Model code will be discussed.

plume contamination↗

The effects of canard-wing flow-field interactions on longitudinal stability, effective dihedral and potential deep-stall trim

The literature available on high aspect ratio canard configurations shows them to have some unique stability characteristics. Using a generic canard-wing model, the effects of canard-wing flow-field interactions on stability were investigated in the NASA Langley Twelve-Foot Low-Speed Wind Tunnel. Results for the attached flow regime indicate linear interactions shift the neutral point of a canard configuration forward, but the effect of a canard on a wing can change significantly when the flow over the surface begins to separate, even several degrees below stall. The asymmetry of the canard downwash in a sideslip condition can result in an increment in effective dihedral roughly proportional to canard lift coefficient. At very high angles of attack the presence of a wing can cause an incremental normal force on a canard, contributing to the possibility of a deep-stall trim point. This effect is greater for a high canard and less for a low one.

Muchmore, C. B., Jr.↗

Extending A Correction Method for Unsteady Transonic Aerodynamics to Variable Camber Continuous Trailing Edge Flap

In this study, we have extended our correction method for transonic aerodynamics to baseline and variable camber continuous trailing edge (VCCTEF) geometries of super-critical airfoils at cruise conditions. This correction method modifies Theodorsen theory, which is based on incompressible potential flow over a flat plate. The method is partly based on CFD RANS simulations of oscillating NASA generic transport model (GTM) baseline cruise and VCCTEF airfoil geometries in pitch and shows promise for addressing transonic flutter problems of aeroelasticity. In the earlier study1, the Theodorsen functions were modulated with non-linear functions of Mach number, M, and reduced frequency, k, and thus the correction was demonstrated for the symmetrical NACA0012 airfoil. In the present study, the modulation functions are also allowed to account for the effects of camber and thickness of the airfoil. This is demonstrated by applying the correction method to various symmetric NACA00xx airfoil geometries (NACA001, NACA0002, NACA0012, NACA0015 and NACA0018), the NASA GTM super-critical baseline and VCCTEF airfoil geometries, e.g., the circular airfoil, VCCTEF222. First results indicate that the correction method accounts for a large variation of thickness of airfoils and camber. The present correction method will guide the development of a new state space model for the VCCTEF system and eventually a new transfer function that will be incorporated in a new aeroelastic framework leading to an appropriate transonic flutter model for use in the future aircraft systems in development under the NASA Advanced Air Transportation Technologies (AATT) project.

Kaul, Upender K.↗

Unsteady disturbances of streaming motions around bodies

The present investigation of flows around bodies with a stagnation point notes that the vortical part of such flows becomes singular along the entire body surface and its wake; as a result, the potential part of the flow will also be singular along the entire surface. A modified splitting of the velocity field is proposed which encompasses: (1) a vortical part having zero streamwise and normal components along the body surface; (2) an entropy-dependent part; and (3) a regular part similar to the potential part which satisfies a linear inhomogeneous wave equation with a modified source term.

Atassi, H. M.↗

Mode Spectra of Natural Disturbances in a Circular Jet and the Effect of Acoustic Forcing

A modal spectrum technique was used to study coherent instability modes triggered by naturally occurring disturbances in a circular jet. This technique was applied to a high Reynolds number jet for both untripped and tripped nozzle exit boundary layers, with both cases having a core turbulence level of 0.15 percent. The region up to the end of the potential core was dominated by the axisymmetric mode, with the azimuthal modes dominating further downstream. The growth of the azimuthal modes was observed closer to the nozzle exit for the jet with a transitional boundary layer. Whether for locally parallel flow or slowly diverging flow, even at low levels of acoustic forcing, the inviscid linear theory is seen to be inadequate for predicting the amplitude of the forced mode. In contrast, the energy integral approach reasonably predicts the evolution of the forced mode.

Raman, G.↗

Development of an unsteady aerodynamic analysis for finite-deflection subsonic cascades

An unsteady potential flow analysis, which accounts for the effects of blade geometry and steady turning, was developed to predict aerodynamic forces and moments associated with free vibration or flutter phenomena in the fan, compressor, or turbine stages of modern jet engines. Based on the assumption of small amplitude blade motions, the unsteady flow is governed by linear equations with variable coefficients which depend on the underlying steady low. These equations were approximated using difference expressions determined from an implicit least squares development and applicable on arbitrary grids. The resulting linear system of algebraic equations is block tridiagonal, which permits an efficient, direct (i.e., noniterative) solution. The solution procedure was extended to treat blades with rounded or blunt edges at incidence relative to the inlet flow.

Verdon, J. M.↗

Direct numerical simulations of turbulent convection with a variable gravity and Keplerian rotation

Thermal convection was proposed as a possible mechanism for generation and maintenance of turbulence in the inner accretion disk regime of the primordial solar nebula. It is of fundamental interest to design experiments with the basic physical features of the solar nebula conditions cannot be produced in the laboratory, numerical simulations of hydrodynamic flows, which have been very successful in describing aerodynamic flows, can be suitable modified to provide experimental data for solar nebula modelling. The goals are to modify an extant, proven hydrodynamics code with the most important features of the solar nebula and other thin accretion disks: bouyancy terms to generate convection, internal heating representing the release of gravitational potential energy, a variable gravity linearly proportional the the distance from the vertical midplane due to centrifugal balance, rapid rotation with axis aligned with gravity, and Keplerian rotational shear; to determine the effect that these features have on the turbulent convection by introducing them individually and to determine the cumulative nature of the turbulent convection for accretion disk conditions; and to model the convection and the turbulence. In this manner, prior solar nebula models can be tested and their deficiencies rectified.

Cabot, William H.↗

Laboratory and field evaluation of an underwater sand height gage

An underwater sand height gage was investigated. This instrument consisted of two transducers, one screened and one unscreened. Laboratory experimentation included static and dynamic tests with three test sands--fine, medium, and coarse. Field tests were conducted at Rudee Inlet, Virginia Beach, Virginia. Test results showed a linear response to up to 10 inches of sand loading. Deviation observed in identical tests appeared to be due to variation in the density of sand. Density differences reflected varying packing styles which, in turn, were a consequence of grain size and flow regime. Further evaluations of the sand height gage reflect this instrument's potential.

Swift, D. J. P.↗

A finite element algorithm for sound propagation in axisymmetric ducts containing compressible mean flow

The described investigation is concerned with he development of a finite element scheme which can be used in a study of the acoustics of aircraft-engine ducts. In the absence of suitable variational principles for acoustic fluctuations within an aircraft fan engine, an acoustic analysis must proceed directly from the differential equations which describe compressible flow. The derived equations cannot be solved algebraically. The numerical technique used for solving them makes use of a linear rectangular element of a type considered by Zienkiewicz (1971). Attention is given to aspects of element derivation, the global matrix assembly, the solution of the matrix equation, questions of acoustic attenuation, and illustrations of the potential of the current model in duct optimization.

Abrahamson, A. L.↗

Study on the application of NASA energy management techniques for control of a terrestrial solar water heating system

This paper addresses the potential for enhanced solar system performance through sophisticated control of the collector loop flow rate. Computer simulations utilizing the TRNSYS solar energy program were performed to study the relative effect on system performance of eight specific control algorithms. Six of these control algorithms are of the proportional type: two are concave exponentials, two are simple linear functions, and two are convex exponentials. These six functions are typical of what might be expected from future, more advanced, controllers. The other two algorithms are of the on/off type and are thus typical of existing control devices. Results of extensive computer simulations utilizing actual weather data indicate that proportional control does not significantly improve system performance. However, it is shown that thermal stratification in the liquid storage tank may significantly improve performance.

Swanson, T. D.↗

Distribution of wave drag and lift in the vicinity of wing tips at supersonic speeds

The point-source-distribution method of calculating the aerodynamic coefficients of thin wings at supersonic speeds was extended to include the effect of the region between the wing boundary and the foremost Mach wave from the wing leading edge. The effect of this region on the surface velocity potential has been determined by an equivalent function, which is evaluated over a portion of the wing surface. In this manner, the effect of angles of attack and yaw as well as the asymmetry of top and bottom wing surfaces may be calculated. As examples of the method, the pressure distribution on a thin plate wing of rectangular plan form as well as the lift and the drag coefficients as a function of Mach number, angle of attack, and aspect ratio are calculated. The equations for the surface velocity potential of several other plan forms are also included. (author)

AIRFOILS, FLAT-PLATE↗

Unsteady flow through compressor stages

The application of the nonsteady Lax-Wendroff technique to problems with asymptotically periodic solution which offers a potentially powerful method for the investigation of the interaction of rotating and stationary blade rows in turbomachinery is reported. A technique for specifying boundary conditions with phase lag was developed to accomplish this. A complete nonlinear analysis is carried out numerically to determine the entire flow field without recourse to the assumption of small disturbances of linear equations which underlie the previous acoustic theories. The result, obtained for the case of equal number of rotor and stator blades shows that transonic flow can be handled without difficulty. In addition, the program is not limited with regard to blade thickness, camber or loading. Extension of this method to incorporate viscous wakes and to analysis of fully three dimensional configuration is feasible, and would greatly expand its utility in practical applications.

Alzner, E.↗

On the mixing of a rectangular jet

Hot-wire measurements in an incompressible rectangular jet, issuing into a quiet environment at ambient conditions, are presented. A blow-down-type air supply system was used to provide the airflow to a cylindrical settling chamber 1.75 m in length and 0.6 m in diameter. The measurements were made with constant-temperature anemometers in conjunction with linearizers. The two signals from the linearizers were sent through a sum and difference unit which was calibrated from dc to 100 kHz. The distributions of mean velocity and the turbulence shear stresses were measured in the two central planes of the jet stations up to 115 widths downstream of the nozzle exit. Three distinct regions characterized the jet flow field: a potential core origin, a two-dimensional-type region, and an axisymmetric type region. The onset of the second region appeared to be at a location where the shear layers separated by the short dimension of the nozzle meet; and the third region occurred at a downstream location where the two shear layers from the short edges of the nozzle meet. In the central plane, similarity was found both in the mean velocity and shear stress profiles beyond 30 widths downstream of the nozzle exit; profiles of rms velocity showed similarity in the second, but not the third region.

Krothapalli, A.↗

Calculation of vortex flows on complex configurations

The calculation of aerodynamic characteristics of complex configurations having strongly coupled vortex flows is a non-linear problem requiring iterative solution techniques. This paper discusses the use of a low-order panel method as a means of obtaining practical solutions to such problems. The panel method is based on piecewise constant source and doublet quadrilateral panels and uses the internal Dirichlet boundary condition of zero perturbation potential. The problems of predicting vortex/surface interaction and vortex separation are discussed. Some example calculations are included but further test cases have yet to be carried out, in particular for comparisons with experimental data. The problem of convergence on the iterative calculation for the shape of the free vortex sheet is addressed and a preprocessor routine, based on an unsteady, two-dimensional version of the panel method, is put forward as a cost-effective way of generating an initial vortex structure for use as a starting solution for general configurations.

Maskew, B.↗

Water Atomization of Barium Fluoride: Calcium Fluoride for Enhanced Flow Characteristics of PS304 Feedstock Powder Blend

PS304 is a plasma spray deposited solid lubricant coating with feedstock composed of NiCr, Cr2O3, Ag, and BaF2-CaF2 powders. The effects of rounded BaF2-CaF2 particles on the gravity-fed flow characteristics of PS304 feedstock have been investigated. The BaF2-CaF2 powder was fabricated by water atomization using four sets of process parameters. Each of these powders was then characterized by microscopy and classified by screening to obtain 45 to 106 micron particles and added incrementally from 0 to 10 wt% to the other constituents of the PS304 feedstock, namely nichrome, chromia, and silver powders. The relationship between feedstock flow rate, measured with the Hall flowmeter, and concentration of fluorides was found to be linear in each case. The slopes of the lines were between those of the linear relationships previously reported using angular and spherical fluorides and were closer to the relationship predicted using the rule of mixtures. The results offer a fluoride fabrication technique potentially more cost-effective than gas atomization processes or traditional comminution processes.

Stanford, Malcolm K.↗

Assessment of Transition Modeling Capability in OVERFLOW with Emphasis on Swept-Wing Configurations

In preparation for comparisons with data obtained from the recently concluded experiments in the National Transonic Facility at the NASA Langley Research Center on the new common research model with natural laminar flow (CRM-NLF), an assessment of the current transition modeling capability in the NASA OVERFLOW 2.2o code has been carried out. A combination of the available experimental data and linear stability analysis is used to evaluate the accuracy and robustness of these models for selected swept-wing type configurations, with significant crossflow. An additional goal for this work involves providing a comparative assessment of the relevant transition models in the context of a single flow solver and identifying model limitations as well as the potential for future improvements that would help strengthen the physical basis of such transition models. Included in this investigation is an assessment of the sensitivities of the underlying transition models to grid resolution (wall-normal, as well as streamwise and spanwise) and the values of extra input parameters such as the level of surface roughness and freestream turbulence variables. The flow configurations targeted in this assessment include the NASA NLF(2)-0415 swept-wing configuration, the sickle-shaped wing introduced by the Technical University of Braunschweig, and the wing-body configuration of the CRM model from the fourth and fifth AIAA CFD Drag Prediction Workshops.

Freestream velocity↗